Reduced Input
*KEYWORD
*TITLE
Rate Stick in 3D container
$
$===========================================================================
$ Units:
$ rho (density): Kg/m^3
$ V (velocity): m/s
$ P (pressure): Pa
$ T (temparature): K
$============================================================================
*CONTROL_TERMINATION
$ endtim endcyc dtmin endeng endmas
15.e-6 20000
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$ Handle output of state and restart data
$
*DATABASE_BINARY_D3PLOT
$ dt/cycl lcdt beam npltc
5.0e-7 0
*DATABASE_BINARY_D3DUMP
$ dt/cycl lcdt beam npltc
10000
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*INCLUDE
struct_setup.k
$
*DUALCESE_MODEL
dcese_setup.k
*END
$---------------------------------------------------------------------------------------------------------------------------------------------
*KEYWORD
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DUALCESE_CONTROL_SOLVER
$ icese igeom iframe MixSelect
Euler 3D fixed hybrid
*DUALCESE_CONTROL_TIMESTEP
$ iddt cfl dtint
2 .9 9.e-9
*DUALCESE_CONTROL_LIMITER
$ idlmt alfa beta epsr
2 1.0 1.0 .0
$
*DUALCESE_INCLUDE_MODEL
fluid_mesh1.k
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$ Setup the boundary conditions for fluid
$
*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
5
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$ Setup the initial conditions for fluid
$
*DUALCESE_INITIAL_HYBRID
$ z1 ra uic vic wic rho_1 rho_a rho_b
0.9999999 0.0 0.0 0.0 0.0 1.225 1905.0 1905.0
$ pic tic
1.01325e5
$ 1.1236e8
$
*DUALCESE_INITIAL_HYBRID_SET
$ ssid idef_func
62
$ z1 ra uic vic wic rho_1 rho_a rho_b
1.0e-7 1.0 0.0 0.0 0.0 1905.0 1905.0 1905.0
$ pic tic
0.33708e11
*DUALCESE_INITIAL_HYBRID_SET
$ ssid idef_func
63
$ z1 ra uic vic wic rho_1 rho_a rho_b
1.0e-7 1.0 0.0 0.0 0.0 1905.0 1905.0 1905.0
$ pic tic
1.1236e8
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$ Setup fluid properties
$
*DUALCESE_PART_MULTIPHASE
$ pid ReactinID EosSet_id mid FSI-type MOVMSHALG
$ 1 3 11
1 3 11 ibm
*DUALCESE_EOS_SET
$ Setid inert_eosid react_eosid product_eosid
11 5 6 7
*DUALCESE_EOS_IDEAL_GAS
$ eosid cv cp t0 v0
5 713.5 1001.5
*DUALCESE_EOS_JWL
$ eosid A B R1 R2 Gama0 rho0 e0
6 778.1e11 -5.031e9 11.3 1.13 0.8938 1905.0 0.2552e6
$ Cv
1.3055e3
*DUALCESE_EOS_JWL
$ eosid A B R1 R2 Gama0 rho0 e0
714.8105e11 0.6379e11 6.2 2.2 0.5 1905.0 0.3622e7
$ Cv
.524934e3
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*DUALCESE_REACTION_RATE_IG
$ react_id Ign aa bb xx Grow1 cc dd
3 4.0e12 0.22 0.667 7.0 4.5e-24 0.667 1.0
$ yy Grow2 ee gg zz IGmax G1max G2min
3.0 3.0e-4 0.667 0.667 1.0 0.02 0.8 0.8
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$ Handle output of state and restart data
$
*DUALCESE_D3PLOT
density
pressure
velocity
total_energy
internal_energy
temperature
volume_fraction
reactant_mass_fraction
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*END
$---------------------------------------------------------------------------------------------------------------------------------------------
*KEYWORD
*TITLE
struct_shell
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*PART
solid
$ PID SECID MID EOSID HGID GRAV ADPOPT TMID
11 6 5 0 1
*SECTION_SOLID
6 1
*MAT_ELASTIC
5,1.0e+4,2.0e+11,0.1,0.0,0.0,0.0
*HOURGLASS
1 4
$-------------------------------------------------------------------------------
*CONTROL_CONTACT
0.10000
2
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*INCLUDE
struct_mesh1.k
*CONTROL_STRUCTURED
*END