SOLID239
3D 20-Node
Diffusion Solid
SOLID239 Element Description
SOLID239 is a 3D 20-node diffusion element. The element has a single degree of freedom, concentration (CONC), at each node. The 20-node elements have compatible concentration shapes and are well suited to model curved boundaries.
The 20-node diffusion element is applicable to a steady-state or transient diffusion analysis. The element can also account for transport effects from a velocity field. See SOLID239 in the SOLID239 - 3D 20-Node Diffusion Solid in the Theory Reference for more details about this element.
SOLID239 Input Data
The geometry, node locations, and the coordinate system for this element are shown in Figure 239.1: SOLID239 Geometry. The element is defined by 20 nodes and the material properties. A prism-shaped element may be formed by defining duplicate K, L, and S; A and B; and O, P, and W node numbers. A tetrahedral-shaped element and a pyramid-shaped element may also be formed as shown in Figure 239.1: SOLID239 Geometry.
Orthotropic material directions correspond to the element coordinate directions. The element coordinate system orientation is described in Coordinate Systems. Properties not input by default are as described in the Material Reference. When saturated concentration (MP,CSAT) is specified, the meaning of the CONC degree of freedom is normalized concentration. The actual concentration is evaluated at the element centroid and output as SMISC,1. When unspecified, the saturated concentration (MP,CSAT) defaults to 1.0 and the nodal degree of freedom CONC is actual concentration. See Diffusion in the Theory Reference for more information on the diffusion analysis.
Nodal loads are defined with the D (Lab
= CONC) and F (Lab = RATE) commands. The
temperature (for material property evaluation only) body loads may be input based on their value
at the element's nodes or as a single element value (BF,
BFE). In general, unspecified nodal values of temperatures default to the
uniform value specified with the BFUNIF or TUNIF commands.
The transport velocity components (in the global coordinate system) are input with the
BF (Lab = VELO) command.
The material properties (DXX, DYY, DZZ, CSAT) and loads (DFLUX, DGEN, VELO) can be input as numerical values or as functions of primary variables by using tabular input. The possible primary variables are time (TIME), x, y, and z location (X, Y, and Z) in local or global coordinates, temperature (TEMP), and the concentration (CONC) degree of freedom. For more information and examples of using tabular input, see Defining Linear Material Properties Using Tabular Input in the Material Reference, Applying Loads Using Tabular Input in the Basic Analysis Guide, and the *DIM command in the Command Reference.
A summary of the element input is given in "SOLID239 Input Summary". A general description of element input is given in Element Input.
SOLID239 Input Summary
- Nodes
I, J, K, L, M, N, O, P, Q, R, S, T, U, V, W, X, Y, Z, A, B
- Degrees of Freedom
CONC
- Real Constants
None
- Material Properties
MP command: DXX, DYY, DZZ, CSAT
- Surface Loads
- Diffusion Flux (DFLUX) --
Face 1 (J-I-L-K), face 2 (I-J-N-M), face 3 (J-K-O-N), face 4 (K-L-P-O), face 5 (L-I-M-P), face 6 (M-N-O-P)
- Body Loads
- Temperature (TEMP) --
T(I), T(J), ..., T(Z), T(A), T(B)
- Diffusing Substance Generations (DGEN) --
DG(I), DG(J), ..., DG(A), DG(B)
- Transport Velocity --
VELO (input transport velocity components with the BF command)
- Special Features
- KEYOPT(1)
Diffusion damping matrix:
- 0 --
Consistent
- 1 --
Diagonalized
- KEYOPT(2)
Implementation of transport effects:
- 0 --
Matrix implementation; produces an unsymmetric matrix. In a linear analysis, a response to transport effects is achieved after one iteration.
- 1 --
Load vector implementation; produces a symmetric matrix. Requires at least two iterations to achieve a response to transport effects.
SOLID239 Output Data
The solution output associated with the element is in two forms:
Nodal concentrations included in the overall nodal solution
Additional element output as shown in Table 239.1: SOLID239 Element Output Definitions
The element output directions are parallel to the element coordinate system. 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 number refers to a table footnote that describes when the item is conditionally available, and “-” indicates that the item is not available. |
Table 239.1: SOLID239 Element Output Definitions
| Name | Definition | O | R |
|---|---|---|---|
| EL | Element Number | Y | Y |
| NODES | Nodes - I, J, K, L, M, N, O, P, Q, R, S, T, U, V, W, X, Y, Z, A, B | Y | Y |
| MAT | Material number | Y | Y |
| VOLU: | Volume | Y | Y |
| XC, YC, ZC | Location where results are reported | 1 | 2 |
| TEMP | Input temperatures T(I), T(J), ..., T(Z), T(A), T(B) | Y | Y |
| CG:X, Y, Z, SUM | Diffusion gradient components and vector magnitude | - | 1 |
| DF:X, Y, Z, SUM | Diffusion flux components and vector magnitude | - | 1 |
| CONC | Element concentration [3] | - | 1 |
The solution value is output only if calculated (based on input data). The element solution is at the centroid.
Available only at centroid as a *GET item.
Actual concentration obtained by multiplying the saturated concentration (MP,CSAT) and the normalized concentration evaluated at the element centroid. For more information, see Normalized Concentration Approach in the Theory Reference.
Table 239.2: SOLID239 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 239.2: SOLID239 Item and Sequence Numbers:
- Name
Output quantity as defined in Table 239.1: SOLID239 Element Output Definitions.
- Item
Predetermined Item label for ETABLE command.
- E
Sequence number for single-valued or constant element data.
SOLID239 Assumptions and Restrictions
The element must not have a zero volume or a zero length side. This occurs most frequently when the element is not numbered properly. Elements may be numbered either as shown in Figure 239.1: SOLID239 Geometry or in an opposite fashion.
An edge with a removed midside node implies that the concentration varies linearly, rather than parabolically, along that edge. See Quadratic Elements (Midside Nodes) in the Modeling and Meshing Guide for more information on the use of midside nodes.
Degeneration to the form of the pyramid should be used with caution. When degenerated, the element sizes should be small in order to minimize the field gradients. Pyramid elements are best used as filler elements in meshing transition zones.
CSAT cannot be temperature-dependent when used with this element. For a temperature-dependent CSAT, perform a thermal-diffusion analysis using the SOLID226 element with KEYOPT(1) = 100010.