Comparison of Reactive Burn Equilibrium Closure Assumptions in CTH.
Abstract not provided.
Ruggirello, Kevin↗
SEARCH · Engineering Papers
Results for “Comparison”
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.
Comparison of edge turbulence characteristics between DIII-D and C-Mod simulations with XGC1
The physical processes taking place at the separatrix and scrape-off layer regions are crucial for the operation of tokamaks as they govern the interaction of hot plasma with the vessel walls. Numerical modeling of the edge with state-of-the-art codes attempts to elucidate the complex interactions between neoclassical drifts, turbulence, poloidal and parallel flows that control the physical set-up of the SOL region. Here, we present post-processing analysis of simulation results from the gyrokinetic code XGC1, comparing and contrasting edge turbulence characteristics from a simulation of the DIII-D tokamak against a simulation of the Alcator C-Mod tokamak. We find that the equilibrium ExB flux across the separatrix has a similar poloidal pattern in both discharges which can be explained by rB-drifts and trapped ion excursions. Furthermore, collisionality is noted to play a major role in the way that it prevents local charge accumulations from having more global effects in the C-Mod case. In both cases, turbulent electron heat flux is observed to be higher than the ion one and is possibly related to the need of electrons to maintain quasineutrality through the only channel available to them for exiting the confinement. By Fourier analysis, we identify turbulent frequencies and growth rates of the dominant mode in both simulations. In the case of C-Mod, these numbers point to the presence of a drift wave. In the DIII-D case, further linear simulations with the Gene code reveal a trapped electron mode. Furthermore, using a blob detection and tracking tool, we present the amplitude and size distributions of the blobs from both simulations. The amplitude distributions are in qualitative agreement with experimental observations while the size distributions are consistent with the fact that most of the blobs are not connecting to the divertor plates and suggest that they are generated by the shearing of the turbulent modes.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Comparison of modular analytical wake models to the Lillgrund wind plant
Efficient and accurate wake models are required for wind plant performance modeling and the suite of engineering processes that support wind plant layout, control, and monitoring. Although many analytical and engineering wake models with low computational costs have been proposed, their ability to represent the power production of large wind plants in a wide range of atmospheric conditions is not completely understood. The following validation study reviews the underlying theory for analytical wake models, outlines quality control procedures for observational data, and compares model results with observational data from the Lillgrund Wind Plant. Lillgrund makes a valuable case study for wake modeling because of its regular arrangement and the relatively close spacing of constituent wind turbines, which lead to regular and significant wake interactions within the wind plant and the development of deep array flow conditions. Formulations for the velocity deficit, wake-added turbulence, and wake superposition methods are considered in a modular sense, yielding many possible configurations to represent wind turbine wakes, of which seven are examined in detail. Velocity deficit models that account for flow conditions in the near wake are better able to reproduce power production for wind turbines in the transitional region of the wind plant, where wind turbines experience as many as five wakes from upstream turbines. Finally, in the deep array, where power production reaches asymptotic values and wake statistics become quasi-periodic, wake superposition schemes become the largest driver in error reduction; using the linear or maximum wake superposition methods can reduce the relative root mean square error by as much as 40% in the deep array.
17 WIND ENERGY↗
Surface structure of linear nanopores in amorphous silica: Comparison of properties for different pore generation algorithms
In this work, we compare the surface structure of linear nanopores in amorphous silica (a-SiO 2 ) for different versions of “pore drilling” algorithms (where the pores are generated by the removal of atoms from the preformed bulk a-SiO 2 ) and for “cylindrical resist” algorithms (where a-SiO 2 is formed around a cylindrical exclusion region). After adding H to non-bridging O, the former often results in a moderate to high density of surface silanol groups, whereas the latter produces a low density. The silanol surface density for pore drilling can be lowered by a final dehydroxylation step, and that for the cylindrical resist approach can be increased by a final hydroxylation step. In this respect, the two classes of algorithms are complementary. We focus on the characterization of the chemical structure of the pore surface, decomposing the total silanol density into components corresponding to isolated and vicinal mono silanols and geminal silanols. The final dehyroxylation and hydroxylation steps can also be tuned to better align some of these populations with the target experimental values.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Comparison of 1D and 2D particle-in-cell simulations for DC magnetron sputtering discharges
One-dimensional and two-dimensional particle-in-cell (PIC)/Monte Carlo collision simulations are performed for a balanced direct current magnetron sputtering (DCMS) discharge. The plasma parameters obtained by both simulations above the target racetrack share similar features and magnitudes under similar discharge intensities, i.e., similar local discharge current densities. The commonly observed high-energy tail of electron energy distribution in DCMS discharges is reproduced in both simulations. Additionally, the accuracy and applicability of using onedimensional simulations to qualitatively investigate the characteristics of balanced DCMS discharge are confirmed. In terms of capturing the key physical mechanisms in DCMS, one-dimensional PIC simulation is an efficient method, which could largely alleviate the computational expense and preserve the physical fidelity of modeling results.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Self-compensation in heavily Ge doped AlGaN: A comparison to Si doping
Self-compensation in Ge- and Si-doped Al 0.3 Ga 0.7 N has been investigated in terms of the formation of III vacancy and donor-vacancy complexes. Both Ge- and Si-doped AlGaN layers showed a compensation knee behavior with impurity compensation (low doping regime), compensation plateau (medium doping regime), and self-compensation (high doping regime). A maximum free carrier concentration of 4–5 × 10 19 cm –3 was obtained by Ge doping, whereas Si doping resulted in only half of that value, ~2 × 10 19 cm –3 . A DFT calculation with the grand canonical thermodynamics model was developed to support the hypothesis that the difference in self-compensation arises from the difference in the formation energies of the V III -n•donor complexes relative to their onsite configurations. The model suggested that the V III -2•donor and V III -3•donor complexes were responsible for self-compensation for both Ge- and Si-doped AlGaN. However, a lower free carrier concentration in Si-doped samples was due to a high V III -3•Si concentration, resulting from a lower energy of formation of V III -3•Si.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Comparison of local and global gyrokinetic calculations of collisionless zonal flow damping in quasi-symmetric stellarators
The linear collisionless damping of zonal flows is calculated for quasi-symmetric stellarator equilibria in flux-tube, flux-surface, and full-volume geometry. Equilibria are studied from the quasi-helical symmetry configuration of the Helically Symmetric eXperiment (HSX), a broken symmetry configuration of HSX, and the quasi-axial symmetry geometry of the National Compact Stellarator eXperiment (NCSX). Zonal flow oscillations and long-time damping affect the zonal flow evolution, and the zonal flow residual goes to zero for small radial wavenumber. The oscillation frequency and damping rate depend on the bounce-averaged radial particle drift in accordance with theory. While each flux tube on a flux surface is unique, several different flux tubes in HSX or NCSX can reproduce the zonal flow damping from a flux-surface calculation given an adequate parallel extent. The flux-surface or flux-tube calculations can accurately reproduce the full-volume long-time residual for moderate kx, but the oscillation and damping time scales are longer in local representations, particularly for small kx approaching the system size.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
A comparison of Raman and pyrometry dynamic temperature measurements of shocked cyclohexane
Optical pyrometry is a well-developed, broadly applicable method to determine temperature in dynamic compression experiments. However, it measures radiation from only the skin depth of the radiant material, is highly susceptible to spurious backgrounds, and requires an often-unmeasurable free parameter—the sample dynamic emissivity. Raman spectroscopy offers a method to determine temperature directly from fundamental theory, but it is practically limited in applications to a small subset of shock experiments. In this work, we describe an experiment allowing simultaneous application of Raman and pyrometric temperature techniques to benchmark a specific instance of pyrometry. We attempted both measurements on multiple experiments and compare the successful temperature results between the separate but highly repeatable shots. Here, we find that dynamic shock temperatures inferred from Raman spectroscopy and pyrometry generally agree within the uncertainties of each measurement, although those uncertainties are larger than typical for a specially designed pyrometry experiment on a shocked metal sample.
74 ATOMIC AND MOLECULAR PHYSICS↗
Comparison of water nanodroplet properties on different graphite-based substrates
The molecular structure and dynamics of water differ considerably at various interfaces. We compare the interfacial water structure–property relationship on three different carbon substrates, namely, amorphous carbon, compressed expanded natural graphite, and pure graphite by utilizing atomistic molecular dynamics simulations. The effect of different substrates on the structural and dynamical properties of water can readily be observed. Here, the density distributions parallel and normal to the substrates show oblate droplet structures. The normal to the substrate water distribution shows a strong hydration layer at the interface that does not vary with substrates. However, the disparity in the structure and dynamics on three different substrates shows that the surface morphologies of the substrates are critical for determining nanoscale water properties. Furthermore, it is observed that the formation of an interfacial water layer or the hydration layer is a direct consequence of both water “confinement” at the nanoscale and “attraction” between water molecules and the carbon substrates.
42 ENGINEERING↗
High-level coupled-cluster energetics by Monte Carlo sampling and moment expansions: Further details and comparisons
We recently proposed a novel approach to converging electronic energies equivalent to high-level coupled-cluster (CC) computations by combining the deterministic CC(P;Q) formalism with the stochastic configuration interaction (CI) and CC Quantum Monte Carlo (QMC) propagations. This article extends our initial study [J. E. Deustua, J. Shen, and P. Piecuch, Phys. Rev. Lett. 119, 223003 (2017)], which focused on recovering the energies obtained with the CC method with singles, doubles, and triples (CCSDT) using the information extracted from full CI QMC and CCSDT-MC, to the CIQMC approaches truncated at triples and quadruples. It also reports our first semi-stochastic CC(P;Q) calculations aimed at converging the energies that correspond to the CC method with singles, doubles, triples, and quadruples (CCSDTQ). The ability of the semi-stochastic CC(P;Q) formalism to recover the CCSDT and CCSDTQ energies, even when electronic quasi-degeneracies and triply and quadruply excited clusters become substantial, is illustrated by a few numerical examples, including the F–F bond breaking in F 2 , the automerization of cyclobutadiene, and the double dissociation of the water molecule.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Phonon- and defect-limited electron and hole mobility of diamond and cubic boron nitride: A critical comparison
Diamond and cBN are two of the most promising ultra-wide bandgap semiconductors for applications in high-power high-frequency electronic devices. Despite extensive studies on carrier transport in these materials, there are large discrepancies in their reported carrier mobilities. In this work, we investigate the phonon- and dopant-limited electron and hole mobilities of cBN and diamond with atomistic first-principles calculations in order to understand their fundamental upper bounds to carrier transport. Our results show that although the phonon-limited electron mobilities are comparable between cBN and diamond, the hole mobility is significantly lower in cBN due to its heavier hole effective mass. Moreover, although lattice scattering dominates the mobility at low doping, neutral impurity scattering becomes the dominant scattering mechanism at higher dopant concentrations due to the high dopant ionization energies. Our analysis provides critical insights and reveals the intrinsic upper limits to the carrier mobilities of diamond and cBN as a function of doping and temperature for applications in high-power electronic devices.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Designing high-performance superconductors with nanoparticle inclusions: Comparisons to strong pinning theory
One of the most promising routes for achieving high critical currents in superconductors is to incorporate dispersed, non-superconducting nanoparticles to control the dissipative motion of vortices. However, these inclusions reduce the overall superconducting volume and can strain the interlaying superconducting matrix, which can detrimentally reduce T c . Consequently, an optimal balance must be achieved between the nanoparticle density n p and size d. Determining this balance requires garnering a better understanding of vortex–nanoparticle interactions, described by strong pinning theory. Here, we map the dependence of the critical current on nanoparticle size and density in (Y 0.77 , Gd 0.23 )Ba 2 Cu 3 O 7-δ films in magnetic fields of up to 35 T and compare the trends to recent results from time-dependent Ginzburg–Landau simulations. We identify consistency between the field-dependent critical current J c (B) and expectations from strong pinning theory. Specifically, we find that J c ∝ B -α , where α decreases from 0.66 to 0.2 with increasing density of nanoparticles and increases roughly linearly with nanoparticle size d/ξ (normalized to the coherence length). At high fields, the critical current decays faster (~B -1 ), suggesting that each nanoparticle has captured a vortex. When nanoparticles capture more than one vortex, a small, high-field peak is expected in J c (B). Due to a spread in defect sizes, this novel peak effect remains unresolved here. Finally, we reveal that the dependence of the vortex creep rate S on nanoparticle size and density roughly mirrors that of α, and we compare our results to low-T nonlinearities in S(T) that are predicted by strong pinning theory.
36 MATERIALS SCIENCE↗
Comparison of private flux region instability in conventional and super-X divertor configurations
Understanding turbulence in the divertor leg of tokamaks is essential to predicting the heat deposition profile on the divertor plate. This in turn is important for evaluating advanced divertor configurations, such as the super-X divertor. Within the divertor region, the private flux region is of interest because it is relatively unaffected by turbulence extending from the outboard midplane, so instabilities in this region could have a particularly pronounced effect on transport. These instabilities are modeled using the Arbitrary Topology Equation Reader (ArbiTER) eigenvalue code. Eigenmodes are examined further by comparing physics models to determine the fundamental mechanisms behind their formation, and quantifying the effect of individual terms. This analysis is conducted on both conventional and super-X divertors to compare these effects. Here, the resulting analysis reveals the presence of a geodesic curvature driven instability that is significantly more pronounced in the super-X configuration.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
RF sheath induced sputtering on Proto-MPEX part 2: Impurity transport modeling and experimental comparison
The Prototype Material Plasma Exposure eXperiment (Proto-MPEX) is a pulsed linear plasma device that uses a radio frequency (RF) helicon antenna and a helicon ceramic window as the plasma ionization source. Modeling performed within this work on the helicon operations predicts hot spots on the helicon ceramic window due to an RF-induced sheath. The RF sheath potential leads to physical sputtering of the ceramic window predominantly where the helicon antenna sits around the helicon window. The role of oxygen in sputtering within Proto-MPEX is also investigated and is needed to match the experimental data. Experiments on Proto-MPEX show reasonable agreement with the modeled erosion/redeposition pattern on the helicon window, and the modeled impurity flux to the target shows radial transport within the device as expected from classical collisional impurity transport. Temperature screening of impurities is shown to be needed for Proto-MPEX to reduce impurity flux to the target. The modeling approach used in this work can be applied to other plasma RF applications for the determination of impurity production and transport.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Observation of a dipole-bound excited state in 4-ethynylphenoxide and comparison with the quadrupole-bound excited state in the isoelectronic 4-cyanophenoxide
Negative ions do not possess Rydberg states, but can have Rydberg-like nonvalence excited states near the electron detachment threshold, including dipole-bound states (DBSs) and quadrupole-bound states (QBSs). While DBSs have been studied extensively, quadrupole-bound excited states have been more rarely observed. 4-cyanophenoxide (4CP – ) was the first anion observed to possess a quadrupole-bound exited state 20 cm-1 below its detachment threshold. Here we report the observation of a DBS in the isoelectronic 4-ethynylphenoxide anion (4EP – ), providing a rare opportunity to compare the behaviors of a dipole-bound and a quadrupole-bound excited state in a pair of very similar anions. Photodetachment spectroscopy (PDS) of cryogenically-cooled 4EP – reveals a DBS 76 cm -1 below its detachment threshold. Photoelectron spectroscopy (PES) at 266 nm shows that the electronic structure of 4EP – and 4CP – are nearly identical. The observed vibrational features in both the PDS and PES, as well as autodetachment from the nonvalence excited states, are also found to be similar for both anions. However, resonant two-photon detachment (R2PD) from the bound vibrational ground state is observed to be very different for the DBS in 4EP – and the QBS in 4CP – . The R2PD spectra reveal that decays take place from both the DBS and QBS to the respective anion ground electronic states within the 5 ns detachment laser pulse due to internal conversion followed by intramolecular vibrational redistribution (IVR), but the decay mechanisms appear to be very different. In the R2PD spectrum of 4EP – , we observe strong threshold electron signals, which are due to detachment, by the second photon, of highly rotationally excited anions resulted from the decay of the DBS. On the other hand, in the R2PD spectrum of 4CP – we observe well-resolved vibrational peaks due to the three lowest-frequency vibrational modes of 4CP – , which are populated from the decay of the QBS. Here, the different behaviors of the R2PD spectra suggest unexpected differences between the relaxation mechanisms of the dipole-bound and quadrupole-bound excited states.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A comparison between the one- and two-step spin-orbit coupling approaches based on the ab initio Density Matrix Renormalization Group
The efficient and reliable treatment of both spin-orbit coupling (SOC) and electron correlation is essential for understanding f-element chemistry. Here, we analyze two approaches to the problem, the one-step approach where both effects are treated simultaneously, and the two-step state interaction approach. We report an implementation of the ab initio density matrix renormalization group(DMRG) with a one-step treatment of the SOC effect which can be compared to prior two-step treatments on an equal footing. Using a dysprosium octahedral complex and bridged dimer as benchmark systems, we identify characteristics of problems where the one-step approach is beneficial for obtaining the low-energy spectrum.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Model comparisons for two-temperature plasma equations of state
When a plasma is generated in the laboratory, energy is often deposited preferentially into either the electrons or the ions, giving rise to a quasiequilibrium state in which the two species, electrons and ions, are well described by two effective temperatures, T e and T i . Accurate hydrodynamic modeling of such a two-temperature plasma requires an equation of state that captures the relevant many-body physics without assuming a strict local thermodynamic equilibrium. Several models have been proposed within the literature, which extend conventional statistical approaches, each employing a different combination of assumptions for modifying the equilibrium equations. Here, in this work, we compare the predictions for several models, presenting derivations of the internal energy and pressure for each microscopic model within a unified framework so that the assumptions of each model may be more easily compared to one another. We find that for sufficiently weak coupling, all models agree with one another. However, as the coupling strength is increased, the disagreement between the models becomes more pronounced. Moreover, the relative sizes of the corrections predicted by each model depend on which species has the higher temperature, T e >T i vs T e < T i .
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
A comparison of past and present computational methods for shape analysis of double-shell x-ray radiographs
Implosion symmetry is a key requirement in achieving a robust burning plasma in inertial confinement fusion experiments. In double-shell capsule implosions, we are interested in the shape of the inner shell as it pushes on the fuel. Shape analysis is a popular technique for studying said symmetry during implosion. Combinations of filtering and contour-finding algorithms are studied for their promise in reliably recovering Legendre shape coefficients from synthetic radiographs of double-shell capsules with applied levels of noise. A radial lineout max(slope) method when used on an image pre-filtered with non-local means and a variant of the marching squares algorithm are able to recover p 0 , p 2 , and p 4 maxslope Legendre shape coefficients with mean pixel discrepancy errors of 2.81 and 3.06, respectively, for the noisy synthetic radiographs we consider. Here, this improves upon prior radial lineout methods paired with Gaussian filtering, which we show to be unreliable and whose performance is dependent on input parameters that are difficult to estimate.
47 OTHER INSTRUMENTATION↗