Critical Stokes Number for Gas-Solid Flow Erosion of Wind Turbine Airfoil

Wind turbine blades are inevitable to be eroded in wind-sand environment,so it is crucial to identify the flow conditions under which the erosion happens.Here,the effect of the sand diameter on wind turbine airfoil is first investigated.When the sand diameter is less than 3μm,the sands will bypass the airfoil and no erosion occurs.When the sand diameter is larger than 4μm,the sand grains collide with the airfoil and the erosion happens.Thus,there must be a critical sand diameter between 3μm and 4μm,at which the erosion is initiated on the airfoil surface.To find out this critical value,aparticle Stokes number is introduced here.According to the range of the critical sand diameter mentioned above,the critical value of particle Stokes number is reasonably assumed to be between 0.007 8and 0.014.The assumption is subsequently validated by other four factors influecing the erosion,i.e.,the angle of attack,relative thickness of the airfoil,different series airfoil,and inflow velocity.Therefore,the critical range of Stokes number has been confirmed.

Critical transition Reynolds number for plane channel flow

The determination of the critical transition Reynolds number is of practical importance for some engineering problems. However, it is not available with the current theoretical method, and has to rely on experiments. For supersonic/hypersonic boundary layer flows, the experimental method for determination is not feasible either. Therefore, in this paper, a numerical method for the determination of the critical transition Reynolds number for an incompressible plane channel flow is proposed. It is basically aimed to test the feasibility of the method. The proposed method is extended to determine the critical Reynolds number of the supersonic/hypersonic boundary layer flow in the subsequent papers.

Quark Number Susceptibility around the Chiral Critical End Point

We calculate the quark number susceptibility （QNS） around the chiral critical end point （CEP）. The CEP is found to be located at （μc,Tc）= （80 MeV, 148 MeV） where μc and Tc are the critical chemical potential and temperature, respectively. The QNS is found to have the highest and sharpest peak at the CEP. It is also found that, when the chemical potential μ is in the range of 60MeV≤ μ ≤ 110MeV, the QNS near the transition temperature is larger than the free field result, which indicates that the space-like damping mode dominates the degree of freedom of motion near the CEP.

从亚临界到临界雷诺数圆柱绕流和分离泡的大涡模拟

A large eddy simulation of wall-adapting local eddy-viscosity model(LES-WALE)is used to simulate the threedimensional flow around a circular cylinder with a diameter of 0.25 m from sub-critical to super-critical Reynolds numbers at 1×10^(5),2.5×10^(5),and 7.2×10^(5),respectively.The present results such as drag crisis,surface pressure distribution,and Strouhal number are in good agreement with the classical experimental data.When entering the critical region,a small plateau was found on the pressure distribution curves,corresponding to the appearance of laminar separation bubbles,and the separation point is delayed and the recirculation bubbles become narrowed and shortened.The tangential velocity of the cylinder surface changes from positive to negative at the separation point.The instantaneous vorticity and timeaveraging separation bubbles embody an unstable feature.Within the separation bubble,the pressure varies dramatically with time,but not with position.The surface pressure fluctuates greatly after the laminar separation bubble appears,and it is gradually stabilized until the basic pressure is reached.The process of laminar separation,transition from laminar flow to turbulent flow and turbulent reattachment is also shown.The three-dimensional Q criterion of vortex structure and the two-dimensional spanwise vorticity reveal the phenomenon that the wake structure narrows with the increase of the Reynolds number.

Critical Transition Reynolds Number for the Incompressible Flat-plate Boundary Layer as Searched by Numerical Simulation

The critical transition Reynolds number is the lowest value at which the turbulent flow can hold in real flows.The determination of the critical transition Reynolds number not only is a scientific problem,but also is important for some engineering problems.However,there is no available theoretical method to search the critical value.For the hypersonic boundary layer with significant importance for engineering problems,there is no available experimental method to search the critical value so far.Consequently,it is imperative to take numerical method to search it.In this paper,direct numerical simulations(DNS)method is employed to determine the critical transition Reynolds number for the incompressible flat-plate boundary layer.Firstly,under the assumption of parallel flow,the temporal mode DNS is performed to determine the critical value as Re_(xpcr)=43767,which is quite close to the numerical results of other people.Secondly,under the condition of nonparallel flow,the spatial mode DNS is performed to determine the critical transition Reynolds number as Re_(xcr)=3×10^(5),which is well consistent with the experimental results.In principle,the proposed method in this paper can be extended to the supersonic/hypersonic boundary layer,and that problem will be discussed in the subsequent papers.

Viscous Fluid Flowing around the Plate: Turbulence

Diffusion of momentum gives rise to viscosity. This article presents a solution in the explicit form of the equation of the momentum diffusion for a viscous fluid flowing around a plate taking into account deceleration. Three characteristic regions of a viscous flow have been described: the mantle, the body of the boundary layer, the viscous sublayer. In the mantle, the effect of viscosity is significant at a considerable distance from the plate. The momentum diffusion is focused in the body of the boundary layer. The diffusion force that produces the momentum of force giving rise to eddies is localized in the viscous sublayer. At the beginning of the plate, a moment of force twists the liquid along the flow, creating eddies that roll along the plate. For this reason, they are pressed against the surface of the plate. But at some distance from the beginning of the plate, the moment of force changes its orientation to the opposite and twists the vortices in the opposite direction, causing the vortices to roll along the plate against the flow. This causes the liquid to detach from the surface of the plate. This is the beginning of turbulence. The diameter of the vortex produced in the viscous sublayer is small being of the order of the thickness of the viscous sublayer. The vortex possesses a large angular velocity. Due to the momentum diffusion and the effect of the eddies combined in passing along the plate, its diameter increases up to the size of the thickness of the boundary layer and even more, whereas its angular velocity decreases down to the values really observed. The value of the critical Reynolds number of the transition from the laminar flow to the turbulent one has been found, and it agrees with the experimental data. The value of the shear stress produced by the viscous fluid on the plate surface has also been obtained. The way of measurement of the friction coefficient characterizing the effect of the plate on the flow has been proposed. It has been shown that the boundary condition of adhesion to the surface of a body flown around, that is applied in the estimation of viscous flows, contradicts the real processes of the flow.

Orifice plate cavitation mechanism and its influencing factors

The orifice plate energy dissipater is an economic and highly efficient dissipater. However, there is a risk of cavitaion around the orifice plate flow： In order to provide references for engineering practice, we examined the cavitation mechanism around the orifice plate and its influencing factors by utilizing mathematical analysis methods to analyze the flow conditions around the orifice plate in view of gas bubble dynamics. Through the research presented in this paper, the following can be observed： The critical radius and the critical pressure of the gas nucleus in orifice plate flow increase with its initial state parameter r0 ; the development speed of bubbles stabilizes at a certain value after experiencing a peak value and a small valley value; and the orifice plate cavitation is closely related to the distribution of the gas nucleus in flow. For computing the orifice plate cavitation number, we ought to take into account the effects of pressure fluctuation. The development time of the gas nucleus from the initial radius to the critical radius is about 107-10-5 s; therefore, the gas nucleus has sufficient time to develop into bubbles in the negative half-cycle of flow fluctuation. The orifice critical cavitation number is closely related to the orifice plate size, and especially closely related with the ratio of the orifice plate radius to the tunnel radius. The approximate formula for the critical cavitation number of the square orifice plate that only considers the main influencing factor was obtained by model experiments.

Coalescence dynamics of two droplets of different viscosities in T-junction microchannel with a funnel-typed expansion chamber

The coalescence behavior of two droplets with different viscosities in the funnel-typed expansion chamber in T-junction microchannel was investigated experimentally and compared with droplet coalescence of the same viscosity.Four types of coalescence regimes were observed:contact non-coalescence,squeeze non-coalescence,two-droplet coalescence and pinch-off coalescence.For droplet coalescence of different viscosities,the operating range of non-coalescence becomes narrowed compared to the droplet coalescence of same viscosity,and it shrinks with increasing viscosity ratioηof two droplets,indicating that the difference in the viscosity of two droplets is conducive to coalescence,especially when1<η<6.Furthermore,the influences of viscosity ratio and droplet size on the film drainage time(Tdr)and critical capillary number(Ca)c)were studied systematically.It was found that the film drainage time declined with the increase of average droplet size,which abided by power-law relation with the size difference and viscosity ratio of the two droplets:Tdr~(ld)^(0.25±0.04)and Tdr~(η)^(﹣0.1±0.02).For droplet coalescence of same viscosity,the relation of critical capillary number with two-phase viscosity ratio and dimensionless droplet size is Cac=0.48λ^(0.26)l^(﹣2.64),while for droplet coalescence of different viscosities,the scaling of critical capillary number with dimensionless average droplet size,dimensionless droplet size difference and viscosity ratio of two droplets is Cac=0.11η^(﹣0.07)ls^(﹣2.23)l^(0.16)_(d).

Locating the QCD critical end point through peaked baryon number susceptibilities along the freeze-out line

We investigate the baryon number susceptibilities up to fourth order along different freeze-out lines in a holographic QCD model with a critical end point（CEP）, and we propose that the peaked baryon number susceptibilities along the freeze-out line can be used as a clean signature to locate the CEP in the QCD phase diagram.On the temperature and baryon chemical potential plane, the cumulant ratio of the baryon number susceptibilities（up to fourth order） forms a ridge along the phase boundary, and develops a sword-shaped ＂mountain＂ standing upright around the CEP in a narrow and oblate region. The measurement of baryon number susceptibilities from heavy-ion collision experiments is along the freeze-out line. If the freeze-out line crosses the foot of the CEP mountain, then one can observe the peaked baryon number susceptibilities along the freeze-out line, and the kurtosis of the baryon number distributions has the highest magnitude. The data from the first phase of the beam energy scan program at the Relativistic Heavy Ion Collider indicates that there should be a peak of the kurtosis of the baryon number distribution at a collision energy of around 5 Ge V, which suggests that the freeze-out line crosses the foot of the CEP mountain and the summit of the CEP should be located nearby, around a collision energy of 3–7 GeV.

Convective stability analysis of a micropolar fluid layer by variational method

This paper studies Rayleigh-B’enard convection of micropolar fluid layer heated from below with realistic boundary conditions.A specific approach for stability analysis of a convective problem based on variational principle is applied to characterize the Rayleigh number for quite general nature of bounding surfaces.The analysis consists of replacing the set of field equations by a variational principle and the expressions for Rayleigh number are then obtained by using trial function satisfying the essential boundary conditions.Further,the values of the Rayleigh number for particular cases of large and small values of the microrotation coefficient have been obtained.The effects of wave number and micropolar parameter on the Rayleigh numbers for onset of stationary instability for each possible combination of the bounding surfaces are discussed and illustrated graphically.The present analysis establishes that the nature of bounding surfaces combination and microrotation have significant effect on the onset of convection.

The Analysis of Stability of Bingham Fluid Flowing Down an Inclined Plane

In this paper, the stability problem of Bingham fluids flowing down an inclinedplane is studied with respect to two dimensional disturbances. The crilical Reynolodsnumber is given in ihe case of long waves, and the effect of yield stress on stability isanalysed.

Laminar Flame Instability of n-Hexane,n-Octane,and n-Decane in Spherical Expanding Flames

This paper focuses on the laminar flame instability of three high molecular weight n-alkanes,namely n-hexane,n-octane,and n-decane.The experiment was carried out in a constant volume combustion bomb to get the flame images.The critical radius under different conditions was extracted using the image processing program.Combined with the existing critical Peclet number theory,the dominant factors of flame instability under current conditions for three n-alkanes can be figured out.Moreover,the average cell size(equivalent cell radius,R_(cell))was extracted to provide quantitative analysis of the flame cellular structure,based on the method developed in this work.The theoretical R_(cell)were also calculated and compared with the experimental results to validate the proposed method.

EXPERIMENTAL STUDY OF FRICTION FACTOR FOR GROUNDWATER FLOW IN A SINGLE ROUGH FRACTURE

To study the relationships between the friction factor f and the flow type in a single rough fracture, the formulae of f for both unconfined and confined flows are deduced based on previous studies. The relationships between f and the Reynolds number （Re） for different relative roughnesses are investigated experimentally. The Moody-type diagram, based on the deduced formula of f, is also plotted and the hydraulic characteristics of the flow in a rough fracture are analyzed. Results show that the Moody-type diagram of the experiment has a similar distribution to that of the conventional Moody diagram. It is found that the value of f in the experiment is much smaller than that of the conventional Moody diagram and turbulent flow appears easier for rough fractures, which can be explained by the separation phenomenon in boundary layers. The critical Re ranging from 650 to 700 in rough fractures is concluded based on the experimental results. It also can be concluded that the friction factor f is related not only with the Re and the relative roughness but also with the absolute roughness.