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 163 records · Page 9

PPINN: Parareal physics-informed neural network for time-dependent PDEs

Physics-informed neural networks (PINNs) encode physical conservation laws and prior physical knowledge into the neural networks, ensuring the correct physics is represented accurately while alleviating the need for supervised learning to a great degree. While effective for relatively short-term time integration, when long time integration of the time-dependent PDEs is sought, the time–space domain may become arbitrarily large and hence training of the neural network may become prohibitively expensive. To this end, we develop a parareal physics-informed neural network (PPINN), hence decomposing a long-time problem into many independent short-time problems supervised by an inexpensive/fast coarse-grained (CG) solver. In particular, the serial CG solver is designed to provide approximate predictions of the solution at discrete times, while initiate many fine PINNs simultaneously to correct the solution iteratively. There is a two-fold benefit from training PINNs with small-data sets rather than working on a large-data set directly, i.e., training of individual PINNs with small-data is much faster, while training the fine PINNs can be readily parallelized. Consequently, compared to the original PINN approach, the proposed PPINN approach may achieve a significant speed-up for long-time integration of PDEs, assuming that the CG solver is fast and can provide reasonable predictions of the solution, hence aiding the PPINN solution to converge in just a few iterations. To investigate the PPINN performance on solving time-dependent PDEs, we first apply the PPINN to solve the Burgers equation, and subsequently we apply the PPINN to solve a two-dimensional nonlinear diffusion–reaction equation. Furthermore, our results demonstrate that PPINNs converge in a few iterations with significant speed-ups proportional to the number of time-subdomains employed.

42 ENGINEERING↗

Progressive TDI Measurements with the PACE OCI ETU

The Plankton Aerosol Cloud ocean Ecosystem (PACE) Ocean Color Instrument (OCI) has completed the ground test program for its engineering unit (ETU) and testing of the flight unit will begin in the near future. OCIis a grating spectrometer with hyperspectral coverage from about 340 nm to 885 nm with 9 additional filtered channels in the SWIR. Two CCDs are used as detectors for the hyperspectral channels. One important operating mode of the CCDs on OCI is progressive time delay integration (or PTDI). In this mode, the charge in the CCD can be held for multiples of the nominal integration times. A series of these measurements can be made with progressively increasing multiples of the nominal integration time as the instrument scans across a uniform source. Ground testing with this operating mode on OCI ETU has shown promising results. This work will present measurements taken with the PTDI mode and the analysis of OCI ETU linearity and dynamic range

PACE↗

Optimization of the moderators in the STS preliminary design

This report details the results for an optimization of the dimensions of the moderators in the preliminary design of the Spallation Neutron Source Second Target Station (STS). This study uses the optimization algorithms of Dakota and an unstructured mesh model for the moderators in MCNP. More details on the unstructured mesh model and the automated mesh generation can be found in [3]. Parallel to this effort, the same moderator geometries have been optimized using a constructive solid geometry (CSG) MCNP model. More details on this model and its results can be found in [4]. Three optimal designs are selected for each moderator: one that is optimized for maximum peak brightness, one for maximum time-integrated brightness, and one for a combination of peak and time-integrated brightness. The backbone of the optimization work flow is provided by Dakota. For each set of design parameters requested by Dakota, a new solid geometry is automatically built in Creo and SpaceClaim, and subsequently exported to Attila4MC to generate an unstructured mesh geometry for MCNP. After the MCNP calculation is finished, the objective function (e.g., brightness metric) is returned to Dakota. After the new design has been evaluated, a result-file is written, and Dakota proposes the next set of design parameters to be evaluated. The loop continues until a specified convergence criterion has been met. The design parameters of the cylindrical (upper) moderator include the hydrogen radius, the premoderator thickness (top, bottom, radial), the beryllium radius and the horizontal position of the moderator. The crucial design choice is the hydrogen radius. A radius of 62 mm is shown to provide the maximum time-integrated brightness. The maximum peak brightness occurs with a radius of 40 mm. A combined (middle) design, which balances peak and time-integrated brightnesses, is obtained with a hydrogen radius of 50 mm. The premoderator thicknesses and the beryllium radius are slightly larger in the design optimized for time-integrated brightness than in the design optimized for peak brightness. The sensitivity to these two parameters is relatively small close to the optimal configurations. The hydrogen vessel and vacuum vessel wall thicknesses are dependent on the radius of the liquid hydrogen due to structural integrity requirements. The increased wall thicknesses for larger vessels significantly penalize the time-integrated brightness, with the maximum obtainable value reduced by more than 10% relative to earlier studies which used fixed vessel wall thicknesses. The impact of the variable wall thicknesses is much less for the peak brightness and combined brightness designs. The design parameters of the tube (lower) moderator selected for the optimization are the tube length, the annular premoderator thickness, the beryllium radius and the horizontal position of the moderator. The tube length is the crucial parameter and is chosen large (210 mm) and small (125 mm) in the designs optimized for time-integrated and peak brightness respectively. A combined optimal design has a tube length of 170 mm. The premoderator thickness and the beryllium radius are chosen larger in the design optimized for time-integrated brightness.

42 ENGINEERING↗

A comparative assessment of R. M. Young and tipping bucket rain gauges

Rain rates as derived from standard tipping bucket rain gauges have variable integration times corresponding to the interval between bucket tips. For example, the integration time for the Weathertronics rain gauge is given by delta(T) = 15.24/R (min), where R is the rain rate expressed in mm/h and delta(T) is the time between tips expressed in minutes. It is apparent that a rain rate of 1 mm/h has an integration time in excess of 15 minutes. Rain rates larger than 15.24 mm/h will have integration times smaller than 1 minute. The integration time is dictated by the time it takes to fill a small tipping bucket where each tip gives rise to 0.254 mm of rainfall. Hence, a uniform rain rate of 1 mm/h over a 15 minute period will give rise to the same rain rate as 0 mm/h rainfall over the first 14 minutes and 15 mm/h between 14 to 15 minutes from the reference tip. Hence, the rain intensity fluctuations may not be captured with the tipping bucket rain gauge for highly variable rates encompassing lower and higher values over a given integration time. The objective of this effort is to provide an assessment of the features of the R. M. Young capacitive gauge and to compare these features with those of the standard tipping bucket rain gauge. A number of rain rate-time series derived from measurements with approximately co-located gauges are examined.

Goldhirsh, Julius↗

National Ignition Facility Opacity Time Resolved Spectrometer Systems Engineering Final Project

The National Ignition Facility (NIF) is the world’s largest and most energetic laser facility. The NIF system is designed to produce high energy density (temperature and pressure) conditions through the application of its 192 laser beams. One of the users of NIF is the opacity platform developed to study the opacities at temperatures and densities relevant to the solar interior and stellar evolution. The platform was developed to study iron (Fe) opacity at temperatures relevant to the solar interior. The opacity campaign uses spectrometers to gather data. Spectrometers utilize crystals to produce x-ray spectra that are recorded on time-integrated and time-resolved detectors. The opacity spectrometer (OpSpec) currently fielded and in use at NIF uses a time integrated film channel to collect data. The opacity spectrometer time resolved (OpSpecTR) will utilize novel hCMOS detectors to capture time resolved images of spectra of interest. The key stakeholders identified for OpSpecTR included the physicists responsible for OpSpec and OpSpecTR, the Target Area Science and Engineering (TASE) department at NIF, the NIF and Photon Science (NIF & PS) Opacity program, the Nevada National Security Site (NNSS) Physics and Engineering program, the Sandia hCMOS manufacturing and testing program, and the Los Alamos National Laboratory (LANL) program sponsor. The Target and Experimental Operations (TEXOPS) was identified as a key stakeholder because the group includes the individuals that will physically interact with the OpSpecTR system as it participates in NIF experiments. The opacity platform collects data in a unique orientation relative to the existing diagnostics fielded at NIF. The existing infrastructure at NIF uses a diagnostic manipulator (DIM) to insert the diagnostic near the target chamber center to collect data during a NIF shot. Existing diagnostics collect data through the center line of the DIM axis and collect relevant data perpendicular to this axis. The opacity platform requires crystals mounted in a specific orientation which requires data collection parallel to the DIM axis. This deviation from standard NIF practices was a key factor in developing requirements.

42 ENGINEERING↗

Multigrid methods with space–time concurrency

Here, we consider the comparison of multigrid methods for parabolic partial differential equations that allow space–time concurrency. With current trends in computer architectures leading towards systems with more, but not faster, processors, space–time concurrency is crucial for speeding up time-integration simulations. In contrast, traditional time-integration techniques impose serious limitations on parallel performance due to the sequential nature of the time-stepping approach, allowing spatial concurrency only. This paper considers the three basic options of multigrid algorithms on space–time grids that allow parallelism in space and time: coarsening in space and time, semicoarsening in the spatial dimensions, and semicoarsening in the temporal dimension. We develop parallel software and performance models to study the three methods at scales of up to 16K cores and introduce an extension of one of them for handling multistep time integration. We then discuss advantages and disadvantages of the different approaches and their benefit compared to traditional space-parallel algorithms with sequential time stepping on modern architectures.

97 MATHEMATICS AND COMPUTING↗

Parallel multi-time step integration on a transputer system

An algorithm which allows different nodes of the finite element mesh to be integrated with different time steps is presented for second order finite element systems. The implementation of this algorithm on a system of transputer processors is discussed and a numerical example is used to evaluate the efficiency of the algorithm.

Smolinski, P.↗

Extreme and far ultraviolet astronomy from Voyagers 1 and 2

The instrumental characteristics, observational capabilities and scientific results of the Voyager 1 and 2 ultraviolet spectrometers are reviewed. These instruments provide current and ongoing access to low resolution spectra for a wide variety of astronomical sources in the 500 to 1700 A band. Observations of the brightest OB stars and hot subluminous stars as faint as V = 15 mag, are possible. In the EUV, at wavelengths shortward of 900 A, several new sources have been detected and a host of potential sources ruled out. In the far UV, particularly at wavelengths between 900 and 1200 A, Voyager is capable of observing a wide range of stellar and non-stellar sources. Such observations can often provide a valuable complement to IUE and other data sets at longer wavelengths. The Voyager spectrometers have proved remarkably stable photon counting instruments, capable of extremely long integration times. The long integration times, relatively large field of view, and location in the outer solar system also provide an ideal platform for observations of sources of faint diffuse emission, such as nebulae and the general sky background.

Holberg, J. B.↗

Comparative Analysis of Radiometer Systems Using Non-Stationary Processes

Radiometers require periodic calibration to correct for instabilities in the receiver response. Various calibration techniques exist that minimize the effect of instabilities in the receivers. The optimal technique depends upon many parameters. Some parameters are constrained by the particular application and others can be chosen in the system design. For example, the measurement uncertainty may be reduced to the limits of the resolution of the measurement (sensitivity) if periodic absolute calibration can be performed with sufficient frequency. However if the period between calibrations is long, a reference-differencing technique, i.e. Dicke-type design, can yield better performance. The measurement uncertainty not only depends upon the detection scheme but also on the number of pixels between calibrations, the integration time per pixel, integration time per calibration reference measurement, calibration reference temperature, and the brightness temperature of what is being measured. The best scheme for reducing the measurement uncertainty also depends, in large part, on the stability of the receiver electronics. In this presentation a framework for evaluating calibration schemes for a wide range of system architectures is presented. Two methods for treating receiver non-stationarity are compared with radiometer measurements.

Racette, Paul↗

Advanced acousto-optic signal processors

The basic acousto-optic signal processing architectures (spectrum analyzer, space-integrating, time-integrating, and triple product processor) systems and algorithms such as the chirp-Z transform are reviewed. New acousto-optic data processing systems and applications that utilze these basic architectures and new ones are described. These include a matched spatial filter acousto-optic processor, two new hybrid time and space-integrating systems, a triple product processor, and four new matrix-vector iterative feedback systems.

Casasent, D.↗

Scattering using real-time path integrals

Background: Path integrals are a powerful tool for solving problems in quantum theory that are not amenable to a treatment by perturbation theory. Most path integral computations require an analytic continuation to imaginary time. While imaginary time treatments of scattering are possible, imaginary time is not a natural framework for treating scattering problems. More importantly, quantum algorithms for calculating path integrals require real-time evolution. Purpose: Here, we test a recently introduced method for performing direct calculations of scattering observables using real-time path integrals in order to understand the challenges facing real-time path integral calculations of scattering observables. Method: The computations are based on a new interpretation of the path integral as the expectation value of a potential functional on cylinder sets of continuous paths with respect to a complex probability distribution. The method can in principle be applied to arbitrary short-range potentials. Results: The method is applied to compute matrix elements of Møller wave operators applied to narrow wave packets. These are used to calculate half-shell sharp-momentum transition matrix elements for one-dimensional potential scattering. The calculations for half-shell transition operator matrix elements converge to the numerical solution of the Lippmann-Schwinger equation. Conclusions: This work presents a proof in principle that scattering observables can be computed using real-time Feynman path integrals. While the computational method is not efficient, it can be improved. It provides a laboratory for studying quantum computational algorithms that are applicable to scattering problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Fully implicit crystal plasticity models representing orientations with modified Rodrigues parameters

Here, this work describes a crystal plasticity formulation combining several mathematical, numerical, and implementation choices to produce a highly efficient model. Specifically, the key choices in the implementation are (1) representing orientations with modified Rodrigues parameters, (2) implementing a fully coupled implicit time integration for the elastic stretch, the crystal orientations, and the model internal variables, (3) implementing the model in the NEML2 constitutive modeling framework, based on PyTorch, to vectorize the calculations and port the computation to GPUs and other hardware accelerators, and (4) an exact implementation of the consistent tangent matrix, even for arbitrary coupling to other field variables beyond the displacements, like temperature, neutron fluence, etc. The first two features of the model are, to our knowledge, novel. The paper considers each of these choices individually as well as the final model as a whole. This includes a full description of modified Rodrigues parameters, their advantages over other representations of orientations, the mathematical formulae and tools required to implement a model with modified Rodrigues parameters, and a detailed description of the geometry of the space of modified Rodrigues parameters (in an appendix). It also includes a description of a fully implicit time integration scheme for the orientations and the advantages in representing orientations with modified Rodrigues parameters in implementing such a model. The work then assess, via numerical examples, the advantages of fully coupled implicit time integration versus more common decoupled and explicit time integration schemes. These studies demonstrate the computational advantages of fully coupled integration versus other time integration algorithms, though the performance of the competing models depends on the complexity of the underlying single crystal model. The study concludes by demonstrating that the choice of time integration method affects the sharpness of the predicted texture, with explicit methods for integrating the orientations overestimating texture sharpness and implicit methods underestimating texture sharpness.

Crystal plasticity↗

Optimizing observations of drizzle onset with millimeter-wavelength radars

Cloud Doppler radars are increasingly used to study cloud and precipitation microphysical processes. Typical bulk cloud properties such as liquid or ice content are usually derived using the first three standard moments of the radar Doppler spectrum. Recent studies demonstrated the value of higher moments for the reduction of retrieval uncertainties and for providing additional insights into microphysical processes. Large effort has been undertaken, e.g., within the Atmospheric Radiation Measurement (ARM) program to ensure high quality of radar Doppler spectra. However, a systematic approach concerning the accuracy of higher moment estimates and sensitivity to basic radar system settings, such as spectral resolution, integration time and beam width, are still missing. Here In this study, we present an approach on how to optimize radar settings for radar Doppler spectra moments in the specific context of drizzle detection. The process of drizzle development has shown to be particularly sensitive to higher radar moments such as skewness. We collected radar raw data (I/Q time series) from consecutive zenith-pointing observations for two liquid cloud cases observed at the cloud observatory JOYCE in Germany. The I/Q data allowed us to process Doppler spectra and derive their moments using different spectral resolutions and integration times during identical time intervals. This enabled us to study the sensitivity of the spatiotemporal structure of the derived moments to the different radar settings. The observed signatures were further investigated using a radar Doppler forward model which allowed us to compare observed and simulated sensitivities and also to study the impact of additional hardware-dependent parameters such as antenna beam width. For the observed cloud with drizzle onset we found that longer integration times mainly modify spectral width (S w ) and skewness (S k ), leaving other moments mostly unaffected. An integration time of 2 s seems to be an optimal compromise: both observations and simulations revealed that a 10 s integration time – as it is widely used for European cloud radars – leads to a significant turbulence-induced increase of S w and reduction of S k compared to 2 s integration time. This can lead to significantly different microphysical interpretations with respect to drizzle water content and effective radius. A change from 2 s to even shorter integration times (0. 4 s) has much smaller effects on S w and S k . We also find that spectral resolution has a small impact on the moment estimations, and thus on the microphysical interpretation of the drizzle signal. Even the coarsest spectral resolution studied, 0. 08 ms -1 , seems to be appropriate for calculation moments of drizzling clouds. Moreover, simulations provided additional insight into the microphysical interpretation of the skewness signatures observed: in low (high)-turbulence conditions, only drizzle larger than 20 µm (40 µm) can generate S k values above the S k noise level (in our case 0.4). Higher S k values are also obtained in simulations when smaller beam widths are adopted.

54 ENVIRONMENTAL SCIENCES↗

Time-dependent integral transport in one-dimensional infinite media using dimensionless variables and the reduced collision formulation

A new method for calculating time-dependent integral neutron transport single collision fluxes has been developed for one-dimensional infinite media problems using dimensionless variables and the reduced collision formulation. By transforming the integral scalar flux equation into dimensionless variables, one may immediately integrate over the time domain and reduce the problem complexity to integrations strictly over the spatial domain, significantly reducing the computational work required to obtain a solution. The results have been compared to published benchmark solutions for validation. Further, reformulating the scalar flux as a Green’s function has allowed for the determination of an analytical expression for the relative error of any collision level scalar flux to the summed scalar flux. This provides the absolute number of collisions required for convergence at any given mean free time, using previous work in one-dimensional Cartesian transport and derived work herein for the one-dimensional spherical transport.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

On a new algorithm for time step integration of nonlinear systems

A new implicit algorithm for time step integration of finite element structural dynamic equations is presented. Convergence, stability and numerical damping properties are discussed. Due to the way nonlinear structural behavior is taken into account, the algorithm is expected to compare favorably with existing ones. Some simple numerical results are presented. A related explicit algorithm is also derived and shortly discussed.

Anderheggen, E.↗

The Complete Spectral Catalog of Bright BATSE Gamma-Ray Bursts

We present a systematic spectral analysis of 350 bright Gamma-Ray Bursts (GRBs) observed by the Burst and Transient Source Experiment (BATSE; approx. 30 keV - 2 MeV; including 17 short GRBs) with high energy and time resolution. Our sample was selected from the complete set of 2704 BATSE GRBs based on their energy fluence or peak photon flux values to assure good statistics. To obtain well-constrained, model-unbiased spectral parameters, a set of various photon models is used to fit each spectrum, and internal characteristics of each model are also investigated. A thorough analysis has been performed on 342 time-integrated and 8459 time-resolved burst spectra, and the effects of integration times in determining the spectral parameters are explored. The analysis results presented here provide the most detailed perspective of spectral aspects of the GRB prompt emission to date. Using the results, we study correlations among spectral parameters and spectral evolutions. The results of all spectral fits are available electronically in FITS format, from the High-Energy Astrophysics Science Archive Research Center (HEASARC).

Kaneko, Yuki↗

Pulmonary venous flow determinants of left atrial pressure under different loading conditions in a chronic animal model with mitral regurgitation

BACKGROUND: The aim of our study was to quantitatively compare the changes and correlations between pulmonary venous flow variables and mean left atrial pressure (mLAP) under different loading conditions in animals with chronic mitral regurgitation (MR) and without MR. METHODS: A total of 85 hemodynamic conditions were studied in 22 sheep, 12 without MR as control (NO-MR group) and 10 with MR (MR group). We obtained pulmonary venous flow systolic velocity (Sv) and diastolic velocity (Dv), Sv and Dv time integrals, their ratios (Sv/Dv and Sv/Dv time integral), mLAP, left ventricular end-diastolic pressure, and MR stroke volume. We also measured left atrial a, x, v, and y pressures and calculated the difference between v and y pressures. RESULTS: Average MR stroke volume was 10.6 +/- 4.3 mL/beat. There were good correlations between Sv (r = -0.64 and r = -0.59, P <.01), Sv/Dv (r = -0.62 and r = -0.74, P <.01), and mLAP in the MR and NO-MR groups, respectively. Correlations were also observed between Dv time integral (r = 0.61 and r = 0.57, P <.01) and left ventricular end-diastolic pressure in the MR and NO-MR groups. In velocity variables, Sv (r = -0.79, P <.001) was the best predictor of mLAP in both groups. The sensitivity and specificity of Sv = 0 in predicting mLAP 15 mm Hg or greater were 86% and 85%, respectively. CONCLUSION: Pulmonary venous flow variables correlated well with mLAP under altered loading conditions in the MR and NO-MR groups. They may be applied clinically as substitutes for invasively acquired indexes of mLAP to assess left atrial and left ventricular functional status.

Pulmonary Veins/physiopathology↗