Search NASA⌕ Search

SEARCH · Search NASA

Results for “Time 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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Optimal trajectories for aeroassisted, coplanar orbital transfer

Classical and minimax optimal control problems arising in the study of aeroassisted coplanar orbit transfer from a high planetary orbit to a low one are considered. Attention is given to (1) the minimization of the energy required for the maneuver; (2) minimization of the time integral of the heating rate; (3) minimization of the time of flight during the atmospheric portion of the trajectory; (4) maximization of the time of flight during the atmospheric portion of the trajectory; (5) minimization of the time integral of the path inclination; and (6) minimization of the sum of the squares of the entry and exit path inclinations.

Miele, A.↗
NASA NTRS1993

Confidence of the three-point estimator of frequency drift

It was shown that a three-point second difference estimator is nearly optimal for estimating frequency drift in many common atomic oscillators. A formula for the uncertainty of this estimate as a function of the integration time and of the Allan variance associated with this integration time is derived.

Weiss, Marc A.↗
NASA NTRS1990

Discrete integration of continuous Kalman filtering equations for time invariant second-order structural systems

A general form for the first-order representation of the continuous second-order linear structural-dynamics equations is introduced to derive a corresponding form of first-order continuous Kalman filtering equations. Time integration of the resulting equations is carried out via a set of linear multistep integration formulas. It is shown that a judicious combined selection of computational paths and the undetermined matrices introduced in the general form of the first-order linear structural systems leads to a class of second-order discrete Kalman filtering equations involving only symmetric sparse N x N solution matrices.

Park, K. C.↗
DOE OSTI2020

Towards improved speed and accuracy of laser powder bed fusion simulations via representation of multiple time scales

Here, we report the result of investigations into alternative representations of the temporal evolution of the laser powder bed fusion (LPBF) process. In order to provide the reader with sufficient context for the following discussion, the modeling challenge and its underlying cause, as well as current attempts to provide a satisfactory solution, will be discussed. Next, the author will describe two numerical methods (multirate time integration and parallel-in-time) to better represent the temporal scales of the problem. Finally, the results of applying these methods to a two-dimensional version of the heat transfer portion of the LPBF problem will be presented, reporting on both convergence and performance behavior, which indicate the possibility of significant speedup of solution calculations (40–100 times), along with well-characterized effects on solution accuracy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Parallel-in-Time Methods for Method-of-Lines Discretizations of Nonlinear Hyperbolic PDEs and Systems (Final Report)

The work for the subcontract is situated in the area of parallel-in-time integration for hyperbolic partial differential equations (PDEs). Parallel-in-time integration is an active area of research due to its ability to enable faster numerical simulations for applications throughout many areas of science. The work in this subcontract builds on a variety of results that were obtained, as part of the work performed for Subcontract No. B648355, for the Multigrid Reduction-in-Time (MGRIT) method from [1] applied to hyperbolic PDEs. This subcontract extends these results further to more efficient methods and to the case of method-of-lines discretizations for nonlinear hyperbolic PDES and systems of PDEs. The following is a summary of the research performed and results achieved during milestone periods 1, 2 and 3 by the PI (Hans De Sterck) and Postdoctoral Research Associate (Oliver Krzysik), for required tasks 1-4 (as listed in the Statement of Work): Research over the previous year has been split into three main projects: (i) solution of acoustic equation system; (ii) solution of nonlinear scalar hyperbolic PDEs; (iii) solution of nonlinear hyperbolic systems of PDEs.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2005

Simplified Solar Modulation Model of Inner Trapped Belt Proton Flux As a Function of Atmospheric Density

No simple algorithm seems to exist for calculating proton fluxes and lifetimes in the Earth's inner, trapped radiation belt throughout the solar cycle. Most models of the inner trapped belt in use depend upon AP8 which only describes the radiation environment at solar maximum and solar minimum in Cycle 20. One exception is NOAAPRO which incorporates flight data from the TIROS/NOAA polar orbiting spacecraft. The present study discloses yet another, simple formulation for approximating proton fluxes at any time in a given solar cycle, in particular between solar maximum and solar minimum. It is derived from AP8 using a regression algorithm technique from nuclear physics. From flux and its time integral fluence, one can then approximate dose rate and its time integral dose. It has already been published in this journal that the absorbed dose rate, D, in the trapped belts exhibits a power law relationship, D = A(rho)(sup -n), where A is a constant, rho is the atmospheric density, and the index n is weakly dependent upon shielding. However, that method does not work for flux and fluence. Instead, we extend this idea by showing that the power law approximation for flux J is actually bivariant in energy E as well as density rho. The resulting relation is J(E,rho)approx.(sum of)A(E(sup n))rho(sup -n), with A itself a power law in E. This provides another method for calculating approximate proton flux and lifetime at any time in the solar cycle. These in turn can be used to predict the associated dose and dose rate.

Wilson, Thomas L.↗
DOE OSTI2019

The multiple depleted mantle components in the Hawaiian-Emperor chain

Oceanic island basalts are targeted for geochemical study because they provide a direct window into mantle composition and a wealth of information on the dynamics and timescales associated with Earth mixing. Previous studies mainly focused on the shield volcanic stage of oceanic islands and the more fusible, enriched mantle components that are easily distinguished in those basalts. Mantle depleted compositions are typically more difficult to resolve unless large amounts of this material participated in mantle melting (e.g., mid-ocean ridges), or unique processes allow for their compositions to be erupted undiluted, such as very small degrees of melting of a source with minimal fusible enriched components (e.g., rejuvenated basalts) or as xenoliths (e.g., abyssal peridotites). Mantle depleted components, defined here as material with low time-integrated Rb/Sr (low 87 Sr/ 86 Sr) and high time-integrated Sm/Nd and Lu/Hf ratios (high 143 Nd/ 144 Nd and 176 Hf/ 177 Hf) relative to primitive mantle, derive from a potentially very large volume reservoir (up to 80% of the mantle), and therefore need adequate characterization in order estimate the composition of the Earth and mantle-derived melts. Here, this review focuses on mantle depleted compositions in oceanic island basalts using the Hawaiian-Emperor chain as a case study. The Hawaiian-Emperor chain is the ~6000 km long geological record of the deeply sourced Hawaiian mantle plume, active for >81 Myr. Hawaiian volcanism evolves through four volcanic stages as a volcano traverses the Hawaiian plume: alkalic preshield, tholeiitic shield (80–90% volcano volume), alkalic postshield (~1%), and silica undersaturated rejuvenated (<0.1%). We report Pb-Sr-Nd-Hf isotope compositions and trace element concentrations of three rejuvenated Northwest Hawaiian Ridge basalts and compare them to an exhaustive compiled dataset of basalts from the Hawaiian Islands to the Emperor Seamounts. The Northwest Hawaiian Ridge (NWHR) includes 51 volcanoes spanning ~42 m.y. between the bend in the Hawaiian-Emperor chain and the Hawaiian Islands where there is no high-precision isotopic data published on the rejuvenated-stage over ~47% of the chain. NWHR and Hawaiian Island rejuvenated basalts are geochemically similar, indicating a consistent source for rejuvenated volcanism over ~12.5 million years. In contrast, shield-stage basalts from the oldest Emperor Seamounts are more depleted in isotopic composition (i.e., higher 176 Hf/ 177 Hf, and 143 Nd/ 144 Nd with lower 87 Sr/ 86 Sr and 208 Pb*/ 206 Pb*) and trace element concentrations (i.e., much lower concentrations of highly incompatible elements) than all other depleted Hawaiian basalts younger than the bend, including NWHR rejuvenated basalts. The strongly depleted source for the oldest Emperor Seamounts (>70 Ma) was likely related to interaction with the Kula-Pacific-Izanagi mid-ocean ridge spreading system active near the Hawaiian plume in the Late Cretaceous. In contrast, the incompatible trace element ratios of NWHR rejuvenated basalts require a distinct source in the Hawaiian mantle plume that was imprinted by ancient (>1 Ga) partial melting, likely ancient recycled oceanic lithosphere. This review of the geochemistry of Hawaiian depleted components documents the need for the sampling of multiple distinctive depleted compositions, each preferentially melted during specific periods of Hawaiian plume activity. This suggests that the composition of depleted components can evolve during the lifetime of the mantle plume, as observed for enriched components in the Hawaiian mantle plume. Changes in the composition of depleted components are dominantly controlled by the upper mantle tectonic configurations at the time of eruption (i.e., proximity to a mid-ocean ridge), as this effect overwhelms the signal imparted by potentially sampling different lower mantle components through time.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2026

Chasing Gamma-Ray Signals from Binary Neutron Star Coalescences with the Cherenkov Telescope Array: Prospects and Observing Strategies

The detection of gravitational waves (GWs) from a binary neutron star (BNS) merger by Advanced LIGO and Advanced Virgo (GW170817), together with its electromagnetic counterpart, the short gamma-ray burst GRB 170817A, heralded the birth of multimessenger astronomy. The detection of TeV emission from GRBs motivates follow-up observations with the Cherenkov Telescope Array Observatory (CTAO), which is ideal for detecting such signals due to its unprecedented sensitivity, rapid response, and wide-field survey capabilities. The aim of this work is to evaluate GeV–TeV GW follow-up strategies for CTAO using a multistep simulation pipeline and to estimate the expected rate of joint GW–GRB detections during observing run O5. Using a simulated sample of BNS systems with corresponding GW detections, gamma-ray emission is simulated through phenomenological prescriptions based on the observed population of short GRBs, including off-axis jet scenarios. CTAO observations are simulated to account for instrument response, sky tiling strategies, integration times, and varying observing conditions. Strategies with variable and constant integration times are investigated. We find that, via an optimized follow-up strategy, about 5% of simulated GW-associated short GRBs produce GeV–TeV radiation detectable by CTAO. Detectability is strongly influenced by the jet opening angle and viewing angle, suggesting that even rough estimates of the viewing angle in GW alerts could enhance targeting. This framework motivates future follow-ups of GW-detectable events, including neutron star–black hole mergers, and further supports the development of advanced strategies incorporating galaxy distributions and synergies with future detectors such as the Einstein Telescope.

Abe, S. [University of Tokyo] (ORCID:0000000172503↗
NASA NTRS2005

Solar Modulation of Inner Trapped Belt Radiation Flux as a Function of Atmospheric Density

No simple algorithm seems to exist for calculating proton fluxes and lifetimes in the Earth's inner, trapped radiation belt throughout the solar cycle. Most models of the inner trapped belt in use depend upon AP8 which only describes the radiation environment at solar maximum and solar minimum in Cycle 20. One exception is NOAAPRO which incorporates flight data from the TIROS/NOAA polar orbiting spacecraft. The present study discloses yet another, simple formulation for approximating proton fluxes at any time in a given solar cycle, in particular between solar maximum and solar minimum. It is derived from AP8 using a regression algorithm technique from nuclear physics. From flux and its time integral fluence, one can then approximate dose rate and its time integral dose.

Lodhi, M. A. K.↗
DOE OSTI2019

IMEX and exact sequence discretization of the multi-fluid plasma model

Multi-fluid plasma models, where an electron fluid is modeled in addition to multiple ion and neutral species as well as the full set of Maxwell's equations, are useful for representing physics beyond the scope of classic MHD. This advantage presents challenges in appropriately dealing with electron dynamics and electromagnetic behavior characterized by the plasma and cyclotron frequencies and the speed of light. For physical systems, such as those near the MHD asymptotic regime, this requirement drastically increases runtimes for explicit time integration even though resolving fast dynamics may not be critical for accuracy. Implicit time integration methods, with efficient solvers, can help to step over fast time-scales that constrain stability, but do not strongly influence accuracy. As an extension, Implicit-explicit (IMEX) schemes provide an additional mechanism to choose which dynamics are evolved using an expensive implicit solve or resolved using a fast explicit solve. In this study, in addition to IMEX methods we also consider a physics compatible exact sequence spatial discretization. Here, this combines nodal bases (H-Grad) for fluid dynamics with a set of vector bases (H-Curl and H-Div) for Maxwell's equations. This discretization allows for multi-fluid plasma modeling without violating Gauss' laws for the electric and magnetic fields. This initial study presents a discussion of the major elements of this formulation and focuses on demonstrating accuracy in the linear wave regime and in the MHD limit for both a visco-resistive and a dispersive ideal MHD problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

Three-dimensional flux-split Euler schemes involving unstructured dynamic meshes

Improved algorithms for the solution of the 3-D time dependent Euler equations are presented for aerodynamic analysis involving unstructured dynamic meshes. The improvements were developed recently to the spatial and temporal discretizations used by unstructured grid flow solvers. The spatial discretization involves a flux split approach which is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves either an explicit time integration scheme using a multistage Runge-Kutta procedure or an implicit time integration scheme using a Gauss-Seidel relaxation procedure which is computationally efficient for either steady or unsteady flow problems. With the implicit Gauss-Seidel procedure, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady flow results are presented for both the NACA 0012 airfoil and the ONERA M6 wing to demonstrate applications of the new Euler solvers. A description of the Euler solvers is presented along with results and comparisons which assess the capability.

Batina, John T.↗
DOE OSTI2024

From Simple Labels to Time-Use Integrations: Supporting the Spectrum of Qualitative Travel Behavior Data

In transportation research, applications of travel behavior data collection are context-specific and require different types of qualitative inputs. These inputs can be viewed as spanning a spectrum of user burden and data quality, from simple trip labels to complex time-use surveys. However, each currently active smartphone-based travel diary platform appears to only support one type of qualitative input, and the effort required for customization is unclear. In this paper, we characterize the spectrum by defining four canonical use cases: (i) trip labels, (ii) trip questionnaire, (iii) counterfactual trips, and (iv) time-use surveys. We then outline a mechanism for supporting configurable user inputs on the same underlying smartphone-based sensing mechanism and demonstrate that it can support all the use cases without any code changes. We further demonstrate that the flexible data model that underpins this mechanism can enable real-time monitoring and analysis. Finally, we evaluate per-user data collection and engagement metrics for large-scale deployments of three canonical use cases, spanning 10 programs, 435 users, and 251,041 trips, and a maximum duration of 800 days. Future efforts may support additional use cases through an expanded configuration and provide greater insight into user engagement. We hope that these insights enable the research community to look at qualitative inputs through a new lens and experiment with novel use cases to fill in the spectrum.

ADVANCED PROPULSION SYSTEMS,POWER TRANSMISSION AND↗
DOE OSTI2026

Measured Indoor Nitrogen Oxides in Homes with a Child with Asthma and Gas or Induction Electric Cooking

The Cooking Energy and Ventilation Impacts on Children's Asthma (CEVICA) study is a randomized control trial investigating the effects of replacing gas with induction electric ranges in the homes of children with asthma in California's San Joaquin Valley. Indoor air quality parameters and respiratory health indicators were measured over 2-week intensive at Baseline and after consecutive 3-month study phases, with half getting electric cooking at the start of Phase 1 and others in Phase 2. Indoor air measurements included time-integrated NO 2 and NO X by passive sampler and time-resolved NO 2 by electrochemical sensors. In the first 40 homes, data were collected during 61 intensives for gas cooking and 55 for induction. Cooking was identified using temperature sensors above the cooktop. Pollutant events were identified from sharp rises in concentrations Time integrated NOX species were lower with electric cooking, with mean differences of 14.2 ppb for NO 2 (95% CI: 9.6, 18.7; p < 0.001), 55.6 ppb for NO X (95% CI: 30.4, 80.9; p < 0.001), and 41.5 ppb for derived NO (95% CI: 19.7, 63.2; p < 0.001). A paired, within home analysis showed larger reductions in Group 2 (gas to induction) than Group 1 that remained electric across Phase 1 and Phase 2 for NO 2 (mean difference in As of 10.5 ppb; 95% CI: 5.6, 15.4; p < 0.001), NO X (32.3 ppb; 95% CI: 10.0, 54.5; p=0.007), and derived NO (21.8 ppb; 95% CI: 3.6, 40.0; p=0.022). Compared to electric cooking, gas had much higher rates of associated NO 2 events, and larger above-baseline NO 2 peaks. These preliminary results are consistent with prior findings that shifting to induction cooking can substantially lower indoor NO 2 compared with gas cooking.

Fang, Yi↗
NASA NTRS1993

Near-field supersonic flow predictions by an adaptive unstructured tetrahedral grid solver

Applicability of a three-dimensional solution adaptive unstructured tetrahedral Euler flow solver about generic models for near-field sonic boom pressure signature predictions is evaluated. Comparisons of computational and experimental data demonstrates the capability of the method for predicting inviscid solutions useful for high speed calculations about simple 3-D geometries. The approach has promising features and results indicate potential for application to more complex configurations. The mesh generation is based on the advancing front technique, and steady state solutions of the Euler equations are achieved by explicit time integration. Spatial discretization uses the Taylor-Galerkin approach; an alternate time integration, based on the Runge-Kutta method, is also included. The solution-adaptive grid procedure is based on either remeshing or mesh refinement techniques. An alternative geometry-adaptive grid procedure has also been incorporated.

Djomehri, M. J.↗
NASA NTRS2006

On Orbit Receiver Performance Assessment of the Geoscience Laser Altimeter System (GLAS) on ICESAT

The GLAS instrument on the NASA's ICESat mission has provided over a billion measurements of the Earth surface elevation and atmosphere backscattering at both 532 and 1064-nm wavelengths. The receiver performance has stayed nearly unchanged since ICESat launch in January 2003. The altimeter receiver has achieved a less than 3-cm ranging accuracy when excluding the effects of the laser beam pointing angle determination uncertainties. The receiver can also detect surface echoes through clouds of one-way transmission as low as 5%. The 532-nm atmosphere backscattering receiver can measure aerosol and clouds with cross section as low as 1e-7/m.sr with a 1 second integration time and molecular backscattering from upper atmosphere with a 60 second integration time. The 1064-nm atmosphere backscattering receiver can measure aerosol and clouds with a cross section as low as 4e-6/m.sr. This paper gives a detailed assessment of the GLAS receiver performance based on the in-orbit calibration tests.

Sun, Xiaoli↗
DOE OSTI2023

A novel integrated time-resolved array avalanche photodiode detection system for nuclear resonant scattering measurements

Here, the nuclear resonant scattering (NRS) experiment requires photon-counting detectors with high time resolution, short dead time, large dynamic range, low noise, and large detection area. An 8-channel avalanche photodiode (APD) array detector system with high integrity, flexibility, and reliability has been developed to adapt to the demands of NRS experiments. The detector system mainly consists of four key parts: (i) an array-APD sensor, (ii) 8-channel integrated fast preamplifiers, (iii) the time-to-digital converter readout electronics, and (iv) a data acquisition system and EPICS support software. Remarkably, the system exhibits a time resolution of better than 500 ps and has a sufficiently low noise level, allowing for the lowest detection energy threshold of 4 keV. The performance of the new array-APD system as well as its real application in nuclear forward scattering (NFS) and nuclear resonant inelastic x-ray scattering (NRIXS) experiments was tested in two synchrotron facilities. With the new system, the NFS signal very close to the prompt electronic scattering signal can be extracted. Thanks to the customized EPICS-areaDetector-based control software, NRIXS spectra can be readily measured with time and energy information of the NRIXS signal stored in the raw data, which is promising for developing NRIXS data analysis in the time domain. The array-APD detector can be deployed for nuclear resonant scattering experiments at various synchrotron radiation facilities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Power balance on a multibeam laser

Inertial confinement fusion (ICF) cryogenic experiments on the 60-beam OMEGA laser have strict requirements for the laser energy delivered on target to be power balanced in order to maximize target-irradiation uniformity. For OMEGA, this quantity (power balance) is inferred from measurements of the time-integrated energy and time-resolved, spatially integrated temporal profile of each of the 60 beams at the output of the laser. The work presented here proposes a general definition of power balance as measured at the laser output and discusses the conditions that are fundamental to achieving laser power balance. Power balance necessitates equal gain across all stages of amplification, equal net losses across each amplifier stage, equal frequency conversion (from 1053 nm to 351 nm) of all 60 beams, and equal beam path lengths (beam timing). Typical OMEGA ICF laser pulse shapes consist of one or more short (100-ps) “pickets” followed by a shaped “drive” pulse of 1 to 2 ns. For these experiments, power balance is assessed for the pickets and the drive independently, with the ultimate goal of achieving root-mean-square (rms) imbalance across all 60 beams of less than 2% rms on both. Our work presents a comprehensive summary of laser shot campaigns conducted to significantly improve laser power balance from typical rms values of 4.7% and 5.2%, respectively, to the 3% level for both features along with a discussion of future work required to further reduce the rms power imbalance of the laser system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗