Modeling fast diffusion processes in time integration of stiff stochastic differential equations.
Abstract not provided.
Han, Xiaoying↗
SEARCH · Engineering Papers
Results for “time integrator”
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.
Leveraging Mixed Precision in Exponential Time Integration Methods
Explore the source record for details and available documents.
Mathematics and Computing↗
Design and characterization of the time-resolved opacity spectrometer (OpSpecTR) for the NIF iron opacity campaign
A new time-resolved opacity spectrometer (OpSpecTR) is currently under development for the National Ignition Facility (NIF) opacity campaign. Here, the spectrometer utilizes Icarus version 2 (IV2) hybridized complementary metal–oxide–semiconductor sensors to collect gated data at the time of the opacity transmission signal, unlocking the ability to collect higher-temperature measurements on NIF. Experimental conditions to achieve higher temperatures are feasible; however, backgrounds will dominate the data collected by the current time-integrating opacity spectrometer. The shortest available OpSpecTR integration time of ~2 ns is predicted to reduce self-emission and other late-time backgrounds by up to 80%. Initially, three Icarus sensors will be used to collect data in the self-emission, backlighter, and absorption regions of the transmission spectrum, with plans to upgrade to five Daedalus sensors in future implementations with integration times of ~1.3 ns. We present the details of the diagnostic design along with recent characterization results of the IV2 sensors.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Analysis of a Real-Time Separation Assurance System with Integrated Time-in-Trail Spacing
This paper describes the implementation and analysis of an integrated ground-based separation assurance and time-based metering prototype system into the Center-TRACON Automation System. The integration of this new capability accommodates constraints in four-dimensions: position (x-y), altitude, and meter-fix crossing time. Experiments were conducted to evaluate the performance of the integrated system and its ability to handle traffic levels up to twice that of today. Results suggest that the integrated system reduces the number and magnitude of time-in-trail spacing violations. This benefit was achieved without adversely affecting the resolution success rate of the system. Also, the data suggest that the integrated system is relatively insensitive to an increase in traffic of twice the current levels.
Aweiss, Arwa S.↗
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↗
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↗
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↗