Preprint
Article

Analytical and Numerical Investigation of two-dimensional Heat Transfer with Periodic Boundary Conditions

Altmetrics

Downloads

97

Views

25

Comments

0

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

10 November 2023

Posted:

10 November 2023

You are already at the latest version

Alerts
Abstract
Two-dimensional heat diffusion problem with heat source which is quasilinear parabolic problem is examined analytically and numerically. Periodic boundary conditions are used as a boundary conditions. Since the problem is not linear, Picard’s successive approximation theorem is used. Under certain conditions of natural regularity and consistency imposed on the input data, establish the existence, uniqueness and constant dependence of the solution on the data using the generalized Fourier method. As a numerical solution, implicit finite difference scheme is used. The results ob-tained from analytical and the numerical solutions are so close to each other.
Keywords: 
Subject: Engineering  -   Mechanical Engineering

1. Introduction

The investigation of mathematical model for a great number of important applications such as chemical diffusions [1,2] application in the heat conduction problems, application of heat conduction problems [3,4,5], population dynamics [6], thermoelasticity [7], medical science, electrochemistry [8], engineering, and control theory require the analysis of the analyses of two-dimensional parabolic partial differential equations with nonlocal boundary conditions [9,10,11].
Heat diffusion equation is one of the two-dimensional quasilinear parabolic problem. Heat diffusion equation is used for determining temperature distribution in every region of the domain at desired time for conduction heat transfer problems. Heat transfer by conduction consists of a temperature differences in a solid or stationary fluid. Information of temperature distribution for a solid can be used to determine structural integrity through the determination of thermal stresses, expansions, deflections and cracks. Evaluation and propagation of granite thermal stress and thermal crack with time are shown in Figure 1 under different temperatures conditions [12]. The temperature distribution could also be used for optimizing the thickness of an insulating material or for determining the compatibility of the material and the special coatings or adhesives used [13].
Initial and boundary conditions are required for solving the two-dimensional quasilinear parabolic problem such as the heat diffusion equation. Some of boundary conditions for heat diffusion equation, Dirichlet (or first kind) determines the fixed temperature at surfaces. For example, the surface is very close to this when it comes into contact with a melting solid and/or boiling liquid. In both cases, heat transfer occurs at the surface, while the surfaces remain at the temperature of phase changing process [13]. Neumann (or the second type) corresponds to the presence of constant heat flux on the surface. This heat flux is related to the temperature at the surface according to Fourier’s law. It can be accomplished by bonding a thin electric heater to the surface. A special case of this situation corresponds to a perfectly insulated or adiabatic surface. Periodic boundary condition is combined with Dirichlet and Neumann boundary conditions [14], and it is set to isolate repeating temperature distribution in solution domain [15].
A special case of the nonlocal boundary [16] condition, the periodic boundary condition, is used in present investigation. Generally, periodic boundary condition is often used in numerical simulations and mathematical models. Also, the periodic boundary condition results from many important applications in heat transfer, fluid dynamics, life sciences and used lunar theory [14,17,18].
In order to solve two-dimensional quasilinear parabolic problem, there are several numerical methods, such as finite difference methods [19,20], finite element methods [21,22], finite volume methods [23], lattice Boltzmann Method [15,24,25]. Denghan solve the one-dimensional heat diffusion equation numerically with finite difference method [26]. In this study, three different schemes were used, namely, backward time centered space (BTCS) implicit scheme, the implicit Crandll’s method, forward time centered space (FTCS) explicit scheme and the Dufort-Frankel three-level techniques. Crank-Nicolson implicit scheme is used to solve numerically for one-dimensional heat diffusion equation with inverse coefficient by Baglan et al [27]. Kanca and Baglan solve the two-dimensional heat diffusion equation with periodic boundary conditions analytically and numerically [28]. Fourier method and implicit finite difference scheme are used for analytical and numerical solutions. Hamila et al. [29] searched the effect of change in thermal conductivity on several transient heat transfer problems with using lattice Boltzmann method (LBM). Benchmark problems containing conduction and/or radiation with constant thermal conductivity were calculated and simulated. Heat diffusion equation were solved numerically using exclusively the lattice Boltzmann method. Numerical results are so close with available result in the literature.
The finite difference method is one of the existing numerical methods which is considered in present investigation can be applied to solve partial differential equations. Finite difference method is based on the discretization of differential equations by finite difference equations. Finite difference approximations have algebraic forms and relate the value of dependent variable at one point in the solution region, to the values at some neighboring points. Using Taylor series is the most usual way to construct these approximations. The numerical method recommended here is the implicit finite difference method. This method supports second order accurate in the spatial grid sizes and first order time grid size. The explicit finite scheme has restriction of determining time step size due to stability requirements for the numerical solution of the two-dimensional heat diffusion equation.
In the present paper, we prove the existence, uniqueness, and continuous dependence of the solution on the data and obtain the numerical solution by Fourier Method and Picard’s successive approximation of the two-dimensional heat diffusion problem with periodic boundary conditions [14,27,28]. For numerical solution, implicit finite difference approximation is used. Numerical solution has a good agreement with analytical solution.

2. The Problem with Periodic Boundary Conditions

Let’s examine the problem given below;
Preprints 90194 i001
Preprints 90194 i002
Preprints 90194 i003
Preprints 90194 i004
Preprints 90194 i005
where Fourier Coefficients: Preprints 90194 i006 Then we obtain the solution:
Preprints 90194 i007

3. Analytical solution of the Problem

Let's assume the following rules for the functions used in the problem: Preprints 90194 i008 where Preprints 90194 i009
Definition 1. 
Preprints 90194 i010is called Banach (B) norm.
Theorem 1. 
Let the assumptions (A)-(B) be satisfied. Then the problem has a unique solution.
Proof. 
As iteration for the problem Preprints 90194 i011
Applying Cauchy, Bessel, Hölder inequalities, Lipschitz condition, we have
Preprints 90194 i012
From the theorem T ( 0 ) ( τ ) B .
Applying for the step N,
Preprints 90194 i013
We get T ( N ) ( τ ) B since T ( N + 1 ) ( τ ) B .
Let us show that, T ( N + 1 ) ( τ ) is converged for N . Preprints 90194 i014 Applying Cauchy, Bessel, Hölder inequalities, Lipschitz condition, we have
Preprints 90194 i015
where let Preprints 90194 i016.
Applying for the step N,
Preprints 90194 i017
We obtain T ( N + 1 ) ( τ ) T ( N ) ( τ ) .
Let us show that, T ( N + 1 ) ( τ ) T ( τ ) is converged for N .
Applying Cauchy, Bessel, Hölder inequalities, Lipschitz condition, we have
Preprints 90194 i018
From the Gronwall’s inequality,
Preprints 90194 i019
From Eq. 12, we get T ( N + 1 ) ( τ ) T ( τ ) .
To show the uniqueness, we get two solutions the problem { T ( τ ) , U ( τ ) } .
Applying Cauchy, Bessel, Hölder inequalities, Lipschitz condition, we have Preprints 90194 i020 From the Gronwall’s inequality,
Preprints 90194 i021
We get T ( τ ) = U ( τ ) .

4. Stability of Solution

Theorem 2. 
Let the assumptions (A)-(B) be satisfied. Then the problem is constantly dependent on the data.
Proof. 
Preprints 90194 i022 Applying Cauchy, Bessel, Hölder inequalities and Lipschitz condition, we have Preprints 90194 i023 From the Gronwall’s inequality,
Preprints 90194 i024
From Eq. 14, for θ θ ¯ . Then T T ¯ .

5. Numerical Method for Problem

In this section, we use implicit finite-difference approximation of the discretized problem (1)-(4);
Preprints 90194 i025
Preprints 90194 i026
Preprints 90194 i027
Preprints 90194 i028
where the computational domain [ 0 , π ] × [ 0 , π ] × [ 0 , τ ] is discretized as follows: Preprints 90194 i029
Where Δ α = π / n α , Δ β = π / n β and Δ τ = τ / n are the space α direction, space β direction and time steps, respectively. nx, ny and n is are three positive integers. Preprints 90194 i030
In order to define periodic boundary condition (Eq. 17 and 18) for implicit finite difference scheme, one-dimensional schematic figure with numerical meshes is represented in Figure 2. Point-1 and point-nα are boundary points of one-dimensional solution domain. When we determine spatial dicretization with finite difference scheme, the used finite difference formulations are defined as;
Preprints 90194 i031
where r α = Δ τ / Δ α 2 , and the coefficient matrix can be rewrittwen for one-dimensional solution domain.
Preprints 90194 i032
This one-dimensional formulation and coefficient matrix can be extended in two-dimensional solution domain. Two-dimensional schematic figure with numerical meshes (9 finite difference meshes) is represented in Figure 3. Of course, matrix is constructed for 9 finite difference meshes, therefore the matrix size is 9×9. Point-1, 2 and 3 are periodic with 7, 8 and 9, respectively, at β-direction. Point-1, 4 and 7 are periodic with 3, 6 and 9, respectively, at α direction. Coefficient matrix can be written as;
Preprints 90194 i033
where r β = Δ τ / Δ β 2 . If large number of meshes are used, the matrix can be extended in a same way. For example, if we have 20 meshes both α and β directions, our matrix size will be 400×400, total element of matrix is 160 000, and our matrix can be constructed in a same manner to define points of periodic boundary condition and inner (no periodic boundary condition points).
Example 1. 
Consider the heat diffision problem with heat source (1)-(4):Preprints 90194 i034
We determine the time step size is Δ τ = 0.005 . Discretization of spaces or mesh sizes of α and β directions are considered equal ( Δ α = Δ β ). In order to determine, mesh size of α and β directions, grid independency study is done. For grid-independence study, fifteen different grid resolutions are used.
In α and β directions, the same numbers of grid are used because of equal grid size ( Δ α = Δ β ) are considered. For every grid resolutions, we observed area-weighted averaged value of temperature for whole solution domain. Figure 4 shows the grid independent study of area-weighted averaged value of temperature. According to Figure 4, mesh number of 6400 is taken as grid independent mesh number. In this case, eighty grids are used for α and β directions, therefore our mesh size of α and β directions are get as 0.03926875.
Figure 5 shows (a) exact and (b) numerical solutions for temperature with time at 0.1s. The exact and numerical solutions are paraboloid. Maximum and minimum temperature values are observed at corners and at the middle of the solution domain, respectively. Because of periodic boundary conditions at side of solution domain, temperatures are same at side of the solution domain. Also, exact and numerical solutions are so quite similar.
Figure 6 shows (a) absolute error and (b) relative absolute error for temperature with time at 0.1s. Absolute error is the absolute value of the difference between exact solution value and numerical solution. According to Figure 6a, maximum absolute error is shown at center of the solution domain, also at corners of the solution domain, absolute error is higher. At center of solution domain, approximately absolute error of 0.1 is observed. Relative absolute error is the ratio of the absolute error to the exact solution. The magnitude of the absolute error in terms of the exact solution is determined with using relative absolute error. Like in absolute error, maximum relative error which is 0.03 is observed at the center of the solution domain. At corners of the solution domain, relative absolute error is order of 0.01. As the time increases, absolute error and relative absolute error remain in the same order, and these error are in reasonable.
Figure 7 presents the numerical predictions of temperature distributions at time (a) 0.1s, (b) 0.25s and (c) 0.5s. As one can see, temperature distributions at the specific times are paraboloid. The concavity of the parabolite increases with time. The minimum and the maximum values of temperature are appeared at center and corners of the solution domain, respectively, at all times.
In order to see all times in one figure, Figure 8 presents the distribution of temperature with time range from 0s to 0.5s. Paraboloid are not shown here. Maximum and minimum temperature values are observed at the corners and center of the solution domain, respectively, at all times. Also with increasing time, temperature values of corners increase and temperature values of centers decreases.

7. Conclusions

Analytical and numerical investigation of two-dimensional heat diffusion problem with heat source is done. This problem is quasi-linear parabolic problem, and we use initial condition and periodic boundary condition to determine temperature in the solution domain. Because of the problem is nonlinear, Picard’s successive approximation theorem is used. Under certain conditions of natural regularity and consistency imposed upon the input data, we constituie the existence, uniqueness and constant dependence of the solution on the data using the generalized Fourier method. Implicit finite difference scheme is used for the numerical solution. The number of numerical meshes where results do not change is determined according to grid independence study. In the light of analytical and numerical solution, distribution of temperature is paraboloid at the certain time. With increasing time, concavity of parabolic increases. At the time of 0.1s, maximum absolute error and relative error are occurred at the middle of the solution domain, they are, 0.1, 0.03, respectively. Therefore, we can say that analytical solution and numerical solutions are close each other. Also, with increasing time, degrees of absolute and relative absolute errors are same.

Nomenclature

Preprints 90194 i035

References

  1. Permikin, D.V.; Zverev, V.S. Mathematical model on surface reaction diffusion in the presence of chemical reaction. International Journal of Heat and Mass Transfer 2013, 57, 215-221. [CrossRef]
  2. Tredenick, E.C.; Farewell, T.W.; Forster, W.A. Modeling of diffusion of a hydrophilic ionic fertilizer in plant cuticles: surfactant and hygroscopic effects. Front Plant Sci. 2018, 9, 1888. [CrossRef]
  3. Canon, J.R. The solution of the heat equation subject to the specification of energy. Quarterly of Appl. Mathematics 1963, 21(2), 155-160. [CrossRef]
  4. Canon, J.R. Determination of an unknown heat source from overspecified boundary data. SIAM Journal on Numerical Analysis 1968, 5(2), 275-286. [CrossRef]
  5. Zhu, Q.; Tran, H.; Yang, H. Heat conduction: Mathematical modeling and experimental data. J Emerging Investigators 2021, 4, 1-4. [CrossRef]
  6. Venkatesha, P.; Eunice, G.B.; Akshaya, B.; Arya Kumari, S. Mathematical modeling of population growth. Int. J. Sci. Eng Management 2017, 2(11), 117-121.
  7. Avey, M.; Fantuzzi, N.; Sofiyev, A. Mathematical modeling and analytical solution of thermoelectric stability problem of functionally graded nanocomposite cylinders within different theories. Mathematics 2022, 10, 1081.
  8. Turnbull, A.; Ferriss, D.H. Mathematical modelling of the electrochemistry in corrosion fatigue cracks in structural steel cathodically protected in sea water. Corrosion Science 1986, 26(8), 601-628. [CrossRef]
  9. Dehghan, M. Implicit locally one-dimensional methods for two-dimensional diffusion with a non-local boundary condition. Math. Computers Simulation 1999, 49, 331-349. [CrossRef]
  10. Siddique, M. Solving two-dimensional diffusion equations with nonlocal boundary conditions by a special class of Padé approximants. Systemics, Cybernetics and Informatics 2010, 8(4), 23-29.
  11. Mohebbi, A. A numerical algorithm for determination of a control parameter in two-dimensional parabolic inverse problems. ACTA Mathematics Appicatae Sinica, English Series 2015, 31(1), 213-224. [CrossRef]
  12. Li, Q.; Yin, T.; Li, X.; Shu, R. Experimental and numerical investigation on thermal damage of granite subjected to heating and cooling. Mathematics 2021, 9, 3027. [CrossRef]
  13. Bergmann, T.L.; Lavine, A.S.; Incropera, F.P.; Dewitt, D.P. Fundamentals of Heat and Mass Transfer, 7th ed.; John Wiley&Sons: Danver, USA, 2011.
  14. Baglan, I.; Kanca, F. Two-dimensional inverse quasilinear parabolic problems with periodic boundary conditions. Applicable Analysis 2019, 98(8), 1549-1565. https://doi.org/10.1080/00036811.2018.1434149. [CrossRef]
  15. Mohamad, A.A. Lattice Boltzmann Method Fundamentals and Engineering Applications with Computer Codes, 2nd ed.; Springer-Verlag London Ltd.: London, England, 2019.
  16. Afshar, S.; Soltanalizadeh, B. Solution of the two-dimensional second-order diffusion equation with nonlocal boundary conditions. Int. J. Pure Appl. Math. 2014, 94(2), 119-131. https://doi.org/10.12732/ijpam.v94i2.1. [CrossRef]
  17. Rao, Y.; Zhang, Y.; Xu, Y.; Ke, H. Experimental study and numerical analysis of heat transfer enhancement and turbulent flow over shallowly dimples channel surfaces. Int. J. Heat Mass Transf. 2020, 160, 120195. [CrossRef]
  18. Hill, G.W. On the part of the motion of the lunar perigee which is a function of the mean motions of the sun and moon. Acta Math 1886, 8, 1-36. [CrossRef]
  19. Morton, K.W.; Mayers, D.F. Numerical Solution of Partial Differential Equation, 2nd ed.; Cambridge University Press: Cambridge, England, 2005.
  20. Dmitriev, V.G.; Danilin, A.N.; Popova, A.R.; Pshenichnova, N.V. Numerical analysis of deformation characteristics of elastic inhomogeneous rotational shells at arbitrary displacements and rotating angles. Computation 2022, 10, 184. [CrossRef]
  21. Benim, A.C.; Zinser, W. Investigation into finite element analysis of confined turbulent flows using a k-ε model of turbulence. Compt. Methods Appl. Mech. Engineering 1985, 51(1-3), 507-523. [CrossRef]
  22. Benim, A.C.; Finite element analysis of confined turbulent swirling flows. Int. J. Num. Meth. Fluids 1990, 11(6), 697-717.
  23. Versteeg, H.K.; Malalasekera, W. An Introduction to Computational Fluid Dynamics, 2nd ed.; Pearson, Prentice Hall: London, England, 2007.
  24. Chen, S.; Doolen, G.D. Lattice Boltzmann method for fluid flows. Annual Review Fluid Mech. 1998, 30, 329-364. [CrossRef]
  25. Succi, S. The Lattice Boltzmann Equation for fluid Dynamics and Beyond; Oxford University Press: Oxford, England, 200.
  26. Denghan, M. Efficient techniques for the second-order parabolic equation subject to nonlocal specifications . Appl. Num. Math. 2005, 12(1), 39-62.
  27. Baglan, I.; Kanca, F.; Mishra, V.N. Determination of an unknown heat source from integral overdetermination condition. Iran J. Sci. Technol. Trans. Sci. 2018, 42, 1373-1382. [CrossRef]
  28. Kanca, F.; Baglan, I. Solution of the boundary-value problem of heat conduction with parabolic conditions. Ukrainian Math. J. 2020, 72(2), 232-245.
  29. Hamila, R.; Chaabane, R.; Askri, F.; Jemni, A.; Nasrallah, S. B. Lattice Boltzmann method for heat transfer problems with variable thermal conductivity. Int. J. Heat Tech. 2017, 35(2), 313-324. [CrossRef]
Figure 1. Evaluation and propagation of granite thermal stress field and thermal cracks [12].
Figure 1. Evaluation and propagation of granite thermal stress field and thermal cracks [12].
Preprints 90194 g001
Figure 2. One-dimensional schematic figure with numerical meshes.
Figure 2. One-dimensional schematic figure with numerical meshes.
Preprints 90194 g002
Figure 3. Two-dimensional schematic figure with numerical meshes
Figure 3. Two-dimensional schematic figure with numerical meshes
Preprints 90194 g003
Figure 4. Grid independence study of area-weighted averaged value of temperature.
Figure 4. Grid independence study of area-weighted averaged value of temperature.
Preprints 90194 g004
Figure 5. The (a) exact and (b) numerical solutions for T(x, y, 0.1).
Figure 5. The (a) exact and (b) numerical solutions for T(x, y, 0.1).
Preprints 90194 g005
Figure 6. The (a) absolute error and (b) relative absolute error of T(α,β,0.1).
Figure 6. The (a) absolute error and (b) relative absolute error of T(α,β,0.1).
Preprints 90194 g006
Figure 7. Numerical predictions of T distributions with time (a) 0.1s, (b) 0.25s and (c) 0.5s.
Figure 7. Numerical predictions of T distributions with time (a) 0.1s, (b) 0.25s and (c) 0.5s.
Preprints 90194 g007
Figure 8. Distribution of temperature with time range from 0s to 0.5s.
Figure 8. Distribution of temperature with time range from 0s to 0.5s.
Preprints 90194 g008
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