Sod shock tube example
This is a classic 1-D model, introduced by G.A. Sod [2] and its purpose is to verify the ability of the CESE solver to solve fluid dynamics problems with shock wave behavior. [2] G. Sod, A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws, J. Comput. Phys., (1978).
https://www.dynaexamples.com/cese/oricese/sod_shock_tube
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Sod shock tube example
This is a classic 1-D model, introduced by G.A. Sod [2] and its purpose is to verify the ability of the CESE solver to solve fluid dynamics problems with shock wave behavior. [2] G. Sod, A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws, J. Comput. Phys., (1978).
*KEYWORD *TITLE *CONTROL_STRUCTURED *CONTROL_TERMINATION *DATABASE_BINARY_D3DUMP *DATABASE_BINARY_D3PLOT *INCLUDE *CESE_CONTROL_SOLVER *CESE_CONTROL_TIMESTEP *CESE_CONTROL_LIMITER *CESE_BOUNDARY_NON_REFLECTIVE_SET *CESE_BOUNDARY_REFLECTIVE_SET *CESE_INITIAL *CESE_INITIAL_SET *CESE_PART *CESE_MAT_GAS *CESE_EOS_IDEAL_GAS *ELEMENT_SOLID *NODE *SET_SEGMENT *SET_SEGMENT *SET_SEGMENT *SET_SEGMENT *SET_SOLID *END
Animation of the pressure distribution:
$-----------------------------------------------------------------------------
$
$ 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