Moving Wedge Shock Wave Interaction
In this test case, a shock wave (initially located at the left boundary) is moving from left to right while a rigid body structural wedge is moving in the opposite direction. From the simulation results, we can see a complicated flow pattern because of the interactions between the off-body shock produced by the moving wedge, the upcoming shock, and the wedge.
*DUALCESE_MODEL
*DUALCESE_CONTROL_SOLVER
*DUALCESE_CONTROL_TIMESTEP
*DUALCESE_CONTROL_LIMITER
*DUALCESE_INCLUDE_MODEL
*DUALCESE_BOUNDARY_NON_REFLECTIVE_SEGMENT_SET
*DUALCESE_BOUNDARY_REFLECTIVE_SEGMENT_SET
*DUALCESE_INITIAL
*DUALCESE_PART
*DUALCESE_MAT_001
*DUALCESE_EOS_IDEAL_GAS
*DUALCESE_D3PLOT
*CONTROL_TERMINATION
*DATABASE_BINARY_D3PLOT
*DATABASE_BINARY_D3DUMP
*INCLUDE
*CONTROL_TIMESTEP
*DEFINE_CURVE_TITLE
*PART
*SECTION_SHELL_TITLE
*MAT_RIGID_TITLE
*BOUNDARY_PRESCRIBED_MOTION_RIGID_ID
*DEFINE_CURVE
*KEYWORD
*TITLE
fluid_2d_bullet_interaction
*CONTROL_TERMINATION
$ endtim endcyc dtmin endeng endmas
.004 200000
$===========================================================================
$ Units:
$ rho (density): Kg/m^3
$ V (velocity): m/s
$ P (pressure): Pa
$ T (temperature): K
$============================================================================
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ Handle output of state and restart data
*DATABASE_BINARY_D3PLOT
$ dt/cycl lcdt beam npltc
2.e-5 0
*DATABASE_BINARY_D3DUMP
$ dt/cycl lcdt beam npltc
10000
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*INCLUDE
struct_setup.k
$
*DUALCESE_MODEL
dcese_setup.k
*END
$---------------------------------------------------------------------------------------------------------------------------------
*KEYWORD
*INCLUDE
struct_mesh.k
$ Time step control for structure side
*CONTROL_TIMESTEP
$# dtinit tssfac isdo tslimt dt2ms lctm erode ms1st
1.0000E-7 0.900 0 0.000 0.000 0 0
$# dt2msf dt2mslc imscl unused unused rmscl
0.000 0 0 0.000
*DEFINE_CURVE_TITLE
minimum_time
$# lcid sidr sfa sfo offa offo dattyp
2 0 1.000 1.000 0.000 0.000 0
$# a1 o1
0.000 1.000000e-07
100.000000 1.000000e-07
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*PART
$# title
Bullet_part
$# pid secid mid eosid hgid grav adpopt tmid
2 2 1 0 0 0 0 0
*SECTION_SHELL_TITLE
bullet_motion
$# secid elform shrf nip propt qr/irid icomp setyp
2 12 1.000000 2 1 0 0 1
$# t1 t2 t3 t4 nloc marea idof edgset
1.000000 1.000000 1.000000 1.000000 0.000 0.000 0.000 0
*MAT_RIGID_TITLE
bullet_material
$# mid ro e pr n couple m alias
1 7850.0 2.10E+11 0.200 0.000 0.000 0.000
$# cmo con1 con2
1.000 5. 7.
$#lco or a1 a2 a3 v1 v2 v3
0.000 0.000 0.000 0.000 0.000 0.000
$# title
*BOUNDARY_PRESCRIBED_MOTION_RIGID_ID
$# id heading
1bullet_part
$# pid dof vad lcid sf vid death birth
2 1 0 3 1.000000 01.0000E+28 0.000
*DEFINE_CURVE
$# lcid sidr sfa sfo offa offo dattyp
3 0 1.00E-6 1.0 0.000 0.000 0
$# a1 o1
0.0 0.0
1.0 340.0
2.0 500.0
5.0 1000.0
7.5 1250.0
10.0 1500.0
1.0e+05 1500.0
*end
$---------------------------------------------------------------------------------------------------------------------------
*KEYWORD
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*DUALCESE_CONTROL_SOLVER
$ ieqns igeom iframeref
Euler 2D fixed
*DUALCESE_CONTROL_TIMESTEP
$ iddt cfl dtint
2 0.5 0.1e-6
*DUALCESE_CONTROL_LIMITER
$ idlmt alfa beta epsr
0 2.0 1.0 0.05
*DUALCESE_INCLUDE_MODEL
dcese_mesh.k
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ Setup the boundary conditions for fluid
$ open BCs
*DUALCESE_BOUNDARY_NON_REFLECTIVE_SEGMENT_SET
$ ssid
2
$
$ bottom solid BCs
$
*DUALCESE_BOUNDARY_REFLECTIVE_SEGMENT_SET
$ ssid
1
3
4
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ Setup the initial conditions for fluid
$
*DUALCESE_INITIAL
$ uic vic wic rhoic pic tic hic
0.0 0.0 0.0 1.12 101325. 300.0
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ Setup fluid properties
$
*DUALCESE_PART
$ pid mid eosid FSI-type MOVMSHALG
1 7 5 IBM
*DUALCESE_MAT_001
$ mid c1 c2 prnd
7 1.458e-6 110.4 0.72
*DUALCESE_EOS_IDEAL_GAS
$ eosid cv cp t0 v0
5 717.5 1004.5
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ Handle output of state and restart data
*DUALCESE_D3PLOT
density
pressure
temperature
velocity
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*END