
[Sponsors] 
August 16, 2000, 08:08 
spatial discretization

#1 
Guest
Posts: n/a

Many scientific journals require at least second order accurate discretization of the convective terms.
As I read it, from the methodology index, this is true for LUD, MARS, and naturally for QUICK. QUICK may result in solutions outside physical bounds, so I would prefer MARS (or LUD), but no reference is given to these 'nonclassical' schemes, which restricts the 'journalacceptability'. Any suggestion on what scheme to choose? and what arguments to give in an article to justify the choice? Thanks DAN 

August 16, 2000, 09:17 
Re: spatial discretization

#2 
Guest
Posts: n/a

(1). If you can obtain the mesh independent solution from your cfd calculations, then submit the results to a journal which does not require higherorder convection schemes. (2). For coarse mesh solutions, in some cases, higherorder schemes tend to give better solution relatively. But this is not always the case. The mesh independent solution is a must,because the solution is repeatable. (3). A nonrepeatable higherorder scheme solution is not a good solution. (4). My suggestion is: try two schemes and obtain mesh independent solutions. Then validate it to see which solution is better.


August 16, 2000, 09:44 
Re: spatial discretization

#3 
Guest
Posts: n/a

(1) I do not agree. There are good reasons for the journals to require higher order schemes. See e.g.
B. P. Leonard "Comments on the Policy Statement on Numerical Accuracy". J. of Fluids Engineering, vol 115, 1993 pp 339340. One reason is that the 1. order upwind actually solves a 'different' problem than was intended, even for a gridindependent solution. (2) The fact that in some (almost all) cases, higherorder schemes are superiour on coarse meshes is because the truncation error (for a given grid) is smaller for a higher order method. This also means that a gridindependent solution requires significantly more grid points for a 1. order method than for a higher order method! (3) Naturally, the solution must be gridindependent and consistent. (4) The 'best' solution for one setup is not necesserily the 'best' solution for other setups. I would rather know, beforehand, that I have used settings which (propably) ensure reasonable solutions. Regards DAN 

August 16, 2000, 11:43 
Re: spatial discretization

#4 
Guest
Posts: n/a

I suggest you to see all refereces of B. P. Leonard because he's interested a bout this subject for long time.
I suggest you to see this references, you can find many ideas. ** P. Tamamidis, and D. N. Assanis, "Evaluation of Various Highorderaccuracy Schemes with and without Flux Limiters", International Journal for Numerical Methods in Fluids, 16(10), 931948, 1993. ** M. K. Patel, and N. C. Markatos, "An Evaluation of Eight Discretization Schemes for Twodimensional ConvectionDiffusion Equations", International Journal for Numerical Methods in Fluids, 6(3), 129154, 1986. ** M. K. Patel, N. C. Markatos, and M. Cross, "A Critical Evaluation of Seven Discretization Schemes for ConvectionDiffusion Equations", International Journal for Numerical Methods in Fluids, 5(3), 225244, 1985. ** D. L. Roberts, and M. S. Selim, "Comparative Study of Six Explicit and Two Implicit Finite Difference Schemes for Solving Onedimensional Parabolic Partial Differential Equations", International Journal for Numerical Methods in Engineering, 20(5), 817844, 1984. ** A. Rigal, and G. Aleix, "Stability Analysis of Some Finite Difference Schemes for the NavierStokes Equations", International Journal for Numerical Methods in Engineering, 12(9), 13991405, 1978. ** John C. Strikwerda, "Highorderaccurate schemes for incompressible viscous flow", International Journal for Numerical Methods in Fluids, 24(7), 715734, 1997 ** Alexander G. Churbanov, Andrei N. Pavlov, and Peter N. Vabishchevich, "Operatorsplitting methods for the incompressible NavierStokes equations on nonstaggered grids. Part 1: Firstorder schemes", International Journal for Numerical Methods in Fluids, 21(8), 617640, 1995 **Carlos M. Lemos, "Higherorder schemes for free surface flows with arbitrary configurations", International Journal for Numerical Methods in Fluids, 23(6), 545566, 1996. 

August 16, 2000, 12:04 
Re: spatial discretization

#5 
Guest
Posts: n/a

(1). Item3, a good idea. But, just go back and review the papers published in the numerical schemes area, and count the number of papers which actually presented the mesh independent solutions. (2). I am using the upwind methods all the time,(including the commercial codes), it is similar to the velocity measurement, using Pitot probe, hot wire, or laser velocimeter. (3). If you are using a higherorder scheme, and the solution has wiggles in it for the particular problem you are trying to solve, then is it more accurate than the smooth upwind solution?


August 16, 2000, 12:46 
Re: spatial discretization

#6 
Guest
Posts: n/a

In general cases it's the upwind scheme gives good results but development in Taylor series of upwind scheme is only first order, so, if a higherorder scheme is used, and the solution for the specific problem has wiggles it's depended of how kind of scheme we use ; we can for exemple use a very simple and very recommended scheme like "Hybrid" order 1.
For the other schemes with and without Flux Limiters it's much difficult to talk a bout them. Try to see this book, you want to find many exemples in order to undrestand first how to make a comparaison between schemes An Introduction to Computational Fluid Dynamics: The Finite Volume Method By: Versteeg, H.K. and Malalasekera, W. AddisonWesley, 1996 

August 16, 2000, 13:41 
Re: spatial discretization

#7 
Guest
Posts: n/a

(1). You are talking about the schemes, while I am talking about the solutions. (2). Even with the upwind solution as the initial guess, it is still very difficult to get the solution to converge, when I switch to the higherorder scheme. (3). And if I use the higherorder scheme right from the begining, the solution does not converge unless I adjust the time steps or the relaxation factors many times in many steps.


August 16, 2000, 19:10 
Re: spatial discretization

#8 
Guest
Posts: n/a

One of the reasons why the first order scheme is easier to converge is that it is numerically dissipative, thus increasing the artifical (numerical) viscosity in the solution. The second order scheme tends to be dispersive, creating the oscillations that you see in the solution, especially in your pressure or Mach number plots, thus making it harder to converge. This has to do with the mathematics of the numerical solutions, and the previous discussion lists a good reference to explain this. I guess it depends on your problem to which scheme that you want to use...perhaps a hybrid scheme like someone listed before so that you don't lose too much accuracy but you introduce enough "stability" to converge your solution.


August 16, 2000, 20:07 
Re: spatial discretization

#9 
Guest
Posts: n/a

(1). Good suggestions. (2). Now, in order to get higherorder scheme solutions, we have to become an artist first. (If on the top of it, there is uncertainty in turbulence models, and lower order boundary conditions, the the whole thing is really a mess.)(3). I think, if one is interested in publishing papers in journals, it is a good idea to invent higherorder schemes.


August 17, 2000, 06:41 
Re: spatial discretization

#10 
Guest
Posts: n/a

In order to get higherorder scheme solutions, it's important to be an artist in Flux Limiters.


August 29, 2000, 10:46 
Re: spatial discretization

#11 
Guest
Posts: n/a

Dan, I've just started trying to use StarCD and I have some suggestions on getting references. (1)there is a reference for LUD: Wilkes, N.S. and Thompson, C.P. 1983 "An Evaluation of Higher Order Upwind Differencing for Elliptic Flow Problems" CSS 137, AERE, Harwell. (2)SFCD, Gamma differencing, and Blended Differencing are hybrid schemes discussed in Culbert Laney's book "Computational Gasdynamics". go to http://capella.colorado.edu/~laney.(3)MARS looks like one of the higher order reconstructionevolution (in contrast w/ limited) schemes discussed also in Laney's book. (3) another way too look for refs is to search the journals for articles written by others using StarCD including employees of CD which may have the appropriate references. often they'll write a paper as a marketing tool to showcase new improvements. the people at fluent do this all the time I don't know about Cd though. Check CD's website which has links to StarCD user conferences (under news) where you may be able to obtain papers etc. (4) the most unlikely way is to ask the people at CD. they may or may not tell you. I've had other companies send me papers that they've published and others have papers online. CD may have another policy in this area but it doesn't hurt to try.


August 30, 2000, 08:31 
Re: spatial discretization

#12 
Guest
Posts: n/a

Dan,
in addition to what I wrote before you can go to http://www.adapcoonline.com/frontpage.html in the "published works" section you'll see many articles written by StarCD users some of which look (by the journal to which they were submitted) to be useful. you may be able to find one that was submitted to the journal you're interested in. 

August 30, 2000, 08:52 
Re: spatial discretization

#13 
Guest
Posts: n/a

Clifford
Thanks for the reply. I have discussed the subject with the CD people and got the reply on the MARS scheme from Mike Lewis that: "We haven't published work on MARS because it is commercially sensitive. However, it is a TVD scheme, so you could use a reference to van leer and another to our methodology manual. (Other commercial codes have similar MUSCL type schemes)" The van Leer reference is: van Leer, B., Towards the ultimate conservative differencing scheme V: A secondorder sequel to Godunov's method, J. Comp. Phys., 23 (1977), pp101136 So, my solution is to write something like (where cite is a citation): The SIMPLE algorithm was used with a second order MUSCL type scheme with a TVD limiter (\shortcite{vanLeer1977}) for the convective terms (denoted MARS in \shortcite{starcd1999}). 

August 30, 2000, 09:29 
Re: spatial discretization

#14 
Guest
Posts: n/a

when i read up on the MARS scheme in the manual the description looked like some kind of ENO scheme so i guess the description is kinda misleading. As for "commercially sensitive" it's probably some warmed over Van Leer type limited scheme I guess.


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Spatial gradients of scalars  Jaishree  Siemens  4  July 2, 2007 11:02 
temporal and spatial discretization.  pity  Main CFD Forum  0  June 10, 2007 22:53 
Spatial and temporal discretization  Giuseppe  FLUENT  2  January 23, 2006 06:18 
spatial and transient B.C. using UDF  Park  FLUENT  0  April 28, 2005 21:32 
FIDAP spatial accuracy  Ingo Meisel  FLUENT  0  March 8, 2004 11:58 