Reduced Input
$================================================================================
*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
$================================================================================