Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “control variates”

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

Raspberry Pi–powered temperature monitoring of growth chamber microclimates

While controlled environments are desirable for growing and measuring plants, growth chambers and greenhouses typically have microclimates that impact plant growth, development, and stress responses. Furthermore, opening and closing the doors of a controlled environment introduces variation in the environment, especially at temperature extremes, affecting both the measurements and the organisms within. Using multiple temperature data loggers to normalize results can be cost-prohibitive and rarely offers real-time feedback on temperature status. We used low-cost single-board computers, cameras, and temperature sensors to manage and capture growth chamber temperatures while acquiring plant image data. Detailed here are methods to document microclimates within a growth chamber so that data can be normalized to measured temperature information. This protocol describes a low-cost method for automated temperature monitoring, which enables both high-throughput measurements of temperature along with plant growth and stress responses via plant imaging.

Plant Sciences↗
DOE OSTI2024

Lattice QCD calculation of the pion distribution amplitude with domain wall fermions at physical pion mass

We present a direct lattice QCD calculation of the x-dependence of the pion distribution amplitude (DA), which is performed using the quasi-DA in large momentum effective theory on a domain-wall fermion ensemble at physical quark masses and spacing a ≈ 0.084 fm. The bare quais-DA matrix elements are renormalized in the hybrid scheme and matched to $\overline{MS}$ with a subtraction of the leading renormalon in the Wilson-line mass. For the first time, we include threshold resummation in the perturbative matching onto the light-cone DA, which resums the large logarithms in the soft gluon limit at next-to-next-to-leading log. The resummed results show controlled scale-variation uncertainty within the range of momentum fraction x ϵ [0.25, 0.75] at the largest pion momentum P z ≈ 1.85 GeV. In addition, we apply the same analysis to quasi-DAs from a highly-improved-staggered-quark ensemble at physical pion mass and a = 0.076 fm. By comparison we find with 2σ confidence level that the DA obtained from chiral fermions is flatter and lower near x = 0.5.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Understanding spectral dependance of laser-induced damage precursors in dielectric materials (Full report_23-ERD-006)

The performance of high-energy laser systems is constrained by laser-induced damage in dielectric coatings, particularly those containing hafnium oxide (HfO 2 ). While thresholds at fundamental Nd harmonics are well studied, the spectral dependence of damage initiation—especially under dual-wavelength irradiation—remains poorly characterized. This project provides the first systematic investigation of wavelength-dependent laser damage in hafnia coatings, focusing on nanoscale precursors such as craze lines, nanobubbles, stoichiometric variations, nodules, and controlled crystallization. Coatings were fabricated via ion beam sputtering and electron-beam deposition and characterized using spectrophotometry, ellipsometry, AFM, GI-XRD, RBS, PCI absorption, and fs/ns laser damage testing. Results show that craze lines, benign under infrared light, strongly initiate damage under UV due to wavelength-selective field intensification. Substituting xenon for argon suppresses nanobubbles and improves UV thresholds by up to 32%. Oxygen modulation reveals that fully oxidized films maximize UV resistance, though at the cost of porosity and stress in multilayers. HfO 2 –SiO 2 composites resist crystallization and defects while achieving ppm-level absorption and elevated thresholds, whereas full crystallization of HfO 2 enhances LIDT by reducing defects and improving thermal transport. Collectively, these findings link photon energy, defect states, and bandgap collapse, providing a predictive framework for wavelength-dependent laser damage. The outcomes directly inform the design of durable, multi-wavelength coatings for facilities such as NIF, MEC, HAPLS, and DPAL, advancing the readiness of next-generation optics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Resolving Low Cloud Feedbacks Globally With E3SM High‐Res MMF: Agreement With LES but Stronger Shortwave Effects

This study investigates low cloud feedback in a warmer climate using global simulations from the High-Resolution Multi-scale Modeling Framework (HR-MMF), which explicitly simulates small-scale eddies globally. Two 5-year simulations—one with present-day sea surface temperatures (SSTs) and a second with SSTs warmed uniformly by 4 K—reveal a positive global shortwave cloud radiative effect (SWCRE = 0.3 W/m 2 /K), comparable to estimates from CMIP models. As the climate warms, significant reductions in low cloud cover occur over stratocumulus regions. This study is the first attempt to compare HR-MMF results with predictions from idealized large-eddy simulations from the CGILS intercomparison. Despite different underlying assumptions, we find qualitative agreement in SWCRE and inversion height changes between HR-MMF and CGILS predictions. This suggests reasonable credibility for the CGILS framework in predicting cloud responses under the out-of-sample conditions found in HR-MMF. However, the HR-MMF exhibits stronger SWCRE changes than predicted by CGILS. We explore potential causes for this discrepancy, examining variations in cloud-controlling factors (CCFs) and cloud conditions. Our results show a fairly homogeneous SWCRE response, with little systematic variation tied to the variations in CCFs. This reveals a dominant role for SST forcing in modulating SWCRE.

boundary-layer clouds↗
DOE OSTI2026

Strain-stabilized interfacial polarization tunes work function over 1 eV in RuO 2 /TiO 2 heterostructures

Interfacial polarization – charge accumulation at the heterointerface – is a well-established tool in semiconductors, but its influence in metals remains unexplored. Here, we demonstrate that interfacial polarization can robustly modulate surface work function in metallic rutile RuO 2 layers in epitaxial RuO 2 /TiO 2 heterostructures grown by hybrid molecular beam epitaxy. Using multislice electron ptychography, we directly visualize polar displacements of transition metal ions relative to oxygen octahedra near the interface, despite the conductive nature of RuO 2 . This interfacial polarization enables over 1 eV modulation of the RuO 2 work function, controlled by small thickness variations (2-3 nm), as measured by Kelvin probe force microscopy, with a critical thickness of ~4 nm – corresponding to the transition from fully-strained to relaxed film. These results establish interfacial polarization as a powerful route to control electronic properties in metals and have implications for designing tunable electronic, catalytic, and quantum devices through interfacial control in polar metallic systems.

Electronic properties and materials↗
DOE OSTI2026

Satellite Reentry Predictions During Sudden Stratospheric Warmings

Prediction of uncontrolled satellite reentry time and trajectory are essential to avoid damages to human being and properties over land. Reentry time and trajectory are largely controlled by aerodynamic drag and mass neutral density below 200 km. In this paper, using a numerical model of the whole atmosphere, WACCM‐X, and a Precision Orbital Determination software, we demonstrate the major sudden stratospheric warming (SSW) in January 2009, during the deep solar minimum of 2009, may cause 0.2–0.3 day of variation in reentry time. The SSW can alter the mass density and composition mixing ratio in the lower thermosphere below 200 km through changed wave forcing and general circulation in the lower thermosphere. This is comparable to the effect of a 2‐SFU (Solar Flux Unit) solar radiance variation on the neutral density. 0.1 day of reentry time uncertainty corresponds to a global shift in impact location, that is, thousands of kilometers distance. For comparison, a major SSW during solar maximum contributes little to reentry time variation, which is predominantly controlled by the solar and geomagnetic forcing. We suggest a few proxies such as O/N 2 in the lower thermosphere or mean meridional wind to account for the lower atmosphere impacts in an operational model. This needs to be investigated in rigorous future work.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Wind turbine power phase control with DC collection bus for onshore/offshore windfarms

A DC bus collection system for a wind farm reduces the overall required number of converters and minimizes the energy storage system requirements. The DC bus collection system implements a power phasing control method between wind turbines that filters the variations and improves power quality. The phasing control method takes advantage of a novel power packet network concept with nonlinear power flow control design techniques that guarantees both stable and enhanced dynamic performance.

Weaver, Wayne W.↗
DOE OSTI2025

Integrating vehicle trajectory planning and arterial traffic management to facilitate eco-approach and departure deployment

Eco-approach and departure (EAD) enable continuous vehicle motion in urban signalized corridors. Since such a motion can extend to the EAD vehicles’ followers, it makes EAD a promising technology to benefit the traffic flow where automated vehicles and conventional vehicles coexist. Most existing EAD studies envision an ideal setting that neglects real-world operational conditions such as lane changes, multi-movement intersection configuration, partially automated fleet, and/or limited traffic state awareness. This study aims to fill the gap by designing an EAD algorithm considering real-world traffic operation constraints. The proposed algorithm uses a model predictive controller to minimize vehicle speed reduction and variation based on the real-time traffic signal control plan and measured queues at the intersection. The required inputs are readily available at many modern intersections. We observed that the proposed controller’s performance might degrade because of lane-changing maneuvers and lead-left turn traffic signals. These observations motivated our development of a lane change management strategy and a signal control implementation strategy to facilitate the EAD implementation. The lane change management strategies separate the EAD operations and lane-changing maneuvers in time and space. The signal control implementation strategy applies lag-left turn signals to enable EAD operation for both the through and left-turn vehicles. Compared to the non-EAD case, our EAD approach produces 2.5% to 7.8% energy savings while keeping similar intersection mobility. Notably, this approach brings about 2.5% to 3.6% energy savings in a 2% CAV case. This result demonstrates the feasibility of deploying EAD at low connected automated vehicle penetration rates.

Arterial corridor management↗
DOE OSTI2026

Tailored Solvent Treatment for Optimized Production of Upcycled Anodes from End-Of-Life Li-Ion Batteries

Recycling processes for lithium-ion batteries typically overlook graphite because of its lower market value relative to that of transition-metal-containing cathode materials. However, graphite recovered from cycled lithium-ion batteries holds additional engineered value associated with the solid-electrolyte interphase (SEI). The SEI contributes critical electronic passivation of the graphite surface but becomes highly resistive with extended cycling, yielding poor cell performance. In this work, we apply tailored solvent treatment to end-of-life (EOL) graphite anodes to selectively remove adverse SEI components while retaining beneficially passivating species. We evaluate a series of polar protic solvents to achieve targeted removal of SEI components and control selectivity through rational variation in solvent properties. The physiochemical properties of treatment solvents correlate with both the retained SEI composition and the corresponding electrochemical performance of solvent-treated “upcycled” graphite anodes. Within the initial set of solvents evaluated, top-performing candidates show capacity and Coulombic efficiency nearly equivalent to those of an analogous pristine anode, as well as promising electrochemical performance enhancement with regard to irreversible capacity-loss metrics. This study establishes critical design principles for an optimized anode upcycling method that enhances the value of recycled graphite by retaining and upgrading the SEI.

25 ENERGY STORAGE↗
DOE OSTI2026

Analysis of a Heat Balance Experiment in the TREAT Reactor with SPARTA and OpenMC

This article presents the interpretation of a heat balance experiment in TREAT using the SPARTA methodology and an OpenMC model. The study examines the effects of control rods and temperature variations on the effective multiplication factor ($k_{\rm eff}$) and correction factors. Results indicate that both core temperature and rod positions have a linear effect on $k_{\rm eff}$ and correction factors in the upper-west channel. A critical hypothesis of SPARTA was validated: reactivities and correction factors can be obtained by linearly combining individual effects. The previously considered negligible effect of temperature on correction factors was found to be significant, necessitating updates to the SPARTA methodology. By fitting an exponential curve to mean temperature data, SPARTA's accuracy and convergence have been improved. Comparisons between SPARTA and OpenMC showed consistent reactivity removal due to temperature increase. While SPARTA effectively interprets complex transients, simpler methods may suffice when correction factors have minimal influence on reactivity calculations.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Elastic Bayesian Model Calibration

Functional data are ubiquitous in scientific modeling. For instance, quantities of interest are modeled as functions of time, space, energy, density, etc. Uncertainty quantification methods for computer models with functional response have resulted in tools for emulation, sensitivity analysis, and calibration that are widely used. However, many of these tools do not perform well when the computer model’s parameters control both the amplitude variation of the functional output and its alignment (or phase variation). This paper introduces a framework for Bayesian model calibration when the model responses are misaligned functional data. The approach generates two types of data out of the misaligned functional responses: (1) aligned functions so that the amplitude variation is isolated and (2) warping functions that isolate the phase variation. These two types of data are created for the computer simulation data (both of which may be emulated) and the experimental data. The calibration approach uses both types so that it seeks to match both the amplitude and phase of the experimental data. The framework is careful to respect constraints that arise, especially when modeling phase variation, and is framed in a way that it can be done with readily available calibration software. In conclusion, we demonstrate the techniques on two simulated data examples and on two dynamic material science problems: a strength model calibration using flyer plate experiments and an equation of state model calibration using experiments performed on the Sandia National Laboratories’ Z-machine.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Kinetic Pathways to Diblock Copolymer Micelles via Thin Film Dissolution and Cosolvent Variation

We investigate how processing pathways control micelle formation during dissolution of lamellar diblock copolymers using dissipative particle dynamics (DPD) simulations and DPD-selfconsistent field theory (DPD-SCFT). Comparing thin film (TF) and cosolvent-assisted (CS) dissolution reveals fundamentally different micellization mechanisms and final states. In the TF pathway, dissolution proceeds through a metastable cylindrical intermediate that dictates final micelle dimensions. The lamellar melt first transforms into cylindrical micelles, which then fragment into spherical micelles. The final spherical radius is intrinsically linked to the cylinder radius by a constant ratio, Rsphere ≈ 1.4Rcylinder. DPD-SCFT confirms that both cylinders and spheres are near their equilibrium dimensions. In contrast, CS dissolution produces kinetically trapped micelles with sizes significantly below equilibrium. Micellization is hindered by high free energy barriers to single-chain exchange and fusion. Growth occurs primarily via fusion, but increasing corona repulsion suppresses further coarsening, preventing equilibration. These results demonstrate that, when starting from a bulk copolymer phase, the metastable cylindrical intermediate plays a central role in determining final micelle size.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Induced and natural variation affect traits independently in hybrid Populus

Abstract The genetic control of many plant traits can be highly complex. Both allelic variation (sequence change) and dosage variation (copy number change) contribute to a plant's phenotype. While numerous studies have investigated the effect of allelic or dosage variation, very few have documented both within the same system, leaving their relative contribution to phenotypic effects unclear. The Populus genome is highly polymorphic, and poplars are fairly tolerant of gene dosage variation. Here, using a previously established Populus hybrid F1 population, we assessed and compared the effect of natural allelic variation and induced dosage variation on biomass, phenology, and leaf morphology traits. We identified QTLs for many of these traits, but our results indicate limited overlap between the QTLs associated with natural allelic variation and induced dosage variation. Additionally, the integration of data from both allelic and dosage variation identifies a larger set of QTLs that together explain a larger percentage of the phenotypic variance. Finally, our results suggest that the effect of the large indels might mask that of allelic QTLs. Our study helps clarify the relationship between allelic and dosage variation and their effects on quantitative traits.

Guo, Weier (ORCID:0000000251789334)↗
DOE OSTI2024

In‐mold rheology and automated process control for injection molding of recycled polypropylene

Abstract Manufacturing plastic parts with secondary feedstocks has risen to the forefront of importance in recent years. However, the variation in molecular weight and rheology of secondary feedstock can lead to inconsistent part quality. This work evaluates the effectiveness of a novel closed‐loop adaptive process control system that adjusts nozzle pressure in response to in‐mold pressure data. Five different recycled polypropylene blends, with a broad distribution of flow properties, were evaluated to determine the effectiveness of the control system at reducing processing variation. The experimental results show that the process control strategy reduced the variation within the mold, as seen by in‐mold pressure curves and calculated in‐mold viscosity values. Additionally, the parameters that control the automated process adjustments were investigated, showing the importance of optimization. The analysis of the correlation between in‐mold rheology and mechanical properties showed a slight variation in the mechanical properties and parts weight with a coefficient of variation of under 5%. Overall, the results demonstrate the ability of pressure‐controlled molding and automated viscosity adjustment to reduce the variability when molding a secondary feedstock. Highlights Pressure‐controlled injection molding of recycled polypropylene. Automated closed‐loop adaptive process control methodology. Methodology resulted in a reduction in pressure variation during molding. Changes in mechanical properties and in‐mold viscosity were investigated. Results show the potential of pressure‐controlled molding at reducing variation.

Krantz, Joshua↗
DOE OSTI2024

The AGORA High-resolution Galaxy Simulations Comparison Project. IV. Halo and Galaxy Mass Assembly in a Cosmological Zoom-in Simulation at z ≤ 2

In this fourth paper from the AGORA Collaboration, we study the evolution down to redshift z = 2 and below of a set of cosmological zoom-in simulations of a Milky Way mass galaxy by eight of the leading hydrodynamic simulation codes. We also compare this CosmoRun suite of simulations with dark matter-only simulations by the same eight codes. We analyze general properties of the halo and galaxy at z = 4 and 3, and before the last major merger, focusing on the formation of well-defined rotationally supported disks, the mass–metallicity relation, the specific star formation rate, the gas metallicity gradients, and the nonaxisymmetric structures in the stellar disks. Codes generally converge well to the stellar-to-halo mass ratios predicted by semianalytic models at z ~ 2. We see that almost all the hydro codes develop rotationally supported structures at low redshifts. Most agree within 0.5 dex with the observed mass–metallicity relation at high and intermediate redshifts, and reproduce the gas metallicity gradients obtained from analytical models and low-redshift observations. We confirm that the intercode differences in the halo assembly history reported in the first paper of the collaboration also exist in CosmoRun, making the code-to-code comparison more difficult. We show that such differences are mainly due to variations in code-dependent parameters that control the time stepping strategy of the gravity solver. We find that variations in the early stellar feedback can also result in differences in the timing of the low-redshift mergers. All the simulation data down to z = 2 and the auxiliary data will be made publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

How Does Escherichia coli Allocate Proteome?

Microorganisms are shown to actively partition their intracellular resources, such as pro- teins, for growth optimization. Recent experiments have begun to reveal molecular com- ponents unpinning the partition; however, it remains unclear quantitatively how individual parts orchestrate to yield precise resource allocation that is both robust and dynamic. Here we developed a coarse-grained mathematical framework that centers on guanosine pentaphosphate (ppGpp)-mediated regulation, and used it to systematically uncover the design principles of proteome allocation in Escherichia coli. Our results showed that cellular ability of resource partition lies in an ultrasensitive, negative feedback control- ling topology with the ultrasensitivity arising from zero-order amino acid kinetics and the negative feedback from ppGpp-controlled ribosome synthesis. In addition, together with the time-scale separation between slow ribosome kinetics and fast turnovers of ppGpp and amino acids, the network topology confers the organism an optimization mechanism which mimics sliding mode control, a nonlinear optimization strategy that is widely used in man-made systems. We further showed that such a controlling mechanism is robust against parameter variations and molecular fluctuations, and is also efficient for biomass production over time. Furthermore, this work elucidates the fundamental controlling mechanism of E. coli proteome allocation, thereby providing insights into quantitative microbial physiology as well as the design of synthetic gene networks.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

A Comprehensive Strategy for Grid Forming Control in DC Coupled Photovoltaic and Battery Energy Storage Inverters

This paper presents an integrated DC-DC and DCAC grid-forming control strategy for DC-coupled photovoltaic (PV) plus battery energy storage systems, considering the effect of DC link voltage variations caused by direct PV connections. A power reference algorithm determines power distribution between the PV and battery to the grid while observing device power ratings to prevent the over-rating of components and keep the battery's state of charge within an acceptable range. The simulated utility-scale model in MATLAB/Simulink illustrates its ability against extreme phase angle variation contingencies in the grid while controlled through grid-forming control with a fast dynamic on DC link voltage. The simulation results confirm the effectiveness of the proposed control in integrating PV plus battery configurations with grid forming control and maintaining reliable grid operation under severe grid disturbances.

battery, boost, control, energy storage, grid form↗
DOE OSTI2026

Definitive Assessment of the Accuracy, Variationality, and Convergence of Relativistic Coupled Cluster and Density Matrix Renormalization Group in 100-Orbital Space

Accuracy, variationality, and convergence underpin the reliability of modern electronic structure methods, yet definitive benchmarks in the relativistic regime remain elusive due to the absence of numerically exact full configuration interaction (CI) references. Recent algorithmic advances in the CI framework, enabled by the small-tensor-product (STP) decomposition approach, have dramatically extended the tractable size of the configuration space, making numerically exact CI calculations feasible in large active spaces previously beyond reach. In this paper, we employ the recently developed STP-CI framework to perform large-scale numerically exact CI calculations and directly benchmark relativistic coupled cluster and density matrix renormalization group methods. Definitive benchmarking of approximate relativistic electronic structure methods is ensured through the application of the gap theorem, which provides rigorous error bounds on the CI reference and establishes a controlled standard for assessing accuracy, variationality, and convergence.

Chemical calculations↗