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.
https://www.dynaexamples.com/icfd/basics-examples/anisotropic
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
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.
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.
Fluid velocity fringes
*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
$-----------------------------------------------------------------------------
$
$ 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