4.1. Galerkin or Streamline Upwind Petrov-Galerkin (SUPG) Formulation

For current-technology elements, PLANE292 and SOLID278, you can choose between the Galerkin or the Streamline Upwind Petrov-Galerkin (SUPG) Formulation described in this section. Instead, only the Galerkin formulation is available for legacy elements, PLANE55 and SOLID70.

4.1.1. Galerkin Formulation

The stability of the Galerkin formulation is measured in terms of the non-dimensional criteria called the element Peclet number (Gresho 59):

(4–1)

where:

= magnitude of the velocity vector
= element length dimension along the velocity vector direction, which depends on mesh size
= density of the fluid (input as DENS on MP command)
= specific heat of the fluid (input as C on MP command)
      (For PLANE292 and SOLID278, you can use either the MP or the TB command to specify and .)
= equivalent thermal conductivity along the velocity vector direction

The magnitude of the velocity vector is calculated as follows:

(4–2)

where:

= fluid velocity (mass transport) in x direction

= fluid velocity (mass transport) in y direction

= fluid velocity (mass transport) in z direction (relevant only for solid elements)

To specify velocity components in your model, issue the BF command for SOLID278 and PLANE292 or the R command for PLANE55 or SOLID70 (see details in Table 4.1: Specifications for Including Mass Transport in an Analysis).

Elements PLANE55 or SOLID70 are restricted to Pe<1, and if this condition is not met, an error message is printed during element formulation. When this error occurs, you must reduce the element size in the direction of the velocity vector to reduce Pe to less than one and enable the solution or switch to one of the current technology elements, SOLID278 or PLANE292 as described in the next section.

4.1.2. Streamline Upwind Petrov-Galerkin (SUPG) Formulation

There is no restriction on Pe for the current-technology elements, PLANE292 and SOLID278. The Streamline Upwind Petrov-Galerkin (SUPG) formulation is available for these elements, and it smooths oscillations without increasing solution time compared to the Galerkin formulation. You can also add one of three additional discontinuity capturing (DC) terms to the governing equation when using the SUPG formulation to further smooth persistent oscillations if they arise in high Pe conditions:

DC1 applies isotropic consistent diffusion to the flow field and works well for simple unidirectional flow fields, but it is not able to stabilize more complex flows. For example, a rotational flow field will require either the DC2 or DC3 term.
DC2 applies isotropic diffusion to the flow field, and works well for more complex flows, for example rotational flow fields.
DC3 applies anisotropic diffusion as a consistent additional diffusion orthogonal to the flow field. Like DC2, it is useful for complex flow fields.

Since all of the DC terms make the problem nonlinear, their use to stabilize the solution requires an iterative solution and increases solution time.

When using the SUPG formulation, you can activate one of the three DC terms to smooth oscillations if they arise in your solution by setting real constants listed in the table below using the R and RMORE commands. To include a DC term to include in your analysis you must set DC1, DC2, or DC3 to a non-zero value, typically 1.0. You can also modify the thermal gradient threshold value (TGRADMAG), above which the selected DC term will be applied (default is 1.0e-6). If you don't specify any real constants, the default setting is to base the analysis on the SUPG formulation without including any DC terms.

Real Constants Used to Activate the SUPG formulation and one of the three DC terms[a]
No. [b] Name Default Description
2 SUPG 1.0Acts as a multiplier on the stabilizing term of the SUPG formulation and enables you to choose also a DC term to smooth oscillations.

To activate the Galerkin Formulation, use a small number(1.e-12), as an exact zero value reverts back to 1.0 internally and enables SUPG formulation.

3 TGRADMAG 1.0e-6Thermal gradient threshold value, above which the DC term that you select will be applied. If the thermal gradient is less than the value you specify for TGRADMAG, the DC term is not applied, and vice versa, if it is greater, the nonlinear stabilizing DC term will be applied.
5 DC1 [c] 0.0Any non-zero value (typically 1.0) selects the DC1 term to be added to your analysis and acts as a multiplier on this stabilizing term.
6 DC2 [c] 0.0Any non-zero value (typically 1.0) selects the DC2 term to be added to your analysis and acts as a multiplier on this stabilizing term.
7 DC3 [c] 0.0Any non-zero value (typically 1.0) selects the DC3 term to be added to your analysis and acts as a multiplier on this stabilizing term.

[a] Available for current-technology elements PLANE292 and SOLID278.

[b] Refers to the real constant values defined by issuing the R and RMORE commands.

[c] Activating more than one DC term at a time will produce an error message.

The SUPG formulation and DC terms enable a smooth solution for the simulation of heat transfer by mass transport in complex flows with discontinuities. This capability is seen in the example below.

Example 4.1: Simulating Heat Transfer in the Rotating Flow of a Brake Disk

The setup of an analysis simulating the dissipation of heat by the rotating flow of a spinning brake disk is illustrated in Figure 4.1: Simulated Mass Transport of Heat - Brake Disk with Rotating Flow (A:Mesh and Loads, B:Prescribed flow field). A transient heat source is generated by friction when the brake pad is in contact with the brake disk. This is modeled as a heat generation thermal load in the first loadstep. In the second loadstep, the heat generation is removed. The region for heat generation is indicated in the figure by the red patch in A. Since the analysis is a 2D model, the heat loss to the ambient environment is modeled as a convection boundary condition at the circumference of the disk as indicated by the red arrows in A. The rotational flow field created by the spinning disk, shown in B, is the prescribed velocity field input for the simulation of advection heat transfer.

Figure 4.1: Simulated Mass Transport of Heat - Brake Disk with Rotating Flow (A:Mesh and Loads, B:Prescribed flow field)

Simulated Mass Transport of Heat - Brake Disk with Rotating Flow (A:Mesh and Loads, B:Prescribed flow field)


The animation below shows the transient computed temperature field in the rotating disk.

Figure 4.2: Animated Results

Animated Results