zone thermal command

Syntax

zone thermal keyword

Primary keywords:

active    analytical    cmodel    implicit    list    property    property-distribution    time-total    zone-based-temperature

Set parameters for a thermal analysis (only available if model configure thermal has been specified). The following keywords apply:

active b

Set thermal calculations on or off. The default is on after model configure thermal has been specified.

analytical keyword

isothermal-planes <x> <y> <z>

Define isothermal boundary layers. For example, the two planes defined by \(x\) = 0 and \(y\) = 0 are isothermal if the command zone thermal analytical isothermal x y is issued.

A plane cannot simultaneously be an isothermal boundary and a symmetry plane.

property keyword

conductivity f

Set the thermal conductivity used in the analytical solution.

diffusivity f

Set the thermal diffusivity used in the analytical solution.

list

Print to the screen information about the thermal analytical analysis.

solve time f

Solve the analytical solution to a time of f.

source keyword

components keyword <range>

Set general information about the heat source components for sources in the specified range.

component-number i

Specify the component of heat source.

fraction f

Specify the fraction of the initial source strength produced by each component.

decay f

Specify the decay constants of the components.

delete <range>

Delete all sources in the range.

grid keyword

Set up a grid of point sources. Source strengths are specified with the strength keyword and the start time is specified with the time-start keyword. These must be provided after the definition of the grid, or, can be given as a separate command with a range. If not provided, the default start time is 0. The source may be outside the FLAC3D grid.

point-1 v

Specify the top left corner of the grid.

point-2 v

Specify the top right corner of the grid.

point-3 v

Specify the bottom right corner of the grid.

n-12 i

Set the number of points across the grid from left to right.

n-23 i

Set the number of point across the grid from top to bottom.

group s1 <slot s2 >

The grid sources are assigned to group s1 when group is used. Group s1 is assigned to slot s2 if the optional slot keyword follows; it is assigned to the slot named Default if not. Use of the group logic is described in the i Group section.

strength f

Specify the strength of each point source in the grid.

time-start f

Set the start time of all point sources in the grid.

line keyword

Set up a line of point sources. Source strengths are specified with the strength keyword and the start time is specified with the time-start keyword. These must be provided after the definition of the grid, or, can be given as a separate command with a range. If not provided, the default start time is 0. The source may be outside the FLAC3D grid.

point-1 v

Specify the line start point.

point-2 v

Specify the line end point.

number i

Set the number of points along the line.

group s1 <slot s2 >

The line sources are assigned to group s1 when group is used. Group s1 is assigned to slot s2 if the optional slot keyword follows; it is assigned to the slot named Default if not. Use of the group logic is described in the i Group section.

strength f

Specify the strength of each point source in the grid.

time-start f

Set the start time of all point sources in the grid.

list <range>

Print to the screen information about the sources in the specified range.

strength f <range>

Specify the source strength for all sources in the range.

time-start f <range>

Specify the start time for all sources in the range.

symmetry-planes <x> <y> <z>

This command is used to define thermal symmetry planes. For example, the plane defined by \(y\) = 0 is a symmetry plane if the command block thermal symmetry-planes y is issued. Symmetry planes are used to represent adiabatic boundaries.

Note

A plane cannot simultaneously be an isothermal boundary and a symmetry plane.

tolerance f

Set a tolerance for thermal calculations. This tolerance can be thought of as a source radius. If temperatures are calculated at a point closer than f to a gridpoint, the distance from the gridpoint is taken as f. This is necessary because the analytical solution for the temperature due to a point source approaches infinity as the distance from the source approaches zero. (The default is f = 0.1.)

cmodel keyword <range>

Associate a thermal constitutive model with specified zones in the grid. It also allows querying of properties and states of constitutive models.

assign keyword

Assign a thermal model to all zones in the range. This keyword may only be given if the model is configured for thermal calculations (see model configure thermal and Mathematical Model Description).

advection-conduction

Assign advection-conduction model.

anisotropic

Assign thermal anisotropic model.

hydration

Assign thermal hydration model.

isotropic

Assign thermal isotropic model.

null

Assign thermally null model.

list

List the available thermal constitutive models and their properties.

implicit b <keyword>

Turn implicit mode on or off. By default implicit mode is off. Implicit mode requires a manually specified timestep. If no keyword is specified, the implicit logic uses the iterative Jacobi solver. If the Jacobi solver fails to converge to the solution, the logic switches to the preconditioned conjugate gradient solver (see below). The following keywords apply.

solver-jacobi

Use the iterative Jacobi solver based on the local matrix formulation (see Implicit Formulation). Note that specifying a timestep that is too large can cause the Jacobi solver to fail to converge. If this keyword is specified and the solver diverges, the system will exit with an error message. The benefit of using the Jacobi solver compared to other solvers is that it is usually faster for small implicit timesteps (<< the divergence timestep).

solver-direct <keyword>

Use the direct solver similar to matrix inversion (see Advanced implicit solvers). This highly efficient solver is able to handle small to medium size models (< 50-100K zones) and it is unconditionally stable for any timestep. Before this solver is called, the global right-hand-side (RHS) vector and global coefficient matrix (M) are formed for the whole model to incorporate contributions from all relevant gridpoints. Then the solver solves a system of linear equations for unknown global vector of temperature changes. The following keyword is available.

preserve <keyword>

Specify to preserve the global heat coefficients matrix and only update its diagonal containing timestep-dependent terms. In general, the global matrix M must be re-created each time the implicit solver is called as model parameters, such as geometry, physical properties, implicit timestep, etc., may change. However, if the model parameters do not change during the simulation but user wishes to specify a different implicit timestep (e.g., to gradually increase it every certain number of steps), only the diagonal terms of the global matrix must be updated, which is very efficient. To stop this procedure and form a new global heat coefficients matrix, the command zone thermal implicit has to be issued (this resets implicit solvers settings)

solver-pcg <keyword>

Use the preconditioned conjugate gradient iterative solver (see Advanced implicit solvers). This highly efficient solver is able to handle very large models and it is unconditionally stable for any timestep due to a special technique of matching desired solution accuracy to the number of iterations required to reach this accuracy. Before the solver is called, the global right-hand-side (RHS) vector and heat coefficients matrix (M) are created for the whole model to incorporate contributions from all relevant gridpoints. Then the solver analyses matrix properties (eigenvalues, condition number) and determines the required number of iterations, after which the iterative solution procedure starts. The following keywords are available.

tolerance

Specify error tolerance (solution accuracy) that will determine when to stop iterations. The default tolerance value is 1.e-6.

iterations-limit

Specify a limit on the number of iterations. By default, the number of iterations is determined based on the analysis of the global matrix properties, such as min/max eigenvalues and the the condition number. This analysis is usually done rather fast, and it assures that the system can be solved to reach the desired solution accuracy. Specifying a limit on the number of iterations skips the global matrix analysis, which may decrease overall solution time, but the risk of this is that the number of iterations specified may be insufficient to reach the desired accuracy. The pcg solver has general limit of 1000 on the number of iterations.

preserve

Specify to preserve the global heat coefficients matrix and only update its diagonal containing timestep-dependent terms. For a detailed description, see the same option above (for the solver-direct keyword).

list keyword

information

List general information about the zone-specific thermal calculation settings.

property keyword <range>

List the values of all thermal properties matching the property in the range.

property keyword <range>

Assign properties for the thermal constitutive models of all zones in the range. Any number of name value pairs may be assigned. Any constitutive model that has a matching property name will be assigned the value. The properties of the Advection-Conduction model are listed here, of the Anisotropic Heat Conduction model are listed here, of the Isotropic Heat Conduction model are listed here, and of the Hydration model are listed here.

property-distribution s akeyword <range>

Assigns a single property name to the thermal constitutive models of all zones in the range. The property value may be varied in space by using the following optional keyword values.

add

Add a to the property value. This keyword only applies to float type property values.

deviation-gaussian f

A Gaussian distribution is used to assign the value randomly, with a mean value of a and a standard deviation of f. This keyword only applies to float type property values.

deviation-uniform f

A uniform distribution is used to assign the value randomly, with a mean value of a and a standard deviation of f. This keyword only applies to float type property values.

gradient v <origin v >

Apply a gradient to the property value. This keyword only applies to float type property values.

multiply

Multiply the property value by a.

vary v

Apply a linear variation to the property value. This keyword only applies to float type property values.

time-total f

Specifies the accumulated zone thermal time, defined as the sum of all the timesteps over which zone thermal is active. Once set, thermal time will continue to accumulate with subsequent cycles.

zone-based-temperature b

Turns zone-based temperatures on and off. By default temperatures are calculated and stored at gridpoints, and the zone temperature is the average gridpoint value. When this is turned on, the zone temperature is stored in the zone and the gridpoint temperatures are ignored. Temperatures must be set manually using the zone.temp FISH intrinsic.