Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BOUNDARY LAYER SEPARATION”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Imaging gate-tunable Tomonaga–Luttinger liquids in 1H-MoSe 2 mirror twin boundaries

One-dimensional electron systems exhibit fundamentally different properties than higher-dimensional systems. For example, electron-electron interactions in one-dimensional electron systems have been predicted to induce Tomonaga-Luttinger liquid behaviour. Naturally occurring grain boundaries in single-layer transition metal dichalcogenides exhibit one-dimensional conducting channels that have been proposed to host Tomonaga-Luttinger liquids, but charge density wave physics has also been suggested to explain their behaviour. Clear identification of the electronic ground state of this system has been hampered by an inability to electrostatically gate such boundaries and tune their charge carrier concentration. Here we present a scanning tunnelling microscopy and spectroscopy study of gate-tunable mirror twin boundaries in single-layer 1H-MoSe 2 devices. Gating enables scanning tunnelling microscopy and spectroscopy for different mirror twin boundary electron densities, thus allowing precise characterization of electron-electron interaction effects. Visualization of the resulting mirror twin boundary electronic structure allows unambiguous identification of collective density wave excitations having two velocities, in quantitative agreement with the spin-charge separation predicted by finite-length Tomonaga-Luttinger liquid theory.

36 MATERIALS SCIENCE↗
DOE OSTI2025

CFD Modeling of Concentration Polarization in High-Performance Flat-Sheet Modules

Membrane-based gas separation is a promising, regeneration-free route for industrial CO₂ capture, but its performance is often limited by concentration polarization (CP). Classical 1-D models neglect CP and thus overpredict fluxes, especially for today’s high-permeance, high-selectivity membranes. Therefore, we built a 3-D CFD model that resolves local species transport, boundary-layer resistance, permeate-side effects, and pressure drop in flat-sheet modules. Using the Concentration Polarization Coefficient (CPC) and the overall Effective Membrane CP Index (EAC), the model quantifies how geometry and operating conditions shape CP intensity. Simulations show CP becomes pronounced for advanced membranes, with feed-channel height emerging as the dominant geometric lever: shrinking the gap mitigates CP but raises pressure drop, forcing a trade-off between mass-transfer enhancement and energy cost. These insights emphasize the need for CFD-level analysis in next-generation module design, where 1-D treatments can severely underestimate required membrane area and misguide scale-up.

carbon capture↗
DOE OSTI2020

Simulating and Evaluating Atmospheric River-Induced Precipitation Extremes Along the U.S. Pacific Coast: Case Studies From 1980–2017

Atmospheric rivers (ARs) are responsible for a majority of extreme precipitation and floodevents along the U.S. West Coast. To better understand the present-day characteristics of AR-relatedprecipitation extremes, a selection of nine most intense historical AR events during 1980–2017 is simulatedusing a dynamical downscaling modeling framework based on the Weather Research and ForecastingModel. We find that the chosen framework and Weather Research and Forecasting Model configurationreproduces both large-scale atmospheric features—including parent synoptic-scale cyclones—as well asthe filamentary corridors of integrated vapor transport associated with the ARs themselves. The accuracyof simulated extreme precipitation maxima, relative to in situ and interpolated gridded observations,improves notably with increasing model resolution, with improvements as large as 40–60% for fine scale(3 km) relative to coarse-scale (27 km) simulations. A separate set of simulations using smoothedtopography suggests that much of these gains stem from the improved representation of complex terrain.Additionally, using the 12 December 1995 storm in Northern California as an example, we demonstratethat only the highest-resolution simulations resolve important fine-scale features—such as localizedorographically forced vertical motion and powerful near hurricane-force boundary layer winds. Given thedemonstrated ability of a targeted dynamical downscaling framework to capture both local extremeprecipitation and key fine-scale characteristics of the most intense ARs in the historical record, we arguethat such a configuration may be highly conducive to understanding AR-related extremes and associatedchanges in a warming climate

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A Shallow-Water Model for Convective Self-Aggregation

Randomly distributed convective storms can self-aggregate in the absence of large-scale forcings. Here we present a 1D shallow-water model to study the convective self-aggregation. This model simulates the dynamics of the planetary boundary layer and represents convection as a triggered process. Once triggered, convection lasts for finite time and occupies finite length. Here we show that the model can successfully simulate self-aggregation, and that the results are robust to a wide range of parameter values. In the simulations, convection excites gravity waves. The gravity waves then form a standing wave pattern, separating the domain into convectively active and inactive regions. We analyze the available potential energy (APE) budget and show that convection generates APE, providing energy for self-aggregation. By performing dimensional analysis, we develop a scaling theory for the size of convective aggregation, which is set by the gravity wave speed, damping time scale, and number density of convective storms. This paper provides a simple modeling framework to further study convective self-aggregation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Indirect Measurements of the Composition of Ultrafine Particles in the Arctic Late-Winter

In this work, we present indirect measurements of size-resolved ultrafine particle composition conducted during the Ocean-Atmosphere-Sea Ice-Snowpack (OASIS) Campaign in Utqiagvik, Alaska, during March 2009. This study focuses on measurements of size-resolved particle hygroscopicity and volatility measured over two periods of the campaign. During a period that represents background conditions in this location, particle hygroscopic growth factors (HGF) at 90% relative humidity ranged from 1.45 to 1.51, which combined with volatility measurements suggest a mixture of ~30% ammoniated sulfates and ~70% oxidized organics. Two separate regional ultrafine particle growth events were also observed during this campaign. Event 1 coincided with elevated levels of H2SO4 and solar radiation. These particles were highly hygroscopic (HGF = 2.1 for 35 nm particles), but were almost fully volatilized at 160 °C. The air masses associated with both events originated over the Arctic Ocean. Event 1 was influenced by the upper marine boundary layer (200–350 m AGL), while Event 2 spent more time closer to the surface (50–150 m AGL) and over open ocean leads, suggesting marine influence in growth processes. Event 2 particles were slightly less hygroscopic (HGF = 1.94 for 35 nm and 1.67 for 15 nm particles), and similarly volatile. We hypothesize that particles formed during both events contained 60–70% hygroscopic salts by volume, with the balance for Event 1 being sulfates and oxidized organics for Event 2. These observations suggest that primary sea spray may be an important initiator of ultrafine particle formation events in the Arctic late-winter, but a variety of processes may be responsible for condensational growth.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

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↗
DOE OSTI2023

Additive Manufacturing of Coreless Flyback Transformers Using Aerosol Jet Printing and Electrochemical Deposition

An additive manufacturing approach combining aerosol jet printing (AJP) and electrodeposition opens a new pathway to the production of lightweight coreless flyback transformer devices for power electronics. AJP of seed layers with resolution on the order of 30 μm is combined with electrodeposition of Cu and Ni for decreased resistance. This combined approach addresses known shortcomings of AJP and electrodeposition. Nanoparticle inks used in AJP of metals have low conductivity versus bulk materials due to their high grain boundary resistance. There is a lack of readily available high-resolution patterning techniques for electrodeposition outside of expensive clean-room-based lithography techniques. Combining these two techniques enables the patterning of high-resolution, high-conductivity components. Here in this manuscript, we report on the construction of coreless flyback transformers consisting of two-layer primary and two-layer secondary spiral inductors separated by layers of a printed UV-curable dielectric. An input voltage of 17 V at 400 kHz was amplified to an output of 1250 V corresponding to a gain of 73.5. COMSOL modeling at the individual inductor level and at the transformer level was used to compare expected inductance, equivalent series resistance (ESR), and coupling with experimentally measured values.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Simulations of divertor designs that spatially separate power and particle exhaust using mid-leg divertor particle pumping

Predictive design modeling of a Dissipation-Focused Divertor for future operation in DIII-D reveals that increasing the poloidal distance of the pump duct entrance from the target surface along the low-field side divertor baffle increases neutral compression and modifies the spatial distribution of power dissipation. With a divertor pump located mid-leg between the target and the X-point, SOLPS-ITER boundary plasma simulations without drifts predict the formation of a dense neutral cloud near the target with > 30x higher neutral compression in detachment, a more stable detachment front located further from the target, and ~25% lower outer midplane separatrix density required for detachment onset, compared to a pump located in the scrape-off layer at the target surface. Up to 19 MW of power flowing into the divertors is modeled using the following two numerical implementations for particle pumping: a specified fraction of particles incident on variable wall sections of the plasma grid is removed from the computational domain (so-called albedo pumping), and a pump duct is modeled which includes dynamics of kinetic neutrals in the duct. The simulations show that the detachment front is located between the divertor target and the X-point and is relatively stable near the pump entrance, without a strong dependence on gas puff rate or injected power. The mid-leg pump design spatially separates the two primary functions of a divertor (power handling and particle exhaust), with the majority of power dissipation occurring near the target plate and particle exhaust taking place further upstream. The benefit of enhanced dissipation using mid-leg pumping comes at the cost of a higher outer midplane separatrix density for a given amount of particle injection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Growth of organized flow coherent motions within a single-stream shear layer: 4D-PTV measurements

Abstract This study investigates the evolution of a single-stream shear layer (SSSL) originating from a wall boundary layer past a backward-facing step. Utilizing a time-resolved 3D-Particle Tracking Velocimetry (4D-PTV) technique, we track the trajectories of fluorescent particles to gain insight into the flow characteristics of the SSSL. A compact water tunnel facility ( $$\textrm{Re}_\tau =1\,240$$ Re τ = 1 240 ) is fabricated to obtain an SSSL with a perpendicular slow entrainment stream past the separation edge. A hybrid interpolation approach that combines ensemble binning and Gaussian weighting is implemented to derive minimally filtered mean and instantaneous lower- and higher-order flow field parameters. Spanwise-dominant coherent motion accompanied by finer flow scales is observed to grow due to flow entrainment through “nibbling” actions of small-scale vortices, “engulfing” by large-scale vortices, and vortex pairing events. Furthermore, the non-zero-speed stream edge grows relatively faster than the zero-speed stream edge, showing a strong asymmetry in mixing composition across a mixing layer. The SSSL reaches self-similarity at a streamwise distance of $$\approx 55\,\theta _{0}$$ ≈ 55 θ 0 , where $$\theta _0$$ θ 0 is the initial momentum thickness from the separation edge, i.e., considerably shorter than reported in previous studies. A literature comparison of growth rate parameters raises intriguing questions regarding a potential inclusive growth scaling unifying the free shear layers. A turbulent kinetic energy (TKE) budget analysis reveals a negative production region immediately downstream of the separation edge attributed to a large positive streamwise gradient of streamwise velocity. In the self-similar region, the phase-averaged flow mapping demonstrates a larger concentration of turbulence production rate around the outer edges of spanwise vortices, specifically at the intersection of braids and vortices. Furthermore, a spatial separation exists in the regions of peak production and dissipation rates within the vortex core region favoring dissipation. The braids exhibit a larger concentration of turbulence diffusion rates, indicating their function as a conduit for exchanging turbulence between neighboring coherent motions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Positive and Negative Ion CIMS Measurements, SGP, May-June 2024

The negative ion measurements were taken using two LTOF-CIMS with two different nitrate inlets. The data labeled Aerodyne_Inlet used an Aerodyne nitrate inlet. Data labeled PCC used a custom transverse inlet. The positive ion measurements were taken using one LTOF-CIMS with hydronium as an ionization reagent using a custom transverse inlet with an inlet voltage difference of 700 V. See https://doi.org/10.1021/acs.jpca.2c01672 for more information. Peaks were autofit using tofware. Positive ion data is published in the form of peak signal/reagent signal. Data files of peak signal/reagent signal for the two mass spectrometers for the negative ion measurements are separate. Concentrations of H2SO4 are combined into one file. When measurements were overlapping, measurements from the Aerodyne nitrate inlet were used.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

The Cascade Effectiveness of 3-Terminal Tandem Photocathode Architectures as Applied to CO 2 Reduction

Cascade catalysis for photoelectrochemical CO 2 reduction (CO 2 R) decouples the overall reaction into sequential steps occurring on separately optimized catalysts (for example, Ag and Cu) between which an intermediate species such as CO is transferred. A 3-terminal tandem (3TT) photovoltaic architecture advantageously holds two different catalytic regions at different potentials under a single illumination source, but its overall efficiency is low. Using a stochastic reaction-diffusion model, we have examined 3TT photocathode design principles, focusing on the coupling of surface chemistry and transport of CO both inside the boundary layer and outward, into the bulk electrolyte. We find that ensuring that the lateral diffusion distance within the boundary layer is short compared to the boundary layer thickness and controlling bulk flow are key, with interdigitated designs showing an overall conversion efficiency improvement by 2 orders of magnitude compared to the side-by-side case without flow. These findings can be adapted for other cascaded architectures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Spatio-spectral optical fission in time-varying subwavelength layers

Transparent conducting oxides are highly doped semiconductors that exhibit favourable optical features compared with metals, including reduced material losses, tuneable electronic and optical properties, and enhanced damage thresholds. Recently, the photonic community has renewed its attention towards these materials, recognizing their remarkable nonlinear optical properties in the near-infrared spectrum. The exceptionally large and ultrafast change in the refractive index, which can be optically induced in these compounds, extends beyond the boundaries of conventional perturbative analysis and makes this class of materials the closest approximation to a time-varying system. Here we report the spatio-spectral fission of an ultrafast pulse trespassing a thin film of aluminium zinc oxide with a non-stationary refractive index. By applying phase conservation to this time-varying layer, our model can account for both space and time refraction and explain, in quantitative terms, the spatial separation of both spectrum and energy. Our findings represent an example of extreme nonlinear phenomena on subwavelength propagation distances, which provides new insights into transparent conducting oxides’ transient optical properties. This can be critical for the ongoing research on photonic time crystals, on-chip generation of non-classical states of light, integrated optical neural networks, ultrafast beam steering and frequency-division multiplexing.

integrated optics↗
DOE OSTI2022

A scale-wise analysis of intermittent momentum transport in dense canopy flows

We investigate the intermittent dynamics of momentum transport and its underlying time scales in the near-wall region of the neutrally stratified atmospheric boundary layer in the presence of a vegetation canopy. This is achieved through an empirical analysis of the persistence time scales (periods between successive zero-crossings) of momentum flux events, and their connection to the ejection–sweep cycle. Using high-frequency measurements from the GoAmazon campaign, spanning multiple heights within and above a dense canopy, the analysis suggests that, when the persistence time scales ( $t_p$ ) of momentum flux events from four different quadrants are separately normalized by $\varGamma _{w}$ (integral time scale of the vertical velocity), their distributions $P(t_p/\varGamma _{w})$ remain height-invariant. This result points to a persistent memory imposed by canopy-induced coherent structures, and to their role as an efficient momentum-transporting mechanism between the canopy airspace and the region immediately above. Moreover, $P(t_p/\varGamma _{w})$ exhibits a power-law scaling at times $t_{p}<\varGamma _{w}$ , with an exponential tail appearing for $t_{p} \geq \varGamma _{w}$ . By separating the flux events based on $t_p$ , we discover that around 80 % of the momentum is transported through the long-lived events ( $t_{p} \geq \varGamma _{w}$ ) at heights immediately above the canopy, while the short-lived ones ( $t_{p} < \varGamma _{w}$ ) only contribute marginally ( $\approx 20\,\%$ ). To explain the role of instantaneous flux amplitudes in momentum transport, we compare the measurements with newly developed surrogate data and establish that the range of time scales involved with amplitude variations in the fluxes tends to increase as one transitions from within to above the canopy.

Mechanics↗
DOE OSTI2024

Reactive transport modeling of organic carbon degradation in marine methane hydrate systems

Abstract Natural methane hydrate has often been observed in sand layers that contain no particulate organic carbon (POC), but are surrounded by organic-rich, fine-grained marine muds. In this paper, we develop a reactive transport model (RTM) of a microbially-mediated set of POC degradation reactions, including hydrolysis of POC driven by extracellular enzymes, fermentation of the resulting high-molecular weight dissolved organic carbon (HMW-DOC), and methanogenesis that consumes low-molecular weight dissolved organic carbon (LMW-DOC). These processes are mediated by two groups of microbes, fermenters and methanogens that are heterogeneously distributed in different lithologies, with the largest numbers of microbes in the large pores of coarse-grained layers. We find that the RTM can reproduce methane hydrate occurrences observed in two different geological environments, at Walker Ridge Site 313-H (Gulf of Mexico) and IODP Site U1325 (Cascadia Margin). We also find that microbes can degrade POC even if they are physically separated, as extracellular enzymes and DOC can diffuse away from where they are produced by microbes. Microbial activity is highest at relatively early times after burial at shallow depths and near lithological boundaries, where concentration gradients transport solutes to intervals that contain the most microbes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

HarDWR - Water Management Area (WMA) Shapefiles

A dataset within the Harmonized Database of Western U.S. Water Rights (HarDWR). For a detailed description of the database, please see the meta-record v2.0. Changelog v2.0 - No changes v1.0 - Initial public release Description Borders of all Water Management Areas (WMAs) across the 11 western-most states of the coterminous United States are available filtered through a single source. The legal name for this set of boundaries varies state-by-state. The data is provided as two compressed shapefiles. One, stateWMAs, contains data for all 11 states. For 10 of those states, Arizona being the exception, the polygons represent the legal management boundaries used by those states to manage their surface and groundwater resources respectively. WMAs refer to the set of boundaries a particular state uses to manage its water resources. Each set of boundaries was collected from the states individually, and then merged into one spatial layer. The merging process included renaming some columns to enable merging with all other source layers, as well as removing columns deemed not required for followup analysis. The retained columns for each boundary are: basinNum - the state provided unique numerical ID; basinName - the state provided English name of the area, where applicable; state - the state name; and uniID - a unique identifier we created by concatenating the state name, and underscore, and the state numerical ID. Arizona is unique within this collection of states in that surface and groundwater resources are managed using two separate sets of boundaries. During our followup analysis (Grogan et al., in review) we decided to focus on one set of boundaries, those for surface water. This is due to the recommendation of our hydrologists that the surface water boundary set is a more realist representation of how water moves across the landscape, as a few of the groundwater boundaries are based on political and/or economic considerations. Therefore, the Arizona surface WMAs are included within stateWMAs. The Arizona groundwater WMAs are provided as a separate file, azGroundWMAs, as a companion to the first file for completeness and general reference. WMA spatial boundary data sources by state: Arizona: Arizona Surface Water Watersheds; Collected February, 2020; https://gisdata2016-11-18t150447874z-azwater.opendata.arcgis.com/datasets/surface-watershed/explore?location=34.158174%2C-111.970823%2C7.50 Arizona: Arizona Ground Water Basins; Collected February, 2020; https://gisdata2016-11-18t150447874z-azwater.opendata.arcgis.com/datasets/groundwater-basin-2/explore?location=34.158174%2C-111.970823%2C7.50 California: California CalWater 2.2.1; Collected February, 2020; https://www.mlml.calstate.edu/mpsl-mlml/data-center/data-entry-tools/data-tools/gis-shapefile-layers/ Colorado: Colorado Water District Boundaries; Collected February, 2020; https://www.colorado.gov/pacific/cdss/gis-data-category Idaho: Idaho Department of Water Resources (IDWR) Administrative Basins; Collected November, 2015; https://data-idwr.opendata.arcgis.com/datasets/fb0df7d688a04074bad92ca8ef74cc26_4/explore?location=45.018686%2C-113.862284%2C6.93 Montana: Collected June, 2019; Directly contacted Montana Department of Natural Resources and Conservation (DNRC) Office of Information Technology (OIT) Nevada: Nevada State Engineer Admin Basin Boundaries; Collected April, 2020 https://ndwr.maps.arcgis.com/apps/mapviewer/index.html?layers=1364d0c3a0284fa1bcd90f952b2b9f1c New Mexico: New Mexico Office of the State Engineer (OSE) Declared Groundwater Basins; Collected April, 2020 https://geospatialdata-ose.opendata.arcgis.com/datasets/ose-declared-groundwater-basins/explore?location=34.179783%2C-105.996542%2C7.51 Oregon: Oregon Water Resources Department (OWRD) Administrative Basins; Collected February, 2020; https://www.oregon.gov/OWRD/access_Data/Pages/Data.aspx Utah: Utah Adjudication Books; Collected April, 2020; https://opendata.gis.utah.gov/datasets/utahDNR::utah-adjudication-books/explore?location=39.497165%2C-111.587782%2C-1.00 Washington: Washington Water Resource Inventory Areas (WRIA); Collected June, 2017; https://ecology.wa.gov/Research-Data/Data-resources/Geographic-Information-Systems-GIS/Data Wyoming: Wyoming State Engineer's Office Board of Control Water Districts; Collected June, 2019; Directly contacted Wyoming State Engineer's Office

Economics↗
DOE OSTI2023

Atomic Dynamics of Multi‐Interfacial Migration and Transformations

Redox-induced interconversions of metal oxidation states typically result in multiple phase boundaries that separate chemically and structurally distinct oxides and suboxides. Directly probing such multi-interfacial reactions is challenging because of the difficulty in simultaneously resolving the multiple reaction fronts at the atomic scale. Using the example of CuO reduction in H 2 gas, a reaction pathway of CuO → monoclinic m-Cu 4 O 3 → Cu 2 O is demonstrated and identifies interfacial reaction fronts at the atomic scale, where the Cu 2 O/m-Cu 4 O 3 interface shows a diffuse-type interfacial transformation; while the lateral flow of interfacial ledges appears to control the m-Cu 4 O 3 /CuO transformation. Together with atomistic modeling, it is shown that such a multi-interface transformation results from the surface-reaction-induced formation of oxygen vacancies that diffuse into deeper atomic layers, thereby resulting in the formation of the lower oxides of Cu 2 O and m-Cu 4 O 3 , and activate the interfacial transformations. In conclusion, these results demonstrate the lively dynamics at the reaction fronts of the multiple interfaces and have substantial implications for controlling the microstructure and interphase boundaries by coupling the interplay between the surface reaction dynamics and the resulting mass transport and phase evolution in the subsurface and bulk.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Collaborative Research: Advancing Arctic Climate Projection Capability at Seasonal to Decadal Scales (Final Technical Report)

The Regional Arctic System Model (RASM) at process resolving configurations has been used to (i) advance understanding of physical processes and feedbacks involved in Arctic amplification and (ii) understand and potentially reduce uncertainty in prediction of arctic climate change at seasonal to decadal scales. RASM consists the atmosphere (Weather and Research Forecasting model, WRF), ocean (Parallel Ocean Program, POP), sea ice (CICE), land hydrology (Variable Infiltration Capacity model, VIC), river routing scheme (RVIC), marine biogeochemistry components and the coupling framework (CPL7). Its domain is pan-Arctic, with the atmosphere and land components configured on a 50-km or 25-km grid and four configurations of the ocean and sea ice components: 1/12°(~9.3km) or 1/48°(~2.4km) and 45 or 60 vertical layers. These RASM configurations have been motivated by the emerging exascale capability for high performance computing to improve model fidelity. The dynamical downscaling of reanalysis allows comparison of RASM results with observations in place and time to: (i) advance system level understanding of physical processes and coupling involved in an event, (ii) optimize model parameter space, (iii) diagnose and reduce model biases and (iv) produce realistic and consistent across all the components initial conditions for predictions and predictability studies, which are all unique capabilities not available in global Earth System Models (ESMs). An evaluation of RASM 1.0 (Cassano et al. 2017) revealed that it had a cold bias over the oceans and a warm bias over land areas due largely to cloud and radiation biases in the model, with too little cloud cover simulated over land and too much cloud cover simulated over sub-polar oceans. This study has motivated an upgrade to WRF version 3.7.1 in RASM and allowed for the inclusion of the radiative impact of convective clouds. A variety of atmospheric physics parameterizations were evaluated against observations (e.g. data from the Arctic Clouds in Summer Experiment (ACSE); Sedlar et al. 2020) to identify an optimal suite of WRF physics options in RASM. The RASM with the optimized WRF physics were used to study the impact of strong mesoscale winds over the ocean around the southern tip of Greenland (DuVivier and Cassano 2016) and their impact on oceanic convection (DuVivier et al. 2017a). Data from the PolarWinds field campaign were used to evaluate WRF boundary layer physics and resolution impacts on the simulation of a Greenland barrier wind event (DuVivier et al. 2017b). The RVIC streamflow routing model has been implemented in RASM to realistically represent high-resolution streamflow processes (Hamman et al. 2017) and to couple the land buoyancy fluxes to the ocean. The RASM-RVIC high-resolution data set of all coastal freshwater fluxes for the Arctic drainage basin and surrounding areas for 1979-2014 was published as a separate product (https://doi.org/10.5281/zenodo.293037). The fidelity of atmospheric momentum transfer to and the response of polar marine Ekman layer in RASM and Community Earth System Model (CESM) was investigated by Roberts et al. (2015). The increased frequency of oceanic flux exchange in CESM, following the RASM guidance, caused a considerable increase in the median inertial ice speed across the Southern Ocean and parts of the Arctic. A comprehensive evaluation of the RASM1.0 atmosphere-ocean-sea ice-land interface was completed by Brunke et al. (2018). RASM was also demonstrated for its capability to simulate extreme events in agreement with observations in space and time (Lee et al. submitted). In particular, the development of three open water events, known as polynyas, have been simulated north of Greenland in February of 2011, 2017 and 2018, in agreement with satellite observations for the past four decades. The optimized RASM sea ice results have been favorably evaluated against satellite observations and a subset of eleven CMIP6 models (Watts et al. submitted). In a complementary project, Jin et al. (2018) have shown that RASM with higher-resolution and new sea-ice processes contributed to lower model errors in sea-ice conditions, concentrations of nutrients and ice algae, in comparison to results from the coarse-resolution (1°) CESM. In two other complementary studies, RASM results were used (i) to explain areas of concentrated use by bowhead whales, the seasonal progression in the use, and the physical environment within those areas (Citta et al. 2015) and (ii) for a synthesis of fall bowhead whales distribution and migration in the Bering-Chukchi-Beaufort (BCB) Sea to investigate whale movements and feeding to the local ocean hydrography and currents (Citta et al. 2018). However, the multi-decadal output from the CESM Large Ensemble yielded unrealistic forcing. Instead, the shorter NCEP CFSv2 9-month forecasts were successfully tested and afforded an increased ensemble size (~30) to demonstrate gains of dynamical downscaling at sub-seasonal to intra-annual time scales. The improved model physics and coupling among RASM model components have yielded more realistic representation of the sea ice cover and consistent across all model components initial conditions. Consequently, RASM demonstrates significant gains compared to simulation of sea ice in the NCEP reanalysis. In addition, RASM 6-month ensemble forecasts yield very realistic sea ice distribution, which demonstrates both significant gains of dynamical downscaling and the continued impact of the initial conditions on forecasts out to 6 months (https://nps.edu/web/rasm/predictions). A follow up study using RASM for dynamical downscaling of the more realistic CESM initialized Decadal Prediction Large Ensemble output is currently ongoing as part of the DOE RGMA HiLAT-RASM project.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Substrate-Directed Underlayer Growth of Bilayer MoS 2 Revealed by Mo Isotope Labeling

Direct control over the vertical formation sequence and stacking registry in van der Waals (vdW) bilayers is essential for device performance and moiré engineering yet difficult to resolve unambiguously with conventional probes. Here, we use Mo isotope labeling in a two-step chemical vapor deposition process to synthesize bilayer MoS 2 and trace its vertical formation on common substrates. By combining site-selective laser thinning, Raman spectroscopy, time-of-flight secondary ion mass spectrometry, and atomic-resolution scanning transimission electron microscopy (STEM), we find a clear substrate dependence: on SiO 2 /Si, the second layer nucleates and grows beneath the first (underlayer), whereas on sapphire, it forms on top (overlayer). Density functional theory indicates that a larger equilibrium interfacial separation and weaker MoS 2 –substrate interactions on amorphous SiO 2 permit confined interfacial diffusion and underlayer nucleation, whereas stronger interactions and smaller separations on sapphire favor overlayer growth. On SiO 2 , confined epitaxy templates commensurate 2H, 3R, and mixed bilayers, as confirmed by second harmonic generation spectroscopy and STEM. During underlayer coalescence, embedded mirror-twin grain boundaries stitch atomically sharp 2H|3R junctions via alternating 4|8 ring motifs. Molecular-dynamics simulations reveal that these alternating 4|8 motifs accommodate interlayer vdW coupling and locally modulate the stacking registry. These results provide mechanistic insight into confined epitaxial growth and establish isotope labeling as a powerful probe of two-dimensional materials synthesis.

MoS2↗