Experimental and numerical analysis of cylindrical turbulators in a double pipe heat exchanger under turbulent flow conditions

In the present work, the use of cylindrical turbulators in a double pipe heat exchanger has been investigated. Cylindrical fin type of turbulators has been placed circumferentially separated by 90° on the outer side of an inner pipe at a regular pitch. Experimental studies were undertaken for different air flow rates in a turbulent regime whose Reynolds number range between 2500 and 10000. Heat transfer characteristics like Nu and friction factor have been experimentally determined. Parametric studies were conducted by changing the pitch and also the orientation of the turbulators. Nu and friction factor were found to increase as the pitch is reduced. A model with alternatively changed orientation outperformed others by exhibiting highest Nu and reduced friction factor.     

A model for crack formation during active solid pyrolysis of a char-forming solid: crack patterns; surface area generation; volatile mass efflux

A model is developed for the formation and propagation of cracks in a material sample that is heated at its top surface, pyrolyses, and then thermally degrades to form char. In this work the sample is heated uniformly over its entire top surface by a hypothetical flame (a heat source). The pyrolysis mechanism is described by a one-step overall reaction that is dependent nonlinearly on the temperature (Arrhenius form). Stresses develop in response to the thermal degradation of the material by means of a shrinkage strain caused by local mass loss during pyrolysis. When the principal stress exceeds a prescribed threshold value, the material forms a local crack. Cracks are found to generally originate at the surface in response to heating, but occasionally they form in the bulk, away from ever-changing material boundaries. The resulting cracks evolve and form patterns whose characteristics are described. Quantities examined in detail are: the crack spacing in the pyrolysis zone; the crack length evolution; the formation and nature of crack loops which are defined as individual cracks that have joined to form loops that are disconnected from the remaining material; the formation of enhanced pyrolysis area; and the impact of all of the former quantities on mass flux. It is determined that the mass flux from the sample can be greatly enhanced over its nominal (non-cracking) counterpart. The mass efflux profile qualitatively resembles those observed in Cone Calorimeter tests.     

AGE method simulations of a turbulent far‐wake compared to spectral DNS

Turbulent flow simulation methods based on finite differences are attractive for their simplicity, flexibility and efficiency, but not always for accuracy or stability. This paper demonstrates that a good compromise is possible with the advected grid explicit (AGE) method. Starting from the same initial field as a previous spectral DNS, AGE method simulations of a planar turbulent wake were carried out as DNS, and then at three levels of reduced resolution. The latter cases were in a sense large‐eddy simulations (LES), although no specific sub‐grid‐scale model was used. Results for the two DNS methods, including variances and power spectra, were very similar, but the AGE simulation required much less computational effort. Small‐scale information was lost in the reduced resolution runs, but large‐scale mean and instantaneous properties were reproduced quite well, with further large reductions in computational effort. Quality of results becomes more sensitive to the value chosen for one of the AGE method parameters as resolution is reduced, from which it is inferred that the numerical stability procedure controlled by the parameter is acting in part as a sub‐grid‐scale model. Copyright © 2002 John Wiley & Sons, Ltd.     

Compulsorily Phase Locked Interferometer Using Orthogonally Polarized Light Source of a Modulated LD with High Extinction Ratio

This paper describes a compulsorily phase locked differential interferometer using an orthogonally polarized light source of a modulated LD with high extinction ratio to reduce non-linearity of the interferometer caused by polarization cross-talk. The current modulated LD is used as a light source to make the interferometer compact and for the scanning phase of the interferometer. The interferometer is operated compulsorily at the maximum inclination point of the fringe intensity curve by fringe scanning and an electric system. A Wollaston prism of high extinction ratio (50 dB) is used to combine the polarizing beams and to make the polarization cross-talk very small. In one light source the polarized output beams are on the same propagation axis; in the other they have a small crossing angle (2.5 mrad ∼ 10 mrad) to completely exclude non-linearity of the interferometer causded by polarization cross-talk. Using jets of a gas mixture of nitrogen and ethylene, this interferometer was demonstrated to be useful in detecting the photothermal effect of a photothermal velocimeter under phase fluctuation in a turbulent flow.     

Numerical analysis of turbulent structure in compound meandering open channel by algebraic Reynolds stress model

The turbulent flow in a compound meandering channel with a rectangular cross section is one of the most complicated turbulent flows, because the flow behaviour is influenced by several kinds of forces, including centrifugal forces, pressure‐driven forces and shear stresses generated by momentum transfer between the main channel and the flood plain. Numerical analysis has been performed for the fully developed turbulent flow in a compound meandering open‐channel flow using an algebraic Reynolds stress model. The boundary‐fitted coordinate system is introduced as a method for coordinate transformation in order to set the boundary conditions along the complicated shape of the meandering open channel. The turbulence model consists of transport equations for turbulent energy and dissipation, in conjunction with an algebraic stress model based on the Reynolds stress transport equations. With reference to the pressure–strain term, we have made use of a modified pressure–strain term. The boundary condition of the fluctuating vertical velocity is set to zero not only for the free surface, but also for computational grid points next to the free surface, because experimental results have shown that the fluctuating vertical velocity approaches zero near the free surface. In order to examine the validity of the present numerical method and the turbulent model, the calculated results are compared with experimental data measured by laser Doppler anemometer. In addition, the compound meandering open channel is clarified somewhat based on the calculated results. As a result of the analysis, the present algebraic Reynolds stress model is shown to be able to reasonably predict the turbulent flow in a compound meandering open channel. Copyright © 2005 John Wiley & Sons, Ltd.     

CFD在燃煤细粒子凝聚过程中的应用

燃煤产生的细粒子富集了大量的有毒痕量元素,对大气环境和人类健康造成严重危害。本文针对燃煤细粒子的形成过程,在计算流体力学(CFD)的基础上,结合气溶胶动力学理论,模拟了烟气中细粒子在圆柱形流场中由于碰撞而凝聚的过程,并通过计算结果分析细粒子颗粒特征参量(速度、质量、直径和颗粒数目)随流场的变化趋势,为深入研究燃煤细粒子的形成演化机理提供理论基础。     

A numerical and experimental investigation of the interactions between a non‐uniform planar array of incompressible free jets

The interaction between multiple incompressible air jets has been studied numerically and experimentally. The numerical predictions have been first validated using experimental data for a single jet configuration. The spreading features of five unequal jets in the configuration of one larger central jet surrounded by four smaller equi‐distant jets, have been studied, for different lateral spacing ratios of 1.5, 2.0 and 2.5 and a central jet Reynolds number of 1.24×105 (corresponding to a Mach number of 0.16). Flow of five equal jets has also been simulated, for the sake of comparison. The jet interactions commence at an axial distance of about 3–4 diameters and complete by an axial distance of about 10 diameters for the lowest spacing ratio of 1.5. For larger spacing ratios, the length required for the start and completion of jet interaction increase. Peripheral jets bend more towards the central jet and merge at a smaller distance, when their sizes are smaller than that of the central jet. The entrainment ratio for multiple jets is higher than that for a single jet. Excellent agreement is observed between the experimental data and theoretical predictions for both mean flow field and turbulent quantities, at regions away from the jet inlet. The potential core length and initial jet development, however, are not predicted very accurately due to differences in the assumed and actual velocity profiles at the jet inlet. Copyright © 2004 John Wiley & Sons, Ltd.     

A time‐accurate pseudo‐compressibility approach based on an unstructured hybrid finite volume technique applied to unsteady turbulent premixed flame propagation

A preconditioning approach based on the artificial compressibility formulation is extended to solve the governing equations for unsteady turbulent reactive flows with heat release, at low Mach numbers, on an unstructured hybrid grid context. Premixed reactants are considered and a flamelet approach for combustion modelling is adopted using a continuous quenched mean reaction rate. An overlapped cell‐vertex finite volume method is adopted as a discretisation scheme. Artificial dissipation terms for hybrid grids are explicitly added to ensure a stable, discretised set of equations. A second‐order, explicit, hybrid Runge–Kutta scheme is applied for the time marching in pseudo‐time. A time derivative of the dependent variable is added to recover the time accuracy of the preconditioned set of equations. This derivative is discretised by an implicit, second‐order scheme. The resulting scheme is applied to the calculation of an infinite planar (one‐dimensional) turbulent premixed flame propagating freely in reactants whose turbulence is supposed to be frozen, homogeneous and isotropic. The accuracy of the results obtained with the proposed method proves to be excellent when compared to the data available in the literature. Copyright © 2004 John Wiley & Sons, Ltd.     

LES of turbulent heat transfer: proper convection numerical schemes for temperature transport

Large eddy simulations of two basic configurations (decay of isotropic turbulence, and the academic plane channel flow) with heat transfer have been performed comparing several convection numerical schemes, in order to discuss their ability to evaluate temperature fluctuations properly. Results are compared with the available incompressible heat transfer direct numerical simulation data. It is shown that the use of regularizing schemes (such as high order upwind type schemes) for the temperature transport equation in combination with centered schemes for momentum transport equation gives better results than the use of centred schemes for both equations. Copyright © 2004 John Wiley & Sons, Ltd.