圆形直管湍流光滑管区的摩擦因数计算

提出了一个简单适用的显式方程,可用于流体处在直圆形湍流光滑管区时的摩擦因数的计算。在3000≤Re≤108时,该方程的计算结果与卡门-普朗特方程的平均偏差为0.503%,最大偏差不超过0.87%,其计算精度高于其他现有的经验公式。

圆形直管湍流粗糙管区的摩擦因数计算

提出了简单适用的管路摩擦因数λ的计算方法,解决了柯尔布鲁克-怀特(Colebrook-white)公式中隐函数不易计算的问题,该方法计算精度高于其他现有简化公式;在10-8≤ε/d≤0.05、3000≤Re≤108时,该方程的计算结果与原Colebrook-White方程的平均相对偏差为0.3315%,最大偏差不超过2.0%。

搅拌釜内气液动态相对稳定状态及雷诺数计算

设计了一种全新的测量搅拌釜内动态体系温度、压力的设备及方案,以CO2-H2O体系为特征体系,对搅拌状态下气液二相动态相对稳定状态进行了研究。搅拌釜内流体发生湍动,达到相平衡的气液二相的传质极限发生偏移。在不同温度(288.15—307.15 K)和转速(190—520 r/min)下进行实验,温度一定时,动态平衡偏移率随湍动程度的增大而增大。温度为307.15 K、转速为520 r/min时,动态平衡偏移率高达27%。利用动态平衡偏移率反算搅拌釜内的雷诺数ReD,将其与传统计算得出的雷诺数ReT进行比较拟合得到关联式。

Comparison of the Reynolds-averaged Turbulence Models on Single Phase Flow Simulation in Agitated Extraction Columns

The flow field of liquid phase (water) of agitated extraction columns is simulated with the help of computational fluid dynamics (CFD). Four kinds of Reynolds-averaged turbulence models, i.e. the standard k-ε model, the RNG (renormalization group) k-s model, the realizable k-ε model and the Reynolds stress model, are compared in detail in order to judge which is the best model in terms of the accuracy, less CPU time and memory required. The performance of the realizable k-s model is obviously improved by reducing the model constant from C2 = 1.90 to C2 = 1.61. It is concluded that the improved realizable k-e model is the optimal model.

Turbulent boundary layers and hydrodynamic flow analysis of nanofluids over a plate

A numerical analysis of the log-law behavior for the turbulent boundary layer of a wall-bounded flow is performed over a flat plate immersed in three nanofluids(Zn O-water,SiO_(2)-water,TiO_(2)-water).Numerical simulations using CFD code are employed to investigate the boundary layer and the hydrodynamic flow.To validate the current numerical model,measurement points from published works were used,and the compared results were in good compliance.Simulations were carried out for the velocity series of 0.04,0.4 and 4 m/s and nanoparticle concentrations0.1% and 5%.The influence of nanoparticles’ concentration on velocity,temperature profiles,wall shear stress,and turbulent intensity was investigated.The obtained results showed that the viscous sub-layer,the buffer layer,and the loglaw layer along the potential-flow layer could be analyzed based on their curving quality in the regions which have just a single wall distance.It was seen that the viscous sub-layer is the biggest area in comparison with other areas.Alternatively,the section where the temperature changes considerably correspond to the thermal boundary layer’s thickness goes a downward trend when the velocity decreases.The thermal boundary layer gets deep away from the leading edge.However,a rise in the volume fraction of nanoparticles indicated a minor impact on the shear stress developed in the wall.In all cases,the thickness of the boundary layer undergoes a downward trend as the velocity increases,whereas increasing the nanoparticle concentrations would enhance the thickness.More precisely,the log layer is closed with log law,and it is minimal between Y^(+)=50 and Y^(+)=95.The temperature for nanoparticle concentration φ=5%is higher than that for φ=0.1%,in boundary layers,for all studied nanofluids.However,it is established that the behavior is inverted from the value of Y^(+)=1 and the temperature for φ =0.1% is more important than the case of φ =5%.For turbulence intensity peak,this peak exists at Y^(+)=100 for v=4 m/s,Y^(+)=10 for v=0.4 m/s and Y^(+)=8 for v=0.04 m/s.

Comparison of a Full Second-Order Moment Model and an Algebraic Stress Two-Phase Turbulence Model for Simulating Bubble-Liquid Flows in a Bubble Column

A full second-order moment (FSM) model and an algebraic stress (ASM) two-phase turbulence modelare proposed and applied to predict turbulent bubble-liquid flows in a 2D rectangular bubble column. Predictiongives the bubble and liquid velocities, bubble volume fraction, bubble and liquid Reynolds stresses and bubble-liquidvelocity correlation. For predicted two-phase velocities and bubble volume fraction there is only slight differencebetween these two models, and the simulation results using both two models are in good agreement with the particleimage velocimetry (PIV) measurements. Although the predicted two-phase Reynolds stresses using the FSM are insomewhat better agreement with the PIV measurements than those predicted using the ASM, the Reynolds stressespredicted using both two models are in general agreement with the experiments. Therefore, it is suggested to usethe ASM two-phase turbulence model in engineering application for saving the computation time.

Deviation of Carbon Dioxide-Water Gas-Liquid Balance from Thermodynamic Equilibrium in Turbulence h Experiment and Correlation

The carbon dioxide-water system was used to investigate the flowing gas-liquid metastable state. The experiment was carded out in a constant volume vessel with a horizontal circulation pipe and a peristaltic pump forced CO2 saturated water to flow. The temperature and pressure were recorded. The results showed that some CO2 escaped from the water in the flow process and the pressure increased, indicating that the gas-liquid equilibrium was broken. The amount of escaped CO2 varied with flow speed and reached a limit in a few minutes, entitled dy- namic equilibrium. Temperature and liquid movement played the same important role in breaking the phase equilib- rium. Under the experimental conditions, the ratio of the excessive carbon dioxide in the gas phase to its thermody- namic equilibrium amount in the liquid could achieve 15%.

Simulation and analysis of airflow stability during fire in mine belt roadway

According to fluid dynamics analysis during the fire, the criteria k-ε two-equation model for solving three-dimensional turbulence was determined, the pollutants generated in the fire disaster were set by adopting Mixture multiphase flow, and the SIMPLE algorithm was used for solving on the basis of comprehensive consideration on the heat radia- tion and components transmission during fire. By simulating the airflow flowing state inside the tunnel during fire disaster of downward ventilation, drift ventilation, and ascensional ventilation, respectively, with regard to the actual situation of No.l, No.3, and No.5 belt roadway in Kongzhuang Coal Mine, the velocity vector distributions of pollutants under different inlet air volumes were obtained, and the damage degree and influential factors of disaster were also clear, which is helpful to control and avoid disaster during belt roadway fire.

Validation of the RANS-SOM Combustion Model Using Direct Numerical Simulation of Incompressible Turbulent Reacting Flows

The second-order moment combustion model, proposed by the authors is validated using the direct numerical simulation （DNS） of incompressible turbulent reacting channel flows. The instantaneous DNS results show the near-wall strip structures of concentration and temperature fluctuations. The DNS statistical results give the budget of the terms in the correlation equations, showing that the production and dissipation terms are most important. The DNS statistical data are used to validate the closure model in RANS second-order moment （SOM） combustion model. It is found that the simulated diffusion and production terms are in agreement with the DNS data in most flow regions, except in the near-wall region, where the near-wall modification should be made, and the closure model for the dissipation term needs further improvement. The algebraic second-order moment （ASOM） combustion model is well validated by DNS.

Numerical Simulation on Ship Bubbly Wake

Based on a volume of fluid two-phase model imbedded in the general computational fluid dynamics code FLUENT6.3.26, the viscous flow with free surface around a model-scaled KRISO container ship （KCS） was first numerically simulated. Then with a rigid-lid-free-surface method, the underwater flow field was computed based on the mixture muitiphase model to simulate the bubbly wake around the KCS hull. The realizable k-e two-equation turbulence model and Reynolds stress model were used to analyze the effects of turbulence model on the ship bubbly wake. The air entrainment model, which is relative to the normal velocity gradient of the free surface, and the solving method were verified by the qualitatively reasonable computed results.

CFD Simulation of Orifice Flow in Orifice-type Liquid Distributor

In this study,a suitable CFD(computational fluid dynamics)model has been developed to investigate the influence of liquid height on the discharge coefficient of the orifice-type liquid distributors.The orifice flow in different diameters and liquid heights has been realized using the shear stress transport(SST)turbulence model and the Gamma Theta transition(GTT)model.In the ANSYS CFX software,two models are used in conjunction with an automatic wall treatment which allows for a smooth shift from a wall function(WF)to a low turbulent-Re near wall formulation(LTRW).The results of the models coupled with LTRW are closer to the experimental results compared with the models with WF,indicating that LTRW is more appropriate for the prediction of boundary layer characteristics of orifice flow.Simulation results show that the flow conditions of orifices change with the variation of liquid height.With respect to the turbulence in orifice,the SST model coupled with LTRW is recommended.However,with respect to the transition to turbulence in orifice with an increase in liquid height,the predictions of GTT model coupled with LTRW are superior to those obtained using other models.

Experimental and CFD Studies on the Performance of Microfiltration Enhanced by a Turbulence Promoter

This paper reports experimental and computational fluid dynamics(CFD) studies on the performance of microfiltration enhanced by a helical screw insert.The experimental results show that the use of turbulence pro-moter can improve the permeate flux of membrane in the crossflow microfiltration of calcium carbonate suspension,and flux improvement efficiency is strongly influenced by operation conditions.The energy consumption analysis indicates that the enhanced membrane system is more energy saving at higher feed concentrations.To explore the intrinsic mechanism of flux enhancement by a helical screw insert,three-dimensional CFD simulation of fluid flow was implemented.It reveals that hydrodynamic characteristics of fluid flow inside the channel are entirely changed by the turbulence promoter.The rotational flow pattern increases the scouring effect on the tube wall,reducing the particle deposition on the membrane surface.The absence of stagnant regions and high wall shear stress are respon-sible for the enhanced filtration performance.No secondary flow is generated in the channel,owing to the streamline shape of helical screw insert,so that the enhanced performance is achieved at relatively low energy consumption.

A new approach to quantifying vehicle induced turbulence for complex traffic scenarios

Traffic-related pollutants adversely affect air quality, especially in regions near major roadways. The vehicleinduced turbulence(VIT) is a significant factor that controls the initial dilution, dispersion, and ultimately the chemical and physical fate of pollutants by altering the conditions in the microenvironment. This study used a computational fluid dynamics(CFD) software FLUENT to model the vehicle-induced turbulence(VIT) generated on roadways, with a focus on impact of vehicle-vehicle interactions, traffic density and vehicle composition on turbulent kinetic energy(TKE). We show, for the first time, that the overall TKE from multiple vehicles traveling in series can be estimated by superimposing the TKE of each vehicle, without considering the distance between them while the distance is greater than one vehicle length. This finding is particularly significant since it enables a new approach to VIT simulations where the overall TKE is calculated as a function of number of vehicles. We found that the interactions between vehicles traveling next to each other in adjacent lanes are insignificant,regardless the directions of the traffic flow. Consequently, simulations of different traffic scenarios can be substantially simplified by treating two-way traffic as one-way traffic, with less than 5% difference in the overall volume-averaged TKE. We also developed equations that allow the estimation of the overall volume-averaged TKE as a function of the number and the type of vehicles.

Introduction to the Turbulent Flows Theory An Axially-Symmetric Peaceful Flows

This paper introduces to fluid state physics （fluid mechanics） a new interpretation of physical phenomena taking place in a fluid in motion. It introduces the base of a new theory claiming that every flow has its own internal structure of motion, which is definite organization of motion, rather than a ＂molecular chaos＂, known from the fluid statics. The paper introduces the new notion of structures vector fields of power and momentum and shows every Newtonian fluid flows are dual in character. It shows that the flow of Newtonian fluid has a dual character. It demonstrates on models and further in mathematical interpretation of physical phenomena. It introduces, on the one hand, the cycloidal motion model into the fluid mechanics, ad on the other hand an addition to the known, the classical model of Poiseuille laminar motion. The theory of dualism （double nature of physical phenomena） allows the description of selected characteristics of the flow, either by using the theory ofcycloidal motion （semicycloidal）, or by using the supplemented theory of laminar motion. The dualism theory is useful to describe each type of flows both, laminar and turbulent. This paper is only an introduction to the theory. It has been assigned number 1. It has been granted a high priority, since it contains basic concepts that will be used in others, following papers of long cycle.

MHD from a Microscopic Concept and Onset of Turbulence in Hartmann Flow

We derive higher order magneto-hydrodynamic （MHD） equations from a microscopic picture using pro-jection and perturbation formalism. In an application to Hartmann flow we find velocity profiles flattening towards the center at the onset of turbulence in hydrodynamic limit. Comparison with the system under the effect of a uniform magnetic field yields difference in the onset of turbulence consistent with observations, showing that the presence of magnetic field inhibits onset of instability or turbulence. The laminar-turbulent transition is demonstrated in a phase transition plot of the development in time of the relative average velocities vs. Reynolds number showing a sharp increase of the relative average velocity at the transition point as determined by the critical Reynolds number.

Numerical Investigation on Two-dimensional Boundary Layer Flow with Transition

As a basic problem in many engineering applications, transition from laminar to turbulence still remains a difficult problem in computational fluid dynamics （CFD）. A numerical study of one transitional flow in two-dimensional is conducted by Reynolds averaged numerical simulation （RANS） in this paper. Turbulence model plays a significant role in the complex flows＇ simulation, and four advanced turbulence models are evaluated. Numerical solution of frictional resistance coefficient is compared with the measured one in the transitional zone, which indicates that Wilcox （2006） k-ω model with correction is the best candidate. Comparisons of numerical and analytical solutions for dimensionless velocity show that averaged streamwise dimensionless velocity profiles correct the shape rapidly in transitional region. Furthermore, turbulence quantities such as turbulence kinetic energy, eddy viscosity, and Reynolds stress are also studied, which are helpful to learn the transition＇s behavior.

Study on heat transfer for falling liquid film flow with consideration of interfacial evaporation

The interfacial evaporative heat transfer was included in the semi-empirical study of the heat transfer for the falling liquid film flow. The investigations showed that, the inclusion of the interfacial eveiporative heat transfer in the turbulent model would lower the predicted convective heat transfer coefficient. Predictions of the new model resulted in a prominent deviation from that predictions of the normal model in the case of large mass flow rate and low wall heat flux. This deviation will be decreased with increasing wall heat flux, such that it will be asymptotic zero at very high wall heat flux. Predictions of the new model agreed well with the current experimental measurements. This study has verified that the Reynolds number is not the sole crucial parameter for heat transfer of falling liquid film flow, and wall heat flux will be another important independent parameter. This result is consistent with our previous studies.

CFD simulation of bubbly turbulent Tayor–Couette flow

Bubbly gas-liquid Taylor-Couette vortex flow has been the subject of several recent investigations both because of interest in bubble-induced drag reduction and because such devices have potential applications to a variety of chemical and biochemical processing problems. In order to quantitatively describe the hydrodynamics of highly turbulent two phase Taylor-Couette flow, a rigorous two-fluid computational fluid dynamics （CFD） model was developed and compared with previously published experimental data. This model includes a comprehensive description of the constitutive closure for inter-phase forces and turbulence was simulated using both the k- and k-to models. In addition, the mechanism by which the dispersed fluid attains a non-uniform radial and axial distribution is analyzed and the relative importance of various interphase forces is discussed. Lastly the model was validated by comparison of simulation predictions with experimental data, and it is shown that the CFD model correctly predicts phase velocity, velocity fluctuation, and gas distribution, and may provide guidance for reactor design and scale-up.

Analysis of Turbulent Heat Transfer and Fluid Flow in Channels with Various Ribbed Internal Surfaces

This paper presents results of a numerical investigation of heat transfer and flow pattern characteristics of a channel with repeated ribs on one broad wall. Numerical computations are performed for seven ribs placed on the bottom wall of a channel for Reynolds numbers ranging from 10,000 to 30,000. The newly modified ribs (the ones with convex pointing upstream/downstream rib, wedge pointing upstream/downstream rib, concave pointing upstream/downstream rib and also concave-concave rib as well as convex-concave rib), are proposed for simulation with prospect to reduce flow separation and extend reattachment area compared to the unmodified square rib. The numerical results are reported in forms of flow structure, temperature field, turbulent kinetic energy, Nusselt number, friction factor and thermal enhancement factor. The results indicate the rib with concave-concave surfaces efficiently suppresses flow separation bubble in the corner of the rib and induces large recirculation zone over those of the others, hence giving the highest Nusselt number and friction factor. On the other hand, the one with convex-concave surface provides the lowest friction factor with moderate Nusselt number. Due to the prominent effect of its low friction factor, the rib with convex-concave surface offers the highest thermal enhancement factor of 1.19.

Hypersonic flow control of shock wave/turbulent boundary layer interactions using magnetohydrodynamic plasma actuators

The effect of magnetohydrodynamic(MHD)plasma actuators on the control of hypersonic shock wave/turbulent boundary layer interactions is investigated here using Reynolds-averaged Navier-Stokes calculations with low magnetic Reynolds number approximation.A Mach 5 oblique shock/turbulent boundary layer interaction was adopted as the basic configuration in this numerical study in order to assess the effects of flow control using different combinations of magnetic field and plasma.Results show that just the thermal effect of plasma under experimental actuator parameters has no significant impact on the flow field and can therefore be neglected.On the basis of the relative position of control area and separation point,MHD control can be divided into four types and so effects and mechanisms might be different.Amongst these,D-type control leads to the largest reduction in separation length using magnetically-accelerated plasma inside an isobaric dead-air region.A novel parameter for predicting the shock wave/turbulent boundary layer interaction control based on Lorentz force acceleration is then proposed and the controllability of MHD plasma actuators under different MHD interaction parameters is studied.The results of this study will be insightful for the further design of MHD control in hypersonic vehicle inlets.