scholarly journals Numerical Solution of Three-Dimensional Transient Heat Conduction Equation in Cylindrical Coordinates

2022 ◽  
Vol 2022 ◽  
pp. 1-8
Author(s):  
Endalew Getnet Tsega

Heat equation is a partial differential equation used to describe the temperature distribution in a heat-conducting body. The implementation of a numerical solution method for heat equation can vary with the geometry of the body. In this study, a three-dimensional transient heat conduction equation was solved by approximating second-order spatial derivatives by five-point central differences in cylindrical coordinates. The stability condition of the numerical method was discussed. A MATLAB code was developed to implement the numerical method. An example was provided in order to demonstrate the method. The numerical solution by the method was in a good agreement with the exact solution for the example considered. The accuracy of the five-point central difference method was compared with that of the three-point central difference method in solving the heat equation in cylindrical coordinates. The solutions obtained by the numerical method in cylindrical coordinates were displayed in the Cartesian coordinate system graphically. The method requires relatively very small time steps for a given mesh spacing to avoid computational instability. The result of this study can provide insights to use appropriate coordinates and more accurate computational methods in solving physical problems described by partial differential equations.

Filomat ◽  
2021 ◽  
Vol 35 (8) ◽  
pp. 2617-2628
Author(s):  
K.Y. Kung ◽  
Man-Feng Gong ◽  
H.M. Srivastava ◽  
Shy-Der Lin

The principles of superposition and separation of variables are used here in order to investigate the analytical solutions of a certain transient heat conduction equation. The structure of the transient temperature appropriations and the heat-transfer distributions are summed up for a straight mix of the results by means of the Fourier-Bessel arrangement of the exponential type for the investigated partial differential equation.


1987 ◽  
Vol 54 (4) ◽  
pp. 854-860 ◽  
Author(s):  
N. Noda ◽  
F. Ashida

The present paper deals with a transient thermoelastic problem for an axisymmetric transversely isotropic infinite solid with a penny-shaped crack. A finite difference formulation based on the time variable alone is proposed to solve a three-dimensional transient heat conduction equation in an orthotropic medium. Using this formulation, the heat conduction equation reduces to a differential equation with respect to the spatial variables. This formulation is applied to attack the transient thermoelastic problem for an axisymmetric transversely isotropic infinite solid containing a penny-shaped crack subjected to heat absorption and heat exchange through the crack surface. Thus, the thermal stress field is analyzed by means of the transversely isotropic potential function method.


2020 ◽  
Author(s):  
Yozo Mikata

Abstract Peridynamics for transient heat conduction problems in general anisotropic materials is developed. In order to develop a new peridynamic governing equation for heat conduction problems, the microconductivity (or microdiffusivity), which contains equivalent information as the constitutive equation for classical heat conduction, is determined by directly requiring the resulting peridynamic equation to converge to a classical heat conduction equation for anisotropic materials as the generalized material horizon approaches 0. Therefore, the convergence proof is built into the theory from the perspective of the governing equation. For the application of the newly obtained peridynamic governing equation, a time-dependent 3D peridynamic heat equation is analytically solved with two types of heat sources, and the results are discussed. These are believed to be the first exact analytical solutions for peridynamic heat conduction.


2020 ◽  
Vol 5 ◽  
pp. 36
Author(s):  
S. Abbas ◽  
A. A. Khan ◽  
B. Shakia

We study a comparison of serial and parallel solution of 2D-parabolic heat conduction equation using a Crank-Nicolson method with an Alternating Direction Implicit (ADI) scheme. The two-dimensional Heat equation is applied on a thin rectangular aluminum sheet. The forward difference formula is used for time and an averaged second order central difference formula for the derivatives in space to develop the Crank-Nicolson method. FORTRAN serial codes and parallel algorithms using OpenMP are used. Thomas tridigonal algorithm and parallel cyclic reduction methods are employed to solve the tridigonal matrix generated while solving heat equation. This paper emphasize on the run time of both algorithms and their difference. The results are compared and evaluated by creating GNU-plots (Command-line driven graphing utility).


Author(s):  
Almério José Venâncio Pains Soares Pamplona ◽  
Karoliny Freitas Silva ◽  
Cláudio Bucar Filho ◽  
Joel Vasco

Author(s):  
Seiji Nomura ◽  
Kosaku Kurata ◽  
Hiroshi Takamatsu

The irreversible electroporation (IRE) is a novel method to ablate abnormal cells by applying a high voltage between two electrodes that are stuck into abnormal tissues. One of the advantages of the IRE is that the extracellular matrix (ECM) may be kept intact, which is favorable for healing. For a successful IRE, it is therefore important to avoid thermal damage of ECM resulted from the Joule heating within the tissue. A three-dimensional (3-D) analysis was conducted in this study to predict temperature rise during the IRE. The equation of electric field and the heat conduction equation were solved numerically by a finite element method. It was clarified that the highest temperature rise occurred at the base of electrodes adjacent to the insulated surface. The result was significantly different from a two-dimensional (2-D) analysis due to end effects, suggesting that the 3-D analysis is required to determine the optimal condition.


Sign in / Sign up

Export Citation Format

Share Document