|
[Sponsors] |
Adapting PISO for condition div(fracFluid*U) = 0 |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
April 11, 2014, 09:10 |
Adapting PISO for condition div(fracFluid*U) = 0
|
#1 | |
New Member
Jean-François Balfroid
Join Date: Oct 2013
Posts: 2
Rep Power: 0 |
Hello everyone,
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 : Quote:
Thanks a lot :-) Shoji |
||
April 15, 2014, 07:34 |
problem solved
|
#2 | |
New Member
Jean-François Balfroid
Join Date: Oct 2013
Posts: 2
Rep Power: 0 |
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 : Quote:
Shoji |
||
June 30, 2015, 12:08 |
|
#3 |
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 838
Rep Power: 17 |
Hell Jean,
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); Thanks.
__________________
My OpenFOAM algorithm website: http://dyfluid.com By far the largest Chinese CFD-based forum: http://www.cfd-china.com/category/6/openfoam We provide lots of clusters to Chinese customers, and we are considering to do business overseas: http://dyfluid.com/DMCmodel.html |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Non-linearity Pressure Equation -- PISO algorithm | gdeneyer | OpenFOAM Programming & Development | 1 | August 23, 2012 05:19 |
Development of a Low mach PISO solver | nishant_hull | OpenFOAM Programming & Development | 0 | August 25, 2009 12:48 |
SIMPLE and PISO | Aditya | Main CFD Forum | 5 | April 1, 2006 18:52 |
PISO - new questions | Benedikt Flurl | Main CFD Forum | 0 | May 7, 2005 12:02 |
PISO vs. SIMPLE | benedikt flurl | Main CFD Forum | 2 | April 14, 2005 06:54 |