
[Sponsors] 
April 8, 2015, 10:34 
Differencing Scheme for Advection Problem

#1 
New Member
Simon Hahn
Join Date: Aug 2014
Location: Germany
Posts: 6
Rep Power: 11 
Hello everyone,
Again I'm doing my a thesis in numerical fluid mechanics and I have a  transient  advection problem that can be described by the following equation: d (rho*phi) \ dt + d (rho*phi*u) \ dx = q_0 (transient term) + advective/convective term = source term As you can see, the Peclet number, defined as quotient of convective/advective and diffusive transport is "quite" high, due to the lack of diffusive mass transport. By discretising it with FVM it leads to: transient term + [Flux_east * phi_east  Flux_west * phi_west) = source term So far I tried several schemes to discretise the property phi on the cell face(s):  1. order upwind: very stable, but massiv numerical diffusion which increases with increasing number of time/iteration steps; starting with a sharp edge at time and nodal point 1, the edge 'melts away' to a long curve  2. order upwind (LUDS): very stable, less numerical diffusion, but finally the diffusion increases massivly again with increasing number of time/iteration steps; same problem like 1. order upwind  QUICK: leads to oscillations, due to the high PecletNumbers  QUICK in combination with Hayase et. al. and deferred correction: stability of Upwind in the coefficient matrix 'yes' but it does not work: the deferred correction does not converge, consequently the 'upwindsolution' is not that much refined  analytical solution: not applicable because of the transient source term q_0 and transient boundary conditions  other central differencing schemes: not applicable due to form of the governing equation During my internet research I found several hints like using "AUSM" (Advection Upstream Splitting Method) or "Newton Raphson" as direct linearisation of phi. But so far I could not find a really good comprehensible paper/journal article/thesis that clearly describes how to apply it on the property phi_east/west in the discretised equation. I would be very thankful for any help and open for questions if anyone has the same problem or about the listed schemes above! Best regards, Simon! 

April 8, 2015, 11:43 

#2  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,675
Rep Power: 71 
Quote:
1) in your case Peclet is infinity 2) depending on the solution, you must use the integral form for nonregular solution (step solution) 3) It is known that linear schemes greater to first order accuracy are not monotone. In such cases you have to use some care (for example flux limiter) to make the scheme nonlinear even for the linear equation 4) You can try the QUICKEST scheme, very accurate for regular solution. It can be supplied by a limiter for non regular solution 

April 8, 2015, 11:59 

#3 
Member
robo
Join Date: May 2013
Posts: 47
Rep Power: 12 
Laney's text on computational gas dynamics covers this class of problems pretty comprehensively; and Pletcher, Anderson & Tannehill also cover gas dynamic problems, including the use of Ausm and similiar schemes. You're looking for FV methods for the Euler equation, in my experience it is not typical to discuss such problems in terms of a Peclet number. Try searching for Riemann solvers, flux limiters, TVD schemes etc.


Tags 
advection, ausm, newtonraphson 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Problem using Central differencing scheme with LES in CFX  selvam2487  CFX  8  April 11, 2016 14:15 
About the central differencing scheme  James  Main CFD Forum  6  July 21, 2009 18:35 
Differencing Scheme MARS  michael  OpenFOAM Running, Solving & CFD  1  August 10, 2006 11:51 
problem about numerical scheme in LES.  libin  Main CFD Forum  4  July 1, 2004 05:32 
Central Differencing Scheme in Fluent  alice  FLUENT  2  March 21, 2004 08:30 