Basics : Conjugate heat transfer
This LS-DYNA simulation shows a simple conjugate heat transfer problem. The default coupling method between the structure and the fluid is a robust monolithic strong coupling approach which allows for accurate solutions. In order to trigger the conjugate heat transfer solver, the FSI keywords must be defined as well, since the FSI algorithm is needed in order to detect which nodes belong to the fluid/solid interface. Of course, the corresponding solid thermal input deck must be set up.
https://www.dynaexamples.com/icfd/basics-examples/conj
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Basics : Conjugate heat transfer
This LS-DYNA simulation shows a simple conjugate heat transfer problem. The default coupling method between the structure and the fluid is a robust monolithic strong coupling approach which allows for accurate solutions. In order to trigger the conjugate heat transfer solver, the FSI keywords must be defined as well, since the FSI algorithm is needed in order to detect which nodes belong to the fluid/solid interface. Of course, the corresponding solid thermal input deck must be set up.
Fluid velocity fringes
*TITLE *CONTROL_TERMINATION *CONTROL_THERMAL_TIMESTEP *CONTROL_THERMAL_SOLVER *CONTROL_TIMESTEP *CONTROL_SOLUTION *DATABASE_BINARY_D3PLOT *DEFINE_CURVE_TITLE *ELEMENT_SHELL *ICFD_BOUNDARY_CONJ_HEAT *ICFD_BOUNDARY_FREESLIP *ICFD_BOUNDARY_FSI *ICFD_BOUNDARY_PRESCRIBED_VEL *ICFD_BOUNDARY_PRESCRIBED_TEMP *ICFD_BOUNDARY_PRESCRIBED_PRE *ICFD_BOUNDARY_NONSLIP *ICFD_CONTROL_FSI *ICFD_CONTROL_TIME *ICFD_DATABASE_DRAG *ICFD_INITIAL *ICFD_MAT *ICFD_PART *ICFD_PART_VOL *ICFD_SECTION *INCLUDE *INITIAL_TEMPERATURE_SET *KEYWORD *MAT *MESH_BL *MESH_SURFACE_ELEMENT *MESH_SURFACE_NODE *MESH_VOLUME *NODE *PARAMETER *PART *SECTION *END
$-----------------------------------------------------------------------------
$
$ Example provided by 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 R8.0.0 (or higher) with double precision
$X
$X------------------------------------------------------------------------------
$# UNITS: Dimensionless.
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
ICFD Conjugate heat transfer
*INCLUDE
mesh.k
*INCLUDE
struc.k
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ PARAMETERS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PARAMETER
R T_end 40.0
R dt_plot 0.50
$
$--- Fluid
$
R v_inlet 1.0
Rrho_fluid 1.0
R mu_fluid 0.005
R dt_fluid 0.050
RTemp_init 10.000
RTemp_inle 20.000
R Temp_cyl 80.000
R HC_fluid 1000.
R TC_fluid 200.
$
$---+----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
*ICFD_CONTROL_FSI
$# owc
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
$# hc tc
&HC_fluid &tc_fluid
*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
10 1 1
$# spid1 spid2 spid3 spid4
1 2 3 4
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD BOUNDARY/INITIAL CONDITIONS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_BOUNDARY_PRESCRIBED_VEL
$# pid dof vad lcid
1 1 1 1
*ICFD_BOUNDARY_PRESCRIBED_VEL
$# pid dof vad lcid
1 2 1 2
*ICFD_BOUNDARY_PRESCRIBED_TEMP
$# pid lcid
1 3
*ICFD_BOUNDARY_PRESCRIBED_PRE
$# pid lcid
2 2
*ICFD_BOUNDARY_FREESLIP
$# pid
3
*ICFD_BOUNDARY_NONSLIP
$# pid
4
*ICFD_BOUNDARY_FSI
$# pid
4
*ICFD_BOUNDARY_CONJ_HEAT
$# pid
4
*ICFD_INITIAL
$# pid vx vy vz temp
0 &Temp_init
*DEFINE_CURVE_TITLE
Velocity inlet
$# lcid sidr sfa sfo offa offo dattyp
1
$# a1 o1
0 0
5 0
6 &v_inlet
10000.0 &v_inlet
*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
Temp inlet
$# lcid sidr sfa sfo offa offo dattyp
3 &Temp_inle
$# a1 o1
0.0 1.0
10000.0 1.0
*DEFINE_CURVE_TITLE
Temp cylinder
$# lcid sidr sfa sfo offa offo dattyp
4 &Temp_cyl
$# a1 o1
0.0 1.0
10000.0 1.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD MESH KEYWORDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MESH_VOLUME
$# volid
1
$# pid1 pid2 pid3 pid4
1 2 3 4
*MESH_BL
$# pid nelth
4 1
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ DATABASE (OUTPUT) $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_DATABASE_DRAG
$# pid
4
*DATABASE_BINARY_D3PLOT
$# dt
&dt_plot
*END
This LS-DYNA simulation shows a simple conjugate heat transfer problem. The default coupling method between the structure and the fluid is a robust monolithic strong coupling approach which allows for accurate solutions. In order to trigger the conjugate heat transfer solver, the FSI keywords must be defined as well, since the FSI algorithm is needed in order to detect which nodes belong to the fluid/solid interface. Of course, the corresponding solid thermal input deck must be set up.