非线性二维导热反问题的混沌-正则化混合解法

考虑热传导系数随温度变化,建立了非线性二维稳态导热反问题数值计算模型。并把混沌优化方法和梯度正则化方法相结合,构成一种混沌-正则化混合算法求该计算模型的全局解。以热传导系数随温度线性变化为例,由布置在结构边界上的观测点温度信息确定了结构材料热传导系数及其随温度变化规律。结果表明混合算法计算结果与初值无关,具有很好的全局寻优性能,而且计算量远比经典遗传算法和单纯采用混沌优化方法小。     

The estimation of the strength of the heat source in the heat conduction problems

A general method is proposed to determine the strength of the heat source in the Fourier and non-Fourier heat conduction problems. A finite difference method, the concept of the future time and a modified Newton–Raphson method are adopted in the problem. The undetermined heat source at each time step is formulated as an unknown variable in a set of equations from the measured temperature and the calculated temperature. Then, an iterative process is used to solve the set of equations. No selected function is needed to represent the undetermined function in advance. Three examples are used to demonstrate the characteristics of the proposed method. The validity of the proposed method is confirmed by the numerical results. The results show that the proposed method is an accurate and stable method to determine the strength of the heat source in the inverse hyperbolic heat conduction problems. Furthermore, the result shows that more future times are needed in the hyperbolic equation than that of parabolic equation. Moreover, the robustness and the accuracy of the estimated results in the non-Fourier problem are not as well as those of the Fourier problem.     

基于分数阶热传导方程激光加热瞬态温度场研究

基于分数阶Taylor(泰勒)级数展开原理,建立单相延迟一阶分数阶近似方程,获得分数阶热传导方程.针对短脉冲激光加热问题建立分数阶热传导方程组,并运用Laplace(拉普拉斯)变换方法进行求解,给出非Gauss(高斯)时间分布的激光内热源温度场解析解.针对具体算例数值研究温度波传播特性.结果表明热传播速度与分数阶阶次有关,分数阶阶次增加,热传播速度减小,温度变化幅度增加.分数阶方程可以用于描述介于扩散方程和热波方程间的热传输过程,且对热传播机制与分数阶热传导方程中分数阶项的关系做了深入剖析.     

Inverse thermal imaging in materials with nonlinear conductivity by material and shape derivative method

The material and shape derivative method is used for an inverse problem in thermal imaging. The goal is to identify the boundary of unknown inclusions inside an object by applying a heat source and measuring the induced temperature near the boundary of the sample. The problem is studied in the framework of quasilinear elliptic equations. The explicit form is derived of the equations that are satisfied by material and shape derivatives. The existence of weak material derivative is proved. These general findings are demonstrated on the steepest descent optimization procedure. Simulations involving the level set method for tracing the interface are performed for several materials with nonlinear heat conductivity. Copyright © 2011 John Wiley & Sons, Ltd.     

An approximate solution of a Fredholm integral equation of the first kind by the residual method

A Tikhonov finite-dimensional approximation is applied to a Fredholm integral equation of the first kind. This allows using a variational regularization method with a regularization parameter from the residual principle and reducing the problem to a system of linear algebraic equations. The accuracy of the approximate solution is estimated with allowance for the error of the finitedimensional approximation of the problem. The use of this approach is illustrated by solving an inverse boundary value problem for the heat conduction equation.     

Determination of a time-dependent diffusivity from nonlocal conditions

The inverse problem of determining the temperature and the time-dependent thermal diffusivity from various additional nonlocal information is investigated. These nonlocal conditions can come in the form of an internal or boundary energy, or, in the one-dimensional case, as a difference boundary temperature or heat flux so as to ensure the uniqueness of solution for the heat conduction equation with unknown thermal diffusivity coefficient. The Ritz-Galerkin method with satisfier function is employed to solve the inverse problems numerically. Numerical results are presented and discussed.     

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.     

Heat transfer in MHD viscoelastic fluid flow over a stretching sheet with variable thermal conductivity,non-uniform heat source and radiation

An analysis has been carried out to study the magnetohydrodynamic boundary layer flow and heat transfer characteristics of a non-Newtonian viscoelastic fluid over a flat sheet with a linear velocity in the presence of thermal radiation and non-uniform heat source. The thermal conductivity is assumed to vary as a linear function of temperature. The basic equations governing the flow and heat transfer are in the form of partial differential equations, the same have been reduced to a set of non-linear ordinary differential equations by applying suitable similarity transformation. The transformed equations are solved analytically by regular perturbation method. Numerical solution of the problem is also obtained by the efficient shooting method, which agrees well with the analytical solution. The effects of various physical parameters such as viscoelastic parameter, Chandrasekhar number, Prandtl number, variable thermal conductivity parameter, Eckert number, thermal radiation parameter and non-uniform heat source/sink parameters which determine the temperature profiles are shown in several plots and the heat transfer coefficient is tabulated for a range of values of said parameters. Some important findings reported in this work reveals that combined effect of variable thermal conductivity, radiation and non-uniform heat source have significant impact in controlling the rate of heat transfer in the boundary layer region.     

Effects of variable viscosity on non-Darcy MHD free convection along a non-isothermal vertical surface in a thermally stratified porous medium

An analysis is performed for non-Darcy free convection flow of an electrically conducting fluid over an impermeable vertical plate embedded in a thermally stratified, fluid saturated porous medium for the case of power-law surface temperature. The present work examines the effects of non-Darcian flow phenomena, variable viscosity, Hartmann–Darcy number and thermal stratification on free convective transport and demonstrates the variation in heat transfer prediction based on three different flow models. The wall effect on porosity variation is approximated by an exponential function. The effects of thermal dispersion and variable stagnant thermal conductivity are taken into consideration in the energy equation. The resulting non-similar system of equations is solved using a finite difference method. Results are presented for velocity, temperature profiles and local Nusselt number for representative values of different controlling parameters.     

Effects of radiation and variable viscosity on unsteady MHD flow of a rotating fluid from stretching surface in porous medium

This work is focused on the study of unsteady magnetohydrodynamics boundary-layer flow and heat transfer for a viscous laminar incompressible electrically conducting and rotating fluid due to a stretching surface embedded in a saturated porous medium with a temperature-dependent viscosity in the presence of a magnetic field and thermal radiation effects. The fluid viscosity is assumed to vary as an inverse linear function of temperature. The Rosseland diffusion approximation is used to describe the radiative heat flux in the energy equation. With appropriate transformations, the unsteady MHD boundary layer equations are reduced to local nonsimilarity equations. Numerical solutions of these equations are obtained by using the Runge–Kutta integration scheme as well as the local nonsimilarity method with second order truncation. Comparisons with previously published work have been conducted and the results are found to be in excellent agreement. A parametric study of the physical parameters is conducted and a representative set of numerical results for the velocity in primary and secondary flows as well as the local skin-friction coefficients and the local Nusselt number are illustrated graphically to show interesting features of Darcy number, viscosity-variation, magnetic field, rotation of the fluid, and conduction radiation parameters.     

非常量导热系数几例微分方程的推立

在本文里,曾先后假设物体的导热系数是依直线和指数函数空间地起改变,就这样来建立了六个二阶热传导微分方程;又对于变密度、变比热、变导热系数这样的更一般的情况也推立了六个二阶热传导的微分方程.     

A superstatistical model for anomalous heat conduction and diffusion

We propose a superstatistical model for anomalous heat conduction and diffusion, which is formulated by the thermal conductivity distribution, overall temperature and heat flux distributions. Our model obeys Fourier's law and the continuity equation at the individual level. The evolution of the thermal conductivity distribution is described by an advection-diffusion equation. We show that the superstatistical model predict anomalous behaviors including the time-dependent effective thermal conductivity and slow long-time asymptotics. The time-dependence of the effective thermal conductivity is determined by the mean square displacement (MSD), which coincides with existing investigations. The superstatistical structure can also be extended into other non-Fourier models including the Cattaneo and fractional-order heat conduction models.     

Direct and inverse solutions of the two-dimensional hyperbolic heat conduction problems

A sequential method is proposed to estimate boundary condition of the two-dimensional hyperbolic heat conduction problems. An inverse solution is deduced from a finite difference method, the concept of the future time and a modified Newton–Raphson method. The undetermined boundary condition at each time step is denoted as an unknown variable in a set of non-linear equations, which are formulated from the measured temperature and the calculated temperature. Then, an iterative process is used to solve the set of equations. No selected function is needed to represent the undetermined function in advance. The example problem is used to demonstrate the characteristics of the proposed method. In the example, a well-known problem is used to demonstrate the validity of the proposed direct method and then the inverse solutions are evaluated. In the second example, the larger value of the relaxation time is implemented in the direct solutions and the inverse solutions. The close agreement between the exact values and the estimated results is made to confirm the validity and accuracy of the proposed method. The results show that the proposed method is an accurate and stable method to determine the boundary conditions in the two-dimensional inverse hyperbolic heat conduction problems.     

Solving evolutionary-type differential equations and physical problems using the operator method

We present a general operator method based on the advanced technique of the inverse derivative operator for solving a wide range of problems described by some classes of differential equations. We construct and use inverse differential operators to solve several differential equations. We obtain operator identities involving an inverse derivative operator, integral transformations, and generalized forms of orthogonal polynomials and special functions. We present examples of using the operator method to construct solutions of equations containing linear and quadratic forms of a pair of operators satisfying Heisenberg-type relations and solutions of various modifications of partial differential equations of the Fourier heat conduction type, Fokker–Planck type, Black–Scholes type, etc. We demonstrate using the operator technique to solve several physical problems related to the charge motion in quantum mechanics, heat propagation, and the dynamics of the beams in accelerators.     

三维热传导型半导体问题的交替方向有限元方法

1　引　　言三维热传导型半导体器件瞬态问题的数学模型由四个非线性偏微分方程描述[1 ,2 ] ,记 Ω为 Ω=[0 ,1 ] 3的边界 ,三维问题-Δψ =α( p -e+ N( x) ) ,　　 ( x,t)∈Ω× [0 ,T] ,( 1 .1 ) e t= . ( De( x) e-μe( x) e ψ) -R( e,p,T) ,　 ( x,t)∈Ω× ( 0 ,T] ,( 1 .2 ) p t= . ( Dp( x) p +μp( x) p ψ) -R( e,p,T) ,　 ( x,t)∈Ω× ( 0 ,T] ,( 1 .3 )ρ( x) T t-ΔT =[( Dp( x) p +μp( x) p ψ) -( De( x) e-μe( x) e ψ) ] . ψ,　　　　　　 ( x,t)∈Ω× ( 0 ,T] . ( 1 .4 )ψ( x,t) =e( x,t) =p( …     

The effect of variable viscosity on MHD viscoelastic fluid flow and heat transfer over a stretching sheet

An analysis has been carried out to study the momentum and heat transfer characteristics in an incompressible electrically conducting non-Newtonian boundary layer flow of a viscoelastic fluid over a stretching sheet. The partial differential equations governing the flow and heat transfer characteristics are converted into highly non-linear coupled ordinary differential equations by similarity transformations. The effect of variable fluid viscosity, Magnetic parameter, Prandtl number, variable thermal conductivity, heat source/sink parameter and thermal radiation parameter are analyzed for velocity, temperature fields, and wall temperature gradient. The resultant coupled highly non-linear ordinary differential equations are solved numerically by employing a shooting technique with fourth order Runge–Kutta integration scheme. The fluid viscosity and thermal conductivity, respectively, assumed to vary as an inverse and linear function of temperature. The analysis reveals that the wall temperature profile decreases significantly due to increase in magnetic field parameter. Further, it is noticed that the skin friction of the sheet decreases due to increase in the Magnetic parameter of the flow characteristics.     

An inverse problem in estimating simultaneously the effective thermal conductivity and volumetric heat capacity of biological tissue

An inverse problem utilizing the Levenberg–Marquardt method (LMM) is applied in this study to determine simultaneously the unknown spatial-dependent effective thermal conductivity and volumetric heat capacity for a biological tissue based on temperature measurements. The accuracy of this inverse problem is examined by using the simulated exact and inexact temperature measurements in the numerical experiments. A statistical analysis is performed to obtain the 99% confidence bounds for the estimated thermal properties. Results show that good estimation on the spatial-dependent thermal conductivity and volumetric heat capacity can be obtained using the present algorithm for the test cases considered in this study.     

三维热传导型半导体问题的特征有限元方法和分析

本文研究三维热传导型半导体瞬态问题的特征有限元方法及其理论分析,其数学模型是一类非线性偏微分方程的初边值问题,对电子位势方程提出Ｇａｌｅｒｋｉｎ逼近;对电子,空穴浓度方程采用特征有限元逼近;对热传导方程采用对时间向后差分的Ｇａｌｅｒｋｉｎ逼近．应用微分方程先验估计理论和技巧得到了最优阶Ｌ＾２误差估计。     

Analytical solutions to nonlinear equations arising in heat transfer by variational iteration,homotopy perturbation,and Adomian decomposition methods

As thermal conductivity plays an important role on fin efficiency, we tried to solve heat transfer equation with thermal conductivity as a function of temperature. In this research, some new analytical methods called homotopy perturbation method, variational iteration method, and Adomian decomposition method are introduced to be applied to solve the nonlinear heat transfer equations, and also the comparison of the applied methods (together) is shown graphically. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

在模糊随机因素影响下的多层圆筒结构的稳态热传导问题的区间数解法

在多层圆筒结构稳态热传导分析中,根据给定固体壁两侧表面温度总传热量公式,首先推导出当边界温度为随机变量情况下总传热量函数统计参数的均值和方差;然后推导出在导热系数为模糊数,边界温度为随机数下的总传热量的区间表达式．通过比较可以知道由区间数算法得到的区间最大,由概率统计算法得到的区间最小．并给出了两者的相对误差公式．最后引进粗糙集中的上、下近似集,提出用一个参数来统一定义模糊和随机区间进行稳态结构的热传导分析．