Large eddy simulation of the rotation effect on the ocean turbulence kinetic energy budget in the surface mixed layer

A non-hydrostatic, Boussinesq, and three-dimensional large eddy simulation(LES) model was used to study the impact of the Earth's rotation on turbulence and the redistribution of energy in turbulence kinetic energy(TKE) budget. A set of numerical simulations was conducted,(1) with and without rotation,(2) at different latitudes(10°N, 30°N, 45°N, 60°N, and 80°N),(3) with wave breaking and with Langmuir circulation, and(4) under different wind speeds(5, 10, 20, and 30 m/s). The results show that eddy viscosity decreases when rotation is included, indicating that rotation weakens the turbulence strength. The TKE budget become tight with rotation and the effects of rotation grow with latitude. However, rotation become less important under Langmuir circulation since the transport term is strong in the vertical direction. Finally, simulations were conducted based on field data from the Boundary Layer and Air-Sea Transfer Low Wind(CBLAST-Low) experiment. The results, although more complex, are consistent with the results obtained from earlier simulations using ideal numerical conditions.

Turbulence Intensity and Turbulent Kinetic Energy Parameters over a Heterogeneous Terrain of Loess Plateau

A deep understanding of turbulence structure is important for investigating the characteristics of the atmospheric boundary layer,especially over heterogeneous terrain.In the present study,turbulence intensity and turbulent kinetic energy（TKE）parameters are analyzed for different conditions with respect to stability,wind direction and wind speed over a valley region of the Loess Plateau of China during December 2003 and January 2004.The purpose of the study is to examine whether the observed turbulence intensity and TKE parameters satisfy Monin–Obukhov similarity theory（MOST）,and analyze the wind shear effect on,and thermal buoyancy function of,the TKE,despite the terrain heterogeneity.The results demonstrate that the normalized intensity of turbulence follows MOST for all stability in the horizontal and vertical directions,as well as the normalized TKE in the horizontal direction.The shear effect of the wind speed in the Loess Plateau region is strong in winter and could enhance turbulence for all stability conditions.During daytime,the buoyancy and shear effect together constitute the generation of TKE under unstable conditions.At night,the contribution of buoyancy to TKE is relatively small,and mechanical shearing is the main production form of turbulence.

Turbulent Kinetic Energy of Flow during Inhale and Exhale to Characterize the Severity of Obstructive Sleep Apnea Patient

This paper aims to investigate and present the numerical investigation of airflow characteristics using Turbulent Kinetic Energy(TKE)to characterize the upper airway with obstructive sleep apnea(OSA)under inhale and exhale breathing conditions.The importance of TKE under both breathing conditions is that it showan accuratemethod in expressing the severity of flow in sleep disorder.Computational fluid dynamics simulate the upper airway’s airflow via steady-state Reynolds-averaged Navier-Stokes(RANS)with k–ωshear stress transport(SST)turbulencemodel.The three-dimensional(3D)airway model is created based on the CT scan images of an actual patient,meshed with 1.29 million elements using Materialise Interactive Medical Image Control System(MIMICS)and ANSYS software,respectively.High TKE were noticed around the region after the necking(smaller cross-sectional area)during the inhale and exhale breathing.The turbulent kinetic energy could be used as a valuablemeasure to identify the severity of OSA.This study is expected to provide a better understanding and clear visualization of the airflow characteristics during the inhale and exhale breathing in the upper airway of patients for medical practitioners in the OSA research field.

Wave Attenuation and Turbulence Driven by Submerged Vegetation Under Current-Wave Flow

A set of laboratory experiments are carried out to investigate the effect of following/opposing currents on wave attenuation.Rigid vegetation canopies with aligned and staggered configurations were tested under the condition of various regular wave heights and current velocities,with the constant water depth being 0.60 m to create the desired submerged scenarios.Results show that the vegetation-induced wave dissipation is enhanced with the increasing incident wave height.A larger velocity magnititude leads to a greater wave height attenuation for both following and opposing current conditions.Moreover,there is a strong positive linear correlation between the damping coefficientβand the relative wave height H_(0)/h,especially for pure wave conditions.For the velocity profile,the distributions of U_(min)and U_(max)show different patterns under combined wave and current.The time-averaged turbulent kinetic energy(TKE)vary little under pure wave and U_(c)=±0.05 m/s conditions.With the increase of flow velocity amplitude,the time-averaged TKE shows a particularly pronounced increase trend at the top of the canopy.The vegetation drag coefficients are obtained by a calibration approach.The empirical relations of drag coefficient with Reynolds and Keulegane-Carpenter numbers are proposed to further understand the wave-current-vegetation interaction mechanism.

Wavelet spectrum analysis on energy transfer of multi-scale structures in wall turbulence

The streamwise velocity components at different vertical heights in wall turbulence were measured. Wavelet transform was used to study the turbulent energy spectra, indicating that the global spectrum results from the weighted average of Fourier spectrum based on wavelet scales. W＇avelet transform with more vanishing moments can express the declining of turbulent spectrum. The local wavelet spectrum shows that the physical phenomena such as deformation position in the boundary layer, and the or breakup of eddies are related to the vertical energy-containing eddies exist in a multi-scale form. Moreover, the size of these eddies increases with the measured points moving out of the wall. In the buffer region, the small scale energy-containing eddies with higher frequency are excited. In the outer region, the maximal energy is concentrated in the low-frequency large-scale eddies, and the frequency domain of energy-containing eddies becomes narrower.

Estimation of Turbulent Kinetic Energy Dissipation Rate in the Bottom Boundary Layer of the Pearl River Estuary

A structure function approach is applied to estimate the turbulent kinetic energy （TKE） dissipation rate in the bottom boundary layer of the Pearl River Estuary （PRE）. Simultaneous measurements with an acoustic Doppler velocimeter （ADV） supplied independent data for the verification of the structure function method. The results show that, 1） the structure function approach is reliable and successfully applied method to estimate the TKE dissipation rate. The observed dissipation rates range between 8.3 ×10^-4 W/kg and 4.9× 10^-6 W/kg in YM01 and between 3.4×10^-4 W/kg and 4.8×10^-7 W/kg in YM03, respectively, while exhibiting a strong quarter-diurnal variation. 2） The balance between the shear production and viscous dissipation is better achieved in the straight river. This first-order balance is significantly broken in the estuary by non-shear production/dissipation due to wave-induced fluctuations.

Significance of including lid thickness and particle shape factor in numerical modeling for prediction of particle trap efficiency of invert trap

Sediment accumulation on the bed of open sewers and drains reduces hydraulic efficiency and can cause localized flooding.Slotted invert traps installed underneath the bed of open sewers and drains can eliminate sediment build-up by catching sediment load.Previous three-dimensional(3D)computational studies have examined the particle trapping performance of invert traps of different shapes and depths under varied sediment and flow conditions,considering particles as spheres.For two-dimensional and 3D numerical modeling,researchers assumed the lid geometry to be a thin line and a plane,respectively.In this 3D numerical study,the particle trapping efficiency of a slotted irregular hexagonal invert trap fitted at the flume bottom was examined by incorporating the particle shape factor of non-spherical sewage solid particles and the thicknesses of upstream and downstream lids over the trap in the discrete phase model of the ANSYS Fluent 2020 R1 software.The volume of fluid(VOF)and the realizable k-turbulence models were used to predict the velocity field.The two-dimensional particle image velocimetry(PIV)was used to measure the velocity field inside the invert trap.The results showed that the thicknesses of upstream and downstream lids affected the velocity field and turbulent kinetic energy at all flow depths.The joint impact of the particle shape factor and lid thickness on the trap efficiency was significant.When both the lid thickness and particle shape factor were considered in the numerical modeling,trap efficiencies were underestimated,with relative errors of-8.66%to-0.65%in comparison to the experimental values of Mohsin and Kaushal(2017).They were also lower than the values predicted by Mohsin and Kaushal(2017),which showed an overall overestimation with errors of-2.3%to 17.4%.

Research on blade tip clearance cavitation and turbulent kinetic energy characteristics of axial flow pump based on the partially-averaged Navier-Stokes model

To reveal the cavitation forms of tip leakage vortex(TLV)of the axial flow pump and the flow mechanism of the flow field,this research adopts the partially-averaged Navier-Stokes(PANS)model to simulate the cavitation values of an axial flow pump,followed by experimental validation.The experimental result shows that compared with the shear stress transport(SST)k-ωmodel,the PANS model significantly reduces the eddy viscosity of the flow field to make the vortex structure clearer and allow the turbulence scale to be more robustly analyzed.The cavitation area within the axial flow pump mainly comprises of TLV cavitation,clearance cavitation and tip leakage flows combined effect of triangular cloud cavitation formed.The formation and development of cavitation are accompanied by the formation and evolution of vortex,and variations in vortex structure also generate and promote the development of cavitation.In addition,an in-depth analysis of the relationship between the turbulent kinetic energy(TKE)transport equation and cavitation patterns was also conducted,finding that the regions with relatively high TKE are mainly distributed around gas/liquid boundaries with serious cavitation and evident gas-liquid change.This phenomenon is mainly attributed to the combined effect of the pressure action term,stress diffusion term and TKE production term.

Analysis of Turbulence Structure in the Stirred Tank with a Deep Hollow Blade Disc Turbine by Time-resolved PIV

The turbulence structure in the stirred tank with a deep hollow blade(semi-ellispe) disc turbine(HEDT) was investigated by using time-resolved particle image velocimetry(TRPIV) and traditional PIV.In the stirred tank,the turbulence generated by blade passage includes the periodic components and the random turbulent ones.Traditional PIV with angle-resolved measurement and TRPIV with wavelet analysis were both used to obtain the random turbulent kinetic energy as a comparison.The wavelet analysis method was successfully used in this work to separate the random turbulent kinetic energy.The distributions of the periodic kinetic energy and the random turbulent kinetic energy were obtained.In the impeller region,the averaged random turbulent kinetic energy was about 2.6 times of the averaged periodic one.The kinetic energies at different wavelet scales from a6 to d1 were also calculated and compared.TRPIV was used to record the sequence of instantaneous velocity in the impeller stream.The evolution of the impeller stream was observed clearly and the sequence of the vorticity field was also obtained for the identification of vortices.The slope of the energy spectrum was approximately-5/3 in high frequency representing the existence of inertial subrange and some isotropic properties in stirred tank.From the power spectral density(PSD) ,one peak existed evidently,which was located at f0(blade passage frequency) generated by the blade passage.

Analysis on the Interaction between Turbulence and Secondary Circulation of the Surface Layer at Jinta Oasis in Summer

The kinetic energy variations of mean flow and turbulence at three levels in the surface layer were calculated by using eddy covariance data from observations at Jinta oasis in 2005 summer. It is found that when the mean horizontal flow was stronger, the turbulent kinetic energy was increased at all levels, as well as the downward mean wind at the middle level. Since the mean vertical flow on the top and bottom were both negligible at that time, there was a secondary circulation with convergence in the upper half and divergence in the lower half of the column. After consideration of energy conversion, it was found that the interaction between turbulence and the secondary circulation caused the intensification of each other. The interaction reflected positive feedback between turbulence and the vertical shear of the mean flow. Turbulent sensible and latent heat flux anomaly were also analyzed. The results show that in both daytime and at night, when the surface layer turbulence was intensified as a result of strengthened mean flow, the sensible heat flux was decreased while the latent heat flux was increased. Both anomalous fluxes contributed to the cold island effect and the moisture island effect of the oasis.

Rainfall effect on wind waves and the turbulence beneath air-sea interface

Rainfall effects on wind waves and turbulence are investigated through the laboratory experiments in a large wind-wave tank. It is found that the wind waves are damped as a whole at low wind speeds, but are enhanced at high wind speeds. This dual effect of rain on the wind waves increases with the increase of rain rate, while the influence of rainfall-area length is not observable. At the low wind speed, the corresponding turbulence in terms of the turbulent kinetic energy （TKE） dissipation rate is significantly enhanced by rain- fall as the waves are damped severely. At the high wind speed, the augment of the TKE dissipation rate is suppressed while the wind waves are enhanced simultaneously. In the field, however, rainfall usually hin- ders the development of waves. In order to explain this contradiction of rainfall effect on waves, a possibility about energy transfer from turbulence to waves in case of the spectral peak of waves overlapping the inertial subrange of turbulence is assumed. It can be applied to interpret the damping phenomenon of gas trans- fer velocity in the laboratory experiments, and the variation of the TKE dissipation rates near sea surface compared with the law of wall.

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.

Convection:a neglected pathway for downward transfer of wind energy in the oceanic mixed layer

Upper-ocean turbulent mixing plays a vital role in mediating air-sea fluxes and determining mixed-layer properties, but its energy source, especially that near the base of the mixed layer, remains unclear. Here we report a potentially significant yet rarely discussed pathway to turbulent mixing in the convective mixed layer. During convection, as surface fluid drops rapidly in the form of convective plumes, intense turbulence kinetic energy(TKE) generated via surface processes such as wave breaking is advected downward, enhancing TKE and mixing through the layer. The related power, when integrated over the global ocean except near the surface where the direct effect of breaking waves dominates, is estimated at O(1)TW, comparable to that required by maintaining the Meridional Overturning Circulation(MOC). The mechanism in question therefore deserves greater research attention, especially in view of the potential significance of its proper representation in climate models.

Multi-physical Modeling and Adjusting for Ultrasonic Assisted Soft Abrasive Flow Processing

The polishing efficiency of the soft abrasive flow(SAF)method is low,which is not in line with the concept of carbon emission reduction in industrial production.To address the above issue,a two-phase fluid multi-physics modeling method for ultrasonic-assisted SAF processing is proposed.The acoustics-fluid coupling mechanic model based on the realizable k-ε model and Helmholtz equation is built to analyze the cavitation effect.The results show that the pro-posed modeling and solution method oriented to ultrasonic-assisted SAF processing have better revealed the flow field evolution mechanism.The turbulence kinetic energy at different ultrasonic frequencies and amplitudes is stud-ied.Simulation results show that the ultrasonic vibration can induce a cavitation effect in the constrained flow chan-nel and promote the turbulence intensity and uniformity of the abrasive flow.A set of comparative polishing experiments with or without ultrasonic vibration are conducted to explore the performance of the proposed method.It can be found that the ultrasonic-assisted SAF method can improve the machining efficiency and uniformity,to achieve the purpose of carbon emission reduction.The relevant result can offer a helpful reference for the SAF method.

Flow characteristics in curved open channel with partial suspended vegetation

Suspended vegetation in open channels such as natural river,lake,reservoir usually affect the flow structure,causing the change of the water environment,sediment transport,bed deformation.In order to study the water flow behavior in curved open channels under the influence of suspended vegetation,experiments were conducted in a Ushaped flume with 180°bend where the suspended vegetation substituted by cylindrical glass rods were partially placed.The particle image velocimeter(PIV)system was employed to measure the flow velocities of various cases with different vegetation arrangements and roots length.Comparison and analysis were conducted for measured data,such as water level,velocity,head loss,Reynolds stress,and turbulence kinetic energy(TKE)to obtain certain general rules of water flow in open curved channels with suspended vegetation.It can be found from the measured data that the water level and the head loss in the vegetation area are closely related to the arrangements of the suspended vegetation.Furthermore,the drag force of the vegetation can not only lead to smaller increments of longitudinal TKE above the vegetation tail than that below it,but also lead to the flow velocities in vegetation area much smaller than areas without vegetation.In addition,suspended vegetation weakens the vortex strength near the water surface and increases the vortex strength below the vegetation tail,and the flow velocity in the vegetation area has a close relationship with the suspended vegetation length in the water.Therefore,it can be concluded that the reasonable arrangements of suspended vegetation group in curved open channel can protect the concave bank from being scoured,and protect the convex bank from being deposited.

Reconstructing the upper ocean thermal profiles using one-dimensional numerical model

The observation data for 5 d at a station in the South China Sea is presented. After brief anMysis of the wind speed, air temperature from the ship-borne meteorological instruments and temperature and salinity profiles from the CTD （conductivity, temperature, depth recorder） data, the authors find that the CTD casts are too sparse for us to understand the diurnal evolution of the thermal structure in the upper ocean. A one-dimensional （1D） numericM code based on Mellor-Yamada turbulence closure model is used to reconstruct the upper ocean thermal structure, utilizing the atmospheric forcing data from ship-borne weather station. The simulation results show good agreement with the observational data; the significance of breaking waves is also briefly discussed. The evolution of turbulence kinetic energy （TKE） and the contribution from shear production and buoy- ancy production are discussed respectively. Finally, several possible factors which might influence the numerical results are briefly analyzed.

Simulation of the ocean surface mixed layer under the wave breaking

A one-dimensional mixed-layer model, including a Mellor- Yamada level 2.5 turbulence closure scheme, was implemented to investi- gate the dynamical and thermal structures of the ocean surface mixed layer in the northern South China Sea. The turbulent kinetic ener- gy released through wave breaking was incorporated into the model as a source of energy at the ocean surface, and the influence of the breaking waves on the mixed layer was studied. The numerical simulations show that the simulated SST is overestimated in summer without the breaking waves. However, the cooler SST is simulated when the effect of the breaking waves is considered, the corre- sponding discrepancy with the observed data decreases up to 20% and the MLD calculated averagely deepens 3.8 m. Owing to the wave-enhanced turbulence mixing in the summertime, the stratification at the bottom of the mixed layer was modified and the tempera- ture gradient spread throughout the whole thermocline compared with the concentrated distribution without wave breaking.

Structures of convection and turbulent kinetic energy in boundary layer over the southeastern edge of the Tibetan Plateau

Based on a comprehensive analysis on Sonic Anemometer and gradient data, wind profile radar(WPR) and GPS sounding data of March–August 2008 from the boundary layer(BL) tower observation system at Dali on the southeastern edge of Tibetan Plateau(TP), it is found that the strengths of turbulent kinetic energy(TKE), buoyancy term and shear term depend on vegetation cover in association with local stability and thermodynamic condition. Strong kinetic turbulence appears when near surface layer in neutral condition with the large contribution from shear term. In an unstable condition within near surface layer, the atmospheric turbulent motion is mainly thermal turbulence, as buoyancy term is obviously larger than shear term. Under a stable condition the intermittent turbulence is accompanied by weak shear and buoyancy term, and TKE is significantly less than neutral or instable condition. The study also presents that the buoyancy term contribution at Nyingchi station in the southern slopes of the TP large topography in spring is significantly larger than that at Dali over the southeastern TP edge, reflecting that the thermal turbulence makes an important contribution to convection activity in the southern slopes of TP. Dali station is located in complex terrain with mountain and valley leading to larger kinetic turbulence. From the perspective of interaction of turbulence-convection in different scales, the study revealed that the height of convective boundary layer(CBL) could reach up to 1500–2000 m. TKE, shear term, and buoyancy term in near surface layer have the notable correlations with BL height and local vertical motion. The daytime thermodynamic turbulence effect of heat flux and buoyancy term has an obvious impact on the height of CBL, whereas mechanical turbulence only exerts a less impact. Mechanical turbulence in near surface layer has a significant impact on vertical motion especially in the forenoon with impacting height of 2500–3000 m. The peaks in diurnal variations of shear term and buoyancy term correspond to the high instable periods, especially in summer forenoon. Our observation analysis characterized the convection activity triggered by TKE source and their interaction in the southeastern TP edge.

Large Eddy Simulation and Study of the Urban Boundary Layer

Based on a pseudo-spectral large eddy simulation (LES) model, an LES model with an anisotropy turbulent kinetic energy (TKE) closure model and an explicit multi-stage third-order Runge-Kutta scheme is established. The modeling and analysis show that the LES model can simulate the planetary boundary layer (PBL) with a uniform underlying surface under various stratifications very well. Then, similar to the description of a forest canopy, the drag term on momentum and the production term of TKE by subgrid city buildings are introduced into the LES equations to account for the area-averaged effect of the subgrid urban canopy elements and to simulate the meteorological fields of the urban boundary layer (UBL). Numerical experiments and comparison analysis show that: (1) the result from the LES of the UBL with a proposed formula for the drag coefficient is consistent and comparable with that from wind tunnel experiments and an urban subdomain scale model; (2) due to the effect of urban buildings, the wind velocity near the canopy is decreased, turbulence is intensified, TKE, variance, and momentum flux are increased, the momentum and heat flux at the top of the PBL are increased, and the development of the PBL is quickened; (3) the height of the roughness sublayer (RS) of the actual city buildings is the maximum building height (1.5-3 times the mean building height), and a constant flux layer (CFL) exists in the lower part of the UBL.

CFD modeling to study the effect of particle size on dispersion in 20l explosion chamber: An overview

Mine disasters occur predominantly due to methane or coal dust explosion or a combination of both.Among the top ten worst coal mine disasters in India, nine are due to coal dust explosion. The current paper describes a general overview of the parameters causing dispersion leading to coal dust explosion,and computational fluid dynamics(CFD) simulation study to observe the effects of particle size on dispersion in Indian coal mines. Turbulent kinetic energy(TKE) and velocity vector path of dust-air mixture and dust-free air were simulated to understand their effects on coal dust dispersion. The TKE contours and velocity vector paths for dust-free air were uniform and symmetrical due to resistance-free path available. The TKE contours and velocity vector paths for dust-air mixture shows the asymmetrical distribution of contours, due to entrainment of air with dust particles. Vortices were observed in velocity vector paths which gradually diminish on increment of time sequence. These vortices are dead centres where velocity and coal dust particles concentration are both zero.