Advanced : Pillar submerged in deep sea
In this example, a partially submerged pillar vibrates due to the periodic incoming waves. Strong FSI coupling is mandatory since the densities of the solid and the fluid are very close.
https://www.dynaexamples.com/icfd/advanced-examples/pillar
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Advanced : Pillar submerged in deep sea
In this example, a partially submerged pillar vibrates due to the periodic incoming waves. Strong FSI coupling is mandatory since the densities of the solid and the fluid are very close.
Fluid pressure fringes
*KEYWORD *TITLE *CONTROL_IMPLICIT_GENERAL *CONTROL_IMPLICIT_DYNAMICS *CONTROL_IMPLICIT_SOLUTION *CONTROL_TERMINATION *CONTROL_TIMESTEP *DATABASE_BINARY_D3PLOT *DEFINE_CURVE *ELEMENT_SOLID *ICFD_BOUNDARY_FREESLIP *ICFD_BOUNDARY_NONSLIP *ICFD_BOUNDARY_FSWAVE *ICFD_CONTROL_OUTPUT *ICFD_CONTROL_TIME *ICFD_MAT *ICFD_PART *ICFD_PART_VOL *ICFD_SECTION *ICFD_SOLVER_TOL_MOM *ICFD_SOLVER_TOL_PRE *INCLUDE *LOAD_BODY *MAT_ELASTIC *MESH_SURFACE_ELEMENT *MESH_SURFACE_NODE *MESH_VOLUME *PART *SECTION_SOLID *END
$----------------------------------------------------------------------------- $ $ 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 R10.1 (or higher) with double precision $X $X------------------------------------------------------------------------------ $# UNITS: (kg/m/s) $X------------------------------------------------------------------------------ $X *KEYWORD *TITLE ICFD Partially submerged pillar $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ $ $ PARAMETERS $ $ $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *PARAMETER R T_end 40. R dt_plot 0.25 $ $--- Fluid $ Rrho_fluid 1000 R mu_fluid 0.001 R dt_fluid 0.035 R grav 9.81 R fsi_bt 0.07 *INCLUDE mesh.k *INCLUDE struc.k $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ $ $ ICFD CONTROL CARDS $ $ $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *ICFD_CONTROL_TIME $# ttm dt &T_end &dt_fluid *ICFD_CONTROL_FSI $# owc bt 0 &fsi_bt *ICFD_CONTROL_OUTPUT $# msglv 3 *ICFD_CONTROL_TURBULENCE 2 0.13 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ $ $ ICFD SOLVER TOLERANCES $ $ $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *ICFD_SOLVER_TOL_MOM 1.e-6 1.e-6 *ICFD_SOLVER_TOL_PRE 1.e-6 1.e-6 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ $ $ ICFD PARTS/ SECTION/ MATERIAL $ $ $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *ICFD_PART $# pid secid mid 1 1 4 *ICFD_PART $# pid secid mid 2 1 4 *ICFD_PART $# pid secid mid 3 1 4 *ICFD_PART $# pid secid mid 4 1 1 *ICFD_PART $# pid secid mid 5 1 4 *ICFD_PART $# pid secid mid 6 1 1 *ICFD_PART $# pid secid mid 7 1 1 *ICFD_PART $# pid secid mid 8 1 1 *ICFD_PART $# pid secid mid 9 1 1 *ICFD_PART $# pid secid mid 10 1 4 *ICFD_PART $# pid secid mid 11 1 4 *ICFD_PART_VOL 12 1 1 2, 4, 6, 7, 8, 9 *ICFD_PART_VOL 13 1 4 1, 2, 3, 5, 10, 11 *ICFD_MAT $# mid flg ro vis 1 0 *ICFD_MAT $# mid flg ro vis 4 1&rho_fluid &mu_fluid $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ $ $ ICFD BOUNDARY/INITIAL CONDITIONS $ $ $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *ICFD_BOUNDARY_FREESLIP 1 *ICFD_BOUNDARY_FREESLIP 3 *ICFD_BOUNDARY_NONSLIP 4 *ICFD_BOUNDARY_FSI 4 *ICFD_BOUNDARY_NONSLIP 5 *ICFD_BOUNDARY_FSI 5 *ICFD_BOUNDARY_FREESLIP 6 *ICFD_BOUNDARY_FSWAVE 7 2 2. 0.35 4. 20. *ICFD_BOUNDARY_FREESLIP 8 *ICFD_BOUNDARY_FREESLIP 9 *ICFD_BOUNDARY_FREESLIP 10 *ICFD_BOUNDARY_FSWAVE 11 2 2. 0.35 4. 20. *LOAD_BODY_Z $# lcid sf 100 1 *DEFINE_CURVE_TITLE Gravity force $# lcid sidr sfa sfo offa offo dattyp 100 &grav $# a1 o1 0.0 1.0 10000.0 1.0 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ $ $ ICFD MESH KEYWORDS $ $ $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *MESH_VOLUME 20 1, 3, 4, 5, 6, 7, 8, 9 10, 11 *MESH_INTERF 20 2 *MESH_BL_SYM 1 *MESH_BL 4,1 *MESH_BL 5,1 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ $ $ DATABASE (OUTPUT) $ $ $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *ICFD_DATABASE_DRAG 5 *ICFD_DATABASE_DRAG 4 *DATABASE_BINARY_D3PLOT &dt_plot *DEFINE_CURVE 1 , 0, 1. , 1. , 0. , 0. 0.,5 1e3,5 *DEFINE_CURVE 2 , 0, 1. , 1. , 0. , 0. 0.,0. 1e3,0. *END
In this example, a partially submerged pillar vibrates due to the periodic incoming waves. Strong FSI coupling is mandatory since the densities of the solid and the fluid are very close.