Search

[Display omitted] The friction and interfacial nanostructure of a water-in-surface-active ionic liquid mixture, 1.6 M 1-butyl-3-methylimidazolium 1,4-bis-2-ethylhexylsulfosuccinate ([BMIm][AOT]), can be tuned by applying potential on Au(1 1 1) and stainless steel. Atomic force microscopy (AFM) was used to examine the friction and interfacial nanostructure of 1.6 M [BMIm][AOT] on Au(1 1 1) and stainless steel at different potentials. Superlubricity (vanishing friction) is observed for both surfaces at OCP+1.0 V up to a surface-dependent critical normal force due to [AOT]− bilayers adsorbing strongly to the positively charged surface thus allowing AFM tip to slide over solution-facing hydrated anion charged groups. High-resolution AFM imaging reveals ripple-like features within near-surface layers, with the smallest amplitudes at OCP+1 V, indicating the highest structural stability and resistance to thermal fluctuations due to highly ordered boundary [AOT]− bilayers templating robust near-surface layers. Exceeding the critical normal force at OCP+1.0 V causes the AFM tip to penetrate the hydrated [AOT]− layer and slide over alkyl chains, increasing friction. At OCP and OCP-1.0 V, higher friction correlates with more pronounced ripples, attributed to the rougher templating [BMIm]+ boundary layer. Kinetic experiments show that switching from OCP-1.0 V to OCP+1.0 V achieves superlubricity within 15 s, enabling real-time friction control. [ABSTRACT FROM AUTHOR]

The turbulent boundary layer in a viscous incompressible fluid developing longitudinally past the surface of a thin cone or cylinder at a finite distance from the laminar–turbulent transition zone is studied. The characteristic Reynolds number, determined from the external flow velocity and the length of the body, is assumed to be large, and the thickness of the boundary layer is small and comparable to the radius of the body. The asymptotic method of multiple scales is used to find solutions to the Navier–Stokes equations. Instead of the traditional decomposition of the solution into time-averaged values and their fluctuations, the velocities and pressure are expressed as an asymptotic series consisting of steady and perturbed terms. As a result, the viscous steady flow ('secondary') that arises in the boundary layer as a mandatory component of fast turbulent fluctuations was described. Analytical and numerical solutions for the radial steady velocity are presented, describing the self-induced suction of fluid from the external flow into the boundary layer. Further analytical solutions are obtained for the longitudinal and circumferential velocities, which differ markedly from the laminar regime. The solutions found are somewhat similar to the degenerate (one-dimensional) case of self-sustaining longitudinal thin structures in turbulent shear flows. A qualitative comparison with direct numerical simulations is presented. [ABSTRACT FROM AUTHOR]

Dispersion within idealized urban environments was studied in a simulated neutrally buoyant, 1:200 scale boundary layer with the Meteorological Wind Tunnel at the EPA’s Fluid Modeling Facility. The measurements are used to offer a baseline of performance for the mechanical turbulence formulation and concentration predictions of AERMOD, the EPA’s preferred Gaussian dispersion model. Scaled meteorological conditions and dispersion characteristics were studied for both uniform and nonuniform building arrays oriented at 0° and 30° with respect to the flow and were compared to baseline, “rural”, measurements without the presence of buildings. Particle image velocimetry (PIV) measured velocity and shear stress profiles within each model configuration, whereas hydrocarbon analyzers (HCAs) measured ethane concentrations at defined points throughout the model. Four source locations were examined for each building array, with two in the urban core and two in a street canyon, each with a source within and above the building canopy. Experimental profiles, regardless of their shape, were fitted to Gaussian profiles to determine lateral and vertical plume spread and shift from the wind tunnel centerline. These parameters were compared against a no-building reference case. Concentration predictions using the formulations in AERMOD are computed for 3 variations of modeled velocity profiles for each source, using factor of 2 (FAC2) and fractional bias (FB) as the governing model evaluation parameters. The two urban configurations were found to decrease the FAC2 performance by 34.1% and 30.1% from the no-building reference for the uniform and nonuniform cases, respectively, while producing modeled concentrations of only 48.1% and 62.4% of the 10 highest observed concentrations. These results encouraged simple first-order corrections to improve model performance with an emphasis on predicting maximum concentrations for regulatory purposes. These corrections proved successful for the uniform cases, mitigating FB, and improving the FAC2 percentage by 11.4% with more mixed results in nonuniform configurations, highlighting the difficulty in applying uniformly derived parameterizations in realistic, nonuniform environments. [ABSTRACT FROM AUTHOR]

This study addresses time-dependent multiple-scale reaction-convection-diffusion initial boundary value systems characterized by strong coupling in the reaction matrix and weak coupling in the convection terms for a locally optimal accurate solution. The discrete problem, which typically loses its tridiagonal structure, expands its bandwidth to four in such coupled systems, resulting in a substantial computational load. Our objective is to mitigate this computational burden through a splitting approach that transforms the non-tridiagonal matrix into a tridiagonal form while maintaining the consistency, local optimal accuracy in space, and stability of the numerical scheme. We employ equidistributed non-uniform grids, guided by a carefully chosen monitor function, to approximate the continuous space domain. The discretization strategy targets local optimal linear accuracy across space and time on the domain's interior points. In addition, we have also provided the global convergence analysis of the present splitting approach, mathematically. The mathematical evidence is also obtained from the numerical experiments by comparing the splitting approach (either diagonal or triangular forms) of the reaction matrix to its coupled form. The results strongly confirm the effectiveness of this approach in delivering uniform linear accuracy, based on the present problem discretizations while significantly reducing the computational costs. [ABSTRACT FROM AUTHOR]

This paper numerically solves Fredholm integro-differential equations with small parameters and integral boundary conditions. The solution of these equations has a boundary layer at the right boundary. A central difference scheme approximates the second-order derivative, a backward difference (upwind scheme) approximates the first-order derivative, and the trapezoidal rule is used for the integral term with a Shishkin mesh. It is shown that theoretically, the proposed scheme is uniformly convergent with almost first-order convergence. Further to improve the order of convergence from first order to second order, we use the post-processing and the hybrid scheme. Two numerical examples are computed to support the theoretical results. [ABSTRACT FROM AUTHOR]

As Arctic sea ice rapidly melts, trans‐Arctic shipping routes are emerging with significant economic value and potentially reshaping global shipping patterns. The safety of ships navigating along these routes will be severely threatened by low visibility caused by Arctic sea fog. Here, we utilize ship and coastal site observations spanning 1979–2023 to investigate the climatological features and mechanisms of Arctic sea fog from June to September. We find that the Northern Sea Route (NSR, along the Russian coast) experiences a mean fog frequency of 20.0%, significantly higher than the 11.5% fog frequency observed along the Northwest Passage (NWP, along the Canadian coast). Fog frequency shows null correlation with sea ice but a strong negative correlation (−0.80) with surface air temperature. The lower mean surface air temperature in the NSR region enhances relative humidity and atmospheric stability, which are conducive to fog formation and maintenance. We further reveal that fog formation mechanisms differ between the NSR and NWP. Along the NSR, fog is usually advection fog, associated with upper boundary layer warming and moistening; in contrast, along the NWP, fog is typically radiation fog linked to near‐surface cooling. Our study highlights the distinct fog frequency and dominant mechanisms along the different trans‐Arctic shipping routes, thereby establishing observed reference for the model validation of simulating Arctic sea fog. Plain Language Summary: The rapid retreat of Arctic sea ice, driven by climate change, has led to a gradual increase in the area of open water in the Arctic. Arctic shipping routes will soon be navigable in summer, allowing for a nearly one‐third reduction in distance from the Far East to western Europe, compared with traditional routes via the Suez Canal. In recent years, it has gradually been realized that the frequent occurrence of sea fog in the Arctic summer poses a threat to the safety and speed of Arctic shipping routes. Nevertheless, due to the sparsity of Arctic observations, the understanding of Arctic sea fog remains relatively preliminary. This paper uses multisource observation data to comprehensively investigate the climatological characteristics and formation mechanisms of Arctic sea fog, particularly highlighting their differences along the shipping routes. This work deepens our understanding of Arctic sea fog and provides a scientific observation guidance for enhancing the safety of Arctic navigation. Key Points: Fog frequency along the Northern Sea Route (NSR) is twice of that along the Northwest Passage, dominated by air temperature rather than sea iceFogs along the NSR are associated with warming and moistening at upper boundary layer, which are mostly advection fogFogs along the Northwest Passage are linked to near‐surface cooling, which are mostly radiation fog [ABSTRACT FROM AUTHOR]

Based on the capillary bundle model, combined with the relationship between the boundary layer thickness and the pressure gradient, a modified low-velocity nonlinear flow model for offshore low permeability reservoir was established. By using the fluid-solid coupling model and the 3D fractured well black oil model based on embedded mesh, the stability of numerical simulation was compared with that of traditional and our proposed low-velocity nonlinear flow model. Finally, the well spacing of five points well pattern in offshore low permeability reservoir was optimized. The results show that compared with previous studies, our proposed nonlinear flow model is more reasonable. The accuracy of fitting with experimental data is higher. Compared with the traditional flow model, our proposed low-velocity nonlinear flow model can reduce the number of nonlinear iterations by 85. 7% in a numerical simulation model of depletion production. At the same time, in an injection-production numerical model of two fracturing wells, our proposed nonlinear flow model makes it stable convergent from non-convergent. With the increase of well spacing, the cumulative oil production of well group increases first and then decreases. There is an optimal well spacing. In addition, the optimal well spacing increases with increasing matrix permeability. This method provides an efficient tool for accurately optimizing well pattern parameters of offshore low permeability reservoir. [ABSTRACT FROM AUTHOR]

This paper is devoted to the study of the vanishing shear viscosity limit and strong boundary layer problem for the compressible, viscous, and heat-conducting planar MHD equations. The main aim is to obtain a sharp convergence rate which is usually connected to the boundary layer thickness. However, The convergence rate would be possibly slowed down due to the presence of the strong boundary layer effect and the interactions among the magnetic field, temperature, and fluids through not only the velocity equations but also the strongly nonlinear terms in the temperature equation. Our main strategy is to construct some new functions via asymptotic matching method which can cancel some quantities decaying in a lower speed. It leads to a sharp L ∞ convergence rate as the shear viscosity vanishes for global-in-time solution with arbitrarily large initial data. [ABSTRACT FROM AUTHOR]

Attempts to mitigate the computational cost of fully resolved large-eddy simulation (LES) in the near-wall region include both the hybrid Reynolds-averaged Navier–Stokes/LES (HRL) and wall-modeled LES (WMLES) approaches. This paper presents an LES wall treatment method that combines key attributes of the two, in which the boundary layer mesh is sized in the streamwise and spanwise directions comparable to WMLES, and the wall-normal mesh is comparable to a RANS simulation without wall functions. A mixing length model is used to prescribe an eddy viscosity in the near-wall region, with the mixing length scale limited based on local mesh size. The RANS and LES regions are smoothly blended using the dynamic hybrid RANS-LES (DHRL) framework. The results are presented for the turbulent channel flow at two Reynolds numbers, and comparison to the DNS results shows that the mean and fluctuating quantities are reasonably well predicted with no apparent log-layer mismatch. A detailed near-wall meshing strategy for the proposed method is presented, and estimates indicate that it can be implemented with approximately twice the number of grid points as traditional WMLES, while avoiding the difficulties associated with analytical or numerical wall functions and modified wall boundary conditions. [ABSTRACT FROM AUTHOR]

In this study, five different time-dependent incompressible non-Newtonian boundary layer models in two dimensions are investigated with the self-similar Ansatz, including external magnetic field effects. The power-law, the Casson fluid, the Oldroyd-B model, the Walter fluid B model, and the Williamson fluid are analyzed. For the first two models, analytical results are given for the velocity and pressure distributions, which can be expressed by different types of hypergeometric functions. Depending on the parameters involved in the analytical solutions of the nonlinear ordinary differential equation obtained by the similarity transformation, a vast range of solution types is presented. It turned out that the last three models lack self-similar symmetry; therefore, no analytic solutions can be derived. [ABSTRACT FROM AUTHOR]

The influence of real gas effects and a turbulent boundary layer on shock wave attenuation in the expansion tube is studied by numerically solving the axisymmetric compressible Navier–Stokes equations with an adaptive mesh refinement technique. Numerical simulation results reveal that the ideal gas assumption is not applicable to the expansion tube, and the turbulent boundary layer plays a major role in decreasing the shock wave speed in the acceleration tube of the expansion tube. Shock wave attenuation is attributed to the turbulent boundary layer decreasing the pressure behind the shock wave. The numerical simulations that include the real gas effects and the development of turbulent boundary layers qualitatively agree with analytical solutions in the shock tube, and they show good agreement with the experimental results, especially for the shock speed in the acceleration tube of the expansion tube. Both effects should be considered in the numerical simulation model aimed to support experiments in expansion tubes. [ABSTRACT FROM AUTHOR]

This study applies dynamic mode decomposition (DMD) to three‐dimensional simulation results of urban heat island circulation (UHIC, which is horizontal circulation) and thermals (vertical convections). The aim of this study is to revisit how these phenomena coexist based on the characteristics of temporal changes in the flow field. We used DMD to obtain the dominant spatial patterns and information on temporal changes. One of the modes of horizontal wind, which does not change temporally (no oscillation or amplification), exhibits a spatial UHIC pattern. The unique feature of this UHIC mode is that there are small‐scale striated structures (150–200 m) and large‐scale convergence. The other modes are time‐varying (oscillating and decaying) and represent smaller spatial‐scale phenomena (150–250 m), such as thermals. The frequency of each mode takes various values, some of which are lower than the lifetime of thermals in accordance with the Deardorff convective scale (~10 min). These low‐frequency modes showed striated structures similar to that observed in the UHIC modes. These results suggest that UHIC and thermals deform each other through components that vary in long temporal scales. [ABSTRACT FROM AUTHOR]

Low-level moisture is an important ingredient for forecasting severe storms, especially over the Great Plains where severe storms often develop along the dryline. Although ground-based observation systems such as lidar or radiosondes on weather balloons provide accurate information on low-level moisture, data are provided at limited locations, and the low temporal resolution of radiosondes makes it difficult to track a rapidly developing dryline. Geostationary satellites provide high spatial and temporal observation, but the channels of current geostationary satellites are mostly sensitive to water vapor at mid- to upper levels. However, the split window difference (SWD) between the "clean" window channel (10.3 μm) and the "dirty" window channel (12.3 μm) is commonly used for estimating low-level water vapor. However, this estimation is complicated by surface temperature contributions, dependence on the lapse rate, and nonlinear relationships between SWD and moisture. This study applies machine learning techniques to infer boundary layer precipitable water (BLPW) from Geostationary Operational Environmental Satellite (GOES) Advanced Baseline Imager (ABI) data. Since there are few observations that cover wide regions for training convolutional neural networks, especially for the atmosphere above the surface, High-Resolution Rapid Refresh (HRRR) model outputs are used as the truth for training. Statistical results show that using ABI channel 13, 14, and 15 brightness temperatures, skin temperature, solar zenith angle, and GOES total precipitable water product as additional inputs shows the lowest mean-square error (MSE). Case study results show that the model developed in this study is good at capturing strong moisture gradients and depicting a dryline. Significance Statement: Knowledge of low-level moisture is important for monitoring drylines and moisture convergence preceding convective initiation. Traditional observations lack the spatial resolution, spatial coverage, or temporal update frequency needed by forecasters. The latest generation of Geostationary Operational Environmental Satellite (GOES), the GOES-R series, provides high spatial resolution and rapid temporal updates. However, the estimation of low-level moisture using the Advanced Baseline Imager (ABI) is complicated by dependencies on other variables and nonlinear relationships. Machine learning is shown to provide good estimates of low-level moisture from ABI observations. This application of machine learning to GOES-R can improve the ability to monitor low-level moisture and forecast convective initiation. [ABSTRACT FROM AUTHOR]

This study investigates the various uses of density of motile microorganisms in the context of the flow of a binary base fluid with double diffusion past a vertical surface. The research aims to comprehend the interactions between motile microorganisms and the fluid dynamics, as well as the heat and mass transport mechanisms in this system. The analysis involves mathematically constructing the governing equations, transforming them into dimensionless nonlinear ordinary differential equations using similarity transformations, and numerically solving them using the MATLAB bvp4c solver. An analysis of the influence of several parameters on the profiles of velocity, temperature, concentration, nanoparticle concentration, and density of motile microorganisms is conducted using graphical representation. The findings demonstrate that boosting the thermophoresis parameter intensifies the temperature profile. In addition, an increase in the nanofluid Schmidt number results in a larger concentration of nanoparticles, whereas a higher bioconvection Lewis number reduces the density of the motile microorganism profile. These findings may find use in biomedical engineering as well as industrial processes that include enhancing the efficiency of mass transfer and bioconvection. Numeric simulation prophesies 99.9% for both shear stress and heat transfer rate intensification for Prandtl values are noticed. [ABSTRACT FROM AUTHOR]

The heat transfer in power‐law fluids across three corrugated circular cylinders placed in a triangular pitch arrangement is studied computationally in a confined channel. Continuity, momentum, and energy balance equations were solved using ANSYS FLUENT (Version 18.0). The flow is assumed to be steady, incompressible, two‐dimensional, and laminar. A square domain of side 300Dh is selected after a detailed domain study. An optimized grid with 98,187 cells is used in the study. The convergence criteria of 10−7 for the continuity, x‐momentum, and y‐momentum balances and 10−12 for the energy equation were used. Constant density and non‐Newtonian power‐law viscosity modules were used. The diffusive term is discretized using a central difference scheme. Convective terms are discretized using the Second‐Order Upwind scheme. Pressure–velocity coupling between continuity and momentum equations was implemented using the semi‐implicit method for pressure‐linked equation scheme. Streamlines show wake development behind the cylinders, which is very dominant at large ReN and n. Isotherm contours are cramped at higher values of ReN and PrN, implying higher heat transfer. Global parameters, like, Cd and Nu, are computed for the wide ranges of controlling dimensionless parameters, such as power‐law index (0.3 ≤ n ≤ 1.5), Reynolds (0.1 ≤ ReN ≤ 40), and Prandtl (0.72 ≤ PrN ≤ 500) numbers. The NuLocal plot attains a pitch near the corrugation of the surface due to abrupt changes in velocity and temperature gradients. Nu increases with ReN and/or PrN and decreases with n under ot herwise identical situations. Nu is correlated with pertinent parameters, namely, ReN, PrN, and n. [ABSTRACT FROM AUTHOR]

To improve the aerodynamic efficiency of the aircraft, the article focuses on optimizing the L/D ratio through postponement of flow separation. Elevating the L/D ratio leads to diminished drag and delays stall occurrence at high angles of attack (α). The wing geometry adopts NACA 6-digit airfoils, specifically NACA 63418 and NACA 63415. Various aerodynamic devices are explored for their ability to defer flow separation, with passive aerodynamic devices being the prevalent choice. The primary goal of the research is to integrate rotational thin wire elements into the aerodynamic components of the wing, aiming to diminish drag, amplify lift, and enhance overall aerodynamic performance. The analysis spans a range of α, from 0° to 20°. The study encompasses two scenarios: one without the incorporation of aerodynamic devices and the other with their implementation. Comparative analyses of increases in CL and reductions in CD, with notable improvements observed at a 15o angle of attack, 28.47% increase in CL and 15.07% decrease in CD. Furthermore, the improved performance in the CL/CD, which increases substantially in stall conditions, thereby demonstrating the potential of these aerodynamic modifications to enhance overall aircraft performance. [ABSTRACT FROM AUTHOR]

The present research enhances the efficiency of an airfoil section from the DTU-10MW Horizontal Axis Wind Turbine (HAWT) via Active Flow Control (AFC) implementation and when using synthetic jets (SJ). The flow around two airfoil sections cut along the wind turbine blade and for a wind speed of 10 m/s is initially simulated using the CFD-2D-RANS-K ω -SST turbulence model, from where the time-averaged boundary layer separation point and the associated vortex shedding frequency are obtained. On a second stage of the paper, and considering one of the two airfoil sections, the boundary layer separation point previously determined is used to locate the SJ groove as well as the groove width; the three remaining AFC parameters, momentum coefficient, jet inclination angle, and jet pulsating frequency, are parametrically optimized. Thanks to the energy assessment presented in the final part of the paper, the study shows that a considerable power increase of the airfoil section can be obtained when attaching the former separated boundary layer. The extension of the optimization process to the rest of the blade sections where the boundary layer is separated would lead to an efficiency increase of the HAWT. The Reynolds numbers associated to the respective airfoil sections analyzed in the present manuscript are R e = 14.088 × 10 6 and R e = 14.877 × 10 6 , the characteristic length being the corresponding chord length for each airfoil. [ABSTRACT FROM AUTHOR]

On 8 April 2024, a rare total solar eclipse (TSE) passed over western New York State (NYS), the first since 1925 and the last one until 2079. The NYS Mesonet (NYSM) consisting of 126 weather stations with 55 on the totality path provides unprecedented surface, profile, and flux data and camera images during the TSE. Here we use NYSM observations to characterize the TSE's impacts at the surface, in the planetary boundary layer (PBL), and on surface fluxes and CO2 concentrations. The TSE‐induced peak surface cooling occurs 17 min after the totality and is 2.8°C on average with a maximum of 6.8°C. It results in night‐like surface inversion, calm winds, and reduced vertical motion and mixing, leading to the shallowing of the PBL and its moistening. Surface sensible, latent and ground heat fluxes all decrease whereas near‐surface CO2 concentration rises as photosynthesis slows down. Plain Language Summary: On 8 April 2024, a rare total solar eclipse (TSE) passed over western New York State (NYS), the first one since 1925 and the last one until 2079. The entire NYS witnessed at least 88% obscuration at the peak of the eclipse. It provides an excellent opportunity to study the impacts of the TSE. The NYS Mesonet (NYSM), an advanced statewide weather network, has 55 stations on the totality path and provides unprecedented measurements of surface meteorological variables, atmospheric vertical profiles, the heat exchange between the atmosphere and the surface and carbon dioxide (CO2) concentration. It enables one to study the TSE in greater details on a regional scale for the first time. This study found that the moon shadow cools the surface by as much as 6.8°C and creates a surface inversion layer. The cooling calms down winds and vertical mixing, leading to less escape of the water vapor and moistening of the air. It also reduces the heat exchange between the surface and the air. Without sunlight, the photosynthesis shuts down, causing a robust rise in near‐surface CO2 concentration. One‐minute camera images provide a fantastic view of the darkening of the sky during the TSE. Key Points: The New York State Mesonet provided unprecedented surface, profile, flux and image data during the 8 April 2024 total solar eclipse across New York StateThe eclipse resulted in significant cooling and moistening near the surface and in the boundary layer, leading to a surface inversion layerIt also weakened surface winds, turbulent mixing, heat fluxes, but caused a robust rise in near‐surface CO2 concentrations [ABSTRACT FROM AUTHOR]

The diffusion process of the nanofibers in a nanofluid displays a combination of diffusion and wave characteristics. The wave-like behavior is attributed to the large aspect ratio of the nanofiber, which can be modeled as a bead-spring chain system, exhibiting wave-like properties. In order to describe this mixed diffusion-wave phenomenon, a fractional diffusion-wave equation is proposed, wherein multiple time fractional derivatives are employed. By appropriately regulating the fractional time terms, the model can be transformed into either the traditional diffusion equation or the traditional wave equation, as required. A numerical schedule is developed through the implementation of suitable time and space discretization, and a stability analysis is conducted. The numerical results substantiate the reliability of the numerical schedule and its applicability to other fractional differential equations. [ABSTRACT FROM AUTHOR]

The article presents the results of the computational and theoretical analysis of the mechanisms of formation of a vortex flow during flame propagation in a channel filled with a combustible gas mixture. It is shown that at the initial stage of combustion development the flow formed in the near-wall region can be described in accordance with the Blasius boundary layer theory. It is established that the nonstationary nature of flame evolution leads to disturbance of the boundary layer and that under certain conditions the development of the boundary layer instability determines the formation of vortex structures ahead of the flame front. In such a case, depending on the composition of the combustible gas mixture, stable flow regimes without the formation of vortex structures and regimes with flow stabilization are possible. [ABSTRACT FROM AUTHOR]

SUMMARY: Understanding the relationship between wind speed and gas exchange in plants is a longstanding challenge. Our aim was to investigate the impact of wind speed on maximum rates of gas exchange and the kinetics of stomatal responses. We conducted experiments in different angiosperm and fern species using an infrared gas analyzer equipped with a controlled leaf fan, enabling precise control of the boundary layer conductance. We first showed that the chamber was adequately mixed even at extremely low wind speed (<0.005 m s−1) and evaluated the link between fan speed, wind speed, and boundary layer conductance. We observed that higher wind speeds led to increased gas exchange of both water vapor and CO₂, primarily due to the increase in boundary layer conductance. This increase in transpiration subsequently reduced epidermal pressure, leading to stomatal opening. We documented that stomatal opening in response to light was 2.5 times faster at a wind speed of 2 m s−1 compared to minimal wind speed in Vicia faba, while epidermal peels in a buffer with no transpiration exhibited a similar opening rate. The increase in stomatal conductance under high wind was also observed in four angiosperm species under field conditions, but it was not observed in Boston fern (Nephrolepis exaltata), which lacks epidermal mechanical advantage. Our findings highlight the significant impact of boundary layer conductance on determining gas exchange rates and the kinetics of gas exchange responses to environmental changes. Our study shows that wind speed significantly influences stomatal conductance and kinetics, highlighting the importance of accounting for wind effects when measuring leaves using gas exchange systems. These findings enhance our understanding of the mechanisms controlling stomatal conductance and plant responses to environmental changes. [ABSTRACT FROM AUTHOR]

Persistent cold air pools (CAPs) trap pollutants in valleys for extended periods, leading to reduced visibility and increased air pollution within these valleys. The structure of the persistent cold air pool that occurred in the Lanzhou Valley in December 2016 was simulated using different Planetary Boundary Layer (PBL) scenarios of the Weather Research and Forecasting (WRF) model, and the simulation of the persistent cold air pool was further optimized in these PBL scenarios. The simulation results indicated that weather‐scale dry subsidence and nighttime ground radiation cooling were significant factors contributing to the accumulation of persistent CAPs and pollutants in the Lanzhou Valley. In contrast, convective lifting from the ground led to the dissipation of persistent CAPs and a reduction in pollution within the valley. During persistent CAPs, the PM2.5 concentration and valley heat deficit (Q) were 66.7% and 62% higher, respectively, than during non‐CAP. In the original MYNN scheme, the average PBL height, double turbulent kinetic energy (QKE), and turbulence length scale during persistent CAPs decreased by 30.79%, 50.5%, and 34.4%, respectively, compared to non‐CAP. Compared with the original MYNN scheme, the optimized MYNN scheme shows a significant improvement in the turbulence simulation during the sustained CAPs, resulting in a more stable atmosphere. The PBL height during the sustained CAPs is reduced by 28 m, the diurnal turbulence length scale is reduced by 31.62%, the stability parameter is reduced by 39%, the diurnal mean QKE is reduced by 27.45%, and the QKE impact height is reduced by 100–400 m. [ABSTRACT FROM AUTHOR]

This study utilizes the wealth of observational data collected during the recent Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) drift experiment to constrain and evaluate close to two‐hundred daily Large‐Eddy Simulations (LES) of Arctic boundary layers and clouds at high resolutions. A standardized approach is adopted to tightly integrate field measurements into the experimental configuration. Covering the full drift represents a step forward from single‐case LES studies, and allows for a robust assessment of model performance against independent data under a range of atmospheric conditions. A homogeneously forced domain is simulated in a Lagrangian frame of reference, initialized with radiosonde and value‐added cloud profiles. Prescribed boundary conditions include various measured surface characteristics. Time‐constant composite forcing is applied, primarily consisting of subsidence rates sampled from reanalysis data. The simulations run for 3 hours, allowing turbulence and clouds to spin up while still facilitating direct comparison to MOSAiC data. Key aspects such as the vertical thermodynamic structure, cloud properties, and surface energy fluxes are well reproduced and maintained. The model captures the bimodal distribution of atmospheric states that is typical of Arctic climate. Selected days are investigated more closely to assess the model's skill in maintaining the observed boundary layer structure. The sensitivity to various aspects of the experimental configuration and model physics is tested. The model input and output are available to the scientific community, supplementing the MOSAiC data archive. The close agreement with observed meteorology justifies the use of LES for gaining further insight into Arctic boundary layer processes and their role in Arctic climate change. Plain Language Summary: The Arctic is one of the regions most affected by global climate change, warming up to four times as fast as the rest of the globe. It is also a particularly inaccessible region to conduct measurements. Fortunately, between 2019 and 2020 the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) campaign collected an unprecedented amount of data in the Arctic. In this study, numerous of these measurements are incorporated into high‐resolution computer simulations of the lowest part of the Arctic atmosphere. This simulation data complements and contextualizes the observations and enables insight into complex physical processes, for example, cloud formation, cloud ice production, or turbulent mixing. The Arctic is an extreme place, and models often struggle to represent the atmosphere accurately. Therefore, the main achievement of this study is to successfully simulate 190 atmospheric situations as measured during the campaign. The generated data set performs well when compared to independent observations. Single cases deliver information about individual atmospheric conditions, and the collection gives insight into how key climate variables behaved throughout the MOSAiC year. Key Points: A standardized LES setup based on campaign data is developed with an aim to supplement the local measurements during the Multidisciplinary drifting Observatory for the Study of Arctic Climate driftIndependent drift‐long statistics on key aspects of the surface energy budget, thermodynamic structure, and clouds are reproducedSensitivity tests indicate microphysics, ice‐radiation interaction and surface representation are critical for successful daily simulations [ABSTRACT FROM AUTHOR]

In the realm of marine science and engineering, hydrofoils play a pivotal role in the efficiency and performance of marine turbines and water-jet pumps. In this investigation, the boundary layer characteristics of an NACA0009 hydrofoil with a blunt trailing edge are focused on. The effectiveness of both the two-equation gamma theta (γ-Reθt) transition model and the one-equation intermittency (γ) transition model in forecasting boundary layer behavior is evaluated. When considering natural transition, these two models outperform the shear stress transport two-equation (SST k-ω) turbulence model, notably enhancing the accuracy of predicting boundary layer flow distribution for chord-length Reynolds numbers (ReL) below 1.6 × 106. However, as ReL increases, both transition models deviate from experimental values, particularly when ReL is greater than 2 × 106. The results indicate that the laminar separation bubble (LSB) is sensitive to changes in angles of attack (AOA) and ReL, with its formation observed at AOA greater than 2°. The dimensions of the LSB, including the initiation and reattachment points, are found to contract as ReL increases while maintaining a constant AOA. Conversely, an increase in AOA at similar ReL values leads to a reduced size of the LSB. The findings are essential for the design and performance optimization of water-jet pumps, particularly in predicting and flow separation and transition phenomena. [ABSTRACT FROM AUTHOR]

This article aims at the construction and analysis of a computational method for a class of singularly perturbed nonlinear differential equations of convection diffusion type with Cauchy data. A computational method, which is a combination of finite difference operators and a Shishkin mesh, is constructed to solve the class of problems. An algorithm involving the continuation technique together with the suggested computational method is developed to solve the class of problems numerically. The newly developed computational method is proved to be essentially first-order parameter-independent convergent. Computational experiments included support the theoretical results. [ABSTRACT FROM AUTHOR]

The wind profile radar data observed from the western mountainous and eastern plain areas of Chengdu from 1 Jan to 31 Dec 2022 are used to analyse the regional differences in the vertical structure and daily variation trend of the boundary layer wind field. The results indicate that: (1) There are regional differences in the prevailing wind direction of the lower boundary layer in Chengdu, with north-north-east winds prevailing in the western mountains area and northeast winds prevailing in the eastern plains, and without seasonal difference. While there is no regional difference in the prevailing wind direction between the east and west regions in the middle and upper boundary layer, both occur alternately with northeast and southwest winds that are consistent with the direction of the western mountain range. (2) The variation trend of the average wind speed profile in the boundary layer between the eastern and western areas of Chengdu is consistent, and the horizontal wind speed in the eastern plain area of the same detection height layer is greater than that in the western mountains area. (3) The local mountain-valley breeze is significant in the lower boundary layer of the western mountains area; the valley breeze with southeast winds prevails during the day; and the mountain breeze with northwest winds prevails at night. (4) There are no regional differences in the diurnal variation characteristics of horizontal wind speed at the same detection height layer, and it is a “single peak and single valley” type, with a peak during the day and a valley at night. (5) The differences in the wind field characteristics of the lower boundary layer between the east and west regions of Chengdu are mainly caused by the local complex terrain of the western mountains area.

To improve the aerodynamic efficiency of the aircraft, the article focuses on optimizing the L/D ratio through postponement of flow separation. Elevating the L/D ratio leads to diminished drag and delays stall occurrence at high angles of attack (α). The wing geometry adopts NACA 6-digit airfoils, specifically NACA 63418 and NACA 63415. Various aerodynamic devices are explored for their ability to defer flow separation, with passive aerodynamic devices being the prevalent choice. The primary goal of the research is to integrate rotational thin wire elements into the aerodynamic components of the wing, aiming to diminish drag, amplify lift, and enhance overall aerodynamic performance. The analysis spans a range of α, from 0o to 20o. The study encompasses two scenarios: one without the incorporation of aerodynamic devices and the other with their implementation. Comparative analyses of increases in CL and reductions in CD, with notable improvements observed at a 15o angle of attack, 28.47% increase in CL and 15.07% decrease in CD. Furthermore, the improved performance in the CL/CD, which increases substantially in stall conditions, thereby demonstrating the potential of these aerodynamic modifications to enhance overall aircraft performance.

A common denominator between conical-symmetry and conical shock interaction is the spanwise pressure gradient, which perform more non-uniformity and its interaction flow is more complicated than the spanwise-homogeneous planar shock wave. Recent advances in conical-symmetry and conical shock interactions with turbulent boundary layer are reviewed in specific areas: (i) quasi-conical swept interactions due to compression ramps and sharp fins, (ii) impinging conical shock wave with interactions of plate wall, (iii) laminar double cone interactions with consideration of real-gas effects. Substantial success has been achieved in describing the phenomena of the time averaged and instantaneous flow features and the low-frequency unsteadiness, including correlations and coherent structures in the separation bubble, through complementary experimental and numerical studies of swept shock interactions. All available observations are here scrutinized to infer underlying mechanisms of interactions in conical flow, and provide theoretical foundation and hints for fluidic control techniques. Comparison with high-fidelity direct numerical simulations is used to quantified the uncertainty of RANS turbulence models in complex interactions. Regarding heat transfer, extensive studies of hypersonic flow over double cone geometries have shown that those can be predicted with reasonable accuracy, even in the presence of high-temperature effects.

Abstract Neglecting the velocity reduction in the induction zone of wind turbines can lead to overestimates of power production predictions, and, thus, of the annual energy production for a wind farm. An experimental study on the induction zone associated with wind turbine operations is performed in the boundary‐layer test section of the BLAST wind tunnel at UT Dallas using stereo particle image velocimetry. This experiment provides a detailed quantification of the wind speed decrease associated with the induction zone for two different incoming flows, namely, uniform flow and boundary layer flow. Operations of wind turbines in different regions of the power curve are modeled in the wind tunnel environment with two porous disks with a solidity of 50.4% and 32.3%, which correspond to thrust coefficients of 0.71 and 0.63, respectively. The porous disks are designed to approximate the wake velocity profiles previously measured for utility‐scale wind turbines through scanning wind LiDARs. The results show that the streamwise velocity at one rotor diameter upwind of both disks decreases 1% more for the boundary layer flow than for the uniform flow. Further, the effect of shear in front of the disk with a higher thrust coefficient can be observed until 1.75 rotor diameter upwind of the disk, whereas for the disk with a lower thrust coefficient, the effect of shear becomes negligible at 1.25 rotor diameter upwind. It is found that at one rotor diameter upwind, for both incoming flows, the disk having a higher thrust coefficient has 2% more velocity reduction than the lower‐thrust‐coefficient disk. The results suggest that the variability in wind shear and rotor thrust coefficient, which is encountered during typical operations of wind turbines, should be considered for the development of improved models for predictions of the rotor induction zone, the respective cumulative effects in the presence of multiple turbines, namely, wind farm blockage, and more accurate predictions of wind farm power capture.

AbstractWhen a fluid flows over a solid surface, it creates a thin boundary layer where the flow velocity is influenced by the surface through viscosity, and can transition from laminar to turbulent at sufficiently high speeds. Understanding and forecasting the fluid dynamics under these conditions is one of the most challenging scientific problems in fluid dynamics. It is therefore of high interest to formulate models able to capture the nonlinear spatio-temporal velocity structure as well as produce forecasts in a computationally efficient manner. Traditional statistical approaches are limited in their ability to produce timely forecasts of complex, nonlinear spatio-temporal structures which are at the same time able to incorporate the underlying flow physics. In this work, we propose a model to accurately forecast boundary layer velocities with a deep double reservoir computing network which is capable of capturing the complex, nonlinear dynamics of the boundary layer while at the same time incorporating physical constraints via a penalty obtained by a Partial Differential Equation (PDE). Simulation studies on a one-dimensional viscous fluid demonstrate how the proposed model is able to produce accurate forecasts while simultaneously accounting for energy loss. The application focuses on boundary layer data in a water tunnel with a PDE penalty derived from an appropriate simplification of the Navier-Stokes equations, showing improved forecasting by the proposed approach in terms of mass conservation and variability of velocity fluctuation against non-physics-informed methods. Supplementary materials for this article are available online, including a standardized description of the materials available for reproducing the work. [ABSTRACT FROM AUTHOR]

The application of nanofluids has exploded in recent decades to improve the local number, mean Nusselt number, and rate of heat transfer. However, boundary layer equations of nanofluid across a flat plate with radiation have not been studied, and therefore this paper studies them mathematically for the first time. For water-based copper and aluminum oxide nanofluids, a similarity solution is presented in this study, and the subsequent system of the ordinary differential equation (ODE) is numerically solved by the Runge–Kutta method in MATLAB. Two different hydraulic boundary conditions are used in the simulations. In the first, the flow across a moving plate and the direction of the flow are analyzed, while in the second, the flow over a nonlinearly moving plate in a still fluid is investigated. The nanoparticle's boundary layer thickness is found less than the thermal and hydraulics boundary layers. The local Nusselt number and friction factor of both the nanofluids are calculated and compared with the base fluid. The results demonstrate that the friction coefficient is high and the Nusselt number is low for nanoparticles with a high volume fraction. It also revealed that the friction factor for water–aluminum oxide is 16% greater than that for the water–CuO whereas the local Nusselt number for water–aluminum oxide is only 5% more than that for the water–CuO. [ABSTRACT FROM AUTHOR]

We introduce a novel TRUNC finite element in n$$ n $$ dimensions, encompassing the traditional TRUNC triangle as a particular instance. By establishing the weak continuity identity, we identify it as crucial for error estimate. This element is utilized to approximate a modified Poisson equation defined on a convex polytope, originating from the nonlocal electrostatics model. We have substantiated a uniform error estimate and conducted numerical tests on both the smooth solution and the solution with a sharp boundary layer, which align with the theoretical predictions. [ABSTRACT FROM AUTHOR]

Several formulae of the coefficient of friction have emerged after that of Nikuradse such as the Colebrook-White formula adopted for the calculation of the coefficient of friction in penstocks. However, this requires enormous computing resources because it suffers from an implicit nature. Other authors have proposed explicit formulae for this in order to reduce calculation times and make applications easy. Among these authors, the most used are the formulae of Swamee-Jain and Haaland. In this work, we set out to study the evolution of friction on the top, bottom and side walls of a penstock using these two formulas. The objective is to find which one best characterizes the friction and gives results close to the implicit Colebrook-White formula. To carry out this work, the study was made by a numerical approach in FLUENT. It appears that the Swamee-Jain correlation gives values closer to that of Colebrook-White along the walls. On the side walls, Haaland's formula better describes the constancy of friction with a much larger range of values than Swamee-Jain. On the upper and lower walls, the friction has a linear character. [ABSTRACT FROM AUTHOR]

The turbulent transport of momentum, heat, and moisture can impact tropical cyclone intensity. However, representing subgrid‐scale turbulence accurately in numerical weather prediction models is challenging due to a lack of observational data. To address this issue, a case study of Hurricane Maria was conducted to analyse the influence of different free tropospheric turbulence parametrisations on sheared tropical cyclones. The study used the current Met Office Unified Model (MetUM) parametrisation, as well as a parametrisation scheme with significantly reduced free tropospheric mixing length. Convection‐permitting ensemble simulations were performed for both mixing schemes at two initialisation times (four 18‐member ensembles in total), revealing an improvement in the intensity forecasts of Hurricane Maria when the mixing length was decreased in the free troposphere. By implementing this change, the less diffuse simulations presented a drier mid‐level. The resolved downward transport of drier air from the mid‐levels into the inflow layer (so‐called "downdraft ventilation") was thus more effective in reducing the storm's intensity. In contrast to earlier studies, where decreasing the diffusivity in the boundary layer intensified the storm, we show that decreasing the free tropospheric diffusivity can weaken the storm by enhancing shear‐related weakening processes. While this study was performed using the MetUM, the findings highlight the general importance of considering turbulence parametrisation, and show that changes in diffusivity can have different impacts on storm intensity depending on the environment and where the changes are applied. [ABSTRACT FROM AUTHOR]

Cold fog refers to a type of fog that forms when the temperature is below 0°C. It can be composed of liquid, ice, and mixed‐phase fog particles. Cold fog happens frequently over mountainous terrain in the cold season, but it is difficult to predict. Using observations from the Cold Fog Amongst Complex Terrain (CFACT) field campaign conducted in Heber Valley, Utah, in the western United States during January and February of 2022, this study investigates the meteorological conditions in the surface and boundary layers that support the formation of wintertime ephemeral cold fog in a local area of small‐scale mountain valleys. It is found that fog formation is susceptible to subtleties in forcing conditions and is supported by several factors: (1) established high pressure over the Great Basin with associated local clear skies, calm winds, and a stable boundary layer; (2) near‐surface inversion with saturation near the surface and strong moisture gradient in the boundary layer; (3) warm (above‐freezing) daytime air temperature with a large diurnal range, accompanied with warm soil temperatures during the daytime; (4) a period of increased turbulence kinetic energy (above 0.5 m2·s−2), followed by calm conditions throughout the fog's duration; and (5) supersaturation with respect to ice. Then, the field observations and identified supporting factors for fog formation were utilized to evaluate high‐resolution (˜400 m horizontal grid spacing) Weather Research and Forecasting (WRF) model simulations. Results show that the WRF model accurately simulates the mesoscale conditions facilitating cold‐fog formation but misses some critical surface and atmospheric boundary conditions. The overall results from this paper indicate that these identified factors that support fog formation are vital to accurately forecasting cold‐fog events. At the same time, they are also critical fields for the NWP model validation. [ABSTRACT FROM AUTHOR]

The new scheme for deriving the near‐surface wind profiles discussed in Ma (in review) is applied to an Australian Bureau of Meteorology operational convective scale model over various domains. Both the new and conventional schemes' diagnostic 10‐m winds are then verified against Australia‐wide automatic weather station observations. Analyses of bulk statistics reveal that the new scheme's 10‐m wind forecasts have generally better accuracy than the current conventional scheme with a consistent reduction of biases over all domains. A widely recognised diurnal bias pattern of surface wind speed over the land is substantially reduced, and the inclusion of Ekman spiral effect on the 10‐m wind marginally improves statistics of the wind direction during the nighttime. The new scheme introduces no systemic bias, given the histogram of a bulk mean bias is analogiased to a Gaussian distribution, and moves the distribution of diagnostic wind speed closer to that observed. [ABSTRACT FROM AUTHOR]

Winds are predicted on the discrete grid of numerical weather and climate models. Winds distribute nonlinearly on the height in the near‐surface layer, and a 10 m wind prediction within the layer is often diagnosed upon the Monin–Obukhov similarity theory flux–profile relationship determined from winds at the lowest grid level, the near‐surface atmospheric stability, and surface properties, which leads to concerns that systemic biases may be introduced to the diagnosed wind. Algorithms are proposed to derive near‐surface wind profiles from the grid‐based numerical model forecasts at multiple model levels under the framework of momentum conservations with an implicit solution, associated with simple logarithmic plus linear interpolation in exceptional exemptional conditions. The diagnosed wind profile coheres to the model prediction at the grid level and exhibits differences from the profile using the conventional scheme in the quasi‐steady thermal stratification and non‐steady transitional conditions, retreating to the same logarithmic profile in the neutral condition. [ABSTRACT FROM AUTHOR]

This paper is devoted to the numerical analysis of the spontaneous acceleration of a slow flame in a semi-closed channel. In particular, the flow development in the channel ahead of the propagating flame is analyzed. The applied detailed numerical model allows the clear observation of all features intrinsic to the reacting flow evolution in the channel, including the formation of perturbations on the scale of the boundary layer and their further development. In all considered cases, perturbations of the boundary layer emerge in the early stages of flame acceleration and decay afterward. The flow stabilizes more rapidly in a narrow channel, where the velocity profile is close to the Poiseuille profile. At the same time, the compression waves generated in the reaction zone travel along the channel. The interaction between compression waves in the area of combustion products can lead to the formation of shock waves. The effect of shock waves on the flow in the fresh mixture causes an increase in the flame area and a corresponding flame acceleration. In addition, shock waves trigger boundary-layer instability in wide channels. The perturbations of the boundary layer grow and evolve into vortexes, while further vortex–flame interaction leads to significant flame acceleration. [ABSTRACT FROM AUTHOR]

Out of the 45 radars composing the Terminal Doppler Weather Radar (TDWR) network, 21 are located in areas of complex terrain. Their mission to observe low-level wind shear at major airports prone to strong shear-induced accidents puts them in an ideal position to fill critical boundary layer observation gaps within the NEXRAD network in these regions. Retrievals such as velocity azimuth display and velocity volume processing (VVP) are used to create time-height profiles of the boundary layer from radar conical scans, but assume that a wide area around the radar is horizontally homogeneous. This assumption is rarely met in regions of complex terrain. This paper introduces a VVP retrieval with limited radius to make these profiling techniques informative for flows affected by topography. These retrievals can be applied to any operational radar to help examine critical boundary layer processes. VVP retrievals were derived from the TDWR for Salt Lake City International Airport, TSLC, during a summertime high ozone period. These observations highlighted thermally driven circulations and variations in boundary layer depth at high vertical and temporal resolution and provided insight on their influence on air quality. Significance Statement: Residents in many urban areas of the United States are exposed to elevated ozone concentrations during the summer months. In complex terrain, thermally driven circulations and terrain-forced flows affect chemical processes by modulating mixing and transport. A novel technique to monitor local boundary layer conditions on small horizontal length scales was applied to data from the Terminal Doppler Weather Radar located near Salt Lake City International Airport during a multiday high ozone event, and effects of these flows on ozone concentrations are illustrated. This technique can be applied to other operational weather radars to create long-term and real-time records of near-surface processes at high vertical and temporal resolution. [ABSTRACT FROM AUTHOR]

The present paper gives an overview of the recent modelling activities under NATO-STO AVT-349, focussed on the understanding and modelling of boundary layers for incompressible, high-Reynolds-number flows subject to non-equilibrium conditions such as strong pressure gradients, three-dimensionality, and surface roughness and heterogeneity. For this, we consider simpler cases where the above flow conditions are present separately or in a reduced number of combinations. First, we focus on the effect of roughness on the outer flow and the problems associated to its characterisation and prediction, with a particular emphasis on the conditions necessary for outer-layer similarity to hold. We then focus on how the presence of adverse and favourable pressure gradients affects the effect of roughness, and to what extent the figures used to quantify it are still useful under such conditions. We also consider the effect of surface heterogeneity, the shortcomings when modelling it and how these can be addressed. We then focus on the effect on the outer layer of pressure gradients and non-equilibrium conditions, to what extent similarity holds in those conditions, and how RANS models perform for such flows, identifying routes for their improvement to handle pressure gradients and non-equilibrium. We also discuss the use of data-driven and machine-aided methods in closure models. [ABSTRACT FROM AUTHOR]

The coral-seawater interface is an important, highly dynamic microenvironment for reef-building corals. Also known as the concentration boundary layer (CBL), it is a thin layer of seawater bordering the coral surface that dictates the biochemical exchange between the coral colony and bulk seawater. The CBL is thus a key feature that modulates coral metabolism. However, CBL variation among small-polyped coral species remains largely unknown. Therefore, we recorded over 100 profiles of dissolved O2 concentration using microsensors to characterize CBL traits (thickness, surface O2 concentration, and flux) of three small-polyped branching coral species, Acropora cytherea, Pocillopora verrucosa, and Porites cylindrica. Measurements were conducted during light and darkness combined with low or moderate water flow (2 and 6 cm s−1). We found that CBL traits differed among species. CBL thickness was lowest in A. cytherea, while P. verrucosa showed the largest depletion of surface O2 in dark and highest dark flux. In addition, we found that O2 concentration gradients in the CBL occurred with three main profile shapes: diffusive, S-shaped, and complex. While diffusive profiles were the most common profile type, S-shaped and complex profiles were more frequent in P. verrucosa and P. cylindrica, respectively, and prevailed under low flow. Furthermore, profile types differed in CBL thickness and flux. Finally, low flow thickened CBLs, enhanced changes in surface O2 concentration, and reduced flux, compared to moderate flow. Overall, our findings reveal CBL variability among small-polyped branching corals and help understand CBL dynamics in response to changes in light and water flow conditions. [ABSTRACT FROM AUTHOR]