shallowWaterFoam using PIMPLE without tmp and h.relax()
 

dear foamers.

i am grateful that people are interested in shallowWaterFoam or the PIMPLE algo. but in OF-2.3.0, i got some questions about the shallowWaterFoam implemented with the PIMPLE.

(1) why the momentumPredictor hUEqn is not declared by tmp<fvVectorMatrix> and later cleared if pimple.nCorrPISO() <=1, as the UEqn in pimpleFoam?

(2) why the water depth h is not relaxed after the hEqn is solved, as the pEqn in pimpleFoam? can h.relax() be added and what's the effect on the result?

(3) since there is a temporal derivative in the hEqn, what makes it different with the general pressure equation (say pEqn in pimpleFoam). what does the hEqn.flux() do and why phi = phiHbyA + hEqn.flux() makes the fluxes conservative? or what is the conservative face flux for the hEqn with a temporal derivative?

if the above questions are related to (or depend on) what i intend to solve, please give me some hints if there is a general manner.

thanks in advance.

best,

karelke

shallowWaterFoam using PIMPLE without tmp and h.relax()
 

i've been struggling on these problems for several days. and i guess the PIMPLE may be some tricky for me.

so if anyone also works on the PIMPLE or the shallowWaterFoam, please comment, please help.

any advices or hints will make my day. thanks a lot.

/karelke