x
Our website uses cookies. By using the website you agree ot its use. More information can be found in our privacy policy.

Reduced Input

$-----------------------------------------------------------------------------
$
$ 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 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