Flow Dynamics of a Spiral-groove Dry-gas Seal

The dry-gas seal has been widely used in different industries. With increased spin speed of the rotator shaft, turbulence occurs in the gas film between the stator and rotor seal faces. For the micro-scale flow in the gas film and grooves, turbulence can change the pressure distribution of the gas film. Hence, the seal performance is influenced. However, turbulence effects and methods for their evaluation are not considered in the existing industrial designs of dry-gas seal. The present paper numerically obtains the turbulent flow fields of a spiral-groove dry-gas seal to analyze turbulence effects on seal performance. The direct numerical simulation （DNS） and Reynolds-averaged Navier-Stokes （RANS） methods are utilized to predict the velocity field properties in the grooves and gas film. The key performance parameter, open force, is obtained by integrating the pressure distribution, and the obtained result is in good agreement with the experimental data of other researchers. Very large velocity gradients are found in the sealing gas film because of the geometrical effects of the grooves. Considering turbulence effects, the calculation results show that both the gas film pressure and open force decrease. The RANS method underestimates the performance, compared with the DNS. The solution of the conventional Reynolds lubrication equation without turbulence effects suffers from significant calculation errors and a small application scope. The present study helps elucidate the physical mechanism of the hydrodynamic effects of grooves for improving and optimizing the industrial design or seal face pattern of a dry-gas seal.

Boundary-layer disturbances subjected to free-stream turbulence and simulation on bypass transition

The phenomena associated with the entrainment of free-stream turbulence （FST） into boundary-layer flows are relevant for a number of subjects. It has been be- lieved that the continuous spectra of the Orr-Sommerfeld （O-S）/Squire equations describe the entrainment process, and thus they are used to specify the inlet condition in simulation of bypass transition. However, Dong and Wu （Dong, M. and Wu, X. On continuous spectra of the Orr-Sommerfeld/Squire equations and entrainment of free-stream vortical disturbances. Journal of Fluid Mechanics, 732, 616-659 （2013）） pointed out that continuous spectra exhibit several non-physical features due to neglecting the non-parallelism. They further proposed a large-Reynolds-number asymptotic approach, and showed that the non-parallelism is a leading-order effect even for the short-wavelength disturbance, for which the response concentrates in the edge layer. In this paper, the asymptotic solution is verified numerically by studying its evolution in incompressible boundary layers. It is found that the numerical results can be accurately predicted by the asymptotic solution, implying that the latter is adequate for moderate Reynolds numbers. By introducing a series of such solutions as the inflow perturbations, the bypass transition is investigated via the direct numerical simulation （DNS）. The transition processes, including the evolution of streaks, the amplification of secondary-instability modes, and the emergence of turbulent spots, agree with the experimental observations.

Numerical simulation on evolution of subharmonic low-speed streaks in minimal channel turbulent flow

The evolution of low-speed streaks in the turbulent boundary layer of the minimum channel flow unit at a low Reynolds number is simulated by the direct numer- ical simulation （DNS） based on the standard Fourier-Chebyshev spectral method. The subharmonic sinuous （SS） mode for two spanwise-aligned low-speed streaks is excited by imposing the initial perturbations. The possibilities and the physical realities of the turbulent sustaining in the minimal channel unit are examined. Based on such a flow field environment, the evolution of the low-speed streaks during a cycle of turbulent sus- taining, including lift-up, oscillation, and breakdown, is investigated. The development of streamwise vortices and the dynamics of vortex structures are examined. The results show that the vortices generated from the same streak are staggered along the streamwise direction, while the vortices induced by different streaks tilt toward the normal direction due to the mutual induction effect. It is the spatial variations of the streamwise vortices that cause the lift-up of the streaks. By resolving the transport dynamics of enstrophy, the strength of the vortices is found to continuously grow in the logarithmic layer through the vortex stretching mechanism during the evolution of streaks. The enhancement of the vortices contributes to the spanwise oscillation and the following breakdown of the low-speed streaks.

Self-consistent parabolized stability equation(PSE) method for compressible boundary layer

The parabolized stability equation （PSE） method has been proven to be a useful and convenient tool for the investigation of the stability and transition problems of boundary layers. However, in its original formulation, for nonlinear problems, the complex wave number of each Fourier mode is determined by the so-called phase-locked rule, which results in non-self-consistency in the wave numbers. In this paper, a modification is proposed to make it self-consistent. The main idea is that, instead of allowing wave numbers to be complex, all wave numbers are kept real, and the growth or decay of each mode is simply manifested in the growth or decay of the modulus of its shape function. The validity of the new formulation is illustrated by comparing the results with those from the corresponding direct numerical simulation （DNS） as applied to a problem of compressible boundary layer with Mach number 6.

Effect of particles on turbulent thermal field of channel flow with different Prandtl numbers

The direct numerical simulation （DNS） of heat transfer in a fully developed non-isothermal particle-laden turbulent channel flow is performed. The focus of this paper is on the modulation of the particles on turbulent thermal statistics in the particle-laden flow with three Prandtl numbers （Pτ = 0.71, 1.5, and 3.0） and a shear Reynolds number （Reτ = 180）. Some typical thermal statistics, including normalized mean temperature and their fluctuations, turbulent heat fluxes, Nusselt number and so on, are analyzed. The results show that the particles have less effects on turbulent thermal fields with the increase of Prandtl number. Two reasons can explain this. First, the correlation between fluid thermal field and velocity field decreases as the Prandtl number increases, and the modulation of turbulent velocity field induced by the particles has less influence on the turbulent thermal field. Second, the heat exchange between turbulence and particles decreases for the particle-laden flow with the larger Prandtl number, and the thermal feedback of the particles to turbulence becomes weak.

Two-phase micro-and macro-time scales in particle-laden turbulent channel flows

The micro- and macro-time scales in two-phase turbulent channel flows are investigated using the direct nu- merical simulation and the Lagrangian particle trajectory methods for the fluid- and the particle-phases, respectively. Lagrangian and Eulerian time scales of both phases are cal- culated using velocity correlation functions. Due to flow anisotropy, micro-time scales are not the same with the theo- retical estimations in large Reynolds number （isotropic） tur- bulence. Lagrangian macro-time scales of particle-phase and of fluid-phase seen by particles are both dependent on particle Stokes number. The fluid-phase Lagrangian inte- gral time scales increase with distance from the wall, longer than those time scales seen by particles. The Eulerian inte- gral macro-time scales increase in near-wall regions but de- crease in out-layer regions. The moving Eulerian time scales are also investigated and compared with Lagrangian integral time scales, and in good agreement with previous measure- ments and numerical predictions. For the fluid particles the micro Eulerian time scales are longer than the Lagrangian ones in the near wall regions, while away from the walls the micro Lagrangian time scales are longer. The Lagrangian integral time scales are longer than the Eulerian ones. The results are useful for further understanding two-phase flow physics and especially for constructing accurate prediction models of inertial particle dispersion.

Near-wall behaviors of oblique-shock-wave/turbulent-boundary-layer interactions

A direct numerical simulation （DNS） on an oblique shock wave with an incident angle of 33.2° impinging on a Mach 2.25 supersonic turbulent boundary layer is performed. The numerical results are confirmed to be of high accuracy by comparison with the reference data. Particular efforts have been made on the investigation of the near-wall behaviors in the interaction region, where the pressure gradient is so significant that a certain separation zone emerges. It is found that, the traditional linear and loga- rithmic laws, which describe the mean-velocity profiles in the viscous and meso sublayers, respectively, cease to be valid in the neighborhood of the interaction region, and two new laws of the wall are proposed by elevating the pressure gradient to the leading order. The new laws are inspired by the analysis on the incompressible separation flows, while the compressibility is additionally taken into account. It is verified by the DNS results that the new laws are adequate to reproduce the mean-velocity profiles both inside and outside the interaction region. Moreover, the normalization adopted in the new laws is able to regularize the Reynolds stress into an almost universal distribution even with a salient adverse pressure gradient （APG）.

Rotation invariant constitutive relation for Reynolds stress structure parameter

A new Reynolds stress constitutive formula is constructed using the first- order statistics of turbulent fluctuations instead of the mean strain rate. It includes zero empirical coefficients. The formula is validated with the direct numerical simulation （DNS） data of turbulent channel flow at Reτ=180. The Reynolds stresses given by the proposed formula agree very well with the DNS results. The good agreement persists even after the multi-angle rotation of the coordinate system, indicating the rotation in- variance of the formula. The autocorrelation of the fluctuating velocity rather than the mean strain rate is close to the essence of the Reynolds stress.

Numerical investigation of turbulent channel flow controlled by spatially oscillating spanwise Lorentz force

A formulation of the skin-friction drag related to the Reynolds shear stress in a turbulent channel flow is derived. A direct numerical simulation （DNS） of the turbulent control is performed by imposing the spatially oscillating spanwise Lorentz force. Under the action of the Lorentz force with several proper control parameters, only the periodi- cally well-organized streamwise vortices are finally observed in the near-wall region. The Reynolds shear stress decreases dramatically, especially in the near-wall area, resulting in a drag reduction.

Lagrangian time scales and its relationship to Eulerian equivalents in turbulent channel flow

Lagrangian and Eulerian time scales were obtained from the direct numerical simulation of turbulent channel flow at two Reynolds numbers based on the friction velocity and channel half-height, Rer= 80, 100. The Lagrangian integral time scales and time microscales were compared to their Eulerian equivalents. It is found that the ratio of Lagrangian to TL Eulerian integral time scales is given by TE/TiE= 1 ＋ 0.1y＋ for y＋ ≤ 10, and that the ratios between the Lagrangian to theEulerian time microscales are almost the same irrespective of the components. Those increase with y＋ are approximated by ≈ 2.75 - 1.75 exp （-v＋/a） . These results also show that these expressions are independent of the Reynolds number.

轴对称射流火焰的数值模拟及稳定性分析

燃烧过程是能源利用过程中不可缺少的环节,而燃烧不稳定现象则是燃烧过程优化应用中潜在的威胁,其导致的热释放率和压力波动会对燃烧器造成极大危害.燃烧不稳定性的产生机理研究和有效控制具有十分重要的意义。近年来很多学者围绕燃烧不稳定性做了大量的理论和实验研究,但绝大部分研究针对预混火焰而不是扩散火焰.人们认为热声的耦合作用是燃烧不稳定性产生的根本原因,瑞利(Rayleigh)准则提供了判断燃烧不稳定性的一个重要判据。本文通过轴对称射流火焰的直接数值模拟(DNS)探究了压力扰动和速度扰动对燃烧不稳定性的影响,并且利用Rayleigh准则对结果进行了分析。结果表明,通过适当地改变射流入口压力和速度脉动,能够使声能的增长受到抑制,从而为控制燃烧不稳定性提供有效的途径。

DIRECT NUMERICAL SIMULATION OF TURBULENT HEAT TRANSFER IN A WALL-NORMAL ROTATING CHANNEL FLOW

Direct Nmerical Simulation （DNS） of turbulent heat transfer in a wall-normal rotating channel flow has been carried out for the rotation number Nr from 0 to 0.1, the Reynolds number 194 based on the friction velocity of non ro taring case and the half-height of the channel, and the Prandtl number 1. The objective of this study is to reveal the effects of rotation on the characteristics of turbulent flow and heat transfer. Based on the present calculated results, two typical rotation regimes are identified. When 0 〈 Nr 〈 0.06, turbu lence and thermal statistics correlated with the spanwise veloc ity fluctuation are enhanced since the shear rate of spanwise mean flow induced by Coriolis force increases; however, the other statistics are suppressed. When Nr 〉 0.06, turbulence and thermal statistics are suppressed significantly because the Coriolis force effect plays as a dominated role in the rotating flow. Remarkable change of the direction of near wall streak structures based on the velocity and temperature fluctuations is identified.

TURBULENT OPEN CHANNEL FLOW SUBJECTED TO THE CONTROL OF A SPANWISE TRAVELING WAVE

Turbulent open channel flows subjected to the control of a spanwise traveling wave have been investigated by means of Direct Numerical Simulation （DNS）. The objective of this study is to reveal the response of the near-wall and surface-influenced turbulence to the spanwise traveling wave control. Three typical frequencies of the spanwise traveling wave, i.e., high-, middle- and low-frequency, corresponding to the exciting periods at 25, 50 and 100, are considered to study the turbulence dynamics in the wall and surface regions. To elucidate the behaviors of turbulence statistics, some typical quantities, including the mean velocity, velocity fluctuations and the structures of turbulence fluctuations, are exhibited and analyzed.

Crystal Structures, Thermal Behavior Analysis and Thermal Safety of Two Energetic Salts of Hydrazinyl-1,2,4,5-tetrazine with 3,5-Dinitrosalicylic Acid

Two energetic salts, DPHT.DNS-H20（1） and DHT.2DNS.2H20（2）[DPHT=3-（3,5-dimethyl-lH-pyrazol- 1-yl）-6-hydrazinyl-1,2,4,5-tetrazine; DHT=3,6-dihydrazinyl-l,2,4,5-tetrazine], were synthesized from S-tetrazine with 3,5-dinitrosalicylic acid（DNS）. Compounds 1 and 2 were structurally characterized by elemental analysis, infrared spectroscopy, and single-crystal X-ray diffraction. The thermal behavior of the title compounds was studied by differential scanning calorimetry（DSC） and thermogravimetry（TG）. The non-isothermal decomposition kinetics of compound 2 were investigated. The self-accelerating decomposition temperature, thermal ignition temperature, and critical temperatures of thermal explosion were obtained to evaluate the thermal safety of compound 2. The results show compounds 1 and 2 decompose at 150.8 and 179.2℃, respectively. The TSADT and Tb of compound 2 are higher than those of DHT, which indicates compound 2 is a potential candidate for energetic materials that have good thermal stability. Keywords Tetrazine compound; Dinitrosalicylic acid（DNS）; Crystal structure; Thermal behavior; Thermal safety

The direct numerical simulation of pipe flow

The conservative difference scheme and the third-order Runge-Kutta scheme in combination with the the Crank-Nicholson scheme are used to directly simulate the flow field in a pipe with the Reynolds number of 2 600. The flow field, including the velocity distribution and the turbulence intensity, is obtained by the direct numerical simulation. From the calculated results, the ratio of the linear average velocity along the ultrasonic propagation path to the profile average velocity on the pipe cross-section is also obtained in an ultrasonic flow meter. It is concluded that the direct numerical simulation method can be used to study the ratio of the profile-linear average velocity at low Reynolds number conditions in the transition region and to improve the measurement accuracy of the ultrasonic flow meter.

DIRECT NUMERICAL SIMULATIONS OF TURBULENT CHANNEL FLOWS WITH CONSIDERATION OF THE BUOYANCY EFFECT OF THE BUBBLE PHASE

Turbulent channel flows with consideration of the buoyancy effect of the bubble phase is investigated by means of the Direct Numerical Simulation （DNS）. This two-phase system is solved by a two-way coupling Lagrangian-Eulerian approach. The Reynolds number based on the friction velocity and the half-width of the channel is 194, and the gravitational acceleration varies from -0.5 to 0.5, ranging from the upflow to the downflow cases. This study aims to reveal the influence of buoyancy on the turbulence behavior and the bubble motion. Some typical statistical quantities, including the averaged velocities and velocity fluctuations for the fluid and bubble phases, as well as the flow structures of the turbulence fluctuations, are analyzed.

Dynamics and geometry of developing planar jets based on the invariants of the velocity gradient tensor

Based on the direct numerical simulation （DNS）, the developing planar jets under different initial conditions, e.g., the con- ditions of the exit Reynolds number and the exit mean velocity profile, are investigated. We mainly focus on the characteristics of the invariants of the velocity gradient tensor, which provides insights into the evolution of the dynamics and the geometry of the planar jets along with the flow transition. The results show that the initial flow near the jet exit is strongly predominated by the dissipation over the enstrophy, the flow transition is accompanied by a severe rotation and straining of the flow elements, where the vortex structure evolves faster than the fluid element deformation, in the fully-developed state, the irrotational dissipation is dominant and the most probable geometry of the fluid elements should remain between the biaxial stretching and the axisymmetric stretching. In addition, with a small exit Re and a parabolic profile for the exit mean streamwise velocity, the decay of the mean flow field and the magnitude of the turbulent variables will be strengthened in the process of the flow transition, however, a large exit Re will promote the flow transition to the fully-developed state. The cross-impact between the exit Re and the exit mean velocity profile is also observed in the present study.