Basics : Cylinder flow FSI
This LS-DYNA simulation is a extension to the Basics : Cylinder flow where a structural part is added and the FSI problem is solved. The structural part is fixed in all direction except the vertical direction, the fluid deck is the same and the only modification is the time step to save computational time.
https://www.dynaexamples.com/icfd/basics-examples/Basics_Cylinder_flow_FSI
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Basics : Cylinder flow FSI
This LS-DYNA simulation is a extension to the Basics : Cylinder flow where a structural part is added and the FSI problem is solved. The structural part is fixed in all direction except the vertical direction, the fluid deck is the same and the only modification is the time step to save computational time.
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ $ FSI deck $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *KEYWORD *INCLUDE *ICFD_BOUNDARY_FSI *ICFD_CONTROL_FSI *END $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ $ Solid deck $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *KEYWORD *PARAMETER *TITLE *BOUNDARY_SPC_SET *CONTROL_IMPLICIT_DYNAMICS *CONTROL_IMPLICIT_GENERAL *CONTROL_IMPLICIT_SOLUTION *CONTROL_TERMINATION *ELEMENT_SHELL *INCLUDE *PART *SECTION_SHELL *SET_NODE_LIST *MAT_RIGID *NODE *END $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ $ Fluid deck $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *KEYWORD *TITLE *INCLUDE *DATABASE_BINARY_D3PLOT *DEFINE_CURVE_TITLE *DEFINE_CURVE_TITLE *ICFD_BOUNDARY_NONSLIP *ICFD_BOUNDARY_PRESCRIBED_VEL *ICFD_BOUNDARY_PRESCRIBED_PRE *ICFD_BOUNDARY_PRESCRIBED_TEMP *ICFD_CONTROL_OUTPUT *ICFD_CONTROL_TIME *ICFD_DATABASE_TIMESTEP *ICFD_DATABASE_TEMP *ICFD_INITIAL *ICFD_MAT *ICFD_PART *ICFD_PART_VOL *ICFD_SECTION *MESH_BL *MESH_BL_SYM *MESH_VOLUME *MESH_NODE *MESH_SURFACE_ELEMENT *END
$-----------------------------------------------------------------------------
$
$ Example provided by Marcus Timgren (DYNAmore Nordic)
$
$ 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 main_fsi.k and select memory=5m.
$X Requires LS-DYNA MPP R8.0.0 (or higher) with double precision.
$X------------------------------------------------------------------------------
$X------------------------------------------------------------------------------
$X
$X
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ main_fsi.k $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$-------------------------------------------------------------------------------
$
*KEYWORD
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ INCLUDE SOLID AND FLUID DECK
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*INCLUDE
struc.k
*INCLUDE
i.k
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ BOUNDARY CONDITIONS
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_BOUNDARY_FSI
$# pid
4
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD CONTROL CARDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_CONTROL_FSI
$# owc bt dt idc lcidsf
0 0.01.00000E28 0.25 0
*END
$
$
$-------------------------------------------------------------------------------
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ Solid Deck $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$-------------------------------------------------------------------------------
*KEYWORD
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ PARAMETERS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PARAMETER
$# prmr1 val1 prmr2 val2 prmr3 val3 prmr4 val4
R t_ends100.0
$--- Struc
$
Rrho_struc1.20
R e_struc2e11
R nu_struc0.3
Rdtm_struc0.05
$--- Struc
$
$--- Struc
$
$--- Struc
$
*INCLUDE
mesh_struc.k
*INCLUDE
boundary_condition.k
*INCLUDE
implicit.k
$--- Struc
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ CONTROL CARDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*CONTROL_TERMINATION
$# endtim endcyc dtmin endeng endmas
&t_ends 0 0.0 0.01.000000E8
*MAT_RIGID
$# mid ro e pr n couple m alias
1&rho_struc&e_struc &nu_struc 0.0 0.0 0.0
$# cmo con1 con2
0.0 0 0
$#lco or a1 a2 a3 v1 v2 v3
0.0 0.0 0.0 0.0 0.0 0.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ PARTS/ SECTION/ MATERIAL $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PART
$# title
$# pid secid mid eosid hgid grav adpopt tmid
1 1 1 0 0 0 0 0
*SECTION_SHELL
$# secid elform shrf nip propt qr/irid icomp setyp
1 12 1.0 2 1.0 0 0 1
$# t1 t2 t3 t4 nloc marea idof edgset
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0
*END
$
$-------------------------------------------------------------------------------
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ Fluid Deck $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$-------------------------------------------------------------------------------
*KEYWORD
*TITLE
$# title
ICFD Cylindrical Flow
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ PARAMETERS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PARAMETER
R t_end100.0
R dt_plot1.00
$
$--- Fluid
$
R v_inlet1.0
Rrho_fluid1.0
R mu_fluid0.005
R dt_fluid0.050
$
$--- Fluid
*TITLE
$# title
ICFD Cylinder Flow
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ INCLUDE SURFACE MESH
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*INCLUDE
mesh.k
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ DATABASE CARD FOR RESULTS OUTPUT
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*DATABASE_BINARY_D3PLOT
$# dt lcdt beam npltc psetid
&dt_plot 0 0 0 0
$# ioopt
0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ LOAD CURVES FOR THE BOUNDARY CONDITIONS
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*DEFINE_CURVE_TITLE
Velocity inlet
$# lcid sidr sfa sfo offa offo dattyp lcint
1 0 1.0&v_inlet 0.0 0.0 0 0
$# a1 o1
0.0 1.0
10000.0 1.0
*DEFINE_CURVE_TITLE
Pressure outlet
$# lcid sidr sfa sfo offa offo dattyp lcint
2 0 1.0 1.0 0.0 0.0 0 0
$# a1 o1
0.0 0.0
10000.0 0.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD BOUNDARY/INITIAL CONDITIONS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_BOUNDARY_FREESLIP
$# pid
3
*ICFD_BOUNDARY_NONSLIP
$# pid
4
*ICFD_BOUNDARY_PRESCRIBED_VEL
$# pid dof vad lcid sf vid death birth
1 1 1 1 1.0 01.00000E28 0.0
*ICFD_BOUNDARY_PRESCRIBED_VEL
$# pid dof vad lcid sf vid death birth
1 2 1 2 1.0 01.00000E28 0.0
*ICFD_BOUNDARY_PRESCRIBED_PRE
$# pid lcid sf death birth
2 2 1.01.00000E28 0.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD CONTROL CARDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_CONTROL_TIME
$# ttm dt cfl lcidsf dtmin dtmax
&t_end &dt_fluid 1.0 0 0.0 0.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ DATABASE (OUTPUT) $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_DATABASE_DRAG
$# pid cpid dtout
4 0 0.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD PARTS/ SECTION/ MATERIAL $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_MAT
$# mid flg ro vis st
1 1&rho_fluid&mu_fluid 0.0
*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
$# pid secid mid
4 1 1
*ICFD_PART_VOL
$# pid secid mid
10 1 1
$# spid1 spid2 spid3 spid4 spid5 spid6 spid7 spid8
1 2 3 4 0 0 0 0
*ICFD_SECTION
$# sid
1
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD MESH KEYWORDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MESH_BL
$# pid nelth blth blfe blst
4 2 0.0 0.0 0
*MESH_VOLUME
$# volid
1
$# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8
1 2 3 4 0 0 0 0
*END
This tutorial explains the structural and FSI part, for the fluid part check the Basics: Cylinder flow case.
Watch the tutorial on YouTube