A meshless procedure for transient heat conduction in functionally graded materials

A meshless numerical procedure is developed for analyzing the transient heat conduction problem in non-homogeneous functionally graded materials. In the proposed method the time derivative of temperature is approximate by the finite difference. At each time step the original nonlinear boundary value problem is transform into a hierarchy of non-homogeneous linear problem by used the homotopy analysis method. In this method a sought solution is presented by using a finite expansion in Taylor series, which consecutive elements are solutions of series linear value problems defining differential deformations. Each of linear boundary value problems with the corresponding boundary conditions is solved by using the method of fundamental solutions and radial basis functions which are used for interpolation of the inhomogeneous term. The accuracy of the obtained approximate solution is controlled by the number of components of the Taylor series, while the convergence of the process is monitored by an additional parameter of the method. Numerical experiments demonstrate the efficiency and accuracy of the present scheme in the solution of the heat conduction problem in nonlinear functionally graded materials. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Time-space boundary elements for transient heat conduction problems

This paper is concerned with a generalized time-space boundary element formulation for transient heat conduction problems in anisotropic media. A weighted residual form of the governing equation is used to obtain the boundary integral equation in terms of the fundamental solution. The resulting boundary integral equation is discretized by means of a wide variety of boundary elements from constant-elements to higher-order isoparametric elements located both in time and space.     

Efficient long-time computations of time-domain boundary integrals for 2D and dissipative wave equation

Linear hyperbolic partial differential equations in a homogeneous medium, e.g., the wave equation describing the propagation and scattering of acoustic waves, can be reformulated as time-domain boundary integral equations. We propose an efficient implementation of a numerical discretization of such equations when the strong Huygens’ principle does not hold.For the numerical discretization, we make use of convolution quadrature in time and standard Galerkin boundary element method in space. The quadrature in time results in a discrete convolution of weights W_j with the boundary density evaluated at equally spaced time points. If the strong Huygens’ principle holds, W_j converge to 0 exponentially quickly for large enough j. If the strong Huygens’ principle does not hold, e.g., in even space dimensions or when some damping is present, the weights are never zero, thereby presenting a difficulty for efficient numerical computation.In this paper we prove that the kernels of the convolution weights approximate in a certain sense the time domain fundamental solution and that the same holds if both are differentiated in space. The tails of the fundamental solution being very smooth, this implies that the tails of the weights are smooth and can efficiently be interpolated. Further, we hint on the possibility to apply the fast and oblivious convolution quadrature algorithm of Schädle et al. to further reduce memory requirements for long-time computation. We discuss the efficient implementation of the whole numerical scheme and present numerical experiments.     

Solution of a three-dimensional boundary-value problem for a fractional differential heat conduction equation

We consider a three-dimensional axisymmetric problem for a differential heat conduction equation with fractional time derivatives. Using the method of homogeneous solutions and integral transformations, we obtain an asymptotic and a numerical solution of the problem. The results of calculations are presented.     

移动边界半线性热传导方程始边值问题分析近似解

本文应用Fourier方法求得移动边界非齐次线性热传导方程始边值问题解及半线性方程问题分析近似解     

Coupled fluid-solid thermal interaction modeling for efficient transient simulation of biphasic water-steam energy systems

This paper proposes a fluid-solid coupled finite element formulation for the transient simulation of water-steam energy systems with phase change due to boiling and condensation. As it is commonly assumed in the study of thermal systems, the transient effects considered are exclusively originated by heat transfer processes. A homogeneous mixture model is adopted for the analysis of biphasic flow, resulting in a nonlinear transient advection-diffusion-reaction energy equation and an integral form for mass conservation in the fluid, coupled to the linear transient heat conduction equation for the solid. The conservation equations are approximated applying a stabilized Petrov-Galerkin FEM formulation, providing a set of coupled nonlinear equations for mass and energy conservation. This numerical model, combined with experimental heat transfer coefficients, provides a comprehensive simulation tool for the coupled analysis of boiling and condensation processes. For the treatment of enthalpy discontinuities traveling with the flow, a novel explicit-implicit time integration method based on Crank-Nicolson scheme is proposed, analyzing its accuracy and stability properties. To reduce problem size and enhance numerical efficiency, a modal superposition method with balanced truncation is applied to the solid equations. Finally, different example problems are solved to demonstrate the capabilities, flexibility and accuracy of the proposed formulation.     

A numerical technique for solving IHCPs using Tikhonov regularization method

This study is intended to provide a numerical algorithm for solving a one-dimensional inverse heat conduction problem. The given heat conduction equation, the boundary conditions, and the initial condition are presented in a dimensionless form. The numerical approach is developed based on the use of the solution to the auxiliary problem as a basis function. To regularize the resultant ill-conditioned linear system of equations, we apply the Tikhonov regularization method to obtain the stable numerical approximation to the solution.     

A BEM for transient thermoelastic fracture analysis of FGMs

A boundary element method for the transient thermoelastic fracture analysis in isotropic, continuously non-homogeneous and linear elastic functionally graded materials subjected to a thermal shock is presented. The material parameters are assumed to be continuous functions of the Cartesian coordinates. Laplace-domain fundamental solutions of linear coupled thermoelasticity for infinite, isotropic, homogeneous and linear elastic solids are applied to derive the boundary-domain integral equation formulation. The numerical implementation is performed by using a collocation method for the spatial discretization. Numerical results for the dynamic stress intensity factors are presented and discussed. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The method of differential quadrature for transient nonlinear diffusion

Numerical solutions to transient nonlinear diffusion problems are obtained by the method of differential quadrature. The accuracy of the solutions is inferred by comparison with the analytical solution for the linear case. The particular problems associated with general boundary conditions are handled by the use of integral methods. Examples from heat diffusion include composite media and radiation-enhanced conduction.     

Transient elastodynamic three-dimensional problems in cracked bodies

The general formulation of the transient elastodynamic second boundary value problem in an isotropic linear elastic body with a crack of arbitrary shape by combining the boundary integral equation method and the Laplace transform with respect to time is presented in this paper. Both finite and infinite elastic bodies are considered. A numerical solution of the transformed boundary integral equations is proposed.     

基于时域精细积分算法的瞬态传热多宗量反演

基于有限元法和精细积分算法,提出了一种求解瞬态热传导多宗量反演问题的新方法．采用有限元法和精细积分算法分别对空间、时间变量进行离散,可以得到正演问题高精度的半解析数值模型,由此建立了多宗量反演的计算模式,并给出敏度分析的计算公式．对一维和二维的热物性参数、热源项、边界条件等进行了单宗量和多宗量的反演求解,初步考虑了初值和噪音等对反演结果的影响,数值算例验证了该方法的有效性．     

On the coefficient inverse problem of heat conduction in a degenerating domain

ABSTRACT

In the paper, we consider a coefficient inverse problem for the heat equation in a degenerating angular domain. It has been shown that the inverse problem for the homogeneous heat equation with homogeneous boundary conditions has a nontrivial solution up to a constant factor consistent with the integral condition. Moreover, the solution of the considered inverse problem is found in explicit form. In conclusion, statements of possible generalizations and the results of numerical calculations are given.     

Boundary element method and wave equation

In this paper the boundary integral expression for a one-dimensional wave equation with homogeneous boundary conditions is developed. This is done using the time dependent fundamental solution of the corresponding hyperbolic partial differential equation. The integral expression developed is a generalized function with the same form as the well-known D'Alembert formula. The derivatives of the solution and some useful invariants on the characteristics of the partial differential equation are also calculated.The boundary element method is applied to find the numerical solution. The results show excellent agreement with analytical solutions.A multi-step procedure for large time steps which can be used in the boundary element method is also described.In addition, the way in which boundary conditions are introduced during the time dependent process is explained in detail. In the Appendix the main properties of Dirac's delta function and the Heaviside unit step function are described.     

Convergence of a numerical scheme of the type of the vortex loop method on a closed surface with an approximation to the surface shape

We consider a linear integral equation with a hypersingular integral treated in the sense of the Hadamard finite value. This equation arises when solving the Neumann boundary value problem for the Laplace equation with the use of the representation of the solution in the form of a double layer potential. We study the case in which an exterior or interior boundary value problem is solved in a domain whose boundary is a smooth closed surface and the integral equation is written out on that surface. For the numerical solution of the integral equation, the surface is approximated by spatial polygons whose vertices lie on the surface. We construct a numerical scheme for solving the integral equation on the basis of such an approximation to the surface with the use of quadrature formulas of the type of the method of discrete singularities with regularization. We prove that the numerical solutions converge to the exact solution of the hypersingular integral equation uniformly on the grid.     

Iterative solution of the generalized Love's problem in anisotropic media

An integral equation method is used to show the well posedness of the generalized Love's problem for an elastic anisotropic layer superimposed to a homogeneous substrate. The solution is derived using an iterative technique and the dispersion equation is obtained by imposing the pertinent boundary conditions. Inhomogeneous layers with different profiles of the material parameters are considered as examples and numerical results are given. The corresponding generalized problem is discussed for Lamb's modes.     

On the convergence of a numerical method for a hypersingular integral equation on a closed surface

We consider a linear integral equation with a hypersingular integral treated in the sense of the Hadamard finite value. This equation arises in the solution of the Neumann boundary value problem for the Laplace equation with a representation of a solution in the form of a double-layer potential. We consider the case in which the interior or exterior boundary value problem is solved in a domain; whose boundary is a smooth closed surface, and an integral equation is written out on that surface. For the integral operator in that equation, we suggest quadrature formulas like the method of vortical frames with a regularization, which provides its approximation on the entire surface for the use of a nonstructured partition. We construct a numerical scheme for the integral equation on the basis of suggested quadrature formulas, prove an estimate for the norm of the inverse matrix of the related system of linear equations and the uniform convergence of numerical solutions to the exact solution of the hypersingular integral equation on the grid.     

二维瞬态热传导的PDDO分析北大核心CSCD

采用近场动力学微分算子(peridynamic differential operator, PDDO)理论求解了二维瞬态热传导问题.将热传导方程和边界条件由其局部微分形式重构为非局部积分形式,引入Lagrange乘数法,采用变分原理的概念,建立了二维瞬态热传导问题的非局部分析模型.通过误差与收敛性分析,与其他数值方法计算结果进行比较,验证了本模型的准确性.在此基础上,将本模型应用于计算不规则边界板和内部含微缺陷(裂纹和圆孔)板的二维瞬态热传导问题.结果表明该方法计算精度高、适用范围广、具有较好的收敛性,为计算二维瞬态热传导问题提供了新的思路.     

非均匀介质GL型热弹性问题解的存在唯一性

We study a new model named the Green-Lindsay type therm-elastic model for nonhomogeneous media that consists of a system of dynamic thermoelasticity equations of displacement and dynamic heat conduction equation. We construct the model based on the classical GL-model for homogeneous material. This system is coupled dynamic problem and the displacement field and heat field must be solved at the same time. By using Fadeo-Galerkin method, we proved that the problem we proposed exist unique weak solution under some regular assumption.     

A numerical study of parabolic problems with nonlinear boundary conditions

In this paper, we consider an initial‐boundary value problem for a parabolic equation with nonlinear boundary conditions. The solution to the problem can be expressed as a convolution integral of a Green's function and two unknown functions. We change the problem to a system of two nonlinear Volterra integral equations of convolution type. By using an explicit procedure on the basis of Sinc‐function properties, the resulting integral equations are replaced by a system of nonlinear algebraic equations, whose solution yields an accurate approximate solution to the parabolic problem. Some examples are considered to illustrate the ability of the proposed method. Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical Solution of the Retrospective Inverse Problem of Heat Conduction with the Help of the Poisson Integral

We consider the retrospective inverse problem that consists in determining the initial solution of the one-dimensional heat conduction equation with a given condition at the final instant of time. The solution of the problem is given in the form of the Poisson integral and is numerically realized by means of a quadrature formula leading to a system of linear algebraic equations with dense matrix. The results of numerical experiments are presented and show the efficiency of the numerical method including the case of the final condition with random errors.