- •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
coupled or connected to the fluid-flow in some way. Process or component cooling are classic examples. The interface includes added functionality for calculating the added dispersion of heat transfer due to turbulence. This is represented by one of the Kays-Crawford or Extended Kays-Crawford Turbulence heat transport models, of by including your own turbulent Prandtl number.
There are four Conjugate Heat Transfer, Turbulent Flow interfaces, and each use the Reynolds-Averaged Navier-Stokes (RANS) equations, solving for the mean velocity field and pressure, along with the k-e model. See Table 13-1 and The Non-Isothermal Flow and Conjugate Heat Transfer, Turbulent Flow Interfaces for details.
Coupling to Other Physics Interfaces
Often, you may be simulating applications that couple heat transfer in turbulent flow to another type of phenomenon described in another physics interface. This can include chemical reactions and mass transport, as covered by the physics interfaces in the Chemical Species Transport branch.
Furthermore, the Chemical Reaction Engineering Module includes, not only support for setting up and simulating chemical reactions, but also for simulating reaction kinetics through the temperature-dependent Arrhenius Expression and Mass Action Law. This interface also includes support for including and calculating thermodynamic data as temperature-dependent expressions, for both reaction kinetics as well as fluid-flow.
In addition, the Heat Transfer Module also includes more detailed descriptions and tools for simulating energy transport, such as surface-to-surface and participating media radiation.
Location of Other Heat Transfer Documentation
Heat transfer through conduction and convection (both non-isothermal flow and conjugate heat transfer) in solid and free media is supported by physics interfaces shipped with the basic COMSOL Multiphysics license. The Joule Heating Interface is described in the COMSOL Multiphysics User’s Guide.
If you are using the Heat Transfer Module with enhanced fluid-flow interface features, due to the presence of the CFD Module, see also the Heat Transfer Module documentation for additional information.
To locate and search all the documentation, in COMSOL, select Help>Documentation from the main menu and either enter a search term or look under Heat Transfer Module in the documentation tree.
406 | 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
T h e H e a t T r a n s f e r I n t e r f a c e s
The following sections list all the physics interfaces and the features associated with them under the Heat Transfer branch. The descriptions follow a structured order as defined by the order in the branch. Because many of the interfaces are integrated with each other, some features described also cross reference to other interfaces. At the end of this section is a summary of the theory that goes towards deriving the physics interfaces under the Heat Transfer branch.
|
• Selecting the Right Interface |
See Also |
• The Heat Transfer Interface Options |
|
|
|
|
The Heat Transfer interfaces model heat transfer by conduction and convection. Surface-to-ambient radiation effects around edges and boundaries can also be included. The interfaces are suitable for modeling heat transfer in solids and fluids and in porous media. The interfaces are available in 1D, 2D, and 3D and for axisymmetric models with cylindrical coordinates in 1D and 2D. The default dependent variable is the temperature, T.
These key topics in this section:
•Accessing the Heat Transfer Interfaces via the Model Wizard
•Heat Transfer in Solids
•Heat Transfer in Fluids
•Boundary Conditions for the Heat Transfer Interfaces
Accessing the Heat Transfer Interfaces via the Model Wizard
There are Heat Transfer interfaces displayed in the Model Builder with the same name but with different icons and default models. After selecting a Heat Transfer interface in the Model Wizard, default settings are added under the main node. For example, if
Heat Transfer in Solids () is selected, a Heat Transfer node is added with a default Heat Transfer in Solids model. If Heat Transfer in Fluids () is selected, a Heat Transfer
T H E H E A T T R A N S F E R I N T E R F A C E S | 407
in Fluids model is added instead, but the parent nodes are both called Heat Transfer. Any interface based on the main Heat Transfer feature has the suffix ht. Select:
•Heat Transfer in Solids (ht) () to model mainly heat transfer in solid materials. A default Heat Transfer in Solids model is added, but all functionality for including fluid
domains is also available.
•Heat Transfer in Fluids (ht) () to model mainly heat transfer in fluid materials. A default Heat Transfer in Fluids model is added, but all functionality for including solid
domains is also available.
•Heat Transfer in Porous Media (ht) () to model mainly heat transfer in porous materials. The Porous Matrix and Heat Transfer in Fluids models are added, but all
functionality for including solid domains is also available.
Select another interface as required. Select:
•Joule Heating (jh) (), found under the Electromagnetic Heating subbranch (), to combine all features from the Electric Currents interface with the Heat Transfer
interface for modeling Joule heating (also called resistive heating or ohmic heating). See The Joule Heating Interface in the COMSOL Multiphysics User’s Guide.
•Conjugate Heat Transfer (), Laminar Flow (nitf) (), and Turbulent Flow (nitf)
() to use a predefined multiphysics coupling consisting of a Single-Phase Flow interface, using a compressible formulation, in combination with a Heat Transfer
interface. See Non-Isothermal Flow Branch.
The Non-Isothermal Flow (), Laminar Flow (nitf) (), and Turbulent Flow (nitf) () interfaces found under the Fluid Flow branch are identical to the Conjugate Heat Transfer interfaces. The only difference is that Fluid is selected as the default model. If Heat transfer in
solids is selected as the default model, the interface changes to a
Note
Conjugate Heat Transfer interface. To change the default model, select the Heat Transfer interface node and locate the Physical Model section in the Settings window.
408 | 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
T h e H e a t T r a n s f e r I n t e r f a c e
The Heat Transfer (ht) interface is available in many forms (see Accessing the Heat Transfer Interfaces via the Model Wizard) and each one has the equations, boundary conditions, and sources for modeling conductive and convective heat transfer, and solving for the temperature.
When this interface is added, default nodes are added to the Model Builder based on the selection made in the Model Wizard—Heat Transfer in Solids or Heat Transfer in Fluids,
Thermal Insulation (the default boundary condition), and Initial Values.
Depending on the version of the Heat Transfer interface selected, these
default nodes may be different.
Note
Right-click the Heat Transfer node to add other features that implement, for example, boundary conditions and sources.
I N T E R F A C E I D E N T I F I E R
The interface identifier is a text string that can be used to reference the respective physics interface if appropriate. Such situations could occur when coupling this interface to another physics interface, or when trying to identify and use variables defined by this physics interface, which is used to reach the fields and variables in expressions, for example. It can be changed to any unique string in the Identifier field.
The default identifier (for the first interface in the model) is ht.
D O M A I N S E L E C T I O N
The default setting is to include All domains in the model to define heat transfer and a temperature field. To choose specific domains, select Manual from the Selection list.
T H E H E A T T R A N S F E R I N T E R F A C E | 409