CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Compressible air flow with large-constant pressure at inlet

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes
  • 1 Post By Alczem

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 4, 2023, 13:28
Default Compressible air flow with large-constant pressure at inlet
  #1
Member
 
Morteza
Join Date: Jan 2018
Posts: 30
Rep Power: 8
mortezahdr is on a distinguished road
Hello all,
I am trying to simulate the air flow entering a tube shown in the following image. The tube has some small holes from which the flow exits, but since the flow field outside the tube is needed, computational domain is much larger than the tube and the outlet is located far from it. The tube boundary is called airbar and the boundaries not labled in the image are called walls.

The only numeric boundary condition I have is 70 psi at inlet (=482600 pa) and because of that, I am using rhoPimpleFoam, as it obviously is a compressible flow.
The problem is that the solution does not become stable and if I use a table to gradually increase the inlet pressure from environment (101325 pa) to 482600 pa, it helps stabilizing the solution, but it gets a very small time step of order of 10^-7 or 10^-8, which is useless.
Is there any solution (schemes/turbulence model for example) to increase the time step? Or do you think it is possible at all to solve a flow with such a large inlet pressure with OpenFOAMor not? I use LES turbulence model and the average y+ on the tube is around 1. Boundary conditions for U and p are shown below:

U:

boundaryField
{
inlet
{
type pressureDirectedInletVelocity;
phi phi;
rho rho;
inletDirection uniform (0 0 -1);
value uniform (0 0 -1);
}

outlet
{
type zeroGradient;
}

airbar
{
type noSlip;
}

walls
{
type noSlip;
}

}

p:
dimensions [1 -1 -2 0 0 0 0];

internalField uniform 101325;

boundaryField
{
inlet
{
type uniformTotalPressure;
p0 tableFile;
file "data_p";
gamma 1.4;
}

outlet
{
type fixedValue;
value uniform 101325;
}

airbar
{
type zeroGradient;
}

walls
{
type zeroGradient;
}
}

Last edited by mortezahdr; April 4, 2023 at 13:32. Reason: Change image, add more text
mortezahdr is offline   Reply With Quote

Old   April 6, 2023, 04:13
Default
  #2
Senior Member
 
Join Date: Dec 2021
Posts: 211
Rep Power: 5
Alczem is on a distinguished road
Hey,


Can you do a steady state simulation with rhoSimpleFoam? It will be faster to achieve a final state. And you could use this to initialize a transient simulation with an already developped flow.


Otherwise, if the steadystate solver is unstable, you might want to give Local Time Stepping a try, by using a ddt scheme like localEuler or SLTS.


If you need a full transient simulation, I don't think you have a choice about the small timesteps are you sure your velocities are correct? You could also do a steady state simulation of the interior of the tube and use this solution to do a transient simulation of the exterior.


Can you coarsen your mesh somehow? It will increase the timestep. Can you share a picture of your mesh near the holes?



Lats thing, you might achieve a larger timestep using PIMPLE and several outer iterations to use a larger Courant number.
hogsonik likes this.
Alczem is offline   Reply With Quote

Reply

Tags
compressible flow, openfoam, pressure-driven flow, rhopimplefoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Cyclic boundary for inlet and outlet with constant pressure gradient to drive flow mldughi OpenFOAM 1 August 1, 2020 20:24
Pressure Outlet Targeted Mass Flow Rate LuckyTran FLUENT 1 November 23, 2016 10:40
Maintaining Static Pressure at Fluid Flow Inlet cdevalve FLUENT 3 January 14, 2012 00:11
flow simulation across a small fan jane luo Main CFD Forum 15 April 12, 2004 17:49
what the result is negatif pressure at inlet chong chee nan FLUENT 0 December 29, 2001 05:13


All times are GMT -4. The time now is 08:45.