CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

CFL-Free numerical methods

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 18, 2008, 05:01
Default CFL-Free numerical methods
  #1
joy
Guest
 
Posts: n/a
I want to use CFL-Free numerical methods. I did google it and come across interesting staff but still I could not understand it.

Can you explain how can we use it, is it for all type of PDEs. what is its relation to centeral methods....

  Reply With Quote

Old   August 18, 2008, 09:35
Default Re: CFL-Free numerical methods
  #2
Jed
Guest
 
Posts: n/a
Central versus upwinded is a stability issue separate from CFL. To avoid CFL, you need to use an implicit method in which case you can choose arbitrarily long time steps. Of course, if the time steps are much longer than CFL, you may not resolve temporal dynamics of the system.
  Reply With Quote

Old   August 18, 2008, 10:15
Default Re: CFL-Free numerical methods
  #3
joy
Guest
 
Posts: n/a
any good start then ?
  Reply With Quote

Old   August 18, 2008, 13:09
Default Re: CFL-Free numerical methods
  #4
Jed
Guest
 
Posts: n/a
Suppose you have an advection-diffusion equation

du/dt + a grad(u) - e div(grad(u)) = 0 in Omega u = 0 on Gamma

where a is the velocity field, e is viscosity. More generally this is written

du/dt + F(u) = 0

Discretize in time using any unconditionally stable ODE method. With backward Euler we get

G(u^{n+1}) = u^{n+1} - u^n + dt F(u^{n+1}) = 0

This is a linear (or nonlinear) system for u^{n+1}. In the general case, it will require a Newton (or Picard) iteration to solve. Let

J(u^n) = (dG/du)(u^n)

be the Jacobian (discretized Frechet derivative) and solve solve linear systems

J(u^n + v_k) v_{k+1} = -G(u^n + v_k)

starting from v_0=0 until convergence. Then

u^{n+1} = u^n + v_k

and move on to the next time step. Note that you only need one Newton step when the problem is linear, hence

u^{n+1} = u^n + v_1.
  Reply With Quote

Old   August 19, 2008, 05:39
Default Re: CFL-Free numerical methods
  #5
joy
Guest
 
Posts: n/a
Many thanks for your response and time, very much appreciated.

Am doing it now and I will see if I can get done with fortran. Just need to know then about the nature of the PDEs i.e. are we really need to confirm hyperbolicity of the equations and the Jacobian or just consider any system and solve it e.g. Euler equations.

What do you think then?

  Reply With Quote

Old   August 19, 2008, 06:31
Default Re: CFL-Free numerical methods
  #6
Jed
Guest
 
Posts: n/a
In advection dominated systems, you will need some form of upwinding to stabilize the advection term. You need to know what boundary conditions are necessary and sufficient (this is related to hyperbolicity which is a concept which is only really meaningful in the linear case).

If you are solving constrained equations (incompressible) using an iterative method, it will require extra work due to it being a saddle point problem. For incompressible Euler equations and modest to high Reynolds number Navier-Stokes, projection methods are appropriate for enforcing incompressibility.

Good luck.
  Reply With Quote

Old   August 19, 2008, 09:53
Default Re: CFL-Free numerical methods
  #7
vadim
Guest
 
Posts: n/a
CFL is related to convection/advection, so CFL free method is which has not time step limitation due to stability.

CFL free methods:

implicit treatment of convection can be stable (but handling nonlinearity is bit difficult and make method inefficient)

beter methods are semi lagrangian methods, e.g. CIP

  Reply With Quote

Old   August 20, 2008, 06:40
Default Re: CFL-Free numerical methods
  #8
student
Guest
 
Posts: n/a
what is CIP ?
  Reply With Quote

Old   August 20, 2008, 07:03
Default Re: CFL-Free numerical methods
  #9
vadim
Guest
 
Posts: n/a
check Yabe or Xiao publication in JCP
  Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Numerical schemes for free surface flows (VOF) botp OpenFOAM 2 March 11, 2011 16:27
numerical jacobian, matrix free methods Nikhil Main CFD Forum 1 December 5, 2008 04:51
How can I find a free Fortran77 numerical library? Alberto Main CFD Forum 3 January 21, 2008 01:56
1D numerical methods compared Andrei Main CFD Forum 0 October 14, 2002 04:22
help:spectral methods & divergence free functionsn D. Puigjaner Main CFD Forum 1 August 28, 2000 11:06


All times are GMT -4. The time now is 20:39.