Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “evolving 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Asynchronicity in opposed-piston RCMs: Does it matter?

Rapid Compression Machines (RCMs) are widely utilized to study combustion phenomena at engine-relevant conditions, and significant efforts are typically made to create a quiescent environment, particularly for investigations of autoignition chemistry. Opposed-piston configurations can be advantageous due to shorter compression times and reduced surface area to volume ratios. Each side must be actuated simultaneously, but this can be challenging in practice. These devices, like most RCMs, utilize hydraulics for actuation, speed control and arrestation of the piston at the end of the stroke; there is no mechanical control or linkage of the two piston trajectories. To quantify the magnitudes and effects of piston asynchronous behavior, this work employs both detailed experimental measurements and, for the first time, high-fidelity, Direct Numerical Simulation (DNS). The boundary conditions are carefully considered applying insight from high-resolution linear variable differential transformer (LVDT) measurements of the piston trajectory and a zero-dimensional kinematics model of the piston-shaft assembly. Sufficient resolution in the piston crevice region is used. The complicated fluid dynamical behavior that can evolve during piston compression and the ensuing delay processes due to offset timings from t offset = 0-10 ms is elucidated. It is found that near t offset = 6 ms and beyond, the boundary layer on the face of the first-seating piston can be sufficiently perturbed, due initially to reemergence of gas from the crevice of the firstseating piston, so that the adiabatic core can become degraded at long ignition delay times. Substantial mixing of colder gas into the interior of the reaction chamber can alter the measurements, similar to effects previously observed for improper piston crevice configuration. In conclusion, experimental techniques to mitigate asynchronous behavior are discussed and demonstrated.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Numerical algorithms for water waves with background flow over obstacles and topography

Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

Ambrose, David M.↗
DOE OSTI2021

Boundary layer description of directional polymer crystallisation

Nearly fifty years ago Lovinger and Gryte suggested that the directional crystallization of a polymer was analogous to the quiescent isothermal crystallization experiment but at a supercooling where the crystal growth velocity was equal to the velocity of the moving front. Our experiments showed that this equivalence holds in a detailed manner at low directional velocities. To understand the underlying physics of these situations, we modeled the motion of a crystallization front in a liquid where the left side boundary is suddenly lowered below the melting point (Stefan's problem) but with the modification that the crystallization kinetics follow a version of the Avrami model. Our numerical results surprisingly showed that the results of the polymer analog track with the Stefan results which were derived for a simple liquid that crystallizes completely at its melting point; in particular, the position of the crystal growth-front evolved with time exactly as in the Stefan problem. The numerical solution also showed that the temperature in the immediate vicinity of the growth-front decreased with increasing front velocity, which is in line with Lovinger and Gryte's ansatz. To provide a clear theoretical understanding of these numerical results we derive a boundary layer solution to the governing coupled differential equations of the polymer problem. The analytical results are in agreement with our observations from experiments and numerical computations but show that this equivalence between the small molecule and polymer analog only holds in the limit where the crystallization enthalpy is much larger than the rate at which heat is conducted away in the polymer. In particular, in the context of the temperature profile, the enthalpy generated by the crystallisation process which is spread out over a narrow spatial region can be approximated as a point source whose location and temperature correspond to the Lovinger–Gryte ansatz.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Aitken Mode Aerosols Buffer Decoupled Mid‐Latitude Boundary Layer Clouds Against Precipitation Depletion

Aerosol-cloud-precipitation interactions are a leading source of uncertainty in estimating climate sensitivity. Remote marine boundary layers where accumulation mode (~100–400 nm diameter) aerosol concentrations are relatively low are very susceptible to aerosol changes. These regions also experience heightened Aitken mode aerosol (~10–100 nm) concentrations associated with ocean biology. Aitken aerosols may significantly influence cloud properties and evolution by replenishing cloud condensation nuclei and droplet number lost through precipitation (i.e., Aitken buffering). We use a large-eddy simulation with an Aitken-mode enabled microphysics scheme to examine the role of Aitken buffering in a mid-latitude decoupled boundary layer cloud regime observed on 15 July 2017 during the Aerosol and Cloud Experiments in the Eastern North Atlantic flight campaign: cumulus rising into stratocumulus under elevated Aitken concentrations (~100–200 mg -1 ). In situ measurements are used to constrain and evaluate this case study. Our simulation accurately captures observed aerosol-cloud-precipitation interactions and reveals time-evolving processes driving regime development and evolution. Aitken activation into the accumulation mode in the cumulus layer provides a reservoir for turbulence and convection to carry accumulation aerosols into the drizzling stratocumulus layer above. Further Aitken activation occurs aloft in the stratocumulus layer. Together, these activation events buffer this cloud regime against precipitation removal, reducing cloud break-up and associated increases in heterogeneity. We examine cloud evolution sensitivity to initial aerosol conditions. With halved accumulation number, Aitken aerosols restore accumulation concentrations, maintain droplet number similar to original values, and prevent cloud break-up. Without Aitken aerosols, precipitation-driven cloud break-up occurs rapidly. In this regime, Aitken buffering sustains brighter, more homogeneous clouds for longer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

TRACER-Uncrewed Aircraft System (TRACER-UAS) Field Campaign Report

Convective processes drive the development of clouds across a variety of scales, thereby playing a central role in governing weather and climate globally. These processes help to set the stage for circulations that redistribute heat and moisture throughout the lower atmosphere, and simultaneously develop storms that supply precipitation to regions that depend on such rain events for supporting a variety of societal needs. At the same time, extreme precipitation events can result in localized flooding of low-lying areas, posing hazards to transportation and regional inhabitants, and potentially damaging structures and infrastructure. As a result, there is an inherent need to better understand the drivers and modulators of clouds and precipitation and leverage such understanding to make projections of how these important systems may evolve with a changing climate. With nearly three billion of Earth’s human inhabitants living within 200 km of a shoreline, development of such understanding is particularly important in coastal regimes. Central drivers of convective cloud development and life cycle in coastal environments include synoptic forcing, local meso- and microscale circulation regimes, microphysical evolution, and aerosol properties. Atmospheric boundary layer (ABL) development and its evolution can help drive the formation and organization of convective clouds, which then impart their own controls through dynamic processes related to latent heating and cooling (including cold-pool formation and propagation), three-dimensional radiative effects, and background environmental conditions (Fan et al. 2016). Such clouds have been contributing to record precipitation events in southeast Texas in recent years, relative to the last several decades of data from regional precipitation gauges (Fagnant et al. 2020). The greater Houston urban area has experienced particularly large increases in the frequency of extreme rainfall events, with enhanced intensity of extreme events (stretching of distribution tails) having potentially significant implications for watershed floodplain planning and mapping in greater Houston and Harris County. While the most extreme precipitation events were associated with tropical cyclones (e.g., Tropical Storm Allison and Hurricane Harvey; Blood 2014, Kao et al. 2019), strong thunderstorms can also deliver large amounts of precipitation over the area over relatively short periods.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

How Well Does the DOE Global Storm Resolving Model Simulate Clouds and Precipitation Over the Amazon?

This study assesses a 40-day 3.25-km global simulation of the Simple Cloud-Resolving E3SM Model (SCREAMv0) using high-resolution ground-based observations from the Atmospheric Radiation Measurement (ARM) Green Ocean Amazon (GoAmazon) field campaign. SCREAMv0 reasonably captures the diurnal timing of boundary layer clouds yet underestimates the boundary layer cloud fraction and mid-level congestus. SCREAMv0 well replicates the precipitation diurnal cycle, however it exhibits biases in the precipitation cluster size distribution compared to scanning radar observations. Specifically, SCREAMv0 overproduces clusters smaller than 128 km, and does not form enough large clusters. Such biases suggest an inhibition of convective upscale growth, preventing isolated deep convective clusters from evolving into larger mesoscale systems. This model bias is partially attributed to the misrepresentation of land-atmosphere coupling. This study highlights the potential use of high-resolution ground-based observations to diagnose convective processes in global storm resolving model simulations, identify key model deficiencies, and guide future process-oriented model sensitivity tests and detailed analyses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Measurement-driven large-eddy simulations of a diurnal cycle during a wake-steering field campaign

Abstract. High-fidelity flow modeling with data assimilation enables accurate representation of the wind farm operating environment under realistic, nonstationary atmospheric conditions. Capturing the temporal evolution of the turbulent atmospheric boundary layer is critical to understanding the behavior of wind turbines under operating conditions with simultaneously varying inflow and control inputs. This paper has three parts: the identification of a case study during a field evaluation of wake steering; the development of a tailored mesoscale-to-microscale coupling strategy that resolved local flow conditions within a large-eddy simulation (LES), using observations that did not completely capture the wind and temperature fields throughout the simulation domain; and the application of this coupling strategy to validate high-fidelity aeroelastic predictions of turbine performance and wake interactions with and without wake steering. The case study spans 4.5 h after midnight local time, during which wake steering was toggled on and off five times, achieving yaw offset angles ranging from 0 to 17°. To resolve nonstationary nighttime conditions that exhibited shear instabilities, the turbulence field was evolved starting from the diurnal cycle of the previous day. These background conditions were then used to drive wind farm simulations with two different models: an LES with actuator disk turbines and a steady-state engineering wake model. Subsequent analysis identified two representative periods during which the up- and downstream turbines were most nearly aligned with the mean wind direction and had observed yaw offsets of 0 and 15°. Both periods corresponded to partial waking on the downstream turbine, which had errors in the LES-predicted power of 4 % and 6 %, with and without wake steering. The LES was also able to capture conditions during which an upstream turbine wake induced a speedup at a downstream turbine and increased power production by up to 13 %.

17 WIND ENERGY↗
DOE OSTI2026

Scaling kinetic Monte-Carlo simulations of grain growth with combined convolutional and graph neural networks

Graph neural networks (GNN) have emerged as a promising machine learning method for microstructure simulations such as grain growth. However, accurate modeling of realistic grain boundary networks requires large simulation cells, which GNN has difficulty scaling up to. To alleviate the computational costs and memory footprint of GNN, we suggest a hybrid architecture combining a convolutional neural network (CNN) based bijective autoencoder to compress the spatial dimensions, and a GNN that evolves the microstructure in the latent space of reduced spatial sizes. Our results demonstrate that the new design significantly reduces computational costs with using fewer message passing layer (from 12 down to 3) compared with GNN alone. The reduction in computational cost becomes more pronounced as the spatial size increases, indicating strong computational scalability. For the largest mesh evaluated (160 3 ), our method reduces memory usage and runtime in inference by 117× and 115×, respectively, compared with GNN-only baseline. More importantly, it shows higher accuracy and stronger spatiotemporal capability than the GNN-only baseline, especially in long-term testing. Such combination of scalability and accuracy is essential for simulating realistic material microstructures over extended time scales. The improvements can be attributed to the bijective autoencoder’s ability to compress information losslessly from spatial domain into a high dimensional feature space, thereby producing more expressive latent features for the GNN to learn from, while also contributing its own spatiotemporal modeling capability. Training data are generated from stochastic grain growth simulations, providing realistic variability for learning robust microstructure evolution. Comprehensive system validation confirms that the model is accurate, robust, and scalable.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Tethered balloon-borne profile measurements of atmospheric properties in the cloudy atmospheric boundary layer over the Arctic sea ice during MOSAiC: Overview and first results

The tethered balloon-borne measurement system BELUGA (Balloon-bornE moduLar Utility for profilinG the lower Atmosphere) was deployed over the Arctic sea ice for 4 weeks in summer 2020 as part of the Multidisciplinary drifting Observatory for the Study of Arctic Climate expedition. Using BELUGA, vertical profiles of dynamic, thermodynamic, aerosol particle, cloud, radiation, and turbulence properties were measured from the ground up to a height of 1,500 m. BELUGA was operated during an anomalously warm period with frequent liquid water clouds and variable sea ice conditions. Three case studies of liquid water phase, single-layer clouds observed on 3 days (July 13, 23, and 24, 2020) are discussed to show the potential of the collected data set to comprehensively investigate cloud properties determining cloud evolution in the inner Arctic over sea ice. Simulated back-trajectories show that the observed clouds have evolved within 3 different air masses (“aged Arctic,” “advected over sea ice,” and “advected over open ocean”), which left distinct fingerprints in the cloud properties. Strong cloud top radiative cooling rates agree with simulated results of previous studies. The weak warming at cloud base is mostly driven by the vertical temperature profile between the surface and cloud base. In-cloud turbulence induced by the cloud top cooling was similar in strength compared to former studies. From the extent of the mixing layer, it is speculated that the overall cloud cooling is stronger and thus faster in the warm oceanic air mass. Larger aerosol particle number concentrations and larger sizes were observed in the air mass advected over the sea ice and in the air mass advected over the open ocean.

Arctic↗
DOE OSTI2021

Large-Eddy Simulations of Marine Boundary Layer Clouds Associated with Cold-Air Outbreaks during the ACTIVATE Campaign. Part I: Case Setup and Sensitivities to Large-Scale Forcings

Large-eddy simulation (LES) is able to capture key boundary layer (BL) turbulence and cloud processes. Yet, large-scale forcing and surface turbulent fluxes of sensible and latent heat are often poorly prescribed for LESs. We derive these quantities from measurements and reanalysis obtained for two cold-air outbreak (CAO) events during Phase I of the Aerosol Cloud Meteorology Interactions over the Western Atlantic Experiment (ACTIVATE) in February–March 2020. We study the two contrasting CAO cases by performing LES and test the sensitivity of BL structure and clouds to large-scale forcings and turbulent heat fluxes. Profiles of atmospheric state and large-scale divergence and surface turbulent heat fluxes obtained from ERA5 data agree reasonably well with those derived from ACTIVATE field measurements for both cases at the sampling time and location. Therefore, we adopt the time-evolving heat fluxes, wind, and advective tendencies profiles from ERA5 data to drive the LES. We find that large-scale thermodynamic advective tendencies and wind relaxations are important for the LES to capture the evolving observed BL meteorological states characterized by the hourly ERA5 data and validated by the observations. We show that the divergence (or vertical velocity) is important in regulating the BL growth driven by surface heat fluxes in LESs. The evolution of liquid water path is largely affected by the evolution of surface heat fluxes. The liquid water path simulated in LES agrees reasonably well with the ACTIVATE measurements. This study paves the path to investigate aerosol–cloud–meteorology interactions using LES informed and evaluated by ACTIVATE field measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The subsurface structure of abraded Al–Zn–Mg–Cu alloy

Surface abrasion has been shown to introduce heavy shear stress and produce a few hundred nanometers thick altered surface layer (ASL) on Al alloys. Such ASLs exhibit microstructural changes at room temperature and different corrosion resistance than the underlying substrate. Here, we report how sample dimensions and temperature affect the deformed subsurface microstructure evolution on an abraded Al–Zn–Mg–Cu alloy using transmission electron microscopy (TEM) and atom probe tomography. The ASL microstructure evolves differently in the bulk sample than in the thin TEM specimen. In the freshly abraded ASL, the original η'/η precipitates and grain boundaries were deformed toward the abrasion direction. After 7-day natural aging of the freshly abraded bulk sample, well-defined subgrains formed locally in the ASL, and unusual Al 2 Cu (θ) and Mg x Zn y phases precipitated at the subgrain boundaries. However, during natural aging of this TEM specimen taken from the abraded and aged bulk samples, much faster changes occurred, including the continued formation of subgrains and θ particles in the entire ASL, and the diffusion of Zn atoms out of Mg x Zn y particles to form a Zn particle. On the other hand, no obvious microstructural changes occurred in the ASL during 3-month natural aging of the TEM specimen immediately taken from the freshly abraded sample. Finally, this work advances the understanding of how Al alloy microstructures respond during and following any shear deformation, and complex and unexpected effects of sample dimensions on the behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Seismic response of vertical dry storage casks under three-dimensional earthquake motions

Ensuring the long-term seismic safety of dry storage casks (DSCs) is becoming increasingly critical as these systems evolve from temporary to de facto permanent repositories for spent nuclear fuels. Traditional seismic soil–structure interaction (SSI) assessment methods use one-dimensional deconvolution or simplified boundary conditions to model incident waves. Although computationally appealing, simplifying assumptions may alter the seismic risk by neglecting the full complexity of three-dimensional (3D) wave propagation effects. To address this challenge, this paper introduces a novel high-fidelity computational framework that leverages the Domain Reduction Method (DRM) with perfectly matched layers (PML) to accurately transfer complex, 3D seismic wavefields from regional-scale fault-rupture simulations into local-scale finite element models of DSCs. Using broadband, physics-based ground motions from a generic M w 7.0 strike-slip event, both single-cask and multi-cask configurations were investigated under near- and far-field conditions. Emphasis is placed on capturing complex SSI, spatial variability in the ground motion, and nonlinear phenomena such as cask rocking and sliding. Numerical results demonstrate that near-field conditions, where forward directivity and fling-step effects dominate, lead to significantly higher DSC rocking and sliding. Far-field cases, by contrast, generally exhibit modest responses. Incorporating SSI tends to amplify or alter DSC response spectra and introduce response variability, which underscores the need for site-specific evaluations and robust modeling approaches to ensure the seismic integrity of DSCs in interim spent fuel storage installations.

Das, Tonmoy↗
DOE OSTI2025

Staircases of passive and active scalar concentration in cellular flow

This paper develops a unified model for staircase formation in both passive and active scalar systems, building upon prior numerical studies by offering new heuristic and physical insights. While prior studies primarily reported numerical results, they did not explore the underlying unifying physics that governs both types of scalar transport; this work addresses that gap by identifying shared mechanisms across both cases. Results of studies of passive and active scalar staircase formation in cellular flows are presented. Staircase formation in cellular flows occurs due to the interplay of fast mixing within cells and slow transport across the inter-cell boundary. The cell boundary emerges as a de facto transport barrier. Special attention is focused on the effects of cellular fluctuations and noise upon staircase structure. A forced, fluctuating vortex array model is used to drive the underlying flow structure. Cellular Peclet number and staircase profile curvature are identified as figures-of-merit to quantify the resiliency of layering. These are related to simple, multi-scatterer scalar random walk models. Results for Peclet number and curvature scaling with flow excitation are presented. We also study staircases of magnetic potential evolving in two-dimensional magnetohydrodynamics as examples of layering of active scalar concentration. Formation of magnetic potential staircases is indeed observed. Flux expulsion inhibits the intercellular transport of magnetic potential and strengthens staircase barriers. Magnetic staircases can be supported against resistive decay by magnetic potential noise forcing. Implications for staircase formation in magnetic confinement experiments are discussed.

Control theory↗
DOE OSTI2026

Systematic analysis of melt pool dynamics in laser processing of mixed powder feedstocks

Functionally graded materials (FGMs) fabricated via additive manufacturing of blended powders offer the potential to spatially tailor properties for new technologies, such as fusion first-wall systems, turbine blades, and spacecraft. However, processing these materials is difficult due to the multiplicity of processing parameters to optimize, all of which must be changed as substrate material, powder feedstock compositions, and melt pool dynamics evolve. Here, this work systematically evaluates the qualitative and quantitative effects of these variables on the melt pool size, shape, composition, and particle distribution in an exemplar Ti-Ta system, and connects the experimental results to Marangoni flow behavior and phenomena observed in other systems. Increasing laser power linearly increases melt pool size and layer thickness, driving engineering considerations such as part/geometrical tolerances. Decreasing laser velocity changes the melt pool shape from lenticular to convex and reduces chemical homogeneity due to extreme thermal and compositional gradients between the melt pool center and boundaries. Thermophysical property differences between the powder feedstock and substrate material, as well as the directionality of the gradient, affect dilution and melt pool dynamics, which in turn affect the melt pool boundary characteristics, shape, and uniformity. Mixed powder feedstocks of intermediate compositions do not behave according to linear interpolations between single-material endpoints, instead building taller and wider melt pools. As such, it is recommended to quantify process maps for at least one intermediate composition in the FGM or multi-material system of interest to ensure optimized processing parameters, predictable melt pool sizes and shapes, and compositional and spatial precision.

Dissimilar↗
DOE OSTI2025

Reconstruction framework advancements to support streaming for the ePIC detector at the EIC

The ePIC collaboration adopted the JANA2 framework to manage its reconstruction algorithms. This framework has since evolved substantially in response to ePIC’s needs. There have been three main design drivers: integrating cleanly with the Podio-based data models and other layers of the key4hep stack, enabling external configuration of existing components, and supporting timeframe splitting for streaming readout. The result is a unified component model featuring a new declarative interface for specifying inputs, outputs, parameters, services, and resources. This interface enables the user to instantiate, configure, and wire components via an external file. One critical new addition to the component model is a hierarchical decomposition of data boundaries into levels such as Run, Timeframe, PhysicsEvent, and Subevent. Two new component abstractions, Folder and Unfolder, are introduced in order to traverse this hierarchy, e.g. by splitting or merging. The pre-existing components can now operate at different event levels, and JANA2 will automatically construct the corresponding parallel processing topology. This means that a user may write an algorithm once, and configure it at runtime to operate on timeframes or on physics events. Overall, these changes mean that the user requires less knowledge about the framework internals, obtains greater flexibility with configuration, and gains the ability to reuse the existing abstractions in new streaming contexts.

Brei, Nathan [Thomas Jefferson National Accelerato↗
DOE OSTI2021

Idealized Large-Eddy Simulations of Stratocumulus Advecting over Cold Water. Part I: Boundary Layer Decoupling

Abstract We explore the decoupling physics of a stratocumulus-topped boundary layer (STBL) moving over cooler water, a situation mimicking warm-air advection (WADV). We simulate an initially well-mixed STBL over a doubly periodic domain with the sea surface temperature decreasing linearly over time using the System for Atmospheric Modeling large-eddy model. Due to the surface cooling, the STBL becomes increasingly stably stratified, manifested as a near-surface temperature inversion topped by a well-mixed cloud-containing layer. Unlike the stably stratified STBL in cold-air advection (CADV) that is characterized by cumulus coupling, the stratocumulus deck in the WADV is unambiguously decoupled from the sea surface, manifested as weakly negative buoyancy flux throughout the subcloud layer. Without the influxes of buoyancy from the surface, the convective circulation in the well-mixed cloud-containing layer is driven by cloud-top radiative cooling. In such a regime, the downdrafts propel the circulation, in contrast to that in CADV regime for which the cumulus updrafts play a more determinant role. Such a contrast in convection regime explains the difference in many aspects of the STBLs including the entrainment rate, cloud homogeneity, vertical exchanges of heat and moisture, and lifetime of the stratocumulus deck, with the last being subject to a more thorough investigation in Part II. Finally, we investigate under what conditions a secondary stratus near the surface (or fog) can form in the WADV. We found that weaker subsidence favors the formation of fog whereas a more rapid surface cooling rate does not. Significant Statement The low-lying blanketlike clouds, called stratocumulus (Sc), reflect much incoming sunlight, substantially modulating Earth’s temperature. While much is known about how the Sc evolves when it moves over warmer water, few studies examine the opposite situation of Sc moving over colder water. We used a high-resolution numerical model to simulate such a case. When moving over cold water, the Sc becomes unambiguously decoupled from the water surface, distinctive from its warm counterpart in which the Sc interacts with the water surface via intermittent cauliflower-like clouds called cumulus clouds. Such decoupling influences many aspects of the Sc–sea surface system, which combine to alter the ability of the Sc to reflect sunlight, thereby influencing the climate. This work laid the foundation for future work that quantifies the contribution of such a decoupled Sc regime to Earth’s radiative budget and climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Grand-scale Atmospheric Imaging Apparatus (GAIA) and Wind LiDAR Multi-scale Measurements in the Atmospheric Surface Layer

Understanding the organization and dynamics of turbulence structures in the atmospheric surface layer (ASL) is important for fundamental and applied research in different fields, including weather prediction, snow settling, particle and pollutant transport, and wind energy. The main challenges associated with probing and modeling turbulence in the ASL are: i) the broad range of turbulent scales associated with the different eddies present in high Reynolds-number boundary layers ranging from the viscous scale ($\mathcal{O}$(mm)) up to large energy-containing structures ($\mathcal{O}$(km)); ii) the non-stationarity of the wind conditions and the variability associated with the daily cycle of the atmospheric stability; iii) the interactions among eddies of different sizes populating different layers of the ASL, which contribute to momentum, energy, and scalar turbulent fluxes. Creative and innovative measurement techniques are required to probe near-surface turbulence by generating spatio-temporally-resolved data in the proximity of the ground and, at the same time, covering the entire ASL height with large enough streamwise extent to characterize the dynamics of larger eddies evolving aloft. To this aim, the U.S. National Science Foundation sponsored the development of the Grand-scale Atmospheric Imaging Apparatus (GAIA) enabling super-large snow particle image velocimetry (SLPIV) in the near-surface region of the ASL. This inaugural version of GAIA provides a comprehensive measuring system by coupling SLPIV and two scanning Doppler LiDARs to probe the ASL at an unprecedented resolution. A field campaign performed in 2021–2022 and its preliminary results are presented herein elucidating new research opportunities enabled by the GAIA measuring system.

17 WIND ENERGY↗
DOE OSTI2021

Is 3D/2D Passivation a Secret to Success for Polycrystalline Thin-Film Solar Cells?

Three leading thin-film photovoltaic (PV) technologies - cadmium telluride (CdTe), CuIn1-xGaxSe2 (CIGS), and perovskite solar cells (PSCs) - are all polycrystalline, but otherwise appear to have little in common. A comprehensive examination of these technologies, however, reveals a common theme: the formation of two-dimensional (2D) van der Waals materials at three-dimensional (3D) absorber interfaces and grain boundaries. In CdTe, the 2D compound is CdCl2; in CIGS, it is XInSe2 (X= K, Rb, Cs) with X depending on the heavy-alkali post-deposition treatment used; and in lead halide PSCs, PbI2 forms naturally, but many new, more stable, 2D perovskites have also been incorporated. Generally, these 2D interfacial materials are present not by design, but instead have evolved from their 3D counterparts during standard device processing. Here, new data, together with evidence compiled from the literature, are presented to illustrate both the existence of 3D/2D interfaces in CdTe, CIGS, and PSCs, and their correlation with improved passivation and device performance. This suggests that 3D/2D passivation may be a heretofore unappreciated key to successful polycrystalline thin-film PV. Finally, the desired attributes of successful low-dimensional layers are presented with rational design strategies for next generation polycrystalline solar cells.

3D/2D↗