
[Sponsors] 
when to use totally explicit schemes or dual timestepping in transient flow? 

LinkBack  Thread Tools  Search this Thread  Display Modes 
March 5, 2016, 06:40 
when to use totally explicit schemes or dual timestepping in transient flow?

#1 
New Member
dengli
Join Date: Mar 2015
Posts: 22
Rep Power: 10 
Hi,guys.In the simulation of a transient phenomenon, we can choose explicit time marching schemes or implicit time stepping(dual timestepping), but which schemes was better for the phenomenon? I saw some tips about this, for example,when the time scale are comparable to the spatial scales over the eigenvalue,i.e. when the CFL number dictated by the physics is of the order of unity. so what is the time scale? and the spatial scales is the length of computational domain or the grid length? and the eigenvalue of the problem is the u+c ??
thanks 

March 6, 2016, 12:50 

#2 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,399
Rep Power: 67 
I think that for transient simulation you should take care of the accuracy... implicit or explicit schemes is a secondary issue


March 8, 2016, 07:29 

#3 
Senior Member
Michael Prinkey
Join Date: Mar 2009
Location: Pittsburgh PA
Posts: 363
Rep Power: 24 
While I agree with Prof Denaro, I would mention that understanding all of the timescales in your problem is really critical. These do not need to necessarily involve complex computations either. Just perterm dimensional analysis will give you a lot of guidance. Length scales are the smallest cell length and the domain length. Velocities are typical and maximum. Timescales (for transient simulations) are the timestep and the final simulation time. If you juggle those in as representative numbers in your transport equations, you will see timescale emerge. Generally, you can balance the transient term against each term in the transport equation (turn off the other effects and look at the time scale for each phenomena).
Once you have those timescales, you can see which effects group together. In the incompressible limit, you will see "flow" timescales that are relatively large compared to acoustic timescales. That leads us to chose pressureimplicit schemes like Fractional Step/Projection methods that implcitly treats the pressure/mass conservation and removes that acoustic timescale from the stability analysis but (as Prof Denaro notes) NOT from accuracy considerations. Similarly, with reacting flows, the timescale of reactions can be very small relative to flow/acoustic timescales. That naturally leads to socalled stiff combustion solvers where computational cells are treated as independent reactors and time integrated with stiff ODE solvers to the next fluid timestep. IMO, understanding the various timescales is critical for simulations that involve transport more complicated than singlephase cold nonreacting flow. This is the only way to both make a solver that captures all important phenomena while not wasting lots of CPU time resolving irrelevant details. 

March 8, 2016, 09:02 

#4 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,399
Rep Power: 67 
from my experience, using explicit timemarching schemes fulfilling the numerical stability constraints are also a requirement for a good physical resolution of the relevant scales of the flow (apart form formulation where the scales are filtered out and then modelled).


March 8, 2016, 09:59 

#5  
New Member
dengli
Join Date: Mar 2015
Posts: 22
Rep Power: 10 
Quote:


March 8, 2016, 10:14 

#6  
Senior Member
Michael Prinkey
Join Date: Mar 2009
Location: Pittsburgh PA
Posts: 363
Rep Power: 24 
Quote:


March 8, 2016, 11:03 

#7  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,399
Rep Power: 67 
Quote:
Using the incompressible model, your governing equations will be not well suitable...Often the acoustic problem is decoupled from the flow solutions and obeys to a different set of equations solved after the incompressible solution is obtained. Alternatively, you can solve the full compressible problem if your Mach is not too low. 

March 9, 2016, 19:43 

#8  
New Member
dengli
Join Date: Mar 2015
Posts: 22
Rep Power: 10 
Quote:
So my opinion is ,in a transient flow simulation, we can only focus on the specific time scale we want to know, and ignore the other time scales?( fulfill the stability criterion) Thanks. 

Tags 
time integration, transient problem 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
High Courant Number @ icoFoam  Artex85  OpenFOAM Running, Solving & CFD  11  February 16, 2017 13:40 
mixerVesselAMI2D's mass is not balancing  sharonyue  OpenFOAM Running, Solving & CFD  6  June 10, 2013 09:34 
Low Mixing time Problem  Mavier  CFX  5  April 29, 2013 00:00 
Convergence criterion when using dual time stepping  abonfi  SU2  7  April 18, 2013 13:05 
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3  bookie56  OpenFOAM Installation  8  August 13, 2011 04:03 