Conversion Between Material and Spatial Frames
This section explains how the user inputs are converted between material and spatial frames. The conversion depends on the dimension of the variables (scalars, vectors, or tensors) and on the density order.
As described in the previous paragraph Material and Spatial Frames, lowercase letters are used to denote the spatial frame coordinates while uppercase letters denote the material frame coordinates. In the followings, a physical quantity A will be referred to as A(x, y, z) in the spatial frame and to as A(X, Y, Z) in the material frame.
Most of the equations solved by the heat transfer interfaces are written in the spatial frame. When an input is specified in the material frame, conversion is necessary to deduce A(x, y, z) from A(X, Y, Z).
Density, Heat Source, Heat Flux
In heat transfer, the following variables are relative scalars of weight one (also called scalar densities):
•
Mass density, ρ (SI unit: kg/m3),
•
Heat source, Q0 (SI unit: W/m3),
•
Production/absorption coefficient, qs (SI unit: W/(m3·K)),
•
Heat flux, q0 (SI unit: W/m2),
•
Heat transfer coefficient, h (SI unit: W/(m2·K)).
For all these variables, the conversion between the material and the spatial frame follows the relation:
This way, the integral of volumetric quantities over the domain, such as the mass density, is invariant between frames:
In these equalities, Ω0 and Ω denote the same domain but represented in the material frame or the spatial frame, respectively. As expected, the same mass is found by integrating ρ(X, Y, Z) over the domain in the material frame or by integrating ρ(x, y, z) over the domain in the spatial frame. The same invariance principle applies to quantities per unit area, in particular for heat flux and heat transfer coefficients:
Here, ∂Ω0 and ∂Ω are the boundaries of the same domain in material and spatial frames, respectively.
Velocity
The relationship between the velocity vectors in the material and spatial frames, u(X, Y, Z) and u(x, y, z), is
This is directly deduced from the differential relation of Equation 4-211.
Thermal Conductivity
Thermal conductivity, k, is a tensor density. The relationship between the value on the spatial frame and the material frame is:
With this relation, and recalling that
and that the unit normal vector n verify that n(X, Y, Z) and FTn(x, y, z) are collinear, the total conductive heat flux through a boundary, computed in both frames according to the integrals below, gives the same result:
Here, ∂Ω0 and ∂Ω are the boundaries of the same domain in material and spatial frames, respectively.
When the material frame differs from the spatial frame, the domain features representing solid domains give the option to specify the thermal conductivity either on the spatial or the material frame. See Heat Conduction, Solid for details about this option.
Thermal Conductivity of a Layer
The same transformations are applied to thermal conductivity but with different transformation matrices. The deformation gradient tensor depends on the layer type:
•
When the layer is resistive, the deformation gradient tensor Fxdim is equal to the deformation gradient tensor F defined in Equation 4-210.
•
When the layer is conductive, the deformation gradient tensor Ft is defined using tangential derivatives as follows:
where xTX corresponds to the tangential derivative x with respect to X, and so on.
•
When the layer is an extra dimension, the deformation gradient tensor Fxdim is defined as follows:
where xTX corresponds to the tangential derivative x with respect to X, and so on. The (nx, ny, nz) vector corresponds to the normal vector in the spatial frame, and the (nX, nY, nZ) vector corresponds to the normal vector in the material frame.
Time Derivative
Partial differential equations often involve time derivatives of physical quantities such as temperature or internal energy in heat transfer. The variations of such derivatives during an elementary time step are studied for the same elementary volume that could be subjected to spatial transformations. The material derivative, denoted d ⁄ dt, is the derivation operator used in such cases. The following relation defines the material derivative in the spatial frame.
The right-hand side of this relation shows a new term u ⋅ ∇ corresponding to convection in the case of fluids or convected quantity by translational motion in the case of a solid.
•
About Frames in the COMSOL Multiphysics Reference Manual.
•
Handling Frames in Heat Transfer