|
[Sponsors] |
Stability of 2nd ordr one sided upwind 1st order in time |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
December 8, 2014, 19:06 |
Stability of 2nd ordr one sided upwind 1st order in time
|
#1 |
Member
Ardalan
Join Date: Feb 2010
Location: Baku
Posts: 50
Rep Power: 16 |
As mentioned in chapter 21 of Hirsch's book:
"second-order one-sided differences in space with a first-order difference in time (explicit), is unconditionally unstable" This can be proved for linear convection equation. I want to know that whether this is true for Euler equations or not. I get converged solutions for one dimensional nozzle flow with shock by applying the above mentioned method! However some overshoots happen at shock region. Also, when I implement limiters (minmod) to my scheme, solution convergence trend stops at about L2norm(p(n+1)-p(n))=5e-5. I can't understand what is happening. |
|
December 15, 2014, 20:29 |
|
#2 |
Member
Ardalan
Join Date: Feb 2010
Location: Baku
Posts: 50
Rep Power: 16 |
no answer
?? |
|
December 17, 2014, 18:24 |
|
#3 |
Senior Member
Join Date: Oct 2011
Posts: 242
Rep Power: 16 |
Hi Ardalan,
Maybe if you use sufficiently coarse grids, then the dissipation is sufficient to damp the enherent instability of the scheme. Try using much finer grids, I expect some troubles. I think the overshoots you observe are only the beginning of a blow up with a finer mesh. Concerning the limiter, I am not suprised you experience convergence problems to reach steady state. This is caused by the non-differentiability of minmod function at 1. Some remedy was proposed by Venkatakrishnan to improve convergence to steady state. They use a smooth limiter function, differentiable. If you look for it on the web you will find many articles showing nice convergence graphs using smooth limiters whereas non-diff limiters such like minmod stall Hope it helps |
|
December 17, 2014, 20:46 |
|
#4 |
Member
Ardalan
Join Date: Feb 2010
Location: Baku
Posts: 50
Rep Power: 16 |
Thanks for your useful reply.
Now I made my code second order both in time and space. Of course. convergence stall happens again. The fact that makes me surprised, is that by increasing cfl number, convergence stalls at lower residuals. For example, at cfl=0.1, convergence stops at about 1e-5 but for cfl=0.5, convergence stops at about 1e-8 and for cfl=0.9, convergence stops at about 1e-12. I'm using minmod limiter. Is this reasonable? |
|
December 18, 2014, 04:09 |
|
#5 |
Senior Member
Join Date: Oct 2011
Posts: 242
Rep Power: 16 |
Hello,
I do not know for sure about the behavior of residual as a function of CFL. My opinion is that using a lower CFL then dt is smaller and so the time truncation error. Then your scheme dissipates more with a higher truncation error and reaches convergence more easily. I really suggest you to try smooth limiters, if your scheme is well implemented it should be very easy to plug one or another limiter. On the wikipedia page you have several examples. I did not do the tests myself, but take for instance the ospre or van albada 1 limiters, it might help to reach convergence. Cheers, François |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
1st or 2nd order gradients approach? | fireflies | Main CFD Forum | 2 | September 18, 2014 08:07 |
AMI interDyMFoam for mixer nu problem | danny123 | OpenFOAM Programming & Development | 8 | September 6, 2013 02:34 |
mixerVesselAMI2D's mass is not balancing | sharonyue | OpenFOAM Running, Solving & CFD | 6 | June 10, 2013 09:34 |
Changing LimitedLinear to blend with 2nd order upwind instead of 1st order upwind | stevenvanharen | OpenFOAM Programming & Development | 0 | April 11, 2011 05:54 |
Multicomponent fluid | Andrea | CFX | 2 | October 11, 2004 05:12 |