Engineering Papers⌕ Search

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.

At least 181 records · Page 10

Three-Dimensional Simulations of Electron Beams Focused by Periodic Permanent Magnets

A fully three-dimensional (3D) model of an electron beam focused by a periodic permanent magnet (PPM) stack has been developed. First, the simulation code MAFIA was used to model a PPM stack using the magnetostatic solver. The exact geometry of the magnetic focusing structure was modeled; thus, no approximations were made regarding the off-axis fields. The fields from the static solver were loaded into the 3D particle-in-cell (PIC) solver of MAFIA where fully 3D behavior of the beam was simulated in the magnetic focusing field. The PIC solver computes the time-integration of electromagnetic fields simultaneously with the time integration of the equations of motion of charged particles that move under the influence of those fields. Fields caused by those moving charges are also taken into account; thus, effects like space charge and magnetic forces between particles are fully simulated. The electron beam is simulated by a number of macro-particles. These macro-particles represent a given charge Q amounting to that of several million electrons in order to conserve computational time and memory. Particle motion is unrestricted, so particle trajectories can cross paths and move in three dimensions under the influence of 3D electric and magnetic fields. Correspondingly, there is no limit on the initial current density distribution of the electron beam, nor its density distribution at any time during the simulation. Simulation results including beam current density, percent ripple and percent transmission will be presented, and the effects current, magnetic focusing strength and thermal velocities have on beam behavior will be demonstrated using 3D movies showing the evolution of beam characteristics in time and space. Unlike typical beam optics models, this 3D model allows simulation of asymmetric designs such as non- circularly symmetric electrostatic or magnetic focusing as well as the inclusion of input/output couplers.

Kory, Carol L.↗

Response of joint dominated space structures

An approximate method is developed for estimating the transient response of nonlinear systems in terms of linearized modes of response. Its advantages are that it is computationally more efficient than the time integration method and that it is possible to view the design problem in the more traditional physical terms of modal response. The major drawback of the approximate method is loss of accuracy. It seems that both approximate methods and time integration have their roles in design. Approximate methods provide efficient tools for performing parametric studies and they supply physical insights into how to optimize system performance that are not easily inferred from strictly numerical methods. Time integration provides a method for assessing the accuracy of the approximate solution for key simulations and for fine tuning the final design. In the procedure presented the nonlinear system is approximated by an equivalent linear system in which the system parameters are constant over the range of transient response.

Source record↗

The Miniaturized Moessbauer Spectrometers MIMOS II on MER: Four Years of Operation - A Summary

The two Miniaturized Moessbauer Spectrometers (MIMOS II) on board the two Mars Exploration Rovers Spirit and Opportunity have now been collecting important scientific data for more than four years. The spectrometers provide information about Fe-bearing mineral phases and determine Fe oxidation states. The total amount of targets analized exceeds 600, the total integration time exceeds 260 days for both rovers. Since landing, more than five half-lives of the Co-57 MB sources have past (intensity at the time of landing approx. 150 mCi). Current integration times are about 50 hours in order to achieve reasonable statistics as opposed to 8 hours at the beginning of the mission. In total, 13 different mineral phases were detected: Olivine, pyroxene, hematite, magnetite and nanophase ferric oxide were detected at both landing sites. At Gusev, ilmenite, goethite, a ferric sulfate phase and a yet unassigned phase (in the rock Fuzzy Smith) were detected. At Meridiani, jarosite, metallic iron in meteoritic samples (kamacite), troilite, and an unassigned ferric phase were detected. Jarosite and goethite are of special interest, as these minerals are indicators for water activity. In this abstract, an overview of Moessbauer results will be given, with a focus on data obtained since the last martian winter. The MER mission has proven that Moessbauer spectroscopy is a valuable tool for the in situ exploration of extraterrestrial bodies and for the study of Febearing samples. The experience gained through the MER mission makes MIMOS II a obvious choice for future missions to Mars and other targets. Currently, MIMOS II is on the scientific payload of two approved future missions: Phobos Grunt (Russian Space Agency; 2009) and ExoMars (European Space Agency; 2013).

Fleischer, I.↗

Connecting Femtosecond Transient Absorption Microscopy with Spatially Coregistered Time Averaged Optical Imaging Modalities

Multimodal all-optical imaging involving coregistered femtosecond transient absorption microscopy (TAM), time-integrated photoluminescence (PL), and steady-state modalities such as confocal reflectance and transmission offers an appealing approach to gain a comprehensive understanding of complex electronic excited-state phenomena in spatially heterogeneous systems. A unique combination of these modalities allows us to unravel not only the competing electronic excited-state dynamical processes but also the underlying morphological information with simultaneous high temporal and spatial resolution. However, correlating the various images obtained from time-resolved and time-independent modalities is generally nontrivial and particularly challenging when the electronic dynamics under study evolve in both time and space. Here, we demonstrate a new approach for rationally correlating time-resolved microscopy with coregistered time-integrated or steady-state modalities. Specifically, our approach involves an extended global lifetime analysis of the time-resolved microscopic data set to separate distinct dynamical processes taking place on commensurate time scales, and the resulting decay-associated amplitude maps (DAAMs) were applied to explore correlations with the images acquired using time-independent modalities. The feasibility of our approach was validated through analyzing a multimodal data set acquired from a thin film of chloride-containing mixed lead halide perovskites (CH 3 NH 3 PbI 3–x Cl x ) using femtosecond transient absorption, time-integrated PL, and confocal reflectance microscopies. Analysis of the results obtained enable us to gain new insight into the complex ultrafast relaxation dynamics in this highly heterogeneous system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fast Multigrid Reduction-in-Time for Advection via Modified Semi-Lagrangian Coarse-Grid Operators

Many iterative parallel-in-time algorithms have been shown to be highly efficient for diffusion-dominated partial differential equations (PDEs) but are inefficient or even divergent when applied to advection-dominated PDEs. We consider the application of the multigrid reduction-in-time (MGRIT) algorithm to linear advection PDEs. Here, the key to efficient time integration with this method is using a coarse-grid operator that provides a sufficiently accurate approximation to the so-called ideal coarse-grid operator. For certain classes of semi-Lagrangian discretizations, we present a novel semi-Lagrangian-based coarse-grid operator that leads to fast and scalable multilevel time integration of linear advection PDEs. The coarse-grid operator is composed of a semi-Lagrangian discretization followed by a correction term, with the correction designed so that the leading-order truncation error of the composite operator is approximately equal to that of the ideal coarse-grid operator. Parallel results show substantial speed-ups over sequential time integration for variable-wave-speed advection problems in one and two spatial dimensions, and using high-order discretizations up to order five. The proposed approach establishes the first practical method that provides small and scalable MGRIT iteration counts for advection problems.

97 MATHEMATICS AND COMPUTING↗

A correlative study of magnetic flux transfer in the magnetosphere

The applied magnetic flux of the southward component of the IMF in GSM coordinates (Bz-) which impinges upon the sunward magnetopause is compared to the time integral of the auroral AL index during 56 intervals within a 3-month period in 1969 when interplanetary records from Heos 1 and Explorer 35 were available. The periods of magnetic activity were those for which Bz was greater than 0 and AL was less than 20gamma at the beginning and end of the interval. It was found that for these intervals, the time integral of the AL index was proportional to the applied magnetic flux with a correlation coefficient of 0.94. In addition, the empirical relationships between magnetic flux transfer, applied southward flux, and the time integral of AL arrived at by Holzer and Slavin (1978) on the basis of expansion and contractions of the forward magnetosphere observed with OGO 5 are reexamined and improved.

Holzer, R. E.↗

Low Frequency 1/f Noise on QWIPs, nBn, and Superlattice Focal Plane Array

Noise Equivalent Difference Temperature (NE△T) is a standard performance metric for most infrared focal plane array (FPA) systems. The frequency bandwidth and range associated with NE△T is at high frequency and it does not describe the long time noise behavior or the very low frequency noise of the FPA. Very low frequency noise measurement requires data capture that takes longer time duration and sampling interval. This study investigates low frequency noise in QWIP, nBn and LWIR n-type Complementary Barrier Infrared Detector (CBIRD) FPAs. The corner frequencies are extracted from the power spectral density (PSD) as function of frequency. The peak wavelength, quantum efficiency of QWIP detector are 6.2 µm and 2 %, respectively. QWIP FPA has a mean NE△T ~ 25 mK at an operating temperature of 65 K and an integration time of 16 msec. The mean QWIP PSD plot shows a corner frequency of < 0.5 mHz. The nBn FPA with 4 µm cut off and quantum efficiency of 67 % has NE△T ~ 15.6 mK at an integration time ~ 7.52 msec and an operating temperature of 120 K. The nBn has corner frequency of > 50 mHz. Lastly the n-type CBIRD FPA with a 50% cutoff at 8.8 µm and quantum efficiency of ~ 50 % has an NE△T ~ 18.6 mK at an integration time of 1.86 msec and operating temperature 120 K. The superlattice FPA has a corner frequency ~ 10 mHz. The investigation of the tail on the NE△T histogram reveals that its origin is not accounted for entirely by the high noise current, but also it needs the inclusion of lower responsivity for some pixels.

Pepper, B.↗

Instruction manual, optical effects module electronic controller and processor, model OEMCP

The OEM-1 electronic module is discussed; it is comprised of four subsystems: the signal processing and display; the stepper motor controls; the chopper controls; and the dc-dc invertor. The OEM-1 module controls the sample wheel so that the relative transmittance of the samples can be compared to the clear aperture position. The 3-1/2 digit digital voltmeter displays the clear aperture signal level as well as the ratio of the remaining sample positions relative to the clear aperture position. The sample wheel position is decoded so that the signals and ratios can be correlated to the data. The OEM is automatically reset to the I sub o on initial turn-on and can be reset to the '0' position by actuating a front panel switch. The sample wheel can be interrupted to change samples or induce a longer integration time if desired by a front panel command. Integration times from 1 - 50 seconds are provided at the front panel, and BCD data for external interfacing is provided.

Source record↗

Local oscillator distribution using a geostationary satellite

A satellite communication system suitable for distribution of local oscillator reference signals for a widely spaced microwave array has been developed and tested experimentally. The system uses a round-trip correction method of the satellite This experiment was carried out using Telstar-5, a commercial Ku-band geostationary satellite. For this initial experiment, both earth stations were located at the same site to facilitate direct comparison of the received signals. The local oscillator reference frequency was chosen to be 300MHz and was sent as the difference between two Ku-band tones. The residual error after applying the round trip correction has been measured to be better than 3psec for integration times ranging from 1 to 2000 seconds. For integration times greater then 500 seconds, the system outperforms a pair of hydrogen masers with the limitation believed to be ground-based equipment phase stability. The idea of distributing local oscillators using a geostationary satellite is not new; several researchers experimented with this technique in the eighties, but the achieved accuracy was 3 to 100 times worse than the present results. Since substantially and the performance of various components has improved. An important factor is the leasing of small amounts of satellite communication bandwidth. We lease three 100kHz bands at approximately one hundredth the cost of a full 36 MHz transponder. Further tests of the system using terminal separated by large distances and comparison tests with two hydrogen masers and radio interferometry is needed.

hydrogen maser↗

Secure boot, trusted boot and remote attestation for ARM TrustZone-based IoT Nodes

With the extensive application of IoT techniques, IoT devices have become ubiquitous in daily lives. Meanwhile, attacks against IoT devices have emerged to compromise IoT devices by tampering with system pre-installed programs or injecting new malware. To mitigate these attacks, integrity enforcement of IoT systems has been proposed. The integrity of an IoT device system includes load-time integrity and runtime integrity. In this paper, we design an IoT system based on ARM TrustZone to enforce the system integrity. First, we establish the root of trust and propose a hybrid booting approach consisting of both secure boot and trusted boot to enforce the system load-time integrity. Second, we investigate a paging-based process integrity measurement method to measure the NW processes and conduct remote attestation based on the measurement results ensuring the NW runtime process integrity. We implement an IoT prototype system on a NXP i.MX6Q SABRE SD development board to assess its feasibility. Finally, real-world experiment results demonstrate that our prototype introduces negligible performance overhead to the original system.

97 MATHEMATICS AND COMPUTING↗

Isotopic analysis of sub-nanogram neodymium loads using new ATONA™ amplifiers

The long-term gains of the ATONA™ amplifiers are shown to have less than 1 ppm variability. Repeat measurements of the JNdi-1 reference material demonstrate the ability of the ATONA™ amplification technology to make measurements of 143 Nd/ 144 Nd ratios with 23 ppm external reproducibility on 100-pg loads. The effect of increasing integration time on analytical reproducibility is also displayed as increasing integration time from 10 to 30s reduced the external measurement uncertainty from 37 to 23 ppm. These measurements represent an improvement of more than a factor of 3 in external measurement reproducibility relative to previously published 143 Nd/ 144 Nd measurements of 100-pg loads. Furthermore, this new technology will allow for the measurement of smaller samples for precise isotope ratios and open new avenues of research in the geological and nuclear forensic communities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

MEUMAPPS (C++ Version)

Many materials, metal alloys in particular, have features on the on micrometer or nanometer scale that have a large impact on the properties of the material. These features are known as the microstructure of the material. Understanding why and how the microstructure forms in a material is of fundamental scientific interest as well as of significant technological interest. The capability to predict microstructure evolution in a material allows the intentional design of microstructures and hence the intentional design of material properties. The phase-field method is one of the leading methods for predicting microstructure evolution. One of the most significant problems for phase-field models is their computational expense. Even limited phase-field simulations can easily require thousands of CPU core-hours to complete, which significantly limits their use. This code provides both a general framework for creating scalable, GPU-accelerated phase-field model applications as well as several applications themselves. The code is capable of using hundreds of GPUs efficiently, which greatly reduces the time required to perform simulations. The code is written with an emphasis on performance portability, that is the ability for the code to run efficiently on a number of different computing architectures without modification of the source code. The performance portability of this code is primarily enabled through the use of two libraries, Kokkos (performance portable data structures and execution patterns) and heFFTe (performance portable distributed 3D fast Fourier transforms). The code consists of a core library, applications, and tests. The core library includes shared functionality between applications. This includes interfaces with fast Fourier transform (FFT) libraries such as heFFTe, data structures based on Kokkos, file input and output capabilities, and a solver for infinitesimal strain mechanical equilibrium problems. Five applications are included in the code. The flagship application is the MEUMAPPS-SS application, which implements the Kim-Kim-Suzuki phase-field model for precipitation for an arbitrary number of phases and components in a metal alloy. Five simpler applications are also included that solve the Eshelby inclusion problem, Allen-Cahn equation, the coupled Allen-Cahn and diffusion equations, and the Cahn-Hilliard equation. The code includes two applications to solve the Cahn-Hilliard equation, one with constant-step-size first-order time integration and the second with adaptive high-order time integration.

DeWitt, Stephen [Oak Ridge National Lab. (ORNL), O↗

A rationale for human operator pulsive control behavior

When performing tracking tasks which involve demanding controlled elements such as those with K/s-squared dynamics, the human operator often develops discrete or pulsive control outputs. A dual-loop model of the human operator is discussed, the dominant adaptive feature of which is the explicit appearance of an internal model of the manipulator-controlled element dynamics in an inner feedback loop. Using this model, a rationale for pulsive control behavior is offered which is based upon the assumption that the human attempts to reduce the computational burden associated with time integration of sensory inputs. It is shown that such time integration is a natural consequence of having an internal representation of the K/s-squared-controlled element dynamics in the dual-loop model. A digital simulation is discussed in which a modified form of the dual-loop model is shown to be capable of producing pulsive control behavior qualitively comparable to that obtained in experiment.

Hess, R. A.↗

Optimal trajectories for aeroassisted, coplanar orbital transfer

Classical and minimax optimal control problems arising in the study of aeroassisted coplanar orbit transfer from a high planetary orbit to a low one are considered. Attention is given to (1) the minimization of the energy required for the maneuver; (2) minimization of the time integral of the heating rate; (3) minimization of the time of flight during the atmospheric portion of the trajectory; (4) maximization of the time of flight during the atmospheric portion of the trajectory; (5) minimization of the time integral of the path inclination; and (6) minimization of the sum of the squares of the entry and exit path inclinations.

Miele, A.↗

Confidence of the three-point estimator of frequency drift

It was shown that a three-point second difference estimator is nearly optimal for estimating frequency drift in many common atomic oscillators. A formula for the uncertainty of this estimate as a function of the integration time and of the Allan variance associated with this integration time is derived.

Weiss, Marc A.↗

Discrete integration of continuous Kalman filtering equations for time invariant second-order structural systems

A general form for the first-order representation of the continuous second-order linear structural-dynamics equations is introduced to derive a corresponding form of first-order continuous Kalman filtering equations. Time integration of the resulting equations is carried out via a set of linear multistep integration formulas. It is shown that a judicious combined selection of computational paths and the undetermined matrices introduced in the general form of the first-order linear structural systems leads to a class of second-order discrete Kalman filtering equations involving only symmetric sparse N x N solution matrices.

Park, K. C.↗

Towards improved speed and accuracy of laser powder bed fusion simulations via representation of multiple time scales

Here, we report the result of investigations into alternative representations of the temporal evolution of the laser powder bed fusion (LPBF) process. In order to provide the reader with sufficient context for the following discussion, the modeling challenge and its underlying cause, as well as current attempts to provide a satisfactory solution, will be discussed. Next, the author will describe two numerical methods (multirate time integration and parallel-in-time) to better represent the temporal scales of the problem. Finally, the results of applying these methods to a two-dimensional version of the heat transfer portion of the LPBF problem will be presented, reporting on both convergence and performance behavior, which indicate the possibility of significant speedup of solution calculations (40–100 times), along with well-characterized effects on solution accuracy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Parallel-in-Time Methods for Method-of-Lines Discretizations of Nonlinear Hyperbolic PDEs and Systems (Final Report)

The work for the subcontract is situated in the area of parallel-in-time integration for hyperbolic partial differential equations (PDEs). Parallel-in-time integration is an active area of research due to its ability to enable faster numerical simulations for applications throughout many areas of science. The work in this subcontract builds on a variety of results that were obtained, as part of the work performed for Subcontract No. B648355, for the Multigrid Reduction-in-Time (MGRIT) method from [1] applied to hyperbolic PDEs. This subcontract extends these results further to more efficient methods and to the case of method-of-lines discretizations for nonlinear hyperbolic PDES and systems of PDEs. The following is a summary of the research performed and results achieved during milestone periods 1, 2 and 3 by the PI (Hans De Sterck) and Postdoctoral Research Associate (Oliver Krzysik), for required tasks 1-4 (as listed in the Statement of Work): Research over the previous year has been split into three main projects: (i) solution of acoustic equation system; (ii) solution of nonlinear scalar hyperbolic PDEs; (iii) solution of nonlinear hyperbolic systems of PDEs.

97 MATHEMATICS AND COMPUTING↗