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

Air blast with spherical mesh - standing low pressure wave

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

LinkBack Thread Tools Search this Thread Display Modes
Old   October 23, 2018, 16:01
Default Air blast with spherical mesh - standing low pressure wave
New Member
Andy S.
Join Date: Jun 2018
Posts: 17
Rep Power: 7
TheMadHungarian is on a distinguished road
Hi all,

Here's an issue my colleagues and I are stuggling for the last few weeks: we are simulating the progress of the blast shockwave in air at certain altitude.

The mesh is set up as a 1/8th of a complete sphere and each flat "side" is a symmetryPlane. The high p/T/U originates from the inside smaller spherical surface (called "blast") while the outer larger surface of the sphere is atmosphere (a waveTransmissive patch).
The boundary conditions are prescribed per the Friedlander equation (although negative T and p is meaningless as both T and p are absolute, so they donot dip below ambient values). In essence it means the T and p on the "blast" surface are described as time-dependent conditions: type uniformFixedValue; uniformValue table (.....

What we see during simulation (ParaFoam - Filters - Data Analysis - Plot over Line) is that the pressure wave is moving outward and decreasing...this is what you'd expect. But we also see a standing wave of low pressure located close to the inner sphere. That low pressure does not seem right since it just "sits" there, not moving, not decreasing even with 0.5 seconds simulation time. I dont think it makes sense at all since the pressure would equalize after a while.

This is a quick moving event, the blast wave comes out at about Mach 3, temps are 5000K at t=0 and decreasing as well as the blast pressure is 450 psi (3.1e6 Pascal) then decreasing.

We use rhoCentralFoam with laminar flow. The medium is air.

One problem is that since the velocity U at t=0 is defined as surfaceNormalFixedValue - this does not allow to describe the U over time although we could do that with T and p.

I have the case tar.gz attached...any input or insight would be useful. What's missing? What's not set up properly?
Are we not using the right kind of BC for T, p and/or U ?

Attached Files
File Type: gz blast_test.tar.gz (2.1 KB, 11 views)
TheMadHungarian is offline   Reply With Quote

Old   May 1, 2020, 01:23
Default blast wave simulation
Join Date: Feb 2015
Location: New York
Posts: 73
Rep Power: 11
opedrofunk is on a distinguished road
Hi Andy,

Came across this post. Interesting idea. There has been substantial progress recently extending OpenFOAM (e.g. the blastFoam solver) for high-explosive and air blast simulation ( If you and your group are working in this area, it might be of interest.

Best regards,
opedrofunk is offline   Reply With Quote


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
Pressure for wave height in an Water Wave Tank Cluain CFX 8 December 6, 2021 03:58
Pressure wave reflection at wall boundaries gotang CFX 3 July 18, 2018 19:52
air bubble is disappear increasing time using vof xujjun CFX 9 June 9, 2009 07:59
How to apply negtive pressure to outlet bioman66 CFX 5 June 3, 2006 01:40
Gas pressure question Dan Moskal Main CFD Forum 0 October 24, 2002 22:02

All times are GMT -4. The time now is 22:09.