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 study aims to highlight the behaviour of one-dimensional and two-dimensional fin models under the natural room conditions, considering the different values of dimensionless Biot number (Bi). The effect of convection and radiation on the heat transfer process has also been demonstrated using the meshless local Petrov–Galerkin (MLPG) approach.

It is true that MLPG method is time-consuming and expensive in terms of man-hours, as it is in the developing stage, but with the advent of computationally fast new-generation computers, there is a big possibility of the development of MLPG software, which will not only reduce the computational time and cost but also enhance the accuracy and precision in the results. Bi values of 0.01 and 0.10 have been taken for the experimental investigation of one-dimensional and two-dimensional rectangular fin models. The numerical simulation results obtained by the analytical method, benchmark numerical method and the MLPG method for both the models have been compared with that of the experimental investigation results for validation and found to be in good agreement. Performance of the fin has also been demonstrated.

The experimental and numerical investigations have been conducted for one-dimensional and two-dimensional linear and nonlinear fin models of rectangular shape. MLPG is used as a potential numerical method. Effect of radiation is also, implemented successfully. Results are found to be in good agreement with analytical solution, when one-dimensional steady problem is solved; however, two-dimensional results obtained by the MLPG method are compared with that of the finite element method and found that the proposed method is as accurate as the established method. It is also found that for higher Bi, the one-dimensional model is not appropriate, as it does not demonstrate the appreciated error; hence, a two-dimensional model is required to predict the performance of a fin. Radiative fin illustrates more heat transfer than the pure convective fin. The performance parameters show that as the Bi increases, the performance of fin decreases because of high thermal resistance.

Though, best of the efforts have been put to showcase the behaviour of one-dimensional and two-dimensional fins under nonlinear conditions, at different Bi values, yet lot more is to be demonstrated. Nonlinearity, in the present paper, is exhibited by using the thermal and material properties as the function of temperature, but can be further demonstrated with their dependency on the area. Additionally, this paper can be made more elaborative by extending the research for transient problems, with different fin profiles. Natural convection model is adopted in the present study but it can also be studied by using forced convection model.

Fins are the most commonly used medium to enhance heat transfer from a hot primary surface. Heat transfer in its natural condition is nonlinear and hence been demonstrated. The outcome is practically viable, as it is applicable at large to the broad areas like automobile, aerospace and electronic and electrical devices.

As per the literature survey, lot of work has been done on fins using different numerical methods; but to the best of authors’ knowledge, this study is first in the area of nonlinear heat transfer of fins using dimensionless Bi by the truly meshfree MLPG method.

The purpose of this paper is to address a novel method for solving parabolic partial differential equations (PDEs) in general, wherein the heat conduction equation constitutes an important particular case. The new method, appropriately named the Improved Transversal Method of Lines (ITMOL), is inspired in the Transversal Method of Lines (TMOL), with strong insight from the method of separation of variables.

The essence of ITMOL revolves around an exponential variation of the dependent variable in the parabolic PDE for the evaluation of the time derivative. As will be demonstrated later, this key step is responsible for improving the accuracy of ITMOL over its predecessor TMOL. Throughout the paper, the theoretical properties of ITMOL, such as consistency, stability, convergence and accuracy are analyzed in depth. In addition, ITMOL has proven to be unconditionally stable in the Fourier sense.

In a case study, the 1-D heat conduction equation for a large plate with symmetric Dirichlet boundary conditions is transformed into a nonlinear ordinary differential equation by means of ITMOL. The numerical solution of the resulting differential equation is straightforward and brings forth a nearly zero truncation error over the entire time domain, which is practically nonexistent.

Accurate levels of the analytical/numerical solution of the 1-D heat conduction equation by ITMOL are easily established in the entire time domain.

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.

The purpose of this paper is to develop a hybrid‐Trefftz (HT) finite element model (FEM) for simulating heat conduction in nonlinear functionally graded materials (FGMs) which can effectively handle continuously varying properties within an element.

In the proposed model, a T‐complete set of homogeneous solutions is first derived and used to represent the intra‐element temperature fields. As a result, the graded properties of the FGMs are naturally reflected by using the newly developed Trefftz functions (T‐complete functions in some literature) to model the intra‐element fields. The derivation of the Trefftz functions is carried out by means of the well‐known Kirchhoff transformation in conjunction with various variable transformations.

The study shows that, in contrast to the conventional FEM, the HT‐FEM is an accurate numerical scheme for FGMs in terms of the number of unknowns and is insensitive to mesh distortion. The method also performs very well in terms of numerical accuracy and can converge to the analytical solution when the number of elements is increased.

The value of this paper is twofold: a T‐complete set of homogeneous solutions for nonlinear FMGs has been derived and used to represent the intra‐element temperature; and the corresponding variational functional and the associated algorithm has been constructed.

A rigorous Finite Element (FE) formulation based on an enthalpy technique is developed for solving coupled nonlinear heat conduction/mass diffusion problems with phase change. The FE formulation consists of a fully coupled heat conduction and solute diffusion formulation, with solid‐liquid phase change, where the effects of pressure and convection are neglected. A full enthalpy method is employed eliminating singularities which result from abrupt changes in heat capacity at the phase interfaces. The FE formulation is based on the fixed grid technique where the elements are two dimensional, four noded quadrilaterals with the primary variables being enthalpy and average solute concentration. Temperature and solid mass fraction are calculated on a local level at each integration point of an element. A fully consistent Newton‐Raphson method is used to solve the global coupled equations and an Euler backward difference scheme is used for the temporal discretization. The solution of the enthalpy‐temperature relationship is carried out at the integration points using a Newton‐Raphson method. A secant method employing the regula falsi technique takes into account sudden jumps or sharp changes in the enthalpy‐temperature behaviour which occur at the phase zone interfaces. The Euler backward difference integration rule is used to calculate the solid mass fraction and its derivatives. A practical example is analysed and results are presented.

The purpose of this study is to compare interpolation algorithms and deep neural networks for inverse transfer problems with linear and nonlinear behaviour.

A series of runs were conducted for a canonical test problem. These were used as databases or “learning sets” for both interpolation algorithms and deep neural networks. A second set of runs was conducted to test the prediction accuracy of both approaches.

The results indicate that interpolation algorithms outperform deep neural networks in accuracy for linear heat conduction, while the reverse is true for nonlinear heat conduction problems. For heat convection problems, both methods offer similar levels of accuracy.

This is the first time such a comparison has been made.

The discretisation scheme presented for an unstructured triangular grid, is able to reflect the physical properties of stationary, nonlinear heat conduction even in the numerical case and leads to a stabilisation of iterative solution techniques for nonlinearities.

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.

A new way of solving the steady‐state coupled radiative‐conductive problem in semi‐transparent media is proposed. An angular discretization technique is applied in order to express the radiative transfer equation (RTE) in an inhomogeneous system of linear differential equations associated with Dirichlet boundary conditions. The system is solved by a direct method, after diagonalizing the characteristic matrix of the medium. The RTE is coupled with the nonlinear heat conduction equation. A simulation of a real semi‐transparent medium composed of silica fibers is illustrated. Comparison with results of other methods validates the new model. Moreover, the general scheme is easy to code and fast. The algorithm proved to be robust and stable.