Unsteady thermal entrance heat transfer in laminar pipe flows with step change in ambient temperature

A fully implicit upwind finite difference numerical scheme has been proposed to investigate the characteristics of thermal entrance heat transfer in laminar pipe flows subject to a step change in ambient temperature. In order to demonstrate the results more clearly, a modified Nusselt number is introduced. The unsteady axial variations of modified Nusselt number, bulk fluid temperature, and wall temperature and the transient temperature profiles at certain axial locations are presented graphically for various outside heat transfer coefficients. The effects of the outside heat transfer coefficient on the heat transport processes in the flow are examined in detail. The results can be comprehensively explained by the interaction between the upstream convective heat transfer and the diffusion heat transfer in the radial direction. Steady state is reached when the axial convection balances the radial diffusion.     

Direct numerical simulation of turbulent heat transfer in a wall-normal rotating channel flow

Investigations into the characteristics of turbulent heat transfer and coherent flow structures in a plane-channel subjected to wall-normal system rotation are conducted using direct numerical simulation (DNS). In order to investigate the influence of system rotation on the temperature field, a wide range of rotation numbers are tested, with the flow pattern transitioning from being fully turbulent to being quasilaminar, and eventually, fully laminar. In response to the Coriolis force, secondary flows appear as large vortical structures, which interact intensely with the wall shear layers and have a significant impact on the distribution of turbulence kinetic energy (TKE), turbulence scalar energy (TSE), temperature statistics, and turbulent heat fluxes. The characteristic length scales of turbulence structures responsible for the transport of TSE are the largest at the quasilaminar state, which demands a very large computational domain in order to capture the two-dimensional spectra of temperature fluctuations. The effects of the Coriolis force on the turbulent transport processes of the temperature variance and turbulent heat fluxes are thoroughly examined in terms of their respective budget balances.     

A physically based flux limiter for QUICK calculations of advective scalar transport

Transient, advective transport of a contaminant into a clean domain will exhibit a moving sharp front that separates contaminated and clean regions. Due to ‘numerical diffusion’—the combined effects of ‘cross‐wind diffusion’ and ‘artificial dispersion’—a numerical solution based on a first‐order (upwind) treatment will smear out the sharp front. The use of higher‐order schemes, e.g. QUICK (quadratic upwinding) reduces the smearing but can introduce non‐physical oscillations in the solution. A common approach to reduce numerical diffusion without oscillations is to use a scheme that blends low‐order and high‐order approximations of the advective transport. Typically, the blending is based on a parameter that measures the local monotonicity in the predicted scalar field. In this paper, an alternative approach is proposed for use in scalar transport problems where physical bounds C_Low?C?C_High on the scalar are known a priori. For this class of problems, the proposed scheme switches from a QUICK approximation to an upwind approximation whenever the predicted upwind nodal value falls outside of the physical range [C_Low, C_High]. On two‐dimensional steady‐state and one‐dimensional transient test problems predictions obtained with the proposed scheme are essentially indistinguishable from those obtained with monotonic flux‐limiter schemes. An analysis of the modified equation explains the observed performance of first‐ and second‐order time‐stepping schemes in predicting the advective transport of a step. In application to the transient two‐dimensional problem of contaminate transport into a streambed, predictions obtained with the proposed flux‐limiter scheme agree with those obtained with a scheme from the literature. Copyright © 2007 John Wiley & Sons, Ltd.     

Variable viscosity in peristalsis of Sisko fluid

This paper aims to examine variable viscosity effects on peristalsis of Sisko fluids in a curved channel with compliant characteristics. Viscous dissipation in a heat transfer is studied. The resulting problems are solved using perturbation and numerical schemes to show qualitatively similar responses for velocity and temperature. A streamline phenomenon is also considered.     

相变传热问题的灵敏度分析与优化设计方法

研究了相变传热问题的优化设计及其灵敏度分析方法. 在有限元-时间差分和等效热容 法求解相变温度场的基础上，提出了相变温度场对设计变量一阶灵敏度的计算方法，给出直 接法和伴随法两种计算格式并分析了它们的特点，建立了相变温度场优化的模型和算法，在有限元分析与优化设计软件JIFEX中实现了该方法. 数值算例表明了灵敏度计算的精度和优 化方法的有效性.     

An efficient 3D particle transport model for use in stratified flow

Random‐walk models are a versatile tool for modelling dispersion of both passive and active tracers in turbulent flow. The physical and mathematical foundations of stochastic Lagrangian models of turbulent diffusion have become more and more solid over the years. An important aspect of these types of models that has not received much attention is the behaviour of the particles near boundaries. Often, a simple stochastic, numerical scheme is used. Because turbulent mixing in the vertical direction is much more complicated than in the two horizontal directions, it is in the vertical direction that a simple numerical scheme, such as the Euler scheme, may cause problems. In this paper our main goal is the development of an efficient 3D particle transport model that can be used in stratified flow. For this type of situation the vertical direction is of special interest. First, a closer look is taken at some considerations that should be regarded when choosing a numerical scheme. Specifically schemes are investigated that can be used in the vertical direction, where the diffusion coefficient is varying in that direction. Experiments are performed regarding the accuracy of different numerical schemes in various situations. The behaviour of the particles near an impermeable layer interface is investigated. The stochastic Heun and Runge–Kutta schemes turn out to be very attractive for this type of model. For the simulation of the transport of various physical quantities, such as salinity, heat, silt, oxygen, or bacteria, different types of models are available. In this case we will take a closer look at the modelling of the transport of pollutants from point sources (either instantaneous or continuous transport). For this purpose a 3D particle transport model has been developed that is especially suited for stratified situations such as can be found in estuaries. The main idea is to use a simple numerical scheme for the horizontal directions and a higher‐order method for the vertical direction. The results play an important role in making specific choices for this type of particle transport model. Copyright © 2005 John Wiley & Sons, Ltd.     

Temperature and related statistics in turbulent jet flames

This paper describes the use of digitally-compensated thermocouples to characterise the time-resolved behaviour of the temperature field in turbulent non-premixed jet flames for Re?2×10⁴, including the transport of heat and the related temperature dissipation. The experimental procedures used are analysed, including the thermocouple time constant determination and the numerical derivative compensation method, and assessments of accuracy are provided for the cross-correlation techniques used in order to estimate high-order temperature statistics in the flames studied. The results quantify the influence of the spatial resolution of the measuring systems on the accuracy of correlation values, regarding the characteristics scales of the flames studied, and demonstrate the applicability of fine-wire thermocouples to characterise the turbulent transport and the dissipation of temperature in non-premixed jet flames.     

Optimization of convective fin systems: a holistic approach

 A numerical analysis of natural convection heat transfer and entropy generation from an array of vertical fins, standing on a horizontal duct, with turbulent fluid flow inside, has been carried out. The analysis takes into account the variation of base temperature along the duct, traditionally ignored by most studies on such problems. One-dimensional fin equation is solved using a second order finite difference scheme for each of the fins in the system and this, in conjunction with the use of turbulent flow correlations for duct, is used to obtain the temperature distribution along the duct. The influence of the geometric and thermal parameters, which are normally employed in the design of a thermal system, has been studied. Correlations are developed for (i) the total heat transfer rate per unit mass of the fin system (ii) total entropy generation rate and (iii) fin height, as a function of the geometric parameters of the fin system. Optimal dimensions of the fin system for (i) maximum heat transfer rate per unit mass and (ii) minimum total entropy generation rate are obtained using Genetic Algorithm. As expected, these optima do not match. An approach to a `holistic' design that takes into account both these criteria has also been presented. Received on 22 February 2001 / Published online: 29 November 2001     

Convective drying analysis of three-dimensional porous solid by mass lumping finite element technique

A numerical analysis of convective drying of a 3D porous solid of brick material is carried out using the finite element method and mass lumping technique. The energy equation and moisture transport equations for the porous solid are derived based on continuum approach following Whitaker’s theory of drying. The governing equations are solved using the Galerkin’s weighted residual method, which convert the governing equations into discretized form of matrix equations. The resulting capacitance matrices are made diagonal matrices by following the classical row-sum mass lumping technique. Hence with the use of the Eulerian time marching scheme, the final equations are reduced to simple algebraic equations, which can be solved directly without using an equation solver. The proposed numerical scheme is initially validated with experimental results for 1D drying problem and then tested by application to convective drying of 3D porous solid of brick material for four different aspect ratios obtained by varying the cross section of the solid. The mass lumping technique could correctly predict the wet bulb temperature of the solid under evaporative drying conditions. A parametric study carried out for three different values of convective heat transfer coefficients, 15, 30 and 45 W/m² K shows an increased drying rate with increase in area of cross section and convective heat transfer coefficient. The proposed numerical scheme could correctly predict the drying behavior shown in the form of temperature and moisture evolutions.     

Natural convection and radiation heat transfer of an externally-finned tube vertically placed in a chamber

A three-dimensional numerical study was made to investigate effects of fin angle, fin surface emissivity, and tube wall temperature on heat transfer enhancement for a longitudinal externally-finned tube placed vertically in a small chamber. The numerical model was first validated through comparison with experimental measurements and the appropriateness of general boundary conditions was examined. The numerical results show that the mean Nusselt number increases with Rayleigh number for all the fin angles investigated. The maximum heat transfer rate per mass occurs when the fin angle is about 60° for fin surface emissivity between 0.7 and 0.8 and 55° when the surface emissivity increases to 0.9. With increasing tube wall temperature, both the natural convection and radiation heat transfer are enhanced, but the fraction of radiation heat transfer decreases in the temperature range studied. Radiation fraction increases with increasing fin surface emissivity. Both convection and radiation heat transfer modes are important.     

Numerical investigation of premixed combustion in a porous burner with integrated heat exchanger

In this paper, we perform a numerical analysis of a two-dimensional axisymmetric problem arising in premixed combustion in a porous burner with integrated heat exchanger. The physical domain consists of two zones, porous and heat exchanger zones. Two dimensional Navier–Stokes equations, gas and solid energy equations, and chemical species transport equations are solved and heat release is described by a multistep kinetics mechanism. The solid matrix is modeled as a gray medium, and the finite volume method is used to solve the radiative transfer equation to calculate the local radiation source/sink in the solid phase energy equation. Special attention is given to model heat transfer between the hot gas and the heat exchanger tube. Thus, the corresponding terms are added to the energy equations of the flow and the solid matrix. Gas and solid temperature profiles and species mole fractions on the burner centerline, predicted 2D temperature fields, species concentrations and streamlines are presented. Calculated results for temperature profiles are compared to experimental data. It is shown that there is good agreement between the numerical solutions and the experimental data and it is concluded that the developed numerical program is an excellent tool to investigate combustion in porous burner.     

An internal penalty discontinuous Galerkin method for simulating conjugate heat transfer in a closed cavity

Using the discontinuous Galerkin (DG) method for conjugate heat transfer problems can provide improved accuracy close to the fluid‐solid interface, localizing the data exchange process, which may further enhance the convergence and stability of the entire computation. This paper presents a framework for the simulation of conjugate heat transfer problems using DG methods on unstructured grids. Based on an existing DG solver for the incompressible Navier‐Stokes equation, the fluid advection‐diffusion equation, Boussinesq term, and solid heat equation are introduced using an explicit DG formulation. A Dirichlet‐Neumann partitioning strategy has been implemented to achieve the data exchange process via the numerical flux of interface quadrature points in the fluid‐solid interface. Formal h and p convergence studies employing the method of manufactured solutions demonstrate that the expected order of accuracy is achieved. The algorithm is then further validated against 3 existing benchmark cases, including a thermally driven cavity, conjugate thermally driven cavity, and a thermally driven cavity with conducting solid, at Rayleigh numbers from 1000 to 100 000. The computational effort is documented in detail demonstrating clearly that, for all cases, the highest‐order accurate algorithm has several magnitudes lower error than first‐ or second‐order schemes for a given computational effort.     

Higher-order accurate implicit time integration schemes for transport problems

The present paper is concerned with the numerical solution of transient transport problems by means of spatial and temporal discretization methods. The generalized initial boundary value problem of various nonlinear transport phenomena like heat transfer or mass transport is discretized in space by p-finite elements. After finite element discretization, the resulting first-order semidiscrete balance has to be solved with respect to time. Next to the classical generalized-α integration method predicated on the Newmark approach and the evaluation at a generalized midpoint also implicit Runge–Kutta time integration schemes, are presented. Both families of finite difference-based integration schemes are derived for general first-order problems. In contrast to the above-mentioned algorithms, temporal discontinuous and continuous Galerkin methods evaluate the balance equation not at a selected time instant within the timestep, but in an integral sense over the whole time step interval. Therefore, the underlying semidiscrete balance and the continuity of the primary variables are weakly formulated within time steps and between time steps, respectively. Continuous Galerkin methods are obtained by the strong enforcement of the continuity condition as special cases. The introduction of a natural time coordinate allows for the application of standard higher-order temporal shape functions of the p-Lagrange type and the well-known Gau?–Legendre quadrature of associated time integrals. It is shown that arbitrary order accurate integration schemes can be developed within the framework of the proposed temporal p-Galerkin methods. Selected benchmark analyses of calcium diffusion demonstrate the properties of all three methods with respect to non-smooth initial or boundary conditions. Furthermore, the robustness of the present time integration schemes is also demonstrated for the highly nonlinear reaction–diffusion problem of calcium leaching, including the pronounced changes of the reaction term and non-smooth changes of Dirichlet boundary conditions of calcium dissolution.     

Numerical and experimental investigation of convective drying in unsaturated porous media with bound water

The heat and mass transfer in an unsaturated wet cylindrical porous bed packed with quartz particles was investigated theoretically for relatively low convective drying rates. Local thermodynamic equilibrium was assumed in the mathematical model describing the multi-phase flow in the unsaturated porous media using the energy and mass conservation equations to describe the heat and mass transfer during the drying. The drying model included convection and capillary transport of the free water, diffusion of bound water, and convection and diffusion of the gas. The numerical results indicated that the drying process could be divided into three periods, the temperature rise period, the constant drying rate period and the decreasing drying rate period. The numerical results agreed well with the experimental data verifying that the mathematical model can evaluate the drying performance of porous media for low drying rates. The effects of drying conditions such as the ambient temperature, the relative humidity, and the velocity of the drying air, on the drying process were evaluated by numerical solution.     

A high-order discontinuous Galerkin method for simulating incompressible fluid-thermal-structural interaction problems

This paper presents a framework for the application of the discontinuous Galerkin(DG) finite element method to the multi-physics simulation of the solid thermal deformation interacting with incompressible flow problems in two-dimensions. Recent applications of the DG method are primarily for thermoelastic problems in a solid domain or fluid-structure interaction problems without heat transfer. Based on a recently published conjugate heat transfer solver, the incompressible Navier-Stokes equation, the fluid advection-diffusion equation, the Boussinesq term, the solid heat equation and the solid linear elastic equation are solved using an explicit DG formulation. A Dirichlet-Neumann partitioning strategy has been implemented to achieve the data exchange process via the numerical flux computed at interface quadrature points in the fluid-solid interface. Formal hp convergence studies employing the method of manufactured solutions demonstrate that the expected order of accuracy is achieved for each solver. The algorithm is then further validated against several existing benchmark cases including the in-plane loaded square, the Timoshenko Beam, the laminated beam subject to thermal-loads and the lid-driven cavity with a flexible bottom wall. The computational effort demonstrates that for all cases the highest order accurate algorithm has several magnitudes lower error than the second-order schemes for a given computational effort. It is a strong justification for the development of such high order discretisations. The solver can be employed to predict thermal deformation of a structure due to convective and conductive heat transfer at low Mach, such as chip deformation on a printed circuit board, wave-guide structure optimization, thermoelectric cooler simulation, and optics mounting method verification.     

Approximate analysis for transient heat conduction with radiation in a slab

In an attempt to minimize the numerical computations associated with the solution of transient heat conduction with radiation in a slab, a perturbation type of analysis is being applied to the temperature field and radiation heat flux simultaneously. The resulting partial differential equations for the perturbation functions for the temperature are solved in explicit forms by use of the energy integral methods, while the radiation heat flux is determined by an appropriate scheme of approximating the temperature distribution in the slab. Included in the analysis are the effects of the parameters: the optical thickness, the ratio of conduction transport to radiation and the wall emissivity. It is found that, in a wide range of these governing parameters, the results compare very favorably with those obtained by the numerical solution of the formulated integro-differential equation. With the present analysis, the temperature, conduction and radiation heat fluxes can be predicted without resorting to lengthy numerical analysis.     

Characteristics of air and water velocity fields during natural convection

Air and water velocity fields have been simulated during natural convection, using a two-dimensional volume of fluid (VOF) model. The results have shown that during unstable thermal stratification, the root-mean-square (RMS) airside velocities are an order of magnitude higher than the RMS waterside velocities, whereas, during the stable thermal stratification, the velocity magnitudes are comparable for air and water sides. Furthermore, the magnitude of the air velocity changed more rapidly with the change in the bulk air–water temperature difference than the water velocity, indicating that the air velocities are more sensitive to the bulk air and water temperature difference than the water velocities. A physical model of the heat and mass transfer across the air–water interface is defined. According to this model, the vortices on the air and water sides play an important role in enhancing the heat and mass transfer. Due to the significance of the flow velocities in the transport process, it has been proposed that the correlations for the heat and mass transfer during natural convection should be improved by incorporating the flow velocity as a parameter.     

Numerical analysis of mixed convection in three‐dimensional rectangular channel with flush‐mounted heat sources based on field synergy principle

In this paper the fluid flow and heat transfer characteristics of mixed convection in three‐dimensional rectangular channel with four heat sources are investigated numerically. The SIMPLEC algorithm is applied to deal with the coupling between pressure and velocity, and a new high‐order stability‐guaranteed second‐order difference (SGSD) scheme is adopted to discretize the convection term. The influence of four parameters is studied: Richardson number, heat source distribution, channel height and inclination angle. The numerical results are analysed from the viewpoint of the field synergy principle, which says that the enhanced convective heat transfer is related not only to the velocity field and temperature field, but also to the synergy between them. It is found that the effects of the four parameters on the thermal performance can all be explained with the field synergy principle. To obtain better electronic cooling, the synergy between the velocity and temperature gradient should be increased when other conditions are unchanged. Copyright © 2006 John Wiley & Sons, Ltd.     

3D numerical simulation on fluid flow and heat transfer characteristics in multistage heat exchanger with slit fins

In this paper, a numerical investigation is performed for three-stage heat exchangers with plain plate fins and slit fins respectively, with a three-dimensional laminar conjugated model. The tubes are arranged in a staggered way, and heat conduction in fins is considered. In order to save the computer resource and speed up the numerical simulation, the numerical modeling is carried out stage by stage. In order to avoid the large pressure drop penalty in enhancing heat transfer, a slit fin is presented with the strip arrangement of “front coarse and rear dense” along the flow direction. The numerical simulation shows that, compared to the plain plate fin heat exchanger, the increase in the heat transfer in the slit fin heat exchanger is higher than that of the pressure drop, which proves the excellent performance of this slit fin. The fluid flow and heat transfer performance along the stages is also provided.     

Effects of Chemical Reaction and Double Dispersion on Non-Darcy Free Convection Heat and Mass Transfer

In this article, the effects of chemical reaction and double dispersion on non-Darcy free convection heat and mass transfer from semi-infinite, impermeable vertical wall in a fluid saturated porous medium are investigated. The Forchheimer extension (non-Darcy term) is considered in the flow equations, while the chemical reaction power–law term is considered in the concentration equation. The first order chemical reaction (n = 1) was used as an example of calculations. The Darcy and non-Darcy flow, temperature and concentration fields in this study are observed to be governed by complex interactions among dispersion and natural convection mechanisms. The governing set of partial differential equations were non-dimensionalized and reduced to a set of ordinary differential equations for which Runge–Kutta-based numerical technique were implemented. Numerical results for the detail of the velocity, temperature, and concentration profiles as well as heat transfer rates (Nusselt number) and mass transfer rates (Sherwood number) are presented in graphs.