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

Differencing Scheme for Advection Problem

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 8, 2015, 09:34
Default Differencing Scheme for Advection Problem
  #1
New Member
 
Simon Hahn
Join Date: Aug 2014
Location: Germany
Posts: 6
Rep Power: 12
Simon Smokes is on a distinguished road
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 Peclet-Numbers

- 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 'upwind-solution' 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!
Simon Smokes is offline   Reply With Quote

Old   April 8, 2015, 10:43
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,839
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Simon Smokes View Post
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 Peclet-Numbers

- 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 'upwind-solution' 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!

1) in your case Peclet is infinity
2) depending on the solution, you must use the integral form for non-regular 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 non-linear 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
FMDenaro is offline   Reply With Quote

Old   April 8, 2015, 10:59
Default
  #3
Member
 
robo
Join Date: May 2013
Posts: 47
Rep Power: 13
robo is on a distinguished road
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.
robo is offline   Reply With Quote

Reply

Tags
advection, ausm, newton-raphson

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
Problem using Central differencing scheme with LES in CFX selvam2487 CFX 8 April 11, 2016 13:15
About the central differencing scheme James Main CFD Forum 6 July 21, 2009 17:35
Differencing Scheme MARS michael OpenFOAM Running, Solving & CFD 1 August 10, 2006 10:51
problem about numerical scheme in LES. libin Main CFD Forum 4 July 1, 2004 04:32
Central Differencing Scheme in Fluent alice FLUENT 2 March 21, 2004 07:30


All times are GMT -4. The time now is 02:01.