This LS-DYNA simulation is a pipe junction with the heat equation activated, it has two inlets and one outlet. This model demonstrates how to use the boundary layer mesh generation and also how to output results in a Paraview format. The geometry is the same as in the pipe junction tutorial in ANSA a pre-processor from BETA.
*KEYWORD *TITLE *INCLUDE *DATABASE_BINARY_D3PLOT *DEFINE_CURVE_TITLE *DEFINE_CURVE_TITLE *ICFD_BOUNDARY_NONSLIP *ICFD_BOUNDARY_PRESCRIBED_VEL *ICFD_BOUNDARY_PRESCRIBED_PRE *ICFD_BOUNDARY_PRESCRIBED_TEMP *ICFD_CONTROL_OUTPUT *ICFD_CONTROL_TIME *ICFD_DATABASE_TIMESTEP *ICFD_DATABASE_TEMP *ICFD_INITIAL *ICFD_MAT *ICFD_PART *ICFD_PART_VOL *ICFD_SECTION *MESH_BL *MESH_BL_SYM *MESH_VOLUME *MESH_NODE *MESH_SURFACE_ELEMENT *END
$-----------------------------------------------------------------------------
$
$ 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
