A space marching method for multidimensional transient inverse heat conduction problems

A new method for the solution of multidimensional heat conduction problems is formulated. The developed space marching method allows to determine quickly and exactly unsteady temperature distributions in the construction elements of irregular geometry. The method is especially appropriate for determining transient temperature distribution in thick-wall pressure components based on temperature measurements at the outer surface. Two examples are included to demonstrate the capabilities of the new approach. Received on 28 October 1998     

Experimental verification of space marching methods for solving inverse heat conduction problems

Two space marching methods for solving one-dimensional nonlinear inverse heat conduction problems developed earlier are validated using data obtained during emergency cooling simulation of PWR pressure vessels. The model of pressure vessel was preheated at 350 ^∘C and then cooled by injection of cold water at 20 ^∘C. Based on the vessel temperature histories measured at the distance of 50 mm from the inner surface and at the outside surface the transient temperature distribution in the pressure vessel was obtained. The comparison of the results obtained by both space marching methods with experimental data shows high accuracy of the space marching methods. Received on 3 May 1999     

New analytic solutions to 2D transient heat conduction problems with/without heat sources in the symplectic space

The two-dimensional (2D) transient heat conduction problems with/without heat sources in a rectangular domain under different combinations of temperature and heat flux boundary conditions are studied by a novel symplectic superposition method (SSM). The solution process is within the Hamiltonian system framework such that the mathematical procedures in the symplectic space can be implemented, which provides an exceptional direct rigorous derivation without any assumptions or predetermination of the solution forms compared with the conventional inverse/semi-inverse methods. The distinctive advantage of the SSM offers an access to new analytic heat conduction solutions. The results obtained by the SSM agree well with those obtained from the finite element method (FEM), which confirms the accuracy of the SSM.     

PRECISE INTEGRAL ALGORITHM BASED SOLUTION FOR TRANSIENT INVERSE HEAT CONDUCTION PROBLEMS WITH MULTI-VARIABLES

IntroductionIHCPs (InverseHeatConductionProblems)arecloselyassociatedwithmanyengineeringaspects,andwelldocumentedintheliteratures,coveringtheidentificationsofthermalparameters[1,2 ],boundaryshapes[3],boundaryconditions[4 ]andsource_relatedterms[5 ,6 ]etc .Howeveritseemsthatonlylittleworkisdirectlyconcernedwithmulti_variablesidentificationsbyauthors’knowledge.Tsengetal.presentedanapproachtodeterminingtwokindsofvariables[7],butonlygavefewnumericalexamplestodeterminethemsimultaneously .Thesol…     

Green’s function solution for transient heat conduction in annular fin during solidification of phase change material

The Green’s function method is applied for the transient temperature of an annular fin when a phase change material (PCM) solidifies on it. The solidification of the PCMs takes place in a cylindrical shell storage. The thickness of the solid PCM on the fin varies with time and is obtained by the Megerlin method. The models are found with the Bessel equation to form an analytical solution. Three different kinds of boundary conditions are investigated. The comparison between analytical and numerical solutions is given. The results demonstrate that the significant accuracy is obtained for the temperature distribution for the fin in all cases.     

Asymptotic averaging technique for heat conduction problems with phase transitions in layered media

Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 36, No. 5, pp. 155–163, September–October, 1995.     

Immersed boundary conditions method for heat conduction problems in slots with time‐dependent geometry

A grid‐less, fully implicit, spectrally accurate algorithm for solving three‐dimensional, both stationary and time‐dependent, heat conduction problems in slots formed by either fixed or time‐dependent boundaries has been developed. The algorithm is based on the concept of immersed boundary conditions (IBC), where the physical domain is immersed within the computational domain and the boundary conditions take the form of internal constraints. The IBC method avoids the need to construct adaptive, time‐dependent grids resulting in the reduction of the required computational resources and, at the same time, maintaining accurate information about the location of the boundaries. The algorithm is spectrally accurate in space and capable of delivering first‐, second‐, third‐ and fourth‐order accuracy in time. Given a potentially large size of the resultant linear algebraic system, various methods that take advantage of the special structure of the coefficient matrix have been explored in search for an efficient solver, including a specialized direct solver as well as serial and parallel iterative solvers. The specialized direct solver has been found to be the most efficient from the viewpoints of the speed of the computations and the memory requirements. Copyright © 2010 John Wiley & Sons, Ltd.     

Comparison of four different versions of the variable metric method for solving inverse heat conduction problems

In this study, four different versions of the variable metric method (VMM) are investigated in solving standard one-dimensional inverse heat conduction problems in order to evaluate their efficiency and accuracy. These versions include Davidon–Fletcher–Powell (DFP), Broydon–Fletcher–Goldfarb–Shanno (BFGS), Symmetric Rank-one (SR1), and Biggs formula of the VMM. These investigations are carried out using temperature data obtained from numerical simulations.     

Non-linear boundary and eigenvalue problems for the emden-fowler equations by group methods

Two pragmatic boundary value and eigenvalue problems of the Emden-Fowler equation $\text{(t}{}^{\mathrm{\alpha }}\text{u′)′ + λt}{}^{\mathrm{\beta }}\text{?(u) = 0,?(u) = u}{}^{\mathrm{\gamma }}$ and e^u are studied using the simple one parameter group properties. In all cases boundary value problems are converted into initial value problems using the property of the invariance group. With $\text{?(u) = u}{}^{\mathrm{\gamma }}$ an eigenvalue problem is detailed and calculations presented.     

On the boundary value problems for ordinary differential equations with turning points

In this paper,we consider the boundary value problems of the formsy″-f(x,ε)y′ g(x,ε)=0 (-a≤x≤b,0≤ε《1 )y(-a)=a,y(b)=βwhere f(x,0)has several and multiple zeros on the interval[-a,b].The conditions forexhibiting boundary and interior layers are given,and the corresponding asymptoticexpansions of solutions are constructed.     

Estimate of the transient conduction of heat in materials with linear thermal properties based on the solution for constant properties

A method of analysis is described which yields quasianalytical solutions for one and multidimensional unsteady heat conduction problems with linearly dependent thermal properties, such as thermal conductivity and volumetric specific heat. The method accomodates rather general thermal boundary conditions including arbitrary variations in surface temperature or in surface heat flux or a convective exchange with a fluid having even varying temperature. Once the solution for the identical problem but with constant properties has been developed, its practical realization is rather direct, being facilitated by a reduced number of iterations. The four applied examples given in this work show that a wide variety of nonlinear heat conduction problems can be tackled by this procedure without much difficulty. These simple solutions compare favorably with more laborious results reported in the archival heat transfer literature.

---

Berechnung nichtstationärer Wärmeleitvorgänge mit linear temperaturabhängigen Stoffwerten aus der Lösung für konstante Stoffwerte

Zusammenfassung Es werden quasi-analytische Lösungen für ein- und mehrdimensionale nichtstationäre Wärmeleitprobleme mit linear temperaturabhängigen Stoffwerten, wie Wärmeleitfähigkeit und volumetrische Wärmekapazität, mitgeteilt. Die Methode gilt für recht allgemeine Randbedingungen wie beliebige Veränderungen der Oberflächentemperatur, der Wärmestromdichte oder auch konvektiven Wärmeaustausch mit veränderlicher Fluidtemperatur. Ist die Lösung für das identische Problem mit konstanten Stoffwerten bekannt, kann die Methode direkt mit einer begrenzten Zahl von Iterationen angewandt werden. Die vier hier mitgeteilten Beispiele zeigen, daß eine große Zahl nichtlinearer Wärmeleitprobleme auf diese Weise ohne Schwierigkeit angepackt werden können. Die einfachen Lösungen stimmen befriedigend mit komplizierteren Ergebnissen aus der Literatur überein.

---

Nomenclature a side of square bar - B _i0 reference Biot number,hR/k₀ - B _i0 ^T transformed Biot number, equation (16) - c geometric parameter, equation (8) - h convective coefficient - k thermal conductivity - k ₀ value ofk atT ₀ - K dimensionless thermal conductivity,k/k ₀ - K _i value ofK at _i - K _i+1 value ofK at _i+1 - m _k slope of theK- line, equation (3) - m _s slope of theS- line, equation (4) - R characteristic length - s volumetric specific heat - s ₀ value of s at T₀ - S dimensionless volumetric specific heat, s/s₀ - S _i value ofS at _i - S _i+1 value of S at _i+1 - t time - T temperature - T ₀ reference temperature - x, y cartesian coordinates - X, Y dimensionless cartesian coordinates,x/a andy/a - thermal diffusivity - _k transformed time, equation (11) - _s transformed time, equation (37) - _k dimensionless time for variable conductivity, equation (8) - _s dimensionless time for variable specific heat, equation (34) - dimensionless temperature,T/T ₀ - dimensionless coordinate,r/R - ₀ value of at T₀ - _i lower value of the interval (_i, _i+1) - _i+1 upper value of the interval (_i, _i+1     

The numerical method for analysing the transient temperature field and transient phase transformation of metal materials in heat treating

In this paper,the Kirchhoffs transformation is popularized to the nonlinear heat conduction problem which the heat conductivity can be expressd as a multinomial of temperature firstly,the boundary condition of heat conduction problem is determined by analytics.Secondly,the incubation peroid superposition and the linear combination law is employed to simulate the transient phasses transformation in the process of heat treatment of materials.That the begin time of phase transformation,the type of phase transformation and the amount of phase constitution is determined simply.Finally,the three-dimension Dual Reciprocity Boundary Element Method is usedto analysis the total process of various heat treatment of component,the results of numerical calculation of examples show that the method provided in this paper is effectivce.     

On the heat conduction problem and a design principle for plasma arc furnace

In this paper a plane heat conduction problem with variable coefficients of heat conductivity K(T) is analysed with given electric power supplied to the plasma arc. The governing equation for unknown temperature distribution is a nonlinear one with a function as its nonhomogeneous term. To make the problem attractable by the method of separation of variables, a set of transformation of governing equation is introduced. An explicit simple formula is found for the efficiency of the furnace , depends linearly on r_o, the nondimensional distance between the arc and surface of melted material, as well as on another nondimensional quantity Q, we describe the above in detail in this paper. This relationship holds for r_o<0.4 and gives a good guidance for the design of furnace.Institute of Mechanics, Academia Sinica     

Modeling for heat and mass transfer with phase change in porous wick of CPL evaporator

A mathematic model is developed to describe heat and mass transfer with phase change in the porous wick of evaporator of capillary pumped loop (CPL). This model with six field variables, including temperature, liquid content, pressure, liquid velocity, vapor velocity and phase-change rate, is closed mathematically with additional pressure relationships introduced. The present model is suitable to the numerical computation, as the established equations become comparatively easy to solve, which is applied to CPL evaporator. The numerical results are obtained and the parameter effects on evaporator are discussed. The study demonstrates that instead of an evaporative interface, there exists an unsaturated two-phase zone between the vapor-saturated zone and the liquid-saturated zone in the wick of CPL evaporator.     

Multiple spatial and temporal scales method for numerical simulation of non-classical heat conduction problems: one dimensional case

A multiple spatial and temporal scales method is studied to simulate numerically the phenomenon of non-Fourier heat conduction in periodic heterogeneous materials. The model developed is based on the higher-order homogenization theory with multiple spatial and temporal scales in one dimensional case. The amplified spatial scale and the reduced temporal scale are introduced respectively to account for the fluctuations of non-Fourier heat conduction due to material heterogeneity and non-local effect of the homogenized solution. By separating the governing equations into various scales, the different orders of homogenized non-Fourier heat conduction equations are obtained. The reduced time dependence is thus eliminated and the fourth-order governing differential equations are derived. To avoid the necessity of C¹ continuous finite element implementation, a C⁰ continuous mixed finite element approximation scheme is put forward. Numerical results are shown to demonstrate the efficiency and validity of the proposed method.     

On the boundary value problems for a class of ordinary differential equations with turning points

In this paper we study the boundary value problems for a class of ordinary differential equations with turning points by the method of multiple scales. The paradox in [1] and the variational approach in [2] are avoided. The uniformly valid asymptotic approximations of solutions have been constructed. We also study the case which does not exhibit resonance.     

Effects of nonmonotonicity and change in sign of the energy flux in the transient regime for the Couette problem with heat transfer

On the basis of the numerical solution of the Boltzmann equation the Couette flow is investigated. The effect of the upper plate velocity on the value and position of the extremum of the dependence of the energy flux to the lower plate on the plate temperature ratio in the free-molecular regime is considered. It is established that in the case of sufficiently small surface temperature ratio the dependence of the energy flux to the lower plate on the Knudsen number Kn has a local extremum. It is found that there is a plate temperature ratio range on which with increase in Kn the energy flux to the lower plate reverses direction. The Kn dependence of its absolute value then has a local maximum.