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

WENO - Basic Question

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 31, 2006, 06:20
Default WENO - Basic Question
  #1
Graham
Guest
 
Posts: n/a
I am studying the WENO scheme as presented in the D S Balsara and C-W Shu paper: "Mononicity Preserving Essentially Non-oscillatory 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 step-by-step explanation.

Any help would be much appreciated.
  Reply With Quote

Old   October 31, 2006, 09:51
Default Re: WENO - Basic Question
  #2
saygin
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.

  Reply With Quote

Old   October 31, 2006, 10:21
Default Re: WENO - Basic Question
  #3
Graham
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 step-by-step 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
  Reply With Quote

Old   October 31, 2006, 14:35
Default Re: WENO - Basic Question
  #4
saygin
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 non-uniform 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 L2-norm 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 non-uniform grids to solve convection-dominated partial differential equations" J.Smit et al. Chemical Engineering Science 60 (2005) 2609-2619
  Reply With Quote

Old   October 31, 2006, 18:52
Default Re: WENO - Basic Question
  #5
Graham
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
  Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 21:34.