Entire Boundary Layer Research Articles

Inconsistent 3‐D Structures and Sources of Sulfate Ammonium and Nitrate Ammonium Aerosols During Cold Front Episodes

AbstractSince the distinct thermostability difference of sulfate ammonium and nitrate ammonium aerosols, their distributions, evolutions and sources could be unpredictable on a long‐range transport condition. Here, we highlighted the 3‐D structures and sources of SO42−, NO3− and NH4+ (SNA) during two cold front episodes in east China. Cold fronts effectively uplift and transport PM2.5 and its precursors from upstream sources to the Yangtze River Delta (YRD). Specifically, in the YRD, surface SO42− is mostly imported from the upstreams, accounting for ∼48%, significantly higher than the contribution from the YRD itself (∼29%). In contrast, NH4NO3 is thermally unstable and more easily lost in the warmer and lower boundary layer (BL) ahead of cold front. Consequently, only 20% of the total NO3− originates from upstreams, while the YRD contributes 28%. In the upper BL, the contribution of SO42− from upstreams remain high (∼49%), with only 18% originating from the YRD. However, due to the intense thermostability of NH4NO3 in colder and wetter air, the YRD’s contribution to NO3− is 27%, and upstreams contribute 20%. The physical processes exert relatively consistent effects on variations of PM2.5 and SNA concentrations. The aerosol chemical process (AERO) of (NH4)2SO4 consistently contributes positively throughout the entire BL. Conversely, the temperature‐sensitive NH4NO3 undergoes repeated dissociation/condensation and deposition, causing positive AERO contributions in upper BL and negative contributions in lower BL. Results indicate that one difference in physicochemical property of species could induce their distinct distributions and sources in large scale, and should be considered in regional air pollution control.

Embedded shear layers in turbulent boundary layers of a NACA0012 airfoil at high angles of attack

An investigation of turbulent boundary layers (TBLs) is presented for a NACA0012 airfoil at angles of attack 9 and 12 deg. Wall-resolved large eddy simulations (LES) are conducted for a freestream Mach number M=0.2 and chord-based Reynolds number Re=4×105, where the boundary layers are tripped near the airfoil leading edge on the suction side. For the angles of attack analyzed, mild, moderate and strong adverse pressure gradients (APGs) develop over the airfoil. Despite the strong APGs, the mean flow remains attached along the entire airfoil suction side. Similarly to other APG-TBLs investigated in the literature, a secondary peak appears in the Reynolds stress and turbulence production profiles. This secondary peak arises in the outer layer and, for strong APGs, it may overcome the first peak typically observed in the inner layer. The analysis of the turbulence production shows that other components of the production tensor become important in the outer layer besides the shear term. For moderate and strong APGs, the mean velocity profiles depict three inflexion points, the third being unstable under inviscid stability criteria. In this context, an embedded shear layer develops along the outer region of the TBL leading to the formation of two-dimensional rollers typical of a Kelvin–Helmholtz instability which are captured by a spectral proper orthogonal decomposition (SPOD) analysis. The most energetic SPOD spatial modes of the tangential velocity show that streaks form along the airfoil suction side and, as the APG becomes stronger, they grow along the spanwise and wall-normal directions, having a spatial support along the entire boundary layer.

Sensitivity of boundary layer features to depth-dependent baroclinic pressure gradient and turbulent mixing in an ocean of finite depth

The present numerical study builds on Ekman (1905)’s work in surface boundary layer and extends the boundary value problem to overcome some of its limitations. Previous studies addressed model’s limitations by assuming that deviations from observations are usually ascribed to different eddy viscosity shapes, but seldom to the presence of baroclinic pressure gradients and shallow seas, which are the mainstays of this work. Improved solutions in the ocean boundary layer are obtained considering both depth-dependent wind-induced eddy viscosity and horizontal density gradients, ranging from well-mixed to highly-stratified conditions in a finite-depth ocean. High-order numerical solutions extend those in previous analytical and numerical works in the literature and widens the parameter space analyzed. Remarkably, the current profiles are obtained without ambiguity as a truly superposition of a geostrophic and a ageostrophic terms. Results indicate that, for a vertically-uniform eddy viscosity without density gradients and in shallow waters, currents are practically aligned with wind. As depth increases, misalignment between currents and wind increases and the complexity of the vertical structure increases. At large depths, Ekman’s values are attained, i.e., deflection angles relative to wind direction, θW, are θS−θW=−45∘ at the surface, where the current is maximum, and θT−θW=−90° for the depth-integrated transport (negative for deflections to the right in the Northern Hemisphere). These features remain regardless of the magnitude of the eddy-viscosity. For non-uniform eddy viscosity, θS−θW decreases from −45∘ up to −90∘ from low to high stratification level, respectively, whereas θT−θW is rather insensitive (θT−θW≈−90∘). Contrary to wind effects, the presence of the only vertically-uniform density gradient forcing, with constant eddy viscosity, deflects the surface angle θS−θD=+45∘ relative to the density-gradient direction, θD, in deep waters. Maximum currents no longer occur at the surface in this case. For non-uniform density-gradient profiles, current magnitude decreases overall while θS−θD≈+45∘ as long as the gradient affects the entire surface boundary layer. The deflection angle θT−θD∼+95∘ remains less sensitive to changes in density-gradient profile.

Vertical Evolution of Ozone Formation Sensitivity Based on Synchronous Vertical Observations of Ozone and Proxies for Its Precursors: Implications for Ozone Pollution Prevention Strategies.

Photochemical ozone (O3) formation in the atmospheric boundary layer occurs at both the surface and elevated altitudes. Therefore, the O3 formation sensitivity is needed to be evaluated at different altitudes before formulating an effective O3 pollution prevention and control strategy. Herein, we explore the vertical evolution of O3 formation sensitivity via synchronous observations of the vertical profiles of O3 and proxies for its precursors, formaldehyde (HCHO) and nitrogen dioxide (NO2), using multi-axis differential optical absorption spectroscopy (MAX-DOAS) in urban areas of the Beijing-Tianjin-Hebei (BTH), Yangtze River Delta (YRD), and Pearl River Delta (PRD) regions in China. The sensitivity thresholds indicated by the HCHO/NO2 ratio (FNR) varied with altitude. The VOC-limited regime dominated at the ground level, whereas the contribution of the NOx-limited regime increased with altitude, particularly on heavily polluted days. The NOx-limited and transition regimes played more important roles throughout the entire boundary layer than at the surface. The feasibility of extreme NOx reduction to mitigate the extent of the O3 pollution was evaluated using the FNR-O3 curve. Based on the surface sensitivity, the critical NOx reduction percentage for the transition from a VOC-limited to a NOx-limited regime is 45-72%, which will decrease to 27-61% when vertical evolution is considered. With the combined effects of clean air action and carbon neutrality, O3 pollution in the YRD and PRD regions will transition to the NOx-limited regime before 2030 and be mitigated with further NOx reduction.

Dynamical Importance of the Trade Wind Inversion in Suppressing the Southeast Pacific ITCZ

AbstractSea surface temperature (SST) gradients are a primary driver of low‐level wind convergence in the east Pacific Inter‐Tropical Convergence Zone (ITCZ) through their hydrostatic relationship to the surface pressure gradient force (PGF). However, the surface PGF may not always align with SST gradients due to variations in boundary layer temperature gradients with height, that is, the boundary layer contribution to the surface PGF. In this study, we investigate the observed northern hemisphere position of the east Pacific ITCZ using a slab boundary layer model (SBLM) driven by different approximations of the boundary layer virtual temperature field. SBLM simulations using the entire boundary layer virtual temperature profile produce a realistic northern hemisphere ITCZ. However, SST‐only simulations produce excessive equatorial divergence and southern hemisphere convergence, resulting in a latitudinally confined double ITCZ‐like structure. Observed virtual temperature gradients highlight the importance of northward temperature gradients strengthening with height from the equator to 15°S below the trade wind inversion (TWI). Our interpretation is that the equatorial cold tongue induces relatively weak high surface pressure and double ITCZ‐like convergence because the resulting layer of cold air is shallow. Concurrently, relatively strong high surface pressure spreads out in the southern hemisphere due to interactions between stratocumulus clouds and the ocean surface. Together, the equatorial cold tongue and the TWI/stratocumulus clouds enable a more northern hemisphere dominant ITCZ. Thus, we provide evidence of a dynamical link between the equatorial cold tongue, low clouds, and double ITCZs, which continue to be problematic in Earth system models.

Scale effect of micro ribs on the turbulent transport in an internal cooling channel

Owing to the limited supply and pressure margin in the air system, a cooling technique providing efficient heat transfer with lower flow loss is highly desirable for gas turbine blades. Microscale ribs have promised to be a potential cooling candidate. In this work, large eddy simulations are implemented to reveal the scale effect of micro ribs on the near-wall turbulent transport in a cooling channel. Considering a mechanistic study and practical applications, both single-rib and rib-array arrangements are studied with a wide range of dimensionless viscous-scaled rib heights involving the entire boundary layer. The results indicate that the rib-induced destruction and regeneration of coherent structures are, respectively, responsible for the weakened momentum transport and enhanced heat transport in the near-wall region. Using tiny ribs, regenerated quasi-streamwise vortices are mainly located in the buffer layer. The resulting turbulence burst greatly enhances wall heat transfer while keeping a lower flow loss due to the weak form drag. Regenerated hairpin vortices using tall ribs are activated in the log-law layer and intensively interact with mainstream. Along with improved wall heat transfer, the significant form drag results in a remarkably high flow loss. Accordingly, heat transfer and flow loss show different dependencies on the rib height, which contributes to an optimum height interval of ribs (e+ = 20–40) located in the high buffer and low log-law layer for maximizing the overall performance. Furthermore, for the rib-array scheme, adequate inter-rib spacing is essential to achieve turbulence regeneration for enhancing near-wall heat transport.

Experimental investigation of a grid plasma jet array in a turbulent boundary layer for skin-friction drag reduction

Toward turbulent skin-friction drag reduction, a novel layout of plasma actuator with a grid anode is devised (grid edge length: L), capable of producing an array of wall-normal plasma jets. The characteristics of this plasma jet array and its interaction with a turbulent boundary layer are investigated experimentally with a particle imaging velocimetry. Results show that the quiescent flow field of the plasma actuator is dominated by a standing vortex ring attached at the inner sides of the square grid, a wall-normal jet flow issued from the grid center, and a downwash flow between adjacent grids. When a tandem array of 11 plasma jets are issued into the cross flow, an equilibrium stage is reached after the third jet. In this stage, the main body of the wall-normal jet bends noticeably to the cross flow, and its leeward side hosts a reverse flow zone, extending downstream to form a slender low-speed wedge (LSW). Two vortical structures are prominent: the streamwise counter-rotating vortex pair residing in the two sides of the jet body, and the arch-shape negative spanwise vortex situated on the LSW root. In the grid-middle plane, the production of turbulence is enhanced across the entire boundary layer, and a second production peak is identified at y+=26. Reduction of the spanwise-averaged wall shear stress is achieved downstream of plasma actuation at x/L≥19.4, and part of the drag reduction fruit earned by ejection is offset by the spanwise transportation of high-velocity fluids toward the middle plane.

Investigating Reynolds number effects in turbulent concentric coaxial pipe flow using stochastic one‐dimensional turbulence modeling

AbstractThe present study numerically investigates turbulent momentum transfer in concentric coaxial (annular) pipe flow with small radius ratios . To model the flow, a stochastic one‐dimensional turbulence (ODT) model formulated for cylindrical geometry is used that provides full‐scale resolution along a representative radial coordinate. The present investigation extends the model validation by Tsai et al. (PAMM, 22:e202200272, 2023), to radius ratios smaller than 0.1 and addresses boundary layers with strong spanwise curvature effects. The focus is on the assessment and analysis of statistical flow features in the vicinity of the inner cylinder wall, particularly in cases with small radius ratios. Following Boersma & Breugem (Flow Turbul. Combust., 86:113–127, 2011), classical boundary‐layer and mixing‐length theory is utilized to analyze the model predictions. The results demonstrate that the ODT model captures leading‐order curvature and mixing‐length effects by its physics‐compatible construction. Utilizing the model for extrapolation to high Reynolds numbers inaccessible to conventional high‐fidelity numerical approaches shows that curvature effects persist and nonlocally affect the entire boundary layer. The model results provide support for a spanwise‐curvature‐modified wall function.

Effects of cooling roofs on mitigating the urban heat island and human thermal stress in the Pearl River Delta, China

The Weather Research and Forecasting model was employed to simulate the effects of cooling roofs on the urban heat island effect and human thermal stress in the Pearl River Delta region, China, in summer (June–August) 2014. It was found that cooling roofs can reduce the 2-m air and ground-surface temperature, wind speed, and planetary boundary layer height, while increasing the specific and relative humidity, in urban areas. The effects can influence the entire urban boundary layer, especially in commercial/industrial areas in the daytime. Cooling roofs can effectively mitigate the urban heat island effect, and the mitigation of white roofs is more significant than that of green roofs. The mitigation of the urban heat island effect is stronger in commercial/industrial areas than in low-intensity residential areas. The city-wide 2-m and ground-surface urban heat island intensity induced by white roofs can reduce by 1.5 °C and 4.8 °C at noon, while the corresponding decreases induced by green roofs can reach 0.7 °C and 1.9 °C, respectively. White roofs improve human thermal comfort both in the daytime and at nighttime. With white roofs, three city-wide thermal-stress indices (wet-bulb globe temperature, apparent temperature, and humidity index) can be reduced at noon by 0.8 °C, 1.3 °C and 1.5 °C, respectively. Green roofs mostly improve human thermal comfort at nighttime but even deteriorate it in the daytime. The above three indices decrease by less than 0.5 °C in different urban categories at nighttime with green roofs.

On the Composite Velocity Profile in Zero Pressure Gradient Turbulent Boundary Layer: Comparison with DNS Datasets

Data obtained by direct numerical simulations (DNS) of the Zero-Pressure-Gradient Turbulent Boundary Layer were analyzed and compared to a mathematical model of the mean velocity profile (MVP) in the range 1000 ≤ Reθ ≤ 6500. The mathematical model is based on the superposition of an accurate description of the inner law and Coles’ wake function with appropriately chosen parameters. It is found that there is excellent agreement between the mathematical model and the DNS data in the inner layer when the Reynolds number based on momentum thickness, Reθ, is greater than 1000. Furthermore, there is very good agreement over the entire boundary layer thickness, when Reθ is greater than 2000. The diagnostic functions Ξ and Γ based on DNS data are examined and their characteristics are discussed in relation to the existence of a logarithmic layer or a power law behavior of the MVP. The diagnostic functions predicted by the mathematical model are also presented.

Air–Sea Enthalpy and Momentum Exchange Coefficients from GPS Dropsonde Measurements in Hurricane Conditions

The intensity of tropical cyclones is highly dependent on air–sea enthalpy and momentum exchange. At extreme wind speeds, the values of the enthalpy, CK, and momentum, CD, exchange coefficients are characterized by high uncertainty. The present study aims to expand the previously used algorithm for CD retrieval to obtain the values of CK from wind speed measurements and the enthalpy profiles obtained from NOAA GPS dropsondes in hurricane conditions. This algorithm uses concepts from technical hydrodynamics, describing turbulent boundary layers on flat plates and pipes. According to this approach, the velocity (and enthalpy) defect profiles are self-similar in the entire boundary layer, including the layer of constant fluxes and the “wake” part, where the airflow adapts to the undisturbed flow region. By using the self-similarity property, the parameters of the constant flow layer (the roughness parameter, friction velocity, and the enthalpy and exchange coefficients CD and CK) could be obtained from measurements in the “wake” part for wind speeds from 20 m/s to 72 m/s. The estimates of the CK/CD ratio revealed values of 0.7 and 0.96 (depending on the self-similar approximation limits), and the results suggest that there are slight variations with the wind speed.

Grand challenges in the design, manufacture, and operation of future wind turbine systems

Abstract. Wind energy is foundational for achieving 100 % renewable electricity production, and significant innovation is required as the grid expands and accommodates hybrid plant systems, energy-intensive products such as fuels, and a transitioning transportation sector. The sizable investments required for wind power plant development and integration make the financial and operational risks of change very high in all applications but especially offshore. Dependence on a high level of modeling and simulation accuracy to mitigate risk and ensure operational performance is essential. Therefore, the modeling chain from the large-scale inflow down to the material microstructure, and all the steps in between, needs to predict how the wind turbine system will respond and perform to allow innovative solutions to enter commercial application. Critical unknowns in the design, manufacturing, and operability of future turbine and plant systems are articulated, and recommendations for research action are laid out. This article focuses on the many unknowns that affect the ability to push the frontiers in the design of turbine and plant systems. Modern turbine rotors operate through the entire atmospheric boundary layer, outside the bounds of historic design assumptions, which requires reassessing design processes and approaches. Traditional aerodynamics and aeroelastic modeling approaches are pressing against the limits of applicability for the size and flexibility of future architectures and flow physics fundamentals. Offshore wind turbines have additional motion and hydrodynamic load drivers that are formidable modeling challenges. Uncertainty in turbine wakes complicates structural loading and energy production estimates, both around a single plant and for downstream plants, which requires innovation in plant operations and flow control to achieve full energy capture and load alleviation potential. Opportunities in co-design can bring controls upstream into design optimization if captured in design-level models of the physical phenomena. It is a research challenge to integrate improved materials into the manufacture of ever-larger components while maintaining quality and reducing cost. High-performance computing used in high-fidelity, physics-resolving simulations offer opportunities to improve design tools through artificial intelligence and machine learning, but even the high-fidelity tools are yet to be fully validated. Finally, key actions needed to continue the progress of wind energy technology toward even lower cost and greater functionality are recommended.

Current shear and turbulence during a near-inertial wave

Surface currents and turbulent mixing were observed during a near-inertial wave (NIW) using an accousting doppler current profiler (ADCP) and satellite-tracked drifters. Drifter trajectories sampled at three depth levels show characteristics of an Ekman solution superposed with the NIW. Velocity and dissipation estimates from the ADCP reveal strong shear with a distinct constant flux layer in between the roughness length and a critical depth at 4m. Below, a shear free slab layer performing an inertial oscillation is observed. Dissipation, as estimated from the vertical beam of the ADCP, peaks in the wave-enhanced friction layer when the current opposes the wind and wave direction. Below the constant flux layer, maximum turbulence is observed when the NIW is in a phase that is in opposite direction to the time-averaged current. During this phase, currents at various depths rapidly realign in the entire boundary layer.

Thermosolutal Marangoni convection on a permeable flat surface with aligned magnetic field and thermophoretic effect: an OHAM solution

In this venture, the challenge to study the mathematical model and its approximate solution for boundary layer thermosolutal Marangoni convection on a permeable flat surface with thermophoretic velocity and aligned magnetic field has been undertaken. Using a special type of transformation, governing PDEs are converted into ODEs. A semi-analytical method known as optimal homotopy analysis method(OHAM) has been used to develop solutions for proposed mathematical model. The novelty of the study is in the exploration of simultaneous impacts of thermophoretic velocity and inclined magnetic field on flow and heat-mass transport while mass suction/injection is using OHAM. Solutions obtained by OHAM are compared with other published literature and it is revealed that our numerical values are improved by approximately up to 7%. The obtained solutions for various important parameters are displayed in the form of graphs and tables. The study explores the variable behaviors of velocity, temperature, and concentration for growths of Marangoni number, magnetic parameter, and thermophoretic parameter, while there exists mass suction or injection. It reveals that in both cases, i.e. suction and injection cases, velocity grows and temperature and concentration reduce with increasing magnetic parameters and inclination angle of aligned magnetic field and these variations are more prominent in case injection. If the magnetic field is orthogonal to the flat surface, then corresponding Lorentz force produces maximum resistance for the transport phenomenon. The spreading effect of induced surface velocity due to Marangoni convection is observed in the entire boundary layer region in case of mass suction, while for injection case it is witnessed only near the surface. Whereas, for higher value of thermophoretic velocity, the concentration exhibits a reducing trend. For rising values of Marangoni number, wall velocity, Nusselt number, and Sherwood number significantly grow, while magnetic field has contrary impacts.

Dynamic Mechanisms Associated with the Structure and Evolution of Roll Vortices and Coherent Turbulence in the Hurricane Boundary Layer: A Large Eddy Simulation During the Landfall of Hurricane Harvey

Roll vortices are a series of large-scale turbulent eddies that nearly align with the mean wind direction and prevail in the hurricane boundary layer. In this study, the one-way nested WRF-LES model simulation results from Li et al. (J Atmos Sci 78(6):1847–1867, https://doi.org/10.1175/JAS-D-20-0270.1, 2021) are used to examine the structure and generation mechanism of roll vortices and associated coherent turbulence in the hurricane boundary layer during the landfall of Hurricane Harvey from 00 UTC 25 to 18 UTC 27 August 2017. Results indicate that roll vortices prevail in the hurricane boundary layer. The intense roll vortices and associated large turbulent eddies above them (at a height of ~ 200 to 3000 m) accumulate within a hurricane radius of 20–40 km. Their intensity is proportional to hurricane intensity during the simulation period. Before and during hurricane landfall, strong inflow convergence leads to horizontal advection of roll vortices throughout the entire hurricane boundary layer. Combined with the strong wind shear, the strongest roll vortices and associated large turbulent eddies are generated near the eyewall with suitable thermodynamic (Richardson number at around − 0.2 to 0.2) and dynamic conditions (strong negative inflow wind shear). After landfall, the decayed inflow weakens the inflow convergence and quickly reduces the strong roll vortices and associated large turbulent eddies. Diagnosis of vertical turbulent kinetic energy indicates that atmospheric pressure perturbation, caused by horizontal convergence, transfers the horizontal component of turbulence to the vertical component with a mean wavelength of about 1 km. The buoyancy term is weak and negative, and the large turbulent eddies are suppressed.

Large-Scale Coherent Turbulence Structures in the Atmospheric Boundary Layer over Flat Terrain

Abstract We investigate characteristics of large-scale coherent motions in the atmospheric boundary layer using field measurements made with two long-range scanning wind lidars. The joint scans provide quasi-instantaneous wind fields over a domain of ∼50 km2, at two heights above flat but partially forested terrain. Along with the two-dimensional wind fields, two-point statistics and spectra are used to identify and characterize the scales, shape, and anisotropy of coherent structures—as well as their influence on wind field homogeneity. For moderate to high wind speeds in near-neutral conditions, most of the observed structures correspond to narrow streaks of low streamwise momentum near the surface, extending several hundred meters in the streamwise direction; these are associated with positive vertical velocity ejections. For unstable conditions and moderate winds, these structures become large-scale rolls, with longitudinal extent exceeding the measuring domain (>∼5 km); they dominate the conventional surface-layer structures in terms of both physical scale and relative size of velocity-component variances, appearing as quasi-two-dimensional structures throughout the entire boundary layer. The observations shown here are consistent with numerical simulations of atmospheric flows, field observations, and laboratory experiments under similar conditions. Significance Statement Coherent structures have attracted the interest of researchers for decades, being viewed as the closest to “order” that we can find within the chaos of turbulence. In the turbulent atmospheric boundary layer, micro- and mesoscale coherent structures come in many shapes and sizes, such as convective cells, rolls, or streaks. In this study we used dual lidars (remote sensing measurements), developing analysis of their tandem usage to characterize in detail some of the large-scale coherent structures generated over flat terrain. This allowed us to better understand the mechanisms that generate such structures and describe their influence on the morphology of the turbulent atmospheric boundary layer across a good deal of its depth.

Effects of thermophoresis on high-pressure binary-species boundary layers with uniform and non-uniform compositions

Direct numerical simulations of high-pressure binary-species temporal boundary layers are performed to investigate the flow physics for three situations: (1) uniform and equal composition, (2) uniform but unequal compositions and (3) non-uniform composition. Both colder- and hotter-wall situations compared with the free stream are simulated. The working fluid is a nitrogen/methane mixture. The analysis is performed at a case-specific self-similar state. Even when the initial composition is uniform, the methane mean mass fraction decreases near the colder wall, whereas it increases near the hotter wall and the mass fraction fluctuates in the entire boundary layer. Analysis of the species-mass diffusion balance and flow structures reveal that both mass-fraction variation and fluctuations are induced by the Soret effect. When the initial composition is non-uniform and the wall is colder, the methane mean mass fraction monotonically increases from the wall akin to its initial profile. However, when the wall is hotter the mean mass fraction decreases near the wall in contrast to its initial profile, a fact traced through the species-mass diffusion balance to the Soret effect being large and enriching methane near the wall. In contrast, the direction of the Soret flux is opposite for the colder wall, thus keeping the methane concentration small. Although the initial magnitude of the difference between the wall and free-stream temperature is the same in all cases, the situation is not symmetric between colder- and hotter-wall cases; the flow structure exhibits much smaller scales when the wall is hotter than when the wall is colder.

Modelling the downstream development of a turbulent boundary layer following a step change of roughness

In this study, we develop an analytical model to predict the turbulent boundary layer downstream of a step-change in the surface roughness where upstream flow conditions are given. We first revisit the classical model of Elliott (Trans. Am. Geophys. Union, vol. 39, 1958, pp. 1048–1054), who modelled the velocity distribution within and above the internal layer with a simple piecewise logarithmic profile, and evolved the velocity profile using the streamwise momentum equation. Elliott's model was originally developed for an atmospheric surface layer, and to make the model applicable to a spatially developing turbulent boundary layer with finite thickness, we propose a number of more physical refinements, including adding a wake function to the velocity profile, considering the growth of the entire boundary layer in the streamwise direction, and using a more realistic shear stress profile in the momentum equation. In particular, we implement the blending model (Li et al., J. Fluid Mech., vol. 923, 2021, p. A18) to account for the deviation of the mean flow within the internal layer from a canonical velocity profile based on the local wall condition. These refinements lead to improved agreement between the prediction and the measurement, especially in the vicinity of the rough-to-smooth change.

Applications of wide-ranging PIV measurements for various turbulent statistics in artificial atmospheric turbulent flow in a wind tunnel

Particle image velocimetry (PIV) is currently the only method which can experimentally quantify the spatial distributions of velocity with high sampling frequencies without causing intrusions into the airflow of the measuring devices. To extend the use of PIV systems for determining the distributions of velocity for wind environmental and engineering purposes, this study showed various applications of wide-ranging PIV systems for an artificially created turbulent flow at the wind-tunnel scale. Various turbulent statistics, which are commonly employed for evaluating the wind structures and environments, were compared among three conditions of the measurement areas covered by PIV systems and a hot-wire anemometry (HWA). The largest measurement covering 500 × 500 mm2 can capture the turbulent motions extending to 40% of the entire boundary layer. In addition, the discrepancies on two-length scales could be discussed using spatial and temporal correlation functions owing to the wide-ranging PIV. Moreover, the limitations of the wide-ranging PIV systems were clarified. While the turbulent statistics compared among three PIV conditions and HWA demonstrated adequate agreement, the PIV with the largest measurement area failed to capture the small-scale turbulent structures due to the coarse grid resolutions. Moreover, the measurement resolutions of PIV systems demonstrated sensitivity in determining the spectrum and probability density functions, whereas a wide-ranging PIV system can cover wide-ranging measurement areas.

Estimation of turbulence parameters in the atmospheric boundary layer of the Bohai Sea, China, by coherent Doppler lidar and mesoscale model.

Obtaining turbulence parameters in the marine atmospheric boundary layer (MABL) is limited by the observation environment and cost. Therefore, estimating based on the weather forecast model or combining the model output with limited observations is a more flexible choice. We conducted cruise observation experiments in the Bohai Sea, China, from May 17 to June 4, 2021. On the basis of the wind profile observed by the coherent Doppler lidar and the temperature, as well as pressure profiles output by the Weather Research and Forecasting (WRF) model, we implemented the Tatarskii turbulence model to estimate the refractive index structure constant C n2 in the atmospheric boundary layer of the Bohai Sea under clear sky. The temporal and spatial variations of turbulence in the Bohai Sea atmospheric boundary layer are studied by combining the vertical velocity variance σ w2, skewness Ske and kurtosis Kur. The performance of simulated C n2 and meteorological parameters in the WRF in the atmospheric boundary layer at the Bohai Sea is evaluated through the experimental measurements of UAV-borne (unmanned aerial vehicle) radiosonde and lidar. Finally, we give the model of the C n2 variation with height in the atmospheric boundary layer at the Bohai Sea. The results show that WRF can better simulate C n2 in most cases. The bias between the measured and simulated C n2 is within one order of magnitude, and the root mean square error ( RMSE ) is within two orders of magnitude. Due to the potential uncertainty of the WRF, the RMSE between the measured and simulated wind speed is 4 ms-1 to 6 ms-1, which is almost two times of the result in previous studies on the underlying land surface. The overall changes of C n2 and σ w2 are similar when the turbulence is well mixed and developed, which shows the consistency in both of optical and dynamics turbulence. But this consistency is not absolute. The temperature difference between the sea surface and the atmosphere leads to the widespread existence of an inversion layer from the sea surface to hundreds of meters in the Bohai Sea. The suppression of the inversion layer weakens the near sea surface turbulence. There is an enhancement of turbulence intensity below the inversion layer and a decrease from the upper inversion layer to top of the boundary layer among the entire boundary layer, also, the position of the inflection point is determined by the height of top of the inversion layer. The main results of this study are the reference significance for further understanding the development and change characteristics of turbulence in the MABL.