Intermediate : Thermal Pipe
This example features a thermal flow in a pipe with a prescribed temperature or heat flux. The heat transfer coefficient is extracted and displayed. The velocity profile is laminar.
https://www.dynaexamples.com/icfd/intermediate-examples/thermalpipe
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Intermediate : Thermal Pipe
This example features a thermal flow in a pipe with a prescribed temperature or heat flux. The heat transfer coefficient is extracted and displayed. The velocity profile is laminar.
Fluid temperature fringes
*KEYWORD *TITLE *DATABASE_BINARY_D3PLOT *DEFINE_CURVE_TITLE *DEFINE_FUNCTION *ICFD_BOUNDARY_NONSLIP
*ICFD_BOUNDARY_PRESCRIBED_VEL
*ICFD_BOUNDARY_PRESCRIBED_PRE
*ICFD_CONTROL_MESH *ICFD_CONTROL_OUTPUT *ICFD_CONTROL_TIME
*ICFD_DATABASE_FLUX
*ICFD_DATABASE_TEMP
*ICFD_DATABASE_HTC *ICFD_MAT *ICFD_PART *ICFD_PART_VOL *ICFD_SECTION *INCLUDE
*MESH_BL
*MESH_BL_SYM *MESH_INTERF *MESH_SURFACE_ELEMENT *MESH_SURFACE_NODE *MESH_VOLUME *PARAMETER *END
$-----------------------------------------------------------------------------
$
$ Example provided by Iñaki and Rodrigo (LSTC)
$
$ E-Mail: info@dynamore.de
$ Web: http://www.dynamore.de
$
$ Copyright © 1992-2010 Livermore Software Technology Corporation
$ 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 R11. (or higher) with double precision
$X
$X------------------------------------------------------------------------------
$# UNITS: Dimensionless.
$X------------------------------------------------------------------------------
$X
*KEYWORD
*title
ICFD Thermal Pipe flow
*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.05
$
$--- Fluid
$
R v_inlet 1.0
R T_inlet 283.
R T_init 333.
R T_wall 333.
Rrho_fluid 1000.
R mu_fluid 1.
R hc_fluid 10.
R tc_fluid 2.
R dt_fluid 0.001
$---+----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_MESH
$# mgsf
1.2
*ICFD_CONTROL_OUTPUT
3
*DATABASE_BINARY_D3PLOT
0.1 0 0 0 0
0
*ICFD_CONTROL_MESH
1.2
$---+----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_fluid &tc_fluid
*ICFD_PART
1 1 1
*ICFD_PART
2 1 1
*ICFD_PART
3 1 1
*ICFD_PART
4 1 1
*ICFD_PART_VOL
5 1 1
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_NONSLIP
3
*ICFD_BOUNDARY_NONSLIP
4
*ICFD_BOUNDARY_PRESCRIBED_VEL
1 3 1 1 1.0 01.00000E28 0.0
*ICFD_BOUNDARY_PRESCRIBED_PRE
2 2 1.01.00000E28 0.0
*ICFD_BOUNDARY_PRESCRIBED_TEMP
1 3 1.01.00000E28 0.0
*ICFD_BOUNDARY_PRESCRIBED_TEMP
3 4 1.01.00000E28 0.0
*ICFD_BOUNDARY_PRESCRIBED_TEMP
4 4 1.01.00000E28 0.0
$*ICFD_BOUNDARY_FLUX_TEMP
$ 3 4 1.01.00000E28 0.0
$*ICFD_BOUNDARY_FLUX_TEMP
$ 4 4 1.01.00000E28 0.0
*DEFINE_CURVE_TITLE
inlet velocity
1 0 1.0 &v_inlet 0.0 0.0 0 0
0.0 1.0
10.0 1.0
*DEFINE_CURVE_TITLE
outlet pressure
2 0 1.0 1.0 0.0 0.0 0 0
0.0 0.0
10.0 0.0
*DEFINE_CURVE_TITLE
inlet temperature
3 0 1.0 &T_inlet 0.0 0.0 0 0
0.0 1.0
10.0 1.0
*DEFINE_CURVE_TITLE
wall temperature
4 0 1.0 &T_wall 0.0 0.0 0 0
0.0 333.0
10.0 333.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD MESH KEYWORDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MESH_BL
3 4 1.20 0.0001 3
*MESH_BL
4 4 1.20 0.0001 3
*MESH_BL_SYM
1 2 0 0 0 0 0 0
*MESH_VOLUME
5
1 2 3 4
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ DATABASE (OUTPUT) $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*DATABASE_BINARY_D3PLOT
&dt_plot
*ICFD_DATABASE_FLUX
1 &dt_plot
*ICFD_DATABASE_FLUX
2 &dt_plot
*ICFD_DATABASE_TEMP
3 &dt_plot
*ICFD_DATABASE_TEMP
4 &dt_plot
*ICFD_DATABASE_HTC
1 0 0.0 0 0 0 0 &dt_plot
*END
Various fringes
This example features a thermal flow in a pipe with a prescribed temperature or heat flux. The heat transfer coefficient is extracted and displayed. The velocity profile is laminar.
References :
[1] Frank P. Incropera, David P. Dewitt, Theodore L. Bergman, Adrienne S. Lavine, “Principles of Heat and Mass transfer”