Local RBF-DQ method for two-dimensional transient heat conduction problems

The meshless local radial basis function-based differential quadrature (RBF-DQ) method is applied on two-dimensional heat conduction for different irregular geometries. This method is the combination of differential quadrature approximation of derivatives and function approximation of radial basis function. Four different geometries with regular and irregular boundaries are considered, and numerical results are compared with those gained by finite element (FE) solution achieved by COMSOL commercial code. Outcomes prove that current technique is in very good agreement with FEM and this fact that RBF-DQ method is an accurate and flexible method in solution of heat conduction problems.     

Heat conduction analysis for irregular functionally graded material geometries using the meshless weighted least-square method with temperature-dependent material properties

Abstract

This article presents the heat conduction analysis for irregular functionally graded material (FGM) with temperature-dependent material properties. For irregular FGM geometries, the meshless weighted least-square (MWLS) method is easy to model, implement, and interpolate those irregularly distributed field variables. To solve the heat conduction problem coupled with temperature-dependent FGM, the Laplace’s equilibrium equation and boundary condition become nonlinear. Thus, the Kirchhoff transformation is employed to convert the nonlinear problem to linear solution. MWLS method as a pure meshless analysis is then used to solve the linear equation of FGM geometries. Next, the temperature field is obtained by the inverse Kirchhoff transformation. Finally, the accuracy and effectiveness of the method were demonstrated by several numerical cases.     

Analysis of Conduction and Radiation Heat Transfer in a Differentially Heated 2‐D Square Enclosure

Radiative heat transfer with and without conduction in a differentially heated 2‐D square enclosure is analyzed. The enclosure with diffuse gray boundaries contains radiating and/or conducting gray homogeneous medium. Radiatively, the medium is absorbing, emitting and scattering. On the south boundary, four types of discrete heated regions, viz., the full boundary, the left one‐third, left two third and middle one third, are considered. In the absence of conduction, distributions of heat flux along the south boundary are studied for the effect of extinction coefficient. In the presence of conduction, distributions of radiation, conduction and total heat fluxes along the south boundary are analyzed for the effects of extinction coefficient, scattering albedo, conduction–radiation parameter, and south boundary emissivity. Effects of these parameters on centerline temperature distribution are also studied. To assess the performance of three commonly used radiative transfer methods, in all cases, the radiative transfer equation is solved using the discrete ordinate method (DOM), the conventional discrete ordinate method (CDOM) and the finite volume method (FVM). In the combined mode problem, with volumetric radiative information known from one of the three methods, viz., DOM, CDOM, and FVM, the energy equation is solved using the finite difference method (FDM). In all cases, the results from FDM‐DOM, FDM‐CDOM, and FDM‐FVM are in good agreement. Computationally, all three sets of methods are equally efficient.     

On preconditioned BiCGSTAB solver for MLPG method applied to heat conduction in complex geometry

Meshless local Petrov–Galerkin (MLPG) method is a promising meshfree method for continuum problems in complex domains, especially for large deformation, moving boundary and phase change problems. For large-scale problems, iterative methods for solving the discretized equations are more suitable than direct methods. Krylov subspace solvers of conjugate gradient type are the most preferred iterative solvers. The convergence rate of these methods depends on preconditioner used. Recently, proposed schedule relaxation Jacobi (SRJ) method can be used as a stand-alone solver and as a preconditioner. In the present work, the SRJ method is tested as a stand-alone solver and as a preconditioner for BiCGSTAB solver in the MLPG method, and its performance has been compared with successive overrelaxation (k) preconditioner. Two-dimensional linear steady-state heat conduction in complex shape geometry has been used as the model test problem.     

Estimation of the time-dependent heat flux using the temperature distribution at a point by conjugate gradient method

In this paper, the conjugate gradient method coupled with adjoint problem is used in order to solve the inverse heat conduction problem and estimation of the time-dependent heat flux using the temperature distribution at a point. Also, the effects of noisy data and position of measured temperature on final solution are studied. The numerical solution of the governing equations is obtained by employing a finite-difference technique. For solving this problem the general coordinate method is used. We solve the inverse heat conduction problem of estimating the transient heat flux, applied on part of the boundary of an irregular region. The irregular region in the physical domain (r,z) is transformed into a rectangle in the computational domain (ξ,η). The present formulation is general and can be applied to the solution of boundary inverse heat conduction problems over any region that can be mapped into a rectangle. The obtained results for few selected examples show the good accuracy of the presented method. Also the solutions have good stability even if the input data includes noise and that the results are nearly independent of sensor position.     

Computational fluid dynamics (CFD) modeling of heat transfer in a polymeric membrane using finite volume method

The efficiency, robustness and reliability of recent numerical methods for finding solutions to flow problems have given rise to the implementation of computational fluid dynamics (CFD) as a broadly used analysis method for engineering problems like membrane separation system. The CFD modeling in this study observes steady and unsteady (transient) heat flux and temperature profiles in a polymeric (cellulose acetate) membrane. This study is novel due to the implementation of user defined scalar (UDS) diffusion equation by using user-defined functions (UDFs) infinite volume method (FVM). Some details of the FVM used by the solver are carefully discussed when implementing terms in the governing equation and boundary conditions (BC). The contours of temperature due to high-temperature gradient are reported for steady and unsteady problems.     

Matrix free meshless method for transient heat conduction problems

In the paper, the element free Galerkin method (EFGM) is applied to calculate two-dimensional unsteady state heat conduction problems. As is well known, most of the meshless methods have higher computational cost than that of finite element method (FEM). In order to overcome this shortcoming especially for transient heat conduction problems, mass lumping procedure is adopted in EFGM, which can decrease the computational cost evidently. Moreover, this technique which can simplify the solution procedure makes the essential boundary conditions enforced directly. The results obtained by EFGM combining mass lumping technique are compared with those obtained by finite element method as well as analytical solutions, which shows that the solutions of the present method are in good agreement with FEM’s and analytical solutions.     

具有辐射边界的三维非规则域内稳态温度场分析

研究了具有辐射边界的空间非规则域内稳态导热问题，求解方法为在球极坐标系内分离变量，获得级数形式的解后，采用边界离散法确定级数项的待定系数，算例表明，边界离形方法不仅可以解决非正交边界问题，而且也可以处理诸如辐射边界的非线性边值问题。     

Weighted Least-Squares Collocation Method (WLSCM) for 2-D and 3-D Heat Conduction Problems in Irregular Domains

In this article, the weighted least-squares collocation method (WLSCM) is adopted to deal with two- and three-dimensional heat conduction problems in irregular domains. A radial basis function (RBF) is selected to construct the approximation function. To improve the accuracy and stability, some auxiliary points are increased within the domain of interest. Only inner nodes are used to construct the approximation function, and the equilibrium equations are satisfied not only at collocation points but also at auxiliary points, so the equations should be solved in a least-square sense. A 2-D case that has an analytical solution is simulated by the proposed method and the outcome verifies that the present method can obtain desired accuracy and efficiency. Then the current method is adopted to compute one 2-D and two 3-D cases of engineering heat conduction problems in irregular complex domains. The results show that the present method can deal effectively with the heat conduction problems of both 2-D and 3-D irregular domains.     

Local RBF meshless scheme for coupled radiative and conductive heat transfer

ABSTRACT

A local radial basis function meshless (LRBFM) method is developed to solve coupled radiative and conductive heat transfer problems in multidimensional participating media, in which compact support radial basis functions (RBFs) augmented on a polynomial basis are employed to construct the trial function, and the radiative transfer equation (RTE) and energy conservation equation are discretized directly at nodes by the collocation method. LRBFM belongs to a class of truly meshless methods which require no mesh or grid, and can be readily implemented in a set of uniform or irregular node distributions with no node connectivity. Performances of the LRBFM is compared to numerical results reported in the literature via a variety of coupled radiative and conductive heat transfer problems in 1D and 2D geometries. It is demonstrated that the local radial basis function meshless method provides high accuracy and great efficiency to solve coupled radiative and conductive heat transfer problems in multidimensional participating media with uniform and irregular node distribution, especially for coupled heat transfer problems in irregular geometry with Cartesian coordinates. In addition, it is extremely simple to implement.