Cloud dispersion models

A new technique has been developed to compute mean and fluctuating concentrations in complex turbulent flows. An initial distribution of material is discretised into any small clouds which are advected by a combination of the mean flow and large scale turbulence. The turbulence can be simulated either by Kinematic Simulation or by a stochastic model for the motion of each cloud centroid. The clouds also diffuse relative to their centroids; the statistics for this are obtained from a separate calculation of the growth of individual clouds in small scale turbulence, generated by Kinematic Simulation. The ensemble of discrete clouds is periodically rediscretised, to limit the size of the small clouds and prevent overlapping. The model is illustrated with simulations of dispersion in uniform flow and in a coastal flow, and the results are compared with analytic, steady state solutions where available.     

Diffusion of macromolecular solutions in the turbulent boundary layer of a cylindrical pipe — Evolution of wall concentration

Summary In this paper, we are presenting a model of the evolution of the wall concentration of a macromolecular solution (PEO) annularly injected in a cylindrical pipe in a turbulent flow. This model valid for all diffusion zones is based on mathematical and physical considerations and proves to be in good agreement with the experimental data.

---

Zusammenfassung Es wird ein Modell der Entwicklung der Wandkonzentration einer makromolekularen Lösung (PÄO) vorgestellt, die in einem wandnahen Ringspalt in die turbulente Strömung durch ein zylindrisches Rohr injiziert worden ist. Dieses für alle Diffusionszonen gültige Modell basiert auf mathematischen und physikalischen Betrachtungen und erweist sich für die Beschreibung der experimentellen Daten als gut geeignet.

---

C _w wall concentration - C _i initial concentration before injection - L ₀ distance from the slot at which the wall concentration drops toe ^-1 of its original value - L _IT ,L _IF ,L _F characteristic lengths - L _I length scale of the second region - x downstream distance from the source - n _I ,n _T ,n _F characteristic exponents - K ₀,K _I ,K _F characteristic constants - ln natural logarithm - q _i flow rate of injection - Q _T flow rate - C _j =C _{i · q i} /Q _T concentration in homogeneous medium - A, B, C, m constants - p andq annex variables - Re Reynolds number With 7 figures     

A von Neumann–Smagorinsky turbulent transport model for stratified shear flows

A simple subgrid turbulent diffusion model based on an analogy to the von Neumann–Richtmyer artificial viscosity is explored for use in modelling mixing in turbulent stratified shear flow. The model may be more generally applicable to multicomponent turbulent hydrodynamics and to subgrid turbulent transport of momentum, composition and energy. As in the case of the von Neumann artificial viscosity and many subgrid-scale models for large-eddy simulation, the turbulent diffusivity explicitly depends on the grid size and is not based on a quantitative model of the unresolved turbulence. In order to address the issue that it is often not known a priori when and where a flow will become turbulent, the turbulent diffusivity is set to zero when the flow is expected to be stable on the basis of a Richardson/Rayleigh–Taylor stability criterion, in analogy to setting the von Neumann artificial viscosity to zero in expanding flows. One-dimensional predictions of this model applied to a simple shear flow configuration are compared to those obtained using a K–ε model. The density and velocity profiles predicted by both models are shown to be very similar.     

A silicon transducer for the determination of wall-pressure fluctuations in turbulent boundary layers

Small and sensitive silicon sensors for turbulent wall-pressure fluctuation measurements have been designed and fabricated using microelectronic technology. For the detection of the pressure fluctuations piezoresistive gauges are placed on a diaphragm and the resistance of these gauges depends on the stresses in the diaphragm. For the determination of the performance of these pressure transducers comparisons with conventional microphones were carried out in a well-defined two-dimensional boundary layer. Power spectra from the silicon pressure transducer revealed a slope in the inertial sublayer corresponding approximately to the 5/3-law of Kolmogorov, and the normalized wall pressure fluctuations agreed well with other direct measurements.     

An assessment of particle methods for approximating anisotropic dispersion

We derive a smoothed particle hydrodynamics (SPH) approximation for anisotropic dispersion that only depends upon the first derivative of the kernel function and study its numerical properties. In addition, we compare the performance of the newly derived SPH approximation versus an implementation of the particle strength exchange (PSE) method and a standard finite volume method for simulating multiple scenarios defined by different combinations of physical and numerical parameters. We show that, for regularly spaced particles, given an adequate selection of numerical parameters such as kernel function and smoothing length, the new SPH approximation is comparable with the PSE method in terms of convergence and accuracy and similar to the finite volume method. On other hand, the performance of both particle methods (SPH and PSE) decreases as the degree of disorder of the particle increases. However, we demonstrate that in these situations the accuracy and convergence properties of both particle methods can be improved by an adequate choice of some numerical parameters such as kernel core size and kernel function. Copyright © 2012 John Wiley & Sons, Ltd.     

On diffusion velocity and the mass conservation equation for components

The mass migration velocity(absolute velocity)of component i in a multicomponent flow is equal to the convection velocity(frame velocity)plus the diffusion velocity(relativevelocity).The diffusion velocity as well as the corresponding diffusion coefficient depends on how the convection velocity is adopted.In turbulent flow,the mass migration velocity of component i is(?)(mass-weighted time average velocity).The diffusion velocity(-a)consists of turbulent diffusionvelocity(?)and molecular diffusion velocity(?)(?is the simple time average velocity of component i and a is a certain convection velocity).So,the part of turbulent diffusion velocity is independent of what convection velocity is taken.In the mass conservation equation for component i,the expression for the diffusion term on its right-hand side will change when the convection velocity on its left-hand side changes.In turbulent flow,there could be no diffusion terms,or a turbulent diffusion term only,or both the turbulent and molecular diffusion     

APPLICATION OF CONVECTION-DIFFUSION EQUATION TO THEANALYSES OF CONTAMINATION BETWEEN BATCHESIN MULTI-PRODUCTS PIPELINE TRANSPORT

I.IntroductionIntilenearfuture,ourcountry'sproductspipelinetransportwillhavebeendevelopedgreatly.Themeansbywhichmanyrefinedpetroleumproductsaretransportedthroughasinglepipelineiscalledbatchingtransport.Thismeansispossessedofadvantagesthatavarietyofrefinedpetroleumproductscanbetrallsportedthroughasinglepipelineandinvestmentcanbereduced.[-lowever,tile1ila-jorproblemresultedfronlthebatchingtrallsportisthecontamillationbetweenbatcllcs,whichmaydeterioratethequalityofproductstransported.Therefore,i…     

用火焰面模型模拟甲烷/空气湍流射流扩散火焰

以层流对撞扩散火焰为基础，利用层流火焰面模型(laminar flamelet model)的方法生成层流火焰面数据库，分别采用预先设定的几率密度函数(propabality density function, PDF)模型和混合物分数-湍流频率的联合几率密度函数输运模型，将火焰面方法应用于甲烷/空气湍流射流扩散火焰结构的模拟计算中.两个模型的计算结果和实验结果进行了比较和分析.     

Analysis of Streamwise Velocity Fluctuations in Turbulent Pipe Flow with the Use of an Ultrasonic Doppler Flowmeter

Data collected from several studies of experimental and numerical nature in wall-bounded turbulent flows and in particular in internal flows (channel and pipe flows, Mochizuki and Nieuwstadt [1]) at different Reynolds numbers R ⁺(Ru _*/ν), indicate that: (i) the peak of the rms-value (normalized by u _*) of the streamwise velocity fluctuations (σ _u ⁺|_peak) is essentially independent of the Reynolds number, (ii) the position of the rms peak value (y ⁺|_peak) is weakly dependent of the Reynolds number, (iii) the skewness of the streamwise velocity fluctuations (S _u ) is close to zero at the position in which the variance has its peak. A series of measurements of streamwise velocity fluctuations has been performed in turbulent pipe flow with the use of an Ultrasonic Doppler Velocimeter and our results support those reported in [1]. This revised version was published online in July 2006 with corrections to the Cover Date.     

Large eddy simulation of turbulent statistical and transport properties in stably stratified flows

Three dimensional large eddy simulation （LES） is performed in the investigation of stably stratified turbulence with a sharp thermal interface. Main results are focused on the turbulent characteristic scale, statistical properties, transport properties, and temporal and spatial evolution of the scalar field. Results show that the buoyancy scale increases first, and then goes to a certain constant value. The stronger the mean shear, the larger the buoyancy scale. The overturning scale increases with the flow, and the mean shear improves the overturning scale. The flatness factor of temperature departs from the Gaussian distribution in a fairly large region, and its statistical properties are clearly different from those of the velocity fluctuations in strong stratified cases. Turbulent mixing starts from small scale motions, and then extends to large scale motions.     

Schlieren measurement of turbulent structure in a diffusion flame

The regular and random mixing structures in a turbulent diffusion flame were investigated using the quantitative, dynamic crossed-beam schlieren method. Evidence was found close to the nozzle relating to the vortexlike structure of eddies surrounding the central fuel jet flow. The observations also make possible resolution of turbulent intensity, scales, convection, and spectra within the diffusion flame without the use of seeding or intrusion of measuring probes. It is found that length scales and other turbulence parameters in the diffusion flame progressively revert to values similar to those expected and observed in scalar passive mixing as the combustion reaction intensity reduces with axial distance from the nozzle system.     

Fixed point index of ultimately compact set-valued mappings in Hausdorff locally convex spaces and its applications

The author defines a concept of fixed point index of ultimately compact set-valued mappings in Hausdorff locally convex spaces. Using this concept, the author establishes several nonzero fixed point theorems of set-valued -condensing mappings. These theorems extend some known results in [1,2,7,8,9].Projects Supported by the Science Fund of the Chinese Academy of Sciences     

Dispersion modeling by kinematic simulation: Cloud dispersion model

A new technique has been developed to compute mean and fluctuating concentrations in complex turbulent flows (tidal current near a coast and deep ocean). An initial distribution of material is discretized into any small clouds which are advected by a combination of the mean flow and large scale turbulence. The turbulence can be simulated either by kinematic simulation (KS) or direct numerical simulation. The clouds also diffuse relative to their centroids; the statistics for this are obtained from a separate calculation of the growth of individual clouds in small scale turbulence, generated by KS. The ensemble of discrete clouds is periodically re-discretized, to limit the size of the small clouds and prevent overlapping. The model is illustrated with simulations of dispersion in uniform flow, and the results are compared with analytic, steady state solutions. The aim of this study is to understand how pollutants disperses in a turbulent flow through a numerical simulation of fluid particle motion in a random flow field generated by Fourier modes. Although this homogeneous turbulent is rather a “simple” flow, it represents a building block toward understanding pollutant dispersion in more complex flow. The results presented here are preliminary in nature, but we expect that similar qualitative results should be observed in a genuine turbulent flow.     

Experimental study on turbulent features in the negative transport region of asymmetric plane channel flow

Turbulent features of streamwise and vertical components of velocity in the negative transport region of asymmetric plane channel flow have been studied experimentally in details. Experiments show that turbulent fluctuations in negative transport region are suppressed, and their probability distributions are far from Gaussian. Besides, the skewness factors attain their negative maxima at the position of the maximum mean velocity, whereas the flatness factors attain their positive maxima at the same position. The project supported by the National Natural Science Foundation of China (19872043)     

A METHOD OF DISCRETE DISPLACEMENTS FOR MODELLING THE CONTRIBUTION OF FINITE-SIZE ORGANIZED FLUID STRUCTURES TO TRANSPORT PROCESSES IN TURBULENT PIPE FLOW

A method of modelling the contribution of finite-size organized streams and fluid structures to the processes of turbulent transport is presented for the example of developed turbulent pipe flow. The method is applied to construct the turbulent length (L) and eddy viscosity coefficient (ν_t) employed to compute the average characteristics of the flow. The average effects of action of these organized fluid structures and streams are modelled as the final results of discrete displacements of certain model turbulent signals between nodes associated in pairs as well as the results of effective discrete displacements of these pairs. The displacement of information about the organization of two nodes into a pair identifies the displacement of the pair. These nodes constitute a network whose parameters have been established a priori analytically by considering a sequence of model turbulent lengths scaled with their distance to the wall. The model turbulent signals are evaluated at respective discrete nodes with the help of a certain finite difference turbulence model closed by L and ν_t and realized on an appropriate numerical grid. The non-uniform grid spacing has been related unambiguously in a rational way to the sequence of model turbulent lengths. Results elucidating specific features of this discrete modelling, particularly its differences from the continuous approach, are presented. Good agreement of the results with available experimental data is demonstrated. The average characteristics of the flow structure predicted for a wide range of Reynolds number (Re) are unique or bifurcated for particular Re intervals. The latter case suggests the occurrence of switching from one type of flow structure organization to another with the ambient conditions unchanged.     

Transient response of enstrophy transport to opposition control in turbulent channel flow

The transient response of the turbulent enstrophy transport to opposition control in the turbulent channel flow is studied with the aid of direct numerical simulation. It is found that the streamwise enstrophy and the spanwise enstrophy are suppressed by the attenuation of the stretching terms at first, while the vertical enstrophy is reduced by inhibiting the tilt of the mean shear. In the initial period of the control, the streamwise enstrophy evolves much slower than the other two components. The vertical vorticity component exhibits a rapid monotonic decrease and also plays an important role in the attenuation of the other two components.     

Computation of Conditional Average Scalar Dissipation in Turbulent Jet Diffusion Flames

The modelling of conditional scalar dissipation in locally self-similar turbulent reacting jets is considered. The streamwise dependence in the transport equation of the conserved scalar pdf is represented by a function solely dependent on centreline mixture fraction. This procedure provides a simple model suitable for non-homogeneous flows and ensures positive values for conditional scalar dissipation. It has been tested in pure hydrogen-air jet diffusion flames using a Conditional Moment Closure method with detailed 12species, 23 reactions chemistry. The calculations show good agreement of the averaged scalar dissipation with reference values and the model proves to be superior to previous models based on homogeneous flows if the distribution of the conditional scalar dissipation in mixture fraction space is compared with experimental results. A dependence of NO predictions on the model of conditional scalar dissipation can be observed. This revised version was published online in July 2006 with corrections to the Cover Date.     

research on rotational dispersion coefficient of fiber in turbulent shear flow of fiber suspension

The rotational dispersion coefficient of the fiber in the turbulent shear flow of fiber suspension was studied theoretically. The function of correlation moment between the different fluctuating velocity gradients of the flow was built firstly. Then the expres- sion, dependent on the characteristic length, time, velocity and a dimensionless parameter related to the effect of wall, of rotational dispersion coefficient is derived. The derived expression of rotational dispersion coefficient can be employed to the inhomogeneous and non-isotropic turbulent flows. Furthermore it can be expanded to three-dimensional turbulent flows and serves the theoretical basis for solving the turbulent flow of fiber suspension.     

气体示踪法研究燃烧器横向湍流混合扩散特性

用气体示踪法对侧边风燃烧器弱旋流一次风和侧边直流二次风横向湍流混合扩散特性进行了详细研究，讨论了各截面混合物浓度分布规律及混合强度分布规律，分析对燃料着火的影响，同时利用湍流运动方程反算气体湍流扩散系数