Intermediate : Pressure wave in a straight cylindrical vessel
A pressure wave is sent down a cylinder and the deformations of the flexible structure are observed. This is a robustness verification test for strong FSI algorithms and methods.
https://www.dynaexamples.com/icfd/intermediate-examples/pwave
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Intermediate : Pressure wave in a straight cylindrical vessel
A pressure wave is sent down a cylinder and the deformations of the flexible structure are observed. This is a robustness verification test for strong FSI algorithms and methods.
*KEYWORD *TITLE *BOUNDARY_SPC_SET_BIRTH_DEATH *CONTROL_IMPLICIT_DYNAMICS *CONTROL_IMPLICIT_GENERAL *CONTROL_IMPLICIT_SOLUTION *CONTROL_TERMINATION *DATABASE_BINARY_D3PLOT *DEFINE_CURVE_TITLE *ELEMENT_SHELL *ICFD_BOUNDARY_FSI *ICFD_BOUNDARY_NONSLIP *ICFD_BOUNDARY_PRESCRIBED_VEL *ICFD_BOUNDARY_PRESCRIBED_PRE *ICFD_CONTROL_FSI *ICFD_CONTROL_OUTPUT
*ICFD_CONTROL_SURFMESH *ICFD_CONTROL_TIME *ICFD_DATABASE_FLUX *ICFD_MAT *ICFD_PART *ICFD_PART_VOL *ICFD_SECTION *INCLUDE *MAT_ELASTIC *MESH_BL *MESH_BL_SYM *MESH_SURFACE_ELEMENT *MESH_SURFACE_NODE *MESH_VOLUME *NODE *PARAMETER *PART *SECTION_SHELL *SET_NODE_LIST *END
$-----------------------------------------------------------------------------
$
$ Example provided by Iñaki (LSTC)
$
$
$X------------------------------------------------------------------------------
$X
$X 1. Run file as is.
$X Requires LS-DYNA MPP R10.0.0 (or higher) with double precision
$X
$X------------------------------------------------------------------------------
$# UNITS: m/ks/s
*KEYWORD
*include
struc.k
*include
mesh_shell.k
*include
mesh.k
$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ FLUID $
$$$$$$$$$$$$$$$$$$$$$$$$$$$
*ICFD_CONTROL_TIME
0.02, 0.0002
*ICFD_CONTROL_SURFMESH
2
*ICFD_CONTROL_OUTPUT
3
*ICFD_BOUNDARY_PRESCRIBED_PRE
1,1
*ICFD_BOUNDARY_PRESCRIBED_PRE
2,2
*DEFINE_CURVE_TITLE
Pressure inflow
$# lcid sidr sfa sfo offa offo dattyp lcint
1 0 1.0 0.1 0.0 0.0 0 0
$# a1 o1
0.0 13332.0
0.005 13332.0
0.0054 0.0
10.0 0.0
*DEFINE_CURVE_TITLE
Pressure outflow
$# lcid sidr sfa sfo offa offo dattyp lcint
2 0 1.0 1 0.0 0.0 0 0
$# a1 o1
0.0 0
10.0 0
*ICFD_BOUNDARY_NONSLIP
$# pid
3
*ICFD_BOUNDARY_FSI
$# pid
3
*ICFD_PART
$# pid secid mid
1 1 1
*ICFD_PART
$# pid secid mid
2 1 1
*ICFD_PART
$# pid secid mid
3 1 1
*ICFD_PART_VOL
$# pid secid mid
5 1 1
$# spid1 spid2 spid3 spid4 spid5 spid6 spid7 spid8
1 2 3 0 0 0 0 0
*ICFD_MAT
$# mid flg ro vis st
1 1 1000 0.003 0.0
*MESH_VOLUME
10
1,2,3
*MESH_BL
3,2
*MESH_BL_SYM
1
*MESH_BL_SYM
2
*END
References :
L. Formaggia, J.-F. Gerbeau, F. Nobile, and A. Quarteroni. On the coupling of 3D and 1D
Navier-Stokes equations for flow problems in compliant vessels. Comp. Meth. Appl. Mech.
Engrg., 191(6-7):561–582, 2001.
M.A. Fern´andez and M. Moubachir. A Newton method using exact Jacobians for solving
fluid-structure coupling. Comp. & Struct., 83:127–142, 2005.
L. Formaggia, J.-F. Gerbeau, F. Nobile, and A. Quarteroni. On the coupling of 3D and 1D
Navier-Stokes equations for flow problems in compliant vessels. Comp. Meth. Appl. Mech.
Engrg., 191(6-7):561–582, 2001.