Thermal Block

Description

A parametric semi-discretized heat transfer problem with varying heat transfer coefficients, the parameters, on subdomains. This model is also called cookie baking problem, and can be viewed as a flattened 2-D version of the skyscraper problem from high performance computing.

Figure 1: The computational domain and boundaries.

Figure 2: A sample heat distribution at time 1.0 for parameter choice [100, 0.01, 0.001, 0.0001].

Figure 3: Sigma magnitude plot of the single parameter variant.

Modeling

Consider a parameter $\mu \in {[10^{-6},10^2]}^4\subset {\mathbb{ℝ}}^4$ and define the heat conductivity $\sigma (\xi ;\mu )$ as $1.0$ when $\xi \in {\mathrm{\Omega }}_0$ and $\sigma (\xi ;\mu )={\mu }_i$ when $\xi \in {\mathrm{\Omega }}_i$. The heat distribution is governed by the equation:

${\partial }_t\theta (t,\xi ;\mu )+\mathrm{\nabla }\cdot (-\sigma (\xi ;\mu )\mathrm{\nabla }\theta (t,\xi ;\mu ))=0,\text{ for }t\in (0,T),\text{ and }\xi \in \mathrm{\Omega },$

with a heat-inflow condition on the left (Neumann boundary)

$\sigma (\xi ;\mu )\mathrm{\nabla }\theta (t,\xi ;\mu )\cdot n(\xi )=u(t)\text{ for }t\in (0,T),\text{ and }\xi \in {\mathrm{\Gamma }}_{in},$

perfect isolation on the top and bottom (Neumann-zero boundary)

$\sigma (\xi ;\mu )\mathrm{\nabla }\theta (t,\xi ;\mu )\cdot n(\xi )=0\text{ for }t\in (0,T),\text{ and }\xi \in {\mathrm{\Gamma }}_N,$

and fixed temperature on the right (Dirichlet boundary)

$\theta (t,\xi ;\mu )=0\text{ for }t\in (0,T),\text{ and }\xi \in {\mathrm{\Gamma }}_D,$

and initial condition

$\theta (0,\xi ;\mu )=0\text{ for }\xi \in \mathrm{\Omega }.$

Discretization

For the discretization, FEniCS 2019.1 was used on a simplicial grid with first order elements. The mesh is generated from the domain specification using gmsh 3.0.6 with 'clscale' set to $0.1$. The Python-based source code for the discretization can be found at Zenodo.

Origin

This benchmark was developed for the MODRED 2019 proceedings.

Data

The benchmark includes the basic domain description as a gmsh input file, Python scripts for the matrix assembly, simulation in pyMOR and visualization as VTK, together with the matrices both as one combined file ABCE.mat or separate matrix market files for all matrices. The sources and the ABCE.mat are available for download at Zenodo.

Note that the heat transfer coefficients are designed as characteristic functions on the domains, such that the system is only well-posed when all entries in $\mu$ are positive.

Dimensions

System structure:

$\begin{array}{rl}E\dot{x}(t)& =(A_0+{\mu }_1{A}_1+{\mu }_2{A}_2+{\mu }_3{A}_3+{\mu }_4{A}_4)x(t)+Bu(t),\\ y(t)& =Cx(t)\end{array}$

System dimensions:

$E\in {\mathbb{ℝ}}^{N\times N}$, $A_0\in {\mathbb{ℝ}}^{N\times N}$, $A_1\in {\mathbb{ℝ}}^{N\times N}$, $A_2\in {\mathbb{ℝ}}^{N\times N}$, $A_3\in {\mathbb{ℝ}}^{N\times N}$, $A_4\in {\mathbb{ℝ}}^{N\times N}$, $B\in {\mathbb{ℝ}}^{N\times 1}$, $C\in {\mathbb{ℝ}}^{4\times N}$,

where $N=7488$ for the system matrices given in ABCE.mat.

Variants

Besides the full four parameter setup, the model can be used in variations with other numbers of independent parameters. The following two are recommended in the original work and have been investigated in the literature,,,.

Single parameter

The interpretation of the thermal block as the "cookie baking" problem with slight variation in the dough leads to an easy one parameter variant. Here the new single parameter $\hat{\mu }\in [10^{-6},10^2]$ is chosen such that $\mu =\hat{\mu }[0.2,0.4,0.6,0.8].$

Non-parametric

The system can be used as a standard LTI state-space model. It is suggested to use $\mu =\sqrt{10}[0.2,0.4,0.6,0.8]$.

Citation

To cite this benchmark, use the following references:

• For the benchmark itself and its data:

S. Rave and J. Saak, Thermal Block. MORwiki - Model Order Reduction Wiki, 2020. http://modelreduction.org/index.php/Thermal_Block

@MISC{morwiki_thermalblock,
  author =       {Rave, S. and Saak, J.},
  title =        {Thermal Block},
  howpublished = {{MORwiki} -- Model Order Reduction Wiki},
  url =          {http://modelreduction.org/index.php/Thermal_Block},
  year =         2020
}

• For the background on the benchmark:

S. Rave and J. Saak, An Instationary Thermal-Block Benchmark Model for Parametric Model Order Reduction. e-prints 2003.00846, arXiv, math.NA (2020).

@TECHREPORT{morRavS20,
  author =       {Rave, S. and Saak, J.},
  title =        {A Non-Stationary Thermal-Block Benchmark Model for Parametric
                  Model Order Reduction},
  institution =  {arXiv},
  type =         {e-print},
  number =       {2003.00846},
  note =         {math.NA},
  year =         2020,
  url =          {https://arxiv.org/abs/2003.00846}
}

Contact

User:Saak