$-----------------------------------------------------------------------------
$
$ Example provided by Iñaki and Zheng (LSTC)
$
$ E-Mail: info@dynamore.de
$ Web: http://www.dynamore.de
$
$ Copyright, 2015 DYNAmore GmbH
$ Copying for non-commercial usage allowed if
$ copy bears this notice completely.
$
$X------------------------------------------------------------------------------
$X
$X 1. Run file shock_tube.k use LS-DYNA version R9.1 or later.
$X------------------------------------------------------------------------------
$X------------------------------------------------------------------------------
$X
$X
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ shock_tube.k $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$-------------------------------------------------------------------------------
*KEYWORD
*TITLE
$# title
Sod_1D_shock_tube
*CONTROL_STRUCTURED
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ DATABASE (OUTPUT) $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*DATABASE_BINARY_D3DUMP
$ dt/cycl lcdt beam npltc
$# cycl lcdt beam npltc psetid
100000.0 0 0 0 0
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$ Handle output of state and restart data
$
*DATABASE_BINARY_D3PLOT
$ dt/cycl lcdt beam npltc
$# dt lcdt beam npltc psetid
0.1 0 0 0 0
$# ioopt
0
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ Include the mesh $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*INCLUDE
$#
shock_tube_m1h_2d_CESE.k
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ CESE CONTROL CARDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*CONTROL_TERMINATION
$# endtim endcyc dtmin endeng endmas nosol
0.2 900000 0.0 0.01.000000E8 0
*CESE_CONTROL_SOLVER
$# icese iflow igeom iframe
0 1 2 0
*CESE_CONTROL_TIMESTEP
$# iddt cfl dtint
2 0.91.00000E-4
*CESE_CONTROL_LIMITER
$# idlmt alfa beta epsr
0 2.0 1.0 0.5
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ CESE BOUNDARY/INITIAL CONDITIONS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*CESE_BOUNDARY_NON_REFLECTIVE_SET
$# ssid
1
2
*CESE_BOUNDARY_REFLECTIVE_SET
$# ssid
3
5
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$ Setup the IC for fluid
$
*CESE_INITIAL
$# u v w rho p t
0.0 0.0 0.0 0.125 0.1 0.0
$
*CESE_INITIAL_SET
$# esid u v w rho p t
111 0.0 0.0 0.0 1.0 1.0 0.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ CESE PARTS/ EOS/ MATERIAL $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*CESE_PART
$# pid mid eosid
1 0 3
$*CESE_MAT_GAS
$$ mid c1 c2 prnd
$ 4 1.458e-6 110.4 0.72
*CESE_EOS_IDEAL_GAS
$# eosid cv cp
3 717.5 1004.5
*END
