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: Dimensionless
\$X------------------------------------------------------------------------------
\$X
*KEYWORD
*TITLE
CESE Shock Boundary layer interaction
*INCLUDE
mesh.k
\$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
\$                                                                              \$
\$                             PARAMETERS                                       \$
\$                                                                              \$
\$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PARAMETER
R    T_end        16
R  dt_plot       0.5
\$
\$--- Fluid
\$
R M_1           2.00
R P_1      0.1785714
R ro_1           1.0
R Cp           0.625
R Cv        0.446429
R dt_fluid    0.0001
Rcfl_fluid       0.9
Rmu_c1      6.566e-6
Rmu_c2        0.9436
RPrandtl         0.7
\$ The oblique shock relations then give the following parameters for a
\$ Betad=32.6 degree obloque shock wave :
*PARAMETER_EXPRESSION
Rg,Cp/Cv
*PARAMETER_EXPRESSION
Rgplus,g+1
*PARAMETER_EXPRESSION
Rgmin,g-1
*PARAMETER_EXPRESSION
*PARAMETER_EXPRESSION
Rsinb,sin(beta)
*PARAMETER_EXPRESSION
Rcotb,ctn(beta)
*PARAMETER_EXPRESSION
RMn1,M_1*sinb
*PARAMETER_EXPRESSION
Rro_2,ro_1*(gplus*(Mn1*Mn1))/(gmin*(Mn1*Mn1)+2)
*PARAMETER_EXPRESSION
RP_2,P_1+(P_1*2)*(g/gplus)*(Mn1*Mn1-1)
*PARAMETER_EXPRESSION
Rtantheta,2*cotb*((M_1*M_1)*(sinb*sinb)-1)/(2+M_1*M_1*(g+cos(2*beta)))
*PARAMETER_EXPRESSION
Rtheta,atan(tantheta)
*PARAMETER_EXPRESSION
RMn2,sqrt(((Mn1*Mn1)+2/gmin)/(2*(g/gmin)*(Mn1*Mn1)-1))
*PARAMETER_EXPRESSION
RM_2,Mn2/sin(beta-theta)
*PARAMETER_EXPRESSION
RU_2,M_2*sqrt(g*(P_2/ro_2))
*PARAMETER_EXPRESSION
RU_1,M_1*sqrt(g*(P_1/ro_1))
*PARAMETER_EXPRESSION
Ruy_2,-1*(U_2*sin(theta))
*PARAMETER_EXPRESSION
Rux_2,U_2*cos(theta)
\$
\$---+----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         0         2
*CESE_CONTROL_TIMESTEP
\$     iddt       cfl     dtint
2&cfl_fluid &dt_fluid
*CESE_CONTROL_LIMITER
\$    idlmt      alfa      beta      epsr
0       0.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         7         3
*CESE_EOS_IDEAL_GAS
\$    eosid        cv        cp
3       &Cv       &Cp
*CESE_MAT_GAS
\$      mid        c1        c2      prnd
7    &mu_c1    &mu_c2  &prandtl
\$---+----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_1       0.0       0.0     &ro_1      &P_1
*CESE_BOUNDARY_PRESCRIBED_SET
\$     ssid
5
\$   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
&Ux_2     &Uy_2       0.0     &ro_2      &P_2
*CESE_BOUNDARY_NON_REFLECTIVE_SET
\$     ssid
2
*CESE_BOUNDARY_NON_REFLECTIVE_SET
\$     ssid
4
*CESE_BOUNDARY_SOLID_WALL_SET
\$     ssid
3
*CESE_INITIAL
\$      uic       vic       wic     rhoic       pic       tic       hic
&Ux_2     &Uy_2       0.0     &ro_2      &P_2
\$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
\$                                                                              \$
\$                             DATABASE (OUTPUT)                                \$
\$                                                                              \$
\$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*DATABASE_BINARY_D3PLOT
&dt_plot
*END```