
[Sponsors] 
December 13, 2006, 02:22 
Flux splitting for WENO scheme

#1 
Guest
Posts: n/a

Hi,
I have implemented the WENO scheme. For a basic study, I have implemented the Van Leer flux splitting, but, instead of using MUSCL, i have used WENO for projecting the variables on to the face. I am solving a compressible NS equation. For now, I have restricted to Euler equations till I get the schemes and code right. Can anyone suggest me a good flux splitting scheme to apply for WENO scheme (for compressible NS eqn) that is less dissipative. Thanks in advance for your help. 

December 13, 2006, 06:23 
Re: Flux splitting for WENO scheme

#2 
Guest
Posts: n/a

I am not well versed in WENO..but why dont u try Steger warming flux splitting method? i hope its less dissipative that van leer. i found it easy to ger euler solution using steger warming.. things would be easy if u work in the positive mach number domain.


December 13, 2006, 07:14 
Re: Flux splitting for WENO scheme

#3 
Guest
Posts: n/a

Hi DSS,
Do you have the repot "ICASE 9765"? If you don't, I strongly recommend you to have it before implementing WENO schemes. At the end of page 23 of that paper, it says: "We would need the positive and negative fluxes to have as many derivatives as the order of the scheme. This unfortunately rules out many popular fluxes (such as those van Leer and Osher) for high order methods in this framework." I can say that the default flux splitting for finite difference WENO is LaxFriedrichs flux splitting, I don't know any other alternative, but there should be some other of course. However, if you want to use finite volume decomposition you can find much more different numerical flux alternatives. Saygin 

December 13, 2006, 07:41 
Re: Flux splitting for WENO scheme

#4 
Guest
Posts: n/a

Hi Saygin,
I implemented the WENO scheme straight from ChiWang Shu's paper. But, since I had already done Vanleer prior, I just performed a wedding ceremony of both But, as you say, I should have a look at the report. Also, recently, when I referred to a paper by CrnjaricZic titled "On different flux splittings and flux functions in WENO schemes for balance laws", I noticed that they use Marquina's Solver and HLL solver. May be I should try HLL which is supposed to be better. As a second question, is WENO worth the computational cost? Because, when I am to do an LES of a compressible jet for mach no. around 0.7, I believe the computational cost becomes an important factor. Thanks for helping me out, DSS 

December 13, 2006, 16:05 
Re: Flux splitting for WENO scheme

#5 
Guest
Posts: n/a

Hi,
Having already implemented a FVS like Van Leer, you may have a try with AUSM+ and more recent derivatives (AUSMPW+, AUSM+up, etc.). Regards 

December 14, 2006, 08:16 
Re: Flux splitting for WENO scheme

#6 
Guest
Posts: n/a

Hi DSS,
First of all, for some other flux implementations look at the papers of Toro. I guess they improved HLL to HLLC. I used global LaxFriedrichs and LLF flux in my Msc. thesis. It is a dissipative one, however, the simplest of all, thus it gives computational advantage. Does it worth the computational cost? It depends on the problem you are trying to solve. High order schemes are powerful in areas like, complex flow structures, shock interaction with complex flow structures, long time simulations ...etc. In that ICASE paper, there is a comparison between TVD and WENO for a long time simulation on vortex evolution. The power of using a high order scheme can be seen clearly there. So, if you have to solve some problem like that, I think it would be better to sacrifice computing time. Also, you can use finite difference WENO, it needs less computation time compared to finite volume version. It directly uses the flux for reconstruction and no need to integrate over the cell boundary as it is finite difference. However, using finite volume will cost you a lot. Because, say you are using 5th order reconstruction, then you have to find a 5th order flux which is entering the cell. So, for example a Gaussian integration needs to be done on the cell face. In order to "reconstruct the reconstructed value" you have to sweep the direction, which corresponds to the cell faces, with WENO again and may be again. I use two point integration point in the thesis, and that means 6 WENO reconstructions (3 for each direction) for each time step for 2 dimensional (dimensionbydimension) case. In addition, you are using RungeKutta time stepping which has to advance one time step in 3 or 4 or if you like 5 substeps. And I can say that, forget using a desktop computer for finite volume with very fine grid in 2D. As I tried to explain, high order methods cost too much, but I think there is a tendency towards using them. CPU powers are increasing daybyday, and for some time in the future, not decades I think, the attraction of using low order methods will completely diminish. Saygin 

January 10, 2007, 11:59 
Re: Flux splitting for WENO scheme

#7 
Guest
Posts: n/a

Hi,
I have implemented the Finite Difference WENO scheme. But, with Finite Volume, I think I need to understand clearly how the reconstruction is done. Please correct me if I am wrong. Considering a 3D FV WENO reconstruction, 1. Obtain face average values at i+1/2 from cell averages with the standard WENO procedure (same as that of Finite Difference) 2. Obtain a reconstructed value for gaussian point along one auxillary dimension 3. For these values, again obtain a reconstructed value for gaussian points in another auxillary dimension 4. Integrate the face with Gausian Quadrature rule (a simple sum in case of two point Gaussian Quadrature) 5. Use the integrated value for computing the flux term So, in case of 2D, for a flux splitting scheme, with N gaussian points, we have One reconstruction along face normal N Reconstructions along auxillary dimension 1 N Reconstructions along auxillary dimension 2 One Gaussian Quadrature Integration Thus, there are (1+N*Dim) constructions required for a Flux term. As there are left and right terms, the total reconstructions would be 2*(1+N*Dim) Hope my understanding is right... Thanks in Advance, DSS 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
flux splitting  emrah  Main CFD Forum  3  August 30, 2014 08:53 
Flux splitting  Dr B.M. Smith (Smith)  OpenFOAM Running, Solving & CFD  19  January 9, 2013 05:07 
Flux Jacobian for flux splitting  mecheng1600  Main CFD Forum  0  December 5, 2009 21:07 
About the Flux splitting method  maximus  Main CFD Forum  0  February 24, 2006 01:05 
Flux vector splitting of Inflow/outflow BC  Quain Tchew  Main CFD Forum  4  March 25, 2002 20:43 