CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Shock wave reflection

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

LinkBack Thread Tools Display Modes
Old   April 13, 2011, 14:35
Default Shock wave reflection
New Member
Join Date: Sep 2010
Location: Russia , Moscow.
Posts: 14
Rep Power: 9
Zato_Ichi is on a distinguished road
Greetings !

I tried to simulate a 1D shock wave generation in closed domain. Pressure field initially was set up with discontinuous conditions, from funkySetFields dictionary:

field p;
expression "3";
condition "(pos().x)>=(1)";
keepPatches 1;

So pressure equals 1 Pa on left and 3 Pa on right.

Boundary conditions for T and p are zeroGradient on both sizes, and U equal (0 0 0) on both sizes, so domain is closed by walls on right and left.
For calculating i used sonicFoam and rhoCentralFoam. And result of simulation was quite strange. Shock wave appears and propagates normally, but after reflection by the wall velocity U turns into zero, while pressure and temperature experience jumps, as it should be. As i know, in theory velocity must change sign and value (but not become a zero).

Maybe I'm putting wrong boundary conditions. I tried slip condition on velocity, but it didn't help.
Zato_Ichi is offline   Reply With Quote


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
shock tube reflection boundary conditions euler 1d micheletuttafesta Main CFD Forum 1 May 8, 2010 11:01
shock reflection Shuo Main CFD Forum 0 October 2, 2007 00:57
Question about the shock wave nucharin CFX 2 January 27, 2005 05:44
Question about the shock wave in CFX software nucharin Main CFD Forum 1 January 25, 2005 09:26
shock wave propagation Chetan FLUENT 1 April 20, 2004 14:45

All times are GMT -4. The time now is 04:20.