Is it correct to initialize this here and not before?

unicfdlab / hybridCentralSolvers

United collection of hybrid Central solvers - one-phase, two-phase and multicomponent versions

GNU General Public License v3.0

95 stars 59 forks source link

Is it correct to initialize this here and not before? #1

Closed ChrisS85 closed 5 years ago

ChrisS85 commented 7 years ago

Question concerning https://github.com/unicfdlab/hybridCentralSolvers/blob/c3ecfa63041d3ee3ecadf8c49226f3196e2d48e6/OpenFOAM-4.1/pisoCentral/commonIncludes/createCommonCentralFields.H#L122

psiU_nei depends on psi_nei, should it not be affected by this correction? Just wondering while looking through the code...

unicfdlab commented 7 years ago

The function correctCentralACMIInterpolation ensures that weighting field passed as the argument to this function (psi_nei in this part of the code) is not zero because to calculate U_nei we use expression psiU_nei / psi_nei