Next: 2.8 Inquiring about D.S. Up: 2 Manipulation of solutions Prev: 2.6 Computing the stresses (D.S. TAE) Contents

2.7 Computing the flux and temperature (D.S. TAE)

Preprocessor FLUXXX constructs D.S. TAE of flux or temperatures obtained by combining the element flux and temperature arrays of D.S. TAE, constructed in conjunction with the element matrices and RHS, and D.S. B containing the temperatures or flux, obtained during the computation of the solution of the thermal problem under consideration.

The finite elements possible are those for which, in the temperature formulation, we compute a flux element array and those, in the flux formulation, for which we compute a temperature element array.

Module FLUXTH, called by FLUXXX, computes the flux or temperature in each element. For a two-dimensional problem (the visualization of the three-dimensional case has not been implemented), structure TAE resulting from FLUXTH is used to visualize the flux or temperatures (see TRSTXX in part II of this user guide).

Program FLUXTH has the following form:

 

      SUBROUTINE FLUXTH(M,XM,NFBE,NIBE,NFTAE,NITAE,NFNDL1,NINDL1,
     +                  NFTAES,NITAES,NTTAES,INDICA)
C   ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C   MODULE  FLUXTH : COMPUTE THE FLUX OR TEMPERATURES DURING THE
C   --------------   RESOLUTION OF A THERMAL PROBLEM FROM THE ELEMENT
C                    FLUX OR TEMPERATURES AND THE SOLUTION
C   ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

where

• M and XM represent the work array,
• NFBE, NIBE are the file number and level of D.S. B containing the temperatures (NFBE > 0: the file is sequential, NFBE = 0, the D.S. is in main memory),
• NFTAE, NITAE are the file number and level of D.S. TAE containing the element arrays (amongst which are the element flux). NFTAE > 0.
• NFNDL1, NINDL1, used if ND = 0, are the file number and level of the thermal D.S. NDL1 (ND = 0 indicates that the number of degrees of freedom per node is not constant),
• NFTAES, NITAES, NTTAES are the file number, level and the number of arrays to associate to the O.D.S. TAE containing the results element by element (NFTAES > 0),
• INDICA is set to 0 to compute the flux, and 1 to compute to temperatures.

The input files required are NFBE, NFTAE and, in certain cases, NFNDL1. The output file is NFTAES.

There is no input data required specifically for this module except when NTTAES is non-zero. In this case, the following input is required:

Loop I from 1 to NTTAES

• name of array I in 4 characters,
• type of array (1 = integer, 2 = real single precision, 4 = character, 5 = real double precision, ...),
• number of words in the array,
• the array variables,
• comments pertaining to the contents of array I in 18 words of 4 characters.

The output structure, TAE, is described in [MODULEF User Guide - 2].

REMTXX is used to compute the temperatures of mixed elements. REHPXX is used to compute the temperatures of primal hybrid elements (only TRIA HP12 is treated).