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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

The area-time-integral technique to estimate convective rain volumes over areas applied to satellite data - A preliminary investigation

The use of the area-time-integral (ATI) technique, based only on satellite data, to estimate convective rain volume over a moving target is examined. The technique is based on the correlation between the radar echo area coverage integrated over the lifetime of the storm and the radar estimated rain volume. The processing of the GOES and radar data collected in 1981 is described. The radar and satellite parameters for six convective clusters from storm events occurring on June 12 and July 2, 1981 are analyzed and compared in terms of time steps and cluster lifetimes. Rain volume is calculated by first using the regression analysis to generate the regression equation used to obtain the ATI; the ATI versus rain volume relation is then employed to compute rain volume. The data reveal that the ATI technique using satellite data is applicable to the calculation of rain volume.

Doneaud, Andre A.↗
DOE OSTI2026

New Time Integrators and Capabilities in SUNDIALS Versions 6.2.0-7.4.0

SUNDIALS is a well-established numerical library that provides robust and efficient time integrators and nonlinear solvers. This article overviews several significant improvements and new features added over the last 3 years to support scientific simulations run on high-performance computing systems. Notably, three new classes of one-step methods have been implemented: low storage Runge–Kutta, symplectic partitioned Runge–Kutta, and operator splitting. In addition, we describe new timestep adaptivity support for multirate methods, adjoint sensitivity analysis capabilities for explicit Runge–Kutta methods, additional options for Anderson acceleration in nonlinear solvers, and improved error handling and logging.

Computer science↗
DOE OSTI2025

Measurement of the time-integrated C P asymmetry in D 0 → K S 0 K S 0 decays using Belle and Belle II data

We measure the time-integrated C P asymmetry in D 0 → K S 0 K S 0 decays reconstructed in e + e − → c c ¯ events collected by the Belle and Belle II experiments. The corresponding data samples have integrated luminosities of 980 and 428 fb − 1 , respectively. The D 0 decays are required to originate from the D * + → D 0 π + decay, which determines the charm flavor at production time. A control sample of D 0 → K + K − decays is used to correct for production and detection asymmetries. The result, ( − 1.4 ± 1.3 ( stat ) ± 0.1 ( syst ) ) % , is consistent with previous determinations and with C P symmetry. Published by the American Physical Society 2025

Adachi, I. (ORCID:0000000322870173)↗
NASA NTRS1993

Further studies to extend and test the area-time-integral technique applied to satellite data

The principal goal of this project is to establish relationships that would allow application of area-time integral (ATI) calculations based upon satellite data to estimate rainfall volumes. The research has been pursued using two different approaches, which for convenience can be designated as the 'fixed-threshold approach' and the 'variable-threshold approach'. In the former approach, an attempt is made to determine a single temperature threshold in the satellite infrared data that would yield ATI values for identifiable cloud clusters which are most closely related to the corresponding rainfall amounts. Results thus far have indicated that a strong correlation exists between the rain volumes and the satellite ATI values, but the optimum threshold for this relationship seems to differ from one geographic location to another. The difference is probably related to differences in the basic precipitation mechanisms that dominate in the different regions. The average rainfall rate associated with each cloudy pixel is also found to vary across the spectrum of ATI values. Work on the second, or 'variable-threshold', approach for determining the satellite ATI values was essentially suspended during this period due to exhaustion of project funds. Most of the ATI work thus far has dealt with cloud clusters from the Lagrangian or 'floating-target' point of view. For many purposes, however, the Eulerian or 'fixed-target' perspective is more appropriate. For a very large target area encompassing entire cluster life histories, the rain volume-ATI relationship would obviously be the same in either case. The important question for the Eulerian perspective is how small the fixed area can be made while maintaining consistency in that relationship.

Smith, Paul L.↗
NASA NTRS1993

Direct time integration of Maxwell's equations in two-dimensional dielectric waveguides for propagation and scattering of femtosecond electromagnetic solitons

We present what are to our knowledge first-time calculations from vector nonlinear Maxwell's equations of femtosecond soliton propagation and scattering, including carrier waves, in two-dimensional dielectric waveguides. The time integration efficiently implements linear and nonlinear convolutions for the electric polarization, and the nonlinear convolution accounts for two quantum effects, the Kerr and Raman interactions. By retaining the optical carrier, the new method solves for fundamental quantities - optical electric and magnetic fields in space and time - rather than a nonphysical envelope function. It has the potential to provide an unprecedented two- and three-dimensional modeling capability for millimeter-scale integrated-optical circuits with submicrometer engineered inhomogeneities.

Joseph, Rose M.↗
DOE OSTI2024

SUNDIALS time integrators for exascale applications with many independent systems of ordinary differential equations

Many complex systems can be accurately modeled as a set of coupled time-dependent partial differential equations (PDEs). However, solving such equations can be prohibitively expensive, easily taxing the world’s largest supercomputers. One pragmatic strategy for attacking such problems is to split the PDEs into components that can more easily be solved in isolation. This operator splitting approach is used ubiquitously across scientific domains, and in many cases leads to a set of ordinary differential equations (ODEs) that need to be solved as part of a larger “outer-loop” time-stepping approach. The SUNDIALS library provides a plethora of robust time integration algorithms for solving ODEs, and the U.S. Department of Energy Exascale Computing Project (ECP) has supported its extension to applications on exascale-capable computing hardware. In this paper, we highlight some SUNDIALS capabilities and its deployment in combustion and cosmology application codes (Pele and Nyx, respectively) where operator splitting gives rise to numerous, small ODE systems that must be solved concurrently.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1996

Investigation of Adaptive-threshold Approaches for Determining Area-Time Integrals from Satellite Infrared Data to Estimate Convective Rain Volumes

The principal goal of this project is to establish relationships that would allow application of area-time integral (ATI) calculations based upon satellite data to estimate rainfall volumes. The research is being carried out as a collaborative effort between the two participating organizations, with the satellite data analysis to determine values for the ATIs being done primarily by the STC-METSAT scientists and the associated radar data analysis to determine the 'ground-truth' rainfall estimates being done primarily at the South Dakota School of Mines and Technology (SDSM&T). Synthesis of the two separate kinds of data and investigation of the resulting rainfall-versus-ATI relationships is then carried out jointly. The research has been pursued using two different approaches, which for convenience can be designated as the 'fixed-threshold approach' and the 'adaptive-threshold approach'. In the former, an attempt is made to determine a single temperature threshold in the satellite infrared data that would yield ATI values for identifiable cloud clusters which are closely related to the corresponding rainfall amounts as determined by radar. Work on the second, or 'adaptive-threshold', approach for determining the satellite ATI values has explored two avenues: (1) attempt involved choosing IR thresholds to match the satellite ATI values with ones separately calculated from the radar data on a case basis; and (2) an attempt involved a striaghtforward screening analysis to determine the (fixed) offset that would lead to the strongest correlation and lowest standard error of estimate in the relationship between the satellite ATI values and the corresponding rainfall volumes.

Smith, Paul L.↗
DOE OSTI2025

Measurement of the time-integrated 𝐶⁢𝑃 asymmetry in 𝐷 0 → 𝐾$^{0}_{S}$𝐾$^{0}_{S}$ decays using opposite-side flavor tagging at Belle and Belle II

We measure the time-integrated 𝐶⁢𝑃 asymmetry in 𝐷 0 → 𝐾$^{0}_{S}$𝐾$^{0}_{S}$ decays reconstructed in 𝑒 + ⁢𝑒 − → $c\bar{c}$ events collected by the Belle and Belle II experiments. The corresponding data samples have integrated luminosities of 980 and 428 fb −1 , respectively. To infer the flavor of the 𝐷 0 meson, we exploit the correlation between the flavor of the reconstructed decay and the electric charges of particles reconstructed in the rest of the 𝑒 + ⁢𝑒 − → $c\bar{c}$ event. This results in a sample which is independent from any other previously used at Belle or Belle II. The result, 𝐴 𝐶⁢𝑃 ⁡(𝐷 0 → 𝐾$^{0}_{S}$𝐾$^{0}_{S}$)=(1.3±2.0±0.2)%, where the first uncertainty is statistical and the second systematic, is consistent with previous determinations and with 𝐶⁢𝑃 symmetry.

CP violation↗
NASA NTRS1991

Direct time integration of Maxwell's equations in linear dispersive media with absorption for scattering and propagation of femtosecond electromagnetic pulses

The initial results for femtosecond pulse propagation and scattering interactions for a Lorentz medium obtained by a direct time integration of Maxwell's equations are reported. The computational approach provides reflection coefficients accurate to better than 6 parts in 10,000 over the frequency range of dc to 3 x 10 to the 16th Hz for a single 0.2-fs Gaussian pulse incident upon a Lorentz-medium half-space. New results for Sommerfeld and Brillouin precursors are shown and compared with previous analyses. The present approach is robust and permits 2D and 3D electromagnetic pulse propagation directly from the full-vector Maxwell's equations.

Joseph, Rose M.↗
DOE OSTI2025

Time-domain all-frequency stable formulation for low-frequency electromagnetic simulation with Newmark-β time integration

An implicitly Coulomb-gauged A-ϕ formulation has previously been proposed and validated for finite ele- ment simulations of low-frequency and multiscale electromag- netic problems in the frequency domain. This formulation has demonstrated numerical stability across all frequencies, with its accuracy, efficiency, and iterative convergence established in various frequency-domain scenarios. However, direct time- domain computation is often preferable for wideband electro- magnetic problems and is typically indispensable in nonlinear and multiphysics simulations. In this work, the A-ϕ formulation is extended to the time domain. By incorporating the well-known Newmark-β time integration scheme, the proposed formulation is validated through capacitive and inductive test cases. The results confirm the solution’s accuracy and demonstrate the formulation’s stability in the time domain.

Mekonnen, Minyechil↗
NASA NTRS2003

A Time Integration Algorithm Based on the State Transition Matrix for Structures with Time Varying and Nonlinear Properties

A variable order method of integrating the structural dynamics equations that is based on the state transition matrix has been developed. The method has been evaluated for linear time variant and nonlinear systems of equations. When the time variation of the system can be modeled exactly by a polynomial it produces nearly exact solutions for a wide range of time step sizes. Solutions of a model nonlinear dynamic response exhibiting chaotic behavior have been computed. Accuracy of the method has been demonstrated by comparison with solutions obtained by established methods.

Bartels, Robert E.↗
NASA NTRS2002

Higher Order Time Integration Schemes for the Unsteady Navier-Stokes Equations on Unstructured Meshes

The efficiency gains obtained using higher-order implicit Runge-Kutta schemes as compared with the second-order accurate backward difference schemes for the unsteady Navier-Stokes equations are investigated. Three different algorithms for solving the nonlinear system of equations arising at each timestep are presented. The first algorithm (NMG) is a pseudo-time-stepping scheme which employs a non-linear full approximation storage (FAS) agglomeration multigrid method to accelerate convergence. The other two algorithms are based on Inexact Newton's methods. The linear system arising at each Newton step is solved using iterative/Krylov techniques and left preconditioning is used to accelerate convergence of the linear solvers. One of the methods (LMG) uses Richardson's iterative scheme for solving the linear system at each Newton step while the other (PGMRES) uses the Generalized Minimal Residual method. Results demonstrating the relative superiority of these Newton's methods based schemes are presented. Efficiency gains as high as 10 are obtained by combining the higher-order time integration schemes with the more efficient nonlinear solvers.

Jothiprasad, Giridhar↗
DOE OSTI2026

Multibody for Everybody (M4E) - A Linearization Approach to Enable Frequency Domain Analysis, Time Integration and Control Co-Design

1.1 Background/Objectives: Marine energy represents a promising yet underexploited source of power. To increase the harvested power, significant efforts have been made to improve wave energy converter (WEC) modeling capabilities and optimize power take-off (PTO) performance; however, these efforts have often treated WEC dynamics, PTO design, and controller development sequentially. In contrast, control co-design (CCD) is emerging as a promising strategy to address these issues directly, creating a growing need for fast analysis tools suitable for repeated simulation and parametric studies [1]. To support this need, this work presents the Multibody for Everybody (M4E) [2] linearization module, which employs a symbolic toolbox to provide deeper insight of WEC design parameters. The objective is to demonstrate that a minimal-coordinate linearization of articulated WEC dynamics can provide accurate wave response predictions and substantial computational savings relative to nonlinear time-domain simulation, while preserving compatibility with broader wave-energy analysis workflows, enabling CCD. 1.2 Approach/Activities: The proposed approach linearizes the equations of motion, generated by M4E, in minimal coordinates about a selected operating point and combines the resulting system with frequencydomain hydrodynamic terms to incorporate the reduced mass, damping, stiffness, and forcing operators. The linearized model is used for both impedance-based response amplitude operator (RAO) prediction and rapid regular-wave time integration. The methodology is demonstrated on a single-flap device and a FOSWEC configuration, with linearized M4E responses compared against the corresponding nonlinear M4E simulations and WEC-Sim results. Regular-wave time histories, RAO trends, and runtime differences are assessed. The framework is also compatible with broader wave-energy workflows, including coupling to WecOptTool, although that capability is not the focus of this work [3]. 1.3 Results/Lessons: The linearized M4E model reproduces key regularwave response characteristics such as integration and Response Amplitude over multiple frequencies. This module matches nonlinear M4E and WEC-Sim results while substantially reducing integration cost. Thus, the proposed framework can serve as a rapid analysis layer for articulated WEC design, parameter studies, and controls-oriented workflows. The analysis is most appropriate in the near-equilibrium regime, about the linearization point.

16 TIDAL AND WAVE POWER↗
NASA NTRS1989

A raingage, radar and satellite simulation study of the estimation of convective rainfall by area-time integrals

Correlations between area-averaged rainrates and the areal fraction of rainrates exceeding a preset threshold, F(T), are examined using data from a network of 22 raingages located within 120 km of the NOAA/Tropical Ocean-Global Atmosphere (TOGA) C-band meteorological radar at Darwin, Australia. The results show that the area averages of convective rainfall are highly correlated with the fraction F(T). To simulate the use of the relationship between area-averaged rainrates and F(T) to obtain space-time averaged rainrates from a satellite sensor, observations from the NOAA/TOGA radar observations are used to estimate F(T) using a reflectivity threshold. The use of area-time integral methods for inferring area-averaged rainrates from satellites is examined, noting the possible use of the methodology for the Tropical Rainfall Measuring Mission.

Short, David A.↗
NASA NTRS1986

Sequential ranging integration times in the presence of CW interference in the ranging channel

The Deep Space Network (DSN), managed by the Jet Propulsion Laboratory for NASA, is used primarily for communication with interplanetary spacecraft. The high sensitivity required to achieve planetary communications makes the DSN very susceptible to radio-frequency interference (RFI). In this paper, an analytical model is presented of the performance degradation of the DSN sequential ranging subsystem in the presence of downlink CW interference in the ranging channel. A trade-off between the ranging component integration times and the ranging signal-to-noise ratio to achieve a desired level of range measurement accuracy and the probability of error in the code components is also presented. Numerical results presented illustrate the required trade-offs under various interference conditions.

Mathur, Ashok↗
NASA NTRS1997

Massively Parallel and Scalable Implicit Time Integration Algorithms for Structural Dynamics

Explicit codes are often used to simulate the nonlinear dynamics of large-scale structural systems, even for low frequency response, because the storage and CPU requirements entailed by the repeated factorizations traditionally found in implicit codes rapidly overwhelm the available computing resources. With the advent of parallel processing, this trend is accelerating because of the following additional facts: (a) explicit schemes are easier to parallelize than implicit ones, and (b) explicit schemes induce short range interprocessor communications that are relatively inexpensive, while the factorization methods used in most implicit schemes induce long range interprocessor communications that often ruin the sought-after speed-up. However, the time step restriction imposed by the Courant stability condition on all explicit schemes cannot yet be offset by the speed of the currently available parallel hardware. Therefore, it is essential to develop efficient alternatives to direct methods that are also amenable to massively parallel processing because implicit codes using unconditionally stable time-integration algorithms are computationally more efficient when simulating the low-frequency dynamics of aerospace structures.

Farhat, Charbel↗
NASA NTRS2001

Designing Adaptive Low-Dissipative High Order Schemes for Long-Time Integrations

A general framework for the design of adaptive low-dissipative high order schemes is presented. It encompasses a rather complete treatment of the numerical approach based on four integrated design criteria: (1) For stability considerations, condition the governing equations before the application of the appropriate numerical scheme whenever it is possible; (2) For consistency, compatible schemes that possess stability properties, including physical and numerical boundary condition treatments, similar to those of the discrete analogue of the continuum are preferred; (3) For the minimization of numerical dissipation contamination, efficient and adaptive numerical dissipation control to further improve nonlinear stability and accuracy should be used; and (4) For practical considerations, the numerical approach should be efficient and applicable to general geometries, and an efficient and reliable dynamic grid adaptation should be used if necessary. These design criteria are, in general, very useful to a wide spectrum of flow simulations. However, the demand on the overall numerical approach for nonlinear stability and accuracy is much more stringent for long-time integration of complex multiscale viscous shock/shear/turbulence/acoustics interactions and numerical combustion. Robust classical numerical methods for less complex flow physics are not suitable or practical for such applications. The present approach is designed expressly to address such flow problems, especially unsteady flows. The minimization of employing very fine grids to overcome the production of spurious numerical solutions and/or instability due to under-resolved grids is also sought. The incremental studies to illustrate the performance of the approach are summarized. Extensive testing and full implementation of the approach is forthcoming. The results shown so far are very encouraging.

Yee, Helen C.↗
NASA NTRS2005

Designing Adaptive Low-Dissipative High Order Schemes for Long-Time Integrations

A general framework for the design of adaptive low-dissipative high order schemes is presented. It encompasses a rather complete treatment of the numerical approach based on four integrated design criteria: (1) For stability considerations, condition the governing equations before the application of the appropriate numerical scheme whenever it is possible. (2) For consistency, compatible schemes that possess stability properties, including physical and numerical boundary condition treatments, similar to those of the discrete analogue of the continuum are preferred. (3) For the minimization of numerical dissipation contamination, efficient and adaptive numerical dissipation control to further improve nonlinear stability and accuracy should be used. (4) For practical considerations, the numerical approach should be efficient and applicable to general geometries, and an efficient and reliable dynamic grid adaptation should be used if necessary. These design criteria are, in general, very useful to a wide spectrum of flow simulations. However, the demand on the overall numerical approach for non-linear stability and accuracy is much more stringent for long-time integration of complex multiscale viscous shock/shear/turbulence/acoustics interactions and numerical combustion. Robust classical numerical methods for less complex flow physics are not suitable or practical for such applications. The present approach is designed expressly to address such flow problems, especially unsteady flows.

Yee, H. C.↗