并列圆柱绕流的格子Boltzmann数值模拟

采用非均匀不可压格子Boltzmann模型对低雷诺数下并列圆柱绕流进行了数值模拟,给出了数值计算结果,分析了间距g对圆柱尾流及升力、阻力的影响,并在此基础上得到了4种尾迹模式.此外,研究了流场的初始扰动对流动分岔现象的影响,发现在适当的扰动下可以很快得到同步同相的尾流.对Re=160和200下圆柱的升、阻力进行了对比,结果表明升力和阻力受间距g的影响大于雷诺数.     

Experimental and Theoretical Sensitivity Analysis of Turbulent Flow Past a Square Cylinder

The Siemens SGT-800 3^rd generation DLE burner fitted to an atmospheric combustion rig has been numerically investigated. Pure methane and methane enriched by 80 vol% hydrogen flames have been considered. A URANS (Unsteady Reynolds Averaged Navier-Stokes) approach was used in this study along with the k ? ω SST and the k ? ω SST-SAS models for the turbulence transport. The chemistry is coupled to the turbulent flow simulations by the use of a laminar flamelet library combined with a presumed PDF. The effect of the mesh density in the mixing and the flame region and the effect of the turbulence model and reaction rate model constant are first investigated for the methane/air flame case. The results from the k ? ω SST-SAS along with flamelet libraries are shown to be in excellent agreement with experimental data, whereas the k ? ω SST model is too dissipative and cannot capture the unsteady motion of the flame. The k ? ω SST-SAS model is used for simulation of the 80 vol% hydrogen enriched flame case without further adjusting the model constants. The global features of the hydrogen enrichment are very well captured in the simulations using the SST-SAS model. With the hydrogen enrichment the time averaged flame front location moves upstream towards the burner exit nozzle. The results are consistent with the experimental observations. The model captures the three dominant low frequency unsteady motion observed in the experiments, indicating that the URANS/LES hybrid model indeed is capable of capturing complex, time dependent, features such as an interaction between a PVC and the flame front.     

Unsteady Laminar to Turbulent Flow in a Spacer-Filled Channel

A combined numerical and experimental investigation has been carried out to study the flow behaviour in a spacer-filled channel, representative of those used in spiral-wound membrane modules. Direct numerical simulation and particle image velocimetry were used to investigate the fluid flow characteristics inside a 2 × 2 cell at Reynolds numbers that range between 100 and 1000. It was found that the flow in this geometry moves parallel to and also rotates between the spacer filaments and that the rate of rotation increases with Reynolds number. The flow mechanisms, transition process and onset of turbulence in a spacer-filled channel are investigated including the use of the velocity spectra at different Reynolds numbers. It is found that the flow is steady for Re < 200 and oscillatory at Re ～ 250 and increasingly unsteady with further increases in Re before the onset of turbulent flow at Re ～ 1000.     

Large-Eddy Simulation of Turbulent Flow Around a Finite-Height Wall-Mounted Square Cylinder Within a Thin Boundary Layer

The paper describes the results of a computational study of the auto-ignition of a fuel spray under Exhaust Gas Recirculation (EGR) conditions, a technique used to reduce the production of NO_x. Large Eddy Simulation (LES) is performed, and the stochastic field method is used for the solution of the joint sub-grid probability density function (pdf) of the chemical species and energy. The fuel spray is n-heptane, a diesel surrogate and its chemical kinetics are described by a reduced mechanism involving 22 species and 18 reaction steps. The method is applied to a constant volume combustion vessel able to reproduce EGR conditions by the ignition of a hot gas mixture previously introduced into the chamber. Once the prescribed conditions are reached the fuel is then injected. Different EGR conditions in terms of temperature and initial ambient chemical composition are simulated. The results are in good overall agreement with measurements both regarding the ignition delay times and the lift-off heights.     

Analysis of Thermal Flow in a Rotating Porous U-Turn Duct Using Lattice Boltzmann Method

The lattice Boltzmann method is carried out to investigate the heat transfer enhancement in a U-turn duct which is partially filled with a porous media. The porous layer is inserted at the core of the duct and is modeled using the Brinkman–Forchheimer assumptions. In order to validate the results, first a channel flow problem without any porous layer is compared with available data. Second, the porous Couette flow and partially porous channel flow are successfully compared with the studies of other researchers. Then, fluid flow in a clear U-turn duct is studied looking carefully at the velocity, curvature and rotation effects. Finally, the effects of porous layer thickness on the rate of heat transfer and pressure drop are investigated. Parametric studies are conducted to evaluate the effects of various parameters (i.e., Reynolds number, Darcy number, rotation number), highlighting their influences on the thermo-hydrodynamics behavior of the flow. The optimum values of porous layer thickness are presented for specific flow parameters.     

Supercritical Flow Over a Sill in an Open Channel

Experimental data on typical profiles of free surface and channel bottom pressure for a supercritical flow over a sill are reported. This flow is shown to have, along with the known critical depth, two other characteristic depths, one of which is at the channel exit to the atmosphere and the other determines conditions under which the disturbances propagate well upstream of the sill. The experimental data are compared with calculation results based on a mathematical model that incorporates turbulent mixing upon wave breaking.     

Dissipation of Turbulent Kinetic Energy in a Cylinder Wall Junction Flow

The subject of this study is the discussion of the dissipation of turbulent kinetic energy and its Reynolds number scaling in front of a wall-mounted cylinder. We employed highly resolved Large-Eddy Simulation and ensured that the computational grid was fine enough to resolve most of the scales. A perceptible fraction of the total dissipation is modeled. However, this fraction - about one third - is small enough so that the total dissipation suffers only marginally from some potential shortcomings of the turbulence model. Individual terms of the pseudo dissipation tensor and their Reynolds number scaling are discussed and compared. This tensor and thus the turbulent small scale structures are not isotropic at the Reynolds numbers investigated. Furthermore, the near-wall anisotropy under the horseshoe vortex is likely to persist to larger Reynolds numbers as it can be linked to a flapping of the near-wall layer. The turbulent length scale shows a strong spatial variability. In the region of the vortex system in the cylinder front, the distribution reveals a similar shape as the one of the turbulent kinetic energy and its amplitude is in the order of magnitude of the cylinder diameter. In contrast to the region dominated by the approach flow, the turbulent length scale is independent of the Reynolds number in the region dominated by the vortex system. Even though the flow investigated is in non-equilibrium, common a priori estimations and scalings of the Kolmogorov length scale based on macro scales give satisfying results.     

DNS of the Turbulent Channel Flow of a Dilute Polymer Solution

A direct numerical simulation of the turbulent channel flow of a dilute polymer solution has been performed in order to compare its turbulence statistics with those obtained in a Newtonian channel flow. The viscoelastic flow has been simulated by solving the whole set of continuity, momentum and constitutive equations for the six independent components of the extra-stress tensor induced by polymer addition. The Finitely Extensible Nonlinear Elastic dumbbell model was adopted in order to simulate a non-linear modulus of elasticity and a finite extendibility of the polymer macromolecules. Simulations were carried out under the narrow channel assumption at a Reynolds number of 169 based on the channel half height and on the friction velocity; they showed a significant reduction in drag, dependent on the influence of the elastic properties of the chains. A qualitative comparison with experiments at a higher Reynolds number has shown that the model here adopted is capable of reproducing all the main features of the polymer solution flow. Analysis of the turbulence statistics suggests that a dilute polymer solution can affect the intensity of the streamwise vortices, leading to an increase in the spacing between low speed streaks and eventually to a turbulent shear stress reduction.     

Passive Scalar in a Turbulent Channel Flow with Wall Velocity Disturbances

Direct Numerical Simulations (DNS) of a passive scalar in a turbulent channel flow with a normal velocity disturbance on the lower wall are presented for high and low Reynolds numbers. The aim is to reproduce the complex physics of turbulent rough flows without dealing with the geometric complexity. In addition, isothermal walls that cannot be easily assigned in an experiment, are considered. The paper explains the increase of heat transfer through the changes of the velocity and thermal structures. As in real rough flows, the transpiration produces an isotropization of the turbulence near the wall.     

Undular Jump in Open‐Channel Flow Over a Sill

The paper reports experimental data on a steady flow with a free surface over a rectangular sill at the bottom of a rectangular channel. The paper focuses on the parameter range in which a qualitative change in flow pattern accompanied by formation of an undular jump is observed. It is shown that transition from the bottom state the surface state of head and tail conjugation to occurs when the second critical depth is reached over the downstream edge of the sill. Conditions of transition from a free state to a submerged state of flow over the sill are considered; conditions of validity of the assumption that the first critical depth is reached over the sill are discussed.     

Permeability Effects of Turbulent Flow Through a Porous Insert in a Backward-Facing-Step Channel

In the present work, a k– model, based on the work of Lee and Howell (Proceedings of the ASME-JSME Thermal Engineering Hawaii, 1987), is rigorously derived based on time average of spatially averaged Navier–Stokes equations. The model is then employed to solve for a flow in a backward-facing step channel with a porous insert. The numerical solver is modified from the STREAM code (Lien and Leschziner, Comput. Meth. Appl. Mech. Eng. 114 (1994a) 123–148), and it has been validated against the experimental data of Seegmiller and Driver (AIAA Journal 23 (1985) 163–171). The code is then used to perform simulation for cases with a porous insert. The resistance of the porous insert can be altered by changing its permeability (), Forchheimers constant (F), or thickness (b). The goal is to examine the influence of each parameter on the resulting flow and turbulent kinetic energy (k) distributions. It is discovered that, by increasing the resistance of the insert, flow eventually enters a transitional regime towards relaminarization. This is due to the contribution of Darcys and Forchheimers terms in the governing equations, and modifying these two terms changes the levels of P_k and, hence, k and . Generally speaking, lowering or raising F results in a greater suppression of P_k than , causing the flow to relaminarize. Meanwhile, if the pore size is reasonably large to sustain turbulence within the porous media, increasing b reduces but does not eliminate the turbulent activity in the porous insert.     

Numerical Simulation of a Turbulent Flow in a Channel with Surface Mounted Cubes

In this paper we report on a fourth-order, spectro-consistent simulation of a complex turbulent flow. A spatial discretization of a convection-diffusion equation is termed spectro-consistent if the spectral properties of the convective and diffusive operators are preserved, i.e. convection skew-symmetric; diffusion symmetric positive definite. We consider a fully developed flow in a channel, where a matrix of cubes is placed at a wall of the channel. The Reynolds number (based on the channel width and the mean bulk velocity) is equal to Re = 13,000. The three-dimensional flow around the surface mounted cubes has served at a test case at the 6th ERCOFTAC/IAHR/COST workshop on refined flow modeling (Delft, June 1997). Here, mean velocity profiles as well as Reynolds stresses at various locations in the channel have been computed without using any turbulence models. The results agree well with the available experimental data.     

Large-Eddy Simulation of Dispersion from Line Sources in a Turbulent Channel Flow

Large-eddy simulations of the dispersion from scalar line sources at various locations within a fully developed turbulent channel flow at Re = uh/ν = 10,400 are presented. Both mean and fluctuating scalar quantities are compared with those from the single available set of experimental data (Lavertu and Mydlarski, J Fluid Mech 528:135–172, 2005) and differences are highlighted and discussed. The results are also discussed in the context of scalar dispersion in other kinds of turbulent flows, e.g. homogeneous shear-flow. Initial computations at a much lower Reynolds number are also reported and compared with the two available direct numerical simulation data sets.     

Use of a Modified Langevin Equation to Describe Turbulent Dispersion of Fluid Particles in a Channel Flow

A stochastic method to represent the positions and velocities of fluid particles in a nonhomogeneous turbulence was pursued. Spatially varying Lagrangian time scales obtained from direct numerical simulations of turbulent flow in a channel and spatially varying joint Gaussian forcing functions were incorporated into a Langevin equation. The model was tested by comparing calculations of the dispersions and velocities of particles originating from point sources with experiments carried out in a DNS of fully-developed turbulent flows in a channel at Re_τ = 150 and 300. The model captured the dispersions, mean velocities and moments of the velocity fluctuations very well. The use of jointly Gaussian (rather than uncorrelated Gaussian) forcing functions greatly improved the calculation of dispersions and mean velocities in the stream wise direction, as well as turbulence quantities which include stream wise velocity fluctuations. The condition of well-mixedness for the model was verified by considering the behavior of a uniform distribution of sources. This revised version was published online in July 2006 with corrections to the Cover Date.     

Pulsatory Turbulent Compressible Fluid Flow and Propagation of Pressure Waves in a Channel

The properties of the damping coefficient and phase velocity of propagation of small-amplitude pressure waves as functions of the oscillation frequency are investigated for the turbulent flow of a weakly compressible fluid in a circular pipe. The wall friction is found by solving numerically the equation of motion and the relaxation equations for the turbulent shear stress and viscosity which provide the basis for a turbulent transfer model developed for unsteady conditions. The properties are explained in terms of an analysis of the calculated data on turbulent transfer. The results obtained are compared with experiments.     

Turbulent Flow and Dispersion of Inertial Particles in a Confined Jet Issued by a Long Cylindrical Pipe

In this work we examine first the flow field of a confined jet produced by a turbulent flow in a long cylindrical pipe issuing in an abrupt angle diffuser. Second, we examine the dispersion of inertial micro-particles entrained by the turbulent flow. Specifically, we examine how the particle dispersion field evolves in the multiscale flow generated by the interactions between the large-scale structures, which are geometry dependent, with the smaller turbulent scales issued by the pipe which are advected downstream. We use Large-Eddy-Simulation (LES) for the flow field and Lagrangian tracking for particle dispersion. The complex shape of the domain is modelled using the immersed-boundaries method. Fully developed turbulence inlet conditions are derived from an independent LES of a spatially periodic cylindrical pipe flow. The flow field is analyzed in terms of local velocity signals to determine spatial coherence and decay rate of the coherent K–H vortices and to make quantitative comparisons with experimental data on free jets. Particle dispersion is analyzed in terms of statistical quantities and also with reference to the dynamics of the coherent structures. Results show that the particle dynamics is initially dominated by the Kelvin–Helmholtz (K–H) rolls which form at the expansion and only eventually by the advected smaller turbulence scales.     

Passive Scalar Transport in a Turbulent Cylinder Wake in the Presence of a Downstream Cylinder

This work aims to investigate how the presence of a downstream cylinder affects the passive scalar transport in a cylinder wake. The wake was generated by two tandem brass circular cylinders of the same diameter (d). The cylinder centre-to-centre spacing L/d was 1.3, 2.5 and 4.0, respectively, covering the three typical flow regimes of this flow. The upstream cylinder was slightly heated. Measurements were conducted at x/d= 10 and Re (≡ dU _∞/ν, where U _∞ is the free-stream velocity and ν is the kinematic viscosity of fluid) = 7000. A three-wire probe consisting of an X-wire and a cold wire was used to measure the velocity and temperature fluctuations, while an X-wire provided a phase reference. The phase-averaged velocity vectors and vorticity display single vortex street behind the downstream cylinder, irrespective of the flow regimes. However, the detailed flow structure exhibits strong dependence on L/d in terms of the Strouhal number, the vortex strength and its downstream evolution. This naturally affects passive scalar transport. The coherent and incoherent heat flux vectors show significant variation for different L/d.     

Reconstruction of Scalar Source Intensity Based on Sensor Signal in Turbulent Channel Flow

We consider a turbulent channel flow, where a scalar point source with a time-harmonic intensity releases a substance that can be modeled as a passive scalar. With the source location known, our objective is to estimate the time history of the source intensity based on sensor measurements at different locations downstream of the source by adopting an adjoint approach. It is shown that the proposed algorithm reproduces the original coherent sinusoidal wave of the scalar source accurately from the chaotic scalar signals measured by our sensors. By systematically changing the source-sensor distance and the pulsation frequency of the source, we clarify how these two factors affect the estimation accuracy. The proposed scheme is also applicable to estimation with multiple sensors. We demonstrate that increasing the number of sensors improves the estimation greatly when the scalar is released from a source away from the wall, where large-scale eddies dominate the scalar dispersion. In contrast, the estimation performance remains poor even with multiple sensors when the scalar source is located near the wall, where the source information is quickly lost due to the strong turbulence activity and the scalar diffusion in the near-wall region.     

Pore-Scale Simulations of Gas Displacing Liquid in a Homogeneous Pore Network Using the Lattice Boltzmann Method

A lattice Boltzmann high-density-ratio model, which uses diffuse interface theory to describe the interfacial dynamics and was proposed originally by Lee and Liu (J Comput Phys 229:8045–8063, 2010), is extended to simulate immiscible multiphase flows in porous media. A wetting boundary treatment is proposed for concave and convex corners. The capability and accuracy of this model is first validated by simulations of equilibrium contact angle, injection of a non-wetting gas into two parallel capillary tubes, and dynamic capillary intrusion. The model is then used to simulate gas displacement of liquid in a homogenous two-dimensional pore network consisting of uniformly spaced square obstructions. The influence of capillary number (Ca), viscosity ratio ( $M$ M ), surface wettability, and Bond number (Bo) is studied systematically. In the drainage displacement, we have identified three different regimes, namely stable displacement, capillary fingering, and viscous fingering, all of which are strongly dependent upon the capillary number, viscosity ratio, and Bond number. Gas saturation generally increases with an increase in capillary number at breakthrough, whereas a slight decrease occurs when Ca is increased from $8.66\times 10^{-4}$ 8.66 × 10 - 4 to $4.33\times 10^{-3}$ 4.33 × 10 - 3 , which is associated with the viscous instability at high Ca. Increasing the viscosity ratio can enhance stability during displacement, leading to an increase in gas saturation. In the two-dimensional phase diagram, our results show that the viscous fingering regime occupies a zone markedly different from those obtained in previous numerical and experimental studies. When the surface wettability is taken into account, the residual liquid blob decreases in size with the affinity of the displacing gas to the solid surface. Increasing Bo can increase the gas saturation, and stable displacement is observed for $Bo>1$ B o > 1 because the applied gravity has a stabilizing influence on the drainage process.