Approximation of Uncoupled Quasi-Static Thermoelasticity Solutions Based on Gaussians

A fast approximation method to three dimensional equations in quasi-static uncoupled thermoelasticity is proposed. We approximate the density via Gaussian approximating functions introduced in the method approximate approximations. In this way the action of the integral operators on such functions is presented in a simple analytical form. If the density has separated representation, the problem is reduced to the computation of one-dimensional integrals which admit efficient cubature procedures. The comparison of the numerical and exact solution shows that these formulas are accurate and provide the predicted approximation rate 2,4,6\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2,4,6$$\end{document} and 8\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$8$$\end{document}.


Introduction
The equations of thermoelasticity describe the elastic and the thermal behavior of elastic, heat conductive media, in particular the reciprocal actions between elastic stresses and temperature differences. We consider the classical thermoelastic system where the elastic part is the usual second-order one in the space variable. In the static uncoupled thermoelasticity, thermal effects on a body are restricted to strains due to a steady-state temperature distribution. Uncoupled quasi-static thermoelasticity can be employed when slowly varying thermal and mechanical loads are encountered and dissipative effected can be neglected. The equations are a coupling of the equations of elasticity and of the heat equation ([2, p.76], [3]) for (x, t) ∈ R 3 × [0, ∞), together with the corresponding initial and boundary conditions. The set of quantities μ, λ, γ, ρ, κ are positive and 3λ + 2μ > 0. We suppose that g : R 3 → R, F = (F 1 , F 2 , F 3 ) : R 3 × [0, ∞) → R 3 with g, F 1 (·, t), F 2 (·, t), F 3 (·, t) ∈ S (R 3 ). Here S (R 3 ) denotes the Schwartz space of smooth functions whose derivatives (including the function itself) decay at infinity faster than any power. The function T (x, t) is the temperature and the vector u(x, t) = (u 1 (x, t), u 2 (x, t), u 3 (x, t)) is the thermoelastic displacement.
In When the temperature field T is known, the displacement field u = (u 1 , u 2 , u 3 ) is obtained by solving (1.1) where the gradient of T is treated as a body force. The displacement field u = (u 1 , u 2 , u 3 ) with lim |x|→∞ |u(x, t)| = 0, t > 0, can be represented by means of the Kelvin fundamental matrix .
We write and u (2) 3 ) and u (2) 3 ) have the following integral representation by means of the Kelvin fundamental matrix Fast formulas of high order for the approximation of u (1) were obtained in [9]. The goal of this paper is to derive semi-analytic cubature formulas for (u (2) , T ) solutions to (1.7)-(1.2)-(1.3) of an arbitrary high-order which are fast and accurate by using the basis functions introduced in the theory approximate approximations ( [11,12]; see also [15] and the reference therein).
The approximate quasi-interpolant has the form where h and D are positive parameters and η is a smooth and rapidly decaying function which satisfies the moment conditions of order N If we define the Fourier transform of η as containing the remainder of the Taylor expansion of g. The functions x,ν (1.11) are rapidly oscillating multivariate trigonometric series and uniformly in x. Denoting Thus, at any point x we have where ∇ k g denotes the vector of partial derivatives {∂ α g} |α=k . The second term in the right hand side of (1.13) is called the saturation error.
Since η ∈ S (R 3 ) implies ε k (D) → 0 as D → ∞ a proper choice of the parameter D allows to make the terms ε k (D) as small as necessary, for example less than the machine precision. Therefore, the quasi-interpolant M h,D g can behave in numerical computations like a converging approximation process. Similar estimates hold in integral norms. Theorem 1.1. [15, p.42] Suppose that η ∈ S (R 3 ) satisfies the moment condition (1.9). Then for any where the constant c η does not depend on g, h and D.
New classes of cubature formulas for important integral operators of mathematical physics by using approximate approximations were studied in [14]. They are based on replacing the density of the integral operator by its quasi-interpolant where the generating function η is chosen such that the operator applied to it can be computed, analytically or at least efficiently. We choose as basis functions products of Gaussians and special polynomials. The use of the Gaussian functions for the numerical solution of the problems under consideration has the main advantage that the action of the integral operators on such functions may be presented in a simple analytical form.
By combining cubature formulas for volume potentials based on approximate approximations with the strategy of separated representations (cf., e.g. [1]), it is possible to derive a method for approximating volume potentials which is accurate and fast also in the multidimensional case and provides approximation formulas of high order. This procedure was applied successfully for the first time to the integration of the harmonic potential [5]. This approach was extended to the biharmonic [7], elastic and hydrodynamic [9] potentials, and to parabolic problems [6]. New approximation formulas for the solutions of nonstationary Stokes system were obtained in [8]. The static thermoelasticity was considered in [10]. Here we show that the fast method can be applied to uncoupled quasi-static thermoelasticity.
The outline of the paper is the following. In Sect. 2 we describe the fast formulas for the approximation of T obtained in [6]. In Sect. 3 we use the approximants obtained in Sect. 2 to construct approximation formulas for u (2) and give error estimates. In Sect. 4 we provide results of numerical experiments, illustrating that our formulas are accurate and provide the predicted approximation rates 2, 4, 6 and 8.

Approximation of T
Cubature formulas for (1.4) are derived by replacing the density g with the quasi-interpolant (1.8). Then provides an approximation formula for T (x, t).
The cubature error can be estimated by the following.  As basis functions in (1.8) we take the tensor products of univariate basis functions where H k are the Hermite polynomials Using the representation ([15, p.55]) By direct computation, the polynomials Q M satisfy Formula (2.3) easily follows.
Using formula (2.3), we can specify the high order approximation T for the generating function η 2M defined in (2.2). This is a semi-analytic cubature formula for (1. The approximation formulas (2.9) are very efficient if g has a separated representation, i.e. for a given accuracy ε it can be represented as the sum of products of vectors in dimension 1 Then T (hs, t) can be approximated by the sum of products of one-dimensional sums

Approximation of u (2)
In this section we propose formulas for the approximation of where T is given in (1.4).
Integrating by parts and using the relation where we set Since T (y, t) = (Pg)(y, t) we can also write where we denote by L the harmonic potential We use the representation (1.4) to get and we change the order of integration We use the representation ([15, p.128]) Now we replace g in (3.5) by the approximate quasi-interpolant (1.8) and we set In the next theorem we estimate the error of the cubature formula N h,D g. Theorem 3.1. Suppose that η satisfies the moment condition (1.9). Let 1 < p < 3, q = 3p/(3 − p), and let g ∈ W L p (R 3 ) with L > 3/p, L ≥ N . Then there exist two constants c and C such that, for any fixed t > 0 ||(N h,D g)(·, t)) − u (2) The constant c does not depend on g, h and D and C is independent of h. t))(x) and u (2) the norm ||∇u|| Lq is equivalent to the norm ||(−Δ) 1/2 u|| Lq ([13, p.458]) and L is the inverse of the Laplacian, we obtain [15, (6.14)], [16, (2.68)]) we see that for any t > 0 and p ≥ 1. In addition, the saturation error converges to zero with the order O(h N ). In [15, Paragraph 6.2.1] the inequality is proved with a constant c α depending on g and t. This shows that Hence, by Theorem 1.1 the assertion follows.
We assume the basis function (2.2). Keeping in mind (2.5) and (2.7) we have, for b > 0 Substituting in (3.6) we obtain, for k = 1, 2, 3, For example, for M = 1 we get the following formula suitable for fast computation

Implementation and Numerical Experiments
In this section we provide numerical experiments for the approximation of u (2) and T by means of (3.10) and (2.8), respectively.
The quadrature of the one-dimensional integrals which appears in (N

(M )
h,D g) k , k = 1, 2, 3, with certain quadrature weights ω p and nodes τ p leads to the approximation formulas at the point of a uniform grid {hs}

The approximation formulas (N (M )
h,D g) k , k = 1, 2, 3 are very efficient if g has a separated representation (2.10). Then an approximate value of u (2) k (hs, t) can be approximated using only one-dimensional operations as follows   Table 3. Absolute error and rate of convergence for u  with the one-dimensional convolutions We provide results of some experiments which show the accuracy and the convergence order of the method. We compute the solution of (1.7),(1.2),(1.3) with g(x) = e −|x| 2 . The exact solution of (1.2)-(1.3) is given by and, by using (3.2) and We assume κ = 1 and the parameters γ, λ, μ such that c γ,λ+2μ = 1. Following [18] the one-dimensional integrals in (3.10) are transformed to integrals over R with integrands decaying doubly exponentially by making the substitutions with certain positive constants α, β, and the computation is based on the classical trapezoidal rule. Then the tensor product structure of the integrands allows the efficient computation of N (M ) h,D g. In Table 1 we compare the exact values u (2) 1 in (4.2) and the approximates values (N 0.025,2 (e −|·| 2 )) 1 at some grid points x = (x, x, x) and t = 1. In Tables 2 and 3 we report on the absolute errors and approximate rates for the computation of u We have chosen α = 6, β = 5 in the transformation (4.3) and τ = 0.003 with 600 terms in the trapezoidal rule. The numerical results confirm the h 2M convergence of the approximating formula when M = 1, 2, 3, 4. For small h, the 8th-order formula has reached the machine precision.
In the next tables we report on numerical experiments for the approximation of T in (4.1) by means of (2.8).
In Table 4 we compare the values of the exact solution and the approximate solution at some points. The approximations in Table 4 have been computed on a uniform grid with step size h = 0.025 and N = 6.
In Tables 5 and 6 we show that formula (2.8) approximates the exact solution with the predicted approximate orders h 2M with M = 1, 2, 3, 4. For small h, the 6th-order and 8th-order formulas have reached the machine precision.
Funding Information Open access funding provided by Università degli Studi di Roma La Sapienza within the CRUI-CARE Agreement.
Open Access. This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.