Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model Predictive Control”

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 1,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Data-driven Mori–Zwanzig modeling of Lagrangian particle dynamics in turbulent flows

The dynamics of Lagrangian particles in turbulence play a crucial role in mixing, transport, and dispersion in complex flows. Their trajectories exhibit highly nontrivial statistical behavior, motivating the development of surrogate models that can reproduce these trajectories without incurring the high computational cost of direct numerical simulations of the full Eulerian field. This task is particularly challenging because reduced-order models typically lack access to the full set of interactions with the underlying turbulent field. Novel data-driven machine learning techniques can be powerful in capturing and reproducing complex statistics of the reduced-order/surrogate dynamics. In this work, we show how one can learn a surrogate dynamical system that is able to evolve a turbulent Lagrangian trajectory in a way that is point-wise accurate for short-time predictions (with respect to Kolmogorov time) and stable and statistically accurate at long times. This approach is based on the Mori–Zwanzig formalism, which prescribes a mathematical decomposition of the full dynamical system into resolved dynamics that depend on the current state and the past history of a reduced set of observables, and the unresolved orthogonal dynamics due to unresolved degrees of freedom of the initial state. We show how by training this reduced order model on a point-wise error metric on short time-prediction, we are able to correctly learn the dynamics of Lagrangian turbulence, such that also the long-time statistical behavior is stably recovered at test time. This opens up a range of applications, for example, for the control of active Lagrangian agents in turbulence.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Low-Cost Heliostat for High-Flux Small-Area Receivers (Final Technical Report)

This project analyzed a two-stage heliostat concept consisting of a tracking stage and a concentrating stage. The tracking stage uses mirrors mounted on a common drive that move to track the sun. The concentrating stage consists of stationary mirrors that each have a unique angle to direct rays towards a small-area, high-flux, point-focused receiver. By splitting the collection and concentrating process into two stages, multiple small, inexpensive mirrors can share a structure and be controlled by a single drive in the tracking stage. The project effort developed modeling techniques that were specifically relevant to this two-stage heliostat concept. Both field-level and unit-level models were developed. The field-level model does not explicitly consider unit-level losses which are predicted by the unit-level model and then integrated into the field-level model through a correlation referred to as an efficiency modifier. This approach is referred to as the two-model approach; the development and demonstration of this two-model approach for a multi-stage heliostat technology is a key outcome of this work. The field-level model is used to design a field that hits a specific design day power given a set of heliostat design parameters. An oversized field is simulated and then heliostat units are removed based on their annual energy production in order to generate the highest performing field. The field reduction procedure fits a smooth curve fit to annual energy production as a function of position in the field which has the effect of reducing the noise that is otherwise caused by the Monte Carlo ray tracing technique. This approach is referred to as the annual energy fit method and substantially reduces computational run time for a given field level modeling accuracy. The annual energy fit approach enables the selection of a properly sized, high-performing field using orders of magnitude fewer rays than would otherwise be possible and the development of this approach is a second key outcome of this work. These models are used within a genetic optimization algorithm in order to optimize the geometric parameters associated with a heliostat in order to achieve the lowest cost per unit of collected design day power. The cost modeling that underlies the optimization is a simple, scaling type analysis backed up by a much more detailed Design for Manufacture and Assembly (DFMA) analysis. Although the figure of merit used for optimization was not cost per mirror area, this metric is reasonable to use as a means of comparison. The optimally designed 500 kW design has a tracking mirror specific cost of $181.85/m 2 , which is significantly larger than the target value and also larger than the current state of the art. The cost of the torque-tube type linkages contributed substantially to the overall cost. Based on this observation, potentially attractive alternative design configuration utilizing a capstan type actuation system should be investigated. Finally, NREL compared the performance of the two-stage heliostat to the performance of a focused and different sized flat conventional heliostats and showed that, as expected, additional losses versus the convention heliostat caused by a worse cosine efficiency, two stages of reflection, and interstage interactions. The two-stage heliostat requires around 75% more reflective area than a flat 1x1 meter conventional heliostat (similar to a focused heliostat) and 40% more than a flat 2x2 meter conventional heliostat.

14 SOLAR ENERGY↗
DOE OSTI2025

On the emergent scale of bedrock groundwater contribution to headwater mountain streams

We investigated the contribution of bedrock groundwater to streamflow as a function of catchment scale in a headwater stream. Synoptic surveys were conducted during hydrologically important periods of the year using multiple environmental tracers in stream water, soil water, and bedrock groundwater, along a first-order montane stream, in west-central Montana. Sampled analytes included 222 Rn, used to constrain total subsurface flux, and major and minor elements, used in end-member mixing analysis (EMMA) to identify the contributions of soil and bedrock groundwater to the stream. Partitioning between soil-derived and bedrock-derived groundwater was then analyzed as a function of the incremental and accumulated sub-catchment sizes. Radon results indicated that subsurface water contributions accounted for the majority of streamflow at all surveyed times. EMMA results revealed that the bedrock groundwater contribution to streamflow varied between 26% during peak snowmelt and 44% during late summer. Streamflow generation was dominated by soil groundwater contribution along the entire reach, but the bedrock groundwater contribution increased consistently with accumulated sub-catchment size. However, groundwater contributions were not well-correlated with incremental sub-catchment size. The scale at which increased bedrock groundwater discharge can be correlated with sub-catchment size appears to be >1 km 2 for our study. Our results are consistent with a conceptual model where streamflow is predominantly generated by a 3D subsurface nested flow system. Local subsurface heterogeneities control the stream source at local scales but begin to average out at scales >2 km 2 . Our study indicates that, while soil groundwater is the dominant source, bedrock groundwater remains an important and predictable contributor to streamflow throughout the year, even in a snow-dominated, mountainous headwater catchment.

environmental tracers↗
DOE OSTI2020

Power system model parameter conditioning tool

A power system model parameter conditioning tool including a server control processor in communication with phasor measurement unit monitored data records of multiple disturbance events, a model calibration unit providing event screening, power system model simulation, and simultaneous tuning of model parameters. The model calibration performing a simulation using default model parameters, the processor comparing the simulation results to the monitored data. If the prediction is within threshold, then terminating conditioning; else performing parameter identifiability analysis to determine differing effects of various model parameters on power system model accuracy, selecting a parameter set causing a degradation in power system model prediction, and updating the default model parameters corresponding to members of the parameter set with values selected to reduce the degradation. A method and a non-transitory computer readable medium are also disclosed.

Baone, Chaitanya Ashok↗
DOE OSTI2021

Prediction of DIII-D Pedestal Structure From Externally Controllable Parameters

The sharp increase of pressure at the edge of a high confinement mode (H-mode) plasma, the pedestal, strongly impacts overall plasma performance. Predicting the pedestal is a necessity to control and optimize tokamak operations. Here, an experimental data-driven machine learning (ML) approach is presented that predicts the pedestal heights and widths of electron density (n e ) and electron temperature (T e ) profiles as well as the separatrix ne from externally controllable parameters such as the plasma shape, heating method and power, and gas puff rate and integrated gas puff. The OMFIT framework was used with DIII-D data to efficiently, robustly, and automatically build a database of pedestal parameters to train machine learning models. Database creation was enabled by the search engine tool for DIII-D data, TokSearch, which parallelizes data fetching, enabling fast searches through basic signals of thousands of DIII-D shots and selection of relevant time intervals. Principal Component Analysis (PCA) separated the database into three clusters that represent classes of plasma shapes that are regularly used in DIII-D. The most important parameters for setting the pedestal structure were plasma current (I p ), toroidal magnetic field (B Φ ), neutral beam heating power (P NBI ) and shaping quantities. The Deep Jointly Informed Neural Networks (DJINN) algorithm was applied to identify suitable neural network (NN) architectures that appropriately capture the features of the pedestal database. Separate NNs were implemented for each pedestal parameter, and ensembling methods were used to improve the prediction accuracy and allowed estimation of the prediction uncertainty. The pedestal predictions of the test dataset lie within the measurement uncertainties of the pedestal parameters. The NN outperformed simple Linear Regression (LR) analysis, indicating non-linear dependencies in the pedestal structure. The presented achievements illustrate a promising path for future research, using feature extraction to infer experimental trends and thereby improve pedestal models as well as deploying NN for a fast pedestal prediction in DIII-D scenario development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Three-Dimensional CFD Investigation of Pre-Spark Heat Release in a Boosted SI Engine

Low-temperature heat release (LTHR) in spark-ignited internal combustion engines is a critical step toward the occurrence of auto-ignition, which can lead to an undesirable phenomenon known as engine knock. As such, correct predictions of LTHR are of utmost importance to improve the understanding of knock and enable techniques aimed at controlling it. While LTHR is typically obscured by the deflagration following the spark ignition, extremely late ignition timings can lead to LTHR occurrence prior to the spark, i.e., pre-spark heat release (PSHR). In this research, PSHR in a boosted direct-injection SI engine was numerically investigated using three-dimensional computational fluid dynamics (CFD). A hybrid approach was used, based on the G-equation model for representing the turbulent flame front and the multi-zone well-stirred reactor model for tracking the chemical reactions within the unburnt region. A recently developed best practice was also employed which keeps the well-stirred reactor model active throughout the entire simulation. This allowed for correct predictions of the previous cycle trapped residuals which have a considerable effect on the onset of LTHR. Multi-cycle simulations were conducted using Co-Optima alkylate and E30 fuels. The predicted in-cylinder pressure and heat release rate agreed well with the experimental data and served as validation for the CFD model. Following the initial validation, the dynamics of PSHR was discussed and a series of parameters of interest were assessed. First, the effect of exhaust valve temperature was investigated, qualitatively highlighting the importance of boundary conditions uncertainty. Further analyses were carried out on the effects of fuel properties, including laminar flame speed (LFS) and heat of vaporization (HOV). The results indicate that PSHR phasing is slightly advanced with lower LFS as more trapped unburnt fuel is made available for the next cycle. A similar trend was observed with lower HOV as less intense spray cooling led to higher mixture temperatures, i.e., higher mixture reactivity. Finally, a comparison of Co-Optima alkylate and E30 fuels was made using the pressure-temperature trajectory framework. It was shown that the differences between Co-Optima alkylate and E30 in PSHR tendency are correlated with both HOV and chemical effects.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Synthesis of results for Brine Availability Test in Salt (BATS) DECOVALEX-2023 Task E

In the 2023 phase of the international collaborative DECOVALEX modeling project, Task E focused on understanding thermal, hydrological, and mechanical (THM) processes related to predicting brine migration in the excavation damaged zone around a heated excavation in salt. Salt is attractive as a disposal medium for radioactive waste because it is self-healing and is essentially impermeable and non-porous in the far field. Investigation of the short-term, near-field behavior is important for radioactive waste disposal because this early period strongly controls the amount of inflowing brine. Brine leads to corrosion of waste forms and waste packages, and possible dissolution of radionuclides with brine transport being a potential transport vector to the accessible environment.The Task was divided into steps. Step 0 included matching unheated brine inflow data from boreholes at the Waste Isolation Pilot Plant (WIPP) and matching temperature observations during a Brine Availability Test in Salt (BATS) heater test. Step 1 included validation of models against a thermo-poroelastic analytical solution, and two-phase flow around an excavation. Finally, Step 2 required all the individual components covered in steps 0 and 1 to come together to match observed brine inflow behavior during the same BATS heater test.There were a range of approaches from the teams, from mechanistic to prescriptive. Given the uncertainties in the problem, some teams used one- or two-dimensional models of the processes, while other teams included more geometrical complexity in three-dimensional models. Task E was a learning experience for the teams involved, and feedback from the modeling teams has led to changes in follow-on BATS experiments at WIPP. The primary Task E lessons learned were the impact of hydrologic initialization methods (wetting up vs. drying down), the difference between confined and unconfined thermal expansion, and the large changes in permeability associated with heating and cooling.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Probing Excited-State Dynamics of Transmon Ionization

The fidelity and quantum nondemolition character of the dispersive readout in circuit QED are limited by unwanted transitions to highly excited states at specific photon numbers in the readout resonator. This observation can be explained by multiphoton resonances between computational states and highly excited states in strongly driven nonlinear systems, analogous to multiphoton ionization in atoms and molecules. In this work, we utilize the multilevel nature of high-𝐸 𝐽 /𝐸 𝐶 transmons to probe the excited-state dynamics induced by strong drives during readout. With up to ten resolvable states, we quantify the critical photon number of ionization, the resulting state after ionization, and the fraction of the population transferred to highly excited states. Moreover, using pulse shaping to control the photon number in the readout resonator in the high-power regime, we tune the adiabaticity of the transition and verify that transmon ionization is a Landau-Zener-type transition. We further extend these methods to a typical transmon with 𝐸 𝐽 /𝐸 𝐶 ≈ 55 and probe the offset-charge dependence of ionization dynamics in a timed-resolved manner. Our experimental results agree well with the theoretical prediction from a semiclassical driven transmon model and may guide future exploration of strongly driven nonlinear oscillators.

cavity quantum electrodynamics↗
DOE OSTI2026

FAIR Data and Interpretable AI Framework for Architectured Metamaterials (Final Report)

This research program established a transformative framework for the discovery and design of mechanical metamaterials, which are architected structures engineered to control physical phenomena like sound and vibration in ways natural materials cannot. To overcome the traditional reliance on trial-and-error, the project developed an interpretable Artificial Intelligence (AI) framework that moves beyond "black box" models to reveal the specific geometric patterns—such as "unit-cell templates"—that govern a material’s performance. A major breakthrough was the development of a hierarchical design method, which allows a single material to block vibrations across multiple frequency ranges simultaneously by layering patterns at different scales without them interfering with one another. This was further expanded to include irregular, graph-based designs that use spanning tree algorithms to ensure structural connectivity while allowing for customized, direction-dependent properties like stiffness and acoustic impedance. Beyond design, the project addressed the practicalities of real-world production by developing uncertainty quantification techniques that account for manufacturing defects and material variability, reducing the need for expensive physical testing by orders of magnitude. To speed up the discovery process, the team implemented Gaussian Process Regression and other surrogate models that provide accurate performance predictions at a fraction of the traditional computational cost. The AI-generated designs were successfully validated through fabrication of physical samples and wave propagation experiments, confirming their ability to accurately guide or reflect waves as predicted. By contributing these tools and high-quality FAIR benchmark datasets to the wider scientific community, this work provides a scalable foundation for advancing technologies in aerospace vibration control, medical imaging, and noise reduction.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Evaluating the WFIP2 updates to the HRRR model using scanning Doppler lidar measurements in the complex terrain of the Columbia River Basin

The wind-energy (WE) industry relies on numerical weather prediction (NWP) forecast models as foundational or base models for many purposes, including wind-resource assessment and wind-power forecasting. During the Second Wind Forecast Improvement Project (WFIP2) in the Columbia River Basin of Oregon and Washington, a significant effort was made to improve NWP forecasts through focused model development, to include experimental refinements to the High Resolution Rapid Refresh (HRRR) model physics and horizontal grid spacing. In this study, the performance of an experimental version of HRRR that includes these refinements is tested against a control version, which corresponds to that of the operational HRRR run by National Oceanic and Atmospheric Administration/National Centers for Environmental Protection at the outset of WFIP2. Furthermore, the effects of horizontal grid resolution were also tested by comparing wind forecasts from the HRRR (with 3-km grid spacing) with those from a finer-resolution HRRR nest with 750-m grid spacing. Model forecasts are validated against accurate wind-profile measurements by three scanning, pulsed Doppler lidars at sites separated by a total distance of 71 km. Model skill and improvements in model skill, attributable to physics refinements and improved horizontal grid resolution, varied by season, by site, and during periods of atmospheric phenomena relevant to WE. In general, model errors were the largest below 150 m above ground level (AGL). Experimental HRRR refinements tended to reduce the mean absolute error (MAE) and other error metrics for many conditions, but degradation in skill (increased MAE) was noted below 150 m AGL at the two lowest-elevation sites at night. Finer resolution was found to produce the most significant reductions in the error metrics.

17 WIND ENERGY↗
DOE OSTI2022

Progress from ASDEX Upgrade experiments in preparing the physics basis of ITER operation and DEMO scenario development

An overview of recent results obtained at the tokamak ASDEX Upgrade (AUG) is given. A work flow for predictive profile modelling of AUG discharges was established which is able to reproduce experimental H-mode plasma profiles based on engineering parameters only. In the plasma center, theoretical predictions on plasma current redistribution by a dynamo effect were confirmed experimentally. For core transport, the stabilizing effect of fast ion distributions on turbulent transport is shown to be important to explain the core isotope effect and improves the description of hollow low-Z impurity profiles. The L–H power threshold of hydrogen plasmas is not affected by small helium admixtures and it increases continuously from the deuterium to the hydrogen level when the hydrogen concentration is raised from 0 to 100%. One focus of recent campaigns was the search for a fusion relevant integrated plasma scenario without large edge localised modes (ELMs). Results from six different ELM-free confinement regimes are compared with respect to reactor relevance: ELM suppression by magnetic perturbation coils could be attributed to toroidally asymmetric turbulent fluctuations in the vicinity of the separatrix. Stable improved confinement mode plasma phases with a detached inner divertor were obtained using a feedback control of the plasma β. The enhanced D α H-mode regime was extended to higher heating power by feedback controlled radiative cooling with argon. The quasi-coherent exhaust regime was developed into an integrated scenario at high heating power and energy confinement, with a detached divertor and without large ELMs. Small ELMs close to the separatrix lead to peeling-ballooning stability and quasi continuous power exhaust. Helium beam density fluctuation measurements confirm that transport close to the separatrix is important to achieve the different ELM-free regimes. Based on separatrix plasma parameters and interchange-drift-Alfvén turbulence, an analytic model was derived that reproduces the experimentally found important operational boundaries of the density limit and between L- and H-mode confinement. Feedback control for the X-point radiator (XPR) position was established as an important element for divertor detachment control. Stable and detached ELM-free phases with H-mode confinement quality were obtained when the XPR was moved 10 cm above the X-point. Investigations of the plasma in the future flexible snow-flake divertor of AUG by means of first SOLPS-ITER simulations with drifts activated predict beneficial detachment properties and the activation of an additional strike point by the drifts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Automated ICRF heating surrogate modeling via machine learning

This work introduces automated machine learning workflows that address critical bottlenecks in surrogate model development for Ion Cyclotron Range of Frequencies (ICRF) heating applications. The automated framework includes data analysis tools that transform raw datasets into actionable insights in seconds, replacing weeks of manual exploratory effort and ensuring consistent, reproducible dataset characterization. By integrating advanced hyperparameter optimization (HPO) methods including Bayesian optimization via BoTorch and Tree-structured Parzen Estimators (TPE), the framework significantly reduces model development time from weeks to hours, decreasing computational cost and required expertise, while enabling high-accuracy surrogate models. Compared to traditional hyperparameter scanning (HPS) techniques such as methodical, randomized, and grid searches, HPO methods achieve superior convergence and predictive performance, even when compared to already well-tuned reference models. On NSTX High Harmonic Fast Wave (HHFW) heating datasets, both Random Forest Regressor (RFR) and neural network surrogates demonstrate improved accuracy, achieving R 2 values beyond 0.97 and 0.98, respectively. The results show that while HPO gains are modest for robust architectures like RFR, they become essential for more sensitive models such as neural networks, highlighting the trade-offs across optimization strategies. Through automated workflows that eliminate manual hyperparameter tuning and require minimal ML expertise, this work enables widespread adoption of high-fidelity surrogate models across the fusion community for real-time plasma control, uncertainty quantification, rapid experimental scenario development, and integrated system optimization.

Sanchez-Villar, Alvaro [Princeton Plasma Physics L↗
DOE OSTI2023

The physics of desiccation cracks 1: Ductile fracturing and dependence on relative humidity

Ductile deformation is ubiquitously found in the shrinkage of geomaterials. The existence of ductility requires elastoplastic mechanics when analyzing the structure deformation under external stress. Here, in this work, we explore the physics of ductile fracturing based on the results from desiccation experiments and triaxial tests. By using the digital image correlation (DIC) method to generate the strain maps of samples undergoing desiccation cracking under different relative humidities, we obtain results showing the previously postulated Cnoidal Wave ductile failure patterns propagating under atmospheric condition-controlled crack velocities. We then correlate these observations with rate-dependent plasticity models calibrated through triaxial tests undergoing several unloading–reloading cycles and velocity stepping. This work demonstrates the necessity to consider ductility in soil cracking, indicating that the formation of crack patterns in soil desiccation is a slow and predictable process.

58 GEOSCIENCES↗
DOE OSTI2025

Environmental Controls on Water Vapor Deuterium Excess in the Coastal Boundary Layer: An Information Theory Perspective

We use information theory to quantify the environmental controls on water vapor deuterium excess (D-excess) in coastal Southern California from June 2023 through February 2024. Using Shannon entropy, mutual information (MI), and joint mutual information, metrics that capture both linear and nonlinear relationships, we identify the most informative variables and variable combinations governing D-excess across contrasting marine and continental regimes. Relative humidity with respect to sea surface temperature (RHS) is consistently the strongest individual predictor, explaining up to 27% of D-excess variability during marine conditions but only 10% in continental air masses. The Relative humidity(RHS) + sea surface temperature (SST) combination demonstrates synergistic effects, where their joint influence (explaining up to 36% of D-excess variability) exceeds what either variable achieves individually, confirming their coupled influence on deuterium excess. Wind direction complements RHS most effectively during continental conditions. The best three-variable combination (RHS + SST + Planetary Boundary Layer height) explains 38% of D-excess variability in marine air, while no combination exceeds 20% explanatory power during continental periods. Information theory shows that heteroscedasticity in D-excess relationships indicates regime shifts in controlling processes and quantifies fundamental constraints on predictor variables: some environmental factors like surface pressure or water vapor flux contain insufficient information content to explain D-excess variability regardless of their physical relevance. These results highlight the different predictability limits between marine and continental regimes, challenging the adequacy of linear models and providing a rigorous framework for quantifying the information content of isotope-climate relationships with implications for both modern and paleoclimate applications.

information theory↗
DOE OSTI2024

Defect control strategies for Al 1- x Gd x N alloys

Tetrahedrally bonded III-N and related alloys are useful for a wide range of applications from optoelectronics to dielectric electromechanics. Heterostructural AlN-based alloys offer unique properties for piezoelectrics, ferroelectrics, and other emerging applications. Atomic-scale point defects and impurities can strongly affect the functional properties of materials, and therefore, it is crucial to understand the nature of these defects and the mechanisms through which their concentrations may be controlled in AlN-based alloys. In this study, we employ density functional theory with alloy modeling and point defect calculations to investigate native point defects and unintentional impurities in Al 1-x Gd x N alloys. Among the native defects that introduce deep midgap states, nitrogen vacancies (V N ⁠) are predicted to be in the highest concentration, especially under N-poor growth conditions. We predict and experimentally demonstrate that V N formation can be suppressed in thin films through growth in N-rich environments. We also find that Al 1-x Gd x N alloys are prone to high levels of unintentional O incorporation, which indirectly leads to even higher concentrations of deep defects. Growth under N-rich/reducing conditions is predicted to minimize and partially alleviate the effects of O incorporation. The results of this study provide valuable insights into the defect behavior in wurtzite nitride-based alloys, which can guide their design and optimization for various applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Rippled metamaterials with scale-dependent tailorable elasticity

Thermally induced ripples are intrinsic features of nanometer-thick films, atomically thin materials, and cell membranes, significantly affecting their elastic properties. Despite decades of theoretical studies on the mechanics of suspended thermalized sheets, controversy still exists over the impact of these ripples, with conflicting predictions about whether elasticity is scale-dependent or scale-independent. Experimental progress has been hindered so far by the inability to have a platform capable of fully isolating and characterizing the effects of ripples. This knowledge gap limits the fundamental understanding of thin materials and their practical applications. Here, we show that thermal-like static ripples shape thin films into a class of metamaterials with scale-dependent, customizable elasticity. Utilizing a scalable semiconductor manufacturing process, we engineered nanometer-thick films with precisely controlled frozen random ripples, resembling snapshots of thermally fluctuating membranes. Resonant frequency measurements of rippled cantilevers reveal that random ripples effectively renormalize and enhance the average bending rigidity and sample-to-sample variations in a scale-dependent manner, consistent with recent theoretical estimations. The predictive power of the theoretical model, combined with the scalability of the fabrication process, was further exploited to create kirigami architectures with tailored bending rigidity and mechanical metamaterials with delayed buckling instability.

Applied Physical Sciences↗
DOE OSTI2024

Convergence in simulating global soil organic carbon by structurally different models after data assimilation

Abstract Current biogeochemical models produce carbon–climate feedback projections with large uncertainties, often attributed to their structural differences when simulating soil organic carbon (SOC) dynamics worldwide. However, choices of model parameter values that quantify the strength and represent properties of different soil carbon cycle processes could also contribute to model simulation uncertainties. Here, we demonstrate the critical role of using common observational data in reducing model uncertainty in estimates of global SOC storage. Two structurally different models featuring distinctive carbon pools, decomposition kinetics, and carbon transfer pathways simulate opposite global SOC distributions with their customary parameter values yet converge to similar results after being informed by the same global SOC database using a data assimilation approach. The converged spatial SOC simulations result from similar simulations in key model components such as carbon transfer efficiency, baseline decomposition rate, and environmental effects on carbon fluxes by these two models after data assimilation. Moreover, data assimilation results suggest equally effective simulations of SOC using models following either first‐order or Michaelis–Menten kinetics at the global scale. Nevertheless, a wider range of data with high‐quality control and assurance are needed to further constrain SOC dynamics simulations and reduce unconstrained parameters. New sets of data, such as microbial genomics‐function relationships, may also suggest novel structures to account for in future model development. Overall, our results highlight the importance of observational data in informing model development and constraining model predictions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Materials Dynamics

The Materials Dynamics area of leadership focuses on understanding process-structure-properties-performance (PSPP) relationships for the extreme conditions of dynamic loading. This research encompasses controlled synthesis of materials to meet dynamic performance requirements and entails computational coupling across length and time scales for three-dimensional microstructure modeling. For this leadership area, we define dynamic loading as strain rates ≥ 10 3 /s and often, high pressures. A key grand challenge of this area is to predict and measure the evolution of microstructural phases, defect structures, and electronic structure under dynamic conditions while also measuring local temperature to understand transition states. Solving this challenge will require agile, multi-dimensional data analysis and interpretation capability.

36 MATERIALS SCIENCE↗