Improved lumped models for transient heat conduction in a slab with temperature-dependent thermal conductivity

This work reports improved lumped-parameter models for transient heat conduction in a slab with temperature-dependent thermal conductivity. The improved lumped models are obtained through two point Hermite approximations for integrals. For linearly temperature-dependent thermal conductivity, it is shown by comparison with numerical solution of the original distributed parameter model that the higher order lumped model (_H1,₁/_H0,₀$H_{1\text{,}1}/H_{0\text{,}0}$ approximation) yields significant improvement of average temperature prediction over the classical lumped model. A unified Biot number limit depending on a single dimensionless parameter β$\beta$ is given both for cooling and heating processes.     

Spatial-resolution,lumped-capacitance thermal model for cylindrical Li-ion batteries under high Biot number conditions

The time-efficient yet accurate thermal modeling of the battery cells for electric and hybrid electric vehicles is essential improving the performance, safety, and lifetime of the battery system. This paper presents a spatial-resolution, lumped-capacitance (LC) thermal model for cylindrical battery cells under high Biot number (Bi ? 1) conditions where the classical LC thermal model is generally inapplicable because of a significant temperature variation in the cell volume. The spatial-resolution LC model was formulated using zero- and first-order Hermite integral approximations. For model validation, a one-dimensional, transient analytical (exact) solution using Green functions was obtained for a cylindrical Li-ion battery cell with uniform volumetric battery heat generation of Joule and entropic heating under convective cooling boundary conditions. It was found from the comparison of the results that the spatial-resolution LC thermal model can accurately and quickly predicts the cell temperatures (core, skin and area-averaged) under various dynamic battery duty cycles even for high Biot numbers due to highly convective conditions such as liquid cooling.     

A Stefan problem for a non-classical heat equation with a convective condition

We prove the existence and uniqueness, local in time, of the solution of a one-phase Stefan problem for a non-classical heat equation for a semi-infinite material with a convective boundary condition at the fixed face x = 0. Here the heat source depends on the temperature at the fixed face x = 0 that provides a heating or cooling effect depending on the properties of the source term. We use the Friedman-Rubinstein integral representation method and the Banach contraction theorem in order to solve an equivalent system of two Volterra integral equations. We also obtain a comparison result of the solution (the temperature and the free boundary) with respect to the one corresponding with null source term.     

Analysis of flow and thermal field in nanofluid using a single phase thermal dispersion model

Flow and thermal field in nanofluid is analyzed using single phase thermal dispersion model proposed by Xuan and Roetzel [Y. Xuan, W. Roetzel, Conceptions for heat transfer correlation of nanofluids, Int. J. Heat Mass Transfer 43 (2000) 3701–3707]. The non-dimensional form of the transport equations involving the thermal dispersion effect is solved numerically using semi-explicit finite volume solver in a collocated grid. Heat transfer augmentation for copper–water nanofluid is estimated in a thermally driven two-dimensional cavity. The thermo-physical properties of nanofluid are calculated involving contributions due to the base fluid and nanoparticles. The flow and heat transfer process in the cavity is analyzed using different thermo-physical models for the nanofluid available in literature. The influence of controlling parameters on convective recirculation and heat transfer augmentation induced in buoyancy driven cavity is estimated in detail. The controlling parameters considered for this study are Grashof number (10³ < Gr < 10⁵), solid volume fraction (0 < ? < 0.2) and empirical shape factor (0.5 < n < 6). Simulations carried out with various thermo-physical models of the nanofluid show significant influence on thermal boundary layer thickness when the model incorporates the contribution of nanoparticles in the density as well as viscosity of nanofluid. Simulations incorporating the thermal dispersion model show increment in local thermal conductivity at locations with maximum velocity. The suspended particles increase the surface area and the heat transfer capacity of the fluid. As solid volume fraction increases, the effect is more pronounced. The average Nusselt number from the hot wall increases with the solid volume fraction. The boundary surface of nanoparticles and their chaotic movement greatly enhances the fluid heat conduction contribution. Considerable improvement in thermal conductivity is observed as a result of increase in the shape factor.     

Numerical study of melting in an enclosure with discrete protruding heat sources

In order to explore the capability of a solid–liquid phase change material (PCM) for cooling electronic or heat storage applications, melting of a PCM in a vertical rectangular enclosure was studied. Three protruding generating heat sources are attached on one of the vertical walls of the enclosure, and generating heat at a constant and uniform volumetric rate. The horizontal walls are adiabatic. The power generated in heat sources is dissipated in PCM (n-eicosane with the melting temperature, T_m = 36 °C) that filled the rectangular enclosure. The advantage of using PCM is that it is able to absorb high amount of heat generated by heat sources due to its relatively high energy density. To investigate the thermal behaviour and thermal performance of the proposed system, a mathematical model based on the mass, momentum and energy conservation equations was developed. The governing equations are next discretised using a control volume approach in a staggered mesh and a pressure correction equation method is employed for the pressure–velocity coupling. The PCM energy equation is solved using the enthalpy method. The solid regions (wall and heat sources) are treated as fluid regions with infinite viscosity and the thermal coupling between solid and fluid regions is taken into account using the harmonic mean of the thermal conductivity method. The dimensionless independent parameters that govern the thermal behaviour of the system were next identified. After validating the proposed mathematical model against experimental data, a numerical investigation was next conducted in order to examine the thermal behaviour of the system by analyzing the flow structure and the heat transfer during the melting process, for a given values of governing parameters.     

Computation of melting with natural convection inside a rectangular enclosure heated by discrete protruding heat sources

This paper presents the results of a numerical investigation of the heat transfer by natural convection during the melting of a phase change material (PCM, n-eicosane with melting point of 36 °C) contained in a rectangular enclosure. This latest is heated by three discrete protruding heat sources (simulating electronic components) placed on one of its vertical walls. The power generated by heat sources is dissipated in PCM. The advantage of using this cooling scheme is that the PCMs are able to absorb high amount of heat generated by the heat sources, without acting the fan during the charging process (melting of the PCM). The thermal behavior and thermal performance of the proposed PCM based-heat sink are numerically investigated by developing a mathematical model based on the mass, momentum and energy conservation equations. The obtained numerical results show the impact of various key parameters on the cooling capacity of the PCM-based heat sink. Correlations encompassing a wide range of parameters were developed in terms of the dimensionless secured operating time (time required by one of the electronic components before reaching its critical temperature, T_cr ∼ 75 °C) and the corresponding liquid fraction, using the asymptotic computational fluid dynamics (ACFD) technique.     

The study of structure optimization of blast furnace cast steel cooling stave based on heat transfer analysis

The three-dimensional mathematical model of temperature and thermal stress field of cast steel cooling stave in a blast furnace has been modeled. Kinds of the parameters optimization of cast steel cooling stave in a blast furnace are proposed based on the heat transfer analysis. The results indicate that the values of the parameters optimization for a cast steel cooling stave are 200 mm for cooling channels interdistance, 25 mm for inner radius of the water channel, 180 mm for thickness of the cooling stave body, 70 mm for thickness of inlaid brick and 1.5 m/s for speed of cooling water. Reducing the water temperature would be uneconomical. The water temperature can be chosen according to the local conditions. The best choice for lining material is silicon nitrogen bond silicon carbide brick or silicon carbide brick.     

Numerical prediction of buoyancy driven micropolar fluid flow within uniformly heated eccentric annulus

In this paper, the problem of buoyancy driven micropolar fluid flow within an annulus formed between two circular concentric/eccentric tubes has been numerically investigated using Fourier spectral method. The annulus inner wall is uniformly heated and maintained at constant heat flux while the outer wall is cooled and kept at constant temperature. The full governing equations of linear momentum, angular momentum and energy have been solved to give the details of flow and thermal fields. The heat convection process in the annulus is mainly controlled by modified Rayleigh number Ra, Prandtl number Pr, radius ratio Rr, eccentricity, e and material parameters of Micropolar fluid. The material parameters are dimensionless spin gradient viscosity λ, dimensionless micro-inertia density B and dimensionless vortex viscosity D. The study considered a range of modified Ra up to 10⁵ and is carried out at three values of Pr, namely Pr = 0.1, 1.0 and 7.0, and at three values of parameter D, namely, D = 2, 4, 8 while the eccentricity is varied between −0.65 and +0.65. The radius ratio is fixed at 2.6 while the material parameters B and λ are assigned the value of 1. The effect of the controlling parameters on flow and thermal fields has been investigated with emphasis on the effect of these parameters on local and mean inner wall temperatures. The study has shown that for certain controlling parameters the steady mean temperature of inner wall of the annulus is maximum at a certain eccentricity. The study has also shown that as the parameter D increases the steady mean inner wall temperature increases. Moreover, the study has shown that as the Pr increases the mean inner wall temperature decreases.     

Inverse heat conduction problem in a thin circular plate and its thermal deflection

An inverse problem of transient heat conduction in a thin finite circular plate with the given temperature distribution on the interior surface of a thin circular plate being a function of both time and position has been solved with the help of integral transform technique and also determine the thermal deflection on the outer curved surface of a thin circular plate defined as 0 ? r ? a, 0 ? z ? h. The results, obtained in the series form in terms of Bessel’s functions, are illustrated numerically.     

Non-local effects and size-dependent properties in Stefan problems with Newton cooling

We model the growth of a one-dimensional solid by considering a modified Fourier law with a size-dependent effective thermal conductivity and a Newton cooling condition at the interface between the solid and the cold environment. In the limit of a large Biot number, this condition becomes the commonly used fixed-temperature condition. It is shown that in practice the size of this non-dimensional number is very small. We study the effect of a small Biot number on the solidification process with numerical and asymptotic solution methods. The study indicates that non-local effects become less important as the Biot number decreases.     

New types of 3-D systems of quadratic differential equations with chaotic dynamics based on Ricker discrete population model

The wide class of 3-D autonomous systems of quadratic differential equations, in each of which either there is a couple of coexisting limit cycles or there is a couple of coexisting chaotic attractors, is found. In the second case the couple consists of either Lorentz-type attractor and another attractor of a new type or two Lorentz-type attractors. It is shown that the chaotic behavior of any system of the indicated class can be described by the Ricker discrete population model: z_i+1 = z_i exp(r − z_i), r > 0, z_i > 0, i = 0, 1, … . The values of parameters, at which in the 3-D system appears either the couple of limit cycles or the couple of chaotic attractors, or only one limit cycle, or only one sphere-shaped chaotic attractor, are indicated. Examples are given.     

Axi symmetric 2D simulation and numerical heat transfer characteristics for the calibration furnace in a rectangular enclosure

This paper presents axi symmetric 2D numerical investigation of the spherical thermocouple calibration furnace in a rectangular enclosure. The focus is on the flow structure inside the Saturn (a hollow spherical cavity), external flow behavior due to annulus block heating and the surface temperature uniformity. Mesh sensitivity analysis is adopted to extract the mesh with minimum number of nodes but with fast convergent finite element solution. The basic strategy here is that temperature perturbation error at a single point instead of a single element contributed to the total perturbation error qualitatively remains the same. Agreement between numerical simulation results and the experiment results is good with a maximum temperature deviation 10 °C for the cavity temperature 400 °C. Finally, standard numerical temperature uncertainty due to variation in thermal conductivity is computed through the sensitivity coefficient using uncertainty analysis.     

Krylov type subspace methods for matrix polynomials

We consider solving eigenvalue problems or model reduction problems for a quadratic matrix polynomial Iλ² − Aλ − B with large and sparse A and B. We propose new Arnoldi and Lanczos type processes which operate on the same space as A and B live and construct projections of A and B to produce a quadratic matrix polynomial with the coefficient matrices of much smaller size, which is used to approximate the original problem. We shall apply the new processes to solve eigenvalue problems and model reductions of a second order linear input-output system and discuss convergence properties. Our new processes are also extendable to cover a general matrix polynomial of any degree.     

Thermal dispersion driven by the spontaneous imbibition process

In this work, we have theoretically analyzed the thermal dispersion process under the influence of the spontaneous imbibition of a liquid trapped in a capillary element, considering the presence of a uniform temperature gradient. The capillary element is represented by a porous medium which is initially found at temperature T₀ and pressure P₀. Suddenly, the lower part of the porous medium touches a liquid reservoir at temperature T_l and pressure P₀. This contact between both phases, in turn causes spontaneously the imbibition process. Using a one-dimensional formulation of the average conservation laws, we derive the corresponding nondimensional momentum and energy equations. The numerical solutions permit us to evaluate the position and velocity of the imbibition front as well as the temperature profiles and Nusselt numbers. The above results are shown by taking into account the influence of three dimensionless parameters: the ratio of the characteristic thermal time to the characteristic imbibition time, β, the ratio of the hydrostatic head of the imbibed liquid to the characteristic pressure difference for the imbibition front, α, and the ratio of the dispersive thermal diffusivity to the effective thermal diffusivity of the medium, Ω. The predictions show that temperature profiles and the heat transfer process are strongly dependent on thermal dispersion effects, indicating a clear deviation in comparison with the case of Ω = 0 that represents the absence of the thermal dispersion.     

Heating of an electron gas by large temperature gradients and scale thermoelectric effects in semiconductors

We study the temperature fields and thermoelectric effects in semiconductors of monopolar conductivity under conditions of heating of an electron gas by large temperature gradients. We study the influence of surface cooling of the current carriers on the scale effects.Translated fromMatematicheskie Metody i Fiziko-Mekhanicheskie Polya, Issue 31, 1990, pp. 16–21.     

Numerical analysis of the thermal behaviour of a shell-and-tube heat storage unit using phase change materials

This work presents a numerical study of a latent heat storage unit (LHSU) consisting of a shell-and-tube. The shell space is filled with two phase change materials (PCMs), P116 and n-octadecane, with different melting temperatures (50 °C and 27.7 °C, respectively). A heat transfer fluid (HTF: water) flows by forced convection through the inner tube, and transfers the heat to PCMs. In order to compare the thermal performances of the latent heat storage unit using two phase change materials (LHSU2) and a single PCM (LHSU1), a mathematical model based on the conservation energy equations was developed and validated with experimental data. Several numerical investigations were conducted in order to examine the impact of the key parameters: the HTF inlet temperature (ranges from 50 to 60 °C), the mass flow rate of the HTF and the proportion mass of PCMs, on the thermal performances of the latent heat storage units using two PCMs and a single PCM, during charging process (melting). This parametric study provides guidelines for system thermal performance and design optimization.     

Two-worker blocking congestion model with walk speed m in a no-passing circular passage system

This paper introduces a blocking model and closed-form expression of two workers traveling with walk speed m (m = integer) in a no-passing circular-passage system of n stations and assuming n = m + 2, 2m + 2, …. We develop a Discrete-Timed Markov Chain (DTMC) model to capture the workers’ changes of walk, pick, and blocked states, and quantify the throughput loss from blocking congestion by deriving a steady state probability in a closed-form expression. We validate the model with a simulation study. Additional simulation comparisons show that the proposed throughput model gives a good approximation of a general-sized system of n stations (i.e., n > 2), a practical walk speed system of real number m (i.e., m ? 1), and a bucket brigade order picking application.     

Symmetric duality with (p, r) − ρ − (η, θ)-invexity

In this paper we introduce a new type of generalized invex function, called (p, r) − ρ − (η, θ)-invex function and study symmetric duality results under these assumptions. In our study the nonnegative orthants for the constraints are replaced by closed convex cones and their polars. We establish weak and strong duality theorems under (p, r) − ρ − (η, θ)-invexity assumptions for the symmetric dual problems. We also give many examples to justify our results.     

Analyses on 3-D gas flow and heat transfer in ladle furnace lid

In order to verify the reasonableness of off-gas pressure and wall temperature, a mathematical model for gas flow and heat transfer in ladle furnace (LF) lid is developed based on 3-D Navier–Stokes equations and k–ε two equation turbulent models as well as energy conservation equation. The gas velocity vector distribution of skirt clearance between the top edge of ladle and furnace lid and electrode gaps between three graphite electrodes and furnace lid, the gas flow line distribution, pressure and temperature distribution on the furnace lid wall are simulated. Simulation results show that appropriate off-gas pressures are 200 Pa, 200 Pa and 150 Pa when electric arc emerges from molten steel surface and alloy hole is unsealed, electric arc emerges from molten steel surface and alloy hole is closure, electric arc immerges into molten steel surface and alloy hole is closure, respectively. The maximum temperature presents in the middle of LF lid in all of heating conditions, and the temperature value are 563, 603 and 343 K. Finally, the relations between gas volume and off-gas pressure are analyzed in different width of skirt clearance, and some relevant mathematical expressions are obtained. By comparing both simulation results and practical data, the advice on reducing off-gas pressure is proposed, and the maximum temperatures of furnace lid wall have good agreement with actual data.     

A new family of unsteady boundary layers over a stretching surface

In this paper, a new family of unsteady boundary layers over a stretching flat surface was proposed and studied. This new class of unsteady boundary layers involves the flows over a constant speed stretching surface from a slot, and the slot is moving at a certain speed. Depending on the slot moving parameter, the flow can be treated as a stretching sheet problem or a shrinking sheet problem. Both the momentum and thermal boundary layers were studied. Under special conditions, the solutions reduce to the unsteady Rayleigh problem and the steady Sakiadis stretching sheet problem. Solutions only exist for a certain range of the slot moving parameter, α. Two solutions are found for −53.55° < α < −45°. There are also two solution branches for the thermal boundary layers at any given Prandtl number in this range. Compared with the upper solution branch, the lower solution branch leads to simultaneous reduction in wall drag and heat transfer rate. The results also show that the motion of the slot greatly affects the wall drag and heat transfer characteristics near the wall and the temperature and velocity distributions in the fluids.