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

Lax Wendroff Euler 2D in physical domain

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 26, 2012, 09:56
Default Lax Wendroff Euler 2D in physical domain
  #1
New Member
 
Join Date: Dec 2012
Posts: 4
Rep Power: 4
Christian87 is on a distinguished road
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 Lax-Wendroff 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 Lax-Wendroff 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
Attached Images
File Type: png control_volume2.png (18.4 KB, 3 views)
Christian87 is offline   Reply With Quote

Old   December 26, 2012, 12:22
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,651
Rep Power: 23
FMDenaro will become famous soon enough
You must use the Lax-Wendroff 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
FMDenaro is offline   Reply With Quote

Old   December 27, 2012, 02:47
Default
  #3
New Member
 
Join Date: Dec 2012
Posts: 4
Rep Power: 4
Christian87 is on a distinguished road
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...
Christian87 is offline   Reply With Quote

Old   December 27, 2012, 05:51
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,651
Rep Power: 23
FMDenaro will become famous soon enough
Quote:
Originally Posted by Christian87 View Post
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...

you can use the FV flux reconstruction method showed in the book of Peric and Ferziger.
Furthermore, I worked on LW-like 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.
FMDenaro is offline   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
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


All times are GMT -4. The time now is 18:38.