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