x

# Reduced Input

\$-----------------------------------------------------------------------------
\$
\$ Example provided by Iñaki (LSTC)
\$
\$ E-Mail: info@dynamore.de
\$ Web: http://www.dynamore.de
\$
\$ 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 Shock Diffraction around a corner
\$
\$
\$                          Normal shock wave moving
\$                                with velocity W
\$                            relative to laboratory
\$                                      |
\$                                      |
\$       2:  Mass motion                |           1: Stagnant gas ahead
\$               U_p                    |              of the moving shock
\$      relative to laboratory          |
\$           induced by                 |
\$           the moving                 |
\$           shock wave                 |
\$                                      |
\$ -------------------------------------|----------------------------------------
\$
\$ Velocity behind the sock relative to the shock is :
\$  U_2= W-U_p
\$
\$ Normal moving shock wave relations give :
\$
\$   Mach = sqrt [ gamma+1/(2*gamma) (P_2/P_1 -1 ) +1 ]
\$   ro_2= ro_1 * [1+(gamma+1)/(gamma-1)*P_2/P_1]/[(gamma+1)/(gamma-1)+P_2/P_1]
\$   U_p = W*(1-ro_1/ro_2)
\$   Mach = W / a_1
\$
\$ Perfect gaz law gives :
\$
\$   a_1 = sqrt (gamma P_1/ ro_1)
\$
\$ Therefore a moving normal shock wave gives :
\$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
\$                                                                              \$
\$                             PARAMETERS                                       \$
\$                                                                              \$
\$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PARAMETER
R    T_end      0.01
R  dt_plot    0.0002
\$
\$--- Fluid
\$
R Mach          1.20
R P_1        1.000e5
R ro_1         1.189
R Cp          1004.5
R Cv           717.5
R dt_fluid    0.0001
Rcfl_fluid       0.5
*PARAMETER_EXPRESSION
Rg,Cp/Cv
*PARAMETER_EXPRESSION
Ra_1_sq,g*(P_1/ro_1)
*PARAMETER_EXPRESSION
Ra_1,sqrt(a_1_sq)
*PARAMETER_EXPRESSION
RW,a_1*Mach
*PARAMETER_EXPRESSION
RP_2,P_1*(1+(2*g)/(g+1)*(Mach*Mach-1))
*PARAMETER_EXPRESSION
Rro_2,ro_1*(1+((g+1)/(g-1))*(P_2/P_1))/((P_2/P_1)+(g+1)/(g-1))
*PARAMETER_EXPRESSION
RU_p,W*(1-(ro_1/ro_2))
*PARAMETER_EXPRESSION
RU_2,(W-U_p)
*INCLUDE
mesh.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
&T_end
*CESE_CONTROL_SOLVER
\$   iframe     iflow     igeom
0         1         2
*CESE_CONTROL_TIMESTEP
\$     iddt       cfl     dtint
2&cfl_fluid &dt_fluid
*CESE_CONTROL_LIMITER
\$    idlmt      alfa      beta      epsr
0       4.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
1                   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
1
\$   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
&U_P       0.0       0.0     &ro_2      &P_2
*CESE_BOUNDARY_NON_REFLECTIVE_SET
\$     ssid
2
*CESE_BOUNDARY_SOLID_WALL_SET
\$     ssid
3
*CESE_BOUNDARY_SOLID_WALL_SET
4
*CESE_INITIAL
\$      uic       vic       wic     rhoic       pic       tic       hic
0.0       0.0       0.0     &ro_1      &P_1
\$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
\$                                                                              \$
\$                             DATABASE (OUTPUT)                                \$
\$                                                                              \$
\$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*DATABASE_BINARY_D3PLOT
&dt_plot
*END