Modelling of turbulent reacting flows in furnace devices

The spatial combustion of turbulent jets in furnace devices is modelled numerically basing on equations for multi-component turbulent reacting gaseous mixtures. The dependencies are obtained for the influence of the secondary air velocity and composition of gaseous components on torch configuration at a diffusion combustion process. The effect of regime parameters on the increase in torch sizes, which arises at the interaction of secondary air with gaseous components, has been elucidated.     

托卡马克边缘湍流中带状流的负粘滞阻尼

基于Hasegawa-Wakatani湍流模型,数值模拟了托卡马克边缘等离子体中漂移波湍流和相关的反常粒子输运。从等离子体动量守恒方程出发导出了不采用常规的布辛涅斯克近似的带状流方程,论证了大振幅密度扰动和湍性粒子流对激发带状流的贡献可等效地对应于低阶负粘滞阻尼效果。数值模拟表明,大振幅密度扰动的非线性大大增强了带状流饱和振幅,从而有效抑制了湍性粒子输运。研究结果阐明了托卡马克边缘等离子体大振幅密度扰动的非线性对驱动等离子体旋转、动量输运和带状流的重要性。     

�п���˱�Ե�����д�״���ĸ�ճ������

基于Hasegawa-Wakatani湍流模型,数值模拟了托卡马克边缘等离子体中漂移波湍流和相关的反常粒子输运.从等离子体动量守恒方程出发导出了不采用常规的布辛涅斯克近似的带状流方程,论证了大振幅密度扰动和湍性粒子流对激发带状流的贡献可等效地对应于低阶负粘滞阻尼效果.数值模拟表明,大振幅密度扰动的非线性大大增强了带状流饱和振幅,从而有效抑制了湍性粒子输运.研究结果阐明了托卡马克边缘等离子体大振幅密度扰动的非线性对驱动等离子体旋转、动量输运和带状流的重要性.     

Suitability of artificial bulk viscosity for large-eddy simulation of turbulent flows with shocks

The artificial bulk viscosity method to numerically capture shocks is investigated for large-eddy simulation (LES). Different variations of this method are tested on a turbulent flow over a cylinder at Reynolds number of 10,000 and free-stream Mach number of 0.85. The artificial bulk viscosity model by Cook and Cabot, which is parameterized by the strain rate magnitude, is found to provide unnecessary bulk viscosity in turbulent regions away from shocks. While developed turbulent structures are found unaffected, this extra bulk viscosity is shown to significantly damp the sound field. An alternative formulation of the model which is parameterized by the rate of dilatation is proposed. This formulation is shown to avoid the unnecessary bulk viscosity and enhance the sound-prediction capability of the model. It was found that standard LES combined with artificial bulk viscosity is a promising approach for simulation of turbulent flows with shocks. The formulation of the model on curvilinear coordinates is presented in the appendix.     

Large eddy simulation of orientation and rotation of ellipsoidal particles in isotropic turbulent flows

The rotational motion and orientational distribution of ellipsoidal particles in turbulent flows are of significance in environmental and engineering applications. Whereas the translational motion of an ellipsoidal particle is controlled by the turbulent motions at large scales, its rotational motion is determined by the fluid velocity gradient tensor at small scales, which raises a challenge when predicting the rotational dispersion of ellipsoidal particles using large eddy simulation (LES) method due to the lack of subgrid scale (SGS) fluid motions. We report the effects of the SGS fluid motions on the orientational and rotational statistics, such as the alignment between the long axis of ellipsoidal particles and the vorticity, the mean rotational energy at various aspect ratios against those obtained with direct numerical simulation (DNS) and filtered DNS. The performances of a stochastic differential equation (SDE) model for the SGS velocity gradient seen by the particles and the approximate deconvolution method (ADM) for LES are investigated. It is found that the missing SGS fluid motions in LES flow fields have significant effects on the rotational statistics of ellipsoidal particles. Alignment between the particles and the vorticity is weakened; and the rotational energy of the particles is reduced in LES. The SGS-SDE model leads to a large error in predicting the alignment between the particles and the vorticity and over-predicts the rotational energy of rod-like particles. The ADM significantly improves the rotational energy prediction of particles in LES.     

The evolution of velocity profiles and turbulent viscosity in a system of currents with wind-induced and density flows

The results of an in-situ study and mathematical modeling of a system of currents with wind-induced and density flows are presented. The model is based on the obtained dependences of its key parameters on the density-field characteristics, changes in the water level, and the water-reservoir topography. The revealed features of the of shear-velocity profile are indicated and a function of its distribution from the bottom to the open surface of the water is given. The model that was developed is verified by the data from measuring the parameters of currents and composition of water in the Petrozavodsk Bay of Onega Lake in September 2007. The regularities for the evolution of distributions of the coefficient of the turbulent exchange and current velocity over depth and in time are revealed.     

Modelling of turbulent boundary layers on the body of revolution in the presence of large eddy breakup devices

A novel approach was proposed to developing a combined algebraic-differential model of turbulent viscosity, which was then used to modify the traditional method for calculating turbulent boundary layers on bodies of revolution immersed in an incompressible flow. Ample experimental data were invoked to test the developed model through prediction of features inherent in turbulent boundary layers on the body of revolution in the presence of Large Eddy BreakUp (LEBU) devices installed in the vicinity of the streamlined surface. Within the developed approach, we showed it possible to calculate the evolution of local governing characteristics of turbulent flow in a broad range of LEBU parameters.     

Reduction of turbulent transport with zonal flows enhanced in helical systems

Gyrokinetic Vlasov simulations of the ion temperature gradient turbulence are performed in order to investigate effects of helical magnetic configurations on turbulent transport and zonal flows. The obtained results confirm the theoretical prediction that helical configurations optimized for reducing neoclassical ripple transport can simultaneously reduce the turbulent transport with enhancing zonal-flow generation. Stationary zonal-flow structures accompanied with transport reduction are clearly identified by the simulation for the neoclassically optimized helical geometry. The generation of the stationary zonal flow explains a physical mechanism for causing the confinement improvement observed in the inward-shifted plasma in the Large Helical Device [O. Motojima, Nucl. Fusion 43, 1674 (2003)].     

Prediction of wall-pressure fluctuation in turbulent flows with an immersed boundary method

The objective of this paper is to assess the accuracy and efficiency of the immersed boundary (IB) method to predict the wall pressure fluctuations in turbulent flows, where the flow dynamics in the near-wall region is fundamental to correctly predict the overall flow. The present approach achieves sufficient accuracy at the immersed boundary and overcomes deficiencies in previous IB methods by introducing additional constraints – a compatibility for the interpolated velocity boundary condition related to mass conservation and the formal decoupling of the pressure on this surfaces. The immersed boundary-approximated domain method (IB-ADM) developed in the present study satisfies these conditions with an inexpensive computational overhead. The IB-ADM correctly predicts the near-wall velocity, pressure and scalar fields in several example problems, including flows around a very thin solid object for which incorrect results were obtained with previous IB methods. In order to have sufficient near-wall mesh resolution for LES and DNS computations, the present approach uses local mesh refinement. The present method has been also successfully applied to computation of the wall-pressure space–time correlation in DNS of turbulent channel flow on grids not aligned with the boundaries. When applied to a turbulent flow around an airfoil, the computed flow statistics – the mean/RMS flow field and power spectra of the wall pressure – are in good agreement with experiment.     

Modelling of particle flows in a scrape-off layer with limiter biasing

The superconducting tokamak Tore Supra will be equipped with an actively cooled toroidal pump limiter (TPL), in the framework of the CIEL (Composants Internes Et Limiteurs) project, dedicated to plasma facing component design for steady state operation. The TPL is equipped with throats, located only on the high field side, for particle collection allowing the control of plasma density which is essential for long plasma discharges. The present design work of the CIEL includes a biasing system in order to enhance the particle pumping. A fluid model, based on the classical fluid equation, is used to estimate the effects of the electric field on the particle flows in the Scrape-Off Layer (SOL). The modifications of the density, the particle flow (toroidal and poloidal) and the position of the stagnation point are discussed as a function of the bias voltage. The model clearly illustrates the different resulting effects on particle pumping for a divertor and a limiter configuration which are designed respectively for poloidal or parallel particle collection. The model is used to interpret the ALT-biasing experiments recently carried out on TEXTOR-94. The pumping capability is shown to be improved by about (15–20)% for positive biasing while the experimental measurements of parallel Mach number are reproduced as a function of the applied voltage. The e-folding length of the edge density in the SOL is also shown to increase from 1.5 to about 2.0 cm for a biased voltage of −400 to 400 V, respectively, in accordance with the model. Finally, the model is used to extrapolate the TEXTOR-94 results to CIEL suggesting that pumping speed enhancement of 25 to 30% may be obtained. Partner in Trilateral Euregio Cluster Partner in Trilateral Euregio Cluster Presented at the Workshop on Role of Electric Fields in Plasma Confinement and Exhaust, Budapest, 18–19 June 2000.     

Instantaneous three-dimensional visualization of concentration distributions in turbulent flows with crossed-plane laser-induced fluorescence imaging

A laser-based technique for measuring instantaneous three-dimensional species concentration distributions in turbulent flows is presented. The laser beam from a single laser is formed into two crossed light sheets that illuminate the area of interest. The laser-induced fluorescence (LIF) signal emitted from excited species within both planes is detected with a single camera via a mirror arrangement. Image processing enables the reconstruction of the three-dimensional data set in close proximity to the cutting line of the two light sheets. Three-dimensional intensity gradients are computed and compared to the two-dimensional projections obtained from the two directly observed planes. Volume visualization by digital image processing gives unique insight into the three-dimensional structures within the turbulent processes. We apply this technique to measurements of toluene-LIF in a turbulent, non-reactive mixing process of toluene and air and to hydroxyl (OH) LIF in a turbulent methane-air flame upon excitation at 248 nm with a tunable KrF excimer laser. PACS 42.30.Va; 32.50.+d; 42.62.Cf     

The interaction of turbulent mixing region with local density perturbation in a pycnocline

Efficient numerical models were developed for the dynamics of local and turbulent formations in a pycnocline for which the averaged flow is characterized by appearance of internal waves of soliton type. Numerical analysis of the process of initiation, interaction, and subsequent propagation of internal waves generated by two local density perturbations in a pycnocline is carried out. Anisotropic turbulence degeneration in the turbulent mixing region in a pycnocline is numerically modeled. The problem of interaction of the turbulent mixing region with the local density perturbation in a pycnocline is numerically investigated in a wide range of local perturbation parameters.     

Coupling of mixing models with manifold based simplified chemistry in PDF modeling of turbulent reacting flows

For general reacting flows the numerical simulation faces two main challenges. One is the high dimensionality and stiffness of the governing conservation equations due to detailed chemistry, which can be solved by using simplified chemical kinetics. The other one is the difficulty of modeling the coupling of turbulence with thermo-chemical source term. The probability density function (PDF) method allows to calculate turbulent reacting flows by solving the thermal-chemical source term in closed form. Usually, the PDF method for turbulent processes such as mixing processes and the reduction method for chemical kinetics are developed separately. However, coupling of both processes plays an important role for the numerical accuracy. To investigate the importance of coupling between turbulence and simplified chemistry, two different coupling strategies for mixing and reduced chemistry are discussed and tested for the well-known Sandia Flames E and F, in which there is a strong interaction between turbulence and chemical kinetics. The EMST mixing model is chosen for turbulent mixing, while the Reaction-Diffusion Manifolds (REDIMs) is used as simplified chemistry. However, the proposed strategies are also valid for other mixing models and manifold based simplified chemistry.     

Modelling of thermogravitation convection in closed volume with local sources of heat release

A mathematical model is presented for unsteady conjugate convective-conductive heat transfer in a closed volume with local sources of heat release under the conditions of a convective-radiation heat exchange on one of external faces of the solution region. The thermogravitation convection regime was analysed for moderate Grashof numbers. The typical temperature and velocity fields were obtained, and the fields of desired quantities were compared for the planar and three-dimensional models in one of the typical sections of the solution region. The work was supported by the Russian Foundation for Basic Research and Administration of the Tomsk Region (No. 05-02-98006, competition r_ob_a).     

Mathematical simulation of conjugate turbulent natural convection in an enclosure with local heat source

A numerical analysis of turbulent regimes of the natural convection in a closed rectangular region with heat-conducting walls of finite thickness was carried out in the presence of a locally concentrated heat source under the conditions of the radiative-convective heat exchange with the ambient medium on one of the external boundaries. The mathematical model was constructed on the basis of the Reynolds equations in dimensionless variables stream function — vorticity vector — temperature. Special attention was paid to the investigation of the influence of the Grashof number er 10⁸≤Gr<10¹⁰, of the unsteadiness factor 0< τ <1000, and the thermal conductivity ratio λ _2,1 = 5.7·10⁻⁴, 6.8·10⁻⁵ on both the local and integral characteristics of the problem.