SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32

All sources NASA NTRS DOE OSTI

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↗

DOE OSTI2026

Mott vs Kondo: Influence of various density functional based methods on the Ce isostructural phase transition mechanism

The cerium iso-structural phase transition (γ–α) is dominated by f-electron localization changes that results in a magnetic ordering change and a volume collapse. Generally, these physics are difficult to capture with ab initio and first principles methods. However, previous works have shown various methods to be successful in predicting at least some of the physics of the γ to α phase transition. Therefore, here, we perform a broad survey of density functional based methods across three levels of theory and types of functions [GGA (generalized gradient approximation), MetaGGA, and hybrid functionals] and compare the results, focusing on hydrostatic compression across the phase boundary at zero kelvin. For the methods that best reproduce experimental results, we directly probe the predicted mechanisms and frame the results in the Mott/Kondo debate, assessing how the underlying methods and assumptions of different functionals can assess the physical drivers in the phase transition, providing insight into the governing processes of this unique phase transition.

36 MATERIALS SCIENCE↗

DOE OSTI2023

A Better Method to Calculate Fuel Burnup in Pebble Bed Reactors Using Machine Learning

Burnup measurement is an important step in material control and accountancy (MC&A) at nuclear reactors, and may be done by examining gamma spectra of fuel samples. Traditional approaches rely on known correlations to specific photopeaks (e.g. 137 Cs) and operate via a standard linear regression method. However, the quality of these regression methods is limited even in the best case, and is significantly poorer at short fuel cool-down times, due to the elevated radiation background by short life-time isotopes, and self-shielding effect of the fuel. For practical operation of pebble bed reactors (PBRs), quick measurements (in minutes) and short cooling times (in hours) are required from a safety and security perspective. We investigated the efficacy and performance of machine learning (ML) methods to predict the burnup of the pebble fuel from full gamma spectra (rather than specific discrete photopeaks) and found a full-spectrum ML approach to far outperform baseline regression predictions in all measurement and cooling conditions - including in operational-like measurement conditions. We also performed model and data ablation experiments to determine the relative performance impact of our ML methods' capacity to model data nonlinearities and the inherent additional information in full spectra. Applying our ML methods, we found a number of surprising results, including improved accuracy at shorter fuel cooling times (the opposite of the norm), remarkable robustness to spectrum compression (via rebinning), and competitive burnup predictions even when using background signal only (i.e. explicitly omitting known isotope photopeaks).

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

DOE OSTI2021

Adjoint Klein-Nishina Sampling Methods: Efficiency, Speed and Applications

Three new rejection sampling methods for generating samples from the adjoint Klein-Nishina cross section are discussed: the two-branch rejection sampling procedure, the three-branch linear rejection sampling procedure and the three-branch inverse rejection sampling procedure. These methods have all been implemented in the Framework for REsearch in Nuclear ScIence and Engineering (FRENSIE). The efficiency and sample generation rate of each of these methods are evaluated to characterize the methods and to make recommendations regarding their use. The use of these methods in realistic transport simulations is also evaluated by incorporating a scattering function into the sampling process. Furthermore, the results of an infinite medium problem are presented to verify that the sampling procedure can be used in an adjoint Monte Carlo simulation to generate results that are in agreement with an equivalent forward simulation.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

High-accuracy method for modeling nucleation and growth of particles

State-of-the-art numerical models describing the kinetics of aerosol particle nucleation and growth from a cooling vapor primarily use a nodal method, in which particles that are smaller than the critical size are omitted from consideration because they are thermodynamically unfavorable. This omission is based on the assumption that most newly formed particles are above the critical size, so that subcritical-size particles are not important to take into account. Due to the nature of the nodal method, it suffers from numerical diffusion, which can cause an artificial broadening of the cluster size distribution leading to a significant overestimation of the number of large-size particles. To address these issues, we propose a more accurate numerical method that explicitly models particles of all sizes, and uses a special numerical scheme that substantially reduces the numerical diffusion and provides high solution accuracy and numerical stability. We extensively compare this novel method to the commonly used nodal solver of the general dynamic equation (GDE) for particle growth and demonstrate that it offers GDE solutions with higher accuracy with low numerical diffusion. Incorporating small subcritical clusters into the solution is crucial for: 1) more precise determination of the entire particle size distribution function and 2) wider applicability of the model to experimental studies with non-monotonic temperature variations leading to particle evaporation. The computational code implementing this numerical method in Python is available upon request.

42 ENGINEERING↗

DOE OSTI2025

A transient near to far field transformation method and verification benchmarking procedure

The numerical calculation of electromagnetic far fields in the time-domain requires a near to far field transformation (NTFF) method. While time-domain NTFF methods for popular finite-difference time-domain (FDTD) approaches are well established, there is little discourse on NTFF methods for finite-element time-domain (FETD) codes. Here, this work is concerned with the development of an NTFF method for the Empire FETD code, which utilizes curl and divergence conforming elements. This discretization presents a difficulty in obtaining the equivalent electric current for the NTFF. Straightforward finite element interpolation of the fields is shown to give poor accuracy. Alternative interpolation methods are recommended. An expanding magnetic quadrupole pulse benchmark problem, which is fully developed in the appendices, provides the basis for quantitative comparison.

FETD↗

DOE OSTI2021

Source term method for binary neutron stars initial data

The initial condition problem for a binary neutron star system requires a Poisson equation solver for the velocity potential with a Neumann-like boundary condition on the surface of the star. Difficulties that arise in this boundary value problem are: (a) the boundary is not known a priori, but constitutes part of the solution of the problem; (b) various terms become singular at the boundary. In this work, we present a new method to solve the fluid Poisson equation for irrotational/spinning binary neutron stars. The advantage of the new method is that it does not require complex fluid surface fitted coordinates and it can be implemented in a Cartesian grid, which is a standard choice in numerical relativity calculations. This is accomplished by employing the source term method proposed by Towers, where the boundary condition is treated as a jump condition and is incorporated as additional source terms in the Poisson equation, which is then solved iteratively. The issue of singular terms caused by vanishing density on the surface is resolved with an additional separation that shifts the computation boundary to the interior of the star. We present two-dimensional tests to show the convergence of the source term method, and we further apply this solver to a realistic three-dimensional binary neutron star problem. By comparing our solution with the one coming from the initial data solver cocal, we demonstrate agreement to approximately 1%. We report our method can be used in other problems with non-smooth solutions like in magnetized neutron stars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Determining the tilt of the Raman laser beam using an optical method for atom gravimeters

The tilt of a Raman laser beam is a major systematic error in precision gravity measurement using atom interferometry. The conventional approach to evaluating this tilt error involves modulating the direction of the Raman laser beam and conducting time-consuming gravity measurements to identify the error minimum. In this work, we demonstrate a method to expediently determine the tilt of the Raman laser beam by transforming the tilt angle measurement into characterization of parallelism, which integrates the optical method of aligning the laser direction, commonly used in freely falling corner-cube gravimeters, into an atom gravimeter. A position-sensing detector (PSD) is utilized to quantitatively characterize the parallelism between the test beam and the reference beam, thus measuring the tilt precisely and rapidly. After carefully positioning the PSD and calibrating the relationship between the distance measured by the PSD and the tilt angle measured by the tiltmeter, we achieved a statistical uncertainty of less than 30 µrad in the tilt measurement. Furthermore, we compared the results obtained through this optical method with those from the conventional tilt modulation method for gravity measurement. The comparison validates that our optical method can achieve tilt determination with an accuracy level of better than 200 µrad, corresponding to a systematic error of 20 µGal in g measurement. This work has practical implications for real-world applications of atom gravimeters.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2020

Methods for the robust computation of the long-period seismic spectrum of broad-band arrays

SUMMARY We describe array methods to search for low signal-to-noise ratio (SNR) signals in long-period seismic data using Fourier analysis. This is motivated by published results that find evidence of solar free oscillations in the Earth's seismic hum. Previous work used data from only one station. In this paper, we describe methods for computing spectra from array data. Arrays reduce noise level through averaging and provide redundancy that we use to distinguish coherent signal from a random background. We describe two algorithms for calculating a robust spectrum from seismic arrays, an algorithm that automatically removes impulsive transient signals from data, a jackknife method for estimating the variance of the spectrum, and a method for assessing the significance of an entire spectral band. We show examples of their application to data recorded by the Homestake Mine 3-D array in Lead, SD and the Piñon Flats PY array. These are two of the quietest small aperture arrays ever deployed in North America. The underground Homestake data has exceptionally low noise, and the borehole sensors of the PY array also have very low noise, making these arrays well suited to finding very weak signals. We find that our methods remove transient signals effectively from the data so that even low-SNR signals in the seismic background can be found and tested. Additionally, we find that the jackknife variance estimate is comparable to the noise floor, and we present initial evidence for solar g-modes in our data through the T2 test, a multivariate generalization of Student's t-test.

Caton, Ross C.↗

DOE OSTI2021

An efficient LC-MS method for isomer separation and detection of sugars, phosphorylated sugars, and organic acids

Abstract Assessing central carbon metabolism in plants can be challenging due to the dynamic range in pool sizes, with low levels of important phosphorylated sugars relative to more abundant sugars and organic acids. Here, we report a sensitive liquid chromatography–mass spectrometry method for analysing central metabolites on a hybrid column, where both anion-exchange and hydrophilic interaction chromatography (HILIC) ligands are embedded in the stationary phase. The liquid chromatography method was developed for enhanced selectivity of 27 central metabolites in a single run with sensitivity at femtomole levels observed for most phosphorylated sugars. The method resolved phosphorylated hexose, pentose, and triose isomers that are otherwise challenging. Compared with a standard HILIC approach, these metabolites had improved peak areas using our approach due to ion enhancement or low ion suppression in the biological sample matrix. The approach was applied to investigate metabolism in high lipid-producing tobacco leaves that exhibited increased levels of acetyl-CoA, a precursor for oil biosynthesis. The application of the method to isotopologue detection and quantification was considered through evaluating 13C-labeled seeds from Camelina sativa. The method provides a means to analyse intermediates more comprehensively in central metabolism of plant tissues.

Plant Sciences↗

DOE OSTI2022

The Davis–Chandrasekhar–Fermi method revisited

ABSTRACT Despite the rich observational results on interstellar magnetic fields in star-forming regions, it is still unclear how dynamically significant the magnetic fields are at varying physical scales, because direct measurement of the field strength is observationally difficult. The Davis–Chandrasekhar–Fermi (DCF) method has been the most commonly used method to estimate the magnetic field strength from polarization data. It is based on the assumption that gas turbulent motion is the driving source of field distortion via linear Alfvén waves. In this work, using MHD simulations of star-forming clouds, we test the validity of the assumption underlying the DCF method by examining its accuracy in the real 3D space. Our results suggest that the DCF relation between turbulent kinetic energy and magnetic energy fluctuation should be treated as a statistical result instead of a local property. We then develop and investigate several modifications to the original DCF method using synthetic observations, and propose new recipes to improve the accuracy of DCF-derived magnetic field strength. We further note that the biggest uncertainty in the DCF analysis may come from the linewidth measurement instead of the polarization observation, especially since the line-of-sight gas velocity can be used to estimate the gas volume density, another critical parameter in the DCF method.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Efficient method for estimating the time evolution of the proto-neutron star mass and radius from a supernova neutrino signal

In this paper, we present a novel method to estimate the time evolution of the proto-neutron star (PNS) structure from the neutrino signal in a core-collapse supernova (CCSN). Employing recent results from multidimensional CCSN simulations, we delve into a relation between the total emitted neutrino energy (TONE) and PNS mass/radius, and we find that they are strongly correlated with each other. We fit the relation by simple polynomial functions connecting the TONE to the mass and radius of the PNS as a function of time. By combining another fitting function representing the correlation between the TONE and the cumulative number of events at each neutrino observatory, the PNS mass and radius can be retrieved from purely observed neutrino data. We demonstrate retrievals of PNS mass and radius from mock data of the neutrino signal, and we assess the capability of our proposed method. While underlining the limitations of the method, we also discuss the importance of the joint analysis with the gravitational wave signal. This would reduce uncertainties of parameter estimations in our method, and may narrow down the possible neutrino oscillation model. The proposed method is a very easy and inexpensive computation, which will be useful in real data analysis of the CCSN neutrino signal.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Investigating event-shape methods in the search for the chiral magnetic effect in relativistic heavy ion collisions

The chiral magnetic effect (CME) is a phenomenon in which electric charge is separated by a strong magnetic field from local domains of chirality imbalance and parity violation in quantum chromodynamics. The CME-sensitive observable, the charge-dependent three-point azimuthal correlator Δ⁢𝛾 , is contaminated by a major physics background proportional to the particle's elliptic flow anisotropy 𝑣 2 . Event-shape engineering (ESE) binning events in dynamical fluctuations of 𝑣 2 and event-shape selection (ESS) binning events in statistical fluctuations of 𝑣 2 are two methods to search for the CME by projecting Δ⁢𝛾 to the measured anisotropy 𝑣 2 = 0 intercept. Here, we conduct a systematic study of these two methods using physics models as well as toy model simulations. It is observed that the ESE method fulfills the general premise of measuring the CME but is statistically hungry. It is found that the intercept from the ESS method depends on the details of the event content, such as the mixtures of background-contributing sources, because of statistical fluctuations of intertwining variables used in the method, and is thus not practically useful to measure the CME.

Relativistic heavy-ion collisions↗

DOE OSTI2021

Two-nucleon 𝑆-wave interactions at the SU(3) flavor-symmetric point with 𝑚 𝑢⁢𝑑 ≈ 𝑚$^{phys}_𝑠$: A first lattice QCD calculation with the stochastic Laplacian Heaviside method

We report on the first application of the stochastic Laplacian Heaviside method for computing multiparticle interactions with lattice QCD to the two-nucleon system. Like the Laplacian Heaviside method, this method allows for the construction of interpolating operators which can be used to construct a set of positive-definite two-nucleon correlation functions, unlike nearly all other applications of lattice QCD to two nucleons in the literature. It also allows for a variational analysis in which optimal linear combinations of the interpolating operators are formed that couple predominantly to the eigenstates of the system. Utilizing such methods has become of paramount importance to help resolve the discrepancy in the literature on whether two nucleons in either isospin channel form a bound state at pion masses heavier than physical, with the discrepancy persisting even in the SU(3)-flavor-symmetric point with all quark masses near the physical strange quark mass. This is the first in a series of papers aimed at resolving this discrepancy. In the present work, we employ the stochastic Laplacian Heaviside method without a hexaquark operator in the basis at a lattice spacing of 𝑎 ≈0.086 fm, lattice volume of 𝐿 = 48⁢𝑎 ≈ 4.1 fm and pion mass 𝑚 𝜋 ≈ 714 MeV. With this setup, the observed spectrum of two-nucleon energy levels strongly disfavors the presence of a bound state in either the deuteron or dineutron channel.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Applicability of semiclassical methods for modeling laser-enhanced fusion rates in a realistic setting

In the context of the potential laser-induced enhancement to the rates of DHe 3 and DT fusion, we discuss the frequently-used Wentzel-Kramers-Brillouin (WKB) method and the imaginary-time method (ITM). For static external electric fields, we find that these methods predict significant enhancement to the fusion cross section for electric-field strengths > 10 14 V/m, especially at low values (≈ keV) for the enter-of-mass (CoM) energy. When considering dynamic electric fields, this enhancement can be amplified by considering increased photon frequencies. However, we also pro- vide a review of the region of laser-parameter phase space where these semiclassical methods are applicable. Here, we conclude that this allowable region decreases for higher photon frequencies in con-junction with lower values for the electric-field strength, motivating the need for future experiments to test the predictions of these methods and their ranges of validity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Training and projecting: A reduced basis method emulator for many-body physics

Here, we present the reduced basis method as a tool for developing emulators for equations with tun able parameters within the context of the nuclear many-body problem. The method uses a basis expansion informed by a set of solutions for a few values of the model parameters and then projects the equations over a well-chosen low-dimensional subspace. We connect some of the results in the eigenvector continuation literature to the formalism of reduced basis methods and show how these methods can be applied to a broad set of problems. As we illustrate, the possible success of the formalism on such problems can be diagnosed beforehand by a principal component analysis. We apply the reduced basis method to the one-dimensional Gross-Pitaevskii equation with a harmonic trap ping potential and to nuclear density functional theory for 48 Ca, achieving speed-ups of more than x150 and x250, respectively, when compared to traditional solvers. The outstanding performance of the approach, together with its straightforward implementation, show promise for its application to the emulation of computationally demanding calculations, including uncertainty quantification.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Gauge-symmetrization method for energy-momentum tensors in high-order electromagnetic field theories

For electromagnetic field theories, canonical energy-momentum conservation laws can be derived from the underpinning spacetime translation symmetry according to the Noether procedure. However, the canonical energy-momentum tensors (EMTs) are neither symmetric nor gauge-symmetric (gauge invariant). The Belinfante-Rosenfeld (BR) method is a well-known procedure to symmetrize the EMTs, which also renders them gauge symmetric for first-order field theories. High-order electromagnetic field theories appear in the study of gyrokinetic systems for magnetized plasmas and the Podolsky system for the radiation reaction of classical charged particles. For these high-order field theories, gauge-symmetric EMTs are not necessarily symmetric and vice versa. In the present study, we develop a new gauge-symmetrization method for EMTs in high-order electromagnetic field theories. The Noether procedure is carried out using the Faraday tensor $F_{μν}$, instead of the 4-potential $A_{μ}$, to derive a canonical EMT $T^{μν}_{N}$. We show that the gauge-dependent part of $T^{μν}_{N}$ can be removed using the displacement-potential tensor $F^{σμν}$ ≡ $D^{σμ}A^{ν}/4π$, where $D^{σμ}$ is the antisymmetric electric displacement tensor. This method gauge-symmetrizes the EMT without necessarily making it symmetric, which is adequate for applications not involving general relativity. For first-order electromagnetic field theories, such as the standard Maxwell system, $F^{σμν}$ reduces to the familiar BR superpotential $S^{σμν}$, and the method developed can be used as a simpler procedure to calculate $S^{σμν}$ without employing the angular momentum tensor in 4D spacetime. When the electromagnetic system is coupled to classical charged particles, the gauge-symmetrization method for EMTs is shown to be effective as well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity

Here we present a realizability-preserving numerical method for solving a spectral two-moment model to simulate the transport of massless, neutral particles interacting with a steady background material moving with relativistic velocities. The model is obtained as the special relativistic limit of a four-momentum-conservative general relativistic two-moment model. Using a maximum-entropy closure, we solve for the Eulerian-frame energy and momentum. The proposed numerical method is designed to preserve moment realizability, which corresponds to moments defined by a nonnegative phase-space density. The realizability-preserving method is achieved with the following key components: (i) a discontinuous Galerkin phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit time-integration method; (iii) a realizability-preserving conserved to primitive moment solver; (iv) a realizability-preserving implicit collision solver; and (v) a realizability-enforcing limiter. Component (iii) is necessitated by the closure procedure, which closes higher order moments nonlinearly in terms of primitive moments. The nonlinear conserved to primitive and the implicit collision solves are formulated as fixed-point problems, which are solved with custom iterative solvers designed to preserve the realizability of each iterate. With a series of numerical tests, we demonstrate the accuracy and robustness of this discontinuous-Galerkin-implicit-explicit method.

79 ASTRONOMY AND ASTROPHYSICS↗