Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sim”

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 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Dynamics of multilayer Rayleigh–Taylor instability at moderately high Atwood numbers

This paper investigates the multilayer Rayleigh–Taylor instability (RTI) using statistically stationary experiments conducted in a gas tunnel. Employing diagnostics such as particle image velocimetry (PIV) and planar laser induced fluorescence (PLIF), we make simultaneous velocity–density measurements to study how dynamics and mixing are linked in this variable density flow. Experiments are conducted in a newly built, blow-down three-layer gas tunnel facility. Mixing between three gas streams is studied, where the top and bottom streams are comprised of air, and the middle stream is an air–helium mixture. Shear is minimized between these streams by matching their inlet velocities. The four experimental conditions investigated here consist of two different density ratios (Atwood numbers 0.3 and 0.6), each investigated at two instability development times (or equivalently, two streamwise locations), and all experiments are with the same middle stream thickness of 3 cm. The growth of the middle layer is measured using laser-based planar Mie scattering visualization. The mixing width is found to grow linearly with time at late times. Various quantitative measures of molecular mixing indicate a very high degree of molecular mixing at late times in the multilayer RTI flow. The vertical turbulent mass flux $a_y$ is calculated. In addition to mostly negative values of $a_y$ , typical of buoyancy-dominated flows due to negative correlation between velocity and density fluctuations, positive regions are also observed in profiles of $a_y$ due to entrainment and erosion at the lower edge of the mixing region. Global energy budgets are calculated for the multilayer RTI flow at late times and it is found that the majority of potential energy released has been dissipated due to viscous effects, and a large value of mixing efficiency ( $\sim$ 60 %) is observed.

Mechanics↗
DOE OSTI2022

Discovery of the most luminous quasar of the last 9 Gyr

Abstract We report the discovery of a bright ( $g = 14.5$ mag (AB), $K = 11.9$ mag (Vega)) quasar at redshift $z=0.83$ — the optically brightest (unbeamed) quasar at $z>0.4$ . SMSS J114447.77-430859.3, at a Galactic latitude of $b=+18.1^{\circ}$ , was identified by its optical colours from the SkyMapper Southern Survey (SMSS) during a search for symbiotic binary stars. Optical and near-infrared spectroscopy reveals broad Mg ii , H $\unicode{x03B2}$ , H $\unicode{x03B1}$ , and Pa $\unicode{x03B2}$ emission lines, from which we measure a black hole mass of $\log_{10}\! (M_{\mathrm{BH}}/\mathrm{M}_{\odot}) = 9.4 \pm 0.5$ . With its high luminosity, $L_{\mathrm{bol}} = (4.7\pm1.0)\times10^{47}\,\mathrm{erg\,s}^{-1}$ or $M_{i}(z=2) = -29.74$ mag (AB), we estimate an Eddington ratio of $\approx1.4$ . As the most luminous quasar known over the last ${\sim}$ 9 Gyr of cosmic history, having a luminosity $8\times$ greater than 3C 273, the source offers a range of potential follow-up opportunities.

Astronomy & Astrophysics↗
DOE OSTI2017

Particle acceleration in relativistic magnetic flux-merging events

Using analytical and numerical methods (fluid and particle-in-cell simulations) we study a number of model problems involving merger of magnetic flux tubes in relativistic magnetically dominated plasma. Mergers of current-carrying flux tubes (exemplified by the two-dimensional ‘ABC’ structures) and zero-total-current magnetic flux tubes are considered. In all cases regimes of spontaneous and driven evolution are investigated. We identify two stages of particle acceleration during flux mergers: (i) fast explosive prompt X-point collapse and (ii) ensuing island merger. The fastest acceleration occurs during the initial catastrophic X-point collapse, with the reconnection electric field of the order of the magnetic field. During the X-point collapse, particles are accelerated by charge-starved electric fields, which can reach (and even exceed) values of the local magnetic field. The explosive stage of reconnection produces non-thermal power-law tails with slopes that depend on the average magnetization $\unicode[STIX]{x1D70E}$ . For plasma magnetization $\unicode[STIX]{x1D70E}\leqslant 10^{2}$ the spectrum power-law index is $p>2$ ; in this case the maximal energy depends linearly on the size of the reconnecting islands. For higher magnetization, $\unicode[STIX]{x1D70E}\geqslant 10^{2}$ , the spectra are hard, $p<2$ , yet the maximal energy $\unicode[STIX]{x1D6FE}_{\text{max}}$ can still exceed the average magnetic energy per particle, ${\sim}\unicode[STIX]{x1D70E}$ , by orders of magnitude (if $p$ is not too close to unity). The X-point collapse stage is followed by magnetic island merger that dissipates a large fraction of the initial magnetic energy in a regime of forced magnetic reconnection, further accelerating the particles, but proceeds at a slower reconnection rate.

Physics↗
DOE OSTI2018

Diffusive shock re-acceleration

Here, we have performed two-dimensional hybrid simulations of non-relativistic collisionless shocks in the presence of pre-existing energetic particles (‘seeds’); such a study applies, for instance, to the re-acceleration of galactic cosmic rays (CRs) in supernova remnant (SNR) shocks and solar wind energetic particles in heliospheric shocks. Energetic particles can be effectively reflected and accelerated regardless of shock inclination via a process that we call diffusive shock re-acceleration. We determine that re-accelerated seeds can drive the streaming instability in the shock upstream and produce effective magnetic field amplification. This can later trigger the injection of thermal protons even at oblique shocks that ordinarily cannot inject thermal particles. We characterize the current in reflected seeds, finding that it tends to a universal value $J\simeq en_{\text{CR}}v_{\text{sh}}$ , where $en_{\text{CR}}$ is the seed charge density and $v_{\text{sh}}$ is the shock velocity. When applying our results to SNRs, we find that the re-acceleration of galactic CRs can excite the Bell instability to nonlinear levels in less than ${\sim}10~\text{yr}$ , thereby providing a minimum level of magnetic field amplification for any SNR shock. Finally, we discuss the relevance of diffusive shock re-acceleration also for other environments, such as heliospheric shocks, galactic superbubbles and clusters of galaxies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Collisional effects on resonant particles in quasilinear theory

A careful examination of the effects of collisions on resonant wave–particle interactions leads to an alternate interpretation and deeper understanding of the quasilinear operator originally formulated by Kennel & Engelmann ( Phys. Fluids , vol. 9, 1966, pp. 2377–2388) for collisionless, magnetized plasmas, and widely used to model radio frequency heating and current drive. The resonant and nearly resonant particles are particularly sensitive to collisions that scatter them out of and into resonance, as for Landau damping as shown by Johnston ( Phys. Fluids , vol. 14, 1971, pp. 2719–2726) and Auerbach ( Phys. Fluids , vol. 20, 1977, pp. 1836–1844). As a result, the resonant particle–wave interactions occur in the centre of a narrow collisional boundary when the collision frequency $\unicode[STIX]{x1D708}$ is very small compared to the wave frequency $\unicode[STIX]{x1D714}$ . The diffusive nature of the pitch angle scattering combined with the wave–particle resonance condition enhances the collision frequency by $(\unicode[STIX]{x1D714}/\unicode[STIX]{x1D708})^{2/3}\gg 1$ , resulting in an effective resonant particle collisional interaction time of $\unicode[STIX]{x1D70F}_{\text{int}}\sim (\unicode[STIX]{x1D708}/\unicode[STIX]{x1D714})^{2/3}/\unicode[STIX]{x1D708}\ll 1/\unicode[STIX]{x1D708}$ . A collisional boundary layer analysis generalizes the standard quasilinear operator to a form that is fully consistent with Kennel–Englemann, but allows replacing the delta function appearing in the diffusivity with a simple integral (having the appropriate delta function limit) retaining the new physics associated with the narrow boundary layer, while preserving the entropy production principle. The limitations of the collisional boundary layer treatment are also estimated, and indicate that substantial departures from Maxwellian are not permitted.

Physics↗
DOE OSTI2021

Gradient-based optimization of 3D MHD equilibria

Using recently developed adjoint methods for computing the shape derivatives of functions that depend on magnetohydrodynamic (MHD) equilibria (Antonsen et al. , J. Plasma Phys. , vol. 85, issue 2, 2019; Paul et al. , J. Plasma Phys. , vol. 86, issue 1, 2020), we present the first example of analytic gradient-based optimization of fixed-boundary stellarator equilibria. We take advantage of gradient information to optimize figures of merit of relevance for stellarator design, including the rotational transform, magnetic well and quasi-symmetry near the axis. With the application of the adjoint method, we reduce the number of equilibrium evaluations by the dimension of the optimization space ( ${\sim }50\text {--}500$ ) in comparison with a finite-difference gradient-based method. We discuss regularization objectives of relevance for fixed-boundary optimization, including a novel method that prevents self-intersection of the plasma boundary. We present several optimized equilibria, including a vacuum field with very low magnetic shear throughout the volume.

Physics↗
DOE OSTI2022

A learned closure method applied to phase mixing in a turbulent gradient-driven gyrokinetic system in simple geometry

We present a new method for formulating closures that learn from kinetic simulation data. We apply this method to phase mixing in a simple gyrokinetic turbulent system – temperature-gradient-driven turbulence in an unsheared slab. The closure, called the learned multi-mode (LMM) closure, is constructed by, first, extracting an optimal basis from a nonlinear kinetic simulation using singular value decomposition. Subsequent nonlinear fluid simulations are projected onto this basis and the results are used to formulate the closure. We compare the closure with other closures schemes over a broad range of the relevant two-dimensional parameter space (collisionality and gradient drive). We find that the turbulent kinetic system produces phase-mixing rates much lower than the linear expectations, which the LMM closure is capable of capturing. We also compare radial heat fluxes. A Hammett–Perkins closure, generalized to include collisional effects, is quite successful throughout the parameter space, producing ${\sim }14\,\%$ root-mean-square (r.m.s.) error. The LMM closure is also very effective: when trained at three (two) points (in a 35 point parameter grid), the LMM closure produces $8\,\%$ ( $12\,\%$ ) r.m.s. errors. The LMM procedure can be readily generalized to other closure problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Towards fast and accurate predictions of radio frequency power deposition and current profile via data-driven modelling: applications to lower hybrid current drive

Three machine learning techniques (multilayer perceptron, random forest and Gaussian process) provide fast surrogate models for lower hybrid current drive (LHCD) simulations. A single GENRAY/CQL3D simulation without radial diffusion of fast electrons requires several minutes of wall-clock time to complete, which is acceptable for many purposes, but too slow for integrated modelling and real-time control applications. The machine learning models use a database of more than 16 000 GENRAY/CQL3D simulations for training, validation and testing. Latin hypercube sampling methods ensure that the database covers the range of nine input parameters ( $n_{e0}$ , $T_{e0}$ , $I_p$ , $B_t$ , $R_0$ , $n_{\|}$ , $Z_{{\rm eff}}$ , $V_{{\rm loop}}$ and $P_{{\rm LHCD}}$ ) with sufficient density in all regions of parameter space. The surrogate models reduce the inference time from minutes to $\sim$ ms with high accuracy across the input parameter space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The structure of 3-D collisional magnetized bow shocks in pulsed-power-driven plasma flows

Here we investigate three-dimensional (3-D) bow shocks in a highly collisional magnetized aluminium plasma, generated during the ablation phase of an exploding wire array on the MAGPIE facility (1.4 MA, 240 ns). Ablation of plasma from the wire array generates radially diverging, supersonic ($M_S$ $\sim 7$), super-Alfvénic ($M_A>1$) magnetized flows with frozen-in magnetic flux $(R_M\gg1)$. These flows collide with an inductive probe placed in the flow, which serves both as the obstacle that generates the magnetized bow shock, and as a diagnostic of the advected magnetic field. Laser interferometry along two orthogonal lines of sight is used to measure the line-integrated electron density. A detached bow shock forms ahead of the probe, with a larger opening angle in the plane parallel to the magnetic field than in the plane normal to it. Since the resistive diffusion length of the plasma is comparable to the probe size, the magnetic field decouples from the ion fluid at the shock front and generates a hydrodynamic shock, whose structure is determined by the sonic Mach number, rather than the magnetosonic Mach number of the flow. The 3-D simulations performed using the resistive magnetohydrodynamic (MHD) code Gorgon confirm this picture, but under-predict the anisotropy observed in the shape of the experimental bow shock, suggesting that non-MHD mechanisms may be important for modifying the shock structure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Multimode theory of electron hole transverse instability

We present Vlasov–Poisson three-dimensional linear stability analysis of an initially planar electron hole structure, solving for the distribution function by integration along unperturbed orbits. The non-sinusoidal potential perturbation shape (parallel to $B$ ) is expanded in eigenfunctions of the adiabatic Poisson operator, generalizing the prior assumption of a rigid shift of the equilibrium. We show that the shiftmode is then modified by a second discrete mode plus an integral over a continuum of wave-like modes. A rigorous treatment shows that the continuum can be approximated effectively by a single mode that satisfies the external wave dispersion relation, thus making the perturbation a weighted sum of three modes. We find numerically the solution for the complex instability frequency, and the corresponding three mode amplitudes determining the perturbation eigenmode. This multimode analysis refines the accuracy of the prior single-mode results, giving slightly higher growth rates at most parameters, as expected from the extra mode shape freedom. Oscillating modes near stability boundaries have larger mode distortions which help explain particle-in-cell simulations that observe instability up to ${\sim }20$ % beyond the prior shiftmode thresholds, and narrowing of the perturbation. At high magnetic field, the multimode analysis predicts a reduction of the already small growth rate.

Physics↗
DOE OSTI2023

Fast-ion transport in quasisymmetric equilibria in the presence of a resonant Alfvénic perturbation

Significant progress has been made in designing magnetic fields that provide excellent confinement of the guiding-centre trajectories of alpha particles using quasisymmetry (QS). Given the reduction in this transport channel, we assess the impact of resonant Alfvén eigenmodes (AEs) on the guiding-centre motion. The AE amplitudes are chosen to be consistent with experimental measurements and large-scale simulations. We evaluate the drift resonance condition, phase-space island width and island overlap criterion for quasisymmetric configurations. Kinetic Poincaré plots elucidate features of the transport, including stiff transport above a critical perturbation amplitude. Our analysis highlights key departures from the AE-driven transport in tokamaks, such as the avoidance of phase-space island overlap in quasihelical configurations and the enhanced transport due to wide phase-space islands in low magnetic shear configurations. In configurations that are closer to QS, with QS deviations $\delta B/B_0 \lesssim 10^{-3}$ , the transport is primarily driven by the AE, while configurations that are further from QS, $\delta B/B_0 \sim 10^{-2}$ , experience significant transport due to the QS-breaking fields in addition to the AE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A buffer-gas trap for the NEPOMUC positron beam: optimization studies with electrons

Buffer-gas traps (BGTs) use inelastic interactions with nitrogen molecules to capture positrons from a continuous beam. These devices are invaluable for high-resolution studies of matter–antimatter interactions, antihydrogen research and positronium laser spectroscopy. We present a new project with the goal of producing a non-neutral plasma containing ${\sim }10^8$ low-energy positrons by installing a BGT on the NEPOMUC (NEutron induced POsitron source MUniCh) high-intensity positron beam. Details of the BGT are outlined and results are presented from experiments in which an electron beam, with a similar intensity and energy spread to the remoderated NEPOMUC beam, was used to create pulses of non-neutral electron plasma. The device is a vital component of the APEX (A Positron Electron eXperiment) project, which aims to create a low-temperature electron–positron pair plasma.

Physics↗
DOE OSTI2020

Recent Advances in Nuclear Forensic Chemistry

Nuclear forensics is a rapidly growing field with a widespread impact on global nonproliferation. From the beginning, nuclear forensic capabilities have been inextricably bound to existing limitations within chemistry and materials science. As nuclear forensics has continued to develop as a unique field of research, new advances in chemistry have revolutionized our ability to rapidly and quantitatively analyze radioactive materials for attribution. In particular, microanalytical techniques originally developed for nanoscience, such as micro-XRF, SIMS, and SEM-EDX, have enabled spatially resolved chemical analyses to be performed with only nanograms of material. Such techniques have redefined the paradigm for pre- and postdetonation nuclear forensics, as large quantities of material are no longer necessary for full characterization and attribution. This Review explains how modern nuclear forensic techniques have developed from their origins in chemistry and provides real-world examples of how these techniques have been used to characterize radioactive materials, including fallout from nuclear detonations such as the 1945 Trinity test.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2016

In Situ Molecular Imaging of the Biofilm and Its Matrix

Molecular mapping of live biofilms at submicron resolution presents a grand challenge. Here, we present the first chemical mapping results of biofilm extracellular polymeric sub-stance (EPS) components in biofilms using correlative imaging be-tween super resolution florescence microscopy and liquid time-of-flight secondary ion mass spectrometry (ToF-SIMS). Shewanella oneidensis is used as a model organism. Heavy metal anions chro-mate (Cr2O72-) consisting of chromium Cr (VI) was a model envi-ronmental stressor used to treat the biofilms. Of particular interest, biologically relevant water clusters have been first observed in the biofilms. Characteristic fragments of biofilm matrix components such as proteins, polysaccharides, and lipids can be spatially im-aged. Furthermore, characteristic fatty acids (e.g., palmitic acid), quinolone signal, and riboflavin fragments are found to respond af-ter the biofilm is treated with Cr (VI), leading to biofilm dispersion. Significant changes in water clusters and quorum sensing signals indicative of intercellular communication in the aqueous environ-ment are observed, suggesting that they might result in fatty acid synthesis and inhibit riboflavin production. The Cr (VI) reduction seems to follow the Mtr pathway leading to Cr (III) formation. Our approach potentially opens a new avenue for mechanistic insight of microbial community processes and communications using in situ imaging mass spectrometry and superresolution optical micros-copy.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2023

Strain-Mediated Sn Incorporation and Segregation in Compositionally Graded Ge 1– x Sn x Epilayers Grown by MBE at Different Temperatures

We investigated the process of Sn incorporation and surface segregation for compositionally graded Ge 1–x Sn x epilayers grown on high-quality Ge (001) substrates. Here, the growth resulted in pseudomorphic GeSn layers with a ~6% maximal Sn fraction at a constant substrate temperature. The maximal fraction of Sn was increased to 9.0% when the growth temperature was continuously lowered while increasing the Sn flux. The analysis of surface droplets and SIMS profiles of elemental composition give evidence of Sn rejection during the growth, potentially associated with a critical energy of elastic strain. The intentional reduction of the coherent strain by decreasing the Sn flux near the sample surface has been shown to trap a higher fraction of Sn in the Ge 1–x Sn x layer and lower surface segregation. Our results demonstrate that strain relief by misfit dislocations in the compositionally graded layer is inhibited, which leads to Sn segregation. Specifically, the compressive strain in the graded Ge 1–x Sn x epilayer is effectively “zero” near the interface with the Ge substrate and increases up to about –1.5 × 10 –2 near the surface. Thus, although the nucleation of a dislocation may reduce the compressive strain for the top region of the epilayer, it is not beneficial for the bottom region.

42 ENGINEERING↗
DOE OSTI2019

Graph Networks as a Universal Machine Learning Framework for Molecules and Crystals

Graph networks are a new machine learning (ML) paradigm that supports both relational reasoning and combinatorial generalization. Here, we develop universal MatErials Graph Network (MEGNet) models for accurate property prediction in both molecules and crystals. We demonstrate that the MEGNet models outperform prior ML models such as the SchNet in 11 out of 13 properties of the QM9 molecule data set. Similarly, we show that MEGNet models trained on $\sim 60,000$ crystals in the Materials Project substantially outperform prior ML models in the prediction of the formation energies, band gaps and elastic moduli of crystals, achieving better than DFT accuracy over a much larger data set. Here, we present two new strategies to address data limitations common in materials science and chemistry. First, we demonstrate a physically-intuitive approach to unify four separate molecular MEGNet models for the internal energy at 0 K and room temperature, enthalpy and Gibbs free energy into a single free energy MEGNet model by incorporating the temperature, pressure and entropy as global state inputs. Second, we show that the learned element embeddings in MEGNet models encode periodic chemical trends and can be transfer-learned from a property model trained on a larger data set (formation energies) to improve property models with smaller amounts of data (band gaps and elastic moduli).

36 MATERIALS SCIENCE↗
DOE OSTI2016

Determining the syringyl/guaiacyl lignin ratio in the vessel and fiber cell walls of transgenic Populus plants

Observation of the spatial lignin distribution throughout the plant cell wall provides insight into the physicochemical characteristics of lignocellulosic biomass. The distribution of syringyl (S) and guaiacyl (G) lignin in cell walls of a genetically modified Populus deltoides and its corresponding empty vector control were analyzed with time-of-flight secondary ion mass spectrometry (ToF-SIMS) and then mapped to determine the S/G lignin ratio of the sample surface and specific regions of interest (ROIs). The surface characterizations of transgenic cross-sections within 1 cm vertical distance of each other on the stem possess similar S/G lignin ratios. Furthermore, the analysis of the ROIs determined that there was a 50% decrease in the S/G lignin ratio of the transgenic xylem fiber cell walls.

09 BIOMASS FUELS↗
DOE OSTI2023

Molecular Imaging Reveals Two Distinct Mixing States of PM 2.5 Particles Sampled in a Typical Beijing Winter Pollution Case

Mixing states of aerosol particles are crucial for understanding the role of aerosols in influencing air quality and climate. However, a fundamental understanding of the complex mixing states is still lacking because most traditional analysis techniques only reveal bulk chemical and physical properties with limited surface and 3-D information. In this research, 3-D molecular imaging enabled by ToF-SIMS was used to elucidate the mixing states of PM 2.5 samples obtained from a typical Beijing winter haze event. In light pollution cases, a thin organic layer covers separated inorganic particles, while in serious pollution cases, ion exchange and an organic-inorganic mixing surface on large-area particles were observed. Finally, the new results provide key 3-D molecular information of mixing states, which is highly desirable for reducing uncertainty and bias in representing aerosol-cloud interactions in current Earth System Models, and improving the understanding of aerosols on air quality and human health.

3-D molecular imaging↗