$-----------------------------------------------------------------------------
$
$ Example provided by Dr. Nils Karajan and Bruno Boll (DYNAmore)
$
$ 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 R8.0.0 (or higher) with double precision
$X
$X------------------------------------------------------------------------------
$# UNITS: kg / m / s / N / Pa / Nm (J) / Pa*s
$X------------------------------------------------------------------------------
$X
*KEYWORD 100M
*TITLE
$# title
ICFD Tool Cooling
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ CONTROL PARAMETERS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PARAMETER
R T_end 20.0
R dt_plot 0.10
$
$--- Fluid
$
R v_inlet 1.5
R T_inlet 20.0
R dt_fluid 0.002
$
$--- Solid
$
R T_tool 50.0
R T_blnk 200.0
R T_surf 50.0
R dt_thrm 0.002
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ FLUID-STRUCTURE-INTERACTION AND HEAT TRANSFER $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*INCLUDE
main_solid.k
*INCLUDE
main_fluid.k
$
$--- Surface for heat transfer
$
*ICFD_BOUNDARY_CONJ_HEAT
$# pid
43
$
$--- Dummy-FSI to be included, even though physically not needed!
$--- If not included, the conjugate heat transfer does not work
$
*ICFD_CONTROL_FSI
$# owc bt dt
1
$
*ICFD_BOUNDARY_FSI
$# pid
43
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ CONTROL CARDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*CONTROL_TERMINATION
$# endtim endcyc dtmin endeng endmas nosol
&T_end
*CONTROL_TIMESTEP
$# dtinit tssfac isdo tslimt dt2ms lctm erode ms1st
&dt_thrm
*CONTROL_SOLUTION
$# soln nlq isnan lcint
2
*CONTROL_THERMAL_SOLVER
$# atype ptype solver cgtol gpt eqheat fwork sbc
1 11
$# msglvl maxitr abstol reltol omega tsf
*CONTROL_THERMAL_TIMESTEP
$# ts tip its tmin tmax dtemp tscp lcts
1.0&dt_thrm
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ DATABASE (OUTPUT) $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*DATABASE_FORMAT
$# iform ibinary
1
*DATABASE_TPRINT
$# dt binary lcur ioopt
&dt_plot 3
*DATABASE_BINARY_D3PLOT
$# dt lcdt nobeam npltc psetid
&dt_plot
*DATABASE_EXTENT_BINARY
$# neiph neips maxint strflg sigflg epsflg rltflg engflg
2 2 2 2
$# cmpflg ieverp beamip dcomp shge stssz n3thdt ialemat
1
$# nintsld pkp_sen sclp hydro msscl therm intout nodout
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ PARTS/ SECTION/ MATERIAL $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$--- Tool solid
$
*PART
Tool
$# pid secid mid eosid hgid grav adpopt tmid
10001 10001 10001 10001
$
*SECTION_SOLID
$# label elform aet
10001 1
$
*MAT_RIGID
$# mid ro e pr n couple m alias/re
10001 7.85E+03 2.1e+11 0.3
$# cmo con1 con2
1.0 7.0 7.0
$# lco a2 a3 v1 v2 v3
$
*MAT_THERMAL_ISOTROPIC
$# tmid tro tgrlc tgmult tlat hlat
10001
$# hc tc
4.80e+02 5.20e+01
$
$--- Tool skin
$
*PART
Outer skin cold
$# pid secid mid eosid hgid grav adpopt tmid
10003 10003 10001 10001
$
*PART
Outer skin hot
$# pid secid mid eosid hgid grav adpopt tmid
10004 10003 10001 10001
$
*PART
Pipe wall
$# pid secid mid eosid hgid grav adpopt tmid
10005 10003 10001 10001
$
*SECTION_SHELL
$# label elform shrf nip propt qr icomp setyp
10003 2
$# t1 t2 t3 t4 nloc marea idof edgset
0.5
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ BOUNDARY/INITIAL CONDITIONS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*CONSTRAINED_RIGID_BODIES
$# pidm pids iflag
10001 10003
10001 10004
10001 10005
$
$---
$
*INITIAL_TEMPERATURE_SET
$# nsid temp loc
10003&T_tool
$
*SET_NODE_GENERAL_TITLE
All nodes
$# sid da1 da2 da3 da4
10003
$# option
ALL
$
$---
$
*INITIAL_TEMPERATURE_SET
$# nsid temp loc
10004&T_blnk
$
*SET_NODE_GENERAL_TITLE
Outer skin hot
$# sid da1 da2 da3 da4
10004
$# option p1
PART 10004
$
$---
$
*INITIAL_TEMPERATURE_SET
$# nsid temp loc
10005&T_surf
*SET_NODE_GENERAL_TITLE
Pipe wall
$# sid da1 da2 da3 da4
10005
$# option p1
PART 10005
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ CONTROL CARDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_CONTROL_TIME
$# ttm dt cfl
&T_end &dt_fluid
$
*ICFD_CONTROL_OUTPUT
$# msgl outl dtout
2
$
*ICFD_CONTROL_SURFMESH
$# rsrf
1
$
$--- Stabilization of the temperature field to avoid negative temperatures
$
*ICFD_CONTROL_STAB
$# on/off
1
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ TURBULENCE MODELS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$# 0: Variational multiscale model (default)
$# 1: RANS k-eps model
$# 2: LES Smagorinsky sub-grid scale model
$# 3: LES Wall adapting local eddy-viscosity (WALE) model
*ICFD_CONTROL_TURBULENCE
$# tmod
2
$# ce1 ce2 qe qk cu
0.2
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ PARTS/ SECTION/ MATERIAL $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ 41 Inlet
$ 42 Outlet
$ 43 Pipe walls
$
*ICFD_PART
$# pid secid mid
41 1 1
$
*ICFD_PART
$# pid secid mid
42 1 1
$
*ICFD_PART
$# pid secid mid
43 1 1
$
$--- Combine the boundaryies of the fluid domain
$
*ICFD_PART_VOL
$# pid secid mid
44 1 1
$# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8
41 42 43
$
$--- Section needs to be generated but serves no further function at the moment
$
*ICFD_SECTION
$# secid
1
$
$--- Material definition of air
$
*ICFD_MAT
$# Source: www.engineeringtoolbox.com/air-properties-d_156.html @ 20 degrees
$# mid flg ro vis thd
1 1 1.205 15.11e-6 0.000
$# hc tc beta
1005.0 0.0257
$ 1005.0 0.0257 3.43e-3
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ BOUNDARY CONDITIONS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ 41 Inlet
$ 42 Outlet
$ 43 Pipe walls
$
$--- Stick condition for pipe walls
$
*ICFD_BOUNDARY_NONSLIP
$# pid
43
$
$--- Prescribe velocity at inlet
$
*ICFD_BOUNDARY_PRESCRIBED_VEL
$# x-velocity
$# pid dof vad lcid sf vid death birth
41 1 3 1
$
*ICFD_BOUNDARY_PRESCRIBED_VEL
$# y-velocity
$# pid dof vad lcid sf vid death birth
41 2 1 2
$
*ICFD_BOUNDARY_PRESCRIBED_VEL
$# z-velocity
$# pid dof vad lcid sf vid death birth
41 3 1 2
$
$--- Prescribe pressure at outlet
$
*ICFD_BOUNDARY_PRESCRIBED_PRE
$# pid lcid sf death birth
42 2
$
$--- Prescribe temperature at inlet
$
*ICFD_BOUNDARY_PRESCRIBED_TEMP
$# pid lcid sf death birth
41 3 1.0 1.0E+28 0.0
$
$--- Curve definitions for velocity/pressure/temperature
$
*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
$
*DEFINE_CURVE_TITLE
Temperature inlet
$# lcid sidr sfa sfo offa offo dattyp
3 &T_inlet
$# a1 o1
0.0 1.0
10000.0 1.0
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ INITIAL CONDITIONS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_INITIAL
$# pid vx vy vz T p
44 0.0 0.0 0.0&T_inlet
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ MESH / NODES/ ELEMENTS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$--- Surface definition for internal mesh generator
$
*MESH_VOLUME
$# volid
44
$# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8
41 42 43
$
$
$--- Boundary layer definition
$
*MESH_BL
$# pid nelth
43 2
$
*MESH_BL_SYM
$# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8
41 42
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ LSO - OUTPUT (SAMPLE POINTS) $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*LSO_TIME_SEQUENCE
$# solver
ICFD
$# dt lcdt lcopt npltc tbeg tend
0.01
$# domid1 domid2 domid3 domid4 domid5 domid6 domid7 domid8
142
$
$---
$
*LSO_DOMAIN
$# domain_type
thist_point
$# solver_name
ICFD
$# outid refid reduct
142 42
$# variable_name
pressure_point
velocity_point
temperature_point
$
$---
$
*LSO_POINT_SET
42
$ x y z
1.5886 -.4722 1.0965
1.6090 -.4714 1.0966
1.6284 -.4707 1.0967
1.6384 -.4861 1.0956
1.6306 -.5025 1.0935
1.6116 -.5032 1.0935
1.5886 -.5041 1.0934
1.6544 -.4696 1.0968
1.6544 -.5017 1.0936
*END
