Search results

This paper aims to develop an efficient numerical method for nonlinear transient heat conduction problems with local radiation boundary conditions and nonlinear heat sources.

Based on the physical characteristic of the transient heat conduction and the distribution characteristic of the Green’s function, a quasi-superposition principle is presented for the transient heat conduction problems with local nonlinearities. Then, an efficient method is developed, which indicates that the solution of the original nonlinear problem can be derived by solving some nonlinear problems with small structures and a linear problem with the original structure. These problems are independent of each other and can be solved simultaneously by the parallel computing technique.

Within a small time step, the nonlinear thermal loads can only induce significant temperature responses of the regions near the positions of the nonlinear thermal loads, whereas the temperature responses of the remaining regions are very close to zero. According to the above physical characteristic, the original nonlinear problem can be transformed into some nonlinear problems with small structures and a linear problem with the original structure.

An efficient and accurate numerical method is presented for transient heat conduction problems with local nonlinearities, and some numerical examples demonstrate the high efficiency and accuracy of the proposed method.

The purpose of this paper is to present a new approach of meshless local B-spline based finite difference (FD) method for solving two dimensional transient heat conduction problems.

In the present method, any governing equations are discretized by B-spline approximation which is implemented in the spirit of FD technique using a local B-spline collocation scheme. The key aspect of the method is that any derivative is stated as neighbouring nodal values based on B-spline interpolants. The set of neighbouring nodes are allowed to be randomly distributed thus enhanced flexibility in the numerical simulation can be obtained. The method requires no mesh connectivity at all for either field variable approximation or integration. Time integration is performed by using the Crank-Nicolson implicit time stepping technique.

Several heat conduction problems in complex domains which represent for extended surfaces in industrial applications are examined to demonstrate the effectiveness of the present approach. Comparison of the obtained results with solutions from other numerical method available in literature is given. Excellent agreement with reference numerical method has been found.

The method is presented for 2D problems. Nevertheless, it would be also applicable for 3D problems.

A transient two dimensional heat conduction in complex domains which represent for extended surfaces in industrial applications is presented.

The presented new meshless local method is simple and accurate, while it is also suitable for analysis in domains of arbitrary geometries.

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.

Topology optimization is a method used for developing optimized geometric designs by distributing material pixels in a given design space that maximizes a chosen quantity of interest (QoI) subject to constraints. The purpose of this study is to develop a problem-agnostic automatic differentiation (AD) framework to compute sensitivities of the QoI required for density distribution-based topology optimization in an unstructured co-located cell-centered finite volume framework. Using this AD framework, the authors develop and demonstrate the topology optimization procedure for multi-dimensional steady-state heat conduction problems.

Topology optimization is performed using the well-established solid isotropic material with penalization approach. The method of moving asymptotes, a gradient-based optimization algorithm, is used to perform the optimization. The sensitivities of the QoI with respect to design variables, required for optimization algorithm, are computed using a discrete adjoint method with a novel AD library named residual automatic partial differentiator (Rapid).

Topologies that maximize or minimize relevant quantities of interest in heat conduction applications are presented. The efficacy of the technique is demonstrated using a variety of realistic heat transfer applications in both two and three dimensions, in conjugate heat transfer problems with finite conductivity ratios and in non-rectangular/non-cuboidal domains.

In contrast to most published work which has either used finite element methods or Cartesian finite volume methods for transport applications, the topology optimization procedure is developed in a general unstructured finite volume framework. This permits topology optimization for flow and heat transfer applications in complex design domains such as those encountered in industry. In addition, the Rapid library is designed to provide a problem-agnostic pathway to automatically compute all required derivatives to machine accuracy. This obviates the necessity to write new code for finding sensitivities when new physics are added or new cost functions are considered and permits general-purpose implementations of topology optimization for complex industrial applications.

This paper aims to review the existing bond-based peridynamic (PD) and state-based PD heat conduction models, and further propose a refined bond-based PD thermal conduction model by using the PD differential operator.

The general refined bond-based PD is established by replacing the local spatial derivatives in the classical heat conduction equations with their corresponding nonlocal integral forms obtained by the PD differential operator. This modeling approach is representative of the state-based PD models, whereas the resulting governing equations appear as the bond-based PD models.

The refined model can be reduced to the existing bond-based PD heat conduction models by specifying particular influence functions. Also, the refined model does not require any calibration procedure unlike the bond-based PD. A systematic explicit dynamic solver is introduced to validate 1 D, 2 D and 3 D heat conduction in domains with and without a crack subjected to a combination of Dirichlet, Neumann and convection boundary conditions. All of the PD predictions are in excellent agreement with the classical solutions and demonstrate the nonlocal feature and advantage of PD in dealing with heat conduction in discontinuous domains.

The existing PD heat conduction models are reviewed. A refined bond-based PD thermal conduction model by using PD differential operator is proposed and 3 D thermal conduction in intact or cracked structures is simulated.

The purpose of this paper is to frame a dual-phase-lag model using the fractional theory of thermoelasticity with relaxation time. The generalized Fourier law of heat conduction based upon Tzou model that includes temperature gradient, the thermal displacement and two different translations of heat flux vector and temperature gradient has been used to formulate the heat conduction model. The microstructural interactions and corresponding thermal changes have been studied due to the involvement of relaxation time and delay time translations. This results in achieving the finite speed of thermal wave. Classical coupled and generalized thermoelasticity theories are recovered by considering the various special cases for different order of fractional derivatives and two different translations under consideration.

The work presented in this manuscript proposes a dual-phase-lag mathematical model of a thick circular plate in a finite cylindrical domain subjected to axis-symmetric heat flux. The model has been designed in the context of fractional thermoelasticity by considering two successive terms in Taylor’s series expansion of fractional Fourier law of heat conduction in the two different translations of heat flux vector and temperature gradient. The analytical results have been obtained in Laplace transform domain by transforming the original problem into eigenvalue problem using Hankel and Laplace transforms. The numerical inversions of Laplace transforms have been achieved using the Gaver−Stehfast algorithm, and convergence criterion has been discussed. For illustrative purpose, the dual-phase-lag model proposed in this manuscript has been applied to a periodically varying heat source. The numerical results have been depicted graphically and compared with classical, fractional and generalized thermoelasticity for various fractional orders under consideration.

The microstructural interactions and corresponding thermal changes have been studied due to the involvement of relaxation time and delay time translations. This results in achieving the finite speed of thermal wave. Classical coupled and generalized thermoelasticity theories are recovered by considering the various special cases for different order of fractional derivatives and two different translations under consideration. This model has been applied to study the thermal effects in a thick circular plate subjected to a periodically varying heat source.

A dual-phase-lag model can effectively be incorporated to study the transient heat conduction problems for an exponentially decaying pulse boundary heat flux and/or for a short-pulse boundary heat flux in long solid tubes and cylinders. This model is also applicable to study the various effects of the thermal lag ratio and the shift time. These dual-phase-lag models are also practically applicable in the problems of modeling of nanoscale heat transport problems of semiconductor devices and accordingly semiconductors can be classified as per their ability of heat conduction.

To the authors’ knowledge, no one has discussed fractional thermoelastic dual-phase-lag problem associated with relaxation time in a finite cylindrical domain for a thick circular plate subjected to an axis-symmetric heat source. This is the latest and novel contribution to the field of thermal mechanics.

We report on the progress in the development and application of a coupled boundary element/finite volume method temperature‐forward/flux‐back algorithm developed to solve conjugate heat transfer arising in 3D film‐cooled turbine blades. We adopt a loosely coupled strategy where each set of field equations is solved to provide boundary conditions for the other. Iteration is carried out until interfacial continuity of temperature and heat flux is enforced. The NASA‐Glenn explicit finite volume Navier‐Stokes code Glenn‐HT is coupled to a 3D BEM steady‐state heat conduction solver. Results from a CHT simulation of a 3D film‐cooled blade section are compared with those obtained from the standard two temperature model, revealing that a significant difference in the level and distribution of metal temperatures is found between the two. Finally, current developments of an iterative strategy accommodating large numbers of unknowns by a domain decomposition approach is presented. An iterative scheme is developed along with a physically‐based initial guess and a coarse grid solution to provide a good starting point for the iteration. Results from a 3D simulation show the process that converges efficiently and offers substantial computational and storage savings.

This paper aims to propose an efficient and accurate method to analyse the transient heat conduction in a periodic structure with moving heat sources.

The moving heat source is modelled as a localised Gaussian distribution in space. Based on the spatial distribution, the physical feature of transient heat conduction and the periodic property of structure, a special feature of temperature responses caused by the moving heat source is illustrated. Then, combined with the superposition principle of linear system, within a small time-step, computation of results corresponding to the whole structure excited by the Gaussian heat source is transformed into that of some small-scale structures. Lastly, the precise integration method (PIM) is used to solve the temperature responses of each small-scale structure efficiently and accurately.

Within a reasonable time-step, the heat source applied on a unit cell can only cause the temperature responses of a limited number of adjacent unit cells. According to the above feature and the periodic property of a structure, the contributions caused by the moving heat source for the most of time-steps are repeatable, and the temperature responses of the entire periodic structure can be obtained by some small-scale structures.

A novel numerical method is proposed for analysing moving heat source problems, and the numerical examples demonstrate that the proposed method is much more efficient than the traditional methods, even for larger-scale problems and multiple moving heat source problems.

An eigenvalue method is presented for solving the transient heat conduction problem with time‐dependent or time‐independent boundary conditions. The spatial domain is divided into finite elements and at each finite element node, a closed‐form expression for the temperature as a function of time can be obtained. Three test problems which have exact solutions were solved in order to examine the merits of the eigenvalue method. It was found that this method yields accurate results even with a coarse mesh. It provides exact solution in the time domain and therefore has none of the time‐step restrictions of the conventional numerical techniques. The temperature field at any given time can be obtained directly from the initial condition and no time‐marching is necessary. For problems where the steady‐state solution is known, only a few dominant eigenvalues and their corresponding eigenvectors need to be computed. These features lead to great savings in computation time, especially for problems with long time duration. Furthermore, the availability of the closed form expressions for the temperature field makes the present method very attractive for coupled problems such as solid—fluid and thermal—structure interactions.

The purpose of this paper is to derive the Trefftz functions (T-functions) for the Pennes’ equation and for the single-phase-lag (SPL) model (hyperbolic equation) with perfusion and then comparing field of temperature in a flat slab made of skin in the case when perfusion is taken into account, with the situation when a Fourier model is considered. When considering the process of heat conduction in the skin, one needs to take into account the average values of its thermal properties. When in biological bodies relaxation time is of the order of 20 s, the thermal wave propagation appears. The initial-boundary problems for Pennes’ model and SPL with perfusion model are considered to investigate the effect of the finite velocity of heat in the skin, perfusion and thickness of the slab on the rate of the thermal wave attenuation. As a reference model, the solution of the classic Fourier heat transfer equation for the considered problems is calculated. A heat flux has direction perpendicular to the surface of skin, considered as a flat slab. Therefore, the equations depend only on time and one spatial variable.

First of all the T-functions for the Pennes’ equation and for the SPL model with perfusion are derived. Then, an approximate solutions of the problems are expressed in the form of a linear combination of the T-functions. The T-functions satisfy the equation modeling the problem under consideration. Therefore, approximating a solution of a problem with a linear combination of n T-functions one obtains a function that satisfies the equation. The unknown coefficients of the linear combination are obtained as a result of minimization of the functional that describes an inaccuracy of satisfying the initial and boundary conditions in a mean-square sense.

The sets of T-functions for the Pennes’ equation and for the SPL model with perfusion are derived. An infinite set of these functions is a complete set of functions and stands for a base functions layout for the space of solutions for the equation used to generate them. Then, an approximate solutions of the initial-boundary problem have been found and compared to find out the effect of finite velocity of heat in the skin, perfusion and thickness of the slab on the rate of the thermal wave attenuation.

The methods used in the literature to find an approximate solution of any bioheat transfer problems are more complicated than the one used in the presented paper. However, it should be pointed out that there is some limitation concerning the T-function method, namely, the greater number of T-function is used, the greater condition number becomes. This limitation usually can be overcome using symbolic calculations or conducting calculations with a large number of significant digits.

The T-functions for the Pennes’ equation and for the SPL equation with perfusion have been reported in this paper for the first time. In the literature, the T-functions are known for other linear partial differential equations (e.g. harmonic functions for Laplace equation), but for the first time they have been derived for the two aforementioned equations. The results are discussed with respect to practical applications.