- •Why CFD is Important for Modeling
- •How the CFD Module Helps Improve Your Modeling
- •Model Builder Options for Physics Feature Node Settings Windows
- •Where Do I Access the Documentation and Model Library?
- •Typographical Conventions
- •Quick Start Guide
- •Modeling Strategy
- •Geometrical Complexities
- •Material Properties
- •Defining the Physics
- •Meshing
- •The Choice of Solver and Solver Settings
- •Coupling to Other Physics Interfaces
- •Adding a Chemical Species Transport Interface
- •Equation
- •Discretization
- •Transport Feature
- •Migration in Electric Field
- •Reactions
- •Reactions
- •Initial Values
- •Initial Values
- •Boundary Conditions for the Transport of Concentrated Species Interface
- •Mass Fraction
- •Mass Fraction
- •Flux
- •Inflow
- •Inflow
- •No Flux
- •Outflow
- •Flux Discontinuity
- •Flux Discontinuity
- •Symmetry
- •Open Boundary
- •Physical Model
- •Transport Properties
- •Model Inputs
- •Fluid Properties
- •Diffusion
- •Migration in Electric Field
- •Diffusion
- •Model Inputs
- •Density
- •Diffusion
- •Porous Matrix Properties
- •Porous Matrix Properties
- •Initial Values
- •Initial Values
- •Domain Features for the Reacting Flow, Concentrated Species Interface
- •Boundary Conditions for the Reacting Flow, Concentrated Species Interface
- •Reacting Boundary
- •Inward Flux
- •Physical Model
- •Transport Properties
- •Fluid Properties
- •Migration in Electric Field
- •Porous Matrix Properties
- •Initial Values
- •Domain Features for the Reacting Flow, Diluted Species Interface
- •Boundary Conditions for the Reacting Flow, Diluted Species Interface
- •Pair and Point Conditions for the Reacting Flow, Diluted Species Interface
- •Multicomponent Mass Transport
- •Multicomponent Diffusion: Mixture-Average Approximation
- •Multispecies Diffusion: Fick’s Law Approximation
- •Multicomponent Thermal Diffusion
- •References for the Transport of Concentrated Species Interface
- •Domain Equations
- •Combined Boundary Conditions
- •Effective Mass Transport Parameters in Porous Media
- •Selecting the Right Interface
- •The Single-Phase Flow Interface Options
- •Laminar Flow
- •Coupling to Other Physics Interfaces
- •The Laminar Flow Interface
- •Discretization
- •The Creeping Flow Interface
- •Discretization
- •Fluid Properties
- •Fluid Properties
- •Mixing Length Limit
- •Volume Force
- •Volume Force
- •Initial Values
- •Initial Values
- •The Turbulent Flow, Spalart-Allmaras Interface
- •The Rotating Machinery, Laminar Flow Interface
- •Rotating Domain
- •Rotating Domain
- •Initial Values
- •Initial Values
- •Rotating Wall
- •Wall
- •Boundary Condition
- •Interior Wall
- •Boundary Condition
- •Inlet
- •Boundary Condition
- •Velocity
- •Pressure, No Viscous Stress
- •Normal Stress
- •Outlet
- •Boundary Condition
- •Pressure
- •Laminar Outflow
- •No Viscous Stress
- •Vacuum Pump
- •Symmetry
- •Open Boundary
- •Boundary Stress
- •Boundary Condition
- •Periodic Flow Condition
- •Flow Continuity
- •Pressure Point Constraint
- •Non-Newtonian Flow—The Power Law and the Carreau Model
- •Theory for the Pressure, No Viscous Stress Boundary Condition
- •Theory for the Laminar Inflow Condition
- •Theory for the Laminar Outflow Condition
- •Theory for the Slip Velocity Wall Boundary Condition
- •Theory for the Vacuum Pump Outlet Condition
- •Theory for the No Viscous Stress Condition
- •Theory for the Mass Flow Inlet Condition
- •Turbulence Modeling
- •Eddy Viscosity
- •Wall Functions
- •Initial Values
- •Wall Distance
- •Inlet Values for the Turbulence Length Scale and Intensity
- •Initial Values
- •The Spalart-Allmaras Turbulence Model
- •Inlet Values for the Turbulence Length Scale and Intensity
- •Pseudo Time Stepping for Turbulent Flow Models
- •References for the Single-Phase Flow, Turbulent Flow Interfaces
- •Selecting the Right Interface
- •Coupling to Other Physics Interfaces
- •Discretization
- •Fluid-Film Properties
- •Initial Values
- •Initial Values
- •Inlet
- •Outlet
- •Wall
- •Symmetry
- •Discretization
- •Initial Values
- •Initial Values
- •Fluid-Film Properties
- •Border
- •Inlet
- •Outlet
- •Conditions for Film Damping
- •The Reynolds Equation
- •Structural Loads
- •Gas Outflow Conditions
- •Rarefaction and Slip Effects
- •Geometry Orientations
- •References for the Thin-Film Flow Interfaces
- •Selecting the Right Interface
- •The Multiphase Flow Interface Options
- •The Relationship Between the Interfaces
- •Bubbly Flow
- •Coupling to Other Physics Interfaces
- •The Laminar Two-Phase Flow, Level Set Interface
- •Discretization
- •The Laminar Two-Phase Flow, Phase Field Interface
- •Domain Level Settings for the Level Set and Phase Field Interfaces
- •Fluid Properties
- •Mixing Length Limit
- •Initial Values
- •Initial Values
- •Volume Force
- •Volume Force
- •Gravity
- •Boundary Conditions for the Level Set and Phase Field Interfaces
- •Wall
- •Boundary Condition
- •Initial Interface
- •The Turbulent Flow, Two-Phase Flow, Level Set Interface
- •The Turbulent Two-Phase Flow, Phase Field Interface
- •Wall Distance Interface and the Distance Equation
- •Level Set and Phase Field Equations
- •Conservative and Non-Conservative Formulations
- •Phase Initialization
- •Numerical Stabilization
- •References for the Level Set and Phase Field Interfaces
- •Stabilization
- •Discretization
- •Level Set Model
- •Initial Values
- •Initial Values
- •Boundary Conditions for the Level Set Function
- •Inlet
- •Initial Interface
- •No Flow
- •Outlet
- •Symmetry
- •Discretization
- •Initial Values
- •Initial Values
- •Phase Field Model
- •Boundary Conditions for the Phase Field Function
- •Initial Interface
- •Inlet
- •Wetted Wall
- •Wetted Wall
- •Outlet
- •The Level Set Method
- •Conservative and Non-Conservative Form
- •Initializing the Level Set Function
- •Variables For Geometric Properties of the Interface
- •Reference for the Level Set Interface
- •About the Phase Field Method
- •The Equations for the Phase Field Method
- •Conservative and Non-Conservative Forms
- •Additional Sources of Free Energy
- •Variables and Expressions
- •Reference For the Phase Field Interface
- •The Laminar Bubbly Flow Interface
- •Reference Pressure
- •Discretization
- •The Turbulent Bubbly Flow Interface
- •Reference Pressure
- •Discretization
- •Fluid Properties
- •Slip Model
- •Initial Values
- •Initial Values
- •Volume Force
- •Volume Force
- •Gravity
- •Gravity
- •Mass Transfer
- •Mass Transfer
- •Boundary Conditions for the Bubbly Flow Interfaces
- •Wall
- •Liquid Boundary Condition
- •Gas Boundary Condition
- •Inlet
- •Liquid Boundary Condition
- •Gas Boundary Condition
- •Outlet
- •Liquid Boundary Condition
- •Gas Boundary Condition
- •Symmetry
- •Gas Boundary Conditions Equations
- •The Mixture Model, Laminar Flow Interface
- •Stabilization
- •Discretization
- •The Mixture Model, Turbulent Flow Interface
- •Stabilization
- •Mixture Properties
- •Mass Transfer
- •Mass Transfer
- •Initial Values
- •Initial Values
- •Volume Force
- •Volume Force
- •Gravity
- •Gravity
- •Boundary Conditions for the Mixture Model Interfaces
- •Wall
- •Mixture Boundary Condition
- •Dispersed Phase Boundary Condition
- •Inlet
- •Mixture Boundary Condition
- •Dispersed Phase Boundary Condition
- •Outlet
- •Mixture Boundary Condition
- •Symmetry
- •The Bubbly Flow Equations
- •Turbulence Modeling in Bubbly Flow Applications
- •References for the Bubbly Flow Interfaces
- •The Mixture Model Equations
- •Dispersed Phase Boundary Conditions Equations
- •Turbulence Modeling in Mixture Models
- •Slip Velocity Models
- •References for the Mixture Model Interfaces
- •Dispersed Phase
- •Discretization
- •Domain Conditions for the Euler-Euler Model, Laminar Flow Interface
- •Phase Properties
- •Solid Viscosity Model
- •Drag Model
- •Solid Pressure Model
- •Initial Values
- •Boundary, Point, and Pair Conditions for the Euler-Euler Model, Laminar Flow Interface
- •Wall
- •Dispersed Phase Boundary Condition
- •Inlet
- •Two-Phase Inlet Type
- •Continuous Phase
- •Dispersed Phase
- •Outlet
- •Mixture Boundary Condition
- •The Euler-Euler Model Equations
- •References for the Euler-Euler Model, Laminar Flow Interface
- •Selecting the Right Interface
- •The Porous Media Flow Interface Options
- •Coupling to Other Physics Interfaces
- •Discretization
- •Fluid and Matrix Properties
- •Mass Source
- •Mass Source
- •Initial Values
- •Initial Values
- •Boundary Conditions for the Darcy’s Law Interface
- •Pressure
- •Pressure
- •Mass Flux
- •Mass Flux
- •Inflow Boundary
- •Inflow Boundary
- •Symmetry
- •No Flow
- •Discretization
- •Fluid and Matrix Properties
- •Volume Force
- •Volume Force
- •Forchheimer Drag
- •Forchheimer Drag
- •Initial Values
- •Initial Values
- •Mass Source
- •Boundary Conditions for the Brinkman Equations Interface
- •Discretization
- •Fluid Properties
- •Porous Matrix Properties
- •Porous Matrix Properties
- •Forchheimer Drag
- •Forchheimer Drag
- •Volume Force
- •Volume Force
- •Initial Values
- •Initial Values
- •Boundary Conditions for the Free and Porous Media Flow Interface
- •Microfluidic Wall Conditions
- •Boundary Condition
- •Discretization
- •Domain, Boundary, and Pair Conditions for the Two-Phase Darcy’s Law Interface
- •Fluid and Matrix Properties
- •Initial Values
- •Initial Values
- •No Flux
- •Pressure and Saturation
- •Pressure and Saturation
- •Mass Flux
- •Inflow Boundary
- •Inflow Boundary
- •Outflow
- •Pressure
- •Darcy’s Law—Equation Formulation
- •About the Brinkman Equations
- •Brinkman Equations Theory
- •References for the Brinkman Equations Interface
- •Reference for the Free and Porous Media Flow Interface
- •Darcy’s Law—Equation Formulation
- •The High Mach Number Flow, Laminar Flow Interface
- •Surface-to-Surface Radiation
- •Discretization
- •Initial Values
- •Initial Values
- •Shared Interface Features
- •Fluid
- •Dynamic Viscosity
- •Inlet
- •Outlet
- •Consistent Inlet and Outlet Conditions
- •Pseudo Time Stepping for High Mach Number Flow Models
- •References for the High Mach Number Flow Interfaces
- •Selecting the Right Interface
- •The Non-Isothermal Flow Interface Options
- •Coupling to Other Physics Interfaces
- •The Non-Isothermal Flow, Laminar Flow Interface
- •Discretization
- •The Conjugate Heat Transfer, Laminar Flow Interface
- •The Turbulent Flow, Spalart-Allmaras Interface
- •Fluid
- •Dynamic Viscosity
- •Wall
- •Boundary Condition
- •Initial Values
- •Pressure Work
- •Viscous Heating
- •Dynamic Viscosity
- •Turbulent Non-Isothermal Flow Theory
- •References for the Non-Isothermal Flow and Conjugate Heat Transfer Interfaces
- •Selecting the Right Interface
- •The Heat Transfer Interface Options
- •Conjugate Heat Transfer, Laminar Flow
- •Conjugate Heat Transfer, Turbulent Flow
- •Coupling to Other Physics Interfaces
- •Accessing the Heat Transfer Interfaces via the Model Wizard
- •Discretization
- •Heat Transfer in Solids
- •Translational Motion
- •Translational Motion
- •Pressure Work
- •Heat Transfer in Fluids
- •Viscous Heating
- •Dynamic Viscosity
- •Heat Source
- •Heat Source
- •Initial Values
- •Initial Values
- •Boundary Conditions for the Heat Transfer Interfaces
- •Temperature
- •Temperature
- •Thermal Insulation
- •Outflow
- •Symmetry
- •Heat Flux
- •Heat Flux
- •Inflow Heat Flux
- •Inflow Heat Flux
- •Open Boundary
- •Periodic Heat Condition
- •Surface-to-Ambient Radiation
- •Boundary Heat Source
- •Boundary Heat Source
- •Heat Continuity
- •Pair Thin Thermally Resistive Layer
- •Pair Thin Thermally Resistive Layer
- •Thin Thermally Resistive Layer
- •Thin Thermally Resistive Layer
- •Line Heat Source
- •Line Heat Source
- •Point Heat Source
- •Convective Cooling
- •Out-of-Plane Convective Cooling
- •Upside Heat Flux
- •Out-of-Plane Radiation
- •Upside Parameters
- •Out-of-Plane Heat Flux
- •Domain Selection
- •Upside Inward Heat Flux
- •Change Thickness
- •Change Thickness
- •Porous Matrix
- •Heat Transfer in Fluids
- •Thermal Dispersion
- •Dispersivities
- •Heat Source
- •Equation Formulation
- •Activating Out-of-Plane Heat Transfer and Thickness
|
+ |
|
w |
C1/ 2k |
|
|
+ |
= |
10 |
|||||
|
= --- |
-------- |
--- |
-------- |
---- |
-- |
-- |
|
|
------------ |
||||
|
w |
|
|
|
|
|
|
|
|
|
|
w1 |
|
Pr1/ 3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
+ |
|
= |
10 |
|
10 |
|
|
|
|
Pr = |
Cp |
||
|
|
|
---- |
--- |
- |
|
|
---------- |
||||||
|
w2 |
|
|
|
|
PrT |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
||||
= 15Pr |
2/ 3 |
– |
PrT |
1 |
+ ln |
|
1000 |
|
|
|
|
|||
|
--2--- |
--- |
|
Pr------ |
-- |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
T |
|
|
|
where in turn is the thermal conductivity, and is the von Karman constant equal to 0.41.
The computational result should be checked so that the distance between the computational fluid domain and the wall, w, is almost everywhere small compared to any geometrical quantity of interest. The distance w is available as a postprocessing variable on boundaries.
References for the Non-Isothermal Flow and Conjugate Heat Transfer Interfaces
1.D.C. Wilcox, “Turbulence Modeling for CFD,” 2nd ed., DCW Industries, 1998.
2.Jonas Larsson, “Numerical Simulation of Turbulent Flows for Turbine Blade Heat Transfer”, Doctoral Thesis for the Degree of Doctor of Philosophy, Chalmers University of Technology, Sweden, 1998.
3.R.L. Panton, “Incompressible Flow”, 2nd ed., John Wiley & Sons, Inc., 1996.
4.W.M. Kays, “Turbulent Prandtl Number — Where Are We?”, ASME Journal of Heat Transfer, 116, pp. 284–295, 1994.
5.B. Weigand, J.R. Ferguson, and M.E. Crawford, “An Extended Kays and Crawford Turbulent Prandtl Number Model,” International Journal of Heat and Mass Transfer, vol. 40, no. 17, pp. 4191–4196, 1997.
6.D. Lacasse, È. Turgeon, and D. Pelletier, “On the Judicious Use of the k— Model, Wall Functions and Adaptivity”, International Journal of Thermal Sciences, vol. 43, pp. 925–938, 2004.
T H E O R Y F O R T H E N O N - I S O T H E R M A L F L O W A N D C O N J U G A T E H E A T T R A N S F E R I N T E R F A C E S | 399
400 | C H A P T E R 1 2 : N O N - I S O T H E R M A L F L O W B R A N C H
13
H e a t T r a n s f e r B r a n c h
This chapter describes the physics interfaces in the Heat Transfer branch (
) in the Model Wizard. As with all other physical descriptions simulated by COMSOL Multiphysics, any description of heat transfer can be directly coupled to any other physical process. This is particularly relevant for systems based on fluid-flow, as well as mass transfer. The interfaces also allows you to account for heat sources and sinks, such as energy evolving from chemical reactions. The Mechanisms for Modeling Heat Transfer in the CFD Module helps you choose the best one to start with.
In this chapter:
•The Heat Transfer Interfaces and The Heat Transfer Interface
•Out-of-Plane Heat Transfer Features
•The Heat Transfer in Porous Media Interface
•Out-of-Plane Heat Transfer Theory
•Theory for the Heat Transfer in Porous Media Interface
•The Non-Isothermal Flow and Conjugate Heat Transfer, Laminar Flow Interfaces and The Non-Isothermal Flow and Conjugate Heat Transfer, Turbulent Flow Interfaces are described in the fluid flow chapter.
401
T h e M e c h a n i s m s f o r M o d e l i n g H e a t T r a n s f e r i n t h e C F D M o d u l e
Heat Transfer is an important part of simulating fluid flow. Fluid flow has such an important effect on heating and cooling operations that it is difficult to simulate such without a complete description of how the flow is applied. Furthermore, the physical and thermodynamic properties of fluids are highly-affected by temperature such that flow can vary greatly in systems with temperature gradients.
While, the standard COMSOL Multiphysics package includes physics interfaces for simulating heat transfer through conduction and convection, the CFD Module provides extra functionality for simulating heat transfer in turbulent flow through The Non-Isothermal Flow and Conjugate Heat Transfer, Turbulent Flow Interfaces. The module also provides extra functionality for simulating heat transfer in porous media through The Heat Transfer in Porous Media Interface.
In this section:
•Selecting the Right Interface
•The Heat Transfer Interface Options
•Coupling to Other Physics Interfaces
Selecting the Right Interface
The Heat Transfer branch (
) included with this module has a number of subbranches to use to describe momentum transport. One or more of them can be added from the Model Wizard; either by themselves or in collaboration with other physics interfaces such as mass transfer, or even chemical reactions. Being a multiphysics interface, it contains the physics for modeling fluid-flow, which can be laminar, turbulent, and Stokes flow, as well as heat transfer through convection and conduction.
Several of the interfaces vary only by one or two default settings (see Table 13-1), which are selected at the Model Wizard level, or from a check box or drop-down list on the individual interfaces. The interfaces in the Non-Isothermal Flow and Conjugate
402 | C H A P T E R 1 3 : H E A T T R A N S F E R B R A N C H
Heat Flow branches relate to each other based on selections made from these drop-down lists in the Settings window (Figure ).
The Conjugate Heat Transfer interface is, by default in the Heat Transfer branch, set up to model heat transfer throughout a fluid in collaboration with a solid where heat is transferred by conduction. If a liquid regime is chosen in the Default model list, then the interface is renamed Non-Isothermal Flow, which is the same interface as Conjugate Heat Transfer but with different default settings as in Table 13-1.
For the Non-Isothermal Flow and Conjugate Heat Transfer branches, all the interfaces have the same interface identifier (nitf), where the differences are based on the default settings in the Physical Model section on the Settings window, and by what is required to model that type of heat transfer (see Table 13-1). Also, for all the
T H E M E C H A N I S M S F O R M O D E L I N G H E A T TR A N S F E R I N T H E C F D M O D U L E | 403
interfaces, the Neglect initial term (Stokes flow) check box is not selected by default.
TABLE 13-1: THE NON-ISOTHERMAL FLOW PHYSICAL MODEL DEFAULT SETTINGS*
INTERFACE (NITF) |
TURBULENCE |
TURBULENCE |
HEAT |
DEFAULT |
|
MODEL TYPE |
MODEL |
TRANSPORT |
MODEL |
|
|
|
TURBULENCE |
|
|
|
|
MODEL |
|
Non-Isothermal Flow, |
None |
n/a |
n/a |
Fluid |
Laminar Flow |
|
|
|
|
|
|
|
|
|
Non-Isothermal Flow, |
RANS |
k- |
Kays-Crawford |
Fluid |
Turbulent Flow, k- |
|
|
|
|
|
|
|
|
|
Non-Isothermal Flow, |
RANS |
Low |
Kays-Crawford |
Fluid |
Turbulent Flow, Low Re |
|
Reynolds |
|
|
k- |
|
number k- |
|
|
|
|
|
|
|
Non-Isothermal Flow, |
RANS |
Spalart-Allma |
Kays-Crawford |
Fluid |
Turbulent Flow, Spalart |
|
ras |
|
|
Allmaras |
|
|
|
|
|
|
|
|
|
Non-Isothermal Flow, |
RANS |
k- |
Kays-Crawford |
Fluid |
Turbulent Flow, k- |
|
|
|
|
|
|
|
|
|
Conjugate Heat Transfer, |
None |
n/a |
n/a |
Heat transfer |
Laminar Flow |
|
|
|
in solids |
|
|
|
|
|
Conjugate Heat Transfer, |
RANS |
k- |
Kays-Crawford |
Heat transfer |
Turbulent Flow, k- |
|
|
|
in solids |
|
|
|
|
|
Conjugate Heat Transfer, |
RANS |
Low |
Kays-Crawford |
Heat transfer |
Turbulent Flow, Low Re |
|
Reynolds |
|
in solids |
k- |
|
number k- |
|
|
|
|
|
|
|
Conjugate Heat Transfer, |
RANS |
Spalart-Allma |
Kays-Crawford |
Heat transfer |
Turbulent Flow, Spalart |
|
ras |
|
in solids |
Allmaras |
|
|
|
|
|
|
|
|
|
Conjugate Heat Transfer, |
RANS |
k- |
Kays-Crawford |
Heat transfer |
Turbulent Flow, k- |
|
|
|
in solids |
|
|
|
|
|
*For all the interfaces, the Neglect Initial Term (Stokes Flow) check box is not selected by default.
404 | C H A P T E R 1 3 : H E A T T R A N S F E R B R A N C H
