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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

Hybrid-Vlasov Simulation of Soft X-Ray Emissions at the Earth’S Dayside Magnetospheric Boundaries

Solar wind charge exchange produces emissions in the soft X-ray energy range which can enable the study of near-Earth space regions such as the magnetopause, the magnetosheath and the polar cusps by remote sensing techniques. The Solar wind–Magnetosphere–Ionosphere Link Explorer (SMILE) and Lunar Environment heliospheric X-ray Imager (LEXI) missions aim to obtain soft X-ray images of near-Earth space thanks to their Soft X-ray Imager (SXI) instruments. While earlier modeling works have already simulated soft X-ray images as might be obtained by SMILE SXI during its mission, the numerical models used so far are all based on the magnetohydrodynamics description of the space plasma. To investigate the possible signatures of ion-kinetic-scale processes in soft X-ray images, we use for the first time a global hybrid-Vlasov simulation of the geospace from the Vlasiator model. The simulation is driven by fast and tenuous solar wind conditions and purely southward interplanetary magnetic field. We first produce global X-ray images of the dayside near-Earth space by placing a virtual imaging satellite at two different locations, providing meridional and equatorial views. We then analyze regional features present in the images and show that they correspond to signatures in soft X-ray emissions of mirror mode wave structures in the magnetosheath and flux transfer events (FTEs) at the magnetopause. Our results suggest that, although the time scales associated with the motion of those transient phenomena will likely be significantly smaller than the integration time of the SMILE and LEXI imagers, mirror-mode structures and FTEs can cumulatively produce detectable signatures in the soft X-ray images. For instance, a local increase by 30% in the proton density at the dayside magnetopause resulting from the transit of multiple FTEs leads to a 12% enhancement in the line-of-sight- and time integrated soft X-ray emissivity originating from this region. Likewise, a proton density increase by 14% in the magnetosheath associated with mirror-mode structures can result in an enhancement in the soft X-ray signal by 4%. These are likely conservative estimates, given that the solar wind conditions used in the Vlasiator run can be expected to generate weaker soft X-ray emissions than the more common denser solar wind. These results will contribute to the preparatory work for the SMILE and LEXI missions by providing the community with quantitative estimates of the effects of small-scale, transient phenomena occurring on the dayside.

magnetosphere↗
DOE OSTI2024

Adaptive tau-leaping methods for microscopic-lattice kinetic Monte Carlo simulations

Traditional Kinetic Monte Carlo (KMC) approaches, rooted in Gillespie’s stochastic simulation algorithm, become computationally demanding in systems with a large range of timescales. The goal of this work is to propose and study new adaptive lattice-KMC time integration strategies for spatially non-uniform systems. To that end, two novel adaptive tau-leaping methods and their corresponding time integration strategies are developed based on the idea of the “n-fold” direct KMC method. These strategies allow for the simultaneous execution of multiple reactions, advancing time by adaptively selected coarse increments. We present numerical experiments comparing the proposed methods with existing approaches in a catalytic surface kinetics application involving ammonia decomposition.

Bimolecular reactions↗
DOE OSTI2014

Long-time uncertainty propagation using generalized polynomial chaos and flow map composition

We present an efficient and accurate method for long-time uncertainty propagation in dynamical systems. Uncertain initial conditions and parameters are both addressed. The method approximates the intermediate short-time flow maps by spectral polynomial bases, as in the generalized polynomial chaos (gPC) method, and uses flow map composition to construct the long-time flow map. In contrast to the gPC method, this approach has spectral error convergence for both short and long integration times. The short-time flow map is characterized by small stretching and folding of the associated trajectories and hence can be well represented by a relatively low-degree basis. The composition of these low-degree polynomial bases then accurately describes the uncertainty behavior for long integration times. The key to the method is that the degree of the resulting polynomial approximation increases exponentially in the number of time intervals, while the number of polynomial coefficients either remains constant (for an autonomous system) or increases linearly in the number of time intervals (for a non-autonomous system). The findings are illustrated on several numerical examples including a nonlinear ordinary differential equation (ODE) with an uncertain initial condition, a linear ODE with an uncertain model parameter, and a two-dimensional, non-autonomous double gyre flow.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1982

Development of mixed time partition procedures for thermal analysis of structures

The computational methods used to predict and optimize the thermal-structural behavior of aerospace vehicle structures are reviewed. In general, two classes of algorithms, implicit and explicit, are used in transient thermal analysis of structures. Each of these two methods has its own merits. Due to the different time scales of the mechanical and thermal responses, the selection of a time integration method can be a difficult yet critical factor in the efficient solution of such problems. Therefore mixed time integration methods for transient thermal analysis of structures are being developed. This proposed methodology would be readily adaptable to existing computer programs for structural thermal analysis.

Liu, W. K.↗
DOE OSTI2021

On the computational solution of vector-density based continuum dislocation dynamics models: A comparison of two plastic distortion and stress update algorithms

Continuum dislocation dynamics models of mesoscale plasticity consist of dislocation transport-reaction equations coupled with crystal mechanics equations. The coupling between these two sets of equations is such that dislocation transport gives rise to the evolution of plastic distortion (strain), while the evolution of the latter fixes the stress from which the dislocation velocity field is found via a mobility law. Earlier solutions of these equations employed a staggered solution scheme for the two sets of equations in which the plastic distortion was updated via time integration of its rate, as found from Orowan's law. In this work, we show that such a direct time integration scheme can suffer from accumulation of numerical errors. We introduce an alternative scheme based on field dislocation mechanics that ensures consistency between the plastic distortion and the dislocation content in the crystal. The new scheme is based on calculating the compatible and incompatible parts of the plastic distortion separately, and the incompatible part is calculated from the current dislocation density field. Stress field and dislocation transport calculations were implemented within a finite element based discretization of the governing equations, with the crystal mechanics part solved by a conventional Galerkin method and the dislocation transport equations by the least squares method. A simple test was first performed to show the accuracy of the two schemes for updating the plastic distortion, which shows that the solution method based on field dislocation mechanics is more accurate. This method then was used to simulate an austenitic steel crystal under uniaxial loading and multiple slip conditions. By considering dislocation interactions caused by junctions, a hardening rate similar to discrete dislocation dynamics simulation results was obtained. Finally, the simulations show that dislocations exhibit some self-organized structures as the strain is increased.

36 MATERIALS SCIENCE↗
NASA NTRS1986

Star of Condor - A strontium critical velocity experiment, Peru, 1983

'Star of Condor' was a critical velocity experiment using Sr vapor produced in a radial shaped charge, which was carried to 571.11 km altitude on a Taurus-Tomahawk rocket launched from Punto Lobos, Peru, and detonated in the plane of the magnetic field lines so that all ranges of pitch angles from parallel to B to perpendicular to B were covered. Sr has a critical velocity of 3.3 km/s, and from observation, 42.5 percent of the neutral Sr gas had a velocity component perpendicular to B exceeding that value. No Sr ion emissions were detected shortly after the burst with usual TV integration times. However, about 10 min after the detonation a faint field-aligned streak was discovered with long TV integration times. The brightness is estimated as 5 R, which, combined with the streak geometry, implies an ion production of 2.4 x 10 to the 19th ions. This is only 0.0036 percent ionization of the Sr vapor. All the ions could easily have been produced by thermal ionization from the original detonation thermal distribution. The breakup of the Sr gas into small bloblike structures may have allowed the high-energy electrons to escape before an ionization cascade could be produced. For whatever reason, the Alfven mechanism proposed for space plasmas in the absence of laboratory walls did not produce an ionization cascade in the experiment.

Wescott, E. M.↗
DOE OSTI2025

Axisymmetric gyrokinetic simulation of ASDEX-Upgrade scrape-off layer using a conservative implicit BGK collision operator

Collisions play an important role in turbulence and transport of fusion plasmas. For kinetic simulations, as the collisionality increases in the domain of interest, the size of the time step to resolve the collisional physics can become overly restrictive in an explicit time integration scheme, leading to high computational cost. With the aim of overcoming such restriction, we have implemented an implicit Bhatnagar–Gross–Krook (BGK) collision operator for use in the discontinuous Galerkin full-f gyrokinetic solver within the Gkeyll framework, which, when combined with Gkeyll's traditional explicit time integrator for collisionless advection, can significantly increase the time step in gyrokinetic simulations of highly collisional regimes. To ensure conservation of density, momentum, and energy, we utilize an iterative scheme to correct the discretized approximation to the equilibrium Maxwellian distribution to which the BGK collision operator relaxes. We have further generalized the BGK infrastructure, both the implicit scheme and the correction routine, to handle cross-species collisions. This improved implicit and conservative BGK operator is benchmarked against the more accurate but more computationally expensive Lenard–Bernstein–Dougherty (LBD) operator, which has been utilized in prior studies with Gkeyll. The implicit BGK operator enables 2D axisymmetric simulations of the ASDEX-Upgrade scrape-off layer to run 56 times faster to completion than the simulations with the LBD operator, because the BGK operator is more robust and converges at a lower resolution than is required by the LBD operator. Additionally, in this more collisional limit, we demonstrate that the results of our simulations utilizing the implicit BGK operator agreed well with simulations utilizing the more computationally expensive LBD operator.

Gyrokinetic simulations↗
DOE OSTI2024

Collocation methods for nonlinear differential equations on low-rank manifolds

We introduce new methods for integrating nonlinear differential equations on low-rank manifolds. These methods rely on interpolatory projections onto the tangent space, enabling low-rank time integration of vector fields that can be evaluated entry-wise. A key advantage of our approach is that it does not require the vector field to exhibit low-rank structure, thereby overcoming significant limitations of traditional dynamical low-rank methods based on orthogonal projection. To construct the interpolatory projectors, we develop a sparse tensor sampling algorithm based on the discrete empirical interpolation method (DEIM) that parameterizes tensor train manifolds and their tangent spaces with cross interpolation. Using these projectors, we propose two time integration schemes on low-rank tensor train manifolds. The first scheme integrates the solution at selected interpolation indices and constructs the solution with cross interpolation. The second scheme generalizes the well-known orthogonal projector-splitting integrator to interpolatory projectors. We demonstrate the proposed methods with applications to several tensor differential equations arising from the discretization of partial differential equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

A High-Order L-Band HTS Filter for Sensitive Detecting

In this paper, a new high-temperature superconducting (HTS)-based microwave radiometer with an improved sensitivity is presented. The cryogenic receiver front end consists of an HTS filter and a cryogenic low noise amplifier (LNA). The cryogenic receiver front end shows an ultra-low noise figure and can suppress radio frequency interference (RFI) effectively. The proposed HTS filter works at a center frequency of 1.4135 GHz with a bandwidth of 25 MHz. The measured mid-band insertion loss, side band steepness, and out-of-band attenuation of the HTS filter are 0.14 dB, 35 dB/MHz, and 80 dB, respectively. The noise figure of the cryogenic LNA is about 0.27 dB at a temperature of 77 K. Compared with other total power radiometers, the proposed radiometer has a lower receiver noise temperature, which can improve the sensitivity with a short integration time of the satellite-based salinity meter. In addition, since the bandwidth of the salinity meter is fixed and the integration time of satellite-based equipment is limited, such a low receiver noise temperature can increase the flexibility of future satellite payload configuration program.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Parallel-In-Time Magnus Integrators

Magnus integrators are a subset of geometric integration methods for the numerical solution of ordinary differential equations that conserve certain invariants in the numerical solution. Our report explores temporal parallelism of Magnus integrators, particularly in the context of nonlinear problems. The method combines the concurrent computation of matrix commutators and exponentials within a time step with a pipelined iteration applied to multiple time steps in parallel. The accuracy and efficiency of time parallel Magnus methods up to order six are highlighted through numerical examples and demonstrate that critical parallel speedup is possible compared to serial methods.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1991

Application of Krylov exponential propagation to fluid dynamics equations

This paper presents an application of matrix exponentiation via Krylov subspace projection, to the solution of fluid dynamics problems. The main idea is to approximate the operation exp(A)v by means of a projection-like process onto a Krylov subspace. This results in a computation of an exponential matrix vector product similar to the one above but of a much smaller size. Time integration schemes can then be devised to exploit this basic computational kernel. The motivation of this approach is to provide time-integration schemes that are essentially of an explicit nature but which have good stability properties.

Saad, Y.↗
NASA NTRS1991

Application of Krylov exponential propagation to fluid dynamics equations

An application of matrix exponentiation via Krylov subspace projection to the solution of fluid dynamics problems is presented. The main idea is to approximate the operation exp(A)v by means of a projection-like process onto a krylov subspace. This results in a computation of an exponential matrix vector product similar to the one above but of a much smaller size. Time integration schemes can then be devised to exploit this basic computational kernel. The motivation of this approach is to provide time-integration schemes that are essentially of an explicit nature but which have good stability properties.

Saad, Youcef↗
DOE OSTI2014

High repetition rate laser-induced breakdown spectroscopy using acousto-optically gated detection

This contribution introduces a new type of setup for fast sample analysis using laser-induced breakdown spectroscopy (LIBS). The novel design combines a high repetition rate laser (up to 50 kHz) as excitation source and an acousto-optical modulator (AOM) as a fast switch for temporally gating the detection of the emitted light. The plasma radiation is led through the active medium of the AOM where it is diffracted on the transient ultrasonic Bragg grid. The diffracted radiation is detected by a compact Czerny-Turner spectrometer equipped with a CCD line detector. Utilizing the new combination of high repetition rate lasers and AOM gated detection, rapid measurements with total integration times of only 10 ms resulted in a limit of detection (LOD) of 0.13 wt.% for magnesium in aluminum alloys. This short integration time corresponds to 100 analyses/s. Temporal gating of LIP radiation results in improved LODs and consecutively higher sensitivity of the LIBS setup. Therefore, an AOM could be beneficially utilized to temporally detect plasmas induced by high repetition rate lasers. The AOM in combination with miniaturized Czerny-Turner spectrometers equipped with CCD line detectors and small footprint diode pumped solid state lasers results in temporally gateable compact LIBS setups.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Long duration x-ray source development for x-ray diffraction at the National Ignition Facility

We present the results of experiments to produce a 10 ns-long, quasi-monochromatic x-ray source. This effort is needed to support time-resolved x-ray diffraction (XRDt) measurements of phase transitions during laser-driven dynamic compression experiments at the National Ignition Facility. To record XRDt of phase transitions as they occur, we use high-speed (∼1 ns) gated hybrid CMOS detectors, which record multiple frames of data over a timescale of a few to tens of ns. Consequently, to make effective use of these imagers, XRDt needs the x-ray source to be narrow in energy and uniform in time as long as the sensors are active. The x-ray source is produced by a laser irradiated Ge foil. Our results indicate that the x-ray source lasts during the whole duration of the main laser pulse. Both time-resolved and time-integrated spectral data indicate that the line emission is dominated by the He-α complex over higher energy emission lines. Time-integrated spectra agree well with a one-dimensional Cartesian simulation using HYDRA that predicts a conversion efficiency of 0.56% when the incident intensity is 2 × 1015 W/cm2 on a Ge backlighter.

Werellapatha, K. (ORCID:0000000293154054)↗
DOE OSTI2017

Computational singular perturbation analysis of stochastic chemical systems with stiffness

Computational singular perturbation (CSP) is a useful method for analysis, reduction, and time integration of stiff ordinary differential equation systems. It has found dominant utility, in particular, in chemical reaction systems with a large range of time scales at continuum and deterministic level. On the other hand, CSP is not directly applicable to chemical reaction systems at micro or meso-scale, where stochasticity plays an non-negligible role and thus has to be taken into account. In this work we develop a novel stochastic computational singular perturbation (SCSP) analysis and time integration framework, and associated algorithm, that can be used to not only construct accurately and efficiently the numerical solutions to stiff stochastic chemical reaction systems, but also analyze the dynamics of the reduced stochastic reaction systems. Furthermore, the algorithm is illustrated by an application to a benchmark stochastic differential equation model, and numerical experiments are carried out to demonstrate the effectiveness of the construction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Temporal Characterization of Hydrates System Dynamics beneath Seafloor Mounds. Integrating Time-Lapse Electrical Resistivity Methods and In Situ Observations of Multiple Oceanographic Parameters

This study was designed to investigate temporal variations in hydrate system dynamics by measuring changes in volumes of hydrate beneath hydrate-bearing mounds on the continental slope of the northern Gulf of Mexico, the landward extreme of hydrate occurrence in this region. Direct Current Resistivity (DCR) measurements were made contemporaneously with measurements of oceanographic parameters at Woolsey Mound, a carbonate-hydrate complex on the mid-continental slope, where formation and dissociation of hydrates are most vulnerable to variations in oceanographic parameters affected by climate change, and where changes in hydrate stability can readily translate to loss of seafloor stability, impacts to benthic ecosystems, and venting of greenhouse gases to the water-column, and eventually, the atmosphere. We focused our study on hydrate within seafloor mounds because the structurally-focused methane flux at these sites likely causes hydrate formation and dissociation processes to occur at higher rates than at sites where the methane flux is less concentrated and we wanted to maximize our chances of witnessing association/dissociation of hydrates. We selected a particularly well-studied hydrate-bearing seafloor mound near the landward extent of the hydrate stability zone, Woolsey Mound (MC118). This mid-slope site has been studied extensively and the project was able to leverage considerable resources from the team’s research experience at MC118. The site exhibits seafloor features associated with gas expulsion, hydrates have been documented at the seafloor, and changes in the outcropping hydrates have been documented, photographically, to have occurred over a period of months. We conducted observatory-based, in situ measurements to 1) characterize, geophysically, the sub-bottom distribution of hydrate and its temporal variability, and 2) contemporaneously record relevant environmental parameters (temperature, pressure, salinity, turbidity, bottom currents) to detect short-term changes within the hydrates system, identify relationships/impacts of local oceanographic parameters on the hydrates system, and improve our understanding of how seafloor instability is affected by hydrates-driven changes. A 2009 DCR survey of MC118 demonstrated that we could image resistivity anomalies to a depth of 75m below the seafloor in water depths of 1km. We reconfigured this system to operate autonomously on the seafloor in a pre-programmed mode, for periods of months. We designed and built a novel seafloor lander and deployment capability that would allow us to investigate the seafloor at potential deployment sites and deploy instruments only when conditions met our criteria. This lander held the DCR system, controlling computers, and battery power supply, as well as instruments to record oceanographic parameters. During the first of two cruises to the study site, we conducted resistivity surveying, selected a monitoring site, and deployed the instrumented lander and DCR, centered on what appeared to be the most active locations within the site, programmed to collect a DCR profile, weekly. After a 4.5-month residence on the seafloor, the team recovered all equipment. Unfortunately, several equipment failures occurred prior to recovery of the instrument packages. Prior to the failures, however, two resistivity profiles were collected together with oceanographic data. Results show, unequivocally, that significant changes can occur in both hydrate volume and distribution during time periods as brief as one week. Occurrences appear to be controlled by both deep and near-surface structure. Results have been integrated with seismic data from the area and show correspondence in space of hydrate and structures, including faults and gas chimneys.

58 GEOSCIENCES↗
DOE OSTI2015

Emission mechanisms in Al-rich AlGaN/AlN quantum wells assessed by excitation power dependent photoluminescence spectroscopy

The optical properties of Al-rich AlGaN/AlN quantum wells are assessed by excitation-power-dependent time-integrated (TI) and time-resolved (TR) photoluminescence (PL) measurements. Two excitation sources, an optical parametric oscillator and the 4th harmonics of a Ti:sapphire laser, realize a wide range of excited carrier densities between 10{sup 12} and 10{sup 21 }cm{sup −3}. The emission mechanisms change from an exciton to an electron-hole plasma as the excitation power increases. Accordingly, the PL decay time is drastically reduced, and the integrated PL intensities increase in the following order: linearly, super-linearly, linearly again, and sub-linearly. The observed results are well accounted for by rate equations that consider the saturation effect of non-radiative recombination processes. Using both TIPL and TRPL measurements allows the density of non-radiative recombination centers, the internal quantum efficiency, and the radiative recombination coefficient to be reliably extracted.

77 NANOSCIENCE AND NANOTECHNOLOGY↗