SEARCH · Engineering Papers

Results for “Numerical techniques”

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 1,243 records · Page 69

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Improving the numerical integration solution of satellite orbits in the presence of solar radiation pressure using modified back differences

The method of modified back differences, a technique that significantly reduces the numerical integration errors associated with crossing shadow boundaries using a fixed-mesh multistep integrator without a significant increase in computer run time, is presented. While Hubbard's integral approach can produce significant improvements to the trajectory solution, the interpolation method provides the best overall results. It is demonstrated that iterating on the point mass term correction is also important for achieving the best overall results. It is also shown that the method of modified back differences can be implemented with only a small increase in execution time.

Lundberg, J. B.↗

DOE OSTI2020

Understanding Differences in California Climate Projections Produced by Dynamical and Statistical Downscaling

Abstract We compare historical and end‐of‐century temperature and precipitation patterns over California from one dynamically downscaled simulation using the Weather Research and Forecast (WRF) model and two simulations statistically downscaled using Localized Constructed Analogs (LOCA). We uniquely separate causes of differences between dynamically and statistically based future climate projections into differences in historical climate (gridded observations versus regional climate model output) and differences in how these downscaling techniques explicitly handle future climate changes (numerical modeling versus analogs). In these methods, solutions between different downscaling techniques differ more in the future compared to the historical period. Changes projected by LOCA are insensitive to the choice of driving data. Only through dynamical downscaling can we simulate physically consistent regional springtime warming patterns across the Sierra Nevada, while the statistical simulations inherit an unphysical signal from their parent Global Climate Model (GCM) or gridded data. The results of our study clarify why these different techniques produce different outcomes and may also provide guidance on which downscaled products to use for certain impact analyses in California and perhaps other Mediterranean regimes.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1977

Probe corrected far-field reconstruction from measurements on a cylinder: A novel formulation and efficient algorithm

A novel and numerically efficient method of far field evaluation from measurements taken on a cylinder is based on the representation of both the antenna and the probe fields as superpositions of plane waves. A system of two integral equations are established whose unknown functions are the azimuthal and elevation components of the antenna pattern and whose known terms are the set of measurement data taken with two different probes - the second probe in most practical instances being simply the same probe with a different geometrical orientation. The equations express the known data - for each angular position of the antenna under measurement - as the integrals of the products of the corresponding components of the unknown antenna and known probe patterns multiplied by a phase term. The convolutional nature of the integral equations makes their solutions straight-forward. If, as is virtually always the case, the probe is small or of moderate size so that the axis of rotation of the antenna mount is in the far field of the probe, the intervention of asymptotic techniques makes the solution numerically very efficient. The agreement of calculated and experimental patterns is excellent.

Borgiotti, G. V.↗

NASA NTRS1978

Research on an augmented Lagrangian penalty function algorithm for nonlinear programming

The augmented Lagrangian (ALAG) Penalty Function Algorithm for optimizing nonlinear mathematical models is discussed. The mathematical models of interest are deterministic in nature and finite dimensional optimization is assumed. A detailed review of penalty function techniques in general and the ALAG technique in particular is presented. Numerical experiments are conducted utilizing a number of nonlinear optimization problems to identify an efficient ALAG Penalty Function Technique for computer implementation.

Frair, L.↗

NASA NTRS1973

Eigenproblem solution by a combined Sturm sequence and inverse iteration technique.

Description of an efficient and numerically stable algorithm, along with a complete listing of the associated computer program, developed for the accurate computation of specified roots and associated vectors of the eigenvalue problem Aq = lambda Bq with band symmetric A and B, B being also positive-definite. The desired roots are first isolated by the Sturm sequence procedure; then a special variant of the inverse iteration technique is applied for the individual determination of each root along with its vector. The algorithm fully exploits the banded form of relevant matrices, and the associated program written in FORTRAN V for the JPL UNIVAC 1108 computer proves to be most significantly economical in comparison to similar existing procedures. The program may be conveniently utilized for the efficient solution of practical engineering problems, involving free vibration and buckling analysis of structures. Results of such analyses are presented for representative structures.

Gupta, K. K.↗

DOE OSTI2020

The use of seismic spatial gradients in a single layer neural network for seismic source discrimination: proof of concept

This report describes a proof-of-concept method of seismic source discrimination using seismic gradiometry and a common machine learning technique. The tests described here are purely numerical, using synthetic seismic data and well understood mathematical techniques. The primary innovation described here is the application of a richer seismic data set derived from seismic gradiometry. Seismic gradiometry is a method to estimate the time variable spatial gradient of the wavefield to compute various wavefield attributes such as slowness, dynamic strain, and rotational motions. With the addition of these wavefield attributes, we are afforded up to twenty "compo- nents" of time series data measured at a single point on, or in, the Earth. This is in direct contrast to conventional three-component seismic data collected at several locations using a seismic network. Using the gradiometrically-derived wavefield components directly in a single-layer neural network, I show that it is possible to discriminate between three common seismic source types (earthquakes, explosions, and opening fractures) for various noise conditions and gradiometry configurations.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

NASA NTRS1984

The effect of rotation on turbulence structure

Rotating turbulent channel flow has been simulated numerically using the large-eddy-simulation technique. With a three-dimensional, time-dependent, numerical simulation of turbulence, it was possible to obtain detailed flow-field information which had not been available previously from experiments. This information was used to study the structure and statistical properties of the flow with externally imposed rotation. It was found that in agreement with experimental observations, the Coriolis force introduced by the system rotation could stabilize/destabilize turbulence. The large-scale motion observed in the experiment was also found to exist in the computed flow field.

Kim, J.↗

NASA NTRS1979

A computational fluid dynamic technique valid at the centerline for non-axisymmetric problems in cylindrical coordinates

A technique is described for the numerical solution of non-axisymmetric flow problems posed in cylindrical coordinates when the z-axis is included in the flowfield. The highlight of the technique is the manner in which the singularities at the centerline are handled. Specifically, the governing flowfield equations at r = 0 are put in a special form by applying L'Hospital's Rule. The required radial derivatives are evaluated using a one-sided, second-order accurate, first-difference. This leads to a smooth, convergent calculation of the flowfield at the centerline. This appears to be the first generally applicable numerical method for avoiding coordinate system singularities in the context of a finite-difference scheme, and could have application to many nonaxisymmetric flows. The technique is illustrated by specific results for the time-dependent flowfield inside an internal combustion engine.

Griffin, M. D.↗

DOE OSTI2025

Non-propagating structures and propagating waves in solar wind turbulence revealed by simulations and observations

Structures and waves are common features of solar wind turbulence at various scales. The interplay between structures and waves is important for processes such as the turbulent energy cascade, plasma heating, and particle scattering. Our understanding of turbulence has been advanced by not only new space missions and numerical simulations, but also techniques that have been developed to interpret the rapidly growing turbulence data. We review basic models of turbulence with a specific focus on the analysis methods for understanding magnetic structures and waves. MHD and kinetic waves in single-spacecraft time series measurements can be identified through mode decomposition or their characteristic polarization signatures. The structures in this paper are considered as zero-frequency, non-propagating or convected modes embedded in the solar wind. The synergy between observations and simulations is most evident in the application of spatial-temporal analysis to multi-spacecraft observation and turbulence simulations. The spatial-temporal analysis has greatly improved our understanding of structures and waves in turbulence. We conclude by discussing prospects for future research.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

CG-Kit: Code Generation Toolkit for performant and maintainable variants of source code applied to Flash-X hydrodynamics simulations

CG-Kit is a new Code Generation tool-Kit that we have developed as a part of the solution for portability and maintainability for multiphysics computing applications. The development of CG-Kit is rooted in the urgent need created by the shifting landscape of high-performance computing platforms and the algorithmic complexities of a particular large-scale multiphysics application: Flash-X. To efficiently use computing resources on a heterogeneous node, an application must have a map of computation to resources and a mechanism to move the data and computation to the resources according to the map. Most existing performance portability solutions are focussed on abstracting the expression of computations so that a unified source code can be specialized to run on different resources. However, such an approach is insufficient for a code like Flash-X, which has a multitude of code components that can be assembled in various permutations and combinations to form different instances of applications. Similar challenges apply to any code that has composability, where a single specified way of apportioning work among devices may not be optimal. Additionally, use cases arise where the optimal control flow of computation may differ for different devices while the underlying numerics remain identical. This combination leads to unique challenges including handling an existing large code base in Fortran and/or C/C++, subdivision of code into a great variety of units supporting a wide range of physics and numerical methods, different parallelization techniques for distributed and shared memory systems and accelerator devices, and heterogeneity of computing platforms requiring coexisting variants of parallel algorithms. All of these challenges demand that scientific software developers apply existing knowledge about domain applications, algorithms, and computing platforms to determine custom abstractions and granularity for code generation. There is a critical lack of tools to tackle those problems. CG-Kit is designed to fill this gap by providing a user with the ability to express their desired control flow and computation-to-resource map in the form a pseudocode-like recipe. It consists of standalone tools that can be combined into highly specific and, we argue, highly effective portability and maintainability toolchains. Here we present the design of our new tools: parametrized source trees, control flow graphs, and recipes. The tools are implemented in Python. They are agnostic to the programming language of the source code targeted for code generation. In conclusion, we demonstrate the capabilities of the toolkit with two examples, first, multithreaded variants of the basic AXPY operation, and second, variants of parallel algorithms within a hydrodynamics solver, called Spark, from Flash-X that operates on block-structured adaptive meshes.

Algorithmic portability↗

DOE OSTI2020

Elevated temperature mechanical properties of Inconel 617 surface oxide using nanoindentation

Inconel 617 is a principal candidate material for helium gas cooled very-high-temperature reactors with outlet temperatures of 700–950 °C. Recent findings showed that this alloy develops unique surface oxide layers especially at high temperature (HT) helium environment with distinctive wear, friction and contact properties. This study investigates the elevated temperature mechanical properties of Inconel 617 top surface layers aged in HT helium environment. Nanoindentation technique is used to obtain load-displacement graphs of the alloy top surface oxide in temperatures ranging from 25 °C up to 600 °C. In addition, using finite element analysis along with an iterative regression technique, a semi-numerical method is developed to further measure and quantify the material parameters and, in particular, time-independent and creep characteristics of the oxide. Although Young's modulus of the oxide is found to be relatively close to the bulk for the tested temperatures, the yield strength and hardness, in comparison to the bulk material, increase significantly as the material is oxidized after aging. The oxide exhibits significant softening as the temperature increases to 600 °C. Unlike the bulk material, diffusion through the grains is found to be the dominant creep mechanism for the oxide. Considerable difference between the mechanical properties of the oxide and the bulk material shows the need for accurate measurements of near surface mechanical properties, if reliable predictive contact and tribological models are sought at elevated temperatures.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Understanding the enhancement of scintillation light in xenon-doped liquid argon

Measuring the scintillation light in noble gases is an important detection technique in particle physics. Numerous rare event searches like neutrino beam experiments, neutrino-less double beta-decay, and dark matter searches use argon-based detectors. In liquid argon, the light yield can be enhanced by the addition of a small quantity of xenon, where ~10–1000 ppm are added. The general enhancement mechanism and its pathway via an energy transfer between argon and xenon excimers is well known, however the importance of absorption of argon excimer emission by atomic xenon has not been fully appreciated. This absorption significantly reduces the light yield in commercially available argon (extracted from air) which contains trace amounts (~0.1 ppm) of xenon. The addition of a small xenon dopant of ~10 ppm recovers this lost light resulting in an increased light yield over un-doped argon of about a factor of two. Here, in this paper, we introduce a model for the light production in xenon doped argon, including absorption and re-emission, and compare it to the measured time dependence of light emission in xenon-doped argon.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2021

Moiré metrology of energy landscapes in van der Waals heterostructures

The emerging field of twistronics, which harnesses the twist angle between two-dimensional materials, represents a promising route for the design of quantum materials, as the twist-angle-induced superlattices offer means to control topology and strong correlations. At the small twist limit, and particularly under strain, as atomic relaxation prevails, the emergent moiré superlattice encodes elusive insights into the local interlayer interaction. Here we introduce moiré metrology as a combined experiment-theory framework to probe the stacking energy landscape of bilayer structures at the 0.1 meV/atom scale, outperforming the gold-standard of quantum chemistry. Through studying the shapes of moiré domains with numerous nano-imaging techniques, and correlating with multi-scale modelling, we assess and refine first-principle models for the interlayer interaction. We document the prowess of moiré metrology for three representative twisted systems: bilayer graphene, double bilayer graphene and H-stacked MoSe 2 /WSe 2 . Moiré metrology establishes sought after experimental benchmarks for interlayer interaction, thus enabling accurate modelling of twisted multilayers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Coulomb-free 1S0 p − p scattering length from the quasi-free p + d → p + p + n reaction and its relation to universality

Abstract The Coulomb-free 1 S 0 proton-proton ( p - p ) scattering length relies heavily on numerous and distinct theoretical techniques to remove the Coulomb contribution. Here, it has been determined from the half-off-the-energy-shell p - p scattering cross section measured at center-of-mass energies below 1 MeV using the quasi-free p + d → p + p + n reaction. A Bayesian data-fitting approach using the expression of the s-wave nucleon-nucleon scattering cross section returned a p - p scattering length $${a}_{pp}=-18.1{7}_{-0.58}^{+0.52}{| }_{stat}\pm 0.0{1}_{syst}$$ a p p = − 18.1 7 − 0.58 + 0.52 ∣ s t a t ± 0.0 1 s y s t fm and effective range r 0 = 2.80 ± 0.05 s t a t ± 0.001 s y s t fm. A model based on universality concepts has been developed to interpret this result. It accounts for the short-range interaction as a whole, nuclear and residual electromagnetic, according to what the s-wave phase-shift δ does in the description of low-energy nucleon-nucleon scattering data. We conclude that our parameters are representative of the short-range physics and propose to assess the charge symmetry breaking of the short-range interaction instead of just the nuclear interaction. This is consistent with the current understanding that the charge dependence of nuclear forces is due to different masses of up-down quarks and their electromagnetic interactions. This achievement suggests that these properties have a lesser than expected impact in the context of the charge symmetry breaking.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Dynamics of qudit gates and effects of spectator modes on optimal control pulses

Qudit gates for high-dimensional quantum computing can be synthesized with high precision using numerical quantum optimal control techniques. Large circuits are broken down into modules and the tailored pulses for each module can be used as primitives for a qudit compiler. Application of the pulses of each module in the presence of extra modes may decrease their effectiveness due to crosstalk. In this paper, we address this problem by simulating qudit dynamics for circuit quantum electrodynamics systems. As a test case, we take pulses for single-qudit swap gates optimized in isolation and then apply them in the presence of spectator modes each of which are in Fock states. We provide an experimentally relevant scaling formula that can be used as a bound on the fidelity decay. Our results show that frequency shift from spectator mode populations has to be ≲ 0.1 % of the qudit's nonlinearity in order for high-fidelity single-qudit gates to be useful in the presence of occupied spectator modes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Electron- and light-induced stimulated Raman spectroscopy for nanoscale molecular mapping

We propose and theoretically analyze a new vibrational spectroscopy, termed electron- and light-induced stimulated Raman (ELISR) scattering, that combines the high spatial resolution of elec-tron microscopy with the molecular sensitivity of surface-enhanced Raman spectroscopy. WithELISR, electron-beam excitation of plasmonic nanoparticles is utilized as a spectrally-broadband butspatially-confined Stokes beam in the presence of a diffraction-limited pump laser. To characterizethis technique, we develop a numerical model and conduct full-field electromagnetic simulations toinvestigate two distinct nanoparticle geometries, nanorods and nanospheres, coated with a Raman-active material. Our results show the significant (10 6 -10 7 ) stimulated Raman enhancement that isachieved with dual electron and optical excitation of thesenanoparticle geometries. Importantly,the spatial resolution of this vibrational spectroscopy for electron microscopy is solely determinedby the nanoparticle geometry and the plasmon mode volume. Our results highlight the promiseof ELISR for simultaneous high-resolution electron microscopy with sub-diffraction-limited Ramanspectroscopy, complementing advances in superresolutionmicroscopy, correlated light and electronmicroscopy, and vibrational electron energy loss spectroscopy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Visual Analytics of Crosstalk in Quantum Hardware

Crosstalk remains a major obstacle to building scalable and fault-tolerant quantum computers. Conventional diagnostic techniques-often based on numerical simulation or statistical modeling-struggle to scale with hardware complexity and offer limited interpretability. In this work, we present a visual analytics framework for diagnosing qubit crosstalk using lightweight, circuit-based models integrated with an interactive user interface. Our approach quantifies correlations between active and idle qubits under parameterized single- and twoqubit operations, enabling detection of both spatial and gateinduced crosstalk. The system incorporates qubit topology and gate performance data to support sector-based exploration and correlation mapping. This tool assists users in identifying correlated error sources, informing qubit placement strategies, and guiding noise-aware circuit design.

Chae, Junghoon [ORNL] (ORCID:0000000206016746)↗

DOE OSTI2026

The role of equatorial waves in triggering precipitation extremes in the Maritime Continent

This review offers a comprehensive analysis of convectively coupled equatorial waves (CCEWs) and their pivotal role in driving precipitation extremes across the Maritime Continent. It examines the current understanding of CCEWs, evaluates the performance of numerical models and forecasting techniques in predicting these phenomena, and pinpoints critical areas for improvement. The discussion centers on three key types of equatorial waves: equatorial Rossby waves, Kelvin waves, and mixed Rossby–gravity waves. By connecting scientific insights with practical forecasting applications, the review sheds light on the challenges of predicting these waves while identifying opportunities to advance both fundamental knowledge and forecasting accuracy. Designed as an educational resource, it targets operational forecasting centers, meteorologists, and researchers, aiming to enhance the prediction of extreme weather events in the region.

Maritime Continent↗