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

negative pressure: Wall boundary, Mach 3, Euler

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 2, 2007, 03:23
Default negative pressure: Wall boundary, Mach 3, Euler
  #1
DSS
Guest
 
Posts: n/a
Hi,

I tried to simulate a test case with an initial start of Mach 3 flow. One of the case is the forward step, and another is the Ron-Ho-Ni channel.

When I run my code, after a few iterations, it crashes because of negative pressure value (in the vicinity of wall). In some literatures, it was mentioned that I need to do some entropy correction at the corners/singularity at wall. But, I am just wondering if there is an alternate way to do that. My code is an unstructed WENO type code, with NRBCs as proposed by Lele. For now, I use HLL flux for interfaces.

I tried to lower down the CFL to 0.01... but, still not getting that right...

Can anyone let me know where could I have possibly erred, or missed out? Should I do the entropy correction? and if I should do that, how can I generalize the boundary condition? (Say a Mach 5 problem for Scramjet intake)

(My code works well till Mach 2)

Thanks in advance for your help.
  Reply With Quote

Old   April 2, 2007, 16:14
Default Re: negative pressure: Wall boundary, Mach 3, Eule
  #2
saygin
Guest
 
Posts: n/a
Singularity at the corner of the forward facing step will cause some problems at the shock reflection area on the bottom boundary after the step. You have to cure that singularity in order to have a correct shock reflection. However, I'm not sure if that singularity is the main reason for negative pressure or not. But I'm sure that those negative pressure effects should be cured by local order reduction. That negative pressure arises because you are using a wide stencil for high order accuracy, which can cause problems when two shocks presents together in one stencil. Or may arise while you are trying to reconstruct that high order polynomial where differences of values are very high across a discontinuity (shock) while one part near a zero value. The high order polynomial can drop below zero at those circumstances. Try to use a local order reduction at that point. But as far as I remember from my work for Mach-3 ffs with WENO-5th order with LF flux, there was no need for local order reduction. So there may be another problem with the implementation.

Saygin
  Reply With Quote

Old   April 3, 2007, 05:46
Default Re: negative pressure: Wall boundary, Mach 3, Eule
  #3
DSS
Guest
 
Posts: n/a
Hi Saygin,

Can you let me know whether you did interpolation on the primitive variables or conserved variables or other means? I think I messed up in these...

Thanks for your help!

DSS
  Reply With Quote

Old   April 3, 2007, 22:45
Default Re: negative pressure: Wall boundary, Mach 3, Eule
  #4
DSS
Guest
 
Posts: n/a
Hi Saygin,

I figured out the error... I made a goof up when I transited my code from interpolating the conserved variables to interpolating the primitive variables. I missed out a division by density in the energy term. Now, the boundary condition works fine. But, am still curious... I read in literature that one should use charateristic variables in interpolation for improved order of accuracy. Did you try those?

DSS
  Reply With Quote

Old   April 4, 2007, 10:06
Default Re: negative pressure: Wall boundary, Mach 3, Eule
  #5
saygin
Guest
 
Posts: n/a
Hi DSS,

I used conservative variables for reconstruction/interpolation. Do you mean using characteristic space by the "characteristic variables"? If you mean that, yes it is correct. Component-wise decompositions are cheaper to implement and to run but gives oscillations and it is not prefered and I don't prefer also. Using characteristic-wise decomposition is costly compared to component-wise but gives non-oscillatory and elegant results. So in summary what I did was; after a Roe-averaging (for linearization -freezing the variables-) projecting the conservative variables to the characteristic space (by that way the vectors are decoupled), applying reconstructions there, and then projecting back to the physical space.

Saygin
  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
External Radiation Boundary Condition (Two sided wall), Grid Interface CFD XUE FLUENT 0 July 8, 2010 06:49
CFX doesn't continue calculation... mactech001 CFX 6 November 15, 2009 22:25
mass flow in is not equal to mass flow out saii CFX 2 September 18, 2009 08:07
Boundary condition in SIPHONS - negative pressure - FLUENT 0 April 3, 2008 04:51
Trimmed cell and embedded refinement mesh conversion issues michele OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 2 July 15, 2005 04:15


All times are GMT -4. The time now is 14:46.