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

Pressure wave reflection compressibleInterIsoFoam

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

LinkBack Thread Tools Search this Thread Display Modes
Old   April 12, 2022, 05:47
Default Pressure wave reflection compressibleInterIsoFoam
New Member
Dennis Thuy
Join Date: Apr 2022
Posts: 3
Rep Power: 2
dplthuy is on a distinguished road
Hi all,

I want to model the breakup of a liquid column under the influence of a step change in pressure & velocity (shock wave), using the compressibleInterIsoFoam (VOF) solver. The goal of the case is to determine a minimum grid cell size to properly resolve the breakup behavior, by comparison to literature data (

The simulations are two-dimensional. I am not yet taking into account the presence of the liquid (i.e. alpha.liquid = 0 in the whole domain), just considering the propagation of the pressure wave through the domain. The domain is initialised with p=p_rgh=8.9e4 Pa (gravity is neglected), U = (0 0 0).
On the left side of the domain, I create an area with higher pressure, p=p_rgh=1.09e5 Pa, U = (80 0 0), as can be seen in attached picture, for t=0. Left domain boundary is an inlet with these values specified with fixedValue BC. Top and bottom boundaries are zeroGradient for all parameters.
After the pressure wave passes through the domain, I would like it to exit at the right outlet boundary. However, as shown in the picture, the pressure wave reflects back into the domain.

Initially, I used zeroGradient conditions for p, p_rgh, and U at the outlet. Since this resulted in wave reflection, I have also tried the advective and waveTransmissive for p, p_rgh and U. I have also applied pressureInletOutletVelocity or zeroGradient for U in combination with advective or waveTransmissive for p and p_rgh. None of these different conditions made any difference in the wave reflection.

I know that density-based solvers (such as rhoCentralFoam), are better suited for pressure wave modeling as compared to pressure-based solvers such as I am using. Indeed, when using rhoCentralFoam in combination with zeroGradient outlet condition for p and U for my case-setup, there is no wave reflection at the outlet. However, since I want to use this case to determine a proper VOF resolution for droplet breakup, I do not have the option of switching to a density-based solver.

Is it possible at all to avoid these pressure wave reflections when using a pressure-based solver? If so, which combination of boundary conditions is recommendable?

I hope I have provided enough information on my case, if not, please let me know what is missing. Thanks in advance!
Attached Images
File Type: jpg pressureFieldShockWave.jpg (34.3 KB, 13 views)
dplthuy is offline   Reply With Quote

Old   July 28, 2022, 11:53
Julio Pieri
Join Date: Sep 2017
Posts: 75
Rep Power: 7
JulioPieri is on a distinguished road
Try using something different for U. waveTransmissive can be paired with zeroGrad in U, for instance. I'm not sure you can use waveTransmissive for p_rgh: I'm facing this very problem right now.
JulioPieri is offline   Reply With Quote


boundary condition, pressure, shock waves

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
Generate convective pressure fluctuation Bananenflanke CFX 10 May 12, 2021 19:33
How to Install Numerical Wave and Pressure Gauges in the Numerical Wave Tank? mushtak Main CFD Forum 0 November 2, 2020 03:34
CFX Solver stopped with error when requested for backup during solver running Mfaizan CFX 40 May 13, 2016 07:50
Pressure wave pattern over body in steadystate? Chebeba CFX 1 March 16, 2008 03:00
Neumann pressure BC and velocity field Antech Main CFD Forum 0 April 25, 2006 03:15

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