Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TRAJECTORIES”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Fuel Property Effects on Knock Propensity and Thermal Efficiency in a Direct-Injection Spark-Ignition Engine

Engine knock remains one of the major barriers to further improvement in thermal efficiency of Direct-Injection Spark-Ignition (DISI) engines. While Research Octane Number and Motor Octane Number are often used as standard rating methods for knock resistance of fuels, the impacts of other fuel properties on knock propensity in modern engines such as heat of vaporization (HoV) and laminar flame speed (LFS) require better understanding in order to co-optimize fuels and engine designs to achieve higher thermal efficiency and lower CO2 emission. In the present study, computational fluid dynamics (CFD) is used to model a boosted DISI engine with a focus on knock prediction and fuel property effects. A level-set G-equation model is employed to capture turbulent premixed combustion, and is coupled with a transported Livengood-Wu (L-W) integral approach to predict autoignition in the unburnt region. A criterion associated with the L-W integral is developed to accurately predict knock onset and knock-limited spark-advance. This model is then applied to a sensitivity analysis of HoV and LFS on knock tendency and thermal efficiency. The pressure-temperature trajectory framework is applied and extended to study the fuel effects on auto-ignition process in the engine. An existing efficiency-based merit function, which is derived from experiments for boosted SI engines, is evaluated and improved based on the current CFD results.

DISI↗
DOE OSTI2022

Investigation of Stochastic Unit Commitment to Enable Advanced Flexibility Measures for High Shares of Solar PV

As the share of solar photovoltaics (PV) in the power system increases, there is a growing need for flexibility from multiple, possibly interdependent sources to adjust to PV's variability, uncertainty, and diurnal dependence. This paper investigates how stochastic unit commitment leveraging probabilistic solar forecasts can support other flexibility measures under high solar shares. We consider two flexibility measures relevant to day-ahead scheduling: battery energy time-shifting and solar ancillary service provision. Unit commitment and economic dispatch simulations are conducted on a realistic test system based on Texas using day-ahead solar trajectories. The benefits of the two flexibility measures are pronounced when the instantaneous solar share is high, offering cost savings of 10%-20% in the spring. For a Texas-sized system, this translates to hundreds of millions of dollars in cost savings once the installed PV capacity enables instantaneous solar shares regularly exceeding 40%. Using probabilistic forecasts also greatly increases the reliability of upward reserve provision from solar PV, reducing unserved reserves by 50%-100%. Both day-ahead forecast resolution and errors can impact system reliability at high solar shares, but the stochastic formulation has significant value, mitigating reliability impacts on over-forecast days.

ancillary services↗
DOE OSTI2022

Optimal management of electric hotel loads in mild hybrid heavy duty truck

Here, the problem of engine idling for heavy-duty trucks has been under study for decades with Auxiliary Power Units (APUs) and Truck Stop Electrification (TSE) as the most compelling solutions. With the electrification of trucks approaching feasibility in terms of cost-effective technology, hybridization offers another “degree of freedom” to tackle the problem. Here, this work aims at exploiting a battery pack of a 48 V mild-hybrid heavy-duty truck to store sufficient onboard energy for powering the auxiliary loads during the hoteling. This problem is not trivial, as the battery packs typically cannot recover the entire energy required through regeneration alone; hence an optimal energy management strategy needs to be employed to charge the battery through the engine during drive operation. This strategy optimizes powertrain performances among the four modes: (i) Engine of Coasting (EOC), (ii) Regeneration by braking, (iii) Regeneration by engine, and (iv) engine idling. This paper presents the development of a Dynamic Programming (DP) framework that employs a multi-objective cost function to minimize the fuel consumption and maximize the regeneration using the above-mentioned four modes. A typical heavy-duty truck drive cycle is used to represent the drive phase, with mandatory hoteling stops as per regulations. A comprehensive powertrain model is developed using validated components’ model. The DP employs two state variables: battery State-of-Charge (SOC) and engine mode, and three control inputs: (i) the engine ON–OFF state, (ii) clutch engagement state, and (iii) power request at the Electric Machine (EM) for calculating optimal SOC trajectory. The framework also tackles rapid engine ON–OFF scenarios to avoid the challenges associated with DP and the compromises in fuel cost with those approaches. Finally, the effectiveness of the proposed framework is tested for potential fuel savings on two different battery packs by performing the full cycle simulations. The results show 6.47% of fuel consumption reduction as compared to traditional APU-based heavy-duty truck.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Heterogeneity or illusion? Track the carbon Kuznets curve of global residential building operations

Residential buildings, the “last mile” sector in the global decarbonization, have become the most significant uncertain factor hindering carbon neutrality with increasing household energy demand. To track the operational carbon in buildings, this study investigates the carbon Kuznets curve (CKC) and the corresponding decoupling status of residential building operations at four emission scales by using the data of 30 countries from 2000 to 2019. The results show that (1) the CKC model can fit more than half of the samples. Most curves have an inverted U-shape, with 76% of emission per household and 82% of total emissions. (2) In the presence of the CKC, over four-fifths of global residential buildings peak regardless of any emission scale. The analysis denotes that the carbon emissions of developed countries reach their peaks earlier. In the total emissions, the samples’ peaking proportion is 20% and 25% with income per capita < 20,000 United States dollars (USD) and 20,000–40,000 USD, respectively. (3) The Tapio decoupling analysis and the threshold regression effectively verify the robustness and the heterogeneity of CKCs, respectively. Strong decoupling effects of CKCs in most countries are demonstrated at the scales of emission per floor space and the total emissions, and the heterogeneity proves the classic inverted U-shaped relationship between economy and emissions doesn’t exist in all emitters. Overall, this study tracks the historical carbon emission trajectories of residential building operations at a global scale, providing reference for different economies to simulate the dynamic of building carbon emissions along with the economic booming.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Survey-wide asteroid discovery with a high-performance computing enabled non-linear digital tracking framework

Modern astronomical surveys detect asteroids by linking together their appearances across multiple images taken over time. This approach faces limitations in detecting faint asteroids and handling the computational complexity of trajectory linking. Here, we present a novel method that adapts “digital tracking” – traditionally used for short-term linear asteroid motion across images – to work with large-scale synoptic surveys such as the Vera Rubin Observatory Legacy Survey of Space and Time (Rubin/LSST). Our approach combines hundreds of sparse observations of individual asteroids across their non-linear orbital paths to enhance detection sensitivity by several magnitudes. To address the computational challenges of processing massive data sets and dense orbital phase spaces, we developed a specialized high-performance computing architecture. We demonstrate the effectiveness of our method through experiments that take advantage of the extensive computational resources at Lawrence Livermore National Laboratory. This work enables the detection of significantly fainter asteroids in existing and future survey data, potentially increasing the observable asteroid population by orders of magnitude across different orbital families, from near-Earth objects (NEOs) to Kuiper belt objects (KBOs).

Asteroid discovery↗
DOE OSTI2026

Atmospheric pion, kaon, and muon fluxes for sub-orbital experiments

Cosmic rays interacting with the Earth's atmosphere generate extensive air showers, which produce Cherenkov, fluorescence and radio emissions. These emissions are key signatures for detection by ground-based, sub-orbital, and satellite-based telescopes aiming to study high energy cosmic ray and neutrino events. However, detectors operating at ground and balloon altitudes are also exposed to a background of atmospheric charged particles, primarily pions, kaons, and muons, that can mimic or obscure the signals from astrophysical sources. In this work, we use coupled cascade equations to calculate the atmospheric pion, kaon and muon fluxes reaching detectors at various altitudes. Our analysis focuses on energies above 10 GeV, where the influence of the Earth's magnetic field on particle trajectories is minimal. We provide angular and energy-resolved flux estimates and discuss their relevance as background for extensive air shower detection. Furthermore, our results are potentially relevant for interpreting data from current and future balloon-borne experiments such as EUSO-SPB2 and for refining trigger and veto strategies in Cherenkov and fluorescence telescopes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Atmospheric ice nuclei concentration measurements over a high altitude-station in the Western Ghats, India

Mixed-phase and ice clouds play a significant role in modulating the Indian summer monsoon rainfall, and modeling studies suggest that inaccurate representation of ice nucleating particles (INP) is one of the factors, which have resulted in a deficiency in the prediction of cloud properties and hydrological cycle. We do not understand the ice nucleation mechanisms to represent the ice formation in the model, and therefore more INP measurements at various supercooled temperatures are needed. This study reports for the first time real-time INP concentration at a high altitude station in the Western Ghats region of India, and these results will be used to constrain the ice nucleation parameterizations in the model. The variation of INP for any particular cloud condition can have strong microphysical responses of mixed-phase precipitation formation processes. Also, to improve the representation of heterogeneous ice nucleation, INP measurements of various types of aerosols commonly observed over the Indian region are needed. This study presents, the INP concentrations and daily concentration of non-refractory chemical composition (organics, sulfate, nitrate, ammonium, and chloride), cloud condensation nuclei (CCN), and aerosol size distribution. Spectrometer for Ice Nuclei (SPIN) was operated to measure INP concentration at the mountain site from August to December 2018. INP measurements were performed at three different temperatures (-25, -30, and - 34 °C) in the immersion freezing mode relevant for mixed-phase cloud conditions. The average INP concentrations approximately varied from 0.18 to 12.4 L-1, 0.39 to 24 L-1, 1.1 to 40.2 L-1 at -25, -30, and - 34 °C, respectively. These concentrations are within the concentrations reported in the literature from different regions across the globe. The air mass back trajectory analysis showed that continental air mass contains more INP compared to maritime air mass. The information on INP concentration and chemical characteristics of aerosol will help to improve heterogeneous ice nucleation parameterization in numerical models.

Kumar, Anil V.↗
DOE OSTI2020

Intra-annual variations of regional total column ozone, aerosol optical depth, and water vapor from ground-based, satellite-based and model-based observations

The Microtops II ozonometer has been deployed to continuously measure the total column ozone (TCO), pre-cipitable water content (PWC) and aerosol optical thickness (AOT) at Atigre village (16.74° N, 74.37° E, 604 m above sea level, masl) located on the southeastern slope of Indian western Ghats. The Microtops II ozonometer measurements during the period from September 2017 to June 2018 were analyzed along with the retrieved products of Ozone Monitoring Instrument (OMI) and Moderate Resolution Imaging Spectroradiometer (MODIS) onboard NASA's Terra satellites and also the ERA-interim reanalysis model. We found that the TCO (or PWC) data from the OMI (or MODIS) and the ERA reanalysis model data products are in good agreement. The mea-surements showed strong seasonal variability such that TCO was lowest in the winter season but increased in the summer season, and both PWC and AOT were lowest in the post-monsoon season and increased in the summer season. In continuous ERA-interim reanalysis TCO observations, we found the three types of periodicities in the daily mean TCO as well as PWC (i) The weekly (7 days), (ii) quasi-biweekly (14 to 16 days), (iii) Madden Julian Oscillations (MJO) (30–60 days) oscillations. These oscillations depend on the periodic weather changes in the troposphere. We also observe that the average day to day variability in daily mean TCO observations in the observational period is about 1.4%, 1.1%, 2.62% in the Microtops II ozonometer, ERA-interim reanalysis, and OMI data respectively. We found the positive correlation between daily mean AOT and PWC in all seasons (except monsoon), which is maximum in the winter season, showing the hygroscopic nature of aerosols. The sources of water vapor and aerosol at our location are also studied using back-trajectory analysis.

Nade, Dada R.↗
DOE OSTI2021

Dynamic consensus and extended high gain observers as a tool to achieve practical frequency synchronization in power systems under unknown time-varying power demand

In this paper, we consider the frequency synchronization problem in a network of lossless, connected, and network-reduced power system. Frequency synchronization is an important problem in power systems as frequency deviation can lead to degraded power quality, tripping of generators, etc. We present a load-estimator-based consensus algorithm in a network of fourth-order generator models that achieves practical frequency synchronization in the presence of unknown time-varying power demand. We do not assume that the demand is generated by a known model. We follow a multiple time-scale approach for the stability analysis. It is shown that synchronization is achieved for the case of constant power demand, by showing that the system is exponentially stable. Furthermore, for the case of time-varying power demand, we show the error trajectories are ultimately bounded, and the synchronization error is proportional to an upper bound on the derivative of the time-varying power demand. We simulate the controller performance on a four area network which is equivalent to the IEEE New England 39-Bus system.

42 ENGINEERING↗
DOE OSTI2026

Higher-order interaction effects among operating conditions and feedstocks shape reactor microbiomes and fatty acid production profiles

Arrested anaerobic digestion (AAD) offers a promising route for producing fatty acids (FAs) from organic residues, yet optimal conditions for selectively generating medium-chain fatty acids (MCFAs) remain poorly defined. Here, we systematically evaluated the main and interaction effects of pH (5, 7, 9), feedstock (food waste, manure), temperature (35 and 45 °C), and inoculum source on microbiome composition and FA production. Anaerobic digester sludge and a novel bison rumen inoculum were compared. Significant higher-order interactions among operating parameters governed FA profiles and microbiome structure. Butyric acid production was driven by a three-way interaction among pH, feedstock, and temperature (p < 0.001), with maximum concentrations achieved in food waste reactors at pH 5.0 and 35 °C (1.2 ± 0.1 g L −1 with sludge and 1.1 ± 0.3 g L −1 with rumen). MCFA production exhibited significant four-way interactions (p < 0.1 to p < 0.001). At 45 °C and pH 5.0, inoculum source tuned MCFA selectivity: sludge favored pentanoic acid (0.4 ± 0.1 g L −1 ), whereas rumen favored hexanoic and heptanoic acids (up to 0.4 ± 0.2 g L −1 ). Manure reactors produced < 0.2 g L −1 MCFAs under all conditions. Genera, including Megasphaera, Prevotella, and Lactobacillus, were associated with production of specific MCFAs. PICRUSt2-based pathway predictions were consistent with MCFA production patterns and suggested a potential role for lactic acid–driven chain elongation pathways. This study provides insights into how interacting operating conditions shape AAD microbiomes, their FA profiles, and advances the trajectory of research aimed at engineering robust and controllable microbiomes for waste valorization.

09 BIOMASS FUELS↗
DOE OSTI2025

Data-driven equation-free dynamics applied to many-protein complexes: The microtubule tip relaxation

Microtubules (MTs) constitute the largest components of the eukaryotic cytoskeleton and play crucial roles in various cellular processes, including mitosis and intracellular transport. The property allowing MTs to cater to such diverse roles is attributed to dynamic instability, which is coupled to the hydrolysis of GTP (guanosine-5'-triphosphate) to GDP (guanosine-5'-diphosphate) within the β-tubulin monomers. Understanding the equilibrium dynamics and the structural features of both GDP- and GTP-complexed MT tips, especially at an all-atom level, remains challenging for both experimental and computational methods because of their dynamic nature and the prohibitive computational demands of simulating large, many-protein systems. This study employs the “equation-free” multiscale computational method to accelerate the relaxation of all-atom simulations of MT tips toward their putative equilibrium conformation. Using large MT lattice systems (14 protofilaments × 8 heterodimers) comprising ~21-38 million atoms, we applied this multiscale approach to leapfrog through time and nearly double the computational efficiency in realizing relaxed all-atom conformations of GDP- and GTP-complexed MT tips. Commencing from an initial 4 μs unbiased all-atom simulation, we interleave coarse projective “equation-free” jumps with short bursts of all-atom molecular dynamics simulation to realize an additional effective simulation time of 1.875 μs. Our 5.875 μs of effective simulation trajectories for each system expose the subtle yet essential differences in the structures of MT tips as a function of whether β-tubulin monomer is complexed with GDP or GTP, as well as the lateral interactions within the MT tip, offering a refined understanding of features underlying MT dynamic instability. Furthermore, the approach presents a robust and generalizable framework for future explorations of large biomolecular systems at atomic resolution.

Wu, Jiangbo [University of Chicago, IL (United Sta↗
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 OSTI2023

Operator inference with roll outs for learning reduced models from scarce and low-quality data

Data-driven modeling has become a key building block in computational science and engineering. However, data that are available in science and engineering are typically scarce, often polluted with noise and affected by measurement errors and other perturbations, which makes learning the dynamics of systems challenging. Here, in this work, we propose to combine data-driven modeling via operator inference with the dynamic training via roll outs of neural ordinary differential equations. Operator inference with roll outs inherits interpretability, scalability, and structure preservation of traditional operator inference while leveraging the dynamic training via roll outs over multiple time steps to increase stability and robustness for learning from low-quality and noisy data. Numerical experiments with data describing shallow water waves and surface quasi-geostrophic dynamics demonstrate that operator inference with roll outs provides predictive models from training trajectories even if data are sampled sparsely in time and polluted with noise of up to 10%.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A novel method for quantifying irradiation damage in nuclear graphite using Raman spectroscopy

Raman spectroscopy has long been used in studying irradiation damage in carbon/graphite materials. It is however unclear if the measurements from different types of materials are directly comparable. Further, decoupling the contribution from irradiation temperature and fluence to the total damage in nuclear graphite possessing irradiation damage gradient is currently not readily feasible as it requires a complete set of test reactor irradiation experiments over the relevant temperature and fluence range. A novel methodology has therefore been proposed and developed to quantify total irradiation damage evolution at crystal level based on graphite Raman G-band position shift. Specifically, G-band positions derived from a large number of spectra collected from the fractured surface of microfine-grained POCO ZXF-5Q graphite possessing proton irradiation damage gradient were plotted with open literature Raman data on HOPG, BEPO (AXGP graphite), PCEA and IG-110 graphite as a function of dpa. The total damage level within 2σ beam radius in this POCO graphite was estimated to be equivalent to ~2–5 dpa at ~350–370 °C. Derived G-band positions were then mapped to the three-stage amorphization trajectory model indicating beam centre area has entered the second stage, i.e., transitioning from nanocrystalline graphite into amorphous carbon. G-band position relative shift (ΔG) curve with a ‘turn-around’ peak as a function of total damage can be used for in-service lifetime prediction. Importantly, the developed methodology has the potential to ‘unify’ total damage levels across different grades of nuclear graphite caused by ion, neutron and proton irradiation at different temperatures.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Computationally guided synthesis of carbon coated mesoporous silica materials

Mesoporous silica materials (MSMs) are well-suited for biomedical applications due to their unique features, including a large surface area and tunable pore size. To enhance their durability, the small pores in MSMs are filled with carbon precursors and then carbonized to prevent them from interacting with unreacted silicic acid. Here, in this study, we synthesized and healed MSMs using a combination of non-reactive and reactive molecular dynamics (MD) simulations. The non-reactive MD simulation revealed that the self-assembly of Pluronic® L64 polymers in water resulted in nearly 80 % hydrogen bonds between the hydrophilic sections of the micelle and water. In the bond-boosted ReaxFF MD simulations, silicic acid precursors were condensed on the micelle surface, with over 60 % of them leading to the creation of periodic mesoporous silica within the system. Condensation of silicic acid precursors at 300 K with bond-boosting and at 1500 K without it both significantly promoted the polymerization of Si(OH) 4 , with the latter doubling the rate compared to the former. Subsequently, we healed the MSM surface by carbonizing carbon precursors inside an MSM pore. Polyethylene (PE) and high-rank lignite were identified as the most suitable precursors due to their ability to form turbostratic graphene structures. High-rank lignite exhibited the highest carbon conversion ratio to 6-membered rings, closely followed by PE, in the carbon ring formation analysis. Additionally, the production of gases, such as H 2 , increased significantly for PE at both 2200 K and 2600 K, indicating the conversion of a considerable portion of carbon into graphitic or turbostratic structures. The carbonization of PE primarily led to the formation of planar (sp 2 ) structures, while sucrose yielded the least planar structures. Finally, we studied the protective blocking of unreacted silicic acid precursor by considering a PET turbostratic graphene structure in a silica mesopore formed at 2600 K. The trajectory analysis showed that the surface of the silica was effectively coated with PET tar, preventing unreacted silicic acid from interacting with the inner silica pore surface. These findings offer valuable insights into the synthesis and carbonization-based healing processes of MSMs, enhancing their potential for various biomedical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A novel xylosylated fucoglucuronan in Penium reveals structural parallels to rhamnogalacturonan-I and its broad evolutionary footprint in lower plants

Green algae inhabit aquatic environments across the planet and play a crucial role in sustaining the global ecosystem. Ancestors of some Charophytes adapted to terrestrial conditions and eventually evolved into land plants. Extant green algae have inherited traits from their ancestors and evolved into their current morphological and chemical forms, as reflected by their cell walls with distinct shapes and compositions. To illuminate the evolution of plant cell walls and bridge the gap between green algae and land plants, we investigated the charophyte Penium margaritaceum, a close relative of terrestrial plants. We discovered a previously unknown polysaccharide in both its culture medium and cell wall. This polysaccharide, termed xylosylated fucoglucuronan (XFG), possesses a rhamnogalacturonan-I (RG-I)-like backbone composed of repeating [-3-α-Fucp-(1,4)-α-GlcpA-] disaccharides that are extensively xylosylated and acetylated. Surveying approximately 20 non-vascular plants revealed that XFG and RG-I (or related structures) first emerge in certain Chlorophyceae and subsequently co-occur throughout lineages along the evolutionary trajectory to bryophytes, thereby bridging aquatic green algae to early land plants. The striking structural parallels between XFG, RG-I, and ulvan suggest a shared evolutionary origin, offering new insight into how plant cell walls adapted during the transition from marine to freshwater environments and ultimately to land.

Algae↗