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

Differencing schemes for equation with no diffusion

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By FMDenaro

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 11, 2018, 03:06
Default Differencing schemes for equation with no diffusion
  #1
Member
 
Join Date: Jan 2018
Posts: 31
Rep Power: 8
lorbekl is on a distinguished road
Ok, imagine for example a 1D flow for which the momentum equation (or in general any sort of convection only equation) has only the convective term whereas the diffusive term is absent. The Peclet number is then undefined and since the problem is uni-dimensional no numerical diffusion is present.



Am I right in assuming that in such a case differencing schemes like the hybrid, exponential and power law scheme don't really have an advantage over a simple upwind scheme?



I suppose higher order differencing schemes like QUICK would however still apply?
lorbekl is offline   Reply With Quote

Old   October 11, 2018, 08:04
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,882
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by lorbekl View Post
Ok, imagine for example a 1D flow for which the momentum equation (or in general any sort of convection only equation) has only the convective term whereas the diffusive term is absent. The Peclet number is then undefined and since the problem is uni-dimensional no numerical diffusion is present.



Am I right in assuming that in such a case differencing schemes like the hybrid, exponential and power law scheme don't really have an advantage over a simple upwind scheme?



I suppose higher order differencing schemes like QUICK would however still apply?



The best and simple model you can start working on is the linear advection equation


df/dt + u*df/dx=0


You have the exact solution f(x,t)=f(x-u*t,0). When you introduce a discretization, you numerical solution can be seen as an exact (in analytical meaning) of the "modified equation", that is the original PDE plus the terms of the local truncation error. Depending on the chosen discretization, the modified equation can show "numerical diffusion", "numerical dispersion", "violation of monotonicity".

The QUICK discretization is not suitable for unsteady solution, you should consider the QUICKEST scheme. However, they do not ensure monotonicity of the solution and numerical oscillations appears, especially for solution with high gradients.
On the other hand the FTUS is monotone but poorly accurate, producing a lot of numerical diffusion. Using the FTCS still not improve the solution of the problem, the scheme has no stability.
lorbekl likes this.
FMDenaro is online now   Reply With Quote

Old   October 12, 2018, 11:12
Default
  #3
Senior Member
 
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,753
Rep Power: 66
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
You still have numerical diffusion. Is it important that Peclet number is or is not defined?

Quote:
Originally Posted by lorbekl View Post
Am I right in assuming that in such a case differencing schemes like the hybrid, exponential and power law scheme don't really have an advantage over a simple upwind scheme?

These schemes are only applied anyway only to the advective term (i.e. central differencing is used on the diffusive term). What you have accomplished by solving a pure advection equation is you have removed a different term (the diffusion term), which is discretized a different way. Ignoring non-linearities, you haven't changed anything about the way the advection terms are discretized and solved.
LuckyTran is offline   Reply With Quote

Old   October 12, 2018, 15:29
Default
  #4
Member
 
Join Date: Jan 2018
Posts: 31
Rep Power: 8
lorbekl is on a distinguished road
Quote:
Is it important that Peclet number is or is not defined?


Well how I understood these differencing schemes is that for example hybrid scheme switches to upwind scheme when Pe>2 and power law switches to upwind at Pe>10. So if there is no diffusion the Pe number goes to infinity and all these schemes would automatically revert to the upwind scheme would they not?
lorbekl is offline   Reply With Quote

Reply


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
Overview Differencing Schemes Tytator OpenFOAM Pre-Processing 0 November 13, 2014 10:57
Harmonic mean Vs Central Differencing (diffusion term- UDS) bharat.cmeri FLUENT 6 October 15, 2014 00:27
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 07:20
DIFFUSION free schemes question Main CFD Forum 10 June 1, 2004 06:04
differencing schemes for 3-D Convection-diffusion problems Nuray Kayakol Main CFD Forum 20 September 16, 1999 05:16


All times are GMT -4. The time now is 07:19.