
[Sponsors] 
October 31, 2006, 06:20 
WENO  Basic Question

#1 
Guest
Posts: n/a

I am studying the WENO scheme as presented in the D S Balsara and CW Shu paper: "Mononicity Preserving Essentially Nonoscillatory Schemes with Increasingly High Order Accuracy", JCP, 160, (2000).
Most of the ideas are fairly straight forward, and I have been able to successfully implement various order schemes in Matlab that work fine for single 1D hypebolic problems and also that solve the Euler Equations. My question is: How are the smoothness indicator coefficients calculated. I am obviously missing something and would be interested in a stepbystep explanation. Any help would be much appreciated. 

October 31, 2006, 09:51 
Re: WENO  Basic Question

#2 
Guest
Posts: n/a

I'm doing my MSc. thesis on WENO and dealing with it for a long period of time, but I haven't dealed with MPWENO. I skimmed this paper, this variaton of WENO uses the same procedure to do the reconstruction in general.
In order to find nonlinear weights assigned for the candidate stencils which are parts of your high order bigger stencil, linear (optimal) weights and the smoothness indicators for these candidates must be found beforehand. For uniform grids they found in a straightforward way. Optimal weights are some numeric constants, and the smoothness indicators are some combination of the values of the cells which are forming these candidate stencils. The formulations are given at the corresponding papers of Shu. So, you use the values of the cells which are a part of the candidates and find smoothness indicators for each of these candidates, then assign nonlinear weights for them with the help of these indicators and optimal weights, and finally use these candidates to form the higher order reconstruction for your bigger stencil. 

October 31, 2006, 10:21 
Re: WENO  Basic Question

#3 
Guest
Posts: n/a

Hi Saygin,
Thank you for your response. If you have got to grips with WENO, then implementation of MPWENO and ACM are fairly minor steps. Yes they do use approximately the same calculation for the Smootness Indicators (SI). However, my question is in relation to the actual stepbystep calculations for the basic WENO SI'. Have you managed to perform the actual calculations? If so, I would be grateful if you could outline the steps, give me a reference, or email me an attachment that details the calculations. Regards 

October 31, 2006, 14:35 
Re: WENO  Basic Question

#4 
Guest
Posts: n/a

Hi Graham,
The actual calculations of these reconstruction polynomials, linear weights and smoothness indicators are not so easy to do. I did some calculations in order to extend my uniform domain to a nonuniform one. I can say that if you won't use rectangular domain these calculations are useless. Because in order to handle curvilinear grids coordinate transformation has to be done, and this transformation projects your domain to a uniform one. About your question, in order to find those equations that you see in the papers, first of all you have to find Lagrange interpolation polynomials for the order you are using. For example, for 5th order WENO you have to expand there 3rd order polynomials. After that the smoothness indicators can be found by calculating the L2norm of all derivatives of these there 3rd order polynomials. Have a look at this paper, it will help you I think. In addition don't forget to use square of the derivatives, it is not printed on the paper. "Grid adaptation with WENO schemes for nonuniform grids to solve convectiondominated partial differential equations" J.Smit et al. Chemical Engineering Science 60 (2005) 26092619 

October 31, 2006, 18:52 
Re: WENO  Basic Question

#5 
Guest
Posts: n/a

Hi Saygin,
Thanks for the information. I have obtained a copy of J Smit's paper and it looks as if it may give me the information I need. Regards 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
ICEM CFX a basic question  Al Mazdeh  CFX  2  April 5, 2008 11:51 
Basic question: UDF for wall heat flux  Carl  FLUENT  1  August 5, 2006 19:01 
Basic mutli component question  Eduardo Mendoza  FLUENT  0  February 17, 2005 09:15 
LES..Basic Question  NIK  Main CFD Forum  13  May 19, 2004 19:03 
Basic Question on Convergence  Anshul  FLUENT  1  August 2, 2002 02:14 