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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Rythmos: Solution and Analysis Package for Differential-Algebraic and Ordinary-Differential Equations

Time integration is a central component for most transient simulations. It coordinates many of the major parts of a simulation together, e.g., a residual calculation with a transient solver, solution with the output, various operator-split physics, and forward and adjoint solutions for inversion. Even though there is this variety in these transient simulations, there is still a common set of algorithms and procedures to progress transient solutions for ordinary-differential equations (ODEs) and differential-alegbraic equations (DAEs). Rythmos is a collection of these algorithms that can be used for the solution of transient simulations. It provides common time-integration methods, such as Backward and Forward Euler, Explicit and Implicit Runge-Kutta, and Backward-Difference Formulas. It can also provide sensitivities, and adjoint components for transient simulations. Rythmos is a package within Trilinos, and requires some other packages (e.g., Teuchos and Thrya) to provide basic time-integration capabilities. It also can be coupled with several other Trilinos packages to provide additional capabilities (e.g., AztecOO and Belos for linear solutions, and NOX for non-linear solutions). The documentation is broken down into three parts: Theory Manual, User's Manual, and Developer's Guide. The Theory Manual contains the basic theory of the time integrators, the nomenclature and mathematical structure utilized within Rythmos, and verification results demonstrating that the designed order of accuracy is achieved. The User's Manual provides information on how to use the Rythmos, description of input parameters through Teuchos Parameter Lists, and description of convergence test examples. The Developer's Guide is a high-level discussion of the design and structure of Rythmos to provide information to developers for the continued development of capabilities. Details of individual components can be found in the Doxygen webpages.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

SIX NEW MILLISECOND PULSARS FROM ARECIBO SEARCHES OF FERMI GAMMA-RAY SOURCES

We have discovered six radio millisecond pulsars (MSPs) in a search with the Arecibo telescope of 34 unidentified gamma-ray sources from the Fermi Large Area Telescope (LAT) four year point source catalog. Among the 34 sources, we also detected two MSPs previously discovered elsewhere. Each source was observed at a center frequency of 327 MHz, typically at three epochs with individual integration times of 15 minutes. The new MSP spin periods range from 1.99 to 4.66 ms. Five of the six pulsars are in interacting compact binaries (period ≤ 8.1 hr), while the sixth is a more typical neutron star-white dwarf binary with an 83 day orbital period. This is a higher proportion of interacting binaries than for equivalent Fermi-LAT searches elsewhere. The reason is that Arecibo's large gain afforded us the opportunity to limit integration times to 15 minutes, which significantly increased our sensitivity to these highly accelerated systems. Seventeen of the remaining 26 gamma-ray sources are still categorized as strong MSP candidates, and will be re-searched.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Implicit integration methods for dislocation dynamics

In dislocation dynamics simulations, strain hardening simulations require integrating stiff systems of ordinary differential equations in time with expensive force calculations, discontinuous topological events, and rapidly changing problem size. Current solvers in use often result in small time steps and long simulation times. Faster solvers may help dislocation dynamics simulations accumulate plastic strains at strain rates comparable to experimental observations. Here, this paper investigates the viability of high order implicit time integrators and robust nonlinear solvers to reduce simulation run times while maintaining the accuracy of the computed solution. In particular, implicit Runge-Kutta time integrators are explored as a way of providing greater accuracy over a larger time step than is typically done with the standard second-order trapezoidal method. In addition, both accelerated fixed point and Newton's method are investigated to provide fast and effective solves for the nonlinear systems that must be resolved within each time step. Results show that integrators of third order are the most effective, while accelerated fixed point and Newton's method both improve solver performance over the standard fixed point method used for the solution of the nonlinear systems.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Multiple time step integrators in ab initio molecular dynamics

Multiple time-scale algorithms exploit the natural separation of time-scales in chemical systems to greatly accelerate the efficiency of molecular dynamics simulations. Although the utility of these methods in systems where the interactions are described by empirical potentials is now well established, their application to ab initio molecular dynamics calculations has been limited by difficulties associated with splitting the ab initio potential into fast and slowly varying components. Here we present two schemes that enable efficient time-scale separation in ab initio calculations: one based on fragment decomposition and the other on range separation of the Coulomb operator in the electronic Hamiltonian. We demonstrate for both water clusters and a solvated hydroxide ion that multiple time-scale molecular dynamics allows for outer time steps of 2.5 fs, which are as large as those obtained when such schemes are applied to empirical potentials, while still allowing for bonds to be broken and reformed throughout the dynamics. This permits computational speedups of up to 4.4x, compared to standard Born-Oppenheimer ab initio molecular dynamics with a 0.5 fs time step, while maintaining the same energy conservation and accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

Rain rate and modeled fade distributions at 20 GHz and 30 GHz derived from five years of network rain gauge measurements

Five years of rain rate and modeled slant path attenuation distributions at 20 GHz and 30 GHz derived from a network of 10 tipping bucket rain gages was examined. The rain gage network is located within a grid 70 km north-south and 47 km east-west in the Mid-Atlantic coast of the United States in the vicinity of Wallops Island, Virginia. Distributions were derived from the variable integration time data and from one minute averages. It was demonstrated that for realistic fade margins, the variable integration time results are adequate to estimate slant path attenuations at frequencies above 20 GHz using models which require one minute averages. An accurate empirical formula was developed to convert the variable integration time rain rates to one minute averages. Fade distributions at 20 GHz and 30 GHz were derived employing Crane's Global model because it was demonstrated to exhibit excellent accuracy with measured COMSTAR fades at 28.56 GHz.

Goldhirsh, Julius↗
DOE OSTI2022

Validation of energetic particle transport models for time-dependent integrated simulations of burning plasmas (Final ST Report)

The major goal of this project is to deliver a suite of validated energetic particle (EP) modeling tools for time-dependent integrated tokamak simulations (with the code TRANSP) that can be used for modeling plasma discharges featuring a substantial EP content, such as burning plasmas. During FY20, issues were resolved with an energetic particle transport model called TGLF-EP, with improvements made to the TGLF-EP workflow, generation of input files, and interpretation of runs on the NERSC supercomputing cluster. The recipient also gave input on simulations with the code TRANSP to reproduce the reference JET discharge, and assessed the accuracy and uncertainties of the simulated B field, current evolution, thermal plasma profiles, stored energy, etc. Input was also given on interpretation of fast-ion instability calculations and recommendations for modifications to JET afterglow discharge for the upcoming DT campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1985

Optical detection of large meteoroids in space

CCD sensors placed across the focal plane of large Schmidt telescopes have great potential for detecting and measuring the very low flux in space of meteoroids with diameters larger than 1 meter. With the Palomar 'Big Schmidt', a detection rate of 1.4 per hour is obtained for meteoroids between 0.6 and 200 meters in diameter. For the Baker-Nunn 'Satellite Tracking Camera', the corresponding rate is about 0.8 per hour. The key to obtaining such high detection rates derives from approximately setting the sensor integration time equal to the time it takes a meteoroid to cross a pixel field of view. This minimizes signal to noise problems and is accomplished, in practice, by multiple summing of short integration time data records to obtain data records of longer effective integration times.

Zook, H. A.↗
NASA NTRS1975

Design Data Collection with Skylab Microwave Radiometer-Scatterometer S-193, Volume 1

The author has identified the following significant results. Observations with S-193 have provided radar design information for systems to be flown on spacecraft, but only at 13.9 GHz and for land areas over the United States and Brazil plus a few other areas of the world for which this kind of analysis was not made. Observations only extended out to about 50 deg angle of incidence. The value of a sensor with such a gross resolution for most overland resource and status monitoring systems seems marginal, with the possible exception of monitoring soil moisture and major vegetation variations. The complementary nature of the scatterometer and radiometer systems was demonstrated by the correlation analysis. Although radiometers must have spatial resolutions dictated by antenna size, radars can use synthetic aperture techniques to achieve much finer resolutions. Multiplicity of modes in the S-193 sensors complicated both the system development and its employment. An attempt was made in the design of the S-193 to arrange optimum integration times for each angle and type of measurement. This unnecessarily complicated the design of the instrument, since the gains in precision achieved in this way were marginal. Either a software-controllable integration time or a set of only two or three integration times would have been better.

Moore, R. K.↗
NASA NTRS1989

Photon statistical limitations for daytime optical tracking

Tracking of interplanetary spacecraft equipped with optical communication systems by using astrometric instruments is being investigated by JPL. Existing instruments are designed to work at night and, for bright sources, are limited by tropospheric errors. To provide full coverage of the solar system, astrometric tracking instruments must either be capable of daytime operation or be space-based. The integration times necessary for the ground-based daytime photon statistical errors to reach a given accuracy level (5 to 50 nanoradians) were computed for an ideal astrometric instrument. The required photon statistical integration times are found to be shorter than the tropospheric integrations times for the ideal detector. Since the astrometric need not be limited by photon statistics even under daytime conditions, it may be fruitful to investigate instruments for daytime optical tracking.

Folkner, W. M.↗
NASA NTRS2016

Characterizing Rocky and Gaseous Exoplanets with 2 m Class Space-based Coronagraphs

Several concepts now exist for small, space-based missions to directly characterize exoplanets in reflected light. While studies have been performed that investigate the potential detection yields of such missions, little work has been done to understand how instrumental and astrophysical parameters will affect the ability of these missions to obtain spectra that are useful for characterizing their planetary targets. Here, we develop an instrument noise model suitable for studying the spectral characterization potential of a coronagraph-equipped, space-based telescope. We adopt a baseline set of telescope and instrument parameters appropriate for near-future planned missions like WFIRST-AFTA, including a 2 m diameter primary aperture, an operational wavelength range of 0.4–1.0 μm, and an instrument spectral resolution of λ/Δλ =70, and apply our baseline model to a variety of spectral models of different planet types, including Earth twins, Jupiter twins, and warm and cool Jupiters and Neptunes. With our exoplanet spectral models, we explore wavelength-dependent planet–star flux ratios for main-sequence stars of various effective temperatures and discuss how coronagraph inner and outer working angle constraints will influence the potential to study different types of planets. For planets most favorable to spectroscopic characterization—cool Jupiters and Neptunes as well as nearby super-Earths—we study the integration times required to achieve moderate signal-to-noise ratio spectra. We also explore the sensitivity of the integration times required to either detect the bottom or presence of key absorption bands (for methane, water vapor, and molecular oxygen) to coronagraph raw contrast performance, exozodiacal light levels, and the distance to the planetary system. Decreasing detector quantum efficiency at longer visible wavelengths makes the detection of water vapor in the atmospheres of Earth-like planets extremely challenging, and also hinders detections of the 0.89 μm methane band. Additionally, most modeled observations have noise dominated by dark currents, indicating that improving CCD performance could substantially drive down requisite integration times. Finally, we briefly discuss the extension of our models to a more distant future Large UV-Optical-InfraRed (LUVOIR) mission.

Astrobiology↗
DOE OSTI2022

The Schwarz alternating method for transient solid dynamics

Abstract In our earlier work, we formulated the Schwarz alternating method as a means for concurrent multiscale coupling in finite deformation solid mechanics for quasi‐static problems. Herein, we advance this method for the study of transient dynamic multiscale solid mechanics problems where information is exchanged back and forth between small and large scales. The extension to dynamics relies on the notion of a global time stepper. Within each global time step, the subdomains are coupled by the standard Schwarz iterative process. Remarkably, each subdomain can use its own time step or even its own time integrator to advance its solution in time, provided that they synchronize at each global time step. We study the performance of the Schwarz method on several examples designed for this purpose. Our numerical experiments demonstrate that the method is capable of coupling regions with different mesh resolutions, different element types, and different time integration schemes (e.g., implicit and explicit), all without introducing any artifacts that afflict other coupling methods for transient dynamics. Finally, we apply the dynamic Schwarz alternating method to the simulation of a bolted joint subjected to dynamic loading, as a demonstration of the performance of the method in a realistic scenario.

Mota, Alejandro↗
NASA NTRS2010

Adaptive Numerical Algorithms in Space Weather Modeling

Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by different physics in different domains. A multi-physics system can be modeled by a software framework comprising of several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and efficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solar wind Roe Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamics (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a block-adaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit numerical schemes. Depending on the application, we find that different time stepping methods are optimal. Several of the time integration schemes exploit the block-based granularity of the grid structure. The framework and the adaptive algorithms enable physics based space weather modeling and even forecasting.

Toth, Gabor↗
DOE OSTI2024

Spectral Analysis of Regular Material Point Method and its Application to Study High Pressure Reverse Osmosis Membrane Compaction and Embossing

Material Point Method (MPM) is gaining widespread interest in applied continuum mechanics. The fact that all the continuum properties are stored on the particles (or material points) and the governing equations are solved on these material points makes MPM extremely suited to problems involving severe material deformations, such as crack propagation, soil movement, and fluid flows. Despite its popularity, only a few studies have focused on the numerical properties of MPM. This presentation introduces a global spectral analysis of the regular material point method. Contrary to previous studies, the analysis focuses on the numerical properties of the method in the spectral space. The amplification factor is derived as a function of the non- dimensional wave numbers. It provides insights into the stability and dissipative properties of the method for various CFL and Fourier numbers. The effect of the grid shape functions, number of particles per cell and their locations inside the grid cell are also analyzed. The EXAGOOP MPM solver (https://github.com/NREL/Exagoop.git) is developed at the National Renewable Energy Laboratory as a part of the NAWI UHPRO project and is based on the AMReX framework. A single-level, uniform cartesian grid is used as the background mesh, while the particle class in AMReX is used to manage the material point operations. Linear hat and B-splines are used as grid shape functions, while the time integration is performed using explicit Euler time integration. EXAGOOP is both CPU and GPU compatible and has been demonstrated to work well on multiple compute architectures. The performance of EXAGOOP on various computing architectures is presented along with its application to study compaction and embossing of high-pressure reverse osmosis membranes. The MPM solution accurately reproduces the membrane deformation. The deformed pore size and structure simulated using MPM also agree well with experimental SEM images.

material point method↗
DOE OSTI2015

Quantum stopping times stochastic integral in the interacting Fock space

Following the ideas of Hudson [J. Funct. Anal. 34(2), 266-281 (1979)] and Parthasarathy and Sinha [Probab. Theory Relat. Fields 73, 317-349 (1987)], we define a quantum stopping time (QST, for short) τ in the interacting Fock space (IFS, for short), Γ, over L{sup 2}(ℝ{sup +}), which is actually a spectral measure in [0, ∞] such that τ([0, t]) is an adapted process. Motivated by Parthasarathy and Sinha [Probab. Theory Relat. Fields 73, 317-349 (1987)] and Applebaum [J. Funct. Anal. 65, 273-291 (1986)], we also develop a corresponding quantum stopping time stochastic integral (QSTSI, for abbreviations) on the IFS over a subspace of L{sup 2}(ℝ{sup +}) equipped with a filtration. As an application, such integral provides a useful tool for proving that Γ admits a strong factorisation, i.e., Γ = Γ{sub τ]} ⊗ Γ{sub [τ}, where Γ{sub τ]} and Γ{sub [τ} stand for the part “before τ” and the part “after τ,” respectively. Additionally, this integral also gives rise to a natural composition operation among QST to make the space of all QSTs a semigroup.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Dynamic calibration of differential equations using machine learning, with application to turbulence models

We present a methodology for calibration of parametric ordinary and partial differential equation models, using off-the-shelf software for back-propagation in Neural Networks (NN). As a prototypical example, we consider calibration of a Reynolds-averaged Navier-Stokes (RANS) turbulence closure model, against ground truth data from direct numerical simulations (DNS) of two different turbulent flows. Numerical time integration is represented as a custom NN, where only the RANS model parameters are trainable. A loss function is defined to quantify the mismatch between the NN prediction and the ground truth over a predefined, finite time integration window. This loss function is then minimized using a gradient descent method utilizing the back-propagation algorithm. Furthermore, this dynamic approach to training is to be contrasted with a static approach, wherein a least square regression estimate for parameters is obtained in the limit of an infinitesimal time integration window. In a first test of static and dynamic approaches against ground truth data generated by the model, the former proves to be significantly faster and more accurate than the latter at recovering the parameters. When both calibration approaches are tested against DNS data, for which it is known that the model cannot achieve a perfect fit, the static approach yields a good prediction only for short times, while the dynamic approach results in physical and stable predictions over the entire integration window. After optimization of the dynamic approach for time step, spatial resolution, stability, and physics-based constraints, we obtain a 50% improvement of outcomes over those obtained from the existing, manually calibrated set of parameters, demonstrating the merits of this systematic and automated procedure.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

PRISMS-PF: A general framework for phase-field modeling with a matrix-free finite element method

Abstract A new phase-field modeling framework with an emphasis on performance, flexibility, and ease of use is presented. Foremost among the strategies employed to fulfill these objectives are the use of a matrix-free finite element method and a modular, application-centric code structure. This approach is implemented in the new open-source PRISMS-PF framework. Its performance is enabled by the combination of a matrix-free variant of the finite element method with adaptive mesh refinement, explicit time integration, and multilevel parallelism. Benchmark testing with a particle growth problem shows PRISMS-PF with adaptive mesh refinement and higher-order elements to be up to 12 times faster than a finite difference code employing a second-order-accurate spatial discretization and first-order-accurate explicit time integration. Furthermore, for a two-dimensional solidification benchmark problem, the performance of PRISMS-PF meets or exceeds that of phase-field frameworks that focus on implicit/semi-implicit time stepping, even though the benchmark problem’s small computational size reduces the scalability advantage of explicit time-integration schemes. PRISMS-PF supports an arbitrary number of coupled governing equations. The code structure simplifies the modification of these governing equations by separating their definition from the implementation of the numerical methods used to solve them. As part of its modular design, the framework includes functionality for nucleation and polycrystalline systems available in any application to further broaden the phenomena that can be used to study. The versatility of this approach is demonstrated with examples from several common types of phase-field simulations, including coarsening subsequent to spinodal decomposition, solidification, precipitation, grain growth, and corrosion.

36 MATERIALS SCIENCE↗
NASA NTRS1988

Multi-time-step integration using nodal partitioning

An algorithm is presented which integrates different groups of nodes of a finite element mesh with different time steps and different integrators. Since the nodal groups are updated independently no unsymmetric systems need be solved. Stability is demonstrated by showing that an energy norm of the solution decreases after every update if the time step is less than a given critical value. The element eigenvalue inequality theorem is used to give the critical time step in terms of element eigenvalues.

Smolinski, P.↗