Implicit steger warming fvs
 

Hiiii, i am writing code for the implicit steger warming flux vector splitting for 1d cd nozzle case. i am using gmres method to solve Ax = b.
One strange thing is happening with in the code, as all of us know that for implicit there is no restriction on CFL number, First point is that if i start with CFL value 0.1 and max CFL value 2, with incrementation in CFL given by CFL*1.1 my code is converging nicely upto log rss value -7.034528 and results are matching.
But if i am taking max CFL number more than 2 say 4 the its blow up and go for log rss value -1.234 so on or if further increase in max CFL say upto 10 i will get nan after some iteration. I am really getting fed up with it from last one month then one of my friend told that u can put your problem in cfd forum surely you will get nice ideas. All ideas are welcome and thanks in advance plz help me

gmres is not easy to get working for flows with shocks.

From your description I assume you are solving a steady state problem.

You do need to use a smaller cfl to begin with. You can slowly increase this as your residual decreases. So you can try something like



You can also put an upper bound on the cfl.

Note that without good preconditioner, gmres will not work.

Thanx parveen for reply, but for 1d cd nozzle case as i said earlier i am starting with a low value of CFl 0.1 and incrementing with CFL*1.01 .
Shall i have to change the iterative method ?

Hey only, one reply from this forum not expected :mad: I am really in very much trouble help me plzzzzzzzzzzz