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

$ Solid BCs (up & down)

*DUALCESE_BOUNDARY_REFLECTIVE_SEGMENT_SET

$*DUALCESE_BOUNDARY_SOLID_WALL_SEGMENT_SET

$ ssid

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

$ 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

$ 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

$ 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

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

*INCLUDE

struct_mesh1.k

*CONTROL_STRUCTURED

*END