Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “modeling and simulation”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Structural studies of the IFNλ4 receptor complex using cryoEM enabled by protein engineering

Abstract IFNλ4 has posed a conundrum in human immunology since its discovery in 2013, with its expression linked to complications with viral clearance. While genetic and cellular studies revealed the detrimental effects of IFNλ4 expression, extensive structural and functional characterization has been limited by the inability to express and purify the protein, complicating explanations of its paradoxical behavior. In this work, we report a method for robust production of IFNλ4. We then use yeast surface display to affinity-mature IL10Rβ and solve the 72 kilodalton structures of IFNλ4 (3.26 Å) and IFNλ3 (3.00 Å) in complex with their receptors IFNλR1 and IL10Rβ using cryogenic electron microscopy. Comparison of the structures highlights differences in receptor engagement and reveals a distinct 12-degree rotation in overall receptor geometry, providing a potential mechanistic explanation for differences in cell signaling, downstream gene induction, and antiviral activities. Further, we perform a structural analysis using molecular modeling and simulation to identify a unique region of IFNλ4 that, when replaced, enables secretion of the protein from cells. These findings provide a structural and functional understanding of the IFNλ4 protein and enable future comprehensive studies towards correcting IFNλ4 dysfunction in large populations of affected patients.

Science & Technology - Other Topics↗
DOE OSTI2025

Drivers and impacts of sediment deposition in Amazonian floodplains

Abstract The Amazon River carries enormous amounts of sediment from the Andes mountains, much of which is deposited in its floodplains. However, accurate quantification of the sediment sink at fine spatiotemporal scales is still challenging. Here, we present a high-resolution hydrodynamic-sediment model to simulate sediment deposition in a representative Amazon/Solimões floodplain. The process is found to be jointly driven by inundation, suspended sediment concentration in the Amazon River, and floodplain hydrodynamics and only weakly correlated with inundation level. By upscaling the sediment deposition rate (1.33 ± 0.24 kg m −2 yr −1 ), we estimate the trapping of 77.3 ± 13.9 Mt (or 6.1 ± 1%) of the Amazon River sediment by the Amazon/Solimões floodplains every year. Widespread deforestation would reduce the trapping efficiency of the floodplains over time, exacerbating downstream river aggradation. Additionally, we show that the deposition of sediment-associated organic carbon plays a minor role in fueling carbon dioxide and methane emissions in the Amazon.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Spontaneously formed phonon frequency combs in van der Waals solid CrGeTe 3 and CrSiTe 3

Optical phonon engineering through nonlinear effects has been utilized in ultrafast control of material properties. However, nonlinear optical phonons typically exhibit rapid decay due to strong mode-mode couplings, limiting their effectiveness in temperature or frequency sensitive applications. Here we report the observation of long-lived nonlinear optical phonons through the spontaneous formation of phonon frequency combs in the van der Waals material CrXTe 3 (X=Ge, Si) using high-resolution Raman scattering. Unlike conventional optical phonons, the highest A g mode in CrGeTe 3 splits into equidistant, sharp peaks forming a frequency comb that persists for hundreds of oscillations and survives up to 200K. These modes correspond to localized oscillations of Ge 2 Te 6 clusters, isolated from Cr hexagons, behaving as independent quantum oscillators. Introducing a cubic nonlinear term to the harmonic oscillator model, we simulate the phonon time evolution and successfully replicate the observed comb structure. Similar frequency comb behavior is observed in CrSiTe 3 , demonstrating the generalizability of this phenomenon. Our findings demonstrate that Raman scattering effectively probes high-frequency nonlinear phonon modes, offering insight into the generation of long-lived, tunable phonon frequency combs with potential applications in ultrafast material control and phonon-based technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Valence-free open nanoparticle superlattices

A cornerstone of advanced materials design is establishing a framework for assembling nanoparticle superstructures with tailored symmetries. A longstanding challenge has been assembling diamond-like superstructures for photonic devices. Traditionally, such open superstructures require functionalized nanoparticles with directional or anisotropic interactions, reminiscent of valence bonding in a diamond. Here, we present a robust strategy for assembling valence-free nanoparticles into a broad array of cubic superstructures. By grafting nanoparticles with oppositely charged, end-functionalized water-soluble polymers of adjustable molecular weight, we gain control over electrostatic interactions and conformational constraints. This unified approach yields lattices analogous to rock salt, CsCl, zinc-blende, diamond, and the rare simple cubic phase, with tunable lattice constants. Theoretical models and simulations elucidate the underlying interactions, providing a framework for engineering valence-free nanoparticle superlattices.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Vegetation responses to air dryness amplify future land surface warming

Temperature exerts a first-order control on vegetation photosynthesis and transpiration. Yet most studies investigating temperature impacts on plants rely on near-surface air temperature, rather than canopy temperature—the temperature plants actually experience. Because canopy temperature directly regulates ecosystem function, it provides a more accurate measure of vegetation–climate interactions. Combining Earth System Model (ESM) simulations and satellite observations in a dual emergent constraint, here we show that canopy temperature is projected to increase substantially more than air temperature (~0.11-degrees more or a 16% increase in their difference) over the 21st century. The ESM ensemble median fails to capture these stronger increases in the majority of vegetated regions. We find that the largest projected increases in the difference between canopy and air temperature are predicted to occur in regions where elevating moisture stress—particularly rising vapor pressure deficit—increasingly constrains vegetation growth and transpiration. This implies that future warming will impose stronger constraints on plant function than currently estimated. Relying on air temperature alone will therefore lead to systematic underestimation of temperature effects on photosynthesis, vegetation growth, and the land carbon sink. Accurate representation of canopy temperature in ESMs is thus essential to improve projections of ecosystem responses and feedbacks to climate change.

Green, Julia K↗
DOE OSTI2025

Unveiling nano-scale crystal deformation using coherent X-ray dynamical diffraction

Visualization of internal deformation fields in crystalline materials helps bridge the gap between theoretical models and practical applications. Applying Bragg coherent diffraction imaging under X-ray dynamical diffraction conditions provides a promising approach to the longstanding challenge of investigating the deformation fields in micron-sized crystals. Here, we present an automatic differentiation-based reconstruction method that integrates dynamical scattering theory to accurately reconstruct deformation fields in large crystals. Using this forward model, our simulated and experimental results demonstrate that three-dimensional local strain information inside a large crystal can be accurately reconstructed under coherent X-ray dynamical diffraction conditions with Bragg coherent X-ray diffraction imaging. These findings open an avenue for extending the investigation of local deformation fields to microscale crystals while maintaining nanoscale resolution, leveraging the enhanced coherence and brightness of advanced X-ray sources.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Unconventional nonlinear Hall effects in twisted multilayer 2D materials

We present the first investigation of unusual nonlinear Hall effects in twisted multilayer 2D materials. Contrary to expectations, our study shows that these nonlinear effects are not merely extensions of their monolayer counterparts. Instead, we find that stacking order and pairwise interactions between neighboring layers, mediated by Berry curvatures, play a pivotal role in shaping their collective nonlinear optical response. By combining large-scale Real-Time Time-Dependent Density Functional Theory (RT-TDDFT) simulations with model Hamiltonian analyses, we demonstrate a remarkable second-harmonic transverse response in hexagonal boron nitride four-layers, even in cases where the total Berry curvature cancels out. Furthermore, our symmetry analysis of the layered structures provides a simplified framework for predicting nonlinear responses in multilayer materials in general. Our investigation challenges the prevailing understanding of nonlinear optical responses in layered materials and opens new avenues for the design and development of advanced materials with tailored optical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Detecting tropospheric composition and climate responses to US air pollution controls in the context of internally-arising variability

Since the 1970s, air pollutant emissions controls in the United States (US) have lowered concentrations of ozone (O 3 ) and aerosols, which have opposing radiative effects on surface temperature. Using a pair of initial-condition ensembles generated by a fully-coupled chemistry-climate model, we simulate the “world avoided” by US air pollution controls. In this counterfactual world, we find tropospheric column O 3 increases, robust to natural internal variability, extending across the Northern Hemisphere. Robust aerosol increases, dominated by sulfate, remain localized near the US. Wintertime Northwest Atlantic cloud droplet number concentration is particularly sensitive to US aerosol. While an ensemble mean US surface cooling signal (−0.4 °C) implies that aerosol-driven cooling prevails over any O 3 -induced warming, we find that large regional internal variability will confound its detection in any single transient realization. Larger signal-to-noise ratios for composition versus climate variables underscore the greater detectability of emissions-driven changes in tropospheric composition compared to their associated climate impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Deciphering yield modification of hadron-triggered semi-inclusive recoil jets in heavy-ion collisions

Jet quenching is recognized as critical evidence for the existence of the quark-gluon plasma (QGP) and serves as an essential probe to study its transport properties. Measurements of hadron-triggered semi-inclusive recoil jets have gained popularity due to its capability to probe jets over an extended phase space at low transverse momenta (p T ) and large radii. Recent ALICE measurements showed that the I AA , yield ratio of recoil jets between heavy-ion and p+p collisions, rises with jet p T and exceeds unity at high p T , contradicting conventional expectations that jet quenching should result in I AA values less than one. In this contribution, we re-examine the surface bias and study the effects of energy losses for both trigger hadrons and recoil jets on I AA , employing the Linear Boltzmann Transport (LBT) model to simulate jet-medium interactions. Our findings suggest that a large portion of hadrons used for the triggers undergoes substantial energy loss, despite surface bias. In particular, the energy loss of the trigger hadrons elevates the I AA baseline, corresponding to the case of no energy loss for recoil jets, to be greatly larger than unity. This enhancement of the baseline implies that the measured I AA values being larger than unity could still signal jet quenching.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Modeling and observations of North Atlantic cyclones: Implications for U.S. Offshore wind energy

To meet the Biden-Harris administration's goal of deploying 30 GW of offshore wind power by 2030 and 110 GW by 2050, expansion of wind energy into U.S. territorial waters prone to tropical cyclones (TCs) and extratropical cyclones (ETCs) is essential. This requires a deeper understanding of cyclone-related risks and the development of robust, resilient offshore wind energy systems. Here, this paper provides a comprehensive review of state-of-the-science measurement and modeling capabilities for studying TCs and ETCs, and their impacts across various spatial and temporal scales. We explore measurement capabilities for environments influenced by TCs and ETCs, including near-surface and vertical profiles of critical variables that characterize these cyclones. The capabilities and limitations of Earth system and mesoscale models are assessed for their effectiveness in capturing atmosphere–ocean–wave interactions that influence TC/ETC-induced risks under a changing climate. Additionally, we discuss microscale modeling capabilities designed to bridge scale gaps from the weather scale (a few kilometers) to the turbine scale (dozens to a few meters). We also review machine learning (ML)-based, data-driven models for simulating TC/ETC events at both weather and wind turbine scales. Special attention is given to extreme metocean conditions like extreme wind gusts, rapid wind direction changes, and high waves, which pose threats to offshore wind energy infrastructure. Finally, the paper outlines the research challenges and future directions needed to enhance the resilience and design of next-generation offshore wind turbines against extreme weather conditions.

17 WIND ENERGY↗
DOE OSTI2025

High-compression implosions based on high density carbon ablator using modified drive and capsule dopant profiles

Laser-driven inertial fusion experiments have, for the first time, achieved a target gain greater than unity in a laboratory setting [Abu-Shawareb et al., Phys. Rev. Lett. 132, 065102 (2024)]. Despite this breakthrough, the burn-up fraction remains limited to about one-fourth of ideal estimates due to insufficient areal density, highlighting the potential for greater gains through enhanced compression. In our previous work, we demonstrated record-high compression of stagnated fuel in indirectly driven implosions using high-density carbon ablators. This was achieved by combining a continuous ramped pulse drive with a modified ablator dopant profile, which reduced mixing at the fuel–ablator interface and improved stability [Tommasini et al., Phys. Rev. Res. 5, L042034 (2023)]. Based on this foundation, the study presented here investigates the limits of compression achievable by combining the continuous ramped pulse drive with different dopant profiles to further minimize unstable interfaces and gradient discontinuities, thereby reducing fuel–ablator mixing. Our results demonstrate that the continuous ramped pulse consistently outperforms designs based on 3-shock drive pulses across all ablator profiles studied, with compression showing only a relatively modest dependence on dopant configurations that reduce the number of interfaces or eliminate discontinuities in the dopant gradient profile. Sub-scale experiments using the continuous ramped pulse achieved compression levels exceeding those of full-scale “HyE” implosions [Kritcher et al., Phys. Plasmas 28, 072706 (2021)] at similar adiabat, anticipating significant performance gains with increased scale, as supported by models and simulations. These findings underscore the critical role of the continuous ramped pulse in reducing mix and achieving improved compression. They also provide a foundation for future large-scale experiments to test the continuous ramped pulse design on deuterium–tritium fuel in the burn-wave propagation regime, leveraging the most effective combinations of continuous ramped pulse and dopant profiles identified in this study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Crack formation in strained β-(AlxGa 1−x ) 2 O 3 films grown on (010) β-Ga 2 O 3 substrates

The cracking and local strain relaxation in (010) (Al x Ga 1−x ) 2 O 3 films grown on Ga 2 O 3 substrates are assessed in terms of film composition and thickness. We utilize x-ray diffraction and electron microscopy techniques combined with simulation and modeling to investigate that film cracking on curvature (flatness) has a directly proportional relationship with film thickness and/or aluminum content. Cross section transmission electron microscopy reveals that cracks along both the (001) and (100) cleavage planes penetrate into the substrate. The diffuse scattered intensity observed in reciprocal space maps (RSMs) is directly correlated with the tilt that is introduced due to the change in the deformation conditions near the cracks. While asymmetric RSMs show that the layers are fully strained, the diffuse scattering distribution in reciprocal space can be interpreted to show that the cracking relaxes and locally tilts the lattice ∼350 nm from the crack edges, which is consistent with the larger radius of curvature associated with the films with higher crack densities. For example, a 200 nm (Al 0.13 Ga 0.87 ) 2 O 3 thick film has an average inter-crack spacing of 3.3 µm, so most of the epitaxial layer is fully strained except near the cracks where it deforms elastically and is consistent with the gallium oxide Poisson ratio. A reciprocal space model was developed, which imports the strain and tilt distributions (based on finite element modeling) to match the features observed in the experimental maps. We also note that previous studies involving (Al x Ga 1−x ) 2 O 3 films may show evidence of cracking as observed in their symmetric and asymmetric RSMs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Observed declines in upper ocean phosphate-to-nitrate availability

Climate warming is increasing ocean stratification, which in turn should decrease the nutrient flux to the upper ocean. This may slow marine primary productivity, causing cascading effects throughout food webs. However, observing changes in upper ocean nutrients is challenging because surface concentrations are often below detection limits. We show that the nutricline depth, where nutrient concentrations reach well-detected levels, is tied to productivity and upper ocean nutrient availability. Next, we quantify nutricline depths from a global database of observed vertical nitrate and phosphate profiles to assess contemporary trends in global nutrient availability (1972–2022). We find strong evidence that the P-nutricline (phosphacline) is mostly deepening, especially throughout the southern hemisphere, but the N-nutricline (nitracline) remains mostly stable. Earth System Model (ESM) simulations support the hypothesis that reduced iron stress and increased nitrogen fixation buffer the nitracline, but not phosphacline, against increasing stratification. These contemporary trends are expected to continue in the coming decades, leading to increasing phosphorus but not nitrogen stress for marine phytoplankton, with important ramifications for ocean biogeochemistry and food web dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Simulating regional workforce impacts of decarbonizing integrated steelmaking

Global efforts to mitigate climate change are increasing pressure on heavy manufacturing industries to decarbonize production. The iron and steel industry is responsible for 7% of CO 2 emissions globally (2% in the United States) and is often a major employer in the regions where iron and steel is produced. Understanding the future prospects for workers in regions with high CO 2 emitting industries—including impacts of phasing out or evolving such industries—will be critical for informing regional economic and clean energy strategies. We simulate the impact of an “in-place” transition that replaces today’s integrated production with direct reduced iron (DRI) used in electric arc furnaces (EAFs), using Southwest Pennsylvania as an application of our generalizable approach. Our results suggest that the integrated steelmaking workforce today has the skills, knowledge, and abilities (SKAs) to fill over 95% of all jobs required by DRI/EAF facilities, but the number of jobs is only 25% of those at integrated plants. We also find that some occupational groups have greater general transferability into the broader job market, while other groups, such as production workers, are ill-equipped today based on current SKAs to transition out of the iron and steel industry. Our methodology further suggests factors that limit transitions: Around 85% of occupations are more limited by missing skills, while 15% are more limited by insufficient wages. These results may help to improve the design of social policy and the targeting of retraining programs, while the simulation approach can be readily adapted for other regions and industries.

decarbonization↗
DOE OSTI2026

Optimisation of the Kaplan hydropower system via PID 2 and digital twin

Here, this paper proposes a proportional–integral-double–derivative (PID 2 ) optimisation method for the Kaplan hydropower system by building a digital twin. The study first uses one multilayer perceptron (MLP) to model the hydroturbine dynamic and then adopts three connected MLPs to model the generator dynamic, both in an open-loop fashion. Inspired by stochastic distribution control (SDC) theory, we regard the training of the turbine's neural network model as a process control problem, and we propose minimising entropy loss to update the network parameters. The next step is to build the digital twin by connecting the neural network models with a PID 2 controller and a lead-lag exciter and run the whole model in a closed-loop fashion. After that, a binary search approach is applied to optimise the PID 2 parameters based on the obtained digital twin model. The simulation results show that the proposed method can reduce the mean square tracking error by more than 90%. Furthermore, the method is extended to jointly optimise the PID 2 controller and excitation system gains through multiobjective optimisation, leveraging Pareto frontier analysis to balance active power and voltage tracking performance. Simulation results confirm the effectiveness of the proposed method, achieving a 83.46% reduction in relative mean square error of active power, a 47.13% reduction in terminal voltage tracking error, and an 82.78% improvement in the overall scalarized objective.

Hydropower system↗
DOE OSTI2025

Path Toward a Breeding, Proliferation-Resistant, Thermal-Spectrum Molten Salt Reactor

Liquid-fueled, thermal-spectrum molten salt breeder reactors (TS-MSBRs) offer the potential for affordable, safe, inexhaustible energy with minimal potential for nuclear material misuse and without significant actinide waste generation. Realizing the full set of TS-MSBR capabilities is only now becoming possible with the advent of advanced fuel-salt processing techniques, improved materials, and a more detailed understanding of fuel-salt properties. Additionally, modern higher-fidelity modeling and simulation methods enable a more detailed evaluation of TS-MSBR design options. TS-MSBRs, however, remain immature and will require substantial, sustained development resources.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Feasibility Study on Implementing a Staggered-Grid Finite Volume Method for System Analysis Code Development Under the MOOSE Framework

Here, this work summarizes a feasibility study on testing numerical algorithms that are suitable and efficient for advanced system analysis code development under the mutli-physics framework, MOOSE. The key to the test bed is the implementation of high-order one-dimensional staggered-grid finite volume method (SG-FVM), and its direct interaction with the linear/nonlinear solver, PETSc. The test bed utilized a more flexible code structure to enable the finite volume method implementation and direct interacting with the solver package, instead of using the natively supported finite element method by the framework. Using a suite of selected test problems with different problem sizes and levels of complexity, the implemented SG-FVM demonstrated superior performance improvement against a direct finite element method implementation through MOOSE. On two computer systems, the speedup was observed to be significant, with at least one order of magnitude of solving time reduction. For a complex reactor model, transient simulation was performed using the newly developed finite volume method code, the results of which agree very well with the reference results from the finite element method code. Overall, this study demonstrates a successful feasibility study on the proposed numerical algorithms and software structure to support advanced system analysis tool development.

MOOSE↗