In this 2D problem, a left moving shock wave propagates first through the stiffened gas before it hits a stationary bubble, which gains speed and loss its circular shape. After the shock wave hits the bubble, circular pressure waves are created, while a pair of vortices appears on the tail of the bubble, increasing the pattern complexity of the interface. This example is a classic multiphase flow example and is used to test our two-phase multiphase flow solver.

A left moving shock wave propagates first through the stiffened gas before it hits a stationary bubble, which gains speed and loss its circular shape. After the shock wave hits the bubble, circular pressure waves are created, while a pair of vortices appears on the tail of the bubble, increasing the pattern complexity of the interface. This example is a classic multiphase flow example [1], here we use it to test our two-phase multiphase flow solver.

References:

[1] G. Allaire, S. Clerc, S. Kokh, A five-equation model for the simulation of interfaces between compressible fluids, J. Comput. Phys. 181(2) (2002) 577-616

*DUALCESE_MODEL

*DUALCESE_CONTROL_SOLVER

*DUALCESE_CONTROL_TIMESTEP

*DUALCESE_CONTROL_LIMITER

*DUALCESE_INCLUDE_MODEL

*DUALCESE_BOUNDARY_NON_REFLECTIVE_SEGMENT_SET

*DUALCESE_BOUNDARY_REFLECTIVE_SEGMENT_SET

*DUALCESE_INITIAL_TWO-PHASE

*DUALCESE_INITIAL_TWO-PHASE_SET

*DUALCESE_PART_MULTIPHASE

*DUALCESE_EOS_SET

*DUALCESE_EOS_VAN_DER_WAALS_GENERALIZED

*DUALCESE_D3PLOT

*CONTROL_TERMINATION

*DATABASE_BINARY_D3PLOT

*DATABASE_BINARY_D3DUMP

*KEYWORD

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

*DUALCESE_CONTROL_SOLVER

$ icese igeom iframe MixSelect

Euler 2D fixed two-phase

*DUALCESE_CONTROL_TIMESTEP

$ iddt cfl dtint

1 .9 1.e-7

*DUALCESE_CONTROL_LIMITER

$ idlmt alfa beta epsr

2 1.0 1.0 .00

$

*DUALCESE_INCLUDE_MODEL

fluid_m2h.k

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

$ Setup the boundary conditions for fluid

$

$ inlet & Outlet

$

*DUALCESE_BOUNDARY_NON_REFLECTIVE_SEGMENT_SET

$ ssid

1

2

$

$ Solid BCs (up & down)

$

*DUALCESE_BOUNDARY_REFLECTIVE_SEGMENT_SET

$*DUALCESE_BOUNDARY_SOLID_WALL_SEGMENT_SET

$ ssid

3

5

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

$ Setup the initial conditions for fluid

$

*DUALCESE_INITIAL_TWO-PHASE

$ z uic vic wic rho_1 rho_2 pic tic

1.0 0.0 0.0 0.0 1.0e+3 0.0 1.0e+5

$

*DUALCESE_INITIAL_TWO-PHASE_SET

$ ssid

61

$ z uic vic wic rho_1 rho_2 pic tic

1.0 -432.69 0.0 0.0 1.23e+3 0.0 1.0e+9

$

$

*DUALCESE_INITIAL_TWO-PHASE_SET

$ ssid

62

$ z uic vic wic rho_1 rho_2 pic tic

0.0 0.0 0.0 0.0 0.0 1.2 1.0e+5

$

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

$ Setup fluid properties

$

*DUALCESE_PART_MULTIPHASE

$ pid ReactinID EosSet_id mid FSI-type MOVMSHALG

1 11

*DUALCESE_EOS_SET

$ Setid inert_eosid reactant_eosid product_eosid

11 5 6

*DUALCESE_EOS_VAN_DER_WAALS_GENERALIZED

$ eosid a b GA BT

5 0.0 0.0 4.4 6.0e+8

*DUALCESE_EOS_VAN_DER_WAALS_GENERALIZED

$ eosid a b GA BT

6 5.0 1.0e-3 1.4 0.0

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

$ Handle output of state and restart data

$

*DUALCESE_D3PLOT

density

pressure

velocity

total_energy

internal_energy

temperature

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

*END