Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time integrator”

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 1,351 records · Page 75

Predicting Spacecraft Trajectories by the WeavEncke Method

A combination of methods is proposed of predicting spacecraft trajectories that possibly include multiple maneuvers and/or perturbing accelerations, with greater speed, accuracy, and repeatability than were heretofore achievable. The combination is denoted the WeavEncke method because it is based on unpublished studies by Jonathan Weaver of the orbit-prediction formulation of the noted astronomer Johann Franz Encke. Weaver evaluated a number of alternatives that arise within that formulation, arriving at an orbit-predicting algorithm optimized for complex trajectory operations. In the WeavEncke method, Encke's method of prediction of perturbed orbits is enhanced by application of modern numerical methods. Among these methods are efficient Kepler s-equation time-of-flight solutions and self-starting numerical integration with time as the independent variable. Self-starting numerical integration satisfies the requirements for accuracy, reproducibility, and efficiency (and, hence, speed). Self-starting numerical integration also supports fully analytic regulation of integration step sizes, thereby further increasing speed while maintaining accuracy.

Weaver, Jonathan K.↗

Enabling topography-resolving structural dynamic contact simulation

Damping of structures and systems is often dominated by frictional dissipation in connections, the prediction of which remains a longstanding scientific challenge. Previous studies have shown that the actual topography of contact interfaces may have a strong effect, especially in the partial slip/liftoff regime. We recently proposed a multi-scale method, which couples finite element and boundary element modeling. The primary benefit of this approach is that it permits to analyze the effect of the actual contact topography on the dynamics of jointed structures. While this multi-scale modeling method was initially developed for quasi-static analysis, we demonstrate herein how it can be used for time step integration and Harmonic Balance analysis. We cross-verify those fully dynamic analysis methods against each other and quasi-static results, for the S4 Beam benchmark. We compare the multi-scale method against state-of-the-art full-FE analysis, in terms of numerical damping and computational performance. Some discrepancy is found to be of physical origin. Depending on the load history, it is shown that the system settles to a slightly different equilibrium. Finally, transient multi-scale simulations enable the prediction of this interesting phenomenon, for the first time, for a structure with bolted joints.

Frictional-unilateral contact↗

Sparse Cholesky factorization for solving nonlinear PDEs via Gaussian processes

In recent years, there has been widespread adoption of machine learning-based approaches to automate the solving of partial differential equations (PDEs). Among these approaches, Gaussian processes (GPs) and kernel methods have garnered considerable interest due to their flexibility, robust theoretical guarantees, and close ties to traditional methods. They can transform the solving of general nonlinear PDEs into solving quadratic optimization problems with nonlinear, PDE-induced constraints. However, the complexity bottleneck lies in computing with dense kernel matrices obtained from pointwise evaluations of the covariance kernel, and its partial derivatives, a result of the PDE constraint and for which fast algorithms are scarce. The primary goal of this paper is to provide a near-linear complexity algorithm for working with such kernel matrices. We present a sparse Cholesky factorization algorithm for these matrices based on the near-sparsity of the Cholesky factor under a novel ordering of pointwise and derivative measurements. The near-sparsity is rigorously justified by directly connecting the factor to GP regression and exponential decay of basis functions in numerical homogenization. We then employ the Vecchia approximation of GPs, which is optimal in the Kullback-Leibler divergence, to compute the approximate factor. This enables us to compute ϵ-approximate inverse Cholesky factors of the kernel matrices with complexity O(N log d (N/ϵ)) in space and O(N log 2d (N/ϵ)) in time. We integrate sparse Cholesky factorizations into optimization algorithms to obtain fast solvers of the nonlinear PDE. We numerically illustrate our algorithm’s near-linear space/time complexity for a broad class of nonlinear PDEs such as the nonlinear elliptic, Burgers, and Monge-Ampère equations. In summary, we provide a fast, scalable, and accurate method for solving general PDEs with GPs and kernel methods.

97 MATHEMATICS AND COMPUTING↗

Electric propulsion system technology

The work performed on the Ion Propulsion System Technology Task in FY90 is described. The objectives of this work fall under two broad categories. The first of these deals with issues associated with the application of xenon ion thrusters for primary propulsion of planetary spacecraft, and the second with the investigation of technologies which will facilitate the development of larger, higher power ion thrusters to support more advanced mission applications. Most of the effort was devoted to investigation of the critical issues associated with the use of ion thrusters for planetary spacecraft. These issues may be succinctly referred to as life time, system integration, and throttling. Chief among these is the engine life time. If the engines do not have sufficient life to perform the missions of interest, then the other issues become unimportant. Ion engine life time was investigated through two experimental programs: an investigation into the reduction of ion engine internal sputter erosion through the addition of small quantities of nitrogen, and a long duration cathode life test. In addition, a literature review and analysis of accelerator grid erosion were performed. The nitrogen addition tests indicated that the addition of between 0.5 and 1.0 percent of nitrogen by mass to the xenon propellant results in a reduction in the sputter erosion of discharge chamber components by a factor of between 20 and 50, with negligible reduction in thruster performance. The long duration test of a 6.35-mm dia. xenon hollow cathode is still in progress, and has accumulated more than 4,000 hours of operation at an emission current of 25 A at the time of this writing. One of the major system integration issues concerns possible interactions of the ion thruster produced charge exchange plasma with the spacecraft. A computer model originally developed to describe the behavior of mercury ion thruster charge exchange plasmas was resurrected and modified for xenon propellant. This model enables one to calculate the flow direction and local density of the charge exchange plasma, and indicates the degree to which this plasma can flow upstream of the thruster exhaust plane. A continuing effort to investigate the most desirable throttling technique for noble gas ion thrusters concentrated this year on experimentally determining the fixed flow rate throttling range of a 30-cm dia. thruster with a two-grid accelerator system. These experiments demonstrated a throttling capability which covers a 2.8 to 1 variation in input power. This throttling range is 55 percent greater than expected, and is due to better accelerator system performance at low net-to-total voltage ratios than indicated in the literature. To facilitate the development of large, higher power ion thrusters several brief studies were performed. These include the development of a technique which simulates ion thruster operation without beam extraction, the development of an optical technique to measure ion thruster grid distortion due to thermal expansion, tests of a capacitance measurement technique to quantify the accelerator system grid separation, and the development of a segmented thruster geometry which enables near term development of ion thrusters at power levels greater than 100 kW. Finally, a paper detailing the benefits of electric propulsion for the Space Exploration Initiative was written.

Brophy, John R.↗

Quantum dynamics of cold atomic gas with SU(1,1) symmetry

Motivated by recent advances in quantum dynamics, we investigate the dynamics of the system with SU(1,1) symmetry. Instead of performing the time-ordered integral for the evolution operator of the time-dependent Hamiltonian, we show that the time evolution operator can be expressed as an SU(1,1) group element. Since the SU(1,1) group describes the “rotation” on a hyperbolic surface, the dynamics can be visualized on a Poincaré disk, a stereographic projection of the upper hyperboloid. As an example, we present the trajectory of the revival of Bose-Einstein condensation and that of the scale-invariant Fermi gas on the Poincaré disk. Further considering quantum gases in an oscillating lattice, we also study the dynamics of the system with time-dependent single-particle dispersion.

74 ATOMIC AND MOLECULAR PHYSICS↗

Integration of 5G and Time Sensitive Networks in Fossil Energy Generation Systems: A Case Study

Precise timing and data transmission within stringent time constraints are critical for numerous applications, such as robotics, virtual and augmented reality, industrial automation, energy and medical business, and various other sectors. Time-sensitive networks (TSN) and fifth-generation wireless communications (5G) are crucial for industrial communications, enabling convergent communication for various services using a common network core. Applications that are time-sensitive and require deterministic communications with low latency fall into this category, such as generation plant systems. This paper presents a simulated model of 5G-TSN for a fossil power plant based on the wired network parameters implemented in situ. Metrics analysis, comparison and future work are presented. © 2024 IEEE.

01 COAL, LIGNITE, AND PEAT↗

AI-driven neutrino diagnostics and radiation-hard beam instrumentation for next-generation neutrino experiments

The Long Baseline Neutrino Facility (LBNF) at Fermilab will deliver a high-intensity, multi-megawatt neutrino beam to the Deep Underground Neutrino Experiment (DUNE), enabling precision tests of the three-neutrino paradigm, CP violation searches, neutrino mass ordering determination, and supernova neutrino studies. To accelerate DUNE’s physics reach and ensure robust beam operations, we propose an integrated AI-driven framework with real-time diagnostics and radiation-hardened instrumentation. At its core is a Real-Time Beam Integrity Monitor using a physics-informed Digital Twin. By reconstructing pion phase space from muon profiles and exploiting magnetic horn optic linearity, it enables spill-by-spill beam correction and flux stabilization. By using this approach, flux-related systematics could be reduced from 5% to 1%, potentially accelerating the discovery of CP violations by four to six years. Complementing this, a US–Japan R&D effort will deploy a LAPPD-based muon monitor in the NuMI beamline. ToF measurements can be acquired with picosecond precision using this radiation-hard system, enhancing sensitivity to horn chromatic effects. Simulations confirm strong response to these effects. ML models predict beam quality and horn current to sub-percent accuracy from muon data, enhancing anomaly detection and stability. This scalable, AI-enabled strategy improves beam fidelity, reduces systematics, and sets a new standard for high-power accelerator operations.

Ganguly, Sudeshna [Fermilab] (ORCID:00000003163482↗

AI-driven neutrino diagnostics and radiation-hard beam instrumentation for next-generation neutrino experiments

The Long Baseline Neutrino Facility (LBNF) at Fermilab will deliver a high-intensity, multi-megawatt neutrino beam to the Deep Underground Neutrino Experiment (DUNE), enabling precision tests of the three-neutrino paradigm, CP violation searches, neutrino mass ordering determination, and supernova neutrino studies. To accelerate DUNE’s physics reach and ensure robust beam operations, we propose an integrated AI-driven framework with real-time diagnostics and radiation-hardened instrumentation. At its core is a Real-Time Beam Integrity Monitor using a physics-informed Digital Twin. By reconstructing pion phase space from muon profiles and exploiting magnetic horn optic linearity, it enables spill-by-spill beam correction and flux stabilization. By using this approach, flux-related systematics could be reduced from 5% to 1%, potentially accelerating the discovery of CP violations by four to six years. Complementing this, a US–Japan R&D effort will deploy a LAPPD-based muon monitor in the NuMI beamline. ToF measurements can be acquired with picosecond precision using this radiation-hard system, enhancing sensitivity to horn chromatic effects. Simulations confirm strong response to these effects. ML models predict beam quality and horn current to sub-percent accuracy from muon data, enhancing anomaly detection and stability. This scalable, AI-enabled strategy improves beam fidelity, reduces systematics, and sets a new standard for high-power accelerator operations.

Ganguly, Sudeshna [Fermilab] (ORCID:00000003163482↗

Analytical and experimental study of vibrations in a gear transmission

An analytical simulation of the dynamics of a gear transmission system is presented and compared to experimental results from a gear noise test rig at the NASA Lewis Research Center. The analytical procedure developed couples the dynamic behaviors of the rotor-bearing-gear system with the response of the gearbox structure. The modal synthesis method is used in solving the overall dynamics of the system. Locally each rotor-gear stage is modeled as an individual rotor-bearing system using the matrix transfer technique. The dynamics of each individual rotor are coupled with other rotor stages through the nonlinear gear mesh forces and with the gearbox structure through bearing support systems. The modal characteristics of the gearbox structure are evaluated using the finite element procedure. A variable time steping integration routine is used to calculate the overall time transient behavior of the system in modal coordinates. The global dynamic behavior of the system is expressed in a generalized coordinate system. Transient and steady state vibrations of the gearbox system are presented in the time and frequency domains. The vibration characteristics of a simple single mesh gear noise test rig is modeled. The numerical simulations are compared to experimental data measured under typical operating conditions. The comparison of system natural frequencies, peak vibration amplitudes, and gear mesh frequencies are generally in good agreement.

Choy, F. K.↗

Discrete Deterministic and Stochastic Petri Nets

Petri nets augmented with timing specifications gained a wide acceptance in the area of performance and reliability evaluation of complex systems exhibiting concurrency, synchronization, and conflicts. The state space of time-extended Petri nets is mapped onto its basic underlying stochastic process, which can be shown to be Markovian under the assumption of exponentially distributed firing times. The integration of exponentially and non-exponentially distributed timing is still one of the major problems for the analysis and was first attacked for continuous time Petri nets at the cost of structural or analytical restrictions. We propose a discrete deterministic and stochastic Petri net (DDSPN) formalism with no imposed structural or analytical restrictions where transitions can fire either in zero time or according to arbitrary firing times that can be represented as the time to absorption in a finite absorbing discrete time Markov chain (DTMC). Exponentially distributed firing times are then approximated arbitrarily well by geometric distributions. Deterministic firing times are a special case of the geometric distribution. The underlying stochastic process of a DDSPN is then also a DTMC, from which the transient and stationary solution can be obtained by standard techniques. A comprehensive algorithm and some state space reduction techniques for the analysis of DDSPNs are presented comprising the automatic detection of conflicts and confusions, which removes a major obstacle for the analysis of discrete time models.

MATHEMATICAL MODELS↗

Introducing and Evaluating the Climate Hazards Center IMERG with Stations (CHIMES) Timely Station-Enhanced Integrated Multisatellite Retrievals for Global Precipitation Measurement

As human exposure to hydroclimatic extremes increase and the number of in situ precipitation observations declines, precipitation estimates, such as those provided by the Integrated Multisatellite Retrievals for Global Precipitation Measurement (GPM) (IMERG) mission, provide a critical source of information. Here, we present a new gauge-enhanced dataset [the Climate Hazards Center IMERG with Stations (CHIMES)] designed to support global crop and hydrologic modeling and monitoring. CHIMES enhances the IMERG Late Run product using an updated Climate Hazards Center (CHC) high-resolution climatology (CHPclim) and low-latency rain gauge observations. CHPclim differs from other products because it incorporates long-term averages of satellite precipitation, which increases CHPclim ’s fidelity in data-sparse areas with complex terrain. This fidelity translates into performance increases in unbiased IMERGlate data, which we refer to as CHIME. This is augmented with gauge observations to produce CHIMES. The CHC’s curated rain gauge archive contains valuable contributions from many countries. There are two versions of CHIMES: preliminary and final. The final product has more copious and better-curated station data. Every pentad and month, bias-adjusted IMERGlate fields are combined with gauge observations to create pentadal and monthly CHIMESprelim and CHIMESfinal. Comparisons with pentadal, high-quality gridded station data show that IMERG late performs well (r = 0.75), but has some systematic biases which can be reduced. Monthly cross-validation results indicate that unbiasing increases the variance explained from 50% to 63% and decreases the mean absolute error from 48 to 39 mm month −1. Gauge enhancement then increases the variance explained to 75%, reducing the mean absolute error to 27 mm month −1.

Chris C funk↗

Mapping National Forest Aboveground Biomass in Mexico By Integrating GEDI and Landsat Times Series Data

Mexico is one of the countries with great potential for the UN's Reducing Emissions from Deforestation and Forest Degradation (REDD+) program, a key nature-based solution for the forest sector. To monitor carbon stock changes, there is a growing demand for unbiased Monitoring Reporting Verification (MRV) systems to facilitate effective forest management and climate change mitigation strategies. Remote sensing-based national aboveground biomass density (AGBD) estimation over Mexico is scarce and often limited to one-time static mapping, leading to spatiotemporal inconsistency in inputs. As an effort under NASA's Carbon Monitoring System (CMS) program, we have developed a remote sensing-based approach to create consistent historical AGBD maps of Mexico using multi-stream remote sensing data, including spaceborne lidar GEDI and long-term Landsat time series, as well as topographic information. We employ the continuous change detection and classification (CCDC) algorithm for temporal modeling of Landsat surface reflectance, followed by the inference of forest AGBD using a random forest machine learning algorithm with the temporal information of land surface dynamics extracted by the CCDC as input. GEDI provides unprecedented forest structure and AGBD sampling datasets for model training and validation practices. In this presentation, we share the progress made in developing a spatially explicit mapping of historical AGBD changes associated with land surface changes and post-disturbance landscapes.

Taejin Park↗

Integration of real-time measurements and atomistic modeling to license nuclear components

An improved testing and data gathering method is described herein with reference to testing a new fuel, as an exemplary component to be tested for licensure purposes. The method includes generally: generating models for the new system; making samples and getting them accepted for a reactor; and testing the samples in a test or commercial reactor until the exposure time is reached for the expected cycle length of the fuel at final use. The method is preferably done concurrent to submitting a license application for commercial use of the new component to the relevant government body to expedite license testing review.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Flip-flop resolving time test circuit

Integrated circuit (IC) flip-flop resolving time parameters are measured by wafer probing, without need of dicing or bonding, throught the incorporation of test structures on an IC together with the flip-flop to be measured. Several delays that are fabricated as part of the test circuit, including a voltage-controlled delay with a resolution of a few picosecs, are calibrated as part of the test procedure by integrating them into, and out of, the delay path of a ring oscillator. Each of the delay values is calculated by subtracting the period of the ring oscillator with the delay omitted from the period with the delay included. The delay measurement technique is sufficiently general for other applications. The technique is illustrated for the case of the flip-flop parameters of a 5-micron feature size NMOS circuit.

Rosenberger, F.↗