
[Sponsors] 
December 26, 2012, 09:56 
Lax Wendroff Euler 2D in physical domain

#1 
New Member
Join Date: Dec 2012
Posts: 4
Rep Power: 5 
Hello,
I discovered this form recently and have some question about the 2D FVM. Forgive me if this thread already exists but I don't know the correct keyword to look for. I'm currently writing a solver for the 2D Euler equations in curvilinear coordinates and the later aim is to compute the flow around a turbine blade. The whole solver is supposed to work with Cartesian coordinates. I want to use the LaxWendroff method and here is the catch since I want to write it directly in the physical domain. Unfortunately, I've got some problems with the flux computation and I'm some kind of confused how to formulate it properly. Starting from a structured vertex centered grid (red points) I introduce a control (thick black line) volume by averaging with the neighbor points. The edges of this control volume are prescribed by the green points which come from the averaging of its surrounding 4 red neighbors. By applying the divergence theorem I can say that the local temporal change of my conserved variable is equal to the sum of all fluxes through the CV's surface or as equation: u_t =  (F_x + G_y) => u_t*dV =  sum[f*S*nx +g*S*ny] , where S is the cell surface and nx and ny the normal vectors. Ok no problem so far but now comes the problem keeping me busy for almost one week: Following LaxWendroff I add artificial viscosity in form of u_tt and express the second derivative in time by saying u_tt = [ (F_x + G_y) ]_t =[ (A*u_x + B*u_y) ]_t = [ ([A*u_t]_x + [B*u_t]_y ] Here I stuck because I haven't figured out yet how to do that for  sum[f*S*nx +g*S*ny] or vice versa how to formulate u_y and u_x for this CV? I know my question is more a mathematical one but if you could recommend me any papers or have some good hints, I would really appreciate this. Best regards, Christian 

December 26, 2012, 12:22 

#2 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,198
Rep Power: 27 
You must use the LaxWendroff method for conservative laws. Give a look to the book Morton & Meyers, Numerical Solution of Partial Differential Equations, subsection 4.6 is what you are looking for.
Also the book of Leveque can be useful 

December 27, 2012, 02:47 

#3 
New Member
Join Date: Dec 2012
Posts: 4
Rep Power: 5 
Thank you for your post. Yes, I just looked it up but unfortunately it is only for 1D.
I don't know but might it be wiser to switch to metrics here? So I can apply all the scheme without worrying about the orientation... 

December 27, 2012, 05:51 

#4  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,198
Rep Power: 27 
Quote:
you can use the FV flux reconstruction method showed in the book of Peric and Ferziger. Furthermore, I worked on LWlike method on unstructured grids while using the integral formulation, you can see details in Int. J. Numer. Meth. Fluids 2002; 38:1133–1157 (DOI: 10.1002/d.278). Specifically, see how the (16) to (18) are obtained and how the operators act on the fluxes Last edited by FMDenaro; December 27, 2012 at 08:34. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Domain Imbalance  HMR  CFX  3  March 6, 2011 21:10 
Difference Between the Oscillations of Lax Wendroff Method and that of Beam Warming  rizwanherl  OpenFOAM Programming & Development  0  June 8, 2010 16:08 
RPM in Wind Turbine  Pankaj  CFX  9  November 23, 2009 05:05 
air bubble is disappear increasing time using vof  xujjun  CFX  9  June 9, 2009 07:59 
Harmonic Balance (Frequency Domain) Euler CFD  Manish  Main CFD Forum  11  December 1, 2005 19:20 