Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLOTTING”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

MRV Challenge 3: velocity and passive scalar comparison in a 3D turbulent flow

The third iteration of a magnetic resonance velocimetry (MRV) challenge comparison activity consisting of 3D flow measurements for a turbulent flow in a water channel moving past centrally positioned cubic obstacles is reported. In this challenge iteration, MRV measurements are coupled with either temperature or concentration measurements to extend the diagnostic utility. Five research teams from around the world conducted the measurements on a single shared apparatus. The water channel included partial elements along the channel sidewalls that varied in height and precluded easy optical accessibility. In addition to a turbulent mainstream flow, a secondary flow entered the channel between the second and third cubic elements from a square injector hole at the channel bottom wall. The injector jet interacts with the mainstream flow and mixes turbulently as it advects downstream. For the selected flow regime and water solutions used as working fluids, the mixing of the higher temperature or concentration through the secondary flow with the mainstream flow satisfies the same dynamics so that both temperature and concentration fields behave as passive scalars and can be directly compared. The measurements are explained in detail by each participating team, and the results are interpolated onto a common grid and compared using line profiles, contour plots, and isosurfaces.

Benson, Mike [ORNL] (ORCID:000000023210116X)↗
DOE OSTI2025

Soil Particulate Organic Matter is Related to Ericoid Mycorrhizal Shrubs, not Ectomycorrhizal Trees, in a Temperate Forest

Mycorrhizal associations are key drivers of soil biogeochemistry, but previous studies have focused almost exclusively on ectomycorrhizal (EcM) and arbuscular mycorrhizal (AM) associations. Ericoid mycorrhizal (ErM) shrubs frequently occur in forest understories and are expanding in response to disturbance, but are rarely considered in biogeochemical frameworks. We investigated the relationships of understory ErM shrubs and overstory trees on carbon (C) and nitrogen (N) in soil organic matter fractions in a southern Appalachian temperate forest. We sampled the 0–10 cm mineral soil layer from 43 plots at the Coweeta Hydrologic Laboratory, across gradients in overstory EcM dominance and understory ErM shrub biomass. Soil C:N ratios increased with both increasing EcM dominance and increasing ErM shrub biomass. However, total particulate organic matter (POM) C, and the proportion of C and N held in POM increased with increasing ErM shrub biomass, but not with increasing EcM dominance. In contrast, mineral-associated organic matter (MAOM) C and N were negatively associated with EcM dominance, but were not related to ErM shrubs. Our findings suggest that ErM shrubs facilitate POM formation while AM trees promote MAOM formation. Because ErM shrub biomass represents a small fraction of total forest biomass, our work provides evidence that ErM shrubs have an outsized effect on soil organic matter, which advocates for their inclusion in mechanistic studies and biogeochemical frameworks.

Bonilla, Kayla A. [University of Georgia, Athens, ↗
DOE OSTI2024

Empirical Characterization and Modeling of Cohesive – to – Adhesive Shear Fracture Mode Transition due to Increased Adhesive Layer Thicknesses of Fiber Reinforced Composite Single – Lap Joints

Here, to ensure a strong adhesive bond, most standards and adhesive manufacturers specify a maximum adhesive gap of 1 mm when bonding fiber reinforced composite structures. In manufacturing large components, such as joining two halves of wind turbine blades, meeting this gap tolerance specification is impractical; gaps larger than 10 mm are common in large adhesively bonded composite structures using state-of-the-art manufacturing techniques. Currently, there is a lack of fundamental understanding of the failure mechanics of adhesive gaps larger than 3 mm. To create such understanding, glass fiber - acrylic thermoplastic composite panels bonded using different epoxy adhesives within single-lap joint samples with adhesive thicknesses of 0.1 mm, 0.3 mm, 1 mm, 3 mm, 5 mm, and 10 mm were sheared to failure. A transition from cohesive to adhesive failure was observed to occur about 1 mm to 3 mm joint thicknesses. Plotting the shear stress normalized by the ratio of the joint width to thickness as a function of the joint thickness normalized by the joint length is shown to result in the ability to fit simple empirically derived models of the cohesive-to-adhesive failure transition, regardless of the adhesive. Furthermore, using these normalized variables, all the observed cohesively failed specimens collapse to a single master curve, as do the adhesively failed specimens.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Unlocking plant-microbial interactions in deep Mollisols in the Midwestern US: Linking depth gradients in roots, microbial activity, and soil carbon in agroecosystems

Deep-rooted plants may build soil carbon (C) stocks, but most research has focused on shallow soils, leaving gaps in our understanding of how shifts in the balance between decomposition and C inputs drive soil C accumulation with depth. Thus, our objectives were to: (1) link depth gradients in root biomass with microbial activity and soil C stocks down to 1 m, and (2) examine the potential of simple C inputs to prime soil C across depths. To this end, we dug 5 quantitative soil pits in Argiudolls under mature perennial miscanthus plots in the SoyFACE Farm (Champaign-Urbana, IL). We added 13 C labeled glucose to our soils to determine the fate of simple C inputs with depth. We found that fine root biomass, total soil C, mineral-associated organic C (MAOC), particulate organic C (POC), and microbial activity (as measured by potential enzyme activity) declined with depth. POC declined more rapidly than MAOC, resulting in an increase in the ratio of MAOC-to-POC. Root biomass, enzyme activity (either acid phosphatase or n-acetyl-glucosaminadase) activity, and microbial respiration explained 74% and 38% of the variability in soil total C and MAOC, respectively, while POC was dependent on root biomass and microbial respiration (47%). Although the incorporation of simple 13 C inputs into MAOC was similar across depths, these inputs led to greater net MAOC losses in shallow soils than in deeper soils between 50 and 100 cm. The divergent impact of simple C inputs across depths may suggest that MAOC in shallow soils is more susceptible to priming losses, while C inputs into deep soils may instead be more persistent. Collectively, our results suggest that depth gradients in soil C stocks represents a balance between inputs, decomposition, and microbial necromass production and that increases in root C inputs by deep-rooted plants may have the potential to build stable MAOC.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Variation in Biomass Yield and Cell Wall Composition in Switchgrass Natural Variants Under Two Nitrogen Regimes

Switchgrass (Panicum virgatum) is a promising lignocellulosic biofuel crop for which biomass and processing quality are important. Inherent plant variability across genotypes and environments challenges uniformity and product quality. In this study, the impact of nitrogen (N) application on switchgrass yield and quality was examined under field conditions using a highly diverse switchgrass panel over a 4-year period at Knoxville, TN. Overall, biomass production was correlated between low (0 kg of added N/ha) and moderate (135 kg of added N/ha) nitrogen treatments, suggesting that the N impact is largely uniform across the genotypes. Nonetheless, high biomass genotypes were identified with high nitrogen-use efficiency ; biomass was congruent or even higher (up to 9-fold) in the low N treatment. Genotypes were also identified with up to 94% nitrogen-remobilization efficiency. Furthermore, nitrogen application appeared to impact lignin content in whole tillers but was neutral to lignin monomer syringyl-to-guaiacyl (S/G) ratio. The same panel grown under natural conditions in Watkinsville, GA, produced significantly more biomass than the Tennessee panel in the first 3 years, but biomass was similar across both sites and treatments in year 4. Top performing genotypes overlapped between sites by 20-37%. There were low correlations in lignin content in whole tillers across the two field sites, but moderate correlations were observed for S/G ratios. The high yielding genotypes from the low N plot identified in this study can be used in breeding programs and management strategies in switchgrass to evade adverse environmental and economic effects.

09 BIOMASS FUELS↗
DOE OSTI2025

Petrology and geochemistry of the Middle Pennsylvanian (Langsettian) Clintwood coalbed, Pike County, Kentucky

Abstract The Middle Pennsylvanian, Langsettian sub-stage Clintwood coal, along with its correlatives, was one of the more important energy resources in eastern Kentucky. The coal thickens from the NE-SW-trending Belfry anticline in central Pike County to the southeast, towards the Virginia-Kentucky border. Much of the thickness increase is in the bright lithotypes below a dull lithology. The bright lithologies transition from a bright clarain + vitrain + fusain lithology to a clarain + fusain lithology to the southeast. Half of the lithotypes in the three analyzed sections exceed 1000-µg/g Rare earth elements + Y + Sc (REYSc) (ash basis). The La N /Sm N vs. Gd N /Yb N and La N /Sm N vs. Ce N /Ce N * plots suggest that the top lithotype may have had best indication of an oxidizing environment. Positive Gd N /Gd N * and Eu N /Eu N * in some of the lithologies suggests that there may have been a hydrothermal influence in the sediment sources and/or in the diagenesis of the coal. Consideration of La N /Sm N and Zr (µg/g; ash basis), Zr vs. Sr, and V/Cr vs. inertinite/(inertinite + vitrinite) indicates that the dull lithotypes existed in a distinct geochemical environment compared to the brighter lithotypes. The isolation of the dull lithotypes from the bright lithotypes is supported by principal components analysis on inertinite/(inertinite + vitrinite), La N /Sm N , V/Cr, Sr, and either Zr or Ln(100(TiO 2 )/Al 2 O 3 ).

Hower, James C. (ORCID:0000000346942776)↗
DOE OSTI2025

Thermodynamic assessment of the quaternary WTaCrV refractory high entropy alloy as a means to guide experimental approaches

The deployment of fusion energy poses challenges for materials in plasma facing components to withstand high temperatures and thermal gradients, particle implantation and neutron damage. The current material of choice is tungsten, although property degradation limits its consideration in future fusion reactors. Hence, materials with better resistance to harsh environments need to be developed for fusion energy to become a reality. High entropy alloys are being explored as potential candidates with some compositions showing good radiation resistance to defect cluster formation. One of these materials is the WTaCrV system, although only one composition has been tested under ion irradiation. In this work, we study the thermodynamic properties of the entire quaternary alloy composition range. Coupling first principles calculations, cluster expansion approaches, and Monte Carlo methods, we access the free energy functionals, short-range ordering as a function of temperature, and atomic configurations that can be compared to experimental observations. We use this data to inform experiments into compositions with higher propensity to form solid solutions, instead of phase separating. With this formalism we have developed thermodynamic database (TDB) files that can be used to plot quaternary phase diagrams.

Cluster Expansion↗
DOE OSTI2025

Quantifying soil organic matter stock distribution and origin following over a century of maize-based cropping in the former tallgrass prairie region of central USA

Tallgrass prairie conversion to maize-based agriculture in central North America has resulted in substantial loss of soil organic carbon (SOC) in less than two centuries. However, evaluations of how management practices may mitigate SOC losses are generally limited in soil depth and/or duration, missing long-term SOC stock outcomes that manifest over timescales of decades or longer. To address this, we sampled soils in year 145 of the Morrow Plots experiment to (i) evaluate effects of crop rotation and fertility management on SOC stocks and (ii) distinguish prairie- versus maize-derived SOC after continuous maize cropping since 1876 using stable carbon isotope ( 13 C) natural abundance. Soil organic carbon stock by equivalent soil mass (ESM) was + 30.7 Mg C ha −1 (+31.7 %) higher under maize-oat-alfalfa than continuous maize, but similar between maize-soybean and continuous maize. NPK fertilization and manuring did not influence SOC stocks by ESM. Response of SOC stocks at 15 cm depth intervals to NPK fertilization varied by depth and crop rotation, with lower SOC stocks at 30–45 cm under continuous maize and maize-soybean. Maize-derived C ranged 19.5–59.6 % of SOC stock across depths, indicating the majority of SOC was still derived from tallgrass prairie even after 145 years of continuous maize cropping. Our results confirm the potential of diversified crop rotation for minimizing SOC losses relative to tallgrass prairie at the supracentennial scale, and highlight the importance of relic prairie soil organic matter for future crop production in central North America.

crop rotation↗
DOE OSTI2026

The productivity of perennial crops miscanthus and switchgrass is more resistant to vapor pressure deficit stress than maize

Rising atmospheric vapor pressure deficit (D) with warming is an increasingly important driver of crop productivity loss, yet the relative sensitivity of conventional seed crops and alternative perennial crops remains poorly resolved. In this study, we combined a 9-year eddy covariance carbon flux record (2008 – 2016) with 600+ midday leaf water potential (ψ L ) measurements (2024 – 2025) from adjacent maize (Zea mays), miscanthus (Miscanthus x giganteus), and switchgrass (Panicum virgatum) plots to: (1) quantify the limitations of elevated D on crop-specific productivity and (2) evaluate how these responses vary as a function of soil moisture (Θ) status. We found that maize strictly regulated ψ L and its gross primary productivity (GPP) was strongly reduced by increasing D. In contrast, miscanthus and switchgrass allowed larger ψ L declines and sustained higher GPP under similar moisture constraints. D exerted a larger limitation on maize GPP than Θ, but D-driven declines were less severe when accompanied by high Θ. GPP sensitivity to D was also influenced by Θ for miscanthus and switchgrass, but the buffering effect was stronger compared to maize. Specifically, high D paired with high Θ were the most productive conditions for the perennials, reflecting temperature-driven gains in photosynthetic efficiency when evaporative stress was mitigated by sufficient soil water supply. Overall, miscanthus and switchgrass displayed greater resistance to atmospheric drought and more stable GPP across hydroclimate variability than maize. These findings highlight that maize is more vulnerable to future projections of rising D than miscanthus or switchgrass.

drought↗
DOE OSTI2026

Spaced out: An economic framework to explore the impacts of PV panel spacing on large-scale farming in Colorado

CONTEXT Agrivoltaic systems co-locate solar technologies with agricultural operations on an integrated plot of land and potentially provide benefits to both energy and agricultural systems. To date, large-scale (>5-MW) agrivoltaic projects in the United States have been limited to grazing and ecovoltaic applications, raising questions about the impact and scalability of agrivoltaic crop systems. Many agrivoltaic designs raise the height of the solar panels to accommodate agricultural practices while keeping energy density high. However, raising the panels results in increased photovoltaic (PV) development costs, which often are higher than the economic returns of crop production underneath the panels. This leads to unfavorable project economics and the need for other agrivoltaic solutions than raising panels. OBJECTIVE To explore other solutions, we perform an initial feasibility analysis for an agrivoltaic solution that can integrate with large-scale farming practices by increasing the row spacing in between panels. Increased PV row spacing is a low-cost approach for scaling agrivoltaics to accommodate crop production and this spacing can be tailored to required crop equipment for different regions. Increasing row spacing will reduce the power density (PV installed per acre), but in areas that are not land limited, these agrivoltaic designs could be economically feasible. Our analysis establishes a framework for a feasibility analysis for where and with what crops spaced out panel agrivoltaic solutions might be economical. METHODS Using a case study for large-scale agriculture crops in Colorado, we establish a framework for wide-row agrivoltaic economic feasibility analysis. We utilized the System Advisor Model to calculate technoeconomic metrics to compare different row spacing solutions and capture tradeoffs of these system designs. RESULTS AND CONCLUSIONS We find that, in some circumstances, wider row agrivoltaic solutions that allow for continued mechanized crop production can provide economic benefits over a traditional utility-scale PV system. For most crops examined in this analysis, roughly $\$$200/acre in agricultural profit justified spacing out the panels to at least 31.7 ft. to accommodate agrivoltaic configurations versus PV only configurations. Additionally, opportunities for increased agricultural revenue with agrivoltaic systems allow PV project economics to tolerate a larger range of CAPEX variability while remaining economically viable relative to the PV only configurations. SIGNIFICANCE This framework can be adapted for a wide variety of crops and regions and allows for examination of economically favorable sites for future agrivoltaic systems that utilize different configuration and expand opportunities for agrivoltaics.

14 SOLAR ENERGY↗
DOE OSTI2025

Transfer learning-based soybean LAI estimations by integrating PROSAIL, UAV, and PlanetScope imagery

Accurate Leaf Area Index (LAI) estimations at the soybean plot scale is achievable using high-resolution Unmanned Aerial Vehicle (UAV) imagery and field measurement samples. However, the limited coverage of UAV flights restricts large-scale remote sensing monitoring in expansive soybean fields. This study leverages the broad coverage and 3-m resolution of PlanetScope satellite imagery to extend LAI prediction from UAV to satellite scales through transfer learning, using UAV-scale LAI estimates as a benchmark to validate cross-scale consistency. To address this challenge, this study proposed the LAI-TransNet, a two-stage transfer learning framework designed for precise and scalable soybean LAI prediction across large areas, demonstrating its effectiveness in cross-scale monitoring. In Stage 1, a UAV-scale benchmark is established using PROSAIL-simulated UAV reflectance data (UAV-Sim) and field-measured soybean LAI. Traditional machine learning, deep learning, and transfer learning models are trained on a hybrid UAV-Sim and field-measured dataset (UAV-Sim_Measured), with the transfer learning model CNN-TL, fine-tuned using pre-trained weights derived from UAV-Sim, achieving the highest accuracy (R 2 = 0.81, RMSE = 0.64 m 2 /m 2 , rRMSE = 11.5 %). In Stage 2, LAI-TransNet is developed by fine-tuning the CNN-TL model on PlanetScope simulated data (PS-Sim), preprocessed via cross-domain mapping to align UAV and satellite spectral features. Real PlanetScope imagery is corrected for reflectance consistency with reference to UAV imagery spectral profiles. LAI-TransNet outperforms other deep learning models trained directly on PS-Sim (R 2 = 0.69 vs. 0.60–0.63), ensuring robust cross-scale consistency. In conclusion, by bridging UAV and satellite scales, LAI-TransNet enables large-scale soybean LAI monitoring, enhancing precision agriculture management through improved monitoring with the PlanetScope imagery.

Leaf area index (LAI)↗
DOE OSTI2025

SAXS Assistant: Automated SAXS analysis for structural discovery in biologics and polymeric nanoparticles

Small-angle x-ray scattering (SAXS) is a powerful technique for assessing macromolecular structure. High-throughput SAXS is limited by the time-consuming and, at times, subjective nature of SAXS data interpretation. Here, we present SAXS Assistant, a Python-based script that streamlines SAXS data analysis to extract features for machine learning (ML) and key structural parameters, including the Guinier radius of gyration (R g ), pair distance distribution function (PDDF)-derived R g , maximum particle dimension (D max ), and Kratky plots. The script builds upon BioXTAS RAW and validates reliability via Guinier/PDDF R g agreement, an important indicator of well-measured data sets. For assistance in D max estimation, a multilayer perceptron regressor was trained with 1940 data files from the Small Angle Scattering Biological Data Bank. The model achieved a test set performance R 2 = 0.90 and mean absolute error = 11.7 Å. Training exclusively with experimental data translates analyses from researchers, including experts in the field, to the ML model, which helps assess D max estimations from PDDF. Gaussian mixture model clustering was implemented to classify profiles into structural classes based on entries in the Small Angle Scattering Biological Data Bank. Users may therefore assess the similarity between experimental samples and known biomolecular shapes within the mapped repository entries. This probabilistic clustering aids in quantifying information from Kratky and generating shape-descriptive features. SAXS Assistant accelerates SAXS data analysis through enforced quality control, ML-ready outputs, and flags for low-confidence results. In addition to providing the ability to analyze large data sets at high throughput, this tool is versatile and may serve researchers in both biological and synthetic polymer research fields.

36 MATERIALS SCIENCE↗
DOE OSTI2025

SynopFrame: Multiscale time-dependent visual abstraction framework for analyzing DNA nanotechnology simulations

We present an open-source framework, SynopFrame, that allows DNA nanotechnology (DNA-nano) experts to analyze and understand molecular dynamics simulation trajectories of their designs. We use a multiscale multi-dimensional abstraction space, connect the representations to a projected conformational space plot of the structure’s temporal sequence, and thus enable experts to analyze the dynamics of their structural designs and, specifically, failure cases of the assembly. In addition, our time-dependent abstraction representation allows the biologists, for the first time in a smooth and structurally clear way, to identify and observe temporal transitions of a DNA-nano design from one configuration to another, and to highlight important periods of the simulation for further analysis. We realize SynopFrame as a dashboard of the different synchronized 3D spatial and 2D schematic visual representations, with a color overlay to show essential properties such as the status of hydrogen bonds. The linking of the spatial, schematic, and abstract views ensures that users can effectively analyze the high-frequency motion. We also categorize the status of the hydrogen bonds into a new format to allow us to color-encode it and overlay it on the representations. To demonstrate the utility of SynopFrame, we describe example usage scenarios and report user feedback.

Abstraction space↗
DOE OSTI2026

Extending SLUSCHI for Automated Diffusion Calculations

We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Catalytic resonance theory for the kinetics of photon-promoted catalysis

The illumination of catalytic surfaces with a continuous or pulsed stream of photons dynamically modulates surface chemistry for faster rates, non-equilibrium conversion, or product selectivity control. To establish fundamental principles of dynamic photon-modulated catalysis, the photocatalytic conversion of a generic surface reaction was simulated using the kinetic Monte Carlo method to understand the kinetic implications of an independent stream of photons that promotes surface product desorption. The time-averaged photocatalytic rate at differential conditions for varying photon flux and temperatures indicated three kinetic regimes described by product thermal desorption control, surface reaction control, and an intermediate kinetic regime with a zero slope Arrhenius plot, consistent with a degree of rate control dominated by the photon arrival frequency (i.e., per-site photon flux). Here, the maximum photocatalytic rate occurred orders of magnitude above the Sabatier limit at the resonance frequency, identified as the photon arrival frequency matching the surface reaction rate constant.

Canavan, Jesse R. [University of Minnesota, Minnea↗
DOE OSTI2024

HTESP (High-throughput electronic structure package): A package for high-throughput ab initio calculations

High-throughput ab initio calculations are the indispensable parts of data-driven discovery of new materials with desirable properties, as reflected in the establishment of several online material databases. The accumulation of extensive theoretical data through computations enables data-driven discovery by constructing machine learning and artificial intelligence models to predict novel compounds and forecast their properties. Efficient usage and extraction of data from these existing online material databases can accelerate the next stage materials discovery that targets different and more advanced properties, such as electron–phonon coupling for phonon-mediated superconductivity. However, extracting data from these databases, generating tailored input files for different ab initio calculations, performing such calculations, and analyzing new results can be demanding tasks. Here, in this work, we introduce a software package named “HTESP” (High-Throughput Electronic Structure Package) written in Python and Bash languages, which automates the entire workflow including data extraction, input file generation, calculation submission, result collection and plotting. Our HTESP will help speed up future computational materials discovery processes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Phase transformation mechanism in irradiation-induced superlattice formation

Atomic kinetic Monte Carlo simulations were used to model void superlattice formation under irradiation in molybdenum, driven by anisotropic diffusion of self-interstitial atoms. A change in the phase transformation mechanism from nucleation and growth to spinodal decomposition occurred with increasing dose rate, with both mechanisms leading to superlattice formation. Analysis of a rate-theory based analytical model showed that an observed change in the kinetics of vacancy accumulation, the appearance of a region of positive second derivative in the plot of average vacancy concentration versus time, was caused by the onset of spinodal instability. Further, the analytical model showed that for molybdenum and several other metals where void superlattice formation is commonly observed, the phase transformation likely occurs by nucleation and growth. However, nickel may offer the possibility of experimental observation of the transition between phase transformation mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Parallel derivative-free optimization for simulation-based design of behind-the-meter energy systems

In this work, the integrated design and dispatch of behind-the-meter or distributed resources (e.g. stationary battery storage and solar PV generation) is considered. A simulation-based framework is employed, generating high-fidelity results with closed-loop predictive control at a fine resolution, at the expense of high computational cost (several minutes to a few hours per design point). To address this challenge, parallel derivative-free design methods are considered. Four methods are compared, including state-of-the-art surrogate-based methods (Radial-Basis Functions and Gaussian processes) and sampling strategies, an evolutionary-based method, and a simple sequential grid refinement method. As a case study, two types of design problem with increasing complexity are considered, namely, the design of behind-the-meter resources (three design variables) and the inclusion of grid capacity (four design variables). The second yields a constrained design problem for which violations can only be determined after solving the computationally expensive simulation. For the three-dimensional case, all methods present a good performance, achieving a solution within 1% of the optimum after the first iteration, with the sequential grid refinement exhibiting the fastest convergence and achieving the best final objective value. This indicates that the parallel evaluation of multiple sampling points may be more important than the choice of method for small decision spaces. For the four-dimensional constrained case, the Genetic Algorithm presents the best tradeoff between performance and computational effort, while the rough objective function terrain generated by constraint violation penalties reduces the performance of surrogate-based methods. Contour plots with flat regions indicate flexibility in the optimal design and highlight the importance of characterizing the solution space.

24 POWER TRANSMISSION AND DISTRIBUTION↗