This examples features a simple CESE input deck coupled with a structural problem. An incoming pressure wave hits a thin plate. Under the load, the elements get progressively eroded allowing the fluid to flow through.
CESE Pressure fringes
*CESE_BOUNDARY_NON_REFLECTIVE_SET *CESE_BOUNDARY_PRESCRIBED_SET *CESE_BOUNDARY_SOLID_WALL_SET *CESE_CONTROL_LIMITER *CESE_CONTROL_SOLVER *CESE_CONTROL_TIMESTEP *CESE_EOS_IDEAL_GAS *CESE_INITIAL *CESE_PART *CONTROL_TERMINATION *DATABASE_BINARY_D3PLOT *DEFINE_CURVE *ELEMENT_SHELL *ELEMENT_SOLID *KEYWORD *MAT_PLASTIC_KINEMATIC *MAT_RIGID *NODE *PARAMETER *PART *TITLE *SECTION_SHELL *SET_PART_LIST *SET_SEGMENT
$-----------------------------------------------------------------------------
$
$ Example provided by Iñaki (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 as is.
$X Requires LS-DYNA MPP R8.0.0 (or higher) with double precision
$X
$X------------------------------------------------------------------------------
$# UNITS: (kg/m/s)
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
CESE Erosion example
*INCLUDE
mesh.k
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ CESE CONTROL CARDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*CESE_CONTROL_SOLVER
$ iframe iflow igeom
200 1 3
*CESE_CONTROL_TIMESTEP
$ iddt cfl dtint
2&cfl_fluid &dt_fluid
*CESE_CONTROL_LIMITER
$ idlmt alfa beta epsr
0 2.0 1.0 0.5
$
$---+----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
22 3
*CESE_EOS_IDEAL_GAS
$ eosid cv cp
3 &Cv &Cp
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ CESE BOUNDARY/INITIAL CONDITIONS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*CESE_BOUNDARY_PRESCRIBED_SET
$ SSID
101
$ LCID_U LCID_V LCID_W LCID_D LCID_P LCID_T
-1
$ SF_U SF_V SF_W SF_D SF_P SF_T
0.0 0.0 &W &ro_2 &P_2
$
$ OUTFLOW BC (RIGHT)
$
*CESE_BOUNDARY_NON_REFLECTIVE_SET
$ SSID
102
$ SOLID BC (other sides)
*CESE_BOUNDARY_SOLID_WALL_SET
103
104
105
106
*CESE_INITIAL
$ uic vic wic rhoic pic tic hic
0.0 0.0 0.0 &ro_1 &P_1
*END