Turbulent transport modelling of separating and reattaching shear flows

The improvement of capabilities for computer simulation of turbulent recirculating flows was investigated. Attention has been limited to two dimensional flows and principally to statistically stationary motion. Improvement of turbulence modeling explored the treatment of the near wall sublayer and of the exterior fully turbulent region, working within the framework of turbulence closures requiring the solution of transport equations for the turbulence energy and its dissipation rate. The work on the numerical procedure, based on the Gosman-Pun program TEACH, addressed the problems of incorporating the turbulence model as well as the extension to time dependent flows, the incorporation of a third order approximation of convective transport, and the treatment of non-orthogonal boundaries

Combustion of hydrogen-air jets in local chemical equilibrium: A guide to the CHARNAL computer program

A guide to a computer program, written in FORTRAN 4, for predicting the flow properties of turbulent mixing with combustion of a circular jet of hydrogen into a co-flowing stream of air is presented. The program, which is based upon the Imperial College group's PASSA series, solves differential equations for diffusion and dissipation of turbulent kinetic energy and also of the R.M.S. fluctuation of hydrogen concentration. The effective turbulent viscosity for use in the shear stress equation is computed. Chemical equilibrium is assumed throughout the flow

Modelling windage power loss from an enclosed spur gear

Within a gearbox the majority of transmission losses can be attributed to bearing losses, meshing losses, or losses due to windage/churning. In this paper the commercial computational fluid dynamics (CFD) code Fluent 6.2.16 is applied in a two-dimensional study of windage power loss (WPL) from a single spur gear rotating in air. By comparing CFD data to published experimental data appropriate grid density and modelling parameters are identified. The model is used to investigate how peripheral shrouding affects WPL and whether WPL can be reduced through minor modifications to tooth tip geometry. Non-dimensional shroud spacings (ratio of gap to gear PCD) of between 0.005 and 0.05 were investigated at shaft speeds between 5000 and 20 000 r/min. Although CFD data compared reasonably well to experimental data, trends were not reproduced and an optimum shroud could not be identified. A full three-dimensional study is recommended. Modifying the tooth tip by adding a small chamfer on the leading edge reduced WPL by ~6 per cent. A small fillet increased total WPL by a similar amount suggesting that WPL may increase as a gear wears. This preliminary study suggests further work in this area would be beneficial

Phenomenology of Wall Bounded Newtonian Turbulence

We construct a simple analytic model for wall-bounded turbulence, containing only four adjustable parameters. Two of these parameters characterize the viscous dissipation of the components of the Reynolds stress-tensor and other two parameters characterize their nonlinear relaxation. The model offers an analytic description of the profiles of the mean velocity and the correlation functions of velocity fluctuations in the entire boundary region, from the viscous sub-layer, through the buffer layer and further into the log-layer. As a first approximation, we employ the traditional return-to-isotropy hypothesis, which yields a very simple distribution of the turbulent kinetic energy between the velocity components in the log-layer: the streamwise component contains a half of the total energy whereas the wall-normal and the cross-stream components contain a quarter each. In addition, the model predicts a very simple relation between the von-K\'arm\'an slope $\kappa$ and the turbulent velocity in the log-law region $v^+$ (in wall units): $v^+=6 \kappa$. These predictions are in excellent agreement with DNS data and with recent laboratory experiments.Comment: 15 pages, 11 figs, included, PRE, submitte

Journal off Fluids Sneering PSL-An Economical Approach to the Numerical Analysis of Near-Wall, Elliptic Flow

The paper points out that, in the numerical computation of elliptic or three-dimensional turbulent flows, the neglect of pressure-variations across the very thin viscosity affected region near the wall allows a fine-grid analysis of this sublayer without prohibitive penalties in core or computational time. The scheme has been successfully applied to the threedimensional flow around a U-bend

Influence of reaction atmosphere (H2O, N2, H2, CO2, CO) on fluidized-bed fast pyrolysis of biomass using detailed tar vapor chemistry in computational fluid dynamics

Secondary pyrolysis in fluidized bed fast pyrolysis of biomass is the focus of this work. A novel computational fluid dynamics (CFD) model coupled with a comprehensive chemistry scheme (134 species and 4169 reactions, in CHEMKIN format) has been developed to investigate this complex phenomenon. Previous results from a transient three-dimensional model of primary pyrolysis were used for the source terms of primary products in this model. A parametric study of reaction atmospheres (H2O, N2, H2, CO2, CO) has been performed. For the N2 and H2O atmosphere, results of the model compared favorably to experimentally obtained yields after the temperature was adjusted to a value higher than that used in experiments. One notable deviation versus experiments is pyrolytic water yield and yield of higher hydrocarbons. The model suggests a not overly strong impact of the reaction atmosphere. However, both chemical and physical effects were observed. Most notably, effects could be seen on the yield of various compounds, temperature profile throughout the reactor system, residence time, radical concentration, and turbulent intensity. At the investigated temperature (873 K), turbulent intensity appeared to have the strongest influence on liquid yield. With the aid of acceleration techniques, most importantly dimension reduction, chemistry agglomeration, and in-situ tabulation, a converged solution could be obtained within a reasonable time (∼30 h). As such, a new potentially useful method has been suggested for numerical analysis of fast pyrolysis

PDF model based on Langevin equation for polydispersed two-phase flows applied to a bluff-body gas-solid flow,

The aim of the paper is to discuss the main characteristics of a complete theoretical and numerical model for turbulent polydispersed two-phase flows, pointing out some specific issues. The theoretical details of the model have already been presented [Minier and Peirano, Physics Reports, Vol. 352/1-3, 2001 ]. Consequently, the present work is mainly focused on complementary aspects, that are often overlooked and that require particular attention. In particular, the following points are analysed : the necessity to add an extra term in the equation for the velocity of the fluid seen in the case of twoway coupling, the theoretical and numerical evaluations of particle averages and the fulfilment of the particle mass-continuity constraint. The theoretical model is developed within the PDF formalism. The important-physical choice of the state vector variables is first discussed and the model is then expressed as a stochastic differential equation (SDE) written in continuous time (Langevin equations) for the velocity of the fluid seen. The interests and limitations of Langevin equations, compared to the single-phase case, are reviewed. From the numerical point of view, the model corresponds to an hybrid Eulerian/Lagrangian approach where the fluid and particle phases are simulated by different methods. Important aspects of the Monte Carlo particle/mesh numerical method are emphasised. Finally, the complete model is validated and its performance is assessed by simulating a bluff-body case with an important recirculation zone and in which two-way coupling is noticeable.Comment: 23 pages, 10 figure