Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary Layer”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28

Identifying Meteorological Drivers for Errors in Modeled Winds along the Northern California Coast

Abstract An accurate wind resource dataset is required for assessing the potential energy yield of floating offshore wind farms that are expected along the California outer continental shelf. The National Renewable Energy Laboratory has developed and disseminated an updated wind resource dataset offshore of California, using the Weather Research and Forecasting Model, referred to as the CA20 dataset. As compared to buoy lidar measurements that have become available recently, the CA20 dataset showed significant positive biases for 100-m wind speeds along Northern California wind energy lease areas. To investigate the meteorological drivers for the model errors, we first consider two 1-yr simulations run with two different planetary boundary layer (PBL) parameterizations: the Mellor–Yamada–Nakanishi–Niino (MYNN) PBL scheme (the chosen configuration in the CA20 dataset) and the Yonsei University PBL scheme (which significantly reduces the bias in modeled winds). By comparing the 1-yr simulations to the concurrent lidar buoy observations, we find that errors are larger with the MYNN PBL scheme in warm seasons. We then dive deeper into the analysis by running simulations for short-term (3-day) case studies to evaluate the sensitivity of initial/boundary condition forcings on model results. By analyzing the short-term simulations, we find that during synoptic-scale northerly flows driven by the North Pacific high and inland thermal low, a coastal warm bias in the MYNN simulation is mainly responsible for the modeled wind speed bias by altering the boundary layer thermodynamics. The results of our analysis will help guide the creation of an updated version of the CA20 dataset.

17 WIND ENERGY↗

Observed volatilization behavior of silicon carbide in flowing hydrogen above 2000 K

In this work, the intrinsic compatibility of silicon carbide (SiC) and hydrogen (H 2 ) at high temperatures (2000-2473 K) and pressure near one atmosphere was evaluated through a combination of thermodynamic calculations and hot hydrogen exposure testing. Thermodynamic calculations predict the decomposition of SiC in a hydrogen environment to form free silicon (Si) and free carbon (C). Free Si is predicted to vaporize from the surface as a volatile species, while free C may interact with H 2 to form the hydrocarbons CH 4 (T < 2100 K) or C 2 H 2 (T > 2100 K). Coupons of high purity chemical vapor deposition (CVD) β-SiC were exposed to slowly flowing hydrogen at temperatures ranging between 2000 and 2473 K. SiC experienced active attack as the result of H 2 exposure, exhibiting linear weight loss kinetics and an Arrhenius dependence of weight loss on exposure temperature. The linear volatilization constant was experimentally evaluated to correspond with an activation energy of 370 ± 18 kJ/mol. Due to the dependence of observed corrosion rates on gas velocity, corrosion of SiC in flowing H 2 was determined to be governed by external mass transfer of volatile Si species through the boundary layer. Experimentally derived mass losses were in good agreement with mass losses predicted by a boundary layer limited gas diffusion model.

36 MATERIALS SCIENCE↗

Atmospheric Turbulent Intermittency Over the Arctic Sea‐Ice Surface During the MOSAiC Expedition

Turbulent motions in the Arctic stable boundary layer are characterized by intermittency, but they are rarely investigated due to limited observations, in particular over the sea‐ice surface. In the present study, we explore the characteristics of turbulent intermittency over the Arctic sea‐ice surface using data collected during the Multidisciplinary drifting Observation for the Study of Arctic Climate expedition from October 2019 to September 2020. We first develop a new algorithm, which performs well in identifying the spectral gap over the Arctic sea‐ice surface. Then the characteristics of intermittency are investigated. It is found that the strength of intermittency increases under the conditions of light surface wind speed, small surface wind speed gradient, and strong surface air temperature gradient. The momentum flux, sensible heat flux, and latent heat flux calculated by raw eddy‐covariance fluctuations are overestimated by 3%, 10%, and 24%, respectively, because submesoscale motions are included. Furthermore, the characteristics of the atmospheric boundary layer structure under various intermittency conditions reveal that strong low‐level jets are favorable to surface turbulent motions that result in weak intermittency, while strong temperature inversions above the surface layer suppress surface turbulent motions and lead to strong intermittency.

54 ENVIRONMENTAL SCIENCES↗

Understanding Low Level Jets in the US Atlantic Offshore

Low level jets (LLJs) in the atmosphere exhibit a local windspeed maximum in the boundary layer, with positive shear beneath the jet and negative shear above the jet. Wind turbines tend to experience increased loads with varying wake recovery characteristics in the presence of an LLJ, therefore understanding the mechanism and impact of LLJs is crucial to wind energy development. The US Mid-Atlantic offshore region is a huge potential wind energy resource, yet LLJs in this area are poorly understood. In particular, the coastal offshore environment does not exhibit the same diurnal cycle that leads to strong LLJs in the well-characterized Great Plains region. In this study, we use the Weather Research and Forecasting model (WRF) plus lidar buoy data to identify case studies for LLJ events in 2020 in the New York Bight. A reduced order model is presented to explain the onset of these events based on the competing effects of baroclinicity and eddy diffusivity. Finally, using the macro-scale WRF, we drive a micro-scale large eddy simulation (LES) to generate a more detailed characterization of the Marine Boundary Layer during an LLJ. Gaining a better understanding of LLJs and their impacts on offshore wind in the mid-Atlantic is crucial for a transition toward renewable energy.

atmospheric boundary layer↗

Bridging Land-Surface Fluxes and Aerosol Concentrations to Triggering Convective Rainfall

During the period August 24-28, 2016, a delayed component the GO-Amazon Project, we launched nine altitude-controlled free balloons (CMET). Smaller than typical rawinsondes, CMET are equipped with altitude control, global communication via Iridium satellite, and aspirated sensors. The aims of our effort were to examine the interactions among convective boundary layer and dual river breeze circulations near the confluence of the Tapajos and Amazon Rivers in the eastern Basin. The week-long field campaign was timed to examine the reestablishment of the breeze circulations shortly after the passage of a strong instability line on August 22. Nine CMET were launched at the Curua-Una hydroelectric dam (2.8S; 54.3W), timed to encounter the Tapajos river breeze front by late afternoon. Soundings were made to establish the thickness of interface between the easterly trade and westerly Tapajos breeze circulation. Careful use of sounding strategies allowed these free balloons to track along the northerly channeled flow in the lowest 300 m above the River. Following the river encounter, balloons tracked to the west, sounding to describe the diurnal course of boundary layer in the forest west of the Tapajos River. The longest flight traveled more than 770 km over three days and twice rested overnight in the rain forest canopy. Ancillary data from surface climate and flux stations as well as the Santarem radiosonde, satellite images will be used to illustrate how the breeze circulations are seen near the surface and how they were disrupted by larger-scale events. Comparisons with HYSPLIT trajectories illustrate how sensitive real trajectories are to the refraction that the encounter with the breeze effects. A second effort was intensive analyses of data from the eddy flux tower at the Tapajós National Forest near Santarém Pará Brazil, with special emphasis on the role of clouds and canopy presence on light absorbance and photosynthetic activity.

54 ENVIRONMENTAL SCIENCES↗

Fine Aerosol Acidity and Water during Summer in the Eastern North Atlantic

Aerosol pH governs many important atmospheric processes that occur in the marine boundary layer, including regulating halogen and sulfur chemistries, and nutrient fertilization of surface ocean waters. In this study, we investigated the acidity of PM1 over the eastern North Atlantic during the Aerosol and Cloud Experiments in Eastern North Atlantic (ACE-ENA) aircraft campaign. The ISORROPIA-II thermodynamic model was used to predict PM1 pH and water. We first investigated the sensitivities of PM1 pH and water predictions to gas-phase NH3 and HNO3 concentrations. Our sensitivity analysis indicated that even though NH3 and HNO3 were present at very low concentrations in the eastern North Atlantic during the campaign, PM1 pH calculations can still be sensitive to NH3 concentrations. Specifically, NH3 was needed to constrain the pH of populations of PM1 that had low mass concentrations of NH4+ and non-volatile cations (NVCs). We next assumed that gas-phase NH3 and HNO3 concentrations during the campaign were 0.15 and 0.09 µg m−3, respectively, based on previous measurements conducted in the eastern North Atlantic. Using the assumption that PM1 were internally mixed (i.e., bulk PM1), we determined that PM1 pH ranged from 0.3–8.6, with a mean pH of 5.0 ± 2.3. The pH depended on both Hair+ and Wi. Hair+ was controlled primarily by the NVCs/SO42− molar ratio, while Wi was controlled by the SO42− mass concentration and RH. Changes in pH with altitude were driven primarily by changes in SO42−. Since aerosols in marine atmospheres are rarely internally mixed, the scenario where non-sea salt species and sea-salt species were present in two separate aerosol modes in the PM1 (i.e., completely externally mixed) was also considered. Smaller pH values were predicted for the aerosol mode comprised only of non-sea salt species compared to the bulk PM1 (difference of around 1 unit on average). This was due to the exclusion of sea-salt species (especially hygroscopic alkaline NVCs) in this aerosol mode, which led to increases in Hair+ values and decreases in Wi values. This result demonstrated that assumptions of aerosol mixing states can impact aerosol pH predictions substantially, which will have important implications for evaluating the nature and magnitude of pH-dependent atmospheric processes that occur in the marine boundary layer.

54 ENVIRONMENTAL SCIENCES↗

Advancing Aerosol Chemical Characterization and Vertical Profiling over the Southern Great Plains Using Uncrewed Aerial Sampling and Offline Aerosol Mass Spectrometry

Recent advancements in uncrewed aerial systems (UASs) and particulate matter (PM) analytical techniques have provided opportunities for atmospheric research. In this study, we deployed the Department of Energy’s fixed-wing ArcticShark UAS to examine PM 2.5 composition at varying altitudes─within and above the planetary boundary layer (PBL)─over the Southern Great Plains atmospheric observatory (SGP). A total of 22 flights were conducted across March, June, and August 2023. Composite filter samples were collected during each flight and analyzed with offline aerosol mass spectrometry (AMS), complemented by on-board real-time sensors and ground-based instrumentation, to provide a comprehensive view of regional aerosol characteristics. Results show clear vertical and seasonal differences in the aerosol composition. Relative to ground-level measurements, aloft samples exhibited shifts in the distribution of organic and inorganic PM, with the organic composition varying distinctly across seasons. Particulate organic nitrogen (ON) was elevated, with bulk compositions similar in March and June but strongly altered in August, likely driven by biomass burning and enhanced photochemical activity. Combined AMS and chemical ionization mass spectrometry analyses detected amines, amides, and amino acids. PM above the planetary boundary layer was enriched in oxidized organic aerosols, while ground-level PM contained higher nitrate and sulfate. Seasonal differences in aqueous-phase processing were also observed, which were strongest in March during persistent cloud cover and weaker in the drier August period, suggesting a shift from aqueous- to gas-phase SOA formation. In conclusion, these findings highlight the value of UAS in advancing PM measurements and vertical profiling of aerosol composition.

54 ENVIRONMENTAL SCIENCES↗

HTR solver: An open-source exascale-oriented task-based multi-GPU high-order code for hypersonic aerothermodynamics

In this study, the open-source Hypersonics Task-based Research (HTR) solver for hypersonic aerothermodynamics is described. The physical formulation of the code includes thermochemical effects induced by high temperatures (vibrational excitation and chemical dissociation). The HTR solver uses high-order TENO-based spatial discretization on structured grids and efficient time integrators for stiff systems, is highly scalable in GPU-based supercomputers as a result of its implementation in the Regent/Legion stack, and is designed for direct numerical simulations of canonical hypersonic flows at high Reynolds numbers. Additionally, the performance of the HTR solver is tested with benchmark cases including inviscid vortex advection, low- and high-speed laminar boundary layers, inviscid one-dimensional compressible flows in shock tubes, supersonic turbulent channel flows, and hypersonic transitional boundary layers of both calorically perfect gases and dissociating air.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Convectively coupled Rossby–Gravity waves in a field campaign: How they are captured in reanalysis products

Abstract Convectively coupled equatorial waves are a significant source of atmospheric variability in the tropics. Current numerical models continue to struggle in simulating the coupled diabatic heating fields that are responsible for the development and maintenance of these waves. This study investigates how the diabatic fields associated with Mixed Rossby–Gravity waves (MRGs) are represented in four reanalysis products by using a unique observational dataset from the TRMM‐KWAJEX (Tropical Rainfall Measuring Mission‐Kwajalein Experiment) field campaign. These reanalyses include ERA5, Japanese 55‐year Reanalysis (JRA‐55), Climate Forecast System Reanalysis (CFSR), and Modern‐Era Retrospective Analysis for Research and Applications (MERRA). We found that all four reanalyses captured the MRG structures in winds and temperature, and to a lesser degree in the humidity field except in the boundary layer. However, only the ERA5 and MERRA reanalyses captured the gradual rise and succession of the diabatic heating from boundary layer turbulence, shallow convection, cumulus congestus, and deep convection within the waves. ERA5 is the only product that also captured the gradual rise of the subgrid‐scale vertical transport of moist static energy. All reanalysis products underestimated the diabatic heating from cumulus congestus. Results provide observational basis on what aspects of MRG can be trusted and what cannot in the reanalysis products.

Wang, Xiaocong↗

Symmetry-breaking-induced rare fluctuations in a time-delay dynamic system

Inspired by the experimental and numerical findings, we study the dynamic instabilities of two coupled nonlinear delay differential equations that are used to describe the coherent oscillations between the top and bottom boundary layers in turbulent Rayleigh–Bénard convection. By introducing two sensitivity parameters for the instabilities of the top and bottom boundary layers, we find three different types of solutions, namely in-phase single-period oscillations, multi-period oscillations and chaos. The chaos solution contains rare but large amplitude fluctuations. The statistical properties of these fluctuations are consistent with those observed in the experiment for the massive eruption of thermal plumes, which causes random reversals of the large-scale circulation in turbulent Rayleigh–Bénard convection. Finally, our study thus provides new insights into the origin of rare massive eruptions and sudden changes of large-scale flow pattern that are often observed in convection systems of geophysical and astrophysical scales.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Assessing turbulence and mixing parameterizations in the gray-zone of multiscale simulations over mountainous terrain during the METEX21 field experiment

Multiscale numerical weather prediction models transition from mesoscale, where turbulence is fully parameterized, to microscale, where the majority of highly energetic scales of turbulence are resolved. The turbulence gray-zone is situated between these two regimes and multiscale models must downscale through these resolutions. Here, we compare three multiscale simulations which vary by the parameterization used for turbulence and mixing within the gray-zone. The three parameterizations analyzed are the Mellor-Yamada Nakanishi and Niino (MYNN) Level 2.5 planetary boundary layer scheme, the TKE-1.5 large eddy simulation (LES) closure scheme, and a recently developed three-dimensional planetary boundary layer scheme based on the Mellor-Yamada model. The simulation domain includes complex (i.e., mountainous) terrain in Nevada that was instrumented with meteorological towers, profiling and scanning lidars, a tethered balloon, and a surface flux tower. Simulations are compared to each other and to observations, with assessment of model skill at predicting wind speed, wind direction and TKE, and qualitative evaluations of transport and dispersion of smoke from controlled releases. This analysis demonstrates that microscale predictions of transport and dispersion can be significantly influenced by the choice of turbulence and mixing parameterization in the terra incognita , particularly over regions of complex terrain and with strong local forcing. This influence may not be apparent in the analysis of model skill, and motivates future field campaigns involving controlled tracer releases and corresponding modeling studies of the turbulence gray-zone.

54 ENVIRONMENTAL SCIENCES↗

Uniform momentum and temperature zones in unstably stratified turbulent flows

Wall-bounded turbulent flows exhibit a zonal arrangement, in which streamwise velocity organizes into uniform momentum zones (UMZs), separated by thin layers of elevated interfacial shear. While significant research efforts have focused on these structural features in neutrally stratified flows, the effects of unstable thermal stratification on UMZs and on analogous uniform temperature zones (UTZs) have not been considered previously. In this article, statistical properties of UMZs and UTZs are investigated using a suite of large eddy simulations of unstably stratified turbulent channel flow spanning weakly to highly convective conditions. When normalized by the friction velocity and stability-dependent mixing length, the mean velocity gradient based on UMZ interfacial velocity jumps and the vorticity thickness exhibits good collapse for all stabilities, establishing a link between UMZ properties and scaling predictions from Monin–Obukhov similarity theory. A similar relationship is found between UTZ properties and surface-layer scaling of the mean temperature gradient. In the mixed layer, mean UMZ depth is quasi-constant with wall-normal distance, while the deepest UTZs are found in the centre of the boundary layer. These instantaneous structures are found to be linked to the well-mixed velocity and temperature profiles in the convective mixed layer. Conditional averaging indicates that both UMZ and UTZ interfaces are associated with ejections of momentum and warm updrafts below the interface and sweeps of momentum and cool downdrafts above the interface. These results demonstrate a tangible connection between instantaneous structural features, mean properties and scaling laws in unstably stratified flows.

Mechanics↗

Large Methane Emissions From the Pantanal During Rising Water‐Levels Revealed by Regularly Measured Lower Troposphere CH 4 Profiles

Abstract The Pantanal region of Brazil is the largest seasonally flooded tropical grassland and, according to local chamber measurements, a substantial CH 4 source. CH 4 emissions from wetlands have recently become of heightened interest because global atmospheric 13 CH 4 data indicate they may contribute to the resumption of atmospheric CH 4 growth since 2007. We have regularly measured vertical atmospheric profiles for 2 years in the center of the Pantanal with the objectives to obtain an estimate of CH 4 emissions using an atmospheric approach, and provide information about flux seasonality and its relation to controlling factors. Boundary layer‐free troposphere differences observed in the Pantanal are large compared to other wetlands. Total emissions based on a planetary boundary layer budgeting technique are 2.0–2.8 TgCH 4 yr −1 (maximum flux ∼0.4 gCH 4 m −2 d −1 ) while those based on a Bayesian inversion using an atmospheric transport model are ∼3.3 TgCH 4 yr −1 . Compared to recent estimates for Amazonia (∼41 ± 3 TgCH 4 yr −1 , maximum flux ∼0.3 gCH 4 m −2 d −1 ) these emissions are not that large. Our Pantanal data suggest a clear flux seasonality with CH 4 being released in large amounts just after water levels begin to rise again after minimum levels have been reached. CH 4 emissions decline substantially once the maximum water level has been reached. While predictions with prognostic wetland CH 4 emission models agree well with the magnitude of the fluxes, they disagree with the phasing. Our approach shows promise for detecting and understanding longer‐term trends in CH 4 emissions and the potential for future wetlands CH 4 emissions climate feedbacks.

Gloor, M.↗

Challenge Problem 1: Preliminary Results of the Direct Numerical Simulation of Transient Flows

This report presents the first direct numerical simulations (DNS) of transient mixed convection in an idealized downcomer-like channel (Challenge Problem 1, Phase II). Using the GPU-accelerated NekRS solver, we modeled a sudden decay in driving pressure, mimicking loss-of-flow events, and tracked the resulting evolution of Reynolds number, boundary-layer structure, turbulence statistics, and heat-transfer metrics. Key findings include the systematic thickening and eventual asymmetry of velocity and thermal boundary layers under buoyant deceleration; minimal “memory” lag in Reynolds shear stress and TKE profiles when sampled at matching Re, yet clear shifts of peak locations toward the cooled wall; overshoots in transient eddy-viscosity and eddy-diffusivity (and corresponding sub-unity turbulent Prandtl numbers) on the cooled side; and a pronounced transient Nusselt-number enhancement driven by wall-temperature inertia and residual eddy mixing. These effects combined to offer a temporary cooling margin above steady-state predictions during reactor LOF transients. Future work will extend this work to a more complex “Case II” geometry (90° turn + lower plenum) and generate multi-Re/Pr datasets for data-driven turbulence closures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Turbulent Flux Measurements of the Near‐Surface and Residual‐Layer Small Particle Events

Abstract According to recent field studies, almost half of the New Particle Formation (NPF) events occur aloft, in a residual layer, near the top of the boundary layer. Therefore, measurements of the meteorological parameters, precursor gas concentrations, and aerosol loadings conducted at the ground level are often not representative of the conditions where the NPFs take place. This paper presents new measurements obtained during the Turbulent Flux Measurements of the Residual Layer Nucleation Particles, conducted at the Southern Great Plains research site. Vertical turbulent fluxes of 3–10 nm‐sized particles were measured using a sonic anemometer and two condensation particle counters with nominal cutoff diameters of 3 nm and 10 nm mounted at the top of the 10‐m telescoping tower. Aerosol number size distribution (5–300 nm) was determined through the ground‐based Scanning Mobility Particle Sizers. The size selected (15–50 nm) particle hygroscopicity was derived with the Humidified Tandem Differential Mobility Analyzer. The ground level observations were supplemented by vertically‐resolved measurements of horizontal and vertical wind speeds and aerosol backscatter. The data analysis suggests that (a) turbulent flux measurements of 3–10 nm particles can distinguish between near‐surface and residual‐layer small particle events; (b) sub‐50 nm particles had a hygroscopicity value of 0.2, suggesting that organic compounds dominate atmospheric nanoparticle chemical composition at the site; and (c) current methodologies are inadequate for estimating the dry deposition velocity of sub‐10 nm particles because it is not feasible to measure particle concentration very near the surface, in the diffusion sublayer.

54 ENVIRONMENTAL SCIENCES↗

First Assessment of Cloud‐Land Coupling in LASSO Large‐Eddy Simulations

Abstract To enhance our understanding of cloud simulations over land, this study provides the first assessment of coupling between cloud and land surface in the Large‐Eddy Simulation (LES) Atmospheric Radiation Measurement Symbiotic Simulation and Observation (LASSO) activity for the shallow convection scenario. The analysis of observation data reveals a diurnal cycle of cloud‐land coupling, which co‐varies with surface fluxes. However, coupled (or decoupled) cumulus clouds are inadequately simulated, manifesting as a too‐high (or low) occurrence frequency during the afternoon. This discrepancy is mirrored by the overestimated cloud liquid water path and cloud‐top height. These overestimations are linked to the overpredicted boundary‐layer development and the easier trigger of shallow convection misrepresented in LES runs. Our study underscores the need to improve the representations of boundary‐layer processes and cloud‐land interactions within LES to better simulate shallow clouds in the future.

58 GEOSCIENCES↗

2023 National Offshore Wind data set (NOW-23)

The 2023 National Offshore Wind data set (NOW-23) is the latest wind resource data set for offshore regions in the United States, which supersedes, for its offshore component, the Wind Integration National Dataset (WIND) Toolkit, which was published about a decade ago and is currently one of the primary resources for stakeholders conducting wind resource assessments in the continental United States. The NOW-23 data set was produced using the Weather Research and Forecasting Model (WRF) version 4.2.1. A regional approach was used: for each offshore region, the WRF setup was selected based on validation against available observations. The WRF model was initialized with the European Centre for Medium Range Weather Forecasts 5 Reanalysis (ERA-5) data set, using a 6-hour refresh rate. The model is configured with an initial horizontal grid spacing of 6 km and an internal nested domain that refined the spatial resolution to 2 km. The model is run with 61 vertical levels, with 12 levels in the lower 300m of the atmosphere, stretching from 5 m to 45 m in height. The MYNN planetary boundary layer and surface layer schemes were used the North Atlantic, Mid Atlantic, Great Lakes, Hawaii, and North Pacific regions. On the other hand, using the YSU planetary boundary layer and MM5 surface layer schemes resulted in a better skill in the South Atlantic, Gulf of Mexico, and South Pacific regions. A more detailed description of the WRF model setup can be found in the WRF namelist files linked at the bottom of this page. For all regions, the NOW-23 data set coverage starts on January 1, 2000. For Hawaii and the North Pacific regions, NOW-23 goes until December 31, 2019. For the South Pacific region, the model goes until 31 December, 2022. For all other regions, the model covers until December 31, 2020. Outputs are available at 5 minute resolution, and for all regions we have also included output files at hourly resolution. The NOW-23 data are provided here as HDF5 files. Examples of how to use the HSDS Service to Access the NOW-23 files are linked below. A list of the variables included in the NOW-23 files is also linked below. No filters have been applied to the raw WRF output.

Array↗

Impact of large-scale effects on mass transfer and concentration polarization in Reverse Osmosis membrane systems

Here, we present well-resolved computational fluid dynamics simulations of a large-scale reverse osmosis membrane-spacer configuration (1 m). Our computational model solves the flow and transport equations with variable solute-dependent properties. We utilize a high resolution computational mesh to resolve all relevant length scales associated with spacer-induced mixing and thin concentration boundary layers. An important contribution of this work is the development of a modified mass-transfer correlation that accounts for the development of the concentration boundary layer along the channel. A set of 2D axisymmetric simulations were performed for a spiral wound module layer with varying cross-flow conditions and spacer diameters which indicate a significant entrance length effect for concentration profile development at lower flow rates while mixing effects dominate at higher flow rates. The mass-transfer correlations at higher flow rates compare well with published correlations while a surrogate model for Sherwood number was obtained that depends on an additional similarity variable that accounted for entrance length effects at lower flow rates. Finally, a large-scale membrane-spacer design relevant to high-pressure reverse osmosis is studied with a non-uniform arrangement of spacers, which indicate a substantial saving in pressure drop (~40%) compared to traditional uniformly spaced pattern with minor variations (~2%) in concentration polarization, product water quality (~1%) and water recovery (~7%) compared to a uniform spacer pattern.

97 MATHEMATICS AND COMPUTING↗