Basics : Cylinder flow FSI — Welcome to LS-DYNA Examples

Home » Knowledge Base » ICFD / Incompressible CFD » Basics Examples » 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.

Keywords

$---+----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

Animated Result

FSI_animation

 

Reduced Input

$-----------------------------------------------------------------------------
$
$ 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
 

Geometry of fluid domain

Model geometry

Download

Picture Results

FSI_image

Tutorial

This tutorial explains the structural and FSI part, for the fluid part check the Basics: Cylinder flow case.

Watch the tutorial on YouTube