Search results

The purpose of this paper is to address unsteady heat conduction in two subsets of ordinary bodies. One subset consists of a large plane wall, a long cylinder and a sphere in one dimension. The other subset consists of a short cylinder and a large rectangular bar in two dimensions. The prevalent assumptions in the two subsets are: constant initial temperature, uniform surface heat flux and thermo-physical properties invariant with temperature. The engineering applications of the unsteady heat conduction deal with the determination of temperature–time histories in the two subsets using electric resistance heating, radiative heating and fire pool heating.

To this end, a novel numerical procedure named the enhanced method of discretization in time (EMDT) transforms the linear one-dimensional unsteady, heat conduction equations with non-homogeneous boundary conditions into equivalent nonlinear “quasi–steady” heat conduction equations having the time variable embedded as a time parameter. The equivalent nonlinear “quasi–steady” heat conduction equations are solved with a finite difference method.

Based on the numerical computations, it is demonstrated that the approximate temperature–time histories in the simple subset of ordinary bodies (large plane wall, long cylinder and sphere) exhibit a perfect matching over the entire time domain 0 < t < ∞ when compared against the rigorous exact temperature–time histories expressed by classical infinite series. Furthermore, using the method of superposition of solutions in the convoluted subset (short cylinder and large rectangular crossbar), the same level of agreement in the approximate temperature–time histories in the simple subset of ordinary bodies is evident.

The performance of the proposed EMDT coupled with a finite difference method is exhaustively assessed in the solution of the unsteady, one-dimensional heat conduction equations with prescribed surface heat flux for: a subset of one-dimensional bodies (plane wall, long cylinder and spheres) and a subset of two-dimensional bodies (short cylinder and large rectangular bar).

In this work, a modified thermoelastic model of heat conduction, including higher order of time derivative, is constructed by extending the Roychoudhuri model (TPL) (Choudhuri, 2007). In this new model, Fourier’s law of heat conduction is replaced by using Taylor series expansions, including three different phase lags for the heat flux, the thermal displacement and the temperature gradient. The generalized thermoelasticity models of Lord–Shulman (Lord and Shulman, 1967), Green and Naghdi (1991), dual-phase lag (Tzou, 1996) and three-phase lag (TPL) (Choudhuri, 2007) are obtained as special cases. The paper aims to discuss these issues.

The aim of this work is to establish a new generalized mathematical model of thermoelasticity that includes TPL in the vector of heat flux, and in the thermal displacement and temperature gradients extending TPL model (Li et al., 2019e). In this model, Fourier law of heat conduction is replaced by using Taylor series expansions to a modification of the Fourier law with introducing three different phase lags for the heat flux vector, the temperature gradient, and the thermal displacement gradient and keeping terms up with suitable higher orders.

The established high-order three-phase-lag heat conduction model reduces to the previous models of thermoelasticity as special cases.

In this paper, a TPL thermoelastic model is developed by extending the Roychoudhuri (Sherief and Raslan, 2017) model (TPL) considering the Taylor series approximation of the equation of heat conduction. This model is an alternative construction to the TPL model. The new model includes high order of TPL in the vector of heat flux, and in the thermal displacement and temperature gradients.

A finite element method for the analysis of combined radiative and conductive heat transport in a finite axisymmetric configuration is presented. The appropriate integro‐differential governing equations for a grey and non‐scattering medium with grey and diffuse walls are developed and solved for several model problems. We consider axisymmetric, cylindrical geometries with top and bottom boundaries of arbitrary convex shape. The method is accurate for media of any optical thickness and is capable of handling a wide array of axisymmetric geometries and boundary conditions. Several techniques are presented to reduce computational overhead, such as employing a Swartz‐Wendroff approximation and cut‐off criteria for evaluating radiation integrals. The method is successfully tested against several cases from the literature and is applied to some additional example problems to demonstrate its versatility. Solution of a free‐boundary, combined‐mode heat transfer problem representing the solidification of a semitransparent material, the Bridgman growth of an yttrium aluminium garnet (YAG) crystal, demonstrates the utility of this method for analysis of a complex materials processing system. The method is suitable for application to other research areas, such as the study of glass processing and the design of combustion furnace systems.

In this paper, a class of nonlinear heat‐conduction equations is derived. The properties of these heat‐conduction equations are analyzed. It is shown that the solutions of these equations contain singularity. That is, when T = T_m, discontinuity, i.e. blown‐up, occurs to the solutions. The occurrence of the blown‐ups is closely related the abnormal distribution of the initial ground temperatures, and so there might be some connections between blown‐ups and earthquakes.

Thermal conduction anisotropy, which is defined by the dependency of thermal conductivity on direction, is an important parameter in many engineering and research studies such as the design of nuclear waste depositional sites. In this context, the authors aim to investigate the effect of grain shape in thermal conduction anisotropy using pore scale modeling that utilizes real shapes of grains, pores and throats to characterize petrophysical properties of a porous medium.

The authors generalize the swelling circle approach to generate porous media composed of randomly arranged but regularly oriented elliptical grains at various grain ratios and porosities. Unlike previous studies that use fitting parameters to capture the effect of grain–grain thermal contact resistance, the authors apply roughness to grains’ surface. The authors utilize Lattice Boltzmann method to solve steady state heat conduction through medium.

Based on the results, when the temperature field is not parallel to either major or minor axes of grains, the overall heat flux vector makes a “deviation angle” with the temperature field. Deviation angle increases by augmenting the ratio of thermal conductivities of solid to fluid and the aspect ratios of grains. In addition, the authors show that porosity and surface roughness can considerably change the anisotropic properties of a porous medium whose grains are elliptical in shape.

The authors developed an algorithm for generation of non-circular-based porous medium with a novel approach to include grain surface roughness. In previous studies, the effect of grain contacts has been simulated using fitting parameters, whereas in this work, the authors impose the roughness based on the its fractal geometry.

In this study, the purpose was to introduce two‐dimensional hyperbolic heat conduction equations in order to simulate the fast precooling process of a cylindrically shaped food product with internal heat generation. A modified model for internal heat generation due to respiration in the food product was proposed to take the effect of relaxation time into account. The obtained governing equations were solved numerically using an efficient finite difference technique. The influence of Biot number and heat generation parameters on thermal characteristics was examined and discussed. The results based on hyperbolic model were compared with the classical parabolic heat diffusion model. The present numerical code was validated via comparison with analytical solution and a good agreement was found.

The obtained governing equations were solved numerically using an efficient finite difference technique.

The influence of Biot number and heat generation parameters on thermal characteristics was examined and discussed. The results based on hyperbolic model were compared with the classical parabolic heat diffusion model. The present numerical code was validated via comparison with analytical solution and a good agreement was found.

Two‐dimensional analysis of fast precooling of cylindrical food product based on hyperbolic heat conduction model has not been investigated yet.

Gas insulated systems, such as gas insulated lines (GIL), use insulating gas, mostly sulfur hexalfluoride (SF₆), to enable a higher dielectric strength compared to e.g. air. However, under high voltage direct current conditions, charge accumulation and electric field stress may occur, which may lead to partial discharge or system failure. Therefore, numerical simulations are used to design the system and determine the electric field and charge distribution. Although the gas conduction shows a more complex current–voltage characteristic compared to solid insulation, the electric conductivity of the SF₆ gas is set as constant in most works. The purpose of this study is to investigate different approaches to address the conduction in the gas properly for numerical simulations.

In this work, two approaches are investigated to address the conduction in the insulating gas and are compared to each other. One method is an ion-drift-diffusion model, where the conduction in the gas is described by the ion motion in the SF₆ gas. However, this method is computationally expensive. Alternatively, a less complex approach is an electro-thermal model with the application of an electric conductivity model for the SF₆ gas. Measurements show that the electric conductivity in the SF₆ gas has a nonlinear dependency on temperature, electric field and gas pressure. From these measurements, an electric conductivity model was developed. Both methods are compared by simulation results, where different parameters and conditions are considered, to investigate the potential of the electric conductivity model as a computationally less expensive alternative.

The simulation results of both simulation approaches show similar results, proving the electric conductivity for the SF₆ gas as a valid alternative. Using the electro-thermal model approach with the application of the electric conductivity model enables a solution time up to six times faster compared to the ion-drift-diffusion model. The application of the model allows to examine the influence of different parameters such as temperature and gas pressure on the electric field distribution in the GIL, whereas the ion-drift-diffusion model enables to investigate the distribution of homo- and heteropolar charges in the insulation gas.

This work presents numerical simulation models for high voltage direct current GIL, where the conduction in the SF₆ gas is described more precisely compared to a definition of a constant electric conductivity value for the insulation gas. The electric conductivity model for the SF₆ gas allows for consideration of the current–voltage characteristics of the gas, is computationally less expensive compared to an ion-drift diffusion model and needs considerably less solution time.

This investigation deals with the equilibration of heat conduction simulation in a very thin film using molecular dynamics. Two parameters, the positional order parameter and the kinetic H‐function, are employed simultaneously to monitor the evolution to the equilibrium. With the different boundary conditions, material parameters, and molecular lattice configurations, the results of the simulation show that the combination of the two parameters can give a consistent prediction to the approach of the equilibrium. At the equilibrium state, the process of heat conduction in a thin film is studied to understand the macroscopic behaviour from the standpoint of molecular dynamic motions. The method used can be applied to solve other microscopic flow and heat transfer problems using molecular dynamic simulation.

A numerical formulation for solving homogeneous anisotropic heat conduction problems based on the use of an isotropic fundamental solution is presented in detail. The analysis is carried out assuming a generic position of the coordinate axes, which may not coincide with the principal directions of orthotropy of the material. The two primary integral equations of the method are derived from the governing differential equation of the problem. Then, the numerical procedure is developed by rewriting the internal degrees of freedom that arise from the domain discretization in terms of the boundary nodes and solving the resulting system of linear equations for the boundary unknowns only. Special attention is given to the differentiation of singular integrals which yields additional terms as well as to the evaluation of the resulting Cauchy principal value integral. The main feature of the proposed formulation is its generality, which makes possible its direct extension to solve the problem of three‐dimensional heat conduction in anisotropic media and, foremost, to three‐dimensional orthotropic and anisotropic elasticity or elastoplasticity.

The present paper describes the applicability of hybrid transfinite element modelling/analysis formulations for non‐linear heat conduction problems involving phase change. The methodology is based on application of transform approaches and classical Galerkin schemes with finite element formulations to maintain the modelling versatility and numerical features for computational analysis. In addition, in conjunction with the above, the effects due to latent heat are modelled using enthalpy formulations to enable a physically realistic approximation to be effectively dealt computationally for materials exhibiting phase change within a narrow band of temperatures. Pertinent details of the approach and computational scheme adapted are described in technical detail. Numerical test cases of comparative nature are presented to demonstrate the applicability of the proposed formulations for numerical modelling/analysis of non‐linear heat conduction problems involving phase change.