Special methods are available in the ICFD solver that allows user to model problems of 3D fluid infiltration through a porous media matrix. This simplified 2D example features a typical Resin Transfer Molding (RTM) application example. The highly viscous fluid has been considered Newtonian and the influence of temperature has been neglected even if both those factors could also be taken into account if needed.
Fluid velocity fringes
*KEYWORD *TITLE *DATABASE_BINARY_D3PLOT *DEFINE_CURVE_TITLE *ICFD_BOUNDARY_PRESCRIBED_PRE *ICFD_CONTROL_OUTPUT *ICFD_CONTROL_POROUS *ICFD_CONTROL_TIME *ICFD_MAT *ICFD_MODEL_POROUS *ICFD_PART *ICFD_PART_VOL *ICFD_SECTION *INCLUDE *MESH_INTERF *MESH_SURFACE_ELEMENT *MESH_SURFACE_NODE *MESH_VOLUME *PARAMETER *END
$-----------------------------------------------------------------------------
$
$ Example provided by Rodrigo and 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 Dev 116424 (or higher) with double precision
$X
$X------------------------------------------------------------------------------
$# UNITS: (kg/m/s)
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
ICFD RTM example
*INCLUDE
mesh.k
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ PARAMETERS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PARAMETER
R T_end 3.0
R dt_plot 0.02
$
$--- Fluid
$
Rrho_fluid 920
R mu_fluid 0.055
R dt_fluid 0.0
R porosity 0.5
R Perm_kxx 1.6e-10
R Perm_kyy 2.0e-11
R Perm_kzz 0.0
R dt_sf 0.4
$---+----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 &dt_sf
*ICFD_CONTROL_OUTPUT
$# mslv
3
*ICFD_CONTROL_POROUS
1
$---+----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
0 100
*ICFD_MAT
$# mid flg ro vis
2 0
0 100
*ICFD_MODEL_POROUS
$# pmmoid pmid
100 4
$# por per/thx ff/thy thz pvlcidx pvlcidy pvlcidz
&porosity
$# kx' ky' kz'
&Perm_kxx &Perm_kyy &Perm_kzz
$#1-x1-pid 1-y2-pid 1-z 2-x 2-y 2-z
1.0000 0.0000 0 0.0000 1.0000 0.
*ICFD_PART
$# pid secid mid
1 1 1
*ICFD_PART
$# pid secid mid
2 1 1
*ICFD_PART
$# pid secid mid
3 1 2
*ICFD_PART_VOL
$# pid secid mid
10 1 1
$# spid1 spid2 spid3 spid4
1 2
*ICFD_PART_VOL
$# pid secid mid
20 1 2
$# spid1 spid2 spid3 spid4
2 3
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD BOUNDARY/INITIAL CONDITIONS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_BOUNDARY_PRESCRIBED_PRE
$# pid lcid sf death birth
1 1
*ICFD_BOUNDARY_PRESCRIBED_PRE
$# pid lcid sf death birth
3 2
*DEFINE_CURVE
$Pressure inlet
$# lcid sidr sfa sfo offa offo dattyp
1 0 1.000000 1.000000 0.000 0.000 0
$# a1 o1
0.000 1.0e+05
1.0e+10 1.0e+05
*DEFINE_CURVE
$# lcid sidr sfa sfo offa offo dattyp
2 0 1.000000 1.000000 0.000 0.000 0
$# a1 o1
0.000 0
1.0e+10 0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD MESH KEYWORDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MESH_VOLUME
$# volid
20
$# pid1 pid2 pid3 pid4
1 3
*MESH_INTERF
$# volid
20
$# pid1 pid2 pid3 pid4
2
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ DATABASE (OUTPUT) $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*DATABASE_BINARY_D3PLOT
&dt_plot
*END
Special methods are available in the ICFD solver that allows user to model problems of 3D fluid infiltration through a porous media matrix. This simplified 2D example features a typical Resin Transfer Molding (RTM) application example. The highly viscous fluid has been considered Newtonian and the influence of temperature has been neglected even if both those factors could also be taken into account if needed.