Shock Bubble interaction (+FSI) (+3D) — Welcome to LS-DYNA Examples

This FSI case is an expansion of the two-phase shock bubble example and adds a flexible structure (baffle) to the problem. In this test case, a left moving shock wave propagates and first hits a stationary bubble, which gains speed and loses its circular shape. After the shock wave hits the bubble, the created circular pressure waves keep moving forward and hit the structure panel causing its deformation. A 2D version and a 3D version of the problem are provided.

This FSI case is an expansion of the two-phase shock bubble example and add a flexible structure (baffle) to the problem. In this test case, a left moving shock wave propagates and first hits a stationary bubble, which gains speed and loses its circular shape. After the shock wave hits the bubble, the created circular pressure waves keep moving forward and hit the structure panel causing its deformation. A 2D and 3D version of the problem are provided.

References:

[1] G. Allaire, S. Clerc, S. Kokh,  A five-equation model for the simulation of interfaces between compressible fluids,  J. Comput. Phys. 181(2) (2002) 577-616

$================================================================================

*KEYWORD

*TITLE

shock_bubble_baffle_interact_2D

$===========================================================================

$ Units:

$ rho (density):  Kg/m^3

$ V   (velocity): m/s

$ P   (pressure): Pa

$ T   (temparature): K

$============================================================================

$

*CONTROL_TERMINATION

$   endtim    endcyc     dtmin     endeng   endmas

6.0e-4     20000

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

$ Handle output of state and restart data

$

*DATABASE_BINARY_D3PLOT

$  dt/cycl      lcdt      beam     npltc

5.0e-5                   0

*DATABASE_BINARY_D3DUMP

$  dt/cycl      lcdt      beam     npltc

10000

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

*INCLUDE

struct_setup.k

$

*DUALCESE_MODEL

fluid_setup.k

*END

$================================================================================

*KEYWORD

*TITLE

struct_shell

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

*PART

shell

11,6,5,0,0,0,0,0

*SECTION_SHELL

6,12,0.0,0.0,0.0

$ 0.0,0.0,0.0,0.0

0.009,0.009,0.009,0.009

*MAT_ELASTIC

5,5.2e+3,2.0e+13,0.1,0.0,0.0,0.0

$——————————————————————————-

*CONTROL_CONTACT

0.10000

 

 

2

*BOUNDARY_SPC_SET

79,0,1,1,1,1,1,1

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

*INCLUDE

struct_shell3.k

*CONTROL_STRUCTURED

*END

$================================================================================

*KEYWORD

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

*DUALCESE_CONTROL_SOLVER

$    icese      igeom   iframe MixSelect

Euler        2D     fixed two-phase

$     Euler      Axis     fixed two-phase

*DUALCESE_CONTROL_TIMESTEP

$     iddt       cfl     dtint

1       0.5     1.e-7

*DUALCESE_CONTROL_LIMITER

$    idlmt      alfa      beta      epsr

2      2.       1.0      .0

$

*DUALCESE_INCLUDE_MODEL

fluid_quad2.k

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

$ Setup the boundary conditions for fluid

$

$ inlet & Outlet

$

*DUALCESE_BOUNDARY_NON_REFLECTIVE_SEGMENT_SET

$     ssid

$         1

2

$

$ Solid BCs (up & down)

$

*DUALCESE_BOUNDARY_REFLECTIVE_SEGMENT_SET

$     ssid

3

5

1

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

$ Setup the initial conditions for fluid

$

*DUALCESE_INITIAL_TWO-PHASE

$        z      uic       vic       wic      rho_1     rho_2       pic       tic

1.0      0.0       0.0       0.0     1.0e+3       0.0    1.0e+5

$

 

*DUALCESE_INITIAL_TWO-PHASE_SET

$     ssid

61

$        z      uic       vic       wic      rho_1     rho_2       pic       tic

1.0  -432.69       0.0       0.0    1.23e+3       0.0    1.0e+9

$

$

*DUALCESE_INITIAL_TWO-PHASE_SET

$     ssid

62

$        z      uic       vic       wic      rho_1     rho_2       pic       tic

0.0      0.0       0.0       0.0        0.0       1.2    1.0e+5

$

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

$ Setup fluid properties

$

*DUALCESE_PART_MULTIPHASE

$      pid ReactinID EosSet_id       mid  FSI-type MOVMSHALG

1                  11            ibm

*DUALCESE_EOS_SET

$    Setid inert_eosid reactant_eosid product_eosid

11        5         6

*DUALCESE_EOS_VAN_DER_WAALS_GENERALIZED

$    eosid         a         b        GA        BT

5       0.0       0.0       4.4    6.0e+8

*DUALCESE_EOS_VAN_DER_WAALS_GENERALIZED

$    eosid         a         b        GA        BT

6       5.0    1.0e-3       1.4       0.0

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

$ Handle output of state and restart data

$

*DUALCESE_D3PLOT

density

pressure

velocity

total_energy

internal_energy

temperature

volume_fraction

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

*END

$================================================================================

– DUALCESE_example_FSIshockbubble.zip

– DUALCESE_example_3DFSIshockbubble.zip