kEpsilon and pisoFoam: when transport equations for k and eps are solved
Hello everyone!
Looking at pisoFoam solver, I am seeing now that, if I am not mistaken, the transport equations for the turbulent kinetic energy (k) and dissipation rate (epsilon) are solved outside piso algorithm loops, i.e., the sentence turbulence->correct() is called upon only when the velocity and pressure fields have already converged. This implies that the equations are written using the nut field computed in the previous time step, and when using wall functions, the term kp used is also from the previous time step. Is this correct or there is something I am not understanding correctly? Could the transport equations be solved inside the loop and checked for convergence too, as in Verteeg and Malasakera? |
Then go to pimpleFoam and adjust the nOuterCorrectors and the residualControl in the fvSolution
|
Hi,
You should also set turbOnFinalIterOnly to no in PIMPLE dictionary, so turbulent properties transport equations are solved on each outer iteration. Otherwise it will not differ much from pisoFoam. |
Many thanks for your responses, this is clearer now.
However, I have tried that: set nOuterCorrectors to 2 in fvSolutions, and turbOnFinalIterOnly to no in controlDict: Code:
turbOnFinalIterOnly no; Code:
PIMPLE Code:
Courant Number mean: 0.673517 max: 3.12374 P.D.1: ssss comments on adjusting the residualControl in fvSolution. Is there a specific residualControl for this? P.D.2. By the way, another related question. I am trying to turn on the debug switch in pimpleControl.C that activates this piece of code: Code:
if (debug) Code:
defineTypeNameAndDebug(pimpleControl, 0); Many thanks, |
Hi,
1. turbOnFinalIterOnly goes to PIMPLE dictionary in fvSolution file (just like nOuterCorrectors). 2. Either you can use forum search with keyword "residualControl", or "grep -r residualControl *" in tutorial folder. Unfortunately there are lots of examples for SIMPLE family of solvers, yet only one for PIMPLE (in compressible/rhoPimpleFoam/ras/angledDuct/system/fvSolution). Though even the example is from compressible solver, it is quite relevant to your case, you PIMPLE dictionary in fvSolution file should be something like: Code:
PIMPLE Code:
DebugSwitches |
Thanks! I was writting turbOnFinalIterOnly in the wrong place.
I find it strange that the turbulence transport equations are only solved at the end of the loop in pisoFoam, so that the PISO algorithm works with the nut from the previous time step. Of course, it is a matter of modifying it or using pimpleFoam. Thanks! |
All times are GMT -4. The time now is 10:30. |