mavguy |
February 15, 2015 22:03 |
Steady-State solution by time-stepping?
This might be a silly question, but I wrote a CFD code (compressible, unstructured mixed element in 2D and 3D) that I am having trouble converging in 3D for relatively low Mach numbers (~0.1).
I am aware that the N-S and Euler equations become stiff at low Mach numbers. How are these cases typically dealt with? Are solutions to complex 3D problems typically time stepped to a steady-state solution? This will certainly make the linear system easier to solve, but seems rather inefficient. Is seems that straight-up Newton-Krylov would be better.
I'm using GMRES to solve the linear system and have put in a lot of work to get a good preconditioner, to the point that I don't think I can do any better on the preconditioning side of things.
Any advice would be great!
|