Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial optimization”

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

Dynamic Modeling, Parameter Estimation, and Data Reconciliation of a Supercritical Pulverized Coal-Fired Boiler

A high-fidelity, spatially distributed dynamic model of a supercritical pulverized coal-fired boiler was developed for calculating transient thermal profiles of the boiler tubes, which are often not measured or unmeasurable in industrial settings due to the harsh operating conditions. In this work, rigorous models of the subcritical and supercritical steam properties, which can be highly nonlinear, were used to capture the transition between two-phase and one-phase flow, respectively, as supercritical boilers frequently transition through the critical point during load-following operation. The transient boiler model was also used for data reconciliation and parameter estimation for the purpose of model validation against industrial data. An approach was developed to reduce the size of large-scale reconciliation problem by focusing on the implementation of bias terms for the reconciled variables, which can be estimated in place of the reconciled variables themselves. In addition, a functional approximation of the reconciliation problem was defined so that it can be more simply optimized. Combining this functional approximation optimization with reformulation of the reconciliation problem considering measurement biases can significantly improve the tractability of many such large-scale problems. The validated boiler model was then used to study transient responses under load-following operation with specific attention on the calculation of unmeasurable variables.

42 ENGINEERING↗
DOE OSTI2025

Sparse chronology strategy for integrating seasonal energy storage in capacity expansion models

Here, this study develops the sparse chronology method to enhance the representative period framework in capacity expansion models, enabling the effective integration of long-duration energy storage modeling. Traditional representative period methods cannot capture the state of charge of seasonal energy storage systems because they do not establish effective inter-day linkages to connect the state of charge between periods. The sparse chronology approach addresses this limitation by establishing inter-day linkages that allow state of charge to shift inter-seasonally. At the same time, it groups identical representative days into partitions, applying constraints sparsely and implicitly to reduce computational load further. Validation results demonstrate that this method successfully simulates long-duration energy storage patterns, achieving close alignment with a continuous yearly benchmark model, with seasonal trends and state of charge cycles clearly represented. The computational load analysis reveals that the sparse chronology method efficiently applies constraints on maximum and minimum state of charge limits within the representative day framework, eliminating the need for detailed constraints on each individual day. By partitioning representative days and constraining only the start and end of each partition, the method significantly decreases computational requirements. Simulation results show that sparse chronology closely approximates the continuous yearly method's accuracy, even with as few as 20 representative days, achieving correlation values with the benchmark of nearly 0.9 in state of charge plots. Furthermore, it maintains computational efficiency, requiring only 4 % of the solver time compared to the continuous yearly method with 20 representative days. This approach allows capacity expansion models to incorporate long-duration energy storage with high temporal, spatial, and technological resolution, enabling more detailed modeling for large-scale power systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Optimizing the location and configuration of disaster resilience hubs under transportation and electric power network failures

Natural disasters often result in failures of transportation network components and blackouts that imperil the wellbeing of vulnerable populations. In response to these events, resilience hubs have been proposed as a pre-disaster planning strategy to improve access to critical services. This paper introduces an optimization-based approach to locate and configure electric power-generating resilience hubs considering the possibility of failures in transportation and electric power systems. The model's objective is to identify hub locations and configurations that maximize transportation accessibility to the hubs and maximize the satisfaction of basic energy needs through hub-generated electric power. Besides a budget constraint, the model accounts for limits on the levels of hub energy generation vis-à-vis community energy demands, and on the transportation network distance of communities to hubs. Three heuristics are presented for the proposed planning problem. The first heuristic is a genetic algorithm (GA) with problem-specific solution generation procedures. The other two heuristics implement greedy search techniques. Numerical experiments were conducted, using data from rural Puerto Rico, to illustrate the application of the proposed model and heuristics, and examine their performance. In the numerical experiments, the GA heuristic found better solutions than the greedy heuristics. Additionally, design solutions consisting of spatially dispersed hubs with low energy generation capacity were better than solutions with spatially concentrated high-capacity hubs. Lastly, across a wide range of hub demand scenarios, only a small number of candidate hub locations consistently ranked among the best locations for establishing a hub.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Fragmentation analysis of a bar with the Lip-field approach

The Lip-field approach was introduced in Moës and Chevaugeon (2021) as a new way to regularize softening material models. It was tested in 1D quasistatic in Moës and Chevaugeon (2021) and 2D quasistatic in Chevaugeon and Moës (2021): this paper extends it to 1D dynamics, on the challenging problem of dynamic fragmentation. The Lip-field approach formulates the mechanical problem to be solved as an optimization problem, where the incremental potential to be minimized is the non-regularized one. Spurious localization is prevented by imposing a Lipschitz constraint on the damage field. Here, the displacement and damage field at each time step are obtained by a staggered algorithm, that is the displacement field is computed for a fixed damage field, then the damage field is computed for a fixed displacement field. Indeed, these two problems are convex, which is not the case of the global problem where the displacement and damage fields are sought at the same time. The incremental potential is obtained by equivalence with a cohesive zone model, which makes material parameters calibration simple. A non-regularized local damage equivalent to a cohesive zone model is also proposed. It is used as a reference for the Lip-field approach, without the need to implement displacement jumps. These approaches are applied to the brittle fragmentation of a 1D bar with randomly perturbed material properties to accelerate spatial convergence. Both explicit and implicit dynamic implementations are compared. Favorable comparison to several analytical, numerical and experimental references serves to validate the modeling approach.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Spatial Study 2021: Sensor-Based Time Series of Surface Water Temperature, Specific Conductance, Total Dissolved Solids, pH, and Dissolved Oxygen from across Multiple Watersheds in the Yakima River Basin, Washington, USA (v3)

This dataset supports a broader study examining the drivers of spatial variability in sediment respiration rates in the Yakima River Basin. The dataset provides two-hour time series hydrological and water chemistry sensor data, manual chamber open channel respiration data, handheld sensor water chemistry data, river substrate grain size photos, general environmental context photos, and field metadata (including qualitative information on instream and river corridor characteristics) collected during the same two-week period at 47 sites within multiple rivers throughout the Yakima River Basin in Washington, USA. Grain size photos can be used to improve estimates of channel substrate D50 data. Related sample-based water chemistry data are published separately at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1898914.This dataset is comprised of four main folders, one containing three sensor-specific subfolders and the others containing photographs. The SFA_SpatialStudy_2021_SensorData main data folder includes file-level metadata (FLMD), data dictionary (dd), installation methods, field metadata, Ultrameter water chemistry data, field data collection protocols, international generic sample number (IGSN) mapping file, and a readme file. The “Sensor_Manual_Specifications” subfolder contains pdf files from the manufacturer of each sensor with details on the sensor specifications. Each sensor subfolder (BarotrollAtm, MantaRiver, and MinidotManualChamber) contains a sensor data subfolder for timeseries data and a subfolder for plots and summary statistics. The BarotrollAtm Data subfolder contains In Situ Rugged BaroTROLL pressure and temperature data. The MantaRiver Data subfolder contains Eureka Manta+ 35B multisonde temperature, specific conductance, and pH data. The MinidotManualChamber Data subfolder contains PME MiniDOT Logger dissolved oxygen (mg/L and percent saturation) and temperature data. The folder SFA_SpatialStudy_2021_EnvironmentalContextPhotos contains environmental context photographs and videos. The folders SFA_SpatialStudy_2021_SedimentQuadratPhotos_Part1 and SFA_SpatialStudy_2021_SedimentQuadratPhotos_Part2 contain sediment quadrat photographs. All files are .csv, .pdf, .R, .jpg, .jpeg, .mp4, or .mov. This data package was originally published September 2022. It was updated January 2023 (modified files) and June 2024 (new and modified files). See the change history in data package readme for more details.We acknowledge the Yakama Nation as owners and caretakers of the lands where we collected these data. We thank the Confederated Tribes and Bands of the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate sample collection and optimization of data usage according to their values and worldview.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Investigating the effects of cooperative transmission expansion planning on grid performance during heat waves with varying spatial scales

There is growing recognition of the advantages of interregional transmission capacity to decarbonize electricity grids. A less explored benefit is potential performance improvements during extreme weather events. This study examines the impacts of cooperative transmission expansion planning using an advanced modeling chain to simulate power grid operations of the United States Western Interconnection in 2019 and 2059 under different levels of collaboration between transmission planning regions. Two historical heat waves in 2019 with varying geographical coverage are replayed under future climate change in 2059 to assess the transmission cooperation benefits during grid stress. The results show that cooperative transmission planning yields the best outcomes in terms of reducing wholesale electricity prices and minimizing energy outages both for the whole interconnection and individual transmission planning regions. Compared to individual planning, cooperative planning reduces wholesale electricity prices by 64.3 % and interconnection-wide total costs (transmission investments + grid operations) by 34.6 % in 2059. It also helps decrease greenhouse gas emissions by increasing renewable energy utilization. However, the benefits of cooperation diminish during the widespread heat wave when all regions face extreme electricity demand due to higher space cooling needs. Despite this, cooperative transmission planning remains advantageous, particularly for California Independent System Operator with significant diurnal solar generation capacity. This study suggests that cooperation in transmission planning is crucial for reducing costs and increasing reliability both during normal periods and extreme weather events. It highlights the importance of optimizing the strategic investments to mitigate challenges posed by wider-scale extreme weather events of the future.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Addressing Experimental Challenges in Probing Dark Energy with Accuracy and Precision with the Rubin Observatory Legacy Survey of Space & Time (LSST)

With support from this award, the Stanford group under the leadership of PI Patricia Burchat focused on pixel-level algorithms that are key to the investigation of dark energy with the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). The work was carried out through strong engagement in the DOE-supported LSST Dark Energy Science Collaboration (DESC), for which the PI served as Deputy Spokesperson (07/2019 - 07/2021) and Technical Coordinator (07/2021 - 07/2023). Members of the Burchat group are Full DESC members and are among scientists who have been accepted to join the Rubin Observatory System Integration, Test, and Commissioning (SIT-Com) effort. The specific intellectual focus of the group was the optimization of gravitational lensing as a probe of the evolution of the distribution of dark matter in space and time to learn about the nature of dark energy. The Stanford group pursued an integrated research program to understand in detail, and then remediate or calibrate, systematic effects that could limit the sensitivity of the LSST to cosmic shear as a dark energy probe -- particularly those related to our understanding of the atmospheric point spread function (PSF), challenges associated with blended galaxy images, and impacts of residual camera artifacts. Key to this work was developing, enhancing, and validating the fidelity of simulations that play a key role in the development of our understanding of subtle systematic biases and then using these simulations to develop techniques to characterize and mitigate biases, and associated algorithms for implementation in the LSST or DESC science pipelines. The group’s work included contributions to the open-source GalSim simulation package, which is widely used in the community; simulations and studies of spatial correlations of atmospheric PSFs; studies of temporal correlations of the atmospheric PSF; impacts of chromatic PSFs on cosmic shear calibration; continued studies of anisotropic Gaussian process interpolation of PSF models across the focal plane; and measuring and mitigating impacts of residual camera artifacts on cosmic shear. A primary impact of this work is that it increases the scientific discovery potential of the DOE and NSF investments in the Rubin Observatory Legacy Survey of Space \& Time, by addressing challenges in interpreting astronomical images through more accurate descriptions of the blurring of images due to the atmosphere. Specifically, this work impacts our ability to understand the fundamental nature of dark energy.

Burchat, Patricia Rose↗
DOE OSTI2023

Adjoint Waveform Tomography for Crustal and Upper Mantle Structure the Middle East and Southwest Asia for Improved Waveform Simulations Using Openly Available Broadband Data

We present a new model of radially anisotropic seismic wavespeeds for the crust and upper mantle of a broad region of the Middle East and Southwest Asia (MESWA) derived from adjoint waveform tomography. We inverted waveforms from 192 Global Centroid Moment Tensor earthquakes (MW 5.5-7.0) recorded by over 1000 openly available broadband seismic stations from permanent and temporary networks in the region. Spatial coverage of the available data is highly uneven due to earthquakes clustered along plate boundaries and sparse coverage of open seismic networks in the region. We considered three possible starting models: the SPiRaL global model (Simmons et al., 2021); MEC-1 (Kaviani et al., 2020); and CSEM2.0 (Noe et al., 2023). Because the SPiRaL model provides good fits to the observed waveforms measured by the time-bandwidth product of selected windows in several period bands, provides all the necessary parameters and covers the entire domain we used it for the starting model with the period band 50-100 seconds. Inversion iterations proceeded using time-frequency phase misfits in six stages and 54 total iterations reducing the minimum period to 30 seconds. Our final model, MESWA, provides improved waveform fits compared to the starting model for both the data used in the inversion and an independent validation data set of 66 events. Two metrics of waveform fit (the time-frequency phase misfit used in the optimization and normalized L2 misfit) were both reduced by nearly 60% for both data sets and MESWA provides significantly larger misfit reductions relative to the SPiRaL model than the MEC-1 or CSEM models. We also find that MESWA provides a larger time-bandwidth product of selected windows indicating that more information content of the observed waveforms is explained by MESWA than the other models. Our new model reveals tectonic features imaged by other studies and methods but in a new holistic model of shear and compressional wavespeeds (v S and v P , respectively) with anisotropy covering the crust and uppermost mantle of a larger domain. MESWA has smaller scale-length features and tends to sharpen some features relative to the SPiRaL starting model. Examples include: low crustal v S in the TurkishIranian Plateau, Zagros Mountains, Afghan Central Blocks and Sulaiman Fold Belt; low mantle vSfollowing divergent (Gulf of Aden, Red Sea) and transform (Dead Sea Fault) margins of the Arabian Plate; low and high v S in the mantle beneath the Arabian Shield and Platform, respectively. Low vS is imaged below Cenozoic volcanic centers of the Arabian Peninsula, the so-called Mecca-Madina-Nafud (MMN) Line. Positive anisotropy (v SH > v SV ) is inferred for asthenospheric depths across the region except where up/downwelling may influence fabric alignment (e.g. Afar, Red Sea, Arabian Shield). Elevated vS tracks Makran subduction under southeast Iran. MESWA resembles the SPiRaL model in its long-wavelength structure, but enhances shorter wavelengths features on the order of 200 km and smaller. The resulting model could be used for as a starting model for further improvements, say using waveforms from in-country seismic networks that are not openly available or smaller-scale studies targeting shorter period waveforms. The model also could be used for source characterization and moment tensor inversion to improve earthquake hazard studies and nuclear explosion monitoring.

58 GEOSCIENCES↗
DOE OSTI2024

Estimates of Lake Nitrogen, Phosphorus, and Chlorophyll‐ a Concentrations to Characterize Harmful Algal Bloom Risk Across the United States

Abstract Excess nutrient pollution contributes to the formation of harmful algal blooms (HABs) that compromise fisheries and recreation and that can directly endanger human and animal health via cyanotoxins. Efforts to quantify the occurrence, drivers, and severity of HABs across large areas is difficult due to the resource intensive nature of field monitoring of lake nutrient and chlorophyll‐aconcentrations. To better characterize how nutrients interact with other environmental factors to produce algal blooms in freshwater systems, we used spatially explicit and temporally matched climate, landscape, in‐lake characteristic, and nutrient inventory data sets to predict nutrients and chlorophyll‐aacross the conterminous US (CONUS). Using a nested modeling approach, three random forest (RF) models were trained to explain the spatiotemporal variation in total nitrogen (TN), total phosphorus (TP), and chlorophyll‐aconcentrations across US EPA's National Lakes Assessment (n = 2,062). Concentrations of TN and TP were the most important predictors and, with other variables, the RF model accounted for 68% of variation in chlorophyll‐a. We then used these RF models to extrapolate lake TN and TP predictions to lakes without nutrient observations and predict chlorophyll‐afor ∼112,000 lakes across the CONUS. Risk for high chlorophyll‐aconcentrations is highest in the agriculturally dominated Midwest, but other areas of risk emerge in nutrient pollution hot spots across the country. These catchment and lake‐specific results can help managers identify potential nutrient pollution and chlorophyll‐ahot spots that may fuel blooms, prioritize at‐risk lakes for additional monitoring, and optimize management to protect human health and other environmental end goals.

Environmental Sciences & Ecology↗
DOE OSTI2021

The SvFUL2 transcription factor is required for inflorescence determinacy and timely flowering in Setaria viridis

Inflorescence architecture in cereal crops directly impacts yield potential through regulation of seed number and harvesting ability. Extensive architectural diversity found in inflorescences of grass species is due to spatial and temporal activity and determinacy of meristems, which control the number and arrangement of branches and flowers, and underlie plasticity. Timing of the floral transition is also intimately associated with inflorescence development and architecture, yet little is known about the intersecting pathways and how they are rewired during development. Here, we show that a single mutation in a gene encoding an AP1/FUL-like MADS-box transcription factor significantly delays flowering time and disrupts multiple levels of meristem determinacy in panicles of the C 4 model panicoid grass, Setaria viridis. Previous reports of AP1/FUL-like genes in cereals have revealed extensive functional redundancy, and in panicoid grasses, no associated inflorescence phenotypes have been described. In S. viridis, perturbation of SvFul2, both through chemical mutagenesis and gene editing, converted a normally determinate inflorescence habit to an indeterminate one, and also repressed determinacy in axillary branch and floral meristems. Our analysis of gene networks connected to disruption of SvFul2 identified regulatory hubs at the intersection of floral transition and inflorescence determinacy, providing insights into the optimization of cereal crop architecture.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Modeling of crystallization fouling on a horizontal-tube falling-film evaporator for thermal desalination

Scale formation within horizontal-tube falling-film evaporators is an important issue for thermal desalination, because of the deterioration of heat transfer performance and added maintenance costs. Predicting the fouling process is crucial for desalination plant design, operation, and maintenance. In this paper, a model is developed to predict the crystallization fouling process during seawater falling-film flow over a horizontal tube bundle. The scaling is estimated based on deposition theory, with couple heat transfer including the in-tube steam condensation, conduction through the tube wall and the scale layer, and falling-film evaporation. The spatial and temporal variations of temperature, heat transfer coefficient, and scale thickness are predicted, and the effects of variations in process parameters of steam, feed seawater, and tube properties on the scale thickness, evaporation rate, and heat transfer coefficient for falling-film evaporation are analyzed. Comparisons to existing experimental data show the scaling layer thickness falls in the general range of reported measurements, with a slight overestimation due to neglecting the delay of scaling onset under real conditions. The scale layer thickness increases dramatically as the in-tube steam pressure is increased or the seawater flow rate is decreased. In general, the scaling layer becomes thicker on the lower tubes in the tube bundle, due to the heat and mass transfer of the falling film. The effect of tube material on scaling appears mainly dependent on thermal conductivity; thus, using polymer rather than a stainless steel tube decreases scaling and evaporation rate 96% and 88%, respectively. Finally, this work has potential to guide thermal desalination plant design, material selection, and operating parameter optimization.

42 ENGINEERING↗
DOE OSTI2025

Accelerating Embedding Potential Optimization by Reconstructing the Pseudo-Valence Electron Density

Density functional embedding theory (DFET) enables use of electronic structure methods with higher accuracy than density functional theory in a local region, with applications thus far ranging from (photo/electro)catalysis to reactions in solution. DFET partitions a large collection of atoms into smaller groups that interact via a shared embedding (interaction) potential V emb , determined via functional optimization. The optimized effective potential (OEP) process used to optimize V emb is time-consuming and becomes a computational bottleneck due to sharp, oscillating features of V emb near nuclei. Here, similar to pseudopotential theory, by reconstructing electron densities used in the OEP process from smoother pseudo-valence-only (PVO) electron densities as proxies for total densities of the full system and subsystems, we can retain accuracy in the embedded electronic structure calculations while potentially reducing the overhead of V emb construction, within the projector augmented-wave (PAW) formalism. We explore three different chemical reactions as exemplars to test PVO–DFET, namely, H 2 dissociative adsorption on a Cu(111) surface, H 2 O adsorption on a Pt(111) surface, and aqueous [Ca 2+ –SO 4 2– ] ion-pair formation. The PVO approximation works well for all three systems with minimal loss of accuracy (∼10–70 meV error relative to the original exact-derivative (ED) approach) while accelerating V emb generation for the Cu and Pt systems respectively by 20× and 5×. Given proper numerical convergence parameters, the spatial distributions of differences between PVO- and ED-based V emb outside the core regions are small, explaining the exceptional agreement between the two approaches. Finally, we anticipate that this more efficient PVO–DFET approximation will be useful whenever computation of V emb is much more expensive than subsequent embedded high-level electron correlation calculations.

approximation↗
DOE OSTI2022

Evaluating Basin-Scale Forest Adaptation Scenarios: Wildfire, Streamflow, Biomass, and Economic Recovery Synergies and Trade-Offs

Active forest management is applied in many parts of the western United States to reduce wildfire severity, mitigate vulnerability to drought and bark beetle mortality, and more recently, to increase snow retention and late-season streamflow. A rapidly warming climate accelerates the need for these restorative treatments, but the treatment priority among forest patches varies considerably. We simulated four treatment scenarios across the 3,450 km 2 Wenatchee River basin in eastern Washington, United States. We used a decision support tool (DST) to assess trade-offs and synergies within and among treatments on wildfire risk and smoke emissions, water yield and snow retention, biomass production, and economic return. Treatment scenarios emphasized prescribed burning ( BurnOnly ), biomass production ( MaxBiomass ), gap-based thinning to optimize water yield ( IdealWater ), and a principle-based restoration scenario ( RA1 ). Fire hazard, smoke emissions, and biomass production metrics were evaluated across scenarios using the Forest Vegetation Simulator, and water yields were modeled using the Distributed Hydrology Soil Vegetation Model. Simulations were summarized to both patch- (10 1 –10 2 ha) and subwatershed- (10 3 –10 4 ha) scales, and treatment effects were evaluated against an untreated baseline landscape. We used logic models to rank effect sizes by scenario across metrics along a continuum between −1 (no or weak effect) to +1 (large effect). All treatments produced benefits across one or more ecosystem services and led to synergistic benefits to water yield and wildfire hazard reduction. Tradeoffs among resource benefits were clear in wilderness where reliance on prescribed burning without mechanical treatment increased costs and eliminated the potential for biomass recovery. The BurnOnly scenario improved fire risk metrics and streamflow, but effect sizes were lower compared to other treatments. IdealWater showed the strongest benefits overall, demonstrating the ability to capture multiple resource benefits through spatially explicit thinning. Our study provides a framework for integrating strategic and tactical models that evaluate tradeoffs and synergies gained through varied management approaches. We demonstrate the utility of decision support modeling to enhance management synergies across large landscapes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

CO 2 storage site characterization using ensemble-based approaches with deep generative models

Estimating spatially distributed properties such as permeability from available sparse measurements is a great challenge in efficient subsurface CO 2 storage operations. In this paper, a deep generative model that can accurately capture complex subsurface structure is tested with an ensemble-based inversion method for accurate and accelerated characterization of CO 2 storage sites. We chose Wasserstein Generative Adversarial Network with Gradient Penalty (WGAN-GP) for its realistic reservoir property representation and Ensemble Smoother with Multiple Data Assimilation (ES-MDA) for its robust data fitting and uncertainty quantification capability. WGAN-GP are trained to generate high-dimensional permeability fields from a low-dimensional latent space and ES-MDA then updates the latent variables by assimilating available measurements. Several subsurface site characterization examples including Gaussian, channelized, and fractured reservoirs are used to evaluate the accuracy and computational efficiency of the proposed method and the main features of the unknown permeability fields are characterized accurately with reliable uncertainty quantification. Furthermore, the estimation performance is compared with a widely-used variational, i.e., optimization-based, inversion approach, and the proposed approach outperforms the variational inversion method in several benchmark cases. We explain such superior performance by visualizing the objective function in the latent space: because of nonlinear and aggressive dimension reduction via generative modeling, the objective function surface becomes extremely complex while the ensemble approximation can smooth out the multi-modal surface during the minimization. This suggests that the ensemble-based approach works well over the variational approach when combined with deep generative models at the cost of forward model runs unless convergence-ensuring modifications are implemented in the variational inversion.

42 ENGINEERING↗
DOE OSTI2024

ELM-Wet: Inclusion of a wetland landunit with sub-grid representation of eco-hydrological patches and hydrological forcing in E3SM Land Model (ELM)

Wetlands emit the most biogenic methane (CH4) and present the greatest uncertainty in the global CH4 budget. Modeling these emissions is challenging due to the temporal and spatial variability in wetland structure and CH4 flux rates, along with complex interactions among hydrological, ecological, meteorological, and microbial processes that govern CH4 dynamics. To address these issues, we aim to enhance the accuracy of wetland representation in the U.S. Department of Energy’s Exascale Earth System Model (E3SM) Land Model, ELM. This effort led to the development of ELM-Wet, which incorporates a dedicated wetland landunit with subgrid representation of eco-hydrological patch types. We implement wetland-specific hydrology by imposing site-specific constraints on surface water levels, thereby allowing different patches to sustain varying depths of inundation. Additionally, we refined the calculation of aerenchyma transport diffusivity based on observed conductance across different vegetation types. We validated these enhancements through site-specific simulations of a coastal freshwater wetland, the Salvador WMA Freshwater Marsh (Ameriflux, site ID US-LA2), located in the coast of Louisiana (29.85N,90.29W) at an elevation of 0 m. The site was simulated with ELM-Wet and the default version ELMv1. We use Bayesian Optimization to parameterize CO2 and CH4 fluxes. Eddy covariance observations of CO2 and CH4 fluxes from 2012-2013 were used to train the model and data from 2021 were used for validation. In this repository, we include all the output of all the simulations performed using different versions of ELMv1 and ELM-Wet, all input required to the model, and the Matlab scripts we used to processed the output data. ELM-Wet_flmd.csv includes a detailed description of the datasets files.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Coriolis forces modify magnetostatic ponderomotive potentials

It is possible to produce a ponderomotive effect in a plasma system without time-varying fields, if the plasma flows over spatial oscillations in the field. This can be achieved by superimposing a spatially oscillatory perturbation on a guide field, then setting up an electric field perpendicular to the guide field to drive flow over the perturbation. However, subtle distinctions in the structure of the resulting electric field can entirely change the behavior of the resulting ponderomotive force. Previous work has shown that, in slab models, these distinctions can be explained in terms of the polarization of the effective wave that appears in the co-moving frame. Here, we consider what happens to this picture in a cylindrical system, where the transformation to the co-moving (rotating) frame is not inertial. It turns out that the non-inertial nature of this frame transformation can lead to counterintuitive behavior, partly due to the appearance of parallel (magnetic-field-aligned) electric fields in the rotating frame even in cases where none existed in the laboratory frame. Apart from the academic interest of this study, the practical impact lies in being better able to anticipate the antenna configuration on the plasma periphery of a cylindrical plasma that will lead to optimal ponderomotive barrier formation in the interior plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Experimental study of energy-dependent angular broadening of MeV electron beams for high-resolution imaging in thick samples

In scanning transmission electron microscopy (STEM), spatial resolution is primarily influenced by the projected size of the electron probe within the specimen. In thin samples, a large semi-convergence angle enables a tightly focused beam and sub-nanometer resolution. However, in thick specimens, resolution is fundamentally limited by transverse beam broadening from multiple large-angle scattering events—for example, a probe with 10 mrad angular divergence can broaden by ∼100 nm over a 10 μm path. Since this broadening scales inversely with beam energy, MeV-STEM offers a promising route for high-resolution imaging in thick materials. To quantitatively assess this effect, we performed high-precision measurements at UCLA’s PEGASUS beamline, characterizing beam divergence and intensity profiles for 3–8 MeV electrons transmitted through a wedged-silicon sample of varying thickness. Our results reconcile discrepancies among analytical models and validate Monte Carlo simulations. Here, we find that increasing beam energy from 3.0 to 5.8 MeV reduces angular broadening by a factor of 2.6, with diminishing returns observed at 7.6 MeV. These findings provide a quantitative framework for optimizing MeV-STEM parameters in high-resolution imaging of thick biological and microelectronic specimens, and for guiding beam energy selection in other advanced imaging modes beyond STEM.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Neutral Mode Dominates the Forced Global and Regional Surface Temperature Response in the Past and Future

Abstract Using a large suite of Green's function perturbation experiments, we construct the emergent dynamical operator for surface temperature and extract the most excitable mode. The leading mode turns out to be the most dominant mode excited by the climate forcing of doubling CO 2 , alone capturing 56% of the total spatial variances of the surface temperature response. The pattern of the leading mode is partly organized by the atmospheric annular modes in both hemispheres, and further modulated by the feedbacks from ocean circulations. Though derived from a single model, the leading mode can capture the spatio‐temporal variabilities of the global mean surface temperature over the past century in both the CMIP6 model simulations and observations. Moreover, the leading mode is most efficiently excited by radiative forcings from the midlatitude bands centered around 45°N and 25°S, where most human activities take place, suggesting the potency of human influence on the climate change.

54 ENVIRONMENTAL SCIENCES↗