A stability condition for turbulence model： From EMMS model to EMMS-based turbulence model

The closure problem of turbulence is still a challenging issue in turbulence modeling. In this work, a stability condition is used to close turbulence. Specifically, we regard single-phase flow as a mixture of turbulent and non-turbulent fluids, separating the structure of turbulence. Subsequently, according to the picture of the turbulent eddy cascade, the energy contained in turbulent flow is decomposed into different parts and then quantified. A turbulence stability condition, similar to the principle of the energy-minimization multi-scale （EMMS） model for gas-solid systems, is formulated to close the dynamic constraint equa- tions of turbulence, allowing the inhomogeneous structural parameters of turbulence to be optimized. We name this model as the ＂EMMS-based turbulence model＂, and use it to construct the corresponding turbulent viscosity coefficient. To validate the EMMS-based turbulence model, it is used to simulate two classical benchmark problems, lid-driven cavity flow and turbulent flow with forced convection in an empty room, The numerical results show that the EMMS-hased turbulence model improves the accuracy of turbulence modeling due to it considers the principle of compromise in competition between viscosity and inertia.     

A stability condition for turbulence model: From EMMS model to EMMS-based turbulence model

The closure problem of turbulence is still a challenging issue in turbulence modeling. In this work, a stability condition is used to close turbulence. Specifically, we regard single-phase flow as a mixture of turbulent and non-turbulent fluids, separating the structure of turbulence. Subsequently, according to the picture of the turbulent eddy cascade, the energy contained in turbulent flow is decomposed into different parts and then quantified. A turbulence stability condition, similar to the principle of the energy-minimization multi-scale (EMMS) model for gas–solid systems, is formulated to close the dynamic constraint equations of turbulence, allowing the inhomogeneous structural parameters of turbulence to be optimized. We name this model as the “EMMS-based turbulence model”, and use it to construct the corresponding turbulent viscosity coefficient. To validate the EMMS-based turbulence model, it is used to simulate two classical benchmark problems, lid-driven cavity flow and turbulent flow with forced convection in an empty room. The numerical results show that the EMMS-based turbulence model improves the accuracy of turbulence modeling due to it considers the principle of compromise in competition between viscosity and inertia.     

充分发展圆管湍流的实验研究

采用粒子数字图像测速（digital particle image velocimetry,DPIV)和定量流动显示技术（quantitative flow visualization,QFA)对充分发展的圆管湍流进行了研究。测量结果和直接数值模拟（direct numerical simulation,DNS)结果进行了比较，结果表明作者开发的DPIV技术取得了满意的精度。在此基础上对圆管湍流的动力学机理进行了研究，分析了上抛和下扫在湍流生成中的贡献以及流动显示结构内的脉动速度分布，测量结果显示在圆管湍流的近壁区存在横向强脉冲现象和流动显示所能观察到的结构为上抛占主导地位的结构。     

Modeling the hydrodynamics of cocurrent gas–solid downers according to energy-minimization multi-scale theory

Cocurrent gas–solid downer reactors have many applications in industry because they possess the technological advantages of a lower pressure drop, shorter residence time, and less solid backmixing when compared with traditional circulating fluidized bed risers. By introducing the concept of particle clusters explicitly, a one-dimensional model with consideration of the interphase interactions between the fluid and particles at both microscale and mesoscale is formulated for concurrent downward gas–solid flow according to energy-minimization multi-scale (EMMS) theory. A unified stability condition is proposed for the differently developed sections of gas–solid flow according to the principle of the compromise in competition between dominant mechanisms. By optimizing the number density of particle clusters with respect to the stability condition, the formulated model can be numerically solved without introducing cluster-specific empirical correlations. The EMMS-based model predicts well the axial hydrodynamics of cocurrent gas–solid downers and is expected to have a wider range of applications than the existing cluster-based models.     

旋转圆管湍流的大涡模拟数值研究

采用大涡模拟（LES）方法,并结合动力学亚格子尺度应力（SGS）模型,通过数值求解柱坐标系下的滤波Navier-Stokes方程,研究了绕管轴旋转圆管内的湍流流动特性.为验证计算的可靠性,以及动力学SGS模型对于旋转湍流的适用性,将大涡模拟计算所得的结果,与相应的直接模拟（DNS）结果和实验数据进行了对比验证,吻合良好.进一步对旋转圆管湍流的物理机理进行了探讨,研究了湍流特性随旋转速率的变化规律.当旋转速率增加时,湍流流动有层流化的发展趋势.基于湍动能变化的关系,分析了旋转效应对湍流脉动生成的抑制作用.     

Large-Scale Turbulent Vortical Structures inside a Sudden Expansion Cylinder Chamber

A large eddy simulation (LES) is performed for turbulent flow around a bluff body inside a sudden expansion cylinder chamber, a configuration which resembles a premixed gas turbine combustor. To promote turbulent mixing and to accommodate flame stability, a flame holder is installed inside the combustion chamber. The Smagorinsky model and the Lagrangian dynamic subgrid-scale model are employed and tested. The calculated Reynolds number is 5,000 based on the bulk velocity and the diameter of inlet pipe. The simulation code is constructed by using a general coordinate system based on the physical contravariant velocity components. The predicted turbulent statistics are evaluated by comparing with the laser-doppler velocimetry (LDV) measurement data. The agreement of LES with the experimental data is shown to be satisfactory. Emphasis is placed on the time-dependent evolutions of turbulent vortical structures behind the flame holder. The numerical flow visualizations depict the behavior of large-scale vortices. The turbulent behavior behind the flame holder is analyzed by visualizing the sectional views of vortical structure. This revised version was published online in July 2006 with corrections to the Cover Date.     

Shallow cavity flow tone experiments: onset of locked-on states

Fully turbulent inflow past a shallow cavity is investigated for the configuration of an axisymmetric cavity mounted in a pipe. Emphasis is on conditions giving rise to coherent oscillations, which can lead to locked-on states of flow tones in the pipe–cavity system. Unsteady surface pressure measurements are interpreted using three-dimensional representations of amplitude–frequency, and velocity; these representations are constructed for a range of cavity depth. Assessment of these data involves a variety of approaches. Evaluation of pressure gradients on plan views of the three-dimensional representations allows extraction of the frequencies of the instability (Strouhal) modes of the cavity oscillation. These frequency components are correlated with traditional models originally formulated for cavities in a free-stream. In addition, they are normalized using two length scales: inflow boundary-layer thickness and pipe diameter. These scales are consistent with those employed for the hydrodynamic instability of the separated shear layer, and are linked to the large-scale mode of the shear layer oscillation, which occurs at relatively long cavity length. In fact, a simple scaling based on pipe diameter can correlate the frequencies of the dominant peaks over a range of cavity depth.The foregoing considerations provide evidence that pronounced flow tones can be generated from a fully turbulent inflow at very low Mach number, including the limiting case of fully developed turbulent flow in a pipe. These tones can arise even for the extreme case of a cavity having a length over an order of magnitude longer than its depth. Suppression of tones is generally achieved if the cavity is sufficiently shallow.     

Experimental investigation of the coherent structures in a spirally corrugated pipe

Previous numerical and theoretical results (Chen et al., 2019; Liu et al., 2018; Zhao et al., 2019) based on the optimization theory of convective heat transfer reveal that the optimized flow structures in a straight circular pipe enhancing convective heat transfer are multiple longitudinal vortices. This conclusion encourages us to find out whether such flow structures really exist in some enhanced heat transfer pipes by means of advanced experimental techniques. Therefore, a typical enhanced heat transfer pipe was selected, namely a spirally corrugated pipe, and stereoscopic particle image velocimetry (SPIV) was employed to measure its internal instantaneous flow field. Moreover, the proper orthogonal decomposition (POD) method was used to extract the large-scale coherent structures from the measured instantaneous velocity fields. Besides the spirally corrugated pipe, the fully developed turbulent flow in a straight pipe was also analyzed as benchmark of the enhanced heat transfer pipes. The results reveal that longitudinal whirling flow with multi-vortices is formed in both the fully developed turbulent flow field of the straight pipe and the spirally corrugated one. It is thus easy to explain the heat transfer enhancement mechanism of the above flow structures from the perspective of momentum transfer. The flow structures of the fully developed turbulent flow in a straight pipe are quite similar to the optimal flow pattern from the optimization theory. More specifically, multiple longitudinal vortices are spontaneously generated due to turbulence without external heat transfer enhancement techniques. Furthermore, the flow structures similar to multiple longitudinal vortices also exist in the spirally corrugated pipe, although these flow structures deviate from symmetric multiple vortices. Moreover, the flow structures in the spirally corrugated pipe are much more energetic than those in the fully developed turbulent flow in a straight pipe. This is probably the reason why a spirally corrugated pipe can enhance heat transfer compared with a straight circular pipe.     

A numerical study of bingham turbulent flow in sudden-expansion straight circular pipe

I.IntroductionBinghamfluidisonebranchofnon-Newtonianfluid,suchascrudeparaffinoil,highsediment--ladenwaterflow,highconcentrationmudandthelikewhicharetransportedinpipelinesinmanyindustries,soit'sofgreatsignificancetostudytheflowmechanismsofBinghamfluid.Tsaietal.II]studiedthelinkagebetweenBinghamfluidandpluggedflow.Wangetal.I2]measuredtheturbulencestructureofBinghammud.Mengetal.[3]researchedthekineticenergycorrectionfactorofBinghamfluidinacircularpipe.However,thestudyofBinghamfluidsofarisn't…     

A dynamic global-coefficient mixed subgrid-scale model for large-eddy simulation of turbulent flows

A dynamic global-coefficient mixed subgrid-scale eddy-viscosity model for large-eddy simulation of turbulent flows in complex geometries is developed. In the present model, the subgrid-scale stress is decomposed into the modified Leonard stress, cross stress, and subgrid-scale Reynolds stress. The modified Leonard stress is explicitly computed assuming a scale similarity, while the cross stress and the subgrid-scale Reynolds stress are modeled using the global-coefficient eddy-viscosity model. The model coefficient is determined by a dynamic procedure based on the global-equilibrium between the subgrid-scale dissipation and the viscous dissipation. The new model relieves some of the difficulties associated with an eddy-viscosity closure, such as the nonalignment of the principal axes of the subgrid-scale stress tensor and the strain rate tensor and the anisotropy of turbulent flow fields, while, like other dynamic global-coefficient models, it does not require averaging or clipping of the model coefficient for numerical stabilization. The combination of the global-coefficient eddy-viscosity model and a scale-similarity model is demonstrated to produce improved predictions in a number of turbulent flow simulations.     

Anisotropic Turbulence and Coherent Structures in Eccentric Annular Channels

Eccentric annular pipe flows represent an ideal model for investigating inhomogeneous turbulent shear flows, where conditions of turbulence production and transport vary significantly within the cross-section. Moreover recent works have proven that in geometries characterized by the presence of a narrow gap, large-scale coherent structures are present. The eccentric annular channel represents, in the opinion of the present authors, the prototype of these geometries. The aim of the present work is to verify the capability of a numerical methodology to fully reproduce the main features of the flow field in this geometry, to verify and characterize the presence of large-scale coherent structures, to examine their behavior at different Reynolds numbers and eccentricities and to analyze the anisotropy associated to these structures. The numerical approach is based upon LES, boundary fitted coordinates and a fractional step algorithm. A dynamic Sub Grid Scale (SGS) model suited for this numerical environment has been implemented and tested. An additional interest of this work is therefore in the approach employed itself, considering it as a step into the development of an effective LES methodology for flows in complex channel geometries. Agreement with previous experimental and DNS results has been found good overall for the streamwise velocity, shear stress and the rms of the velocity components. The instantaneous flow field presented large-scale coherent structures in the streamwise direction at low Reynolds numbers, while these are absent or less dominant at higher Reynolds and low eccentricity. After Reynolds averaging is performed over a long integration time the existence of secondary flows in the cross session is proven. Their shape is found to be constant over the Reynolds range surveyed, and dependent on the geometric parameters. The effect of secondary flows on anisotropy is studied over an extensive Reynolds range through invariant analysis. Additional insight on the mechanics of turbulence in this geometry is obtained.     

An Investigation of SGS Stress Anisotropy Modeling in Complex Turbulent Flow Fields

An anisotropy-resolving subgrid-scale (SGS) model for large eddy simulation was investigated. Primary attention was given to the predictive performance of the SGS model in the case of complex turbulence with flow impingement and/or flow separation. The SGS model was constructed by combining an isotropic linear eddy-viscosity model with an extra anisotropic term. Since the extra anisotropic term was modeled to prevent undesirable energy transfer between the grid-scale and SGS parts, the model is expected not to seriously affect computational stability. To validate the model performance for complex turbulent flow fields, the SGS model was applied to numerical simulations of a plane impinging jet and 3-D diffuser flow as well as fundamental plane channel flows. The SGS model provided reasonable predictions for these test cases. Furthermore, the predicted SGS stress components were decomposed into linear and anisotropic parts and their roles were investigated in detail. The usefulness of the present anisotropy-resolving SGS model in practical engineering applications was thus described.     

Subgrid-Scale Stress Modelling Based on the Square of the Velocity Gradient Tensor

A new subgrid scale model is proposed for Large Eddy Simulations in complex geometries. This model which is based on the square of the velocity gradient tensor accounts for the effects of both the strain and the rotation rate of the smallest resolved turbulent fluctuations. Moreover it recovers the proper y ³ near-wall scaling for the eddy viscosity without requiring dynamic procedure. It is also shown from a periodic turbulent pipe flow computation that the model can handle transition.     

Frequency response correction of hot-film measurements in a turbulent airflow

Turbulence measurements were performed in a fully developed turbulent pipe flow using a wedged-shaped hot-film probe. The dynamic sensitivity of hot-film probes is frequency dependent and cannot be obtained from static calibration. As it is complicated to make a dynamic calibration, a correction method for the turbulent energy spectrum is used instead. The measurements were compared with hot-wire turbulence measurements and show good agreement. Received: 28 May 1999/Accepted: 24 August 1999     

Numerical prediction of turbulent flows using Reynolds‐averaged Navier–Stokes and large‐eddy simulation with uncertain inflow conditions

Numerous comparisons between Reynolds‐averaged Navier–Stokes (RANS) and large‐eddy simulation (LES) modeling have already been performed for a large variety of turbulent flows in the context of fully deterministic flows, that is, with fixed flow and model parameters. More recently, RANS and LES have been separately assessed in conjunction with stochastic flow and/or model parameters. The present paper performs a comparison of the RANS k ? ε model and the LES dynamic Smagorinsky model for turbulent flow in a pipe geometry subject to uncertain inflow conditions. The influence of the experimental uncertainties on the computed flow is analyzed using a non‐intrusive polynomial chaos approach for two flow configurations (with or without swirl). Measured quantities including an estimation of the measurement error are then compared with the statistical representation (mean value and variance) of their RANS and LES numerical approximations in order to check whether experiment/simulation discrepancies can be explained within the uncertainty inherent to the studied configuration. The statistics of the RANS prediction are found in poor agreement with experimental results when the flow is characterized by a strong swirl, whereas the computationally more expensive LES prediction remains statistically well inside the measurement intervals for the key flow quantities.Copyright © 2012 John Wiley & Sons, Ltd.     

Nonlinear dynamics of cantilevered pipes conveying fluid: Towards a further understanding of the effect of loose constraints

The nonlinear dynamics of a fluid-conveying cantilevered pipe with loose constraints placed somewhere along its length is investigated. The main objective of this study is to determine the effects of several geometrical and physical parameters of the loose constraints on the characteristics and behavior of pipes conveying fluid. Based on the full nonlinear equation of motion, the dynamical behavior of the pipe system is investigated. Phase portraits and bifurcation diagrams are constructed for a selected set of system parameters. Typical results are firstly compared to numerical ones reported previously and excellent agreement is obtained. Then, the threshold flow velocities for several key bifurcations including pitchfork, period doubling, chaos, and sticking behaviors are predicted, showing that in many cases, the gap size, stiffness, and asymmetry of the loose constraints have remarkable effects on the nonlinear responses of the cantilevered pipe conveying fluid. For a pipe system with small/large constraint gap sizes, small constraint stiffness, or large constraint offset, some of the complex dynamical behaviors including chaos and period-doubling bifurcations would disappear, at least in the flow velocity range of interest.     

Turbulent flow in a conical diffuser: A review

This is a review of experimental studies of turbulent flow in a conical diffuser by eight Ph.D. students, eleven M.Sc. students, one M.Eng. student, and myself in the past 29 years. During this time, two conical diffusers were constructed: the first was of cast aluminum construction, and the second was of plastic fabrication. These two diffusers were basically the same in geometry except that the pipe section was constructed as an integral part of the plastic diffuser to avoid the lip at the junction of the inlet pipe and the diffuser. The conical diffuser had a total divergence angle of 8°, an area ratio of 4:1, and an inlet diameter of 0.1016 m (4 in.).

The flow at the inlet of the diffuser was usually fully developed pipe flow, but sometimes it was boundary layer grown on the pipe wall. Hot-wire and pulse-wire anemometry together with computer facilities were used to obtain the results of complex flow present in the conical diffuser. Mean velocity profiles were obtained throughout the diffuser, which in turn were used to obtain strain rates and their principal direction. Turbulence moments up to fourth order were measured. The results were used to assess momentum, turbulent kinetic energy, and shear stress equations. Other features such as instantaneous flow reversals in the wall region, relative strength of large eddies, extra strain rate, and the production of kinetic energy also were investigated to find the dynamical picture in the diffuser flow.     

多体系统接触碰撞问题的牛顿-欧拉线性互补方法

基于非光滑动力学方法的多体系统接触碰撞分析是目前多体系统动力学的研究热点.本文采用牛顿-欧拉方法建立多体系统接触、碰撞问题的动力学模型,给出一种牛顿-欧拉型线性互补公式.该建模方法与目前一般采用的拉格朗日建模方法的不同之处是约束条件中除了库仑摩擦、单边约束之外还含有光滑等式约束.在建立系统动力学模型时,首先解除摩擦约束和单边约束得到原系统对应的基本系统.牛顿-欧拉方法采用最大数目坐标建立基本系统的动力学方程,由于坐标不相互独立,因此基本系统中带有等式约束,其数学模型为一组微分代数方程.借助约束雅可比矩阵,在基本系统微分代数方程中添加摩擦接触和单边约束对应的拉氏乘子,就可以得到系统全局运动的具有变拓扑结构特征的动力学方程,再结合非光滑约束互补条件便可构成完备的系统动力学模型.完备的动力学模型由动力学微分方程以及等式约束和不等式约束组成.线性互补公式采用分块矩阵形式进行推导,简化了推导过程.数值计算采用基于线性互补的时间步进算法.时间步进算法是目前流行的非光滑数值算法,其突出特点是可以免去数值积分中繁琐的事件检测过程,而数值积分过程中通过对线性互补问题的求解可以确定系统的触-离状态.通过对典型的曲柄滑块间隙机构进行数值分析,验证本文方法的有效性.     

Analysis of pipe flow transition. Part II. Energy transfer

The present article describes the results from a study of nonlinear mechanisms at work during the process of transition to turbulence in pipe flows. Using an accurate hybrid finite-difference code for the simulation of unsteady incompressible pipe flow, we have performed a direct numerical simulation designed to model experiments performed by Han, Tumin and Wygnanski [12]. Based on these numerical data, we have conducted a meticulous investigation of the dynamic interactions of the structures and flow modes that can be observed during this process. Based on this study, we can paint a detailed picture of the dynamical interactions of flow structures during both the linear and nonlinear stages of pipe flow transition. While this picture does have some similarities to earlier proposed mechanisms, we find that even for the simple cases considered here the structure of the pertinent interactions is much richer than suggested by these earlier models.     

Transport modeling of sedimenting particles in a turbulent pipe flow using Euler–Lagrange large eddy simulation

A volume-filtered Euler–Lagrange large eddy simulation methodology is used to predict the physics of turbulent liquid–solid slurry flow through a horizontal periodic pipe. A dynamic Smagorinsky model based on Lagrangian averaging is employed to account for the sub-filter scale effects in the liquid phase. A fully conservative immersed boundary method is used to account for the pipe geometry on a uniform cartesian grid. The liquid and solid phases are coupled through volume fraction and momentum exchange terms. Particle–particle and particle–wall collisions are modeled using a soft-sphere approach. Three simulations are performed by varying the superficial liquid velocity to be consistent with the experimental data by Dahl et al. (2003). Depending on the liquid flow rate, a particle bed can form and develop different patterns, which are discussed in light of regime diagrams proposed in the literature. The fluctuation in the height of the liquid-bed interface is characterized to understand the space and time evolution of these patterns. Statistics of engineering interest such as mean velocity, mean concentration, and mean streamwise pressure gradient driving the flow are extracted from the numerical simulations and presented. Sand hold-up calculated from the simulation results suggest that this computational strategy is capable of predicting critical deposition velocity.