31. THERMAL_SOURCE Namelist

31.1. Overview

The THERMAL_SOURCE namelist is used to define external volumetric heat sources (power per unit volume). This source is in addition to any other sources coming from other physics, such as a Joule heat source. Each instance of this namelist defines a source, and the final source is the sum of all such sources (subject to some limitations).

Two forms of sources \(q\) can be defined:

  1. \(q(t,x) = f(t,x,T)\chi_S(x)\), where \(f\) is a user-defined function and \(S\) is a subdomain corresponding to one or more user-specified mesh cell sets. Here \(\chi_S\) is the characteristic function on \(S: \chi_S = 1\) for \(x \in S\) and \(\chi_S = 0\) for \(x \notin S\). Sources of this form can be summed as long as the interiors of their subdomains do not intersect.

  2. \(q(t,x) = A(t)\Sigma_jq_j\chi_j(x)\), where \(A\) is a user-defined time-dependent prefactor, \(q_j\) is a constant source of mesh cell \(j\), and \(\chi_j\) is the characteristic function on cell \(j\). The collection of values {\(q_j\)} is read from a data file.

31.2. THERMAL_SOURCE Namelist Features

Required/Optional : Optional
Single/Multiple Instances: Multiple

31.3. Components

Equation

Description : A unique name used to identify a particular instance of this namelist.
Type : string (31 characters max)
Default : none

cell_set_ids

Description : A list of cell set IDs that define the subdomain where the source is applied.
Type : a list of up to 32 integers
Default : none
Valid values: any valid mesh cell set ID
Note : Different instances of this namelist must apply to disjoint subdomains; overlapping of source functions of this form is not supported. Exodus II mesh element blocks are interpreted by Truchas as cell sets having the same IDs.

source

Description : The constant value of the heat source. To specify a function, use source_func instead.
Physical dimension: \(E/T L^3\)
Type : real
Default : none

source_func

Description : The name of a FUNCTION namelist that defines the source function. That function is expected to be a function of (t,x,y,z,T).
Type : string
Default : none

data_file

Description : The path to the data file. It is expected to be a raw binary file consisting of a sequence of 8-byte floating point values, the number of which equals the number of mesh cells. The order of the values is assumed to correspond to the external ordering of the mesh cells. The file can be created with most any programming language. In Fortran use an unformatted stream access file.
Type : string
Default : none

prefactor

Description : The constant value of the prefactor A. For a function use prefactor_func.
Type : real
Default : none

prefactor_func

Description : The name of a FUNCTION namelist that defines the function that computes the value of the time dependent prefactor \(A(t)\).
Type : string
Default : none