Adapting PISO for condition div(fracFluid*U) = 0
I'm studying the heat transfert and water flow in a nuclear pool.
The model i use for nuclear waste is porosity. So, i have introduced a variable "fracFluid" that represent the fraction of fluid.
So, i have 2 regions in my pool : the porous medium with phi = 0.57, and the water above with phi = 1.
Then, my continuity equation is not anymore div(U) = 0, but div(fracFluid*U) =0.
I've done some modification on my PISO loop, but it doesn't work.
Do you have any ideas why it's not good ?
here is my code :
Thanks a lot :-)
Hi everyone !
I solved my problem. I post the solution since it may be usefull for some of you...
The trick was to make interpolation(fracFluid*rAU) instead of interpolation(fracFluid)*interpolation(rAu) and so on for all variables.
here is the new equation of p_rgh :
Why is it U = HbyA + rAU*fvc::reconstruct((phig - p_rghEqn.flux())/rAUf)
instead of U = HbyA - rAU*fvc::reconstruct((phig - p_rghEqn.flux())/rAUf);
|All times are GMT -4. The time now is 12:21.|