Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical integration”

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

A Novel 'Smart Microchip Proppants' Technology for Precision Diagnostics of Hydraulic Fracture Networks (Edited Final Report)

This project introduces innovative technology to improve subsurface characterization, visualization, and diagnostics of unconventional reservoirs (fossil resources). Through a collaborative effort involving the University of Kansas, UCLA, MicroSilicon Inc., and EOG Resources, the project aims to deliver precision diagnostics for hydraulic fractures using novel high-resolution imaging technology based on smart microchip proppants. Additionally, it seeks to enhance the accuracy and predictability of integrated numerical, and machine-learning modeling techniques for hydraulic fracture characterization and simulation. This groundbreaking technology addresses significant gaps in understanding unconventional and tight reservoir behavior and optimizing well-completion strategies, enabling more cost-efficient recovery of unconventional resources.

02 PETROLEUM↗
DOE OSTI2025

Exact Timestep for a Pairwise Coulomb Collision

Standard numerical integrators work well for many-body Coulomb repulsion problems when the timestep is much shorter than the timescale of relative position changes. However, for ‘hard’ collisions in which two particles have a near miss and exchange a lot of momentum within one timestep, they understandably struggle. This note proposes using the exact solution of Keplerian two-body orbits (usually hyperbolic) to calculate the momentum exchange with other particles: either a selection of the ‘closest’ ones or all of them.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

A Quantum Approach for Implementing Fixed-Point Arithmetic in Solving Ordinary Differential Equations

Differential equations (DEs) serve as fundamental tools in mathematical modeling across scientific disciplines, yet classical numerical solvers face limitations with large-scale or computationally intensive problems. This study explores a quantum-inspired approach to solving DEs, combining quantum- inspired techniques with classical methods. It focuses on fixed- point arithmetic on quantum circuits, utilizing basic quantum gates to manipulate DE solutions. We expand upon the techniques introduced by Zanger et al. [Quantum, 5, 502 (2021)] by offering a precise computation for a fixed-point signed multiplication scheme, while also presenting a quantum circuit capable of executing the fixed-point division algorithm. We demonstrate the feasibility of our approach through the simulation of a linear Ordinary Differential Equation (ODE), where initial conditions and parameters are encoded into quantum circuits using fixed- point representation. By executing sequences of quantum gates mimicking numerical integration steps, we obtain approximate solutions to the ODE with specified fixed-point precision.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Performance Improvements of the Griffin Solvers in FY24

The Griffin code is a MOOSE-based reactor physics application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the Department of Energy Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation Program. This fiscal year, we have made significant efforts to improve the performance of transport solver options and cross-section generation for the efficient use of Griffin in advanced reactor applications. For the HFEM-PN solver, the residual evaluations of HFEM kernels were optimized by utilizing the pre- computed averaged cross sections for individual elements. Numerical integration involving the evaluation of basis functions at quadrature points was bypassed by facilitating precomputed element mass matrices for response matrices. Red-black iterations were improved by introducing a new generalized minimum residual based solver. The memory usage of response matrix storage was significantly reduced by applying basis function rotations on interfaces and calculating volumetric odd-parity moments on the fly. Additionally, the adjoint flux and transient calculation capabilities of the HFEM-PN solver were successfully implemented and verified using the TWIGL benchmark problem. For the DFEM-SN solver, memory footprint and computation time were significantly reduced by not treating angular flux vectors as the MOOSE nonlinear system vectors. Specifically for IQS, scalar adjoint weighting was introduced to further eliminate angular adjoint flux storage in the MOOSE auxiliary system. It was demonstrated through the three-dimensional Advanced Burner Test Reactor core problem that the memory usage for transient calculations with the IQS method was reduced by over 7.5× compared to before the optimizations. For the self-shielding application programming interface, a new double-heterogeneity treatment method, named the Bell Function-Based Analytic Two-Region Slowing Down Method, was developed to efficiently flux-volume homogenize TRISO particles with the matrix. Additionally, optimizations were made to hyper- fine group (HFG) slowing down calculations by pretabulating collision probability coefficients and grouping isotopes, significantly reducing the computational time for calculating scattering sources per HFG. Lastly, the pin power reconstruction module was extended to account for temporal behavior in a microreactor analysis problem, specifically for a control drum transient. Verification tests for each of these improvements demonstrated significant performance enhancements and memory reduction.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

On the Efficient Evaluation of the Exchange Correlation Potential on Graphics Processing Unit Clusters

The predominance of Kohn–Sham density functional theory (KS-DFT) for the theoretical treatment of large experimentally relevant systems in molecular chemistry and materials science relies primarily on the existence of efficient software implementations which are capable of leveraging the latest advances in modern high-performance computing (HPC). With recent trends in HPC leading toward increasing reliance on heterogeneous accelerator-based architectures such as graphics processing units (GPU), existing code bases must embrace these architectural advances to maintain the high levels of performance that have come to be expected for these methods. In this work, we purpose a three-level parallelism scheme for the distributed numerical integration of the exchange-correlation (XC) potential in the Gaussian basis set discretization of the Kohn–Sham equations on large computing clusters consisting of multiple GPUs per compute node. In addition, we purpose and demonstrate the efficacy of the use of batched kernels, including batched level-3 BLAS operations, in achieving high levels of performance on the GPU. We demonstrate the performance and scalability of the implementation of the purposed method in the NWChemEx software package by comparing to the existing scalable CPU XC integration in NWChem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Potential Themis-family Asteroid Contribution to the Jupiter-family Comet Population

Recent dynamical analyses suggest that some Jupiter family comets (JFCs) may originate in the main asteroid belt instead of the outer solar system. This possibility is particularly interesting given evidence that icy main-belt objects are known to be present in the Themis asteroid family. We report results from dynamical analyses specifically investigating the possibility that icy Themis family members could contribute to the observed population of JFCs. Numerical integrations show that such dynamical evolution is indeed possible via a combination of eccentricity excitation apparently driven by the nearby 2:1 mean-motion resonance with Jupiter, gravitational interactions with planets other than Jupiter, and the Yarkovsky effect. We estimate that, at any given time, there may be tens of objects from the Themis family on JFC-like orbits with the potential to mimic active JFCs from the outer solar system, although not all, or even any, may necessarily be observably active. We find that dynamically evolved Themis family objects on JFC-like orbits have semimajor axes between 3.15 and 3.40 au for the vast majority of their time on such orbits, consistent with the strong role that the 2:1 mean-motion resonance with Jupiter likely plays in their dynamical evolution. We conclude that a contribution from the Themis family to the active JFC population is plausible, although further work is needed to better characterize this contribution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Meteoroid Stream of Comet C/1961 T1 (Seki) and Its Relation to the December ρ-Virginids and γ-Sagittariids

We answer the question of whether long-period comet C/1961 T1 (Seki) is the parent body of a meteor shower observed in Earth’s atmosphere. We model discrete parts of a meteoroid stream associated with this comet and follow the orbital evolution of its particles via a numerical integration of their orbits. The gravitational perturbations of major planets and the Poynting–Robertson effect, which influences the motion of the particles, are taken into account. We predict two showers which are related to C/1961 T1. One of these showers was identified with the December ρ-Virginids, #502. This identification confirms an earlier suggested relationship between the shower and the comet. The second predicted shower was identified with the γ-Sagittariids, #657. The radiant areas of both showers are symmetrical in respect to the apex of Earth’s heliocentric motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The Turbulent Stress Spectrum in the Inertial and Subinertial Ranges

For velocity and magnetic fields, the turbulent pressure and, more generally, the squared fields such as the components of the turbulent stress tensor, play important roles in astrophysics. For both one and three dimensions, we derive the equations relating the energy spectra of the fields to the spectra of their squares. We solve the resulting integrals numerically and show that for turbulent energy spectra of Kolmogorov type, the spectral slope of the stress spectrum is also of Kolmogorov type. For shallower turbulence spectra, the slope of the stress spectrum quickly approaches that of white noise, regardless of how blue the spectrum of the field is. For fully helical fields, the stress spectrum is elevated by about a factor of two in the subinertial range, while that in the inertial range remains unchanged. We discuss possible implications for understanding the spectrum of primordial gravitational waves from causally generated magnetic fields during cosmological phase transitions in the early universe. We also discuss potential diagnostic applications to the interstellar medium, where polarization and scintillation measurements characterize the square of the magnetic field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Phase-space Energization of Ions in Oblique Shocks

Abstract Examining energization of kinetic plasmas in phase space is a growing topic of interest, owing to the wealth of data in phase space compared to traditional bulk energization diagnostics. Via the field-particle correlation (FPC) technique and using multiple means of numerically integrating the plasma kinetic equation, we have studied the energization of ions in phase space within oblique collisionless shocks. The perspective afforded to us with this analysis in phase space allows us to characterize distinct populations of energized ions. In particular, we focus on ions that reflect multiple times off the shock front through shock-drift acceleration, and how to distinguish these different reflected populations in phase space using the FPC technique. We further extend our analysis to simulations of three-dimensional shocks undergoing more complicated dynamics, such as shock ripple, to demonstrate the ability to recover the phase-space signatures of this energization process in a more general system. This work thus extends previous applications of the FPC technique to more realistic collisionless shock environments, providing stronger evidence of the technique’s utility for simulation, laboratory, and spacecraft analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1962

Analytical Study of the Tumbling Motions of Vehicles Entering Planetary Atmospheres

The tumbling motion of vehicles entering planetary atmospheres is analyzed. A differential equation governing the tumbling motion, its arrest, and the subsequent oscillatory motion is obtained and identified as the equation for the fifth Painleve transcendant. An approximate analytical solution for the transcendant is derived. Comparisons with results obtained from numerical integration of the exact equations of motion indicate that the solution for the angle-of-attack history is sufficiently accurate to be of practical use.

PLANETARY ENTRY↗
NASA NTRS1968

Computer program determines exact two-sided tolerance limits for normal distributions

Computer program determines by numerical integration the exact statistical two-sided tolerance limits, when the proportion between the limits is at least a specified number. The program is limited to situations in which the underlying probability distribution for the population sampled is the normal distribution with unknown mean and variance.

Friedman, H. A.↗