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
