The following is a list of the QTRAN data arrays to which user-supplied subroutines may require access. The variable dimensions for these arrays (e.g., J1, J2, ..., MAXT, etc.) are stored in the DIMS common block and may be found in the COMMONBLK file.

The following arrays are arranged in their order of common block declarations. The IA common blocks are Integer Array common blocks, while the RA common blocks are Real Arrays.

In the following, the arrays have been grouped according to usage.

IALIAS(N) contains the model’s node number for QTRAN’s internal node number N. The sign of IALIAS(N) determines whether the node’s temperature will be printed out in the QOUTDAT file, with a negative sign indicating that it will be printed out and a positive sign indicating that it will not be printed out.

This whole process of “aliased” node numbers is necessary if a code is to allow nonsequential node numbering in the model. QTRAN essentially looks at the node numbers that are used in the model and then condenses out any node numbers which are not used in the model that lie between 1 and the largest node number specified. This adds some unfortunate complexity to the coding, but it allows more flexibility in modeling without incurring unacceptable performance costs.

A useful QTRAN subroutine for determining a QTRAN internal node number for a given model node number is subroutine GETNUM. The calling sequence for GETNUM is as follows:

CALL GETNUM( YNODE, QNODE, IALIAS, NTEMPS, MATCH )

INTEGER YNODE, QNODE, IALIAS(MAXT), NTEMPS LOGICAL MATCH

and YNODE is the node number from the model, QNODE is the QTRAN internal node number that QTRAN has assigned to YNODE (QNODE is returned by GETNUM), IALIAS is the node alias array, NTEMPS comes from the COMMONBLK file and is the number of nodes which have been defined in the model, and MATCH is returned as .TRUE. or .FALSE., depending on whether GETNUM was able to match YNODE in the IALIAS array.

T1RC(K) and T2RC(K) contain QTRAN’s internal node numbers for conductive resistor K’s nodes 1 and 2. RCPROP(K) contains trains internal MPID number for thermal conductivity for resistor K, and RCAL(K) contains the area/length ratio for conductive resistor K. If solution option 1 (SOL = 1 in the QINDAT file) is being used, GVALCA(K) contains the conductance (kA/L) for resistor K.

RCN(N) contains the number of conductive resistors associated with QTRAN’s internal node N. RCNPNT is a pointer array used by QTRAN to determine which conductive resistors are associated with a particular node, and is used in conjunction with the RCN array. If RCN(1) is 10, the first 10 elements of the RCNPNT are conductive resistors which are associated with QTRAN internal node number 1. The conductive resistors associated with QTRAN internal node number 2 (if any) would be contained in RCNPNT(11) to RCNPNT(10+RCN(2)).

The QTRAN internal node numbers associated with convective resistor K are stored in T1RH(K), T2RH(K), and T3RH(K), and the convective configuration (CFIG value) for convective resistor K is stored in JTYPE(K). The GP values for convective resistor K are stored in GP(K,1) through GP(K,J6) while the QTRAN internal MPID numbers for convective resistor K are stored in JPROP(K,1) through JPROP(K,J2).

RHN(N) is the number of convective resistors associated with node N. The RHNPNT array interacts with the RHN array in a manner similar to the conductive resistors and RCN/RCNPNT interactions.

The QTRAN internal node numbers associated with advective resistor K are stored in T1RF(K) and T2RF(K). The QTRAN internal MPID numbers associated with specific heat and with variable mass flow rate for advective resistor K are stored in ICPFLO(K,1) and ICPFLO(K,2), respectively. The constant mass flow rate (or flow rate multiplier) associated with advective resistor K is stored in RMDOT(K).

The RFN and RFNPNT arrays are analogous to the conductive resistor RCN and RCNPNT arrays.

The QTRAN internal node numbers associated with hydraulic resistor K are stored in T1RP(K) and T2RP(K), and the fluid configuration (FCFIG value) for hydraulic resistor K is stored in PTYPE(K). The GP values for hydraulic resistor K are stored in PGP(K,1) through PGP(K,J46) while the QTRAN internal MPID numbers for hydraulic resistor K are stored in PPROP(K,1) through PPROP(K,J46).

RPN(N) is the number of hydraulic resistors associated with node N. The RPNPNT array interacts with the RPN array in a manner similar to the conductive resistors and RCN/RCNPNT interactions.

The hydraulic network introduces an additional network that is solved as a separate group of nodes with its own solution procedure separate from the thermal network solution. This imposes some restriction on the node definitions. All hydraulic nodes must be specified before the purely thermal nodes are specified, but there is no restriction on how the nodes are to be numbered. As with gaps in the thermal nodes require an aliasing array to relate user nodes numbers to internal node number, the same is needed for the fluid nodes, plus index to relate internal node numbers to the internal thermal node numbers is required.

PALIAS(N) contains the model node number as specified by the user to the internal node number N. PIALAS(N) contains the internal thermal node number corresponding to the internal fluid node number N. HIALAS(N) contains the internal fluid node number corresponding to the internal thermal node number N.

The QTRAN internal node numbers associated with gray radiative resistor K are stored in T1RR(K), T2RR(K), and T3RR(K). The resistor subtype for gray radiative resistor K is stored in IGTYPE(K), while the QTRAN internal MPID number associated with the resistor (for emissivity, transmissivity, or extinction coefficient) is stored in IGPROP(K). The view factor distance (if any) is stored in GDIST(K), the F value is stored in GSFACT(K), and the A value (if any) is stored in GAREA(K).

The RRN and RRNPNT arrays are analogous to the conductive resistor RCN and RCNPNT arrays.

The QTRAN internal node numbers associated with wavelength-dependent radiative resistor K are stored in T1RW(K), T2RW(K), and T3RW(K). The resistor subtype for wavelength-dependent radiative resistor K is stored in IWTYPE(K), while the QTRAN internal MPID number associated with the resistor (for emissivity, transmissivity, or extinction coefficient) is stored in IWPROP(K). The view factor distance (if any) is stored in WDIST(K), the F value is stored in SFACTR(K), and the A value (if any) is stored in WAREA(K). The shortest and longest wavelengths associated with wavelength-dependent radiative resistor K are stored in WSPROP(K) and WLPROP(K), respectively.

The RWN and RWNPNT arrays are analogous to the conductive resistor RCN and RCNPNT arrays.

The QTRAN internal node number for capacitor K is stored in T1CAPS(K), while the QTRAN internal specific heat, density, and phase change MPID numbers are stored in CPROP(K), CRHO(K), and PIDID(K), respectively. The volume associated with capacitor K is stored in CVOL(K).

The CN and CNPNT arrays are analogous to conductive resistor arrays RCN and RCNPNT.

These arrays are used to store material property data if the MPID defined was NOT a phase change data set (IEVAL not set to PH in the QINDAT file). For these MPID sets, QTRAN stores the MPID numbers into the MID array in the order they appear in the QINDAT or MATDAT files. For QTRAN internal nemophilas change MPID number K, the MPID which is defined for this material is stored in MID(K).

The MDATA1 and MDATA2 values for QTRAN internal MPID K are stored in COEFF(K,2) to COEFF(K,J4) for MDATA1 and in EXPO(K,2) to EXPO(K,J4) for the MDATA2 values. COEFF(K,1) is used to store the number of MDATA1/MDATA2 data pairs, while EXPO(K,1) is used to store a numeric key for the evaluation option IEVAL. Specifically:

These arrays are used to store material property data if the MPID defined WAS a phase change data set (IEVAL set to PH in the QINDAT file). For these MPID sets, QTRAN stores your MPID numbers into the PID array in the order they appear in the QINDAT or MATDAT files. For QTRAN internal phase change MPID number K, the MPID which is defined for this material is stored in PID(K).

The MDATA1 and MDATA2 values for QTRAN internal phase change MPID K are stored in PIDSET(K,1) through PIDSET(K,J23). For phase change set K, the number of MDATA1/MDATA2 data pairs is stored in PIDPAR(K).

MFID(K) contains the microfunction ID number (MFID) for the K'th microfunction which is defined in the QINDAT or MICRODAT files.

QCARD(K,1) contains an argument code for the microfunction which specifies whether the microfunction is to be a function of time, temperature, delta temperature, delta T4, or average temperature. Values are 0, 1, 2, 3, and 4, respectively.

QCARD(K,2) stores the microfunction option number as per Microfunction Library (Ch. 10).

QCARD(K,3) contains the number of parameters or data pairs entered for MDATA1 and MDATA2.

QCARD(K,4) contains a row pointer L into the P or TX/TY arrays where the MDATA values for the microfunction are stored.

For the K'th microfunction, L is stored in QCARD(K,4). The parameters MICDAT1-MICDAT4 are stored in P(L,1) through P(L,4) for those microfunctions which use parameters. For tabular microfunctions, the MICDAT1/MICDAT2 ordered pairs are stored in TX(L,2) -TX(L,J10) for MICDAT1 values and TY(L,2)-TY(L,J10) for MICDAT2 values. TX(L,1) contains the number of MDATA data pairs, and TY(L,1) is not used.

For the K'th QMACROfunction, IFLIST(K,1) contains the QTRAN internal node number to which the heat source is assigned.

IFLIST(K,2) contains the number of microfunctions used to build this QMACROfunction.

IFLIST(K,3) contains the T1 QTRAN internal node number.

IFLIST(K,4) contains the T2 QTRAN internal node number.

IFLIST(5)-IFLIST( 4+IFLIST(K,2) ) contains the QTRAN internal microfunction ID's which are used to build this QMACROfunction.

QMFACT(K) contains the scale factor for the K'th QMACROfunction.

QIP and QIPPNT are analogous to the conductive resistor RCN and RCNPNT arrays, respectively.

For the K'th MMACROfunction, IMLIST(K,1) contains the QTRAN internal node number to which the fluid mass flow rate is assigned.

IMLIST(K,2) contains the number of microfunctions used to build this QMACROfunction.

IMLIST(K,3) contains the P1 QTRAN hydraulic internal node number.

IMLIST(K,4) contains the P2 QTRAN hydraulic internal node number.

IMLIST(5)-IPLIST( 4+IMLIST(K,2) ) contains the QTRAN internal microfunction ID's which are used to build this MMACROfunction.

MMFACT(K) contains the scale factor for the K'th MMACROfunction.

MIP and MIPPNT are analogous to the conductive resistor RCN and RCNPNT arrays, respectively.

For the K'th TMACROfunction, ITLIST(K,1) contains the QTRAN internal node number to which the temperature source is assigned.

ITLIST(K,2) contains the number of microfunctions used to build this TMACROfunction.

ITLIST(K,3) contains the T1 QTRAN internal node number.

ITLIST(K,4) contains the T2 QTRAN internal node number.

ITLIST(5)-ITLIST( 4+ITLIST(K,2) ) contains the QTRAN internal microfunction IDs which are used to build this TMACROfunction.

TMFACT(K) contains the scale factor for the K'th TMACROfunction.

For the K'th MMACROfunction, IPLIST(K,1) contains the QTRAN internal node number to which the temperature source is assigned.

IPLIST(K,2) contains the number of microfunctions used to build this MMACROfunction.

IPLIST(K,3) contains the P1 QTRAN internal hydraulic node number.

IPLIST(K,4) contains the P2 QTRAN internal node hydraulic number.

IPLIST(5)-IPLIST( 4+IPLIST(K,2) ) contains the QTRAN internal microfunction ID's which are used to build this MMACROfunction.

PMFACT(K) contains the scale factor for the K'th MMACROfunction.

QBASE(N) is the constant nodal heat source value associated with QTRAN internal node number N. Once read from the QINDAT file in subroutine INPUT4, it is never updated by QTRAN. However, it can be modified by any of the User-Supplied subroutines.

MBASE(N) is the constant nodal net mass flow rate value associated with QTRAN internal node number N. Once read from the QINDAT file in subroutine INPUT4, it is never updated by QTRAN. However, it can be modified by any of the User-Supplied subroutines.

VELOCP is the average velocity at internal hydraulic node. PRHO is the average density at the internal hydraulic node. OPRESS is pressure from the previous interaction at the internal hydraulic node. PRESS is the current pressure at the internal hydraulic node. MDOTND is the net mass flow rate at the internal hydraulic node.

These arrays define parameters that apply to the hydraulic elements rather than results at the node locations. PRHOE is the average fluid density throughout an element. MDOTP is the mass flow rate of fluid flowing in a given hydraulic element. HYCCE is the hydraulic conductance throughout a given element. DIFHED is the difference in static head from one side of a hydraulic element to the other. QMDOTP is the volumetric flow rate of the fluid flowing through a given element.

TFIX(N) contains the fixed, not fixed, or TMACROfunction controlled code for node N. The values are 0: for not fixed, 1 for fixed, and 2 for TMACROfunction controlled.

The QINDAT CFIX information is stored in the FIXNUM, SETTIM, and FIXVAL arrays. For the K'th CFIX data line in the QINDAT file, FIXNUM(K) contains the QTRAN internal node number whose TFIX classification is being changed, SETTIM(K) is the time at which the change is to occur, and FIXVAL(K) is the new TFIX control code (see TFIX above) which is 0, 1, or 2.

PFIX(N) contains the fixed or not fixed code for the internal hydraulic node N. The hydraulic fixed flag operates on the fluid pressure. If a Hydraulic node is fixed, the pressure at that node is fixed.

OTEMPS(K) is the temperature of internal node number K at the beginning of a time step for transient runs. It is not used for steady-state runs.

TEMPS(K) is the temperature of internal node number K for steady-state runs. It is also the temperature at time t+dt (the end of the time step) for transient runs.

The temperatures stored in these arrays are in the temperature scale specified by the QINDAT parameter ICCALC.

TERROR(N,1) through TERROR(N,J27) contains the iterative error for QTRAN internal node N for the last J27 iterations. TERROR(N,1) is the most recent iterative error while TERROR(N,J27) is the least recent iterative error. This array is used in conjunction with the EPSIT2 parameter for determining which nodes should continue to be iterated upon.

QVECT(N) is the net nodal heat flow into QTRAN internal node number N, not including the heat flow to/from the capacitors associated with node N. This value is updated in subroutine CCAPAC.

ALPHA(N) is the explicit stable time step associated with QTRAN internal node number N. It is updated in subroutine CCAPAC.

DTMAXA is an array that is used to store the DTMAX adjustments. See QTRAN Run Control Parameters and Node Number Declarations (Ch. 8). The DTMAXA data pairs from the QINDAT file are stored in the DTMAXA array in the order in which they appear in the QINDAT file. DTMAXA(1,1) contains the number of DTMAXA data pairs in the QINDAT file. DTMAXA(1,2) contains a pointer to the next row of the DTMAXA array that will be used to adjust the DTMAX value. DTMAXA(K,1) contains the new DTMAX values, while DTMAXA(K,2) contains the implementation times.

DTMXHA is an array that contains the time step increment for hydraulic solutions. When time gets to the implementation time defined by DTMAXA(K,2), a new fluid time step DTMAXH is defined by DTMXHA(K).

PRINTA is an array that is used to store the TPRINT adjustments. See QTRAN Run Control Parameters and Node Number Declarations (Ch. 8). The PRINTA data pairs from the QINDAT file are stored in the PRINTA array in the order in which they appear in the QINDAT file. PRINTA(1,1) contains the number of PRINTA data pairs in the QINDAT file. PRINTA(1,2) contains a pointer to the next row of the PRINTA array that will be used to adjust the TPRINT value. PRINTA(K,1) contains the new TPRINT values, while PRINTA(K,2) contains the implementation times.

This array stores the indexes of the node that are defined in the plot output block.

These array are in common block RA43 and contain the node locations in the global coordinate system. The locations are packed in the same order that the nodes were defined in the node definition block.

Note that all three arrays are in the same named common block. This is the same as referencing the values as a doubly dimensioned array XYZCRD( MAXT, 3 ). Since this type of array is loaded by column in Fortran, the unit dimension for the array can not be used to reference the arrays. Either the MAXT value must be input for the specific case being executed, another subroutine called where the MAXT argument is included as one of the parameters passed and included in the declaration statement or only a single dimensioned array is specified and the actual index is calculated. For example, if the single index is used, the index would indicate the x location but MAXT and 2*MAXT would have to be added for the y and z locations respectively.

Another point is that the radial dimension for axisymmetric translation is always in the XRCRDV and the z-axis dimension is always in the YZCRDV arrays regardless of which plane was used to create the model. Thus if a model was created in the zx plane with the z-axis being defined as the radial direction all radial dimension would be in the XRCRDV array and the z-axis dimension which would have been along the x-axis would be in the YZCRDV array.

These arrays contain information related to calculations and printout of the relaxation parameters. The INDRLX, IPRLXC, and IRRLXC arrays are pointers that indicate the type of relaxation applied to a given node and counters which indicate if the relaxation parameter is to be updated. IRLXGR keeps the node IDs that have the greater iterative delta or system error for each type of node. RLXTBS and RLXTBT contain the relaxation parameters that are input through the QINDAT File. The first index is the node type, advection, radiation, etc. and the second index is the type of parameter: maximum value, damping factor, and multiplying factor. RELAXV is the relaxation parameter used by the specified node, EFACTB is the system error multiplier which is the measure of how quick the node is converging. RELAXM is the relaxation multiplier that is applied at each node. This is how under relaxation is applied. RERROR is the iterative delta associated with each node. The last three values are retained in the array with the first index representing the oldest value. RLXGRP is a summary of the maximum values encountered for each node group. This information is saved for print purposes.

The internal node ID contained in array TCPLND is lumped into internal node ID defined at the same position in array TCPLCN. Although the internal energy balanced is performed as though the TCPLND node is the same as the TCPLCN node, its individuality is retained and is updated after each converged iteration. NTCPL is the total number of temperature coupling and is defined in common block IB106. The arrays are originally sized with variable J52 which is defined in the DIMS Common Block.