problem about pressure driven flow
Dear ALL Foamers:
I have a very simple problem but is very difficult for me.
I am simulating a pressure driven compressible flow, ie, the inlet and outlet pressure are given, and velocity is calculated from the inlet and outlet conditions. but from the tutorial cases from compressible solver, velocity is given and pressure is calculated.
Do any of you have any ideas about how to set the boundary condition for pressure driven flow, OR any solver is suitable for this problem?
thank you in advance.
I think velocity changes and pressure changes are related. So in solving for one, you are left with an unknown constant offset. That is why we need to fix a reference pressure when it is not set at a boundary condition, and I think the same is true for velocity. If the "inlet" and "outlet" are zero gradient for velocity, there is a continuum of possible flow solutions that give the same pressure distribution, even with no-slip conditions on the walls. In other words, your solution may "blow up" with the velocity ever increasing from inlet to outlet for no apparent reason.
I had good results with the following settings for compressible, pressure driven flows:
p - TotalPressure
U - pressureInletVelocity
p - fixedValue (static pressure)
U - pressureInletOutletVelocity
I think, every compressible solver should be okay. I would start with rhoSimpleFoam for stationary flows and/or rhoPisoFoam for transient flows.
Are these BCs ok for an incompressible solver as well?
Test case ----> New case
p - zeroGradient ----> totalPressure where p0 = value from test case
U - flowRateInletVelocity ----> pressureInletVelocity
p - fixedValue ----> fixedValue (static pressure)
U - inletOutlet ----> pressureInletOutletVelocit
When I do this, the code bombs very quickly-2 iterations. I'm just starting to look at this more closely but I'm still clueless . The only thing I've tried beyond that is run slip walls on the duct, specify a small Pt, and turn off the turbulence model. I'm trying to make the simplest case possible to figure this out. Any tips would be appreciated.
|All times are GMT -4. The time now is 20:16.|