|
[Sponsors] |
January 20, 2014, 12:53 |
Supersonic Flow Error
|
#1 |
New Member
Join Date: Oct 2013
Posts: 4
Rep Power: 12 |
Hey all,
I'm trying to model a Mach 2 flow over a missile. However I keep getting errors in the middle of my simulation saying "time step reduced in 3 cells due to excessive temperature change". I get a similar error for pressure as well. I attached a pic of my mesh and the error (along with residuals) for reference. I feel this problem is due to wrong boundary conditions and I'm new to supersonic flows in Fluent. Here are my settings: Density solver Transient Inviscid with energy equation on Density as ideal gas Constant specific heat Boundary condition: -On the far boundaries (inlet and far from the missile) chose a pressure farfield setting with Gauge=101325Pa, Mach 2, 300K and the flow direction -On the outlet chose pressure outlet with Gauge=0Pa normal to boundary -Operating pressure=0Pa Solution initialization: compute from inlet Solution Methods: Implicit, AUSM flux, Gradient=Green Gauss Cell Based, Flow=3rd order MUSCL CFL=2 Residuals at 1e-4 Time step=.1 Max Iterations/step=20 Let me know what you think. Thanks. fluent error.jpg mesh.PNG |
|
January 20, 2014, 13:11 |
|
#2 |
Super Moderator
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,399
Rep Power: 46 |
Compressible supersonic flows are not my strong point, but are you sure it is a good idea to set both the gauge pressure and the operating pressure to zero?
The time step size seems may be way to high and I would not start the first simulation with a third order scheme. It might aswell have to do with bad mesh quality, but unfortunately there is nothing to see in the picture you posted. Does the simulation really have to be transient? |
|
January 20, 2014, 13:25 |
|
#3 |
New Member
Join Date: Oct 2013
Posts: 4
Rep Power: 12 |
I read somewhere on these forums that supersonic flows needed 0Pa for the operating pressure. As far as the gauge pressure it's only 0 on the outlet, but 101325Pa in the farfield and inlet, but thinking about that twice I feel that they should be consistent (maybe I'll try both at 101325Pa).
I don't believe is is the mesh quality, at it has roughly 8.5 million cells (although I'm not sure about the orthogonality of the mesh). The mesh radius is roughly 20 diameters away. I chose transient due to that it is slightly easier for Fluent to solve (as opposed to steady where some damping of the solution occurs). I was going to run transient until a steady state is approached. I don't really have a reasoning for the time step size though, I'll try some lower. |
|
January 20, 2014, 13:32 |
|
#4 |
New Member
Join Date: Oct 2013
Posts: 4
Rep Power: 12 |
I should also add that I get this same "time step in 3 cells changes due to excessive temperature change, reducing the CFL to something lower" regardless if I use a steady OR transient solver.
|
|
January 20, 2014, 13:41 |
|
#5 | |
Super Moderator
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,399
Rep Power: 46 |
Quote:
Check the mesh quality explicitly. Operating and gauge pressure set to zero result in an absolute pressure of zero (as far as I know) which could explain one of the error messages you get. |
||
September 14, 2019, 07:52 |
|
#6 |
New Member
Join Date: Apr 2019
Posts: 1
Rep Power: 0 |
Hello, I recently met this problem, very like your first picture, and I also make a model about a Mach 2 flow in the combustion chamber. So, could you please give my some advises?
|
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Building OpenFOAM1.7.0 from source | ata | OpenFOAM Installation | 46 | March 6, 2022 13:21 |
Ansys Fluent 13.0 UDF compilation problem in Window XP (32 bit) | Yogini | Fluent UDF and Scheme Programming | 7 | October 3, 2012 07:24 |
POSDAT problem | piotka | STAR-CD | 4 | June 12, 2009 08:43 |
Problem with compile the setParabolicInlet | ivanyao | OpenFOAM Running, Solving & CFD | 6 | September 5, 2008 20:50 |
DecomposePar links against liblamso0 with OpenMPI | jens_klostermann | OpenFOAM Bugs | 11 | June 28, 2007 17:51 |