SEARCH · Engineering Papers
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.
Motional Stark effect (MSE) data acquired during large fast-ion stabilized sawteeth are critically reexamined. Here, the safety factor at the sawtooth crash changes by Δq ≃ 0.15, much more than any likely errors.
Analysis of the results from our investigations of the temperature dependence of ac susceptibilities (χ' and χ'') and magnetic viscosity in a 7.6 nm thin film of La 0.7 Sr 0.3 MnO 3 grown on SrTiO 3 (001) substrate is presented. The Curie temperature (T C ) of this film is magnetic field (H) dependent, varying as [T C (H) - T C (0)] ~ H n with T C (0) = 270 K and n = 0.15. Temperature dependence of χ' and χ'' measured at ac frequencies f m in the range of 0.1–10 kHz shows a broad peak near 230 K associated with the blocking temperature (T B ) of spin clusters present in the 1.4 nm surface layer and a frequency dependent peak near 270 K associated with T C . The mean relaxation time τ determined from the Cole–Cole analysis of the temperature dependent χ' and χ'' is shown to fit the Vogel–Fulcher law: τ = τ 0 exp[ΔE/k B (T - T 0 )] with T 0 = 245 K, ΔE/k B = 270 K = T C and τ 0 = 1.2 x 10 -9 s. Magnetic viscosity S (measured in H = 0 for a sample cooled in H = 50 Oe) determined from the fit to the magnetization: M (t)= M (0) – S ln t, with time t up to 2 h, shows a peak at 230 K above which M (2 h) switches to negative values for temperatures up to T C (0). It is argued that this negative magnetization results from an interaction between the superparamagnetic spin clusters and the ferromagnetic phase, both being present in the temperature range of T B < T < T C .
Experiments on the 1-MA, 220-ns COBRA generator at Cornell University were conducted to provide detailed measurements of structured cylindrical gas-puff z pinches. In the experiments, a 7 cm diameter triple-nozzle gas valve assembly with concentric outer and inner annular nozzles and a central gas jet initialize the z-pinch load with various working gases, radial density profiles, and externally applied axial magnetic fields. Planar laser-induced fluorescence provides a measure of the initial neutral gas density of the load, while three-frame laser shearing interferometry and multi-frame extreme ultraviolet (XUV) cameras reveal the formation and propagation of a magneto-Rayleigh–Taylor (MRT) unstable shock layer. Implosion trajectories are compared to simple, experimentally informed models and found to be in good agreement. Differences in the structure of the accelerating plasma sheath and evolution of the MRT instability are observed for different gas species and axial magnetic field strengths, correlating with differences in pinch uniformity and x-ray emission. Here, the average instability growth is compared to linear MRT theory predictions using the instantaneous acceleration of the best-fit implosion models and characteristic instability wavelength, with the effective Atwood number and seed perturbation size as fit parameters. For high density argon center jets, ionization prior to the arrival of the imploding plasma sheath suggests a heating mechanism consistent with photoionization by XUV self-emission.
We present the design and characterization of a large-area Cryogenic PhotoDetector designed for active particle identification in rare event searches, such as neutrinoless double beta decay and dark matter experiments. The detector consists of a 45.6 cm 2 surface area by a 1-mm-thick 10.6 g Si wafer. It is instrumented with a distributed network of Quasiparticle-trap-assisted Electrothermal feedback Transition-edge sensors with superconducting critical temperature T c = 41.5 mK to measure a thermal phonons released from interactions with photons.
Here, we performed ab initio molecular dynamics (AIMD) simulations to benchmark bulk liquid structures and to evaluate results from all-atom force field molecular dynamics (FFMD) simulations with the generalized Amber force field (GAFF) for organophosphorus (OP) and organochlorine (OC) compounds. Our work also addresses the current and important topic of force field validation, applied here to a set of nonaqueous organic liquids. Our approach differs from standard treatments, which validate force fields based on thermodynamic data. Utilizing radial distribution functions (RDFs), our results show that GAFF reproduces the AIMD-predicted asymmetric liquid structures moderately well for OP compounds that contain bulky alkyl groups. Among the OCs, RDFs obtained from FFMD overlap well with AIMD results, with some offsets in position and peak structuring. However, re-parameterization of GAFF for some OCs is needed to reproduce fully the liquid structures predicted by AIMD. The offsets between AIMD and FFMD peak positions suggest inconsistencies in the developed force fields, but, in general, GAFF is able to capture short-ranged and long-ranged interactions of OPs and OCs observed in AIMD. Along with the local coordination structure, we also compared enthalpies of vaporization. Overall, calculated bulk properties from FFMD compared reasonably well with experimental values, suggesting that small improvements within the FF should focus on parameters that adjust the bulk liquid structures of these compounds.
The design space for inertial confinement fusion (ICF) experiments is vast, and experiments are extremely expensive. Researchers rely heavily on computer simulations to explore the design space in search of high-performing implosions. However, ICF multiphysics codes must make simplifying assumptions, and thus deviate from experimental measurements for complex implosions. For more effective design and investigation, simulations require input from past experimental data to better predict future performance. In this work, we describe a cognitive simulation method for combining simulation and experimental data into a common, predictive model. This method leverages a machine learning technique called “transfer learning,” the process of taking a model trained to solve one task, and partially retraining it on a sparse dataset to solve a different, but related task. In the context of ICF design, neural network models are trained on large simulation databases and partially retrained on experimental data, producing models that are far more accurate than simulations alone. Here, we demonstrate improved model performance for a range of ICF experiments at the National Ignition Facility and predict the outcome of recent experiments with less than 10% error for several key observables. We discuss how the methods might be used to carry out a data-driven experimental campaign to optimize performance, illustrating the key product—models that become increasingly accurate as data are acquired.
Localized high concentration electrolytes have been proposed as an effective route to construct stable solid–electrolyte interphase (SEI) layers near Li-metal anodes. However, there is still a limited understanding of decomposition mechanisms of electrolyte components during SEI formation. In this work, we investigate reactivities of lithium bis(fluorosulfonyl)imide (LiFSI, salt), 1,2-dimethoxyethane (DME, solvent) and tris(2,2,2-trifluoroethyl)orthoformate (TFEO, diluent) species in DME+TFEO mixed solvents and 1M LiFSI/DME/TFEO solutions. By supplying excess of electrons into the simulation cell, LiFSI is initially reduced via a 4-electron charge transferring reaction yielding F - and N(SO 2 ) 2 3- . The local solvation environment has little effect on the subsequent TFEO reaction, which typically requires 6 |e| to decompose into F - , HCOO - , CH2CF and - OCH 2 CF 3 . Besides, TFEO dehydrogenation reaction mechanism under the attack of anions is also identified. Unlike salt and diluent, DME shows good stability with any excess of electrons. Here, the energetics of most relevant reactions are characterized. Most reactions are thermodynamically favorable with low activation barriers.
Inhomogeneities in a current-carrying conductor promote non-uniform heating and expansion through the complex feedback between current density, electrical resistivity, Ohmic heating, temperature, and hydrodynamics. Three-dimensional-magnetohydrodynamic (3D-MHD) simulations suggest that μm-scale resistive inclusions or voids seed local overheating and through hydrodynamic explosion generate continuously growing craters which become several times larger than the initial perturbation. The ejected mass is the genesis of an electrothermally driven plasma filament which develops at lower current than plasmas on uniform surfaces adjacent to the defect. This result suggests that 1D or even 2D treatments are largely inadequate for detailed prediction of plasma formation. To test computational predictions, z-pinch experiments driven to 1 MA studied ultra-high-purity aluminum rods which were then machined to include pairs of quasi-hemispherical voids or “engineered defects (ED)” on the current-carrying surface. ED are the dominant current-density perturbation and reproducibly drive local overheating which can be compared with 3D-MHD simulation. Data from high-resolution-gated imagers of visible surface emissions confirm many simulation predictions, including the surface topography of local overheating, and the propensity for neighboring ED to prematurely source plasmas which then connect to form a plasma filament. Results also provide conditional support of theory which suggests heating similarity; that is, heating is independent of ED size for geometrically scaled ED.
In this work, we introduce a diagrammatic approach to facilitate the automatic derivation of analytical nuclear gradients for tensor hyper-contraction (THC) based electronic structure methods. The automatically derived gradients are guaranteed to have the same scaling in terms of both operation count and memory footprint as the underlying energy calculations, and the computation of a gradient is roughly three times as costly as the underlying energy. The new diagrammatic approach enables the first cubic scaling implementation of nuclear derivatives for THC tensors fitted in molecular orbital basis (MO-THC). Furthermore, application of this new approach to THC-MP2 analytical gradients leads to an implementation, which is at least four times faster than the previously reported, manually derived implementation. Finally, we apply the new approach to the 14 tensor contraction patterns appearing in the supporting subspace formulation of multireference perturbation theory, laying the foundation for developments of analytical nuclear gradients and nonadiabatic coupling vectors for multi-state CASPT2.
The solid-electrolyte interphase (SEI) layer is a critical constituent of battery technology, which incorporates the use of lithium metals. Since the formation of the SEI is difficult to avoid, the engineering and harnessing of the SEI are absolutely critical to advancing energy storage. One problem is that much fundamental information about SEI properties is lacking due to the difficulty in probing a chemically complex interfacial system. One such property that is currently unknown is the dissolution of the SEI. This process can have significant effects on the stability of the SEI, which is critical to battery performance but is difficult to probe experimentally. Here, we report the use of ab initio computational chemistry simulations to probe the solution state properties of SEI components LiF, Li 2 O, LiOH, and Li 2 CO 3 in order to study their dissolution and other solution-based characteristics. Ab initio molecular dynamics was used to study the solvation structures of the SEI with a combination of radial distribution functions, discrete solvation structure maps, and vibrational density of states, which allows for the determination of free energies. From the change in free energy of dissolution, we determined that LiOH is the most likely component to dissolve in the electrolyte followed by LiF, Li 2 CO 3 , and Li 2 O although none were favored thermodynamically. Furthermore, this indicates that dissolution is not probable, but Li 2 O would make the most stable SEI with regard to dissolution in the electrolyte.
Previously, ELM-free and inter-ELM divertor peak heat flux reduction induced by an edge harmonic oscillation (EHO) was observed in NSTX (K.F. Gan et al., 2017 Nucl. Fusion 57 126053). This paper introduces new analysis of the EHO impact on the divertor heat flux. It was found that enhanced edge turbulence significantly altered the divertor heat flux footprint in NSTX H-modes, relative to ELM-free discharges. When the background edge fluctuation level was low, the EHO significantly increased the heat flux width and decreased the divertor peak heat flux. When the background edge fluctuation level was high, the EHO actually increased the divertor peak heat flux. As a result, it was also found that the heat flux width increased with increasing frequency of the EHO.
Here, the characteristics of the surface waves along the interface between a plasma and a dielectric material have been investigated using kinetic particle-in-cell simulations. A microwave source of GHz frequency has been used to trigger the surface wave in the system. The outcome indicates that the surface wave gets excited along the interface of plasma and the dielectric tube and appears as light and dark patterns in the electric field profiles. The dependency of radiation pressure on the dielectric permittivity and supplied input frequency has been investigated. Further, we assessed the capabilities of neural networks to predict the radiation pressure for a given system. The proposed deep neural network model is aimed at developing accurate and efficient data-driven plasma surface wave devices.
Gas-puff z-pinch implosions are characterized by the formation of a dense annular plasma shell, the sheath, that is driven to the axis by magnetic forces and therefore subject to the magneto-Rayleigh–Taylor instability. Here, the conditions within these sheaths are measured on the 1-MA COBRA generator at Cornell University for various gas species and initial fill densities. The gas-puff loads are initialized by a 7 cm diameter triple-nozzle gas valve assembly with concentric outer and inner annular nozzles and a central gas jet. Thomson scattering and laser interferometry provide spatially resolved flow, temperature, and electron density profiles midway through the implosion, while extreme ultraviolet pinhole cameras record the evolution of the plasma column and photoconducting diodes measure x-ray emission. Analysis of the scattering spectra includes a means of discriminating between thermal and non-thermal broadening to test for the presence of hydrodynamic turbulence. Two types of sheath profiles are observed, those with sharp discontinuities at the leading edge and those with smooth gradients. In both cases, non-thermal broadening is generally peaked at the front of the sheath and exhibits a characteristic decay length that roughly scales with the sheath ion mean free path. Here, we demonstrate that this non-thermal broadening term is inconsistent with laminar velocity gradients and is more consistent with dissipative turbulence driven by unstable plasma waves in a collisionless shock. The resulting differences in sheath profile are then set by the sheath ion collisionality in a manner consistent with recent 1D kinetic simulations.
This paper considers vertical plasma motion resulting from plasma current decay during the disruption event. The presented filament-based model describes the motion in the ideal wall limit as an adiabatically slow evolution of the plasma equilibrium. Here, the equilibrium exhibits a pitchfork bifurcation when the decaying plasma current passes a critical value determined by the external magnetic field. This bifurcation affects the disruption-induced mechanical loads on the first wall.
We present a computational investigation of the dependence of material erosion on the incident ion angle at rough graphite and silicon carbide divertor surfaces. Ion angle distributions (IADs) for D plasmas at NSTX-U and DIII-D divertors were calculated by an equation-of-motion model that traces the ion trajectories in the sheath. Then, the effective sputtering yields and ion shadowed area fractions were calculated by a Monte Carlo micro-patterning and roughness code that applied the calculated IADs to surface topographic data that were obtained from optical confocal microscopy of rough graphite and SiC divertor surfaces from NSTX-U and DIII-D experiments. The calculations found that the effective sputtering yields, the sputtering pattern, and the shadowed area are determined by the detailed surface topology rather than the root mean square roughness R RMS , which represents deviations from a flat surface. The suppression of the effective sputtering yields for rough surfaces compared to the yield for a smooth surface was accounted for by the change of the mean local incident ion angle (LIIA) (θ'). The mean surface inclination angle distribution (SIAD) (δ) was found to be a useful parameter to estimate the LIIA from the calculated IADs. We report global empirical formulas for the mean LIIA and fraction of the area shadowed from the main ions for D plasmas for rough surfaces with B-field incident angles α = 85°–89° as a function of the mean SIAD (δ). We propose the use of the mean LIIA (θ') to estimate the sputtering yield for rough surfaces from the angular dependence of the sputtering yield.
The bond strength and photodissociation dynamics of MgI+ are determined by a combination of theory, photodissociation spectroscopy, and photofragment velocity map imaging. From 17 000 to 21 500 cm−1, the photodissociation spectrum of MgI+ is broad and unstructured; photofragment images in this region show perpendicular anisotropy, which is consistent with absorption to the repulsive wall of the (1) Ω = 1 or (2) Ω = 1 states followed by direct dissociation to ground state products Mg+ (2S) + I (2P3/2). Analysis of photofragment images taken at photon energies near the threshold gives a bond dissociation energy D0(Mg+-I) = 203.0 ± 1.8 kJ/mol (2.10 ± 0.02 eV; 17 000 ± 150 cm−1). At photon energies of 33 000–41 000 cm−1, exclusively I+ fragments are formed. Over most of this region, the formation of I+ is not energetically allowed via one-photon absorption from the ground state of MgI+. Images show the observed product is due to resonance enhanced two-photon dissociation. The photodissociation spectrum from 33 000 to 38 500 cm−1 shows vibrational structure, giving an average excited state vibrational spacing of 227 cm−1. This is consistent with absorption to the (3) Ω = 0+ state from ν = 0, 1 of the (1) Ω = 0+ ground state; from the (3) Ω = 0+ state, absorption of a second photon results in dissociation to Mg* (3P°J) + I+ (3PJ). From 38 500 to 41 000 cm−1, the spectrum is broad and unstructured. We attribute this region of the spectrum to one-photon dissociation of vibrationally hot MgI+ at low energy and ground state MgI+ at higher energy to form Mg (1S) + I+ (3PJ) products.
Inertial confinement fusion capsules fielded at the National Ignition Facility are filled with deuterium and tritium fuel by means of a fill tube. The fill tube introduces a low-density pathway into the fuel region of the capsule that allows high Z contaminant to invade the hot spot during the course of the implosion. A recent series of nominally identical high-yield implosions on the NIF has exhibited significant variability in performance. We evaluate the impact of the fill tube in these implosions computationally to determine whether variations in fill tube geometry could have contributed to this variability. The main contrast between the fill tube geometry in the six shots was the outer diameter of the capsule bore hole, a conical hole into which the fill tube is inserted. In our simulations, the geometry of the bore hole can play a significant role in the development of nonlinear flows seeded by the fill tube. Additionally, we find that the amount of space between the bore hole and the fill tube is the primary factor that determines the amount of contaminant jetted into the hot spot by the fill tube and, in turn, the level of yield reduction due to the fill tube in our simulations. As a consequence, some capsules with 5 μm fill tubes are predicted to outperform capsules with 2 μm fill tubes. We also find that micrometer-scale changes to bore hole size can impact fusion yields by up to four times near the ignition threshold. Nevertheless, simulation trends do not reproduce experimental yield trends, suggesting that the fill tube geometry was not the primary factor contributing to the observed variability in performance and that the fill tube could be masking sensitivity to other asymmetries such as other micrometer-scale capsule defects like voids that were not included in our simulations.
An improvement in confinement associated with the injection of pellets has been observed in TJ-II and other devices. By means of a simple model, we study here the modification of turbulent transport by the injection of pellets in TJ-II and how this modification affects the particle confinement time. The analysis of the results shows the relationship of the confinement improvement with the evolution of the shear flow due to turbulence, especially near low order rational surfaces. This work is based on pellet injections into the neutral beam injection heated phase of TJ-II plasmas. For this, single, double, or triple hydrogen pellet injection are considered and the post-injection evolution of selected plasma parameters are analyzed and compared.