SURF152


3D Thermal Surface Effect

Valid Products: Pro | Premium | Enterprise | PrepPost | Solver | AS add-on

SURF152 Element Description

SURF152 may be used for various load and surface effect applications. It may be overlaid onto an area face of any 3D thermal element. The element is applicable to 3D thermal analyses. Various loads and surface effects may exist simultaneously. See SURF152 in the Mechanical APDL Theory Reference for more details about this element.

Figure 152.1: SURF152 Geometry

SURF152 Geometry

SURF152 Input Data

The geometry, node locations, and the coordinate system for this element are shown in Figure 152.1: SURF152 Geometry. The element is defined by four to ten nodes and by the material properties. The nodes for this element must share the nodes of the underlying solid element. An extra node (away from the base element) may be used for convection or radiation effects. Two extra nodes (away from the base element) may be used to more accurately capture convection effects. A triangular element may be formed by defining duplicate K and L node numbers as described in Degenerated Shape Elements. The element x-axis is parallel to the I-J side of the element.

The mass, volume, and heat-generation calculations use the element thicknesses at node I, J, K, and L (real constants TKI, TKJ, TKK, and TKL, respectively). Thickness TKI defaults to 0.0, and thicknesses TKJ, TKK, and TKL default to TKI. The mass calculation uses the density (material property DENS).

See Element Loading for a description of element loads. Convections or heat fluxes may be input as surface loads on the element.

The convection surface conductivity matrix calculation uses the film coefficient (input on the SFE command with KVAL = 0 and CONV as the label). If the extra node is used, its temperature becomes the bulk temperature. If the extra node is not used, the CONV value input with KVAL = 2 becomes the bulk temperature. The convection surface heat flow vector calculation uses the bulk temperature. On a given face, either a heat flux or a convection may be specified, but not both simultaneously.

For the extra node option (KEYOPT(5) = 1), film effectiveness and free stream temperatures may also be input for convection surface loads (input on the SFE command with the CONV label and KVAL = 3 and 4, respectively). If film effectiveness is input, bulk temperature is ignored.

Setting KEYOPT(7) = 1 multiplies the evaluated film coefficient by the empirical term ITS - TBIn, where TS is the element surface temperature, TB is the fluid bulk temperature, and n is an empirical coefficient (real constant ENN).

Convections and heat fluxes are multiplied by an area to obtain the heat flows. KEYOPT(12) determines whether the current area or the original area is selected for the calculation.

If KEYOPT(5) = 1 and flow information is available from FLUID116 with KEYOPT(2) = 1, the bulk temperature may be adjusted to the adiabatic wall temperature using KEYOPT(6) = 1, real constants OMEG (rotational speed) and NRF (recovery factor), and the logic described in the Mechanical APDL Theory Reference. For this adjustment, the axis of rotation may be defined as the global Cartesian X, Y or Z coordinate axis (KEYOPT(3)). When using the OMEG real constant, you can specify either numerical values or table inputs. If specifying table inputs, enclose the table name in % signs (for example, %tabname%). Rotational speed (OMEG) can vary with time and location. Use the *DIM command to dimension the table and identify the variables. For more information and examples on using table inputs, see Array Parameters of the Ansys Parametric Design Language Guide, Applying Loads Using TABLE Type Array Parameters in the Basic Analysis Guide, and Doing a Thermal Analysis Using Tabular Boundary Conditions in the Thermal Analysis Guide, as well as the description of the *DIM command in the Command Reference.

A film coefficient specified by the SFE command may be modified by activating the user subroutine USERCV with the USRCAL command. USERCV may be used to modify the film coefficient of a surface element with or without an extra node. It may be used if the film coefficient is a function of temperature and/or location.

If the surface element has an extra node (KEYOPT(5) = 1 or 2), the bulk temperature and/or the film coefficient may be redefined in a general way by user-programmable routine USRSURF116. USRSURF116 may be used if the bulk temperature and/or the film coefficient is a function of fluid properties, velocity and/or wall temperature. If a bulk temperature is determined by USRSURF116, it overrides any value specified by SFE or according to KEYOPT(6). Also, if a film coefficient is determined by USRSURF116, it overrides any values specified by SFE or USRCAL, USERCV. USRSURF116 calculation are activated by modifying the USRSURF116 subroutine and creating a custom version of Mechanical APDL; there will be no change in functionality without modifying USRSURF116. For more information, see User-Programmable Features (UPFs) in the Advanced Analysis Guide.

Heat generation rates are input on a per unit volume basis and may be input as an element body load at the nodes, using the BFE command. Element body loads are not applied to other elements connected at the same nodes. The node I heat generation HG(I) defaults to zero. If all other heat generations are unspecified, they default to HG(I). If all corner node heat generations are specified, each midside node heat generation defaults to the average heat generation of its adjacent corner nodes. For any other input heat generation pattern, unspecified heat generations default to zero. The heat generation load vector calculation uses the heat generation rate values.

As an alternative to using the BFE command, you can specify heat generation rates directly at the nodes using the BF command. For more information on body loads, see Body Loads in the Basic Analysis Guide.

SURF152 allows for radiation between the surface and the extra node. The emissivity of the surface (input as material property EMIS for the material number of the element) is used for the radiation surface conductivity matrix. The form factor FORMF and the Stefan-Boltzmann constant SBCONST are also used for the radiation surface conductivity matrix. The form factor can be either input as a real constant (defaults to 1) using KEYOPT(9) = 1 or it can be calculated automatically as a cosine effect using KEYOPT(9) = 2 or 3. For information on how the cosine effect depends on basic element orientation and the extra node location, see the Mechanical APDL Theory Reference. There is no distance effect included in the cosine effect. The Stefan-Boltzmann constant defaults to 0.119x10-10 (Btu/hr*in2* °R4)).

When KEYOPT(4) = 0, an edge with a removed midside node implies that the temperature varies linearly, rather than parabolically, along that edge. See Quadratic Elements (Midside Nodes) in the Modeling and Meshing Guide for more information about the use of midside nodes.

A summary of the element input is given in "SURF152 Input Summary". A general description of element input is given in Element Input.

SURF152 Input Summary

Nodes
I, J, K, L if KEYOPT (4) = 1 and KEYOPT (5) = 0
I, J, K, L, M if KEYOPT (4) = 1 and KEYOPT (5) = 1
I, J, K, L, M, N if KEYOPT (4) = 1 and KEYOPT (5) = 2
I, J, K, L, M, N, O, P if KEYOPT (4) = 0 and KEYOPT (5) = 0
I, J, K, L, M, N, O, P, Q if KEYOPT (4) = 0 and KEYOPT (5) = 1
I, J, K, L, M, N, O, P, Q, R if KEYOPT (4) = 0 and KEYOPT (5) = 2
Degrees of Freedom
KEYOPT(11) SettingDOF for all nodes except extra node(s)DOF for extra node(s) (KEYOPT(5) = 1 or 2)
0TEMPTEMP
1TTOPTEMP
2TBOTTEMP
Real Constants
FORMF, SBCONST, (Blank), OMEG, NRF, VABS,
TKI, TKJ, TKK, TKL, (Blank), (Blank),
ENN, GC, JC
See Table 152.1: SURF152 Real Constants for a description of the real constants
Material Properties
MP command: DENS (density), EMIS (emissivity) if KEYOPT(9) > 0)
Surface Loads
Convections -- 

face 1 (I-J-K-L) if KEYOPT(8) > 1

Heat Fluxes -- 

face 1 (I-J-K-L) if KEYOPT(8) = 1

Body Loads
Heat Generation -- 

HG(I), HG(J), HG(K), HG(L), and, if KEYOPT(4) = 0, HG(M), HG(N), HG(O), HG(P)

Special Features

Birth and death

KEYOPT(1)

Adiabatic wall temperature option:

0, 1, 2 -- 

See Adiabatic Wall Temperature as Bulk Temperature for information on these options.

KEYOPT(2)

Recovery factor (FR) option:

0, 1, or 2 -- 

See Adiabatic Wall Temperature as Bulk Temperature for information on these options.

KEYOPT(3)

Axis of symmetry:

0 -- 

OMEG used about global Cartesian X-axis

1 -- 

OMEG used about global Cartesian Y-axis

2 -- 

OMEG used about global Cartesian Z-axis

KEYOPT(4)

Midside nodes:

0 -- 

Has midside nodes

1 -- 

No midside nodes

KEYOPT(5)

Extra nodes:

0 -- 

No extra nodes. Use this option if the bulk temperature is known.

1 -- 

One extra node (optional if KEYOPT (8) > 1; required if KEYOPT (9) > 0). Valid for convection and radiation calculations. Use this option if the bulk temperature is unknown. The extra node gets the bulk temperature from a FLUID116 element.

2 -- 

Two extra nodes (optional if KEYOPT (8) > 1). Only valid for convection calculations. Use this option if the bulk temperature is unknown. The extra nodes get bulk temperatures from the two nodes of a FLUID116 element. This is generally more accurate than the one extra node option.

KEYOPT(6) (used only if KEYOPT(5) = 1 and KEYOPT(8) > 1)

Use of bulk temperature:

0 -- 

Extra node temperature used as bulk temperature

1 -- 

Adiabatic wall temperature used as bulk temperature

KEYOPT(7)

Empirical term:

0 -- 

Do not multiply film coefficient by empirical term.

1 -- 

Multiply film coefficient by empirical term |TS-TB|n.

KEYOPT(8)

Heat flux and convection loads:

0 -- 

Ignore heat flux and convection surface loads (if any)

1 -- 

Include heat flux, ignore convection

Use the following to include convection (ignore heat flux):

2 -- 

Evaluate film coefficient hf (if any) at average film temperature, (TS +TB)/2

3 -- 

Evaluate hf at element surface temperature, TS

4 -- 

Evaluate hf at fluid bulk temperature, TB

5 -- 

Evaluate hf at differential temperature, | TS - TB |

KEYOPT(9)

Radiation form factor calculation:

0 -- 

Do not include radiation

1 -- 

Use radiation with the form factor real constant

2 -- 

Use radiation with cosine effect calculated as an absolute value (ignore real constant)

3 -- 

Use radiation with cosine effect calculated as zero if negative (ignore real constant)

KEYOPT(11)

Label used for all nodal degrees of freedom of SURF152 when it is overlaid onto SHELL131 or SHELL132 (except for the extra node):

0 -- 

TEMP

1 -- 

TTOP

2 -- 

TBOT

The extra nodes, if requested with KEYOPT(5) > 0, are always TEMP.

KEYOPT(12)

Area to use for the heat-flow calculation:

0 -- 

Current area (default)

1 -- 

Original area

KEYOPT(13)

Film coefficient matrix key:

0 -- 

Program determines whether to use diagonal or consistent film coefficient matrix (default).

1 -- 

Use a diagonal film coefficient matrix.

2 -- 

Use a consistent matrix for the film coefficient.

Table 152.1: SURF152 Real Constants

No.NameDescription
1FORMFForm factor
2SBCONTStefan-Boltzmann constant
3(Blank)--
4OMEGAAngular velocity (KEYOPT(6) = 1)
5NRFRecovery factor
6VABSAbsolute value of fluid velocity (KEYOPT(1) = 0)
7TKIThickness at node I
8TKJThickness at node J (defaults to TKI)
9TKKThickness at node K (defaults to TKI)
10TKLThickness at node L (defaults to TKI)
11-12(Blank)--
13ENNEmpirical coefficient
14GCGravitational constant used for units consistency
15JCJoule constant used to convert work units to heat units

SURF152 Output Data

The solution output associated with the element is in two forms:

Convection heat flux is positive out of the element; applied heat flux is positive into the element. A general description of solution output is given in Solution Output. See the Basic Analysis Guide for ways to view results.

The Element Output Definitions table uses the following notation:

A colon (:) in the Name column indicates that the item can be accessed by the Component Name method (ETABLE, ESOL). The O column indicates the availability of the items in the file jobname.out. The R column indicates the availability of the items in the results file.

In either the O or R columns, “Y” indicates that the item is always available, a letter or number refers to a table footnote that describes when the item is conditionally available, and “-” indicates that the item is not available.

Table 152.2: SURF152 Element Output Definitions

NameDefinitionOR
ELElement NumberYY
SURFACE NODESNodes - I, J, K, LYY
EXTRA NODEExtra node (if present)YY
MATMaterial numberYY
AREASurface areaYY
VOLU:VolumeYY
XC, YC, ZCLocation where results are reportedY7
VN(X, Y, Z)Components of unit vector normal to center of element-Y
DENSITYDensity-1
MASSMass of element-1
HGENHeat generations HG(I), HG(J), HG(K), HG(L), HG(M), HG(N), HG(O), HG(P)2-
HEAT GEN. RATEHeat generation rate over entire element (HGTOT)22
HFLUXInput heat flux at nodes I, J, K, L3-
HEAT FLOW RATEInput heat flux heat flow rate over element surface area (HFCTOT)33
HFILMFilm coefficient at each face node44
TBULKBulk temperature at each face node or temperature of extra node44
TAVGAverage surface temperature44
TAWAdiabatic wall temperature55
RELVELRelative velocity55
SPHTFLSpecific heat of the fluid55
RECFACRecovery factor55
CONV. HEAT RATEConvection heat flow rate over element surface area (HFCTOT)44
CONV. HEAT RATE/AREAAverage convection heat flow rate per unit area4-
EMISSURAverage emissivity of surface (for element material number)66
EMISEXTEmissivity of extra node66
TEMPSURAverage temperature of surface66
TEMPEXTTemperature of extra node66
FORM FACTORAverage form factor of element66
RAD. HEAT RATERadiation heat flow rate over entire element (HRTOT)66
RAD. HEAT RATE/AREAAverage radiation heat flow rate per unit area6-

  1. If dens > 0

  2. If heat generation load is present

  3. If KEYOPT(8) = 1

  4. If KEYOPT(8) > 1

  5. If KEYOPT(6) = 1 and KEYOPT(8) > 1

  6. If KEYOPT(9) > 0

  7. Available only at centroid as a *GET item.

Table 152.3: SURF152 Item and Sequence Numbers lists output available through the ETABLE command using the Sequence Number method. See The General Postprocessor (POST1) in the Basic Analysis Guide and The Item and Sequence Number Table in this reference for more information. The following notation is used in Table 152.3: SURF152 Item and Sequence Numbers:

Name

output quantity as defined in Table 152.2: SURF152 Element Output Definitions

Item

predetermined Item label for ETABLE command

E

sequence number for single-valued or constant element data

I,J,K,L

sequence number for data at nodes I, J, K, L

Table 152.3: SURF152 Item and Sequence Numbers

Output Quantity NameETABLE and ESOL Command Input
ItemE
HGTOTSMISC1
HFCTOTSMISC2
HRTOTSMISC3
AREANMISC1
VNXNMISC2
VNYNMISC3
VNZNMISC4
HFILMNMISC5
TAVGNMISC6
TBULKNMISC7
TAWNMISC8
RELVELNMISC9
SPHTFLNMISC10
RECFACNMISC11
EMISSURNMISC12
EMISEXTNMISC13
TEMPSURNMISC14
TEMPEXTNMISC15
FORM FACTORNMISC16
DENSNMISC17
MASSNMISC18

SURF152 Assumptions and Restrictions

  • The element must not have a zero area.

  • If KEYOPT(9) > 0 (radiation is used):

    • element is nonlinear and requires an iterative solution

    • extra node must be present.

    • if KEYOPT(4) = 0, midside nodes may not be dropped.

  • If real constants TKI, TKJ, TKK, TKL (in-plane thicknesses), and TKPS (out-of-plane thickness) are defined for the element, the element volume is greater than zero. If BF, BFE, or BFUNIF commands are issued under this circumstance, heat generation loads are activated. However, the damping matrix is not activated even though the element volume is greater than zero.

SURF152 Product Restrictions

When used in the product(s) listed below, the stated product-specific restrictions apply to this element in addition to the general assumptions and restrictions given in the previous section.

Ansys Mechanical Pro  —  

  • Birth and death is not available.