Basics : Anisotropic Porous Media flow — Welcome to LS-DYNA Examples

Home » Knowledge Base » ICFD / Incompressible CFD » Basics Examples » Basics : Anisotropic Porous Media flow

This LS-DYNA simulation shows a simple ICFD anisotropic porous media flow input deck. The flow entering the porous media is deflected and follows the orientation of the porous media.

Description

This LS-DYNA simulation shows a simple ICFD anisotropic porous media flow input deck. The flow entering the porous media is deflected and follows the orientation of the porous media, the porous media is located between the two rectangular blocks.

Animated Result

Fluid velocity fringes

ezgif.comvideotogif.gif

 

Keywords

*DATABASE_BINARY_D3PLOT
*DEFINE_CURVE_TITLE
*END
*ICFD_BOUNDARY_FREESLIP
*ICFD_BOUNDARY_PRESCRIBED_PRE
*ICFD_BOUNDARY_PRESCRIBED_VEL
*ICFD_BOUNDARY_NONSLIP
*ICFD_CONTROL_TIME
*ICFD_MAT
*ICFD_MODEL_POROUS
*ICFD_PART
*ICFD_PART_VOL
*ICFD_SECTION
*INCLUDE
*KEYWORD
*MESH_INTERF
*MESH_SURFACE_ELEMENT
*MESH_SURFACE_NODE
*MESH_VOLUME
*PARAMETER
*TITLE

Reduced Input

$-----------------------------------------------------------------------------
$
$ 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 R9.0.0 (or higher) with double precision 
$X
$X------------------------------------------------------------------------------
$# UNITS: Dimensionless. 
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
ICFD Anisotropic Porous media flow
*INCLUDE
mesh.k
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                             PARAMETERS                                       $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PARAMETER
R    T_end      50.0
R  dt_plot      1.00
$
$--- Fluid
$
R  v_inlet       1.0
Rrho_fluid       1.0
R mu_fluid     0.005 
R dt_fluid     0.000
Rpor_fluid       0.8
R ff_fluid       0.1
R kx_fluid         1
R ky_fluid    0.0001
R kz_fluid    0.0001
$
$---+----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
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                       ICFD PARTS/ SECTION/ MATERIAL                          $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_SECTION
$#     sid   
         1
*ICFD_MAT
$#     mid       flg        ro       vis          
         1         1&rho_fluid &mu_fluid 
*ICFD_MAT
$#     mid       flg        ro       vis          
         2         1&rho_fluid &mu_fluid 
$#      hc        tc

$#  nnmoid    pmmoid
                   1
*ICFD_MODEL_POROUS
$#  pmmoid      pmid
         1         4
$#     por       per        ff
&por_fluid           &ff_fluid
$permeability vector in ppal axes (x',y',z')    
$#      kx        ky        kz
 &kx_fluid &ky_fluid &kz_fluid
$#     1-X       1-Y       1-Z       2-X       2-Y       2-Z
   0.93970  -0.34202        0.   0.34202   0.93970        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
$#     pid     secid       mid 
         6         1         2
*ICFD_PART
$#     pid     secid       mid 
         7         1         2
*ICFD_PART_VOL
$#     pid     secid       mid 
        10         1         1
$#   spid1     spid2     spid3     spid4     spid5
         1         2         3         4         6
*ICFD_PART_VOL
$#     pid     secid       mid 
        20         1         2
$#   spid1     spid2     
         6         7
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                    ICFD BOUNDARY/INITIAL CONDITIONS                          $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_BOUNDARY_PRESCRIBED_VEL
$#     pid       dof       vad      lcid       
         1         1         1         1      
*ICFD_BOUNDARY_PRESCRIBED_VEL
$#     pid       dof       vad      lcid       
         1         2         1         2       
*ICFD_BOUNDARY_PRESCRIBED_PRE
$#     pid      lcid        sf     death     birth     
         2         2     
*ICFD_BOUNDARY_FREESLIP
$#     pid   
         3
*DEFINE_CURVE_TITLE
Velocity inlet
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp
         1                    &v_inlet
$#                a1                  o1
                 0.0                 1.0
             10000.0                 1.0
*DEFINE_CURVE_TITLE
Pressure outlet
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp
         2
$#                a1                  o1
                 0.0                 0.0
             10000.0                 0.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                            ICFD MESH KEYWORDS                                $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MESH_VOLUME
$#   volid     
        10
$#    pid1      pid2      pid3      pid4      pid5      
         1         2         3         4         7
*MESH_INTERF
$#    vpid     
        10
$#     pid
         6
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                             DATABASE (OUTPUT)                                $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*DATABASE_BINARY_D3PLOT
$#      dt
  &dt_plot
*END
 

Figures

Cylinder Flow

 

Download