Comparison between exact thermal boundary condition and harmonic mean conductivity condition at the solid–fluid interface for finite thickness shrouded non-isothermal fin array

A numerical study is made comparing the exact thermal boundary condition and a harmonic mean conductivity condition at the solid–fluid interface for a finite thickness shrouded non-isothermal fin array. Results highlight that there exists a significant deviation of pressure drop across the length of the fin for the exact thermal boundary condition, which is as high as 20% as compared to that obtained using the harmonic mean conductivity condition. The exact thermal boundary condition forecasts relatively more non-isothermal fin as compared to a harmonic mean conductivity condition. The greater the fin spacing the larger the non-isothermal behavior of the fin and it also depends upon Grashof number as well as inlet fluid velocity. The larger the Grashof number the greater is the non-isothermal behavior of fin. The greater the inlet velocity, the larger is the non-isothermal behavior of fin. Bulk fluid temperature is over predicted by as much as 13% by the harmonic mean conductivity condition for larger fin spacing with highest Grashof number coupled with larger velocity. This deviation is only 6% for smaller fin spacing. Overall Nusselt number is over predicted for the harmonic mean conductivity condition as compared to exact thermal boundary condition. This over prediction is limited to about 8%.     

Solving two-phase freezing Stefan problems: Stability and monotonicity

The two-phase Stefan problems with phase formation and depletion are special cases of moving boundary problems with interest in science and industry. In this work, we study a solidification problem, introducing a front-fixing transformation. The resulting non-linear partial differential system involves singularities, both at the beginning of the freezing process and when the depletion is complete, that are treated with special attention in the numerical modelling. The problem is decomposed in three stages, in which implicit and explicit finite difference schemes are used. Numerical analysis reveals qualitative properties of the numerical solution spatial monotonicity of both solid and liquid temperatures and the evolution of the solidification front. Numerical experiments illustrate the behaviour of the temperatures profiles with time, as well as the dynamics of the solidification front.     

Heat Conduction and Finite Measures for Transition Times between Steady States

A finite time, the mean action time associated with the conductivetransition from a constant initial temperature to thermal equilibriumat a constant ambient temperature, is related to time lag constants,mean energy residence times, mean first passage times, and Green'sfunction properties for linear equations, and to freezing timesand finite transition times for problems with phase transitions.When the boundary conditions are linear of mixed type and theconductivity is constant, or when they are of Dirichlet typeand the conductivity perhaps temperature dependent, the meanaction time is given by the solution of a linear Poisson problem.It is then easily found and is a useful finite comparative timefor the thermal transition process, giving a measure of itsdependence on size and other geometric factors.     

Free vibration analysis of FGM cylindrical shells surrounded by Pasternak elastic foundation in thermal environment considering fluid-structure interaction

This paper presents an investigation on partially fluid-filled cylindrical shells made of functionally graded materials (FGM) surrounded by elastic foundations (Pasternak elastic foundation) in thermal environment. Material properties are assumed to be temperature dependent and radially variable in terms of volume fraction of ceramic and metal according to a simple power law distribution. The shells are reinforced by stiffeners attached to their inside and outside in which the material properties of shell and the stiffeners are assumed to be continuously graded in the thickness direction. The formulations are derived based on smeared stiffeners technique and classical shell theory using higher-order shear deformation theory which accounts for shear flexibility through shell's thickness. Displacements and rotations of the shell middle surface are approximated by combining polynomial functions in the meridian direction and truncated Fourier series with an appropriate number of harmonic terms in the circumferential direction. The governing equations of liquid motion are derived using a finite strip element formulation of incompressible inviscid potential flow. The dynamic pressure of the fluid is expanded as a power series in the radial direction. Moreover, the quiescent liquid free surface is modeled by concentric annular rings. A detailed numerical study is carried out to investigate the effects of power-law index of functional graded material, fluid depth, stiffeners, boundary conditions, temperature and geometry of the shell on the natural frequency of eccentrically stiffened functionally graded shell surrounded by Pasternak foundations.     

Numerical simulation of thermal problems coupled with magnetohydrodynamic effects in aluminium cell

A system of partial differential equations describing the thermal behavior of aluminium cell coupled with magnetohydrodynamic effects is numerically solved. The thermal model is considered as a two-phases Stefan problem which consists of a non-linear convection–diffusion heat equation with Joule effect as a source. The magnetohydrodynamic fields are governed by Navier–Stokes and by static Maxwell equations. A pseudo-evolutionary scheme (Chernoff) is used to obtain the stationary solution giving the temperature and the frozen layer profile for the simulation of the ledges in the cell. A numerical approximation using a finite element method is formulated to obtain the fluid velocity, electrical potential, magnetic induction and temperature. An iterative algorithm and 3-D numerical results are presented.     

运动热源产生的温度场的有限元分析

本文采用运动元法,考虑了热辐射、热交换和变热传导系数,在较宽的速度范围内研究了运动热源所产生的温度场.给出了在运动坐标和静止坐标系中,不同速度情况下温度场随时间变化的情形,给出了在运动坐标系中不同速度情况下温度场的(定常)等温线分布图.本文还讨论了动态裂纹扩展中裂纹端部过程区的塑性变形所引起的温度场,结果表明,结构钢中过程区的温度一般不会超过1000℃,或1832℉.     

Influence of electromagnetic stirrer position on fluid flow and solidification in continuous casting mold

A computational study of the effect of stirrer position on fluid flow and solidification in a continuous casting billet mold with in-mold electromagnetic stirring has been carried out. The numerical investigation uses a full coupling method in which alternating magnetic field equations are solved simultaneously with the governing equations of fluid flow and heat transfer. An enthalpy-porosity technique is used for the solidification analysis while the magnetohydrodynamics technique is used for studying the fluid flow behavior under the electromagnetic field. The streamline, liquid fraction, and solid shell thickness at the mold wall have been predicted with and without EMS application at different positions along the length of the mold. Recirculation loops are seen to be formed above and below the stirrer position when fluid flow and electromagnetic field equations were solved, without incorporating the solidification model. Application of the solidification model interestingly resulted in the reduction of the size of the recirculation loops formed. The tangential component of velocity of the fluid near the solidification front, stirring intensity and the effective length of stirring below the stirrer decrease as the stirrer position is moved downwards. Significant changes in characteristics of solid shell formation like delay in initiation of solidification at the mold wall and formation of a gap in the re-solidified shell have been observed with change in stirrer position.     

热-压电效应对轴压混合层合圆柱曲板后屈曲的影响

基于Karman-Donnell型非线性壳体方程，给出带压电作动器混合层合圆柱曲板在机械荷载、电荷载和热荷载作用下的后屈曲分析．假定温度场为均匀分布，电场仅有沿板厚方向的分量Ez，且假定材料性能常数与温度和电场的变化无关。将壳体屈曲的边界层理论推广到混合层合圆柱曲板受复合荷载作用的情况.相应的奇异摄动法用于确定圆柱曲板的屈曲荷载和后屈曲平衡路径.分析中同时考虑非线性前屈曲变形和初始几何缺陷的影响．数值算例给出完善和非完善，含整体覆盖或内埋压电作动器正交铺设层合圆柱曲板的后屈曲平衡路径。讨论了温度变化、控制电压、铺层方式、面内边界条件和初始几何缺陷等各种参数变化的影响。     

Modelling ultra-fast nanoparticle melting with the Maxwell–Cattaneo equation

The role of thermal relaxation in nanoparticle melting is studied using a mathematical model based on the Maxwell–Cattaneo equation for heat conduction. The model is formulated in terms of a two-phase Stefan problem. We consider the cases of the temperature profile being continuous or having a jump across the solid–liquid interface. The jump conditions are derived from the sharp-interface limit of a phase-field model that accounts for variations in the thermal properties between the solid and liquid. The Stefan problem is solved using asymptotic and numerical methods. The analysis reveals that the Fourier-based solution can be recovered from the classical limit of zero relaxation time when either boundary condition is used. However, only the jump condition avoids the onset of unphysical “supersonic” melting, where the speed of the melt front exceeds the finite speed of heat propagation. These results conclusively demonstrate that the jump condition, not the continuity condition, is the most suitable for use in models of phase change based on the Maxwell–Cattaneo equation. Numerical investigations show that thermal relaxation can increase the time required to melt a nanoparticle by more than a factor of ten. Thus, thermal relaxation is an important process to include in models of nanoparticle melting and is expected to be relevant in other rapid phase-change processes.     

STUDY OF CRYSTAL GROWTH AND SOLUTE PRECIPITATION THROUGH FRONT TRACKING METHOD

Crystal growth and solute precipitation is a Stefan problem. It is a free boundary problem for a parabolic partial differential equation with a time-dependent phase interface. The velocity of the moving interface between solute and crystal is a local function. The dendritic structure of the crystal interface, which develops dynamically, requires high resolution of the interface geometry. These facts make the Lagrangian front tracking method well suited for the problem. In this paper, we introduce an upgraded version of the front tracking code and its associated algorithms for the numerical study of crystal formation. We compare our results with the smoothed particle hydrodynamics method （SPH） in terms of the crystal fractal dimension with its dependence on the Damkohler number and density ratio.     

A Mixed Analytical/Numerical Method for Velocity and Heat Transfer of Laminar Power-Law Fluids

This paper presents a relatively simple numerical method to investigate the flow and heat transfer of laminar power-law fluids over a semi-infinite plate in the presence of viscous dissipation and anisotropy radiation. On one hand, unlike most classical works, the effects of power-law viscosity on velocity and temperature fields are taken into account when both the dynamic viscosity and the thermal diffusivity vary as a power-law function. On the other hand, boundary layer equations are derived by Taylor expansion, and a mixed analytical/numerical method (a pseudo-similarity method) is proposed to effectively solve the boundary layer equations. This method has been justified by comparing its results with those of the original governing equations obtained by a finite element method. These results agree very well especially when the Reynolds number is large. We also observe that the robustness and accuracy of the algorithm are better when thermal boundary layer is thinner than velocity boundary layer.     

On the role of radiative heat transfer in air gaps in vertical continuous casting

A recent asymptotic thermomechanical model for the formation and evolution of air gaps in vertical continuous casting in round moulds is used as the basis for understanding the importance of radiative heat transfer between the mould wall and the solidified shell. Asymptotic analysis is used to systematically reduce the model to a moving boundary problem with a boundary condition in integro-differential form. In addition, a dimensionless parameter is found which determines whether radiation is significant or not for prescribed process operating conditions. Sample computations are carried out using parameters relevant to the continuous casting of copper, and the results are used to interpret earlier findings in the literature on the interaction between air-gap width and thermal radiation.     

An accurate finite-difference method for ablation-type Stefan problems

A recently derived numerical algorithm for one-dimensional time-dependent Stefan problems is extended for the purpose of solving one-phase ablation-type moving boundary problems; in tandem with the Keller box finite-difference scheme, the so-called boundary immobilization method is used. An important component of the work is the use of variable transformations that must be built into the numerical algorithm in order to preserve second-order accuracy in both time and space. The analysis also determines that the ablation front initially moves as the time raised to the power 3/2; hence, it evolves considerably more slowly than the phase-change front in the classical Stefan problem with isothermal cooling.     

Numerical study of forced and free convective boundary layer flow of a magnetic fluid over a flat plate under the action of a localized magnetic field

The two-dimensional, steady, laminar, forced and free convective boundary layer flow of a magnetic fluid over a semi-infinite vertical plate, under the action of a localized magnetic field, is numerically studied. The magnetic fluid is considered to be water-based with temperature dependent viscosity and thermal conductivity. The study of the boundary layer is separated into two cases. In case I the boundary layer is studied near the leading edge, where it is dominated by the large viscous forces, whereas in case II the boundary layer is studied far from the leading edge of the plate where the effects of buoyancy forces increase. The numerical solution, for these two different cases, is obtained by an efficient numerical technique based on the common finite difference method. Numerical calculations are carried out for the value of Prandl number Pr =  49.832 (water-based magnetic fluid) and for different values of the dimensionless parameters entering into the problem and especially for the magnetic parameter Mn, the viscosity/temperature parameter Θ _r and the thermal/conductivity parameter S*. The analysis of the obtained results show that the flow field is influenced by the application of the magnetic field as well as by the variation of the viscosity and the thermal conductivity of the fluid with temperature. It is hoped that they could be interesting for engineering applications.     

MHD power-law fluid flow and heat transfer over a non-isothermal stretching sheet

This article presents a numerical solution for the magnetohydrodynamic (MHD) non-Newtonian power-law fluid flow over a semi-infinite non-isothermal stretching sheet with internal heat generation/absorption. The flow is caused by linear stretching of a sheet from an impermeable wall. Thermal conductivity is assumed to vary linearly with temperature. The governing partial differential equations of momentum and energy are converted into ordinary differential equations by using a classical similarity transformation along with appropriate boundary conditions. The intricate coupled non-linear boundary value problem has been solved by Keller box method. It is important to note that the momentum and thermal boundary layer thickness decrease with increase in the power-law index in presence/absence of variable thermal conductivity.     

Numerical approach of temperature distribution in a free boundary model for polythermal ice sheets

Summary. A numerical method for solving the thermal subproblem appearing in the modelization of polythermal ice sheets is described. This thermal problem mainly involves three nonlinearities: a reaction term due to the viscous dissipation, a Signorini boundary condition associated to the geothermic flux and an enthalpy term issued from the two phase Stefan formulation of the polythermal regime. The stationary temperature is obtained as the limit of an evolutive problem which is discretized in time with an upwind characteristics scheme and in space with finite elements. The nonlinearities are solved either by Newton-Raphson method or by duality techniques applied to maximal monotone operators. The application of the algorithms provides the dimensionless temperature distribution approximation and allows to identify the cold and temperate ice regions. Received February 1, 1998 / Published online July 28, 1999     

Inverse determination of thermal conductivity using semi-discretization method

In this work a semi-discretization method is presented for the inverse determination of spatially- and temperature-dependent thermal conductivity in a one-dimensional heat conduction domain without internal temperature measurements. The temperature distribution is approximated as a polynomial function of position using boundary data. The derivatives of temperature in the differential heat conduction equation are taken derivative of the approximated temperature function, and the derivative of thermal conductivity is obtained by finite difference technique. The heat conduction equation is then converted into a system of discretized linear equations. The unknown thermal conductivity is estimated by directly solving the linear equations. The numerical procedures do not require prior information of functional form of thermal conductivity. The close agreement between estimated results and exact solutions of the illustrated examples shows the applicability of the proposed method in estimating spatially- and temperature-dependent thermal conductivity in inverse heat conduction problem.     

Unsteady Marangoni convection heat transfer of fractional Maxwell fluid with Cattaneo heat flux

Fractional shear stress and Cattaneo heat flux models are introduced in characterizing unsteady Marangoni convection heat transfer of viscoelastic Maxwell fluid over a flat surface. Governing equations and boundary condition are formulated firstly via the balance between the surface tension and shear stress. Numerical solutions are obtained by new developed numerical technique and some novel phenomena are found. Results shown that the fractional derivative parameters, Marangoni number and power law exponent have significant influence on characteristics velocity and temperature fields. As fractional derivative parameters increase, the temperature profiles rise remarkably and the viscoelastic effects of the fluid enhance with delayed response to surface tension, however the temperature profiles decline significantly with a thinner thickness of thermal boundary layer with the increase of Marangoni number. The average skin friction coefficient increases with the augment of Marangoni number, while the average Nusselt number decreases for larger values of power law exponent.     

Numerical study of forced and free convective boundary layer flow of a magnetic fluid over a flat plate under the action of a localized magnetic field

The two-dimensional, steady, laminar, forced and free convective boundary layer flow of a magnetic fluid over a semi-infinite vertical plate, under the action of a localized magnetic field, is numerically studied. The magnetic fluid is considered to be water-based with temperature dependent viscosity and thermal conductivity. The study of the boundary layer is separated into two cases. In case I the boundary layer is studied near the leading edge, where it is dominated by the large viscous forces, whereas in case II the boundary layer is studied far from the leading edge of the plate where the effects of buoyancy forces increase. The numerical solution, for these two different cases, is obtained by an efficient numerical technique based on the common finite difference method. Numerical calculations are carried out for the value of Prandl number Pr =  49.832 (water-based magnetic fluid) and for different values of the dimensionless parameters entering into the problem and especially for the magnetic parameter Mn, the viscosity/temperature parameter Θ _r and the thermal/conductivity parameter S*. The analysis of the obtained results show that the flow field is influenced by the application of the magnetic field as well as by the variation of the viscosity and the thermal conductivity of the fluid with temperature. It is hoped that they could be interesting for engineering applications.