Thermoelastic state of a half-space having a thermally active crack parallel to its boundary

Using the boundary integral equation method, the problem of stationary heat conduction and thermoelasticity for a semi-infinite body with a crack parallel to its boundary is solved. Temperature or heat flow on the crack is prescribed. The body boundary is heat-insulated or is at zero temperature. The dependence of the stress intensity factor on the depth of occurrence of a circular crack at a constant temperature or under a constant heat flow is studied. In contrast to mechanical loading, thermal loading shows less SIF values than in an infinite body __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 4, pp. 46–54, April 2007.     

一种全耦合多相流分析的并行计算方法1)

研究了孔隙介质中热、水和汽流全耦合分析的并行计算方法.模型中采用了考虑毛细压力关系的修正有效应力概念,并考虑了相变和潜热传递.基本变量为位移、毛细压力、汽压和温度.并行程序是在国家高性能计算中心(北京)的曙光1000A上借助PVM(Parallel Virtual Machine)软件系统实现的,考题显示出较高的并行加速比和效率.     

COUPLING BETWEEN LAMINAR FILM CONDENSATION AND NATURAL CONVECTION ON OPPOSITE SIDES OF A VERTICAL PLATE

Theoretical analyses which incorporate one-dimensional heat conduction along a plate and transverse heat conduction approximations are presented to predict the net heat transfer between laminar film condensation of a saturated vapour on one side of a vertical plate and boundary layer natural convection on the other side. It is assumed that countercurrent boundary layer flows are formed on the two sides. The governing boundary layer equations of this problem and their corresponding boundary conditions are all cast into dimensionless forms by using a non-similarity transformation. Thus the resulting system of equations can be solved by using the local non-similarity method for the boundary layer equations and a finite difference method for the heat conduction equation of the plate. The plate temperature and the heat flux through the plate are repetitively determined until the solutions for each side of the plate match. The predicted results show that the effect of Pr_c is not negligible for larger values of A* (thermal resistance ratio between natural convecti on side and condensing film side) and the approximation of transverse heat conduction overpredicts the plate temperature for lower values of R_t (thermal resistance ratio between plate and condensing film). However, no significant differences are observed between the two different approximations for higher values of R_t. © by 1997 John Wiley & Sons, Ltd.     

Hamiltonian description of the heat conduction

It is shown that the linear boundary value problems of the heat conduction in a homogeneous slab can be mapped on the initial value problem for a Hamiltonian motion whose phase-space trajectories are subject to an additional restriction, the “arrival condition”. The physical consequences of this formal analogy for the macroscopic heat conduction are discussed in detail. Received on 27 April 1998     

Planar solidification of a warm flowing-liquid under different boundary conditions

Planar solidification of a warm flowing liquid with the convective heat transfer to the growing solid layer, has been analysed for the boundary conditions of constant temperature, constant heat flux and convective heat flux at the surface respectively. The mathematical formulation of the problem resulted in a coupled set of two differential equations in temperature and solid thickness as function of position, time and the problem parameters. Analytical expressions for the temperature distribution within the growing solid layer, the rate of solidification and the solidification time are obtained. The perturbation techniques employed here is simple and straight forward in contrast with the earlier techniques. Good agreement between the experimental results and the present solutions is obtained for the convective heat flux boundary condition. The results of this analysis are useful in the design and analysis of experiments dealing with freezing/melting in one dimension. The role of the parameter Stefan number which is small for phase change materials, is discussed in context with the storage of thermal energy.     

Solution of the classical single-phase Stefan problem in a modified formulation for semitransparent media

A numerical simulation of the classical Stefan problem is performed in a modified formulation for a semitransparent gray medium. Temperature fields and the resulting radiation flux field are obtained, and the dynamics of displacement of the interface and the evolution of the temperature rise on the black left boundary of the sample for different values of the emissivity factor of the exposed right boundary are studied.     

Axisymmetric non-Fourier temperature field in a hollow sphere

The non-Fourier axisymmetric (2+1)-dimensional temperature field within a hollow sphere is analytically investigated by the solution of the well-known Cattaneo–Vernotte hyperbolic heat conduction equation. The material is assumed to be homogeneous and isotropic with temperature-independent thermal properties. The method of solution is the standard separation of variables method. General linear time-independent boundary conditions are considered. Ultimately, the presented solution is applied to a (1+1)—as well as a (2+1)—dimensional problem, and their respective non-Fourier thermal behavior is studied. The present solution can be reduced to special cases of interest by choosing appropriate boundary conditions parameters. Dedicated to Prof. Gholamali Atefi, with appreciation and admiration on the occasion of his 65th birthday.     

Transient heat conduction in two-dimensional functionally graded hollow cylinder with finite length

In this paper a thick hollow cylinder with finite length made of two dimensional functionally graded material (2D-FGM) subjected to transient thermal boundary conditions is considered. The volume fraction distribution of materials, geometry and thermal boundary conditions are assumed to be axisymmetric but not uniform along the axial direction. The finite element method with graded material properties within each element is used to model the structure and the Crank–Nicolson finite difference method is implemented to solve time dependent equations of the heat transfer problem. Two-dimensional heat conduction in the cylinder is considered and variation of temperature with time as well as temperature distribution through the cylinder are investigated. Effects of variation of material distribution in two radial and axial directions on the temperature distribution and time response are studied. The achieved results show that using two-dimensional FGM leads to a more flexible design so that transient temperature, maximum amplitude and uniformity of temperature distributions can be modified to achieve required specifications by selecting a suitable material distribution profile in two directions.     

Numerical simulation of coupled heat and mass transfer in wood dried at high temperature

The mutual effect between heat and mass transfer is investigated for wood dried at high temperature. A numerical model of coupled heat and mass transfer under the effect of the pressure gradient is presented. Based on the macroscopic viewpoint of continuum mechanics, the mathematical model with three independent variables (temperature, moisture content and gas pressure) is constructed. Mass transfer in the pores involves a diffusional flow driven by the gradient of moisture content, convectional flow of gaseous mixture governed by the gradient of gas pressure, the Soret effect and phase change of water. Energy gain or loss due to phase change of water is taken as the heat source. Numerical methods, the finite element method and the finite difference method are used to discretize the spatial and time dimension, respectively. A direct iteration method to solve the nonlinear problem without direct evaluation of the tangential matrix is introduced. The local convergence condition based on the contraction–mapping principle is discussed. The mathematical model is applied to a 3-D wood board dried at high temperature with the Neumann boundary conditions for both temperature and moisture content, and the Dirichlet boundary conditions for gas pressure.     

一种全耦合多相流分析的并行计算方法

研究了孔隙介质中热、水和汽流全耦合分析的并行计算方法．模型中采用了考虑毛细压力关系的修正有效应力概念,并考虑了相变和潜热传递．基本变量为位移、毛细压力、汽压和温度．并行程序是在国家高性能计算中心（北京）的曙光１０００Ａ上借助ＰＶＭ（ＰａｒａｌｅｌＶｉｒｔｕａｌＭａｃｈｉｎｅ）软件系统实现的,考题显示出较高的并行加速比和效率     

Ablation of a hydrogen pellet in a light-gas gun bore

An ablation model for a hydrogen pellet accelerating in a thermonuclear facility fuel-injector bore is developed. The model is based on the Lagrange internal-ballistic problem complemented by the Stefan conditions on the ablating surface of the pellet Calculations demonstrate that, during pellet motion in the bore under the action of the gas, the pellet length remains almost unchanged, while due to the heat exchange with the bore walls the pellet diameter decreases by 4–10% St. Petersburg. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 3–11, January–February, 1999.     

A quasi diffusive reflective method for the calculation of radiative heat transfer in semitransparent media between parallel planes with angle dependent reflection behaviour

 The problem of radiation heat transfer in a semitransparent infinite plane parallel layer with purely diffusive reflection, emission and transmission behaviour (referred to as optical boundary conditions), in which the temperature depends only on the coordinate perpendicular to the boundaries, can be solved exactly by using exponential integrals of different orders. On the other hand the attenuation for the reflected, emitted and transmitted heat currents in the case of specular reflecting, emitting and transmitting boundaries has to be described by more complicate functions or the problem has to be treated by different methods as the discrete ordinates method (DOM). Within this work, a new method to describe radiation heat transfer for this geometry especially in the case of strongly angle dependent optical boundary conditions, including total reflection as described by Fresnel's equations, is developed. This method uses discrete angle ranges with averaged optical properties of the boundaries as well as the exponential integral formalism and is therefore referred to as quasi-diffusive reflective method (QDR-method). Exact solutions and solutions by means of the QDR-method including heat transfer by conduction are compared. The results of the QDR-method are sufficiently accurate in the most practical cases and allow calculations with a temperature dependent absorption coefficient in a much easier way. Received on 13 March 2000     

The analyticity of solutions of the Stefan problem

The Stefan problem of a semi-infinite body with arbitrarily prescribed initial and boundary conditions is studied. One of the objectives of the paper is to investigate the analyticity of the solutions. For this purpose, the prescribed initial and boundary conditions are considered to be series of fractional powers of their arguments. It is found that the exact solutions of the problem for various forms of the initial and boundary conditions can be established in series of parabolic cylinder functions and time t. Existence and convergence of the series solutions are studied and proved. The present solutions include the known exact solutions as special cases. On the basis of the present solutions, the question of the analyticity of solutions of the Stefan problem, raised by Rubinstein in his book, can be answered. Conditions for analyticity of the solutions with various initial and boundary conditions are fully discussed.     

Exact Solutions of the Problem of Unsteady Flow of a Viscoplastic Medium in a Circular Pipe

Two one-parameter series of real solutions describing the process of deceleration and acceleration of a viscoplastic medium under the action of a time-varying pressure gradient are obtained. The problem of axisymmetric unsteady viscoplastic flow is reduced to the solution of the Stefan boundary-value problem for the heat conduction equation with a nonlinear condition on the boundary of the quasi-rigid core. By a self-similar change of variables the problem can be reduced to a second-order ordinary differential equation. The solutions of this equation are represented in terms of Bessel and elementary functions. As a result, two one-parameter series of solutions, the first of which describes the acceleration and the second the deceleration of a viscoplastic medium in a pipe under the action of a time-varying pressure gradient are obtained.     

Thermal interaction between free convection and forced convection along a vertical conducting wall

A theoretical study is presented in this paper to investigate the conjugate heat transfer across a vertical finite wall separating two forced and free convection flows at different temperatures. It is assumed that the heat conduction in the wall is only in the transversal direction. We also assume that countercurrent boundary layers are formed on the both sides of the wall. The governing equations of this problem and their corresponding boundary conditions are all cast into a dimensionless form by using a non-similarity transformation. These resulting equations, which are singular at the points ξ_c=0 and 1, are solved numerically using a very efficient singular perturbation method. The effects of the resistance parameters and of the Prandtl numbers on heat transfer characteristics are investigated and presented in a table and ten figures. Received on 8 April 1998     

Freezing around a finite heat sink immersed in an infinite phase change medium

Freezing around a spherical heat sink immersed in an infinite phase change medium — a free boundary problem involving growth and decay of the free boundary — is analysed here. A one-dimensional conduction model is formulated and the resulting partial differential equations are solved by finite difference methods. The energy discharged from the phase change medium during the heat transfer process is analysed for latent heat thermal energy storage applications. Results are presented for a wide range of parameters that are encountered in energy storage devices. The cases of slab/cylindrical heat sink are reexamined for a range of parameters not covered by the earlier investigators     

The effects of casting speed on steel continuous casting process

A three dimensional simulation of molten steel flow, heat transfer and solidification in mold and “secondary cooling zone” of Continuous Casting machine was performed with consideration of standard k−ε model. For this purpose, computational fluid dynamics software, FLUENT was utilized. From the simulation standpoint, the main distinction between this work and preceding ones is that, the phase change process (solidification) and flow (turbulent in mold section and laminar in secondary cooling zone) have been coupled and solved jointly instead of dividing it into “transient heat conduction” and “steady fluid flow” that can lead to more realistic simulation. Determining the appropriate boundary conditions in secondary cooling zone is very complicated because of various forms of heat transfer involved, including natural and forced convection and simultaneous radiation heat transfer. The main objective of this work is to have better understanding of heat transfer and solidification in the continuous casting process. Also, effects of casting speed on heat flux and shell thickness and role of radiation in total heat transfer is discussed.     

An ALE finite element method for a coupled Stefan problem and Navier–Stokes equations with free capillary surface

In this article, an ALE finite element method to simulate the partial melting of a workpiece of metal is presented. The model includes the heat transport in both the solid and liquid part, fluid flow in the liquid phase by the Navier–Stokes equations, tracking of the melt interface solid/liquid by the Stefan condition, treatment of the capillary boundary accounting for surface tension effects and a radiative boundary condition. We show that an accurate treatment of the moving boundaries is crucial to resolve their respective influences on the flow field and thus on the overall energy transport correctly. This is achieved by a mesh‐moving method, which explicitly tracks the phase boundary and makes it possible to use a sharp interface model without singularities in the boundary conditions at the triple junction. A numerical example describing the welding of a thin‐steel wire end by a laser, where all aforementioned effects have to be taken into account, proves the effectiveness of the approach.Copyright © 2012 John Wiley & Sons, Ltd.     

High Péclet number heat exchange between cocurrent streams

Summary  This paper concentrates on the analysis of the heat transfer between two cocurrent laminar flows in parallel channels. For high values of the Péclet number Pe, a boundary layer arises near the wall separating the streams. Matched asymptotic expansions (MAE) are used to obtain approximate solutions. We consider arbitrary inlet temperatures and derive higher-order corrections of the boundary problem. The separating wall is supposed to be sufficiently thin to neglect the heat conduction in it. Analyticity and adiabatic conditions at the outer walls impose restrictions on the inlet temperatures. It turns out, however, that only the inlet temperatures at the wall separating the two fluids enter the leading-order problem. The Nusselt numbers thus calculated are in the leading order proportional to (Pe/x)^1/3, where x is the stream-wise coordinate. An estimate of the thickness of the separating wall to validate the MAE approach is obtained. It is also demonstrated that the MAE analysis is unable to describe the heat exchange of counterflowing fluids. Received 9 June 1999; accepted for publication 17 November 1999