Preprint
Article

Numerical Simulation for Transient Diffusion Convection Reaction Problems of a Class of Anisotropic FGMs

Altmetrics

Downloads

96

Views

25

Comments

0

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

12 April 2023

Posted:

12 April 2023

You are already at the latest version

Alerts
Abstract
In this paper a combined Laplace transform (LT) and boundary element method (BEM) is used to find numerical solutions to problems of anisotropic functionally graded media which are governed by the transient diffusion-convection-reaction equation. First, the variable coefficients governing equation is reduced to a constant coefficients equation. Then, the Laplace-transformed constant coefficients equation is transformed a boundary-only integral equation. Using a BEM, the numerical solutions in the frame of Laplace transform may then be obtained from this integral equation. Then the solutions are inversely transformed numerically using the Stehfest formula. Some problems considered are those of compressible or incompressible flow, and of media which are quadratically, exponentially and trigonometrically graded materials. The results obtained show that the analysis used to transform the variable coefficients equation into the constant coefficients equation is valid, and the mixed LT-BEM is easy to implement for obtaining the numerical solutions. The numerical solutions are verified by showing their accuracy and steady state. For symmetric problems, the symmetry of solutions is also justified. Moreover, the effect of the anisotropy and inhomogeneity of the material on the solutions are also shown, as to suggest that it is important to take the anisotropy and inhomogeneity into account when doing experimental studies.
Keywords: 
Subject: Computer Science and Mathematics  -   Computational Mathematics

1. Introduction

Over the last ten years, functionally graded materials (FGMs) have become a popular research topic, and many studies have been conducted on FGMs for various applications. FGMs are defined by authors as materials that are inhomogeneous and have properties that change spatially in a continuous manner, such as thermal conductivity, hardness, toughness, ductility, and corrosion resistance.
The diffusion convection reaction (DCR) equation has many applications in engineering, medicine, biology, and ecology. Several studies have been conducted to find numerical solutions to the DCR equation. Some of these studies include Fendoglu et al. [1] in 2018, Wang and Ang [2] in 2018, Sheu et al. [3] in 2000, Xu [4] in 2018, and AL-Bayati and Wrobel [5] in 2019, who considered the DCR equation with constant coefficients. Samec and Škerget [6] in 2004, Rocca et al. [7] in 2005, and AL-Bayati and Wrobel [8,9] in 2018 studied the DCR equation with variable velocity. Martinez et al. [10] used nonstandard finite difference schemes based on Green’s function formulations for reaction-diffusion-convection systems in 2013.
This paper is intended to extend the recently published works [11] on the steady DCR equation to the transient DCR equation for anisotropic functionally graded materials of the form
x i d i j x c x , t x j x i v i x c x , t k x c x , t = α x , t c x , t t
The continuously varying coefficients d i j , v i , k , α in (1) respectively represent the anisotropic diffusivity, velocity, decay reaction and change rate coefficients of the medium of interest. Therefore equation (1) is relevant for FGMs. Equation (1) provides a wider class of problems since it applies for anisotropic and inhomogeneous media but nonetheless covers the case of isotropic diffusion taking place when d 11 = d 22 , d 12 = 0 and also the case of homogeneous media which occurs when the coefficients d i j x , v i x , k x and α x , t are constant.

2. The initial boundary value problem

Referred to the Cartesian frame O x 1 x 2 we will consider initial boundary value problems governed by (1) where x = x 1 , x 2 . The coefficient d i j i , j = 1 , 2 is a real positive definite symmetrical matrix. Also, in (1) the summation convention for repeated indices applies, so that explicitly (1) can be written as
x 1 d 11 c x 1 + x 1 d 12 c x 2 + x 2 d 12 c x 1 + x 2 d 22 c x 2 x 1 v 1 c x 2 v 2 c k c = α c t
By knowing the coefficients d i j x , v i x , k x , α x , t we will seek solutions c x , t and its derivatives to (1) which are valid for time interval t 0 and in a region Ω in R 2 with boundary Ω which consists of a finite number of piecewise smooth curves. On Ω 1 the dependent variable c x , t is specified, and
P x , t = d i j x c x , t x i n j
is specified on Ω 2 where Ω = Ω 1 Ω 2 and n = n 1 , n 2 denotes the outward pointing normal to Ω . The initial condition is
c x , 0 = 0

3. The boundary integral equation

We restrict the coefficients d i j , v i , k , α to be of the form
d i j x = d ^ i j g x
v i x = v ^ i g x
k x = k ^ g x
α x , t = α ^ t g x
where g x is a differentiable function, d ^ i j , v ^ i , k ^ are constants and α ^ t is a function of time t. Substitution of (4)-(7) into (1) gives
d ^ i j x i g c x j v ^ i g c x i k ^ g c = α ^ g c t
Assume
c x , t = g 1 / 2 x ψ x , t
therefore use of (4) and (9) in (2) gives
P x , t = P g x ψ x , t + g 1 / 2 x P ψ x , t
where
P g x = d ^ i j g 1 / 2 x x j n i P ψ x , t = d ^ i j ψ x , t x j n i
Moreover, equation (8) can be written as
d ^ i j x i g g 1 / 2 ψ x j v ^ i g 1 / 2 ψ x i k ^ g 1 / 2 ψ = α ^ g g 1 / 2 ψ t
d ^ i j x i g g 1 / 2 ψ x j + ψ g 1 / 2 x j v ^ i g 1 / 2 ψ x i + ψ g 1 / 2 x i k ^ g 1 / 2 ψ = α ^ g g 1 / 2 ψ t
d ^ i j x i g 1 / 2 ψ x j + g ψ g 1 / 2 x j v ^ i g 1 / 2 ψ x i + ψ g 1 / 2 x i k ^ g 1 / 2 ψ = α ^ g 1 / 2 ψ t
Use of the identity
g 1 / 2 x i = g 1 g 1 / 2 x i
implies
d ^ i j x i g 1 / 2 ψ x j ψ g 1 / 2 x j v ^ i g 1 / 2 ψ x i + ψ g 1 / 2 x i k ^ g 1 / 2 ψ = α ^ g 1 / 2 ψ t
Rearranging and neglecting some zero terms gives
g 1 / 2 d ^ i j 2 ψ x i x j v ^ i ψ x i ψ d ^ i j 2 g 1 / 2 x i x j + v ^ i g 1 / 2 x i k ^ g 1 / 2 ψ = α ^ g 1 / 2 ψ t
So that if g satisfies
d ^ i j 2 g 1 / 2 x i x j + v ^ i g 1 / 2 x i λ g 1 / 2 = 0
where λ is a constant, then the transformation (9) brings the variable coefficients equation (1) into a constant coefficients equation
d ^ i j 2 ψ x i x j v ^ i ψ x i λ + k ^ ψ = α ^ ψ t
Taking the Laplace transform of (9), (10), (12) and applying the initial condition (3) we obtain
ψ * x , s = g 1 / 2 x c * x , s
P ψ * x , s = P * x , s + P g x ψ * x , s g 1 / 2 x
d ^ i j 2 ψ * x i x j v ^ i ψ * x i λ + k ^ + s α ^ * ψ * = 0
By using Gauss divergence theorem, equation (15) can be transformed into a boundary integral equation
η ζ ψ * ζ , s = Ω P ψ * x , s Φ x , ζ P x Φ x , ζ + Γ x , ζ ψ * x , s d S x
where
P v x = v ^ i n i x
For 2-D problems the fundamental solutions Φ ( x , ζ ) and Γ ( x , ζ ) for are given as
Φ ( x , ζ ) = ρ i 2 π D exp v ˙ R ˙ 2 D K 0 μ ˙ R ˙ Γ ( x , ζ ) = d ^ i j Φ x , ζ x j n i
where
μ ˙ = v ˙ / 2 D 2 + λ + k ^ + s α ^ * / D D = d ^ 11 + 2 d ^ 12 ρ r + d ^ 22 ρ r 2 + ρ i 2 / 2 R ˙ = x ˙ ζ ˙ x ˙ = x 1 + ρ r x 2 , ρ i x 2 ζ ˙ = ζ 1 + ρ r ζ 2 , ρ i ζ 2 v ˙ = v ^ 1 + ρ r v ^ 2 , ρ i v ^ 2 R ˙ = x 1 + ρ r x 2 ζ 1 ρ r ζ 2 2 + ρ i x 2 ρ i ζ 2 2 v ˙ = v ^ 1 + ρ r v ^ 2 2 + ρ i v ^ 2 2
where ρ r and ρ i are respectively the real and the positive imaginary parts of the complex root ρ of the quadratic equation d ^ 11 + 2 d ^ 12 ρ + d ^ 22 ρ 2 = 0 and K 0 is the modified Bessel function. Use of (13) and (14) in (16) yields
η g 1 / 2 c * = Ω g 1 / 2 Φ P * + P g P v g 1 / 2 Φ g 1 / 2 Γ c * d S
Equation (17) provides a boundary integral equation for determining the numerical solutions of c * and its derivatives at all points of Ω . The derivative solutions c * / ξ 1 and c * / ξ 2 can be determined using the following equations
c * ζ 1 = g 1 / 2 Ω g 1 / 2 Φ ζ 1 P * + P g P v g 1 / 2 Φ ζ 1 g 1 / 2 Γ ζ 1 c * d S c * g 1 / 2 ζ 1 c * ζ 2 = g 1 / 2 Ω g 1 / 2 Φ ζ 2 P * + P g P v g 1 / 2 Φ ζ 2 g 1 / 2 Γ ζ 2 c * d S c * g 1 / 2 ζ 2
Knowing the solutions c * x , s and its derivatives c * / x 1 and c * / x 2 from (17), the numerical Laplace transform inversion technique using the Stehfest formula is then employed to find the values of c x , t and its derivatives c / x 1 and c / x 2 . The Stehfest formula is
c x , t ln 2 t m = 1 N V m c * x , s m c x , t x 1 ln 2 t m = 1 N V m c * x , s m x 1 c x , t x 2 ln 2 t m = 1 N V m c * x , s m x 2
where
s m = ln 2 t m V m = 1 N 2 + m k = m + 1 2 min m , N 2 k N / 2 2 k ! N 2 k ! k ! k 1 ! m k ! 2 k m !
Possible multi-parameter solutions g 1 / 2 x to (11) are
g 1 / 2 x = constant , λ = 0 exp β 0 + β i x i , d ^ i j β i β j + v ^ i β i λ = 0
If the flow is incompressible, that is the divergence of the velocity is zero, then
v i x x i = 0
Therefore the governing equation (1) reduces to
x i d i j x c x , t x j v i x c x , t x i k x c x , t = α x c x , t t
Also, from (5) we obtain
v i x x i = 2 g 1 / 2 x v ^ i g 1 / 2 x x i = 0
so that
v ^ i g 1 / 2 x x i = 0
Therefore equation (11) reduces to
d ^ i j 2 g 1 / 2 x i x j λ g 1 / 2 = 0
Thus, for incompressible flow, possible multi-parameter functions g 1 / 2 x satisfying (20) are
g 1 / 2 x = β 0 + β i x i , λ = 0 cos β 0 + β i x i + sin β 0 + β i x i , d ^ i j β i β j + λ = 0 exp β 0 + β i x i , d ^ i j β i β j λ = 0

4. Numerical examples

We will examine multiple analytical and non-analytical test problems to demonstrate the accuracy and effectiveness of the mixed Laplace transform and boundary element method used in deriving the boundary integral equation (17). We will also analyze the efficiency, accuracy, and consistency of the combined LT-BEM method.
We assume each problem belongs to a system which is valid in given spatial and time domains and governed by equation (1). The system also is assumed to satisfy the initial condition (3) and some boundary conditions as mentioned in Section 2. The characteristics of the system which are represented by the coefficients d i j x , v i x , k x , α x , t in equation (1) are assumed to be of the form (4), (5), (6) and (7). They represents respectively the diffusivity or conductivity, the velocity of flow existing in the system, the reaction coefficient and the change rate of the unknown or dependent variable c x , t .
Standard BEM with constant elements is employed to obtain numerical results. For a simplicity, a unit square depicted in Figure 1 and Figure 11 will be taken as the geometrical domain for all problems. A number of 320 elements of equal length, namely 80 elements on each side of the unit square, are used.
A FORTRAN script is developed to compute the numerical solutions. A subroutine that evaluates the values of the coefficients V m , m = 1 , 2 , , N of the Stehfest formula in (18) for any number N is embedded in the script. Table 1 shows the values of V m for N = 6 , 8 , 10 , 12 resulted from the subroutine.

4.1. A test problem

The problems will consider three types of inhomogeneity functions g x , namely exponential function of the form (19) with compressible flow, and quadratic and trigonometric functions of the form (21) with incompressible flow. For all test problems we take coefficients d ^ i j and k ^
d ^ i j = 1 0.35 0.35 0.25 k ^ = 0.5
and a set of boundary conditions (see Figure 1)
P is given on side AB, BC, CD
c is given on side AD
For each problem, numerical solutions c and the derivatives c / x 1 and c / x 2 are sought at 19 × 19 interior points x 1 , x 2 = . 05 , . 1 , . 15 , , . 9 , . 95 × . 05 , . 1 , . 15 , , . 9 , . 95 and 9 time-steps t = . 0004 , π 8 , π 4 , 3 π 8 , π 2 , 5 π 8 , 3 π 4 , 7 π 8 , π . The value t = . 0004 is the approximating value of t = 0 as the singularity of the Stehfest formula. The individual relative error E I at each interior point and the aggregate relative error E A of the numerical solutions are computed using the formulas
E I = c n , i c a , i c a , i E A = i = 1 19 × 19 c n , i c a , i 2 i = 1 19 × 19 c a , i 2 1 2
where c n and c a are the numerical and analytical solutions c or the derivatives respectively.

4.1.1. Case 1

First, we suppose that the function g x is an exponential function
g 1 / 2 x = exp 1 + 0.15 x 1 0.25 x 2
that is, the material under consideration is an exponentially graded material. We choose
v ^ i = 1 , 1 α ^ t = 0.192625 t
so that the system has a compressible flow. In order for g x to satisfy (19) then λ = 0.088125 . The analytical solution c x , t for this problem is
c x , t = t exp 0.2 x 1 + 0.3 x 2 exp 1 + 0.15 x 1 0.25 x 2 .
Figure 2 (top row) shows the aggregate relative errors E A of the numerical solutions c obtained using N = 6 , 8 , 10 , 12 for the Stehfest formula (18). It indicates convergence of the Stehfest formula when the value of N changes from N = 6 to N = 10 . For this specific case (Case 1) the value of N is optimized at N = 10 . Increasing N to N = 12 does not give more accurate solutions. According to Hassanzadeh and Pooladi-Darvish [12] increasing N will increase the accuracy up to a point, and then the accuracy will decline due to round-off errors. Bottom row of Figure 2 depicts individual relative errors E I for the 19 × 19 interior points at time t = π / 2 (left) and t = π (right) with N = 10 as the optimized value of N. It indicates that the errors E I decreases as t changes from t = π / 2 to t = π . This result agrees with the result of the aggregate relative error E A in the top row of Figure 2.
For the derivative solution c / x 1 , Figure 3 (top row) shows that N = 6 is the optimized value of N for the aggregate relative errors E A . Bottom row of Figure 3 depicts individual relative errors E I with N = 6 . It indicates that the errors E I stay steady as t changes from t = π / 2 to t = π . This result agrees with the result of the aggregate relative error E A in the top row of Figure 3.
Whereas for the derivative solution c / x 2 , Figure 4 (top row) shows that N = 10 is the optimized value of N for the aggregate relative errors E A . Bottom row of Figure 4 depicts individual relative errors E I with N = 10 .

4.1.2. Case 2

Next, we choose an analytical solution
c x , t = exp 0.2 x 1 + 0.3 x 2 1 + 0.15 x 1 0.25 x 2 sin t
Suppose the function g x and the coefficients are
g 1 / 2 x = 1 + 0.15 x 1 0.25 x 2 v ^ i = 1 , 0.6 α ^ t = 0.031 t tan t
Therefore the considered system involves a quadratically graded material with an incompressible flow. From (21) we have the parameter λ = 0 .
Figure 5 (top row) indicates that N = 10 is the optimized value of N for the aggregate relative errors E A of the numerical solutions c. Increasing N to N = 12 gives worse solutions. Bottom row of Figure 5 depicts individual relative errors E I with N = 10 .
N = 10 is also the optimized value of N for the aggregate relative errors E A of the numerical solutions c / x 1 . This results is shown in Figure 6 (top row). Bottom row of Figure 6 depicts individual relative errors E I with N = 10 .
Whereas for the derivative solution c / x 2 , Figure 7 (top row) shows that N = 6 is the optimized value of N for the aggregate relative errors E A . Bottom row of Figure 7 depicts individual relative errors E I with N = 6 .

4.1.3. Case 3

Now, we consider a trigonometrically graded material with a grading function
g 1 / 2 x = cos 1 + 0.15 x 1 0.25 x 2
We choose
v ^ i = 1 , 0.6 , α ^ t = 0.003625 1 exp t
so that the system has an incompressible flow. From (21) we have λ = 0.011875 . The analytical solution c x , t for this problem is
c x , t = exp 0.2 x 1 + 0.3 x 2 1 exp t cos 1 + 0.15 x 1 0.25 x 2
Based on the results in Figure 8, Figure 9 and Figure 10 (top rows) we assume that N = 12 is the optimized value for the aggregate relative errors E A of the solutions c and the derivatives c / x 1 and c / x 2 . The corresponding individual relative errors E I are shown at the bottom row of each figure.

4.2. A problem without analytical solution

Further, we will show that the anisotropy and inhomogeneity of materials give impacts on the solutions. We will use d ^ i j , v ^ i , k ^ , g 1 / 2 x in Case 3 of Problem Section 4.1 for this problem, which are
d ^ i j = 1 0.35 0.35 0.25 v ^ i = 1 , 0.6 k ^ = 0.5 g 1 / 2 x = cos 1 + 0.15 x 1 0.25 x 2
We choose
α ^ t = 1
As we aim to show the impacts of the anisotropy and inhomogeneity of the material, we need to consider the case of homogeneous material and the case of isotropic material. We assume when the material is homogeneous then
g 1 / 2 x = 1
and for an isotropic material is under consideration then
d ^ i j = 1 0 0 1
The boundary conditions are (see Figure 11)
P = 0 on side AB
c = 0 on side BC
P = 0 on side CD
P = P t on side AD
where P t is associated with four cases, namely
Case 1: P t = 1
Case 2: P t = exp t
Case 3: P t = t
Case 4: P t = t / t + 0.01
Figure 12 shows for all cases when the material is isotropic and homogeneous the solutions c 0.5 , 0.3 , t and c 0.5 , 0.7 , t coincide. This is to be expected as the problem is geometically symmetric at x 2 = 0.5 when the material is isotropic and homogeneous. The results in Figure 12 also indicate that the anisotropy and inhomogeneity of the material give effects on the solutions. Moreover, as also expected, the variation of the solution with respect to t mimics the time function P t as the boundary condition on side AD.
Whereas, the results in Figure 13 show that the Case 1 of P ( t ) = 1 and Case 4 of P ( t ) = t / t + 0.01 have the same steady state solution. This is to be expected as P ( t ) = t / t + 0.01 will converge to 1 when t approaches infinity.

5. Conclusion

Several problems for a class of anisotropic FGMs (quadratically, exponentially and trigonometrically graded materials) have been solved using a combined BEM and Laplace transform. From the results of the considered problems in Section 4.1 and Section 4.2, we may conclude that the analysis of reduction to constant coefficients equation (in Section 3) for deriving the boundary-only integral equation (17) is valid, and the BEM and Stehfest formula is appropriate for solving such problems as defined in Section 2. Moreover, the results of the test problem in Section 4.1 show the accuracy of the method, whereas the results of the problem in Section 4.2 exhibit the consistency of the numerical solutions. The effect of the inhomogeneity and anisotropy of materials as well as the obtained steady-state solutions are as expected.

Funding

The APC was funded by Hasanuddin University, Makassar, Indonesia.

Data Availability Statement

The data that supports the findings of this study are available within the article.

Acknowledgments

This work is supported by Hasanuddin University and Ministry of Education, Culture, Research, and Technology of Indonesia.

Conflicts of Interest

The author declares that there is no conflicts of interest for publishing the manuscript.

Abbreviations

The following abbreviations are used in this manuscript:
FGM Functionally Graded Material
BEM Boundary Element Method
LT Laplace Transform
DCR Diffusion Convection Reaction

References

  1. Fendoglu, H.; Bozkaya, C.; Tezer-Sezgin, M. DBEM and DRBEM solutions to 2D transient convection-diffusion-reaction type equations. Eng. Anal. Boundary Elem. 2018, 93, 124–134. [Google Scholar] [CrossRef]
  2. Wang, X.; Ang W-T. A complex variable boundary element method for solving a steady-state advection–diffusion–reaction equation. Appl. Math. Comput. 2018, 321, 731–744. [Google Scholar] [CrossRef]
  3. Sheu, T.W.H.; Wang, S.K.; Lin, R.K. An Implicit Scheme for Solving the Convection–Diffusion–Reaction Equation in Two Dimensions. J. Comput. Phys. 2000, 164, 123–142. [Google Scholar] [CrossRef]
  4. Xu, M. A modified finite volume method for convection-diffusion-reaction problems. Int. J. Heat Mass Transfer 2018, 117, 658–668. [Google Scholar] [CrossRef]
  5. AL-Bayati, S.A.; Wrobel, L.C. Radial integration boundary element method for two-dimensional non-homogeneous convection–diffusion–reaction problems with variable source term. Eng. Anal. Boundary Elem. 2019, 101, 89–101. [Google Scholar] [CrossRef]
  6. Samec, N.; Škerget, L. Integral formulation of a diffusive–convective transport equation for reacting flows. Eng. Anal. Boundary Elem. 2004, 28, 1055–1060. [Google Scholar] [CrossRef]
  7. Rocca, A.L.; Rosales, A.H.; Power, H. Radial basis function Hermite collocation approach for the solution of time dependent convection–diffusion problems. Eng. Anal. Boundary Elem. 2005, 29, 359–370. [Google Scholar] [CrossRef]
  8. AL-Bayati, S.A.; Wrobel, L.C. The dual reciprocity boundary element formulation for convection-diffusion-reaction problems with variable velocity field using different radial basis functions. Int. J. Mech. Sci. 2018, 145, 367–377. [Google Scholar] [CrossRef]
  9. AL-Bayati, S.A.; Wrobel, L.C. A novel dual reciprocity boundary element formulation for two-dimensional transient convection–diffusion–reaction problems with variable velocity. Eng. Anal. Boundary Elem. 2018, 94, 60–68. [Google Scholar] [CrossRef]
  10. Hernandez-Martinez, E.; Puebla, H.; Valdes-Parada, F.; Alvarez-Ramirez, J. Nonstandard finite difference schemes based on Green’s function formulations for reaction–diffusion–convection systems. Chem. Eng. Sci. 2013, 94, 245–255. [Google Scholar] [CrossRef]
  11. Azis, M.I. Standard-BEM solutions to two types of anisotropic-diffusion convection reaction equations with variable coefficients. Eng. Anal. Boundary Elem. 2019, 105, 87–93. [Google Scholar] [CrossRef]
  12. Hassanzadeh, H.; Pooladi-Darvish, M. Comparison of different numerical Laplace inversion methods for engineering applications. Appl. Math. Comput. 2007, 189, 1966–1981. [Google Scholar] [CrossRef]
Figure 1. The boundary conditions for Problem 4.1
Figure 1. The boundary conditions for Problem 4.1
Preprints 70886 g001
Figure 2. Top: The aggregate relative error E A of the numerical solutions c with N = 6 , 8 , 10 , 12 for Case 1 (left) and zoom-in view for N = 8 , 10 (right). Bottom: The individual relative errors E I at t = π / 2 (left) and t = π (right) with N = 10 .
Figure 2. Top: The aggregate relative error E A of the numerical solutions c with N = 6 , 8 , 10 , 12 for Case 1 (left) and zoom-in view for N = 8 , 10 (right). Bottom: The individual relative errors E I at t = π / 2 (left) and t = π (right) with N = 10 .
Preprints 70886 g002
Figure 3. Top: The aggregate relative error E A of the numerical solutions c / x 1 with N = 6 , 8 , 10 , 12 for Case 1 (left) and zoom-in view for N = 6 , 8 , 10 (right). Bottom: The individual relative errors E I at t = π / 2 (left) and t = π (right) with N = 6 .
Figure 3. Top: The aggregate relative error E A of the numerical solutions c / x 1 with N = 6 , 8 , 10 , 12 for Case 1 (left) and zoom-in view for N = 6 , 8 , 10 (right). Bottom: The individual relative errors E I at t = π / 2 (left) and t = π (right) with N = 6 .
Preprints 70886 g003

Figure 4. Top: The aggregate relative error E A of the numerical solutions c / x 2 with N = 6 , 8 , 10 , 12 for Case 1 (left) and zoom-in view for N = 6 , 8 , 10 (right). Bottom: The individual relative errors E I at t = π / 2 (left) and t = π (right) with N = 10 .
Figure 4. Top: The aggregate relative error E A of the numerical solutions c / x 2 with N = 6 , 8 , 10 , 12 for Case 1 (left) and zoom-in view for N = 6 , 8 , 10 (right). Bottom: The individual relative errors E I at t = π / 2 (left) and t = π (right) with N = 10 .
Preprints 70886 g004

Figure 5. Top: The aggregate relative error E A of the numerical solutions c with N = 6 , 8 , 10 , 12 for Case 2 (left) and zoom-in view for N = 10 , 12 (right). Bottom: The individual relative errors E I at t = π / 2 (left) and t = π (right) with N = 10 .
Figure 5. Top: The aggregate relative error E A of the numerical solutions c with N = 6 , 8 , 10 , 12 for Case 2 (left) and zoom-in view for N = 10 , 12 (right). Bottom: The individual relative errors E I at t = π / 2 (left) and t = π (right) with N = 10 .
Preprints 70886 g005

Figure 6. Top: The aggregate relative error E A of the numerical solutions c / x 1 with N = 6 , 8 , 10 , 12 for Case 2 (left) and zoom-in view for N = 10 , 12 (right). Bottom: The individual relative errors E I at t = π / 2 (left) and t = π (right) with N = 10 .
Figure 6. Top: The aggregate relative error E A of the numerical solutions c / x 1 with N = 6 , 8 , 10 , 12 for Case 2 (left) and zoom-in view for N = 10 , 12 (right). Bottom: The individual relative errors E I at t = π / 2 (left) and t = π (right) with N = 10 .
Preprints 70886 g006

Figure 7. Top: The aggregate relative error E A of the numerical solutions c / x 2 with N = 6 , 8 , 10 , 12 for Case 2. Bottom: The individual relative errors E I at t = π / 2 (left) and t = π (right) with N = 6 .
Figure 7. Top: The aggregate relative error E A of the numerical solutions c / x 2 with N = 6 , 8 , 10 , 12 for Case 2. Bottom: The individual relative errors E I at t = π / 2 (left) and t = π (right) with N = 6 .
Preprints 70886 g007

Figure 8. Top: The aggregate relative error E A of the numerical solutions c with N = 6 , 8 , 10 , 12 for Case 3 (left) and zoom-in view for N = 10 , 12 (right). Bottom: The individual relative errors E I at t = π / 2 (left) and t = π (right) with N = 12 .
Figure 8. Top: The aggregate relative error E A of the numerical solutions c with N = 6 , 8 , 10 , 12 for Case 3 (left) and zoom-in view for N = 10 , 12 (right). Bottom: The individual relative errors E I at t = π / 2 (left) and t = π (right) with N = 12 .
Preprints 70886 g008

Figure 9. Top: The aggregate relative error E A of the numerical solutions c / x 1 with N = 6 , 8 , 10 , 12 for Case 3 (left) and zoom-in view for N = 10 , 12 (right). Bottom: The individual relative errors E I at t = π / 2 (left) and t = π (right) with N = 12 .
Figure 9. Top: The aggregate relative error E A of the numerical solutions c / x 1 with N = 6 , 8 , 10 , 12 for Case 3 (left) and zoom-in view for N = 10 , 12 (right). Bottom: The individual relative errors E I at t = π / 2 (left) and t = π (right) with N = 12 .
Preprints 70886 g009

Figure 10. Top: The aggregate relative error E A of the numerical solutions c / x 2 with N = 6 , 8 , 10 , 12 for Case 3 (left) and zoom-in view for N = 10 , 12 (right). Bottom: The individual relative errors E I at t = π / 2 (left) and t = π (right) with N = 12 .
Figure 10. Top: The aggregate relative error E A of the numerical solutions c / x 2 with N = 6 , 8 , 10 , 12 for Case 3 (left) and zoom-in view for N = 10 , 12 (right). Bottom: The individual relative errors E I at t = π / 2 (left) and t = π (right) with N = 12 .
Preprints 70886 g010

Figure 11. The boundary conditions for Problem 4.2.
Figure 11. The boundary conditions for Problem 4.2.
Preprints 70886 g011
Figure 12. Solutions c 0.5 , 0.3 , t and c 0.5 , 0.7 , t for all cases of Problem Section 4.2
Figure 12. Solutions c 0.5 , 0.3 , t and c 0.5 , 0.7 , t for all cases of Problem Section 4.2
Preprints 70886 g012aPreprints 70886 g012b
Figure 13. Solutions c 0.5 , 0.5 , t for Case 1 and 4 of Problem Section 4.2
Figure 13. Solutions c 0.5 , 0.5 , t for Case 1 and 4 of Problem Section 4.2
Preprints 70886 g013
Table 1. Values of V m of the Stehfest formula for N = 6 , 8 , 10 , 12
Table 1. Values of V m of the Stehfest formula for N = 6 , 8 , 10 , 12
V m N = 6 N = 8 N = 10 N = 12
V 1 1 1 / 3 1 / 12 1 / 60
V 2 49 145 / 3 385 / 12 961 / 60
V 3 366 906 1279 1247
V 4 858 16394 / 3 46871 / 3 82663 / 3
V 5 810 43130 / 3 505465 / 6 1579685 / 6
V 6 270 18730 236957.5 1324138.7
V 7 35840 / 3 1127735 / 3 58375583 / 15
V 8 8960 / 3 1020215 / 3 21159859 / 3
V 9 164062.5 8005336.5
V 10 32812.5 5552830.5
V 11 2155507.2
V 12 359251.2
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated