Shock Bubble interaction (+FSI) (+3D)
This FSI case is an expansion of the two-phase shock bubble example and add a flexible structure (baffle) to the problem. In this test case, a left moving shock wave propagates and first hits a stationary bubble, which gains speed and loses its circular shape. After the shock wave hits the bubble, the created circular pressure waves keep moving forward and hit the structure panel causing its deformation. A 2D and 3D version of the problem are provided.
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
*INCLUDE
*PART
*SECTION_SHELL
*MAT_ELASTIC
*CONTROL_CONTACT
*BOUNDARY_SPC_SET
*CONTROL_STRUCTURED
*END
$================================================================================
*KEYWORD
*TITLE
shock_bubble_baffle_interact_2D
$===========================================================================
$ Units:
$ rho (density): Kg/m^3
$ V (velocity): m/s
$ P (pressure): Pa
$ T (temparature): K
$============================================================================
$
*CONTROL_TERMINATION
$ endtim endcyc dtmin endeng endmas
6.0e-4 20000
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$ Handle output of state and restart data
$
*DATABASE_BINARY_D3PLOT
$ dt/cycl lcdt beam npltc
5.0e-5 0
*DATABASE_BINARY_D3DUMP
$ dt/cycl lcdt beam npltc
10000
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*INCLUDE
struct_setup.k
$
*DUALCESE_MODEL
fluid_setup.k
*END
$================================================================================
*KEYWORD
*TITLE
struct_shell
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*PART
shell
11,6,5,0,0,0,0,0
*SECTION_SHELL
6,12,0.0,0.0,0.0
$ 0.0,0.0,0.0,0.0
0.009,0.009,0.009,0.009
*MAT_ELASTIC
5,5.2e+3,2.0e+13,0.1,0.0,0.0,0.0
$-------------------------------------------------------------------------------
*CONTROL_CONTACT
0.10000
2
*BOUNDARY_SPC_SET
79,0,1,1,1,1,1,1
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*INCLUDE
struct_shell3.k
*CONTROL_STRUCTURED
*END
$================================================================================
*KEYWORD
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DUALCESE_CONTROL_SOLVER
$ icese igeom iframe MixSelect
Euler 2D fixed two-phase
$ Euler Axis fixed two-phase
*DUALCESE_CONTROL_TIMESTEP
$ iddt cfl dtint
1 0.5 1.e-7
*DUALCESE_CONTROL_LIMITER
$ idlmt alfa beta epsr
2 2. 1.0 .0
$
*DUALCESE_INCLUDE_MODEL
fluid_quad2.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
$ ssid
3
5
1
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$ 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 ibm
*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
volume_fraction
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*END
$================================================================================