Implementation of a porous domain
I am simulating reactions through a porous domain using porousExplicitSourceReactingParcelFoam. The header of the .C file of this solver reads:
Note: ddtPhiCorr not used here when porous zones are active
- not well defined for porous calculations
I guess this means that the solver is not well defined for porous domains because of the given reason, not the opposite (i.e. ddtPhiCorr is not well defined for porous domains, reason why it is not used here), right?
If so, my question is: why does ddtPhiCorr make the solver ill-defined and how can it be improved? I guess the programmers would have done it if simply using ddtPhiCorr was the answer...
Ok, I figured out what it means by opening the corresponding pEqn.H file.
However, I still do not know what the REASON for ddtPhiCorr not to be well defined is.
Furthermore, IŽd be glad to know if the fact of not using this function leads to the "right" solution or simply avoids the numerical concern but returns a solution that is not as good as it could be.
|All times are GMT -4. The time now is 10:22.|