Encyclopaedia Index

T Tabular output TABULATION TACT,the Cooling-Tower Program TALK MODE TECPLOT TEM1/2 TEMP0/1/2 TEMPERATURE TEMPR1 TEMPR2
TERMINATION OF SWEEPS AND ITERATIONS
TERMS TEST TEXT TFIRST TFRAC THEN THINPLT object THREE-DIMENSIONAL THRME1 THRME2 TICK
TIM/A TIME TIME-dependent boundary conditions TIME-OUT TINY TITLE TITLES of Input Files TLAST TMP1,TMP2,etc/2 TO TOLERANCE TRACE TRACK Transfer grid planes
Transfer Object Transfer Objects

Transience, time-step specification Transient flows, setting of Transport equations Truncation inaccuracies TSTSWP

Turbulence Turbulence energy generation Turbulent kinematic viscosity Turbulent Prandtl number
TURBULENCE MODELS TURBULENCE MODELS for two-phase flow
Turbulence_Model_Constants
TURMOD Tutorials

TWO-EQUATION turbulence models
Two-Fluid turbulence model
TWO-PHASE FLOWS
TWO-Layer KE-EP turbulence model.
TWO-SCALE KE-EP TURBULENCE MODEL.

TYPE settings for PATCHes TYPES of VR objects


T

is a PIL logical variable signifying "TRUE" (equivalent to the .TRUE. of Fortran)


Tabular output

To see how PHOENICS can print selected data into designated files of CSV format, click here.

Tabulation and plotting spot-values and residuals

ITABL)

TALK mode

(see Q1)


TEMP0

----- PIL real; default=0.0; group 9, ----

TEMP0....the negative of the absolute zero of temperature on the currently used temperature scale, i.e. the number which must be added to the current temperature to yield the absolute temperature.

It should be set to 273 when the Celsius scale is in use, and to 492 when the Fahrenheit scale is in use.

It is used for the calculation of enthalpy and specific heat when TEM1 and/or TEM2 is solved for, and also for radiative heat- transfer calculations.


TEMP1

----- PIL integer name; group 7 ----------

TEMP1.... is a name which may be used in GROUND coding to access the temperature of the first phase, when this is deduced from the enthalpy. It is not a PIL variable.

TEMP2

----- PIL integer name; group 7 ----------

TEMP1.... is a name which may be used in GROUND coding to access the temperature of the first phase, when this is deduced from the enthalpy. It is not a PIL variable.


TEMPerature

--------------------------------

Temperature can be treated in PHOENICS as either a derived property, or a solved-for variable.

In the former case, the enthalpy of the phase in question must be a solved-for variable; then the temperature is derived from it by some such relation as:

temperature = const1 * enthalpy + const2.

In the latter, temperature is the represented by the solved-for variables:
TEM1, for phase 1; and
TEM2, for phase 2.

See PHENC entries for: H1, TMP1, TMP1A, TMP1B, TMP1C, TEM1, CP1,

Further relevant information can be found in the entries on: PRPS and PROPS.


TEMPR1

Integer used in GXTEMPR to denote temperature of the first phase .


TEMPR2

Integer used in GXTEMPR to denote temperature of the second phase.


TERMS

----- Command; in group 8 ------------

The command to determine which terms are active in the balance equation for variables solved is:

TERMS(variable index,Y or N,Y or N,.. six times) (if no change desired, enter P for pass)
The six questions answered by the Y's and N's are:

  1. Built-in sources active?
  2. Convection active?
  3. Diffusion active?
  4. Transient term active?
  5. Variable belongs to first phase?
  6. Interphase transport active?

Note on 1.

This may be used to cut out the built-into-EARTH sources; it has no effect on GROUND-set sources. A list of the built-in sources is provided under the SOURCE entry.

Note on 4.

The transient term is automatically inactive when STEADY=T

Note on 5 and 6.

These entries are inactive when ONEPHS is T, in which case all variables are associated with the first phase, except U2, V2, W2, H2, R2 and RS which may be solved when ONEPHS=F only. The variables for which the default 5 entry is N, when ONEPHS=F, are: R2,RS,H2,C2,C4,C6,C8,C10,C12,C14,C16,C18,C20,C22,C24,C26,C28,C30,C32 and C34,

ie those with indices,
10,11,15,17,19,21,23,25,27,29,31,33,35,37,39,41,43,45,47 and 49.


Terms, activation of

(see TERMS command, Group 8)


TEST

------ PIL logical; default=F; group 25 --- -

TEST....activates calls to the SEARCH subroutines at certain locations in the algorithm-controlling subroutines. Set this T when SEARCHing.


TFIRST

---- PIL real; default= 0.0; group 2 --- -

TFIRST....initial value of time in the calculation of a time-dependent process.


TFRAC

----- PIL real array; default= 100*1.0; gr -

See PHENC


TFRAC

TFRAC is a PIL-language array for setting the time-step intervals. There are three ways of using TFRAC to set the time intervals namely:

  1. the direct method
  2. the method of pairs
  3. the use of GRDPWR (qv)

These methods are detailed below.

(a) The direct method.

When successive values of the TFRAC array are set to a series of monotonically increasing values, the time at the end of each time step will be equal to TLAST times the value in the array, minus TFIRST. Thus, if TLAST=10.0 (seconds), TFIRST=0 and LSTEP=4 (implying that four time intervals will be calculated), the following settings of TFRAC have the implications shown:-
TFRAC(1)= 0.1 implies that the first calculation is from 0.0 to 1.0 sec;
TFRAC(2)= 0.25 implies that the second calculation is from 1.0 to 2.5 sec;
TFRAC(3)= 0.75 implies that the third calculation is from 2.5 to 7.5 sec;
TFRAC(4)= 1.42 implies that the fourth calculation is from 7.5 to 14.2 sec.

Although TFRAC is properly interpreted as being a fraction, in that it multiplies TLAST, there is no objection to its taking values in excess of 1.0. provided each successive member of the array exceeds the one before it. If the last value of TFRAC is 1.0 then TLAST is the duration of the run and the initial time TFIRST is added to the time printed in the output.

(b) The method of pairs.

When TFRAC(1) is set to a negative number, PHOENICS-EARTH is informed that the method of pairs is in use. Thereafter, the minus sign is ignored.

In this method, the odd-numbered TFRAC's must be given whole- number values; but the even-numbered TFRAC's can take any positive values. The interpretation by EARTH of the following settings:


TFRAC(1)=-4.0 , TFRAC(2)=0.25 ,
TFRAC(3)= 1.0 , TFRAC(4)=0.137 ,
TFRAC(5)=42.0 , TFRAC(6)=1.2E-4 ,

will be that time will be divided into:- 4 intervals each equal to 0.25 times TLAST, followed by 1 interval equal to 0.137 times TLAST, followed by 42 intervals each equal to 0.00012 times TLAST.

(c) The command GRDPWR

GRDPWR(T,number of intervals, total time, exponent)

automatically sets TLAST, LSTEP and TFRAC. See GRDPWR for information.


THEN

---------- Advanced PIL command --- -

THEN is an element of the PIL IF construct. See the HELP entry on IF for further details.


THINPLT, a PHOENICS-VR object type

The THINPLT object represents a zero-thickness obstacle to flow, which may be porous. It has a nominal thickness and material for heat transfer through the plate. See the description in the PHOENICS_VR Reference Guide, TR326


Three-dimensional grids in PINTO

(see TR218)


THRME1

Integer used in GXTHRMX to denote 1st phase volumetric thermal-expansion coefficient.


THRME2

Integer used in GXTHRMX to denote 2nd phase volumetric thermal-expansion coefficient.


TICK command in AUTOPLOT

TI[CK]

Cause axis tick marks to be drawn on other side of axes. Revert to original by repeating TICK.

See also HELP on: AXES, BOX, FRAME


TIM

TIM is a Fortran real variable used in GROUND. It represents the current time.


TIMA

---------- PIL real; group 13 -----------

TIMA... is used to specify the wave amplitude of a periodic source term. See the help and encyclopaedia entries on ???, and GXTIM for further information.


Time arguments, setting of

(see GRDPWR command, Group 2)


Time-Out feature

PHOENICS permits users to set the maximum clock time of a simulation calculation, by setting the PIL variable MAXSEC to the maximum allowable number of seconds.

It was introduced for the convenience of users of CHAM's pay-by-use Remote Computing Service. However, it may of course be used in stand-alone installations also.

The coding which effects this is in the open-source file:
grex3.htm.

This employs the following call to an in-EARTH sub-routine:

CALL SECONDS(NUMSEC)

which returns the number of seconds of time which have elapsed since the EARTH run started to execute.

In the GREX3 sequence provided by CHAM, NUMSEC is simply compared with ISG20, in order to determine whether execution is to cease. Users who incorporate it in their own coding may well find other uses for it.


Time as a plotting direction in PHOTON

If NZ=1, and STEADY=F, it is possible to dump field values which can be plotted by PHOTON, with the time dimension replacing Z. In order to activate this feature, it is necessary only to set IDISPA > 0 in the Q1 file; then, because of the statement:

IF(.NOT.STEADY.AND.IDISPA.GT.0.AND.PNAM.EQ.' ') CALL GXPARA

in Group 19 Section 8 of GREX3.F, a call is made to subroutine GXPARA. IDISPA dictates the frequency (in terms of time steps) of dumping.

The file containing the dumped information is called PARPHI, if PHIDA=F in PREFIX, and otherwise PARADA.

The first time step at which dumping occurs is IDISPB and the last is IDISPC. However, if these are left as zero, dumping occurs for the whole time range.


Time direction, setting number of regions for

(see NREGT, Group 2)


Time step, setting index number for

(see LSTEP)


Time-marching computations

(see STEADY)


Time-step interval, print-out of

(see NPRMNT)


Time-step intervals, setting of

(see TFRAC)


Time-step specification

(see GROUP 2)


TINY

---------------------------------------

Real flag; value=1.E-20.

TINY....a small real number, defaulted to 1.E-20 which is used in EARTH to guard against division by zero. This value may be too small for machines with short word length, and would normally be set to an appropriate value at the time of installation.

When TINY.LE.1.E-20, underflow trapping occurs at various locations in the calculation. Setting TINY to a higher value therefore saves some time, and is acceptable for most problems on most machines.


TITLE

------ PIL character ---------------------

TITLE....is a special character variable which contains the character string generated by the TEXT command. Its default setting is:

NAME TITLE OF RUN HERE. MAX OF 40 CHARS.

This variable should be treated as a "read only" PIL variable. Only the TEXT command should be used to modify it.


TITLE

-------------------------------------- Photon Help ----

After pressing the button, type in the title of the current frame to be saved in the PHOTON save-file.


TLAST

----- PIL real; default= 1.0; group 2 --- -

TLAST....multiplier of TFRAC. If TLAST is set to GRND, EARTH visits group 2 of GROUND for a setting of the time-step size DT, rather than DT being determined from the settings of TFRAC.

See the Encyclopaedia for full details.


To

----------------------------------------- Photon Help ----

The last GEOMETRY element to be [DELETE]ed or turned ON/OFF.


TOLERANCE

'Tolerance' is a variable, having the dimensions of distance, which the VR Editor uses when creating a grid containing VR-objects which are specified as 'affecting the grid'.

All such objects do so by default.

Three 'tolerance' values can be set, one for each of the Cartesian-coordinate directions. By default, all three are set to 0.001 (m).

'Affecting the grid' means requiring the Satellite to create 'region' boundaries co-incident with the edges of the object's bounding box. The grid can then be specified within each region, either automatically or manually.

If the distance between two region boundaries, or a region boundary and the domain edge, is less than the tolerance in that direction, the second region boundary will not be created.

If it is desired that an object should not create region boundaries, the part of the Q1 file which pertains to it should contain:
> OBJ, GRID, NO
if the object is not to affect the grid in any direction, or
> OBJ, GRID, p,p,p
where 'p' is Y (for Yes) or N (for No) for the X, Y and Z co-ordinate directions respectively.

Tolerance values can be specified by the user either interactively or by editing the RSET(M,......) line in the Q1 file.

Judicious settings for the tolerance can help eliminate very thin 'spaghetti' cells formed when objects almost - but not quite - line up.


TRACE

----------- PIL logical, default F -------

When set equal to T, TRACE switches on a primitive tracing facility in the Q1 file which functions as an aid to debugging complicated Q1 files.

After the command TRACE=T has been read, each line of PIL is written either to the VDU or to the file lupvr (depending on the PHOENICS version) with the message ... 123.45 TRACING ... prefixed to it, before it is interpreted by the Satellite.

The numerical value prefixing the line shows the current machine clock time in seconds. The line-printing procedure is switched off when the command TRACE=F is encountered.


Track dir.

------------------------------- Photon Help ----

The direction in which the streamlines are being tracked; it can be either UP or DOWN.


TRANSFER, a PHOENICS-VR object type

The TRANSFER object transfers sources between calculations. See the description in the PHOENICS_VR Reference Guide, TR326, and click here for more information on transferring sources.


Transfer grid planes

(see GSET(T,...))


Transience, time-step specification

(see GROUP 2, GRDPWR, RSET) (See also TR326 Section 10, Space and Time Grids - Time)


Transient flows, setting of

(see STEADY=F) (See also TR326 Section 10, Space and Time Grids - Time)


Transport equations

The partial differential equations which govern fluid flow and heat transfer, for example those described here.


Truncation inaccuracies

A computer works to only a limited number of significant figures. It is therefore advisable to avoid producing calculated results of the following kinds:

These can be avoided by use of suitable datum levels for the enthalpy or pressure (see PRESS0) respectively, or by non- dimensionalisation.


TSTSWP

---- PIL integer; default=1; group 22 --- -

TSTSWP....dictates frequency (in terms of sweeps, or of hydrodynamic iterations for PARAB) with which sums of the residuals of the equations solved are written to the VDU, for the purpose of on-line monitoring of the convergence path of the solution. Note that the printed quantities are divided by the RESREFs.

If TSTSWP is set to a negative number and USEGRX=T, then EARTH will activate GXMONI. This will open a graphics window and display spot values and residuals both graphically and numerically.

TSTSWP=-1 dictates that changes are displayed every sweep or iteration,
TSTSWP=-2,-3,... means that changes are displayed every 2,3... sweeps or iterations.

An example of what appears on the screen can be seen by clicking here.


TURBULENCE

Flowing fluids are classified as 'turbulent' or 'laminar' according to whether they do or do not exhibit, superimposed upon an evident slowly or scarcely varying pattern of flow, seemingly-random transient deviations in lateral directions.

Present-day computers lack the power which they would need were they required to compute the small-scale details of turbulence; therefore 'models' of turbulence have been provided which permit then to compute some aspects of turbulence, sometimes to be approximately computed.

Effective-viscosity models

An essential supposition of many such models is that the time-averaged behaviour of a turbulent flow can be simulated by acting as though it were laminar, but with an increased 'effective' viscosity.

The problem then becomes one of calculating the increase ratio.

The widely-used so-called 'k-epsilon' model does this by supposing that the 'effective viscosity' depends on the 'turbulence energy', k, and its rate of dissipation, epsilon, for each of which soluble transport equations have been devised.

Population models

Other models seek to compute more than the time-average properties of the turbulence. They recognise that a sampling probe held at a fixed position in the smoke plume, for example, would at different times find pure air, dense smoke, and samples of varying degrees of smokiness; and they attempt to express this fact as a 'population distribution'.

The 'two-fluid' model and 'multi-fluid model', with which PHOENICS is supplied, are of this kind.

Which model should be used? And when?"

These are important questions for CFD users; yet there is no agreement among specialists as to the answers.

A broad-brush answer to which some would subscribe is:

Further information

Click here for a lecture on effective viscosity models, and here for the Encyclopaedia entry on turbulence models in general.

See also PHENC entry Turmod.


TURBULENCE-ENERGY GENERATION RATE, print-out of

This is effected when a variable named GEN1 is stored by means of STORE(GEN1) in the Q1 file, and the default print-out provisions are used. In addition, printout of the total rate of strain is effected when a variable named GENK is stored in the Q1 file.


Turbulent kinematic viscosity

See ENUT and

VISTRB


Turbulent Prandtl number, setting variable of

(see PRT)


TURMOD

---- Command; group 9 ---------------

TURMOD....command to activate a turbulence model. The syntax is:

TURMOD(model name) where model name is one of:-


        LVEL                KLMODL,          KEMODL,

        KEMODL-YAP,         KERNG,           KECHEN,

        KOMODL,             KOMODL-LOWRE,

        TSKEMO,             KECHEN-LOWRE,    KEMODL-LOWRE,

        KEMODL-LOWRE-YAP,   KEMODL-2L,       MIXLEN,

        MIXLEN-RICE,        SGSMOD,

        REYSTRS

WARNING! TURMOD's actions begin by wiping out all prior settings of enut, el1, etc. The TURMOD command must therefore be followed by any desired EL1 settings.

TURMOD is not used for the constant-effective-viscosity model, for which it suffices to set, in the Q1 file:

ENUT = whatever value is chosen.

>>> See the PHENC entry 'Turbulence Models in PHOENICS'.


TWO-EQUATION turbulence models

TWO-EQUATION turbulence models are those which solve two differential equations for the prediction of statistical properties of a turbulent flows, typically the "effective viscosity" and the "dissipation rate".

The most widely known model of this class is the so-called k-epsilon model, first proposed by Harlow and Nakayama.

There are many variants.

See PHENC entry: Turbulence models in PHOENICS


Two-phase flows, specification of

(see ONEPHS=F)


TWO-Layer KE-EP turbulence model.

See PHENC entry Two-layer KE-EP turbulence model


TWO-SCALE KE-EP TURBULENCE MODEL.

See PHENC entry Two-scale KE-EP turbulence model


TYPE settings for PATCHes

TYPE....the second argument of PATCH, used in Groups 11, 13 and 23.

The complete list of TYPEs, with their index equivalents and their significances, is provided below.
GroupTYPEvaluemeaning
11INIVAL26 initial value
11LINVLX27 linear with IX
11LINVLY28 linear with IY
11LINVLZ29 linear with IZ
---------------- --------
23PROFIL24 line-printer profile
23CONTUR25 line-printer contour
------------------------
13 multiplier
13CELL1 1.0
13EAST2 east cell area
13WEST3 west cell area
13NORTH4 north cell area
13SOUTH5 south cell area
13HIGH6 high cell area
13LOW7 low cell area
13VOLUME8 cell volume
13FREEE9 free east area
13FREEN10 free north area
13FREEH11 free high area
13FREEVL12 free volume
13PHASEM13 mass of phase
13RGRAD14 volume-fraction gradient
13OUTFLO15 1.0, or 0.0 if mass flow is inward
13EWALL17 area * gamma /distance
13WWALL18 area * gamma /distance
13NWALL19 area * gamma /distance
13SWALL20 area * gamma /distance
13HWALL21 area * gamma /distance
13LWALL22 area * gamma /distance
13INFLO23 1.0, or 0.0 if mass flow is outward

Notes:

  1. "Multiplier", for group 13 patches (i.e. sources and boundary conditions) signifies that by which CO * ( VAL - phi ) should be multiplied in order to create the required source;
  2. "gamma" signifies the near-wall value of viscosity divided by Prandtl number;
  3. "distance" signifies the distance from the near-wall cell centre to the wall.

wbs