$-----------------------------------------------------------------------------
$
$ 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_pipe.k and select memory=1m.
$X Requires LS-DYNA MPP R9.0.0 (or higher) with double precision.
$X------------------------------------------------------------------------------
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
$# title
ICFD 3D pipe junction
$---+----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
0.05 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
pressure
$# lcid sidr sfa sfo offa offo dattyp lcint
1 0 1.0 1.0 0.0 0.0 0 0
$# a1 o1
0.0 0.0
10.0 0.0
*DEFINE_CURVE_TITLE
velocity
$# lcid sidr sfa sfo offa offo dattyp lcint
2 0 1.0 1.0 0.0 0.0 0 0
$# a1 o1
0.0 1.0
10.0 1.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ INCLUDE SURFACE MESH
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*INCLUDE
$#
pipe_junction_shell_mesh.k
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ BOUNDARY CONDITIONS
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_BOUNDARY_NONSLIP
$# pid
4
$Velocity inlets conditions
*ICFD_BOUNDARY_PRESCRIBED_VEL
$# pid dof vad lcid sf vid death birth
1 1 1 2 -1.0 01.00000E28 0.0
*ICFD_BOUNDARY_PRESCRIBED_VEL
$# pid dof vad lcid sf vid death birth
2 4 1 2 1.0 01.00000E28 0.0
$Outlet zero pressure condition
*ICFD_BOUNDARY_PRESCRIBED_PRE
$# pid lcid sf death birth
3 1 1.01.00000E28 0.0
$Temperature inlets conditions
*ICFD_BOUNDARY_PRESCRIBED_TEMP
$# pid lcid sf death birth
1 2 20.01.00000E28 0.0
*ICFD_BOUNDARY_PRESCRIBED_TEMP
$# pid lcid sf death birth
2 2 60.01.00000E28 0.0
$Temperature condition for the pipe wall
*ICFD_BOUNDARY_PRESCRIBED_TEMP
$# pid lcid sf death birth
4 2 20.01.00000E28 0.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD CONTROL CARDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$Outputs the results in a Paraview format and outputs the mesh in LS-PrePost format.
*ICFD_CONTROL_OUTPUT
$# msgl outl dtout lsppout
4 7 0.05 1
$Controls the termination time and time step
*ICFD_CONTROL_TIME
$# ttm dt cfl lcidsf dtmin dtmax
2.0 0.0 1.0 0 0.0 0.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ ICFD DATABASE CARDS
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$Outputs a file with time step information
*ICFD_DATABASE_TIMESTEP
$# outlv
1
$Outputs a file with temperature for the outlet
*ICFD_DATABASE_TEMP
$# pid
3
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ PART/Initial conditions/Fluid properties
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$Sets the initial temperature to 20 degrees celsius
*ICFD_INITIAL
$# pid vx vy vz t p
0 0.0 0.0 0.0 20.0 0.0
*ICFD_MAT
$# mid flg ro vis st
1 1 998.2 0.001005 0.0
$# hc tc beta prt hcsflcid tcsflcid
4181.0 0.597 0.0 0.85 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
5 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
$
$ MESH
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MESH_BL
$# pid nelth blth blfe blst
4 7 0.0022.00000E-4 2
*MESH_BL_SYM
$# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8
1 2 3 0 0 0 0 0
*MESH_VOLUME
$# volid
5
$# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8
1 2 3 4 0 0 0 0
*END
