Search NASA⌕ Search

SEARCH · Search NASA

Results for “fast time simulation”

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 415 records · Page 23

Plasma flows during the ablation stage of an over-massed pulsed-power-driven exploding planar wire array

We characterize the plasma flows generated during the ablation stage of an over-massed exploding planar wire array, fielded on the COBRA pulsed-power facility (1 MA peak current, 250 ns rise time). The planar wire array is designed to provide a driving magnetic field (80–100 T) and current per wire distribution (about 60 kA), similar to that in a 10 MA cylindrical exploding wire array fielded on the Z machine. Over-massing the arrays enables continuous plasma ablation over the duration of the experiment without implosion. The requirement to over-mass on the Z machine necessitates wires with diameters of 75–100μm, which are thicker than wires usually fielded on wire array experiments. To test ablation with thicker wires, we perform a parametric study by varying the initial wire diameter between 33 and 100 μm. The largest wire diameter (100 μm) array exhibits early closure of the cathode-wire gap, while the gap remains open over the duration of the experiment for wire diameters between 33 and 75 μm. Laser plasma interferometry and time-gated extreme-ultraviolet (XUV) imaging are used to probe the plasma flows ablating from the wires. The plasma flows from the wires converge to generate a pinch, which appears as a fast-moving (V≈100kms−1) column of increased plasma density (n¯e≈2×1018cm−3) and strong XUV emission. Finally, we compare the results with three-dimensional resistive-magnetohydrodynamic (MHD) simulations performed using the code GORGON, the results of which reproduce the dynamics of the experiment reasonably well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Critical Assemblies: Dragon Burst Assembly and Solution Assemblies

This work reviews the historical literature associated with the Dragon experiment and water boiler reactors operated at Los Alamos during the Manhattan Project. Frisch’s invited talk given at the American Nuclear Society’s Fast Burst Reactor Conference held at the University of New Mexico in Albuquerque, New Mexico, in 1969 is quoted. From the literature review, basic models for the Dragon experiment and for a water boiler–type assembly (aqueous homogeneous reactor) were created that can be used for conducting multiphysics simulations for criticality excursion studies. This methodology utilizes the coupled neutronic-hydrodynamic method to perform a time-dependent dynamic simulation of a criticality excursion. MCNP® was utilized to calculate important nuclear kinetic parameters that were incorporated into the models. Simulation results compare reasonably well with historic data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Imaging multiphoton ionization dynamics of CH 3 I at a high repetition rate XUV free-electron laser

XUV multiphoton ionization of molecules is commonly used in free-electron laser experiments to study charge transfer dynamics. However, molecular dissociation and electron dynamics, such as multiple photon absorption, Auger decay, and charge transfer, often happen on competing time scales, and the contributions of individual processes can be difficult to unravel. We experimentally investigate the Coulomb explosion dynamics of methyl iodide upon core–hole ionization of the shallow inner-shell of iodine (4d) and classically simulate the fragmentation by phenomenologically introducing ionization dynamics and charge transfer. Under our experimental conditions with medium fluence and relatively long XUV pulses (~75 fs), we find that fast Auger decay prior to charge transfer significantly contributes to the charging mechanism, leading to a yield enhancement of higher carbon charge states upon molecular dissociation. Furthermore, we argue for the existence of another charging mechanism for the weak fragmentation channels leading to triply charged carbon atoms. This study shows that classical simulations can be a useful tool to guide the quantum mechanical description of the femtosecond dynamics upon multiphoton absorption in molecular systems.

Optics↗

Feasibility of an Active Interrogation System to Classify Waste with He-4 Neutron Spectroscopy

This work investigates a 4 He-detector active interrogation system that leverages neutron spectroscopy to classify nuclear waste streams. MCNP models tested the concept through the simulation of a D-D neutron generator, an array of 4 He detectors, and various waste compositions. The fast-neutron Differential Die-Away signature was augmented with a neutron-energy discrimination signature. This signature isolates induced fission neutrons, the energy of which is greater than that of the D-D monoenergetic spectrum. With the incorporation of this spectroscopic technique, the measurement time decreased by 3–9% (depending on the degree of neutron moderation and absorption presented by the sample), demonstrating how neutron spectroscopy can enhance active interrogation methods. The reduced measurement times would have significant financial and logistical benefits for facilities with large footprints of low-level waste production.

4He detectors↗

Sub-system quantum dynamics using coupled cluster downfolding techniques

In this paper, we discuss extending the sub-system embedding sub-algebra coupled cluster (SES-CC) formalism and the double unitary coupled cluster (DUCC) ansatz to the time domain. As we demonstrated in earlier studies, it is possible, using these formalisms, to calculate the energy of the entire system as an eigenvalue of downfolded/effective Hamiltonian in the active space, that is identifiable with the sub-system of the composite system. In these studies, we demonstrated that downfolded Hamiltonians integrate out Fermionic degrees of freedom that do not correspond to the physics encapsulated by the active space. We extend these results to the time-dependent Schrödinger equation, showing that a similar construct is possible to partition a system into a sub-system that varies slowly in time and a remaining subsystem that corresponds to fast oscillations. This time dependent formalism allows coupled cluster quantum dynamics to be extended to larger systems and for the formulation of novel quantum algorithms based on the quantum Lanczos approach, which have recently been considered in the literature.

coupled cluster, Electron correlation, quantum dyn↗

Benchmarking fast-neutron response of 7 Li-enriched C 7 LYC scintillators: Simulations and experiment

New insights into the fast neutron response of 7 Li-enriched C 7 LYC scintillators have been obtained through detailed Geant4 simulations of 35 Cl(n,p) and 35 Cl(n, α) reactions that are primarily responsible for the fast neutron response of C 7 LYC. The light output of a 1” x1” C 7 LYC detector was simulated using the Birks equation for light quenching in scintillators. The simulations were benchmarked against neutron scattering data from a 12 C target obtained with a white neutron beam at the Los Alamos LANSCE facility and fourteen 1”x1” C 7 LYC detectors. Incident neutron energies extracted from time-of-flight measurements, and appropriate recoil corrections for the different detector angles were used to map the response of C 7 LYC to mono-energetic neutrons. Resonances in the 35 Cl(n,p) reaction cross-section below 3 MeV incident neutron energy were directly observed in C 7 LYC. The experimental data and simulations are compared with ENDF databases, as well as recent direct cross-section measurements of neutrons scattered from a 35 Cl target. In conclusion, the results suggest a scrutiny of the 35 Cl(n,p) and 35 Cl(n, α) cross sections at low energy in the ENDF library.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Topology induced crossover between Langevin, subdiffusion, and Brownian diffusion regimes in supercooled water

Despite extensive studies of supercooled water, it remains challenging to understand its peculiar dynamic anomalous properties. In this work, we integrated full atomistic simulations of supercooled water over the temperature range of room temperature to 200 K using a quantum-mechanics-based polarizable force field with the dressed dynamics method that couples fast collision events and slow reorganization dynamics of hydrogen-bond networks. Our analysis unveils the salient multiscale features in the transient relaxation dynamics of supercooled water. A classical Langevin behavior dominates at fast timescales, while long-time relaxations unveil two different activation barriers in two temperature regions: below and above 230 K. The modulation of the entropy spectrum by temperature is elucidated in terms of a three-state model underlined by the complexity of the water dynamics associated with a topological transition of a strong hydrogen-bond network. This state-dependent network topology is quantitatively characterized by power-law exponents of inverse network connectivity from 200 to 298 K. Finally, this work provides valuable guidance for further studies on the transient relaxation dynamics of supercooled water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Application of Monte Carlo Algorithms to Cardiac Imaging Reconstruction

Monte Carlo algorithms have a growing impact on nuclear medicine reconstruction processes. One ofthe main limitations of myocardial perfusion imaging (MPI) is the effective mitigation of the scattering component,which is particularly challenging in Single Photon Emission Computed Tomography (SPECT). In SPECT,no timing information can be retrieved to locate the primary source photons. Monte Carlo methods allow anevent-by-event simulation of the scattering kinematics, which can be incorporated into a model of the imagingsystem response. This approach was adopted in the late Nineties by several authors, and recently took advantageof the increased computational power made available by high-performance CPUs and GPUs. These recent developmentsenable a fast image reconstruction with improved image quality, compared to deterministic approaches.Deterministic approaches are based on energy-windowing of the detector response, and on the cumulative estimateand subtraction of the scattering component. In this paper, we review the main strategies and algorithms tocorrect the scattering effect in SPECT and focus on Monte Carlo developments, which nowadays allow the threedimensionalreconstruction of SPECT cardiac images in a few seconds.

Pharmacology & Pharmacy↗

Laboratory Study of Magnetic Reconnection in Lunar-relevant Mini-magnetospheres

Abstract Mini-magnetospheres are small ion-scale structures that are well suited to studying kinetic-scale physics of collisionless space plasmas. Such ion-scale magnetospheres can be found on local regions of the Moon, associated with the lunar crustal magnetic field. In this paper, we report on the laboratory experimental study of magnetic reconnection in laser-driven, lunar-like ion-scale magnetospheres on the Large Plasma Device at the University of California, Los Angeles. In the experiment, a high-repetition rate (1 Hz), nanosecond laser is used to drive a fast-moving, collisionless plasma that expands into the field generated by a pulsed magnetic dipole embedded into a background plasma and magnetic field. The high-repetition rate enables the acquisition of time-resolved volumetric data of the magnetic and electric fields to characterize magnetic reconnection and calculate the reconnection rate. We notably observe the formation of Hall fields associated with reconnection. Particle-in-cell simulations reproducing the experimental results were performed to study the microphysics of the interaction. By analyzing the generalized Ohm’s law terms, we find that the electron-only reconnection is driven by kinetic effects through the electron pressure anisotropy. These results are compared to recent satellite measurements that found evidence of magnetic reconnection near the lunar surface.

Astronomy & Astrophysics↗

Quantifying Impacts of Biomass Pelletization on Fast Pyrolysis Using a Single-Particle Reactor, X-ray Computed Tomography, and Computational Modeling

The pore structure and density of lignocellulosic feedstocks dictate intraparticle transport phenomena and thereby play an important role in thermochemical conversion processes such as fast pyrolysis for biofuel and biochemical production. Variations in microstructure are inherent from different biomass species and can be introduced by preprocessing techniques such as cutting and pelletization. Morphological changes also occur during conversion and lead to vastly different pore structures and behavior during pyrolysis, which impact required conversion times and product distributions. The current work presents a comprehensive comparison of fast pyrolysis of neat and pelletized pine feedstocks, which includes single-particle experiments, modeling, and 3D imaging by X-ray computed tomography (XCT). The particle-scale model included anisotropic heat and mass transport in a shrinking particle with pyrolysis reactions based on the CRECK mechanism with boundary conditions informed by reactor-scale simulations of the single-particle reactor. The models were validated by measurements of the temperature and mass loss from single-particle pyrolysis experiments of neat and pelletized pine. Quantitative analysis of XCT geometries revealed that pyrolytic conversion yielded chars with increased porosity and permeability compared to the unpyrolyzed materials, along with decreased tortuosity and anisotropy. Pelletization of the pine feedstock resulted in a much denser, less permeable material, which converted slower and produced more residual char after pyrolysis compared to neat pine. The results from particle modeling revealed that accounting for the dynamic and anisotropic heat and mass transport caused by differences in pore structure is critical to achieving agreement with experimental results. Overall, this study highlights the dramatic differences in conversion behavior imparted by pelletization and the importance of capturing microstructural attributes in computational models to guide the design and optimization of pyrolysis processes for specific biomass feedstocks.

09 BIOMASS FUELS↗

Assessment of the associated particle technique with high-resolution gamma-ray spectroscopy for in-field identification of chemical warfare agents and explosives

In this work, high-resolution prompt γ-ray neutron activation analysis (PGNAA) is performed using the associated particle (AP) technique with a deuterium-tritium neutron generator and high-purity germanium detector. Although the time resolution is inferior compared to similar systems employing fast scintillation γ-ray detectors, the strong background suppression combined with the high γ-ray energy resolution provides an important boost in sensitivity to certain key elements, especially relevant to field measurements of complex chemicals encased within thick layers of metal, such as in the case of recovered chemical warfare (CW) munitions. Results for various CW and explosive simulants contained within mock munitions are presented and compared with results obtained from traditional PGNAA measurements without the aid of the AP technique. Significant improvements in identification and discrimination between certain chemicals are achieved with the high energy resolution AP technique.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Characterization of the Fast-Neutron Irradiator and the Fast-Flux Tube Irradiation Fixtures at the Pennsylvania State Breazeale Reactor

Accurate knowledge of the neutron spectrum at a nuclear research reactor is a prerequisite for planning irradiation experiments, as well as for evaluating irradiation exposure results. The neutron-flux spectrum in the fast-neutron irradiator (FNI) and the fast-flux tube (FFT) irradiation fixtures at the Pennsylvania State Breazeale Reactor (PSBR) were characterized using the multi-foil neutron activation method. These irradiation fixtures make use of graded shielding to produce unique neutron fields. Multiple foil sets were irradiated in the fixtures with different exposure times and reactor powers to understand the stability over a wide range of operating conditions. Measured results were evaluated against a MCNP6 simulation to produce a measurement-informed neutron flux-energy spectrum for each fixture using STAYSL_PNNL. Simulated estimates of the FNI fixture, a newer fixture (~25 years old), demonstrated excellent agreement with measured results; the FFT did not. Thermal neutron measurements from the FFT suggest there is additional thermal leakage not captured in the simulation model. In conclusion, possible explanations for the discrepancy include burn-out or degradation (i.e., micro-cracking, gaps, etc.) in the boral and cadmium liners over the lifetime of the fixture (~40 years old).

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Solid, structured composite neutron detectors with high dynamic range capability

Neutron detectors are essential across disciplines such as fundamental science, nuclear security, safeguards, and civilian applications. While 3 He-filled gas proportional counters have long been revered for their efficacy in detecting thermal neutrons and praised for their efficiency, neutron/gamma discrimination, and stability, the scarcity of 3He has spurred a search for alternatives. Here, we explore a solid structured scintillating particle composite (SPC) consisting of 6 Li-containing scintillating glass particles within an acrylic matrix as a neutron detector for high dynamic range applications. We show for the first time that an SPC neutron detector can boast an intrinsic detection efficiency of 0.261% for pure 252 Cf fission neutrons and an overall neutron detection efficiency of (0.546 ± 0.003)% at the Neutron Free-in-Air facility while being able to function in an intense gamma-ray environment. We also show that the SPC neutron detector supports fast neutron capture times and enables a dual-readout scheme that extends the detector dynamic range to high incident neutron fluxes. A scalable fabrication process allows for tailoring the SPC detector properties to the requirements of specific applications. Good agreement is found between the experimental results taken with a National Institute of Standards and Technology traceable 252 Cf source and the coupled MCNP6 and optical-ray-tracing simulations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Extremely massive disc galaxies in the nearby Universe form through gas-rich minor mergers

ABSTRACT In our hierarchical structure-formation paradigm, the observed morphological evolution of massive galaxies – from rotationally supported discs to dispersion-dominated spheroids – is largely explained via galaxy merging. However, since mergers are likely to destroy discs, and the most massive galaxies have the richest merger histories, it is surprising that any discs exist at all at the highest stellar masses. Recent theoretical work by our group has used a cosmological, hydrodynamical simulation to suggest that extremely massive (M* > 1011.4 M⊙) discs form primarily via minor mergers between spheroids and gas-rich satellites, which create new rotational stellar components and leave discs as remnants. Here, we use UV-optical and H i data of massive galaxies, from the Sloan Digital Sky Survey, Galaxy Evolution Explorer, Dark Energy Camera Legacy Survey (DECaLS), and Arecibo Legacy Fast ALFA surveys, to test these theoretical predictions. Observed massive discs account for ∼13 per cent of massive galaxies, in good agreement with theory (∼11 per cent). ∼64 per cent of the observed massive discs exhibit tidal features, which are likely to indicate recent minor mergers, in the deep DECaLS images (compared to ∼60 per cent in their simulated counterparts). The incidence of these features is at least four times higher than in low-mass discs, suggesting that, as predicted, minor mergers play a significant (and outsized) role in the formation of these systems. The empirical star formation rates agree well with theoretical predictions and, for a small galaxy sample with H i detections, the H i masses and fractions are consistent with the range predicted by the simulation. The good agreement between theory and observations indicates that extremely massive discs are indeed remnants of recent minor mergers between spheroids and gas-rich satellites.

79 ASTRONOMY AND ASTROPHYSICS↗

Flexible Resource Scheduler for FAST-DERMS (FRS-FASTDERMS) v0.9

The Flexible Resource Scheduler is a hierarchical controller that manages the distributed energy resources in a distribution substation or distribution feeder to provide a firm commitment of power flow at the substation or feeder head to be scheduled in transmission-level markets as an aggregated demand resource. It is the reference controller for the FAST-DERMS Architecture, developed in tandem with the architecture under the DOE FAST-DERMS project. It is comprised of a day-ahead stochastic optimization, which schedules substation power flow and reserves, a intra-hour MPC, which generates dispatch base points for DER, and a real-time PID controller maintaining that dispatches DER to maintain the substation power around the base points. The repository also includes a representative aggregator controller, and all of the necessary components to run a simulation using PNNL's GridAPPS-D software with the controller.

MacDonald, Jason [Lawrence Berkeley National Labor↗

Multiscale Approach to Fast ModSim for Laser Processing of Metals for Future Nuclear Deterrence Environments

Predicting performance of parts produced using laser-metal processing remains an out- standing challenge. While many computational models exist, they are generally too computationally expensive to simulate the build of an engineering-scale part. This work develops a reduced order thermal model of a laser-metal system using analytical Green's function solutions to the linear heat equation, representing a step towards achieving a full part performance prediction in an "overnight" time frame. The developed model is able to calculate a thermal history for an example problem 72 times faster than a traditional FEM method. The model parameters are calibrated using a non-linear solution and microstructures and residual stresses calculated and compared to a non-linear case. The calibrated model shows promising agreement with a non-linear solution.

36 MATERIALS SCIENCE↗

GPU-enabled extreme-scale turbulence simulations: Fourier pseudo-spectral algorithms at the exascale using OpenMP offloading

Fourier pseudo-spectral methods for nonlinear partial differential equations are of wide interest in many areas of advanced computational science, including direct numerical simulation of three-dimensional (3-D) turbulence governed by the Navier-Stokes equations in fluid dynamics. This paper presents a new capability for simulating turbulence at a new record resolution up to 35 trillion grid points, on the world's first exascale computer, Frontier, comprising AMD MI250x GPUs with HPE's Slingshot interconnect and operated by the US Department of Energy's Oak Ridge Leadership Computing Facility (OLCF). Key programming strategies designed to take maximum advantage of the machine architecture involve performing almost all computations on the GPU which has the same memory capacity as the CPU, performing all-to-all communication among sets of parallel processes directly on the GPU, and targeting GPUs efficiently using OpenMP offloading for intensive number-crunching including 1-D Fast Fourier Transforms (FFT) performed using AMD ROCm library calls. With 99% of computing power on Frontier being on the GPU, leaving the CPU idle leads to a net performance gain via avoiding the overhead of data movement between host and device except when needed for some I/O purposes. Memory footprint including the size of communication buffers for MPI_ALLTOALL is managed carefully to maximize the largest problem size possible for a given node count. Detailed performance data including separate contributions from different categories of operations to the elapsed wall time per step are reported for five grid resolutions, from 2048 3 on a single node to 32768 3 on 4096 or 8192 nodes out of 9408 on the system. Both 1D and 2D domain decompositions which divide a 3D periodic domain into slabs and pencils respectively are implemented. The present code suite (labeled by the acronym GESTS, GPUs for Extreme Scale Turbulence Simulations) achieves a figure of merit (in grid points per second) exceeding goals set in the Center for Accelerated Application Readiness (CAAR) program for Frontier. The performance attained is highly favorable in both weak scaling and strong scaling, with notable departures only for 2048 3 where communication is entirely intra-node, and for 32768 3 , where a challenge due to small message sizes does arise. Communication performance is addressed further using a lightweight test code that performs all-to-all communication in a manner matching the full turbulence simulation code. Performance at large problem sizes is affected by both small message size due to high node counts as well as dragonfly network topology features on the machine, but is consistent with official expectations of sustained performance on Frontier. Overall, although not perfect, the scalability achieved at the extreme problem size of 32768 3 (and up to 8192 nodes — which corresponds to hardware rated at just under 1 exaflop/sec of theoretical peak computational performance) is arguably better than the scalability observed using prior state-of-the-art algorithms on Frontier's predecessor machine (Summit) at OLCF. New science results for the study of intermittency in turbulence enabled by this code and its extensions are to be reported separately in the near future.

3D fast Fourier transform↗

A first look at AS-ROC: a Si-Ge integrated chip readout for fast timing

Abstract Advances in timing detector technology require new specialized readout electronics. Applications demand below 10 ps time of arrival resolution and low power for a low repetition rate. A possible path to achieve O(10 ps) time resolution is an integrated chip using Silicon Germanium (SiGe) technology. Using DoE SBIR funding, Anadyne, Inc., in collaboration with UC Santa Cruz, has developed a prototype SiGe front-end readout chip optimized for low power and timing resolution. Two versions of the chip were produced with performance in simulation: a more power version with 10 ps resolution at 5 fC with 1.1 mW/channel, and a less power version with 10 ps resolution at 8 fC with 0.6 mW/channel.The chip was produced at Tower Semiconductor with 350 nm technology. The ASIC from the prototype run shows good performance: a rise time of 0.7–1 ns and 25 mV per fC response with RMS noise <1 mV. Simulation and results from the prototype will be reported in this paper.

Instruments & Instrumentation↗