
[Sponsors] 
February 10, 2012, 21:42 
Confusing, but Fun results.

#1 
Member
Join Date: Mar 2009
Posts: 72
Rep Power: 10 
Hello,
Wife and kid sleeping since, well several hours ago (almost 3 AM here). A long week of work and what do I do? I spend my time writing CFD codes, for fun, of course (I use commercial codes only at work) I have written a very basic MACtype incompressible solver that I now benchmark against some different cases. The confusion that I have is the following: I tested to increase the timestep well above the stability limit (both convective and diffusive terms are handled explicitly so the limit should be of the order deltaX^2). Well the code did not blow up until I raised the timestep above convective term stability criteria (approx deltaX). The results were good compared to the following paper for Re 100 and 200, looking at the length of separation it was more or less identical. I might try to go higher in Re tomorrow but then I need to extend the mesh in the flow direction (probably, it seems like how things should work and how they actually work is very different) http://courses.washington.edu/me431/...alyjfm83.pdf I realize that I might have provided too little information to actually give any input or comments, but I try anyway. I can always add information when my head is a bit more clear So, any thoughts? Cheers and good night! On a side note: The code converges nicely with the lower, "correct" timestep as well, but it obviously takes much, much longer.
__________________
"Trying is the first step to failure."  Homer Simpson 

February 10, 2012, 22:41 

#2  
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 699
Rep Power: 19 
Quote:
I guess the reason is that error will build up slowly and it will blow up at some time. How long did you let run. Anyway in MAC type solvers the problem is not low Re but high Re, at high Re it is very sensitive to time stepping. 

February 11, 2012, 06:40 

#3 
Member
Join Date: Mar 2009
Posts: 72
Rep Power: 10 
Hello,
I tried to run it for approx 20000 timesteps and it did not blow up, but produced the correct answer (the same answer was also obtained in far fewer iterations). Is it true that the stability of a method depends on the Reynolds number? I would have guessed that either a method is stable or it is not. Thank you.
__________________
"Trying is the first step to failure."  Homer Simpson 

February 11, 2012, 15:01 

#4 
Member
Join Date: Jul 2011
Posts: 59
Rep Power: 8 
I think the Reynolds number will come into play. Not exactly sure about the NS equations but I know when solving the heat equation the limit is something like alpha*dt/dx^2 < C where C is like 1/2 or 1/4. Anyways the alpha term is related to the density, conductivity, and heat capacity. So I'm thinking the limit for the viscous term in the NS equations would be like dt/(Re*dx^2) < C. You could try a Reynolds number of 1 or 0.1 and I bet it would blow up with a time step that is stable at Re 100.


February 11, 2012, 17:00 

#5 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 3,349
Rep Power: 37 
well, your question needs to be completed by some information ...
 are you using a standard firstorderintime and secondorder central discretization on staggered grid?  what about your test configuration, is a steady or unsteady test?  are you using stretched grid ?  how do you estimate the stability constraint, independently for diffusion and convection part? I would just remark that the linear stability analysis (Von Neumann) on a 3d LINEAR diffusionconvection equation discretized by a FTCS provides a constraint on the time step that simoultaneously involves the Re number and the velocity intensity in terms of the mesh sizes. But...such analysis is not exactly valid for the nonlinear NS equations. MAC is able to filter some spurious modes, therefore the numerical instability can be slow... 

February 11, 2012, 18:56 

#6  
Member
Join Date: Mar 2009
Posts: 72
Rep Power: 10 
Quote:
Unsteady and solve until steadystate. Uniform mesh spacing. I use CFL condition as well as a condition where Reynolds number as well as dX^2 is included. Normally both conditions are evaluated and the lowest is chosen for the timestep. For the test mentioned above I just used the CFL condition. Could you define numerical instability being slow? Do you mean an instability that will appear if we let the solution run indefinitely? Cheers!
__________________
"Trying is the first step to failure."  Homer Simpson 

February 12, 2012, 05:06 

#7  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 3,349
Rep Power: 37 
Quote:
[/QUOTE]I use CFL condition as well as a condition where Reynolds number as well as dX^2 is included. Normally both conditions are evaluated and the lowest is chosen for the timestep. For the test mentioned above I just used the CFL condition. Could you define numerical instability being slow? Do you mean an instability that will appear if we let the solution run indefinitely? Cheers![/QUOTE] There is a stability region that is determined by a simultaneous combination of Re, cfl parameters. I can suggest a reading of my paper Int. J. Numer. Meth. Fluids 2006; 52:393–432, you will find some analysis of the parabolichyperbolic equation discretized with a semiimplicit scheme. But the issue is generally valid also for the FTCS. Numerical instability can be asymptotic if some eigenvalues is just a bit greater than 1. let me know 

February 12, 2012, 05:16 

#8 
Super Moderator

FordPerfect, I have seen similar behaviour as you report. For heat equation, I cannot use time step which is greater than the stability limit. But for convectiondiffusion (NavierStokes), it seems to be not very strict; even a time step based on convection seems to work. My guess is that due to convection, the instabilities/disturbances are convected away and out of the computational domain, and do not get a chance to grow much. This does not happen for heat equation where disturbances will very soon affect all parts of the domain. I would suspect that for convectiondiffusion, it will blow up after a sufficiently long time if your time step does not account for diffusion.


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
OpenFOAM  Validation of Results  Ahmed  OpenFOAM Running, Solving & CFD  9  June 22, 2011 18:59 
Velocity spots in openFoam results  Valle  OpenFOAM Running, Solving & CFD  4  August 19, 2009 05:53 
How to plot a function over a time period?  Cirion0000  CFX  4  July 18, 2009 12:48 
how to view the complete results  shokry  FLOW3D  5  February 3, 2009 14:56 
CFX cylinder or sphere benchmark results  Mel  CFX  1  August 8, 2005 18:47 