Search NASA⌕ Search

SEARCH · Search NASA

Results for “Methods”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A High-Order Eulerian–Lagrangian Runge–Kutta Finite Volume (EL–RK–FV) Method for Scalar Nonlinear Conservation Laws

Abstract We present a class of high-order Eulerian–Lagrangian Runge–Kutta finite volume methods that can numerically solve Burgers’ equation with shock formations, which could be extended to general scalar conservation laws. Eulerian–Lagrangian (EL) and semi-Lagrangian (SL) methods have recently seen increased development and have become a staple for allowing large time-stepping sizes. Yet, maintaining relatively large time-stepping sizes post shock formation remains quite challenging. Our proposed scheme integrates the partial differential equation on a space-time region partitioned by linear approximations to the characteristics determined by the Rankine–Hugoniot jump condition. We trace the characteristics forward in time and present a merging procedure for the mesh cells to handle intersecting characteristics due to shocks. Following this partitioning, we write the equation in a time-differential form and evolve with Runge–Kutta methods in a method-of-lines fashion. High-resolution methods such as ENO and WENO-AO schemes are used for spatial reconstruction. Extension to higher dimensions is done via dimensional splitting. Numerical experiments demonstrate our scheme’s high-order accuracy and ability to sharply capture post-shock solutions with large time-stepping sizes.

Chen, Jiajie↗
DOE OSTI2025

Further Development of the Tamped Richtmyer-Meshkov Instability Method and Application to Molybdenum Dynamic Strength Calibration and Tabulation

The high pressure and high strain rate dynamic strength of Mo is experimentally and computationally investigated in the 3–20 GPa stress, 50–600 C temperature, and 10 5 –10 6 /s strain rate regimes using a modified tamped Richtmyer-Meshkov instability (RMI) method. Modifications to the tamped RMI method include a method to determine loading states during strain, a new strength calibration function based on interface shape, and a robust uncertainty quantification method. These modifications improve fidelity of the tamped RMI method, allowing evaluation of the compensating effects of pressure hardening, strain rate hardening, strain hardening, and thermal softening. The new calibration function based on interface shape is not limited to sinusoidal corrugations and could be applied to additional interface shapes. Plate impact experiments are performed at Argonne National Laboratory’s Advanced Photon Source’s Dynamic Compression Sector operated by Washington State University (DCS), driving a planar shock front through a corrugated Mo-D 2 O or Mo-C 8 F 18 interface, forcing the corrugation to significantly deform. The extent of interfacial deformation, RMI growth, is experimentally observed using X-ray phase contrast imaging at the DCS. RMI jet lengths and jet shapes are extracted from the experimental radiographs, then used to calibrate numerical simulations performed with the Sandia National Laboratories (SNL) hydrocode CTH. Mo yield strength, Y, as a function of shock pressure, P, strain rate, $\dot{\varepsilon }$, accumulated strain, ϵ, relative volumetric compression, RD , and temperature, T , is determined for each impact experiment and presented. The calibrated Mo yield strength values range 1.2–1.8 GPa, with strength generally decreasing as the impact stress increases. This trend is likely caused by thermal softening or strain localization. The tabular yield strength versus loading condition data presented in this paper can be used to fit complex strength models.

Voorhees, T. J. [Sandia National Lab. (SNL-CA), Li↗
DOE OSTI2026

Optimizing bioenergy biofuel harvest: a comparative analysis of stepwise and integrated methods for economic and environmental sustainability

Switchgrass is a promising bioenergy feedstock due to its high biomass yield potential, adaptability to marginal lands, and low carbon intensity for feedstock production. However, accurate cost estimation and assessment of greenhouse gas (GHG) emissions for the energy-intensive harvesting process are essential for evaluating the sustainability of bioenergy. This study provides a comparative analysis of two harvesting methods: the Stepwise Method, which separates operations into multiple stages, and the Integrated Method, which combines mowing and raking into a single pass. The analysis was conducted under four scenarios based on field sizes and biomass yields. Using three years of field-scale switchgrass harvest data from 125 sites, GHG emissions, energy consumption, and harvesting costs were quantified using the GREET model and techno-economic analysis. Additionally, regression analysis identified key climate and operational factors affecting fuel consumption. The Stepwise method was the most cost-effective for large fields with high biomass yield, achieving the lowest harvesting costs ($37.70 per ton). In contrast, the Integrated Method performed better in small fields and low-yield conditions, reducing GHG emissions by 9 % and energy use by 5 %. Regression analysis confirmed that a larger field size reduced fuel consumption, while higher biomass yield and longer operational time increased fuel use. Maximum temperature also contributed to a slight increase in fuel consumption. Furthermore, these results provide actionable insights for optimizing harvesting strategies based on field-specific conditions and operational goals, contributing to the economic and environmental sustainability of bioenergy production.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Dynamic flux surrogate-based partitioned methods for interface problems

Loosely coupled partitioned methods for multiphysics problems treat each subproblem as a separate entity and advance them independently in time. In so doing these methods enable code reuse, increase concurrency and provide a convenient framework for plug-and-play multiphysics simulations. However, mathematically loosely coupled schemes are equivalent to a single step of an iterative solution method, which can compromise their accuracy and stability. We present a new data-driven partitioned method for coupled parametric PDEs that can improve upon the accuracy of traditional loosely coupled methods without incurring a performance penalty. To that end, we replace conventional field transfers across the interface by a surrogate for the dynamics of the interface flux exchanged between the subdomains. To develop this surrogate we apply dynamic mode decomposition to a non-standard staggered-in-time state, comprising the interface flux and small solution patches near the interface. The new approach shifts the main computational burden to an offline training phase, whereas application of the surrogate in the online phase amounts to a single matrix–vector multiplication. In conclusion, we provide stability analysis of the surrogate-based partitioned scheme and include numerical results that demonstrate its potential.

Dynamic mode decomposition (DMD)↗
DOE OSTI2025

Multi-material ALE remap with interface sharpening using high-order matrix-free finite element methods

The arbitrary Lagrangian-Eulerian (ALE) technique involves remapping field quantities from a Lagrangian mesh to an optimized mesh in a conservative, accurate and bounds-preserving manner. For methods based on arbitrary order finite elements, as described in a reference, material volume fractions are advected in pseudo-time using flux-corrected transport (FCT) without any form of interface reconstruction. In practice, this can lead to excessive propagation of small volume fractions throughout the domain. In addition, this method requires assembly of a global advection matrix to compute the bounds-preserving low-order FCT solution. In this work, we introduce a new approach for ALE remap using a high-order matrix-free technique which incorporates a flux modification to sharpen material interfaces in a conservative manner. Our approach begins with computing a bounds-preserving low-order solution to the ALE remap equations at the element level. We then compute a sharp interface solution (not guaranteed to be bounds-preserving) which comes from solving an augmented version of the ALE remap equations with a conservative flux modification which acts to sharpen material volume fractions based on their gradients and transport directions. Using the sharp interface solution, we make global corrections to the bounds-preserving solution while maintaining preservation of bounds. By blending with the sharpened solution at the global level we are able to globally conserve mass without hindering the remap pseudo-time step. This new interface-aware ALE remap method is based entirely on partial assembly techniques where globally assembled matrix operators are no longer needed, resulting in a globally matrix-free FCT method for multi-material, multi-field ALE remap with high performance on GPU architectures. We present results of our new remap method on 1D, 2D and 3D benchmarks and describe the algorithmic tailoring for GPU architectures that was developed.

Vargas, Arturo [Lawrence Livermore National Labora↗
DOE OSTI2025

Performance evaluation of automated data-driven feature extraction and selection methods for practical and scalable building energy consumption prediction models

Here, this study quantifies the impact of automated feature engineering methods (feature extraction and selection) on the quality and accuracy of machine learning models that predict building energy consumption. The case study compares model performance for three main scenarios: baseline (no feature extraction and selection), feature extraction only, and feature extraction combined with feature selection (filter and/or wrapper methods) for fully trained machine learning models for 200 metered/sub-metered energy measurements across 118 real buildings. For consistency, the same machine learning model architecture (a black box deep learning neural network with probabilistic forecast output) was used for all scenarios. Based on results, all feature engineering methods provided noticeable prediction accuracy improvements (e.g., 29%-68% median prediction improvement) compared to baseline scenarios. However, in this application, feature selection methods provide little practical value due to their limited performance gains and high computational cost. Smarter algorithm development supported by better computational environments will be needed before feature selection methods can reliably and efficiently improve predictive model performance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Collocation methods for nonlinear differential equations on low-rank manifolds

We introduce new methods for integrating nonlinear differential equations on low-rank manifolds. These methods rely on interpolatory projections onto the tangent space, enabling low-rank time integration of vector fields that can be evaluated entry-wise. A key advantage of our approach is that it does not require the vector field to exhibit low-rank structure, thereby overcoming significant limitations of traditional dynamical low-rank methods based on orthogonal projection. To construct the interpolatory projectors, we develop a sparse tensor sampling algorithm based on the discrete empirical interpolation method (DEIM) that parameterizes tensor train manifolds and their tangent spaces with cross interpolation. Using these projectors, we propose two time integration schemes on low-rank tensor train manifolds. The first scheme integrates the solution at selected interpolation indices and constructs the solution with cross interpolation. The second scheme generalizes the well-known orthogonal projector-splitting integrator to interpolatory projectors. We demonstrate the proposed methods with applications to several tensor differential equations arising from the discretization of partial differential equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Evaluation of a high-throughput method for processing sponge-stick samples to detect viable, non-spore-forming biothreat agents

After a bioterrorism incident, surface sampling is often used to determine the extent of contamination and exposure, guiding decontamination efforts and decisions for re-occupancy of affected sites. The sponge-stick (SS) is a preferred and commonly used device for sample collection to detect both spore-forming and non-spore-forming biothreat agents from non-porous surfaces. Here, in this study, a recently developed high-throughput method (HTM) for processing SS samples to detect viable Bacillus anthracis spores was adapted for detection of non-spore-forming biothreat agents, Yersinia pestis and Francisella tularensis. The scalable HTM was used to process up to 20 SS samples simultaneously, compared to the current stomacher-based method which processes one SS at a time. Comparisons of the HTM and the stomacher-based method were statistically indistinguishable for most experiments (P > 0.05) with HTM recoveries of 37–60 % for Y. pestis inoculated at 102–103 cells/SS and held 48 h at 4 °C to mimic sample transport/storage. The HTM was integrated with Rapid Viability-Polymerase Chain Reaction (RV-PCR) analysis to detect viable Y. pestis in the presence of particulate contamination (Arizona Test Dust, ATD). This approach detected Y. pestis inoculated at 20 cells/SS and ATD did not impact detection (P > 0.05). F. tularensis showed significantly lower recoveries between no-hold time and 48-h hold time (4 °C, P < 0.05) using the HTM, which further testing showed could be due to toxicity of the neutralizing buffer used for SS pre-wetting. With modifications, this method could enhance throughput capacity while maintaining similar recovery efficiencies to current methods for other non-spore-forming bacterial pathogens.

Biological and medical sciences↗
DOE OSTI2026

A polarimetry-based field-deployable non-interruptive mirror soiling detection method

The soiling level of heliostat mirrors in Concentrated Solar Power (CSP) fields is one of the key factors that significantly influences optical efficiency. State-of-the-art methods of monitoring heliostats soiling levels still face various challenges, including slow speed, labor-intensive operations, resolution and accuracy constraints or interruptions to solar field operations. Here, we present a rapid, cost-effective, and non-intrusive method for mirror soiling detection based on polarimetric imaging, referred to as Polarimetric Imaging-based Mirror Soiling (PIMS). The compact PIMS device is designed for integration with unmanned aerial vehicles (UAVs), enabling rapid, large-area assessments of heliostat mirrors for efficient soiling detection. Our method utilizes the correlation between the Degree of Linear Polarization (DoLP) and surface soiling level based on Mie scattering theory and Monte Carlo simulations. Field deployment of the PIMS method requires minimal device installation, and its UAV-based operation allows for soiling detection without interrupting plant activities. The PIMS method holds the potential for mirror soiling detection across various concentrated solar power (CSP) plants and can be further adapted for other types of solar fields, such as parabolic trough systems.

CSP Field↗
DOE OSTI2024

Quantifying Phospholipids in Organic Samples Using a Hydrophilic Interaction Liquid Chromatography–Inductively Coupled Plasma High-Resolution Mass Spectrometry (HILIC-ICP-HRMS) Method

Here, in this study, a novel method using hydrophilic interaction liquid chromatography (HILIC) coupled with inductively coupled plasma high-resolution mass spectrometry (ICP-HRMS) was introduced for the quantification of phospholipids in oil samples. The method employed a bridged ethyl hybrid (BEH) stationary phase HILIC column with a tetrahydrofuran (THF)/water mobile phase, enhancing the solubility and detection of phospholipids. During the study, a gradient/matrix effect on ICP-HRMS sensitivity was observed and successfully compensated for experimentally, ensuring reliable quantification results. This approach has proven effective for a wide range of different oil samples including vegetable oils, animal fats, and phospholipid supplements. Notably, this method allowed the direct quantification of phospholipids in oil samples, bypassing the need for prior sample preparation methods, such as solid phase extraction (SPE), thereby streamlining the analytical process. The precision, accuracy, and reduced need for extensive sample preparation offered by this method mark a significant advancement in lipids analysis. Its robustness and broad applicability have substantial implications for industries such as food and renewable energy production, where both efficient and accurate lipid identification and quantification are crucial.

09 BIOMASS FUELS↗
DOE OSTI2026

Multireference Methods for Chemistry and Materials Science: Automated Active Spaces, Efficient Dynamic Correlation, and Extended Systems

While multiconfigurational approaches have long been relegated to expert practitioners working on a case-by-case basis, recent developments have increasingly made these methods more routine and applicable to broader sets of systems. This article outlines the state-of-the-art in multiconfigurational approaches, with an emphasis on moving from delicate hand-selected pathways through configuration space toward more robust and efficient approaches to treating a host of challenging chemical systems accurately. First, we overview recent work in automated active-space selection, which has enabled increasingly large-scale applications of multireference methods to modeling vertical excitations and reactivity. Second, we highlight the increasingly efficient methods for recovering correlation energy beyond the active space, as headlined by extensions of pair-density functional theory and its role in accurate and efficient treatment of excited-state dynamics and its utilization to train machine-learned potentials. Finally, we highlight recent efforts to treat extended systems that until recently have lied beyond the traditional limits of active-space methods, giving center stage to product-form wave functions of the localized active space family of methods that allow for the computation of multiconfigurational band structures. These recent advancements point to a broader use of multireference approaches for high-impact chemical and materials science applications.

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

A Unified Analytical Method Greenness Score ( uAMGS ) Quantifies How Microscopic Imaging Is Greener Than Conventional Liquid Chromatography

Green chemistry is a set of principles for assessing, developing, and implementing methods that are safer, more efficient, and less detrimental to the environment. The analytical method greenness score (AMGS) is one of many metrics that attempt to evaluate traditional liquid chromatography (LC) based on the energy consumption of the instrument and the safety, health risks, and environmental impact of the solvents employed. Unfortunately, in practice, the AMGS is primarily focused on traditional separation methods in the pharmaceutical industry and is not amenable to cutting-edge separation science, including miniaturization. To broaden this scope, the unified Analytical Method Greenness Score (uAMGS) is presented here, which clarifies and expands on the underlying mathematics and incorporates both dimensional and uncertainty analysis, enabling its application to a broader range of analytical techniques. The uAMGS is used to compare the greenness of two distinct methods: single-molecule microscopy (SMM) and high-performance liquid chromatography (HPLC), which were used to collect equivalent data. uAMGS determines that SMM is significantly greener than HPLC due primarily to decreased solvent consumption. Overall, the uAMGS should allow chemists ranging from undergraduates to industrial PhDs to assess the greenness of a wide range of separations.

chemical separations↗
DOE OSTI2025

Methods for Carbon Mass Closure in Polyolefin Hydrocracking

Heterogeneous catalytic hydrocracking of polyolefins is a promising approach for the processing of postconsumer plastics, but product quantification methods remain inconsistent across the literature. In systems that generate a large fraction of vapor-phase products, typical product capture methods can result in large carbon balance deficits, exceeding 50%, compromising reported yields and selectivities. Here, we identify the major sources of product loss and develop enhanced capture methods to improve the quantification accuracy. Seven supplemental techniques were evaluated, targeting either increased vapor recovery (by increasing the volatility or system volume) or enhanced retention in the liquid phase (by decreasing volatility). Among these, a flow collection approach using a continuous helium sweep and downstream gas sampling bag capture yielded the highest recovery, achieving a 96 ± 9.2% carbon balance closure. We show that the efficacy of these methods is strongly dependent on product distribution. In general, solvent addition was most effective when condensable species dominate the product distribution, while flow collection was preferred when both condensable species and light gases are present in high concentrations. These results highlight the need for method-specific workup strategies and demonstrate that no single protocol is universally optimal. We provide general guidelines for selecting and implementing robust product capture techniques, enabling accurate yield and selectivity determinations in polyolefin hydrocracking systems.

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

Numerical Investigation of Observational Flux Partitioning Methods for Water Vapor and Carbon Dioxide

Abstract While yearly budgets of CO 2 flux (F c ) and evapotranspiration (ET) above vegetation can be readily obtained from eddy‐covariance measurements, the separate quantification of their soil (respiration and evaporation) and canopy (photosynthesis and transpiration) components remains an elusive yet critical research objective. In this work, we investigate four methods to partition observed total fluxes into soil and plant sources: two new and two existing approaches that are based solely on analysis of conventional high frequency eddy‐covariance (EC) data. The physical validity of the assumptions of all four methods, as well as their performance under different scenarios, are tested with the aid of large‐eddy simulations, which are used to replicate eddy‐covariance field experiments. Our results indicate that canopies with large, exposed soil patches increase the mixing and correlation of scalars; this negatively impacts the performance of the partitioning methods, all of which require some degree of uncorrelatedness between CO 2 and water vapor. In addition, best performances for all partitioning methods were found when all four flux components are non‐negligible, and measurements are collected close to the canopy top. Methods relying on the water‐use efficiency (W) perform better whenWis known a priori, but are shown to be very sensitive to uncertainties in this input variable especially when canopy fluxes dominate. We conclude by showing how the correlation coefficient between CO 2 and water vapor can be used to infer the reliability of differentWparameterizations.

Environmental Sciences & Ecology↗
DOE OSTI2024

Intercomparison of sensible and latent heat flux measurements from combined eddy covariance, energy balance, and Bowen ratio methods above a grassland prairie

We present a comparison of four different methods of measuring sensible (H) and latent (LE) heat fluxes for a year over a mixed grass prairie ecosystem in the Nebraska SandHills [eddy covariance (EC), energy balance/Bowen ratio (EBBR), residual energy (RES), modified Bowen ratio (MBR) methods]. Additionally, we developed a set of quality control criteria for each method and present a simplification to the traditional EBBR setup. Using EC as reference, all methods yielded similar estimates of yearly H (regression slopes (m) ~ 2% from unity; H EC > H EBBR , H RES , and H MBR ). For yearly LE, EBBR and RES yielded similar estimates with EC (m ~ 2% from unity; LE EC < LE EBBR and LE RES ), while a larger bias was found from MBR (m ~ 8% from unity; LE EC > LE MBR ). At shorter time scales (~ hourly), moderate scatter was found about linear regression fits for H between EBBR and EC (R 2 = 0.81), with smaller scatter between RES and MBR, and EC (R 2 = 0.91). For LE, smaller scatter was also measured between EC, and EBBR and RES (R 2 = 0.89 and 0.87, respectively), with the larger scatter between EC and MBR (R 2 = 0.65). This suggests methods other than EC may be well suited to longer-term applications (≥ yearly), but have larger uncertainty on individual measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Accurate energies for ππ* excited states via exchange scaling: the XS-CASSCF method

The state-averaged complete-active space self-consistent field method (SA-CASSCF) is a widely employed electronic structure method used for studying photochemistry and dynamics owing to its ability to provide a reliable description even of complicated cases while still retaining computational efficiency. However, SA-CASSCF suffers from one Achilles heel, related to the description of ionic ππ* excited states, whose energy is often overestimated by 1–2 eV. In light of this challenge, we present the XS-CASSCF method, a new approach based on the idea of exchange scaling (XS) that screens the involved energy terms to improve the excitation energies of singlet ionic ππ* states. First, we illustrate the power of the XS-CASSCF method using hexatriene and para-quinodimethane as examples, showing that it corrects the targeted ionic states while leaving the other states largely unaffected, giving root-mean-square errors (RMSE) below 0.2 eV for the four lowest states in both cases. Subsequently, XS-CASSCF vertical excitation energies are tested against theoretical best estimates for a set of 11 molecules and 56 excited states. XS-CASSCF performs exceptionally well for the ππ* states of hydrocarbons, reducing the RMSE over 21 excitation energies from 0.96 to 0.27 eV. In the challenging subset of molecules with heteroatoms and a larger number of ππ* and nπ* states, we find that improvements can also be obtained, albeit not as pronounced. We conclude with an outlook into more realistic molecular materials focusing on their singlet–triplet (S 1 /T 1 ) gaps, finding that significant improvements can be obtained along the whole range of S 1 /T 1 gaps studied, going from 0.1 eV to more than 1.5 eV. Owing to notable improvements across significant classes of molecules combined with its conceptual simplicity, we believe that XS-CASSCF is a promising addition to the electronic structure toolbox, serving both as a standalone electronic structure method and as a starting point for further correlated treatment.

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

Determining the Minimum Phase‐to‐Phase Gap Distance for Unconventional Transmission Lines Using Numerical Method

With the increasing demand for power transmission, compact, high surge impedance loading (HSIL) high-voltage transmission lines have emerged as a viable solution due to their reduced land acquisition costs and higher power delivery capability. The compactness of a transmission line depends on effective insulation coordination, particularly in determining the phase-to-phase clearance, which is governed by the critical flashover voltage under switching and lightning overvoltage conditions. Traditional methods for phase-to-phase clearance rely on empirical formulas derived from experimental data, which are convenient for conventional high-voltage lines. However, unconventional HSIL lines require a faster and more adaptable evaluation method, as they involve optimized conductor positioning to reduce right-of-way requirements while enhancing natural power loadability. This study presents a simplified numerical approach to determine the minimum phase-to-phase gap, utilizing arc propagation viability curves, and offers an efficient alternative to conventional empirical methods. The proposed method was successfully applied to a 500 kV conventional line as well as 500 and 735 kV unconventional line designs, demonstrating its capability in accurately assessing insulation requirements. Results reveal that the method can support reduced gap clearances while still maintaining reliability, thereby validating its usefulness in optimizing compact transmission line configurations.

Arafat, Easir [The University of Texas at Dallas, ↗
DOE OSTI2025

Phase-based velocity extraction method for photonic Doppler velocimetry with potential higher time resolution

We present an extension of the [Takeda et al., J. Opt. Soc. Am. 72, 156 (1982)] phase extraction method to heterodyne photonic Doppler velocimetry applications. The method yields results equivalent to those obtained by the short-time Fourier transform (STFT), while offering potential improvements in time resolution. Unlike STFT, which relies on window functions, such as the Hamming window, that emphasize central data points and diminish the influence of edges, the extended Takeda method utilizes all data uniformly. This uniform treatment allows for the derivation of empirical equations that directly relate velocity error to the actual time resolution rather than to the local analysis duration. The established equation provides a useful metric for both optimizing hardware configuration and guiding data analysis. Simulation and experimental results confirm that, for a given dataset, specifying a target time resolution yields consistent velocity errors for both methods. These findings underscore the Takeda method’s advantages, particularly its potential higher time resolution and reduced computational burden, making it a valuable tool for high-throughput applications such as laser dynamic compression experiments.

Computer simulation↗