SEARCH · Search NASA
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.
The linearity of the equations of motion of dynamical systems such as spacecraft or manipulators when controlled by acceleration commands is exploited in generating exact fast maneuvers. A simulated robot example thereof is examined, with advantages and limitations of the proposed design approach discussed in its context.
PeleLMeX's unique capabilities are allowing for reactive flow modeling at unprecedented scales and a reasonable time and cost. The code is currently being extended to tackle more practical, design-oriented simulations by implementing Large Eddy Simulation and data-driven chemical models, providing a fast but accurate tool for engineers considering the emergence of GPU-accelerated platforms. These extensions are critical for enabling the physical insight required to design the next generation of combustion devices as a key component of a renewable energy future.
In Fermilab's PIP-II machine protection system, beam loss signals from various detectors are digitized at 125 MS/s. Noise from both high-frequency sources and low-frequency 60 Hz AC power equipment can contaminate the data. To suppress noise across these ranges especially 60 Hz and its harmonics, which overlap with beam loss signal frequencies advanced digital processing beyond standard filtering is required. Several real-time functional blocks were simulated and tested on an FPGA: (1) a dual time-constant discharging integrator filter, (2) a de-ripple baseline extraction and storage block, and (3) a fast-recovery discharging integrator. The nonlinear IIR integrator filter removes high-frequency noise and feeds into the baseline extractor. Upon detecting abrupt beam loss, it switches to a longer time constant to prevent baseline distortion. The de-ripple block calculates a valid baseline by averaging over multiple 60 Hz periods, storing results in a 4096-word FPGA RAM. This baseline is subtracted from raw data before integration by the fast-recovery block, which resets quickly after use. All blocks achieved expected performance.
The quantum simulation of entangled spin systems can play a central role in quantum magnetic materials discovery. Additionally, the simulation of spectroscopic signatures, such as the dynamical structure factor accessed in inelastic neutron scattering (INS), necessitates a long timescale for circuit evolution. This is because the energy resolution is directly related to the time over which the circuit could be meaningfully evolved. However, canonical Trotterization requires deep circuits precluding such long-time evolution—even for a small number of qubits. Here, in this study, we demonstrate “direct” resource efficient fast-forwarding (REFF) measurements with short-depth circuits that can be used to capture longer time dynamics of spin Hamiltonians. We showcase the results of the dynamics of a quantum spin dimer, the basic quantum unit of emergent many-body spin systems, whose density of states we simulate accurately. The long temporal evolution and measurement of the two-spin correlation functions enable the calculation of the dynamical structure factor S(Q = 0, ω) measured in the neutron scattering cross-section. We exhibit the clarity of the triplet gap and the triplet splitting of the quantum dimer with class-leading fidelity that enables comparison to experimental neutron data. Our results on current circuit hardware outline an important workflow to predict and benchmark against the outputs of INS experiments of quantum magnets.
The kinematics of a 1-D electron hole is studied using a novel Particle-In-Cell simulation code. A hole tracking technique enables us to follow the trajectory of a fast-moving solitary hole and study quantitatively hole acceleration and coupling to ions. We observe a transient at the initial stage of hole formation when the hole accelerates to several times the cold-ion sound speed. Artificially imposing slow ion speed changes on a fully formed hole causes its velocity to change even when the ion stream speed in the hole frame greatly exceeds the ion thermal speed, so there are no reflected ions. The behavior that we observe in numerical simulations agrees very well with our analytic theory of hole momentum conservation and the effects of “jetting.”
Model-based optimization offers a systematic approach to advanced scenario planning. In this case, the feedforward-control inputs (actuator trajectories) that are needed to attain and sustain a desired scenario are obtained by solving a nonlinear constrained optimization problem. This class of problems generally minimize a cost function that measures the difference between desired and actual plasma states. Several numerical optimization algorithms, such as sequential quadratic programming, require repeated calculation of the cost function gradients with respect to the input trajectories. Calculating these gradients numerically can be computationally intensive, increasing the time needed to solve the feedforward-control optimization problem. Here, this work introduces a method to analytically calculate these cost function gradients from the current profile evolution model. This can significantly reduce the computational time and allow for fast feedforward-control optimization, which would eventually enable optimal scenario planning between discharges. The performance of the feedforward optimizer with analytical gradients is compared to a traditional optimization algorithm based on numerical gradients for different NSTX-U scenarios. The plasma dynamics in the optimization algorithm are simulated using the Control Oriented Transport SIMulator (COTSIM). Results of the work show that analytical gradients consistently reduce the computation time while achieving trajectories that are comparable to those obtained by traditional optimization algorithms based on numerical gradients.
We simulate the dynamics, including laser cooling, of three-dimensional (3-D) ion crystals confined in a Penning trap using a newly developed molecular dynamics-like code. The numerical integration of the ions’ equations of motion is accelerated using the fast multipole method to calculate the Coulomb interaction between ions, which allows us to efficiently study large ion crystals with thousands of ions. In particular, we show that the simulation time scales linearly with ion number, rather than with the square of the ion number. By treating the ions’ absorption of photons as a Poisson process, we simulate individual photon scattering events to study laser cooling of 3-D ellipsoidal ion crystals. Initial simulations suggest that these crystals can be efficiently cooled to ultracold temperatures, aided by the mixing of the easily cooled axial motional modes with the low frequency planar modes. In our simulations of a spherical crystal of 1000 ions, the planar kinetic energy is cooled to several millikelvin in a few milliseconds while the axial kinetic energy and total potential energy are cooled even further. This suggests that 3-D ion crystals could be well suited as platforms for future quantum science experiments.
Synchrotron micro- and nanoprobe beamlines have demonstrated great potential to advance photovoltaic devices. Most importantly, their small X-ray spot size has enabled the direct correlation of electrical performance with elemental composition at subgrain resolution for a variety of polycrystalline solar cells. Whereas the bulk of most inorganic semiconductors is stable under the high X-ray flux of focused X-ray beams, semiconductors with organic components are prone to a variety of degradation mechanisms. This is particularly critical to evaluate for the emerging organometal halide perovskite solar cells. In this study, we investigate the effects of hard X-rays on the nanoscale performance and elemental distribution of these solar cells. We show that their composition does not change during common operando and in situ measurements at synchrotron nanoprobes. However, we found a significant X-ray-induced electronic degradation of solar cells with methylammonium lead iodide absorbers. Time- and dose-dependent measurements unveiled two characteristic degradation time constants on the order of 12 and 200 s that are independent of the X-ray flux. On the basis of heat and dose simulations, we attribute the fast decay to the dose-driven creation of recombination centers, while the slow decay is compatible with the observation of compositional changes. Finally, we detail how degradation-induced measurement artifacts can be outrun and showcase the high correlation of the X-ray-beam-induced current with the iodine and lead distribution.
Self-assembly of colloidal nanocrystals (NCs) into superlattices (SLs) is an appealing strategy to design hierarchically organized materials with promising functionalities. Mechanistic studies are still needed to uncover the design principles for SL self-assembly, but such studies have been difficult to perform due to the fast time and short length scales of NC systems. To address this challenge, we developed an apparatus to directly measure the evolving phases in situ and in real time of an electrostatically stabilized Au NC solution before, during, and after it is quenched to form SLs using small-angle X-ray scattering. By developing a quantitative model, we fit the time-dependent scattering patterns to obtain the phase diagram of the system and the kinetics of the colloidal and SL phases as a function of varying quench conditions. The extracted phase diagram is consistent with particles whose interactions are short in range relative to their diameter. We find the degree of SL order is primarily determined by fast (subsecond) initial nucleation and growth kinetics, while coarsening at later times depends nonmonotonically on the driving force for self-assembly. We validate these results by direct comparison with simulations and use them to suggest dynamic design principles to optimize the crystallinity within a finite time window. The combination of this measurement methodology, quantitative analysis, and simulation should be generalizable to elucidate and better control the microscopic self-assembly pathways of a wide range of bottom-up assembled systems and architectures.
Herein we present the first experimental study of plasmoid formation in a magnetic reconnection layer undergoing rapid radiative cooling, a regime relevant to extreme astrophysical plasmas. Two exploding aluminum wire arrays, driven by the Z machine, generate a reconnection layer (S L ≈ 120) in which the cooling rate far exceeds the hydrodynamic transit rate (τ hydro =τ cool > 100). The reconnection layer generates a transient burst of >1 keV x-ray emission, consistent with the formation and subsequent rapid cooling of the layer. Time-gated x-ray images show fast-moving (up to 50 km s –1 ) hotspots in the layer, consistent with the presence of plasmoids in 3D resistive magnetohydrodynamic simulations. X-ray spectroscopy shows that these hotspots generate the majority of Al K-shell emission (around 1.6 keV) prior to the onset of cooling, and exhibit temperatures (170 eV) much greater than that of the plasma inflows and the rest of the reconnection layer, thus providing insight into the generation of high-energy radiation in radiatively cooled reconnection events.
Modeling studies of the interaction of a delta wing aircraft with direct lightning strikes were carried out using an approximate scale model of an F-106B. The model, which is three feet in length, is subjected to direct injection of fast current pulses supplied by wires, which simulate the lightning channel and are attached at various locations on the model. Measurements are made of the resulting transient electromagnetic fields using time derivative sensors. The sensor outputs are sampled and digitized by computer. The noise level is reduced by averaging the sensor output from ten input pulses at each sample time. Computer analysis of the measured fields includes Fourier transformation and the computation of transfer functions for the model. Prony analysis is also used to determine the natural frequencies of the model. Comparisons of model natural frequencies extracted by Prony analysis with those for in flight direct strike data usually show lower damping in the in flight case. This is indicative of either a lightning channel with a higher impedance than the wires on the model, only one attachment point, or short streamers instead of a long channel.
Fast, efficient parallel algorithms are presented for discrete event simulations of dynamic channel assignment schemes for wireless cellular communication networks. The driving events are call arrivals and departures, in continuous time, to cells geographically distributed across the service area. A dynamic channel assignment scheme decides which call arrivals to accept, and which channels to allocate to the accepted calls, attempting to minimize call blocking while ensuring co-channel interference is tolerably low. Specifically, the scheme ensures that the same channel is used concurrently at different cells only if the pairwise distances between those cells are sufficiently large. Much of the complexity of the system comes from ensuring this separation. The network is modeled as a system of interacting continuous time automata, each corresponding to a cell. To simulate the model, conservative methods are used; i.e., methods in which no errors occur in the course of the simulation and so no rollback or relaxation is needed. Implemented on a 16K processor MasPar MP-1, an elegant and simple technique provides speedups of about 15 times over an optimized serial simulation running on a high speed workstation. A drawback of this technique, typical of conservative methods, is that processor utilization is rather low. To overcome this, new methods were developed that exploit slackness in event dependencies over short intervals of time, thereby raising the utilization to above 50 percent and the speedup over the optimized serial code to about 120 times.
The Fast Scattering Code (FSC) is a frequency domain noise prediction program developed at the NASA Langley Research Center (LaRC) to simulate the acoustic field produced by the interaction of known, time harmonic incident sound with bodies of arbitrary shape and surface impedance immersed in a potential flow. The code uses the equivalent source method (ESM) to solve an exterior 3-D Helmholtz boundary value problem (BVP) by expanding the scattered acoustic pressure field into a series of point sources distributed on a fictitious surface placed inside the actual scatterer. This work provides additional code validation studies and illustrates the range of code parameters that produce accurate results with minimal computational costs. Systematic noise prediction studies are presented in which monopole generated incident sound is scattered by simple geometric shapes - spheres (acoustically hard and soft surfaces), oblate spheroids, flat disk, and flat plates with various edge topologies. Comparisons between FSC simulations and analytical results and experimental data are presented.
This paper describes the adoption of a Test Driven Development approach and a Continuous Integration System in the development of the Trick Simulation Toolkit, a generic simulation development environment for creating high fidelity training and engineering simulations at the NASA Johnson Space Center and many other NASA facilities. It describes the approach, and the significant benefits seen, such as fast, thorough and clear test feedback every time code is checked into the code repository. It also describes an approach that encourages development of code that is testable and adaptable.
Accurate thermal scattering laws (TSLs) are essential for reliable neutron transport simulations, particularly for systems with low-energy neutrons. Pulsed-neutron die-away (PNDA) experiments provide a highly sensitive platform for TSL validation but face challenges from fast neutron leakage and material constraints in nonhydrogenous targets such as beryllium. This work investigates the use of high-density polyethylene (HDPE) as a premoderator to reduce leakage, improve thermalization, and shorten experimental run times while preserving sensitivity to the Be TSL. Here, three premoderated configurations, encapsulation, slab, and interstitial, were evaluated using MCNP6.3® simulations with ENDF/B-VIII.0 nuclear data. The effects of the Be TSL was quantified by comparing decay constants from simulations with and without 𝑆(𝛼,𝛽) treatments for both beryllium and HDPE. The results show that premoderation enables the use of as little as 5.20% of the beryllium volume required for an unreflected geometry, achieving uncertainties below 0.5% and reducing run time by up to 72.2% in the studied configurations. The interstitial configuration achieved the best balance between low statistical uncertainty and high TSL sensitivity, outperforming the slab and encapsulation designs. These findings demonstrate that carefully optimized premoderation can significantly enhance the efficiency, feasibility, and precision of PNDA experiments for Be TSL validation.
This study develops more computationally efficient versions of the TwO-Moment Aerosol Sectional (TOMAS) microphysics algorithms, collectively called Fast TOMAS. Several methods for speeding up the algorithm were attempted, but only reducing the number of size sections was adopted. Fast TOMAS models, coupled to the GISS GCM II-prime, require a new coagulation algorithm with less restrictive size resolution assumptions but only minor changes in other processes. Fast TOMAS models have been evaluated in a box model against analytical solutions of coagulation and condensation and in a 3-D model against the original TOMAS (TOMAS-30) model. Condensation and coagulation in the Fast TOMAS models agree well with the analytical solution but show slightly more bias than the TOMAS-30 box model. In the 3-D model, errors resulting from decreased size resolution in each process (i.e., emissions, cloud processing wet deposition, microphysics) are quantified in a series of model sensitivity simulations. Errors resulting from lower size resolution in condensation and coagulation, defined as the microphysics error, affect number and mass concentrations by only a few percent. The microphysics error in CN70CN100 (number concentrations of particles larger than 70100 nm diameter), proxies for cloud condensation nuclei, range from 5 to 5 in most regions. The largest errors are associated with decreasing the size resolution in the cloud processing wet deposition calculations, defined as cloud-processing error, and range from 20 to 15 in most regions for CN70CN100 concentrations. Overall, the Fast TOMAS models increase the computational speed by 2 to 3 times with only small numerical errors stemming from condensation and coagulation calculations when compared to TOMAS-30. The faster versions of the TOMAS model allow for the longer, multi-year simulations required to assess aerosol effects on cloud lifetime and precipitation.
We propose a general framework leveraging the galaxy–halo connection to link galaxies observed at different redshift in a statistical way, and use the link to infer the redshift evolution of the galaxy population. Our tests based on hydrodynamic simulations show that our method can accurately recover the stellar mass assembly histories up to z ∼ 3 for present star-forming and quiescent galaxies down to 10 10 h −1 M ⊙ . Applying the method to observational data shows that the stellar mass evolution of the main progenitors of galaxies depends strongly on the properties of descendants, such as stellar mass, halo mass, and star formation states. Galaxies hosted by low-mass groups/haloes at the present time have since z ∼ 1.8 grown their stellar mass ∼2.5 times as fast as those hosted by massive clusters. This dependence on host halo mass becomes much weaker for descendant galaxies with similar star formation states. Star-forming galaxies grow about 2–4 times faster than their quiescent counterparts since z ∼ 1.8. Both TNG and EAGLE simulations overpredict the progenitor stellar mass at z > 1, particularly for low-mass descendants.
The design of a new parallel-processing digital simulator is described. The new simulator has been developed specifically for analysis of wind energy systems in real time. The new processor has been named: the Wind Energy System Time-domain simulator, version 3 (WEST-3). Like previous WEST versions, WEST-3 performs many computations in parallel. The modules in WEST-3 are pure digital processors, however. These digital processors can be programmed individually and operated in concert to achieve real-time simulation of wind turbine systems. Because of this programmability, WEST-3 is very much more flexible and general than its two predecessors. The design features of WEST-3 are described to show how the system produces high-speed solutions of nonlinear time-domain equations. WEST-3 has two very fast Computational Units (CU's) that use minicomputer technology plus special architectural features that make them many times faster than a microcomputer. These CU's are needed to perform the complex computations associated with the wind turbine rotor system in real time. The parallel architecture of the CU causes several tasks to be done in each cycle, including an IO operation and the combination of a multiply, add, and store. The WEST-3 simulator can be expanded at any time for additional computational power. This is possible because the CU's interfaced to each other and to other portions of the simulation using special serial buses. These buses can be 'patched' together in essentially any configuration (in a manner very similar to the programming methods used in analog computation) to balance the input/ output requirements. CU's can be added in any number to share a given computational load. This flexible bus feature is very different from many other parallel processors which usually have a throughput limit because of rigid bus architecture.