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

Euler Equations Boundary Conditions

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 5, 2013, 14:51
Default Euler Equations Boundary Conditions
  #1
New Member
 
Join Date: Dec 2012
Posts: 4
Rep Power: 13
Christian87 is on a distinguished road
Dear CFD Community,

I'm currently trying to fix a problem with my subsonic boundary condtions for the 2d Euler equations. Attached you can see a contour plot of my problem. The problem can be found in

" A. LERAT, C. CORRE. A residual-based compact scheme for the compressible Navier-Stokes equations. Journal of Computational Physics, 170, p. 642-675, 2001. "

I'm using a two step Lax Wendroff scheme with a third smoothing step and the solver is vertex centered=> the boundary match with the physical boundary. Moreover there are fictitious cells around my boundary. My problem is the high entropy production on the lower side after the bump, the too small peak Mach number (Lerat has 0.95) and the deviation at the inlet. The walls are modelled by using fictitious cells where the velocity of the first interior line is mirror-sided and thus enforcing zero mass flux.

For the supersonic case everything is alright but for this case I'm suspecting my boundary condtions for the deviation. The reference is set to freestream and e.g. u8=1 and rho8=1. Based on the freestream Mach number M8 I can get all freestream variabels.

I prescribe the incoming Riemann Invariant R+= 1 + 2/(M8*(gamma-1)) and the entropy s8=1 at the inlet and extrapolate R- from the interior. Based on that I get my actual speed of sound a and the normal velocity.

For the outlet I do it similarly where I extrapolate s and R+ and prescribe the pressure (here freestream pressure p8). This gives me the density by using the entropy relation and so the speed of sound.

I do all of this exactly on the boundary since my second lines (2:i_max-1 & 2:j_max-1) define my boundary and all sourrounding fictitious points are just updated by linear extrapolation except for the fictitious wall points.

What might I do wrongly here?

Thank you very much

P.S. In the plot upper and lower are changed.

Christian
Attached Images
File Type: jpg contour.jpg (18.7 KB, 14 views)
File Type: jpg plot.jpg (30.5 KB, 14 views)
Christian87 is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Implementation of boundary conditions for FVM Tom Main CFD Forum 7 August 26, 2014 06:58
Water subcooled boiling Attesz CFX 7 January 5, 2013 04:32
domain imbalance for enrgy equation happy CFX 14 September 6, 2012 02:54
[Commercial meshers] Trimmed cell and embedded refinement mesh conversion issues michele OpenFOAM Meshing & Mesh Conversion 2 July 15, 2005 05:15
Exit Boundary Conditions N-S equations lavan Main CFD Forum 2 July 2, 2005 04:10


All times are GMT -4. The time now is 08:40.