Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Agree”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Reaction dynamics of O({sup 1}D) + HCOOD/DCOOH investigated with time-resolved Fourier-transform infrared emission spectroscopy

We investigated the reaction dynamics of O({sup 1}D) towards hydrogen atoms of two types in HCOOH. The reaction was initiated on irradiation of a flowing mixture of O{sub 3} and HCOOD or DCOOH at 248 nm. The relative vibration-rotational populations of OH and OD (1 ≦ v ≦ 4, J ≤ 15) states were determined from time-resolved IR emission recorded with a step-scan Fourier-transform spectrometer. In the reaction of O({sup 1}D) + HCOOD, the rotational distribution of product OH is nearly Boltzmann, whereas that of OD is bimodal. The product ratio [OH]/[OD] is 0.16 ± 0.05. In the reaction of O({sup 1}D) + DCOOH, the rotational distribution of product OH is bimodal, but the observed OD lines are too weak to provide reliable intensities. The three observed OH/OD channels agree with three major channels of production predicted with quantum-chemical calculations. In the case of O({sup 1}D) + HCOOD, two intermediates HOC(O)OD and HC(O)OOD are produced in the initial C−H and O−D insertion, respectively. The former undergoes further decomposition of the newly formed OH or the original OD, whereas the latter produces OD via direct decomposition. Decomposition of HOC(O)OD produced OH and OD with similar vibrational excitation, indicating efficient intramolecular vibrational relaxation, IVR. Decomposition of HC(O)OOD produced OD with greater rotational excitation. The predicted [OH]/[OD] ratio is 0.20 for O({sup 1}D) + HCOOD and 4.08 for O({sup 1}D) + DCOOH; the former agrees satisfactorily with experiments. We also observed the v{sub 3} emission from the product CO{sub 2}. This emission band is deconvoluted into two components corresponding to internal energies E = 317 and 96 kJ mol{sup −1} of CO{sub 2}, predicted to be produced via direct dehydration of HOC(O)OH and secondary decomposition of HC(O)O that was produced via decomposition of HC(O)OOH, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Response of the plasma to the size of an anode electrode biased near the plasma potential

As the size of a positively biased electrode increases, the nature of the interface formed between the electrode and the host plasma undergoes a transition from an electron-rich structure (electron sheath) to an intermediate structure containing both ion and electron rich regions (double layer) and ultimately forms an electron-depleted structure (ion sheath). In this study, measurements are performed to further test how the size of an electron-collecting electrode impacts the plasma discharge the electrode is immersed in. This is accomplished using a segmented disk electrode in which individual segments are individually biased to change the effective surface area of the anode. Measurements of bulk plasma parameters such as the collected current density, plasma potential, electron density, electron temperature and optical emission are made as both the size and the bias placed on the electrode are varied. Abrupt transitions in the plasma parameters resulting from changing the electrode surface area are identified in both argon and helium discharges and are compared to the interface transitions predicted by global current balance [S. D. Baalrud, N. Hershkowitz, and B. Longmier, Phys. Plasmas 14, 042109 (2007)]. While the size-dependent transitions in argon agree, the size-dependent transitions observed in helium systematically occur at lower electrode sizes than those nominally derived from prediction. The discrepancy in helium is anticipated to be caused by the finite size of the interface that increases the effective area offered to the plasma for electron loss to the electrode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Note: Calibration of EBT3 radiochromic film for measuring solar ultraviolet radiation

Solar (UVA + UVB) exposure was assessed using the Gafchromic EBT3 film. The coloration change was represented by the net reflective optical density (Net ROD). Through calibrations against a UV-tube lamp, operational relationships were obtained between Net ROD and the (UVA + UVB) exposures (in J cm⁻²p or J m⁻²). The useful range was from ~0.2 to ~30 J cm⁻². The uniformity of UV irradiation was crucial for an accurate calibration. For solar exposures ranging from 2 to 11 J cm⁻², the predicted Net ROD agreed with the recorded values within 9%, while the predicted exposures agreed with the recorded values within 15%.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Lateral migration of a microdroplet under optical forces in a uniform flow

The behavior of a microdroplet in a uniform flow and subjected to a vertical optical force applied by a loosely focused Gaussian laser beam was studied numerically. The lattice Boltzmann method was applied to obtain the two-phase flow field, and the dynamic ray tracing method was adopted to calculate the optical force. The optical forces acting on the spherical droplets agreed well with the analytical values. The numerically predicted droplet migration distances agreed well with the experimentally obtained values. Simulations of the various flow and optical parameters showed that the droplet migration distance nondimensionalized by the droplet radius is proportional to the S number (z{sub d}/r{sub p} = 0.377S), which is the ratio of the optical force to the viscous drag. The effect of the surface tension was also examined. These results indicated that the surface tension influenced the droplet migration distance to a lesser degree than the flow and optical parameters. The results of the present work hold for the refractive indices of the mean fluid and the droplet being 1.33 and 1.59, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Smoothed particle hydrodynamics model for Landau-Lifshitz-Navier-Stokes and advection-diffusion equations

We propose a novel smoothed particle hydrodynamics (SPH) discretization of the fully coupled Landau-Lifshitz-Navier-Stokes (LLNS) and stochastic advection-diffusion equations. The accuracy of the SPH solution of the LLNS equations is demonstrated by comparing the scaling of velocity variance and the self-diffusion coefficient with kinetic temperature and particle mass obtained from the SPH simulations and analytical solutions. The spatial covariance of pressure and velocity fluctuations is found to be in a good agreement with theoretical models. To validate the accuracy of the SPH method for coupled LLNS and advection-diffusion equations, we simulate the interface between two miscible fluids. We study formation of the so-called “giant fluctuations” of the front between light and heavy fluids with and without gravity, where the light fluid lies on the top of the heavy fluid. We find that the power spectra of the simulated concentration field are in good agreement with the experiments and analytical solutions. In the absence of gravity, the power spectra decay as the power −4 of the wavenumber—except for small wavenumbers that diverge from this power law behavior due to the effect of finite domain size. Gravity suppresses the fluctuations, resulting in much weaker dependence of the power spectra on the wavenumber. Finally, the model is used to study the effect of thermal fluctuation on the Rayleigh-Taylor instability, an unstable dynamics of the front between a heavy fluid overlaying a light fluid. The front dynamics is shown to agree well with the analytical solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Analysis of shockless dynamic compression data on solids to multi-megabar pressures: Application to tantalum

Magnetically-driven, planar shockless-compression experiments to multi-megabar pressures were performed on tantalum samples using a stripline target geometry. Free-surface velocity waveforms were measured in 15 cases; nine of these in a dual-sample configuration with two samples of different thicknesses on opposing electrodes, and six in a single-sample configuration with a bare electrode opposite the sample. Details are given on the application of inverse Lagrangian analysis (ILA) to these data, including potential sources of error. The most significant source of systematic error, particularly for single-sample experiments, was found to arise from the pulse-shape dependent free-surface reflected wave interactions with the deviatoric-stress response of tantalum. This could cause local, possibly temporary, unloading of material from a ramp compressed state, and thus multi-value response in wave speed that invalidates the free-surface to in-material velocity mapping step of ILA. By averaging all 15 data sets, a final result for the principal quasi-isentrope of tantalum in stress-strain was obtained to a peak longitudinal stress of 330 GPa with conservative uncertainty bounds of ±4.5% in stress. The result agrees well with a tabular equation of state developed at Los Alamos National Laboratory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Graphics processing unit accelerated three-dimensional model for the simulation of pulsed low-temperature plasmas

A 3-dimensional particle-in-cell/Monte Carlo collision simulation that is fully implemented on a graphics processing unit (GPU) is described and used to determine low-temperature plasma characteristics at high reduced electric field, E/n, in nitrogen gas. Details of implementation on the GPU using the NVIDIA Compute Unified Device Architecture framework are discussed with respect to efficient code execution. The software is capable of tracking around 10 × 10{sup 6} particles with dynamic weighting and a total mesh size larger than 10{sup 8} cells. Verification of the simulation is performed by comparing the electron energy distribution function and plasma transport parameters to known Boltzmann Equation (BE) solvers. Under the assumption of a uniform electric field and neglecting the build-up of positive ion space charge, the simulation agrees well with the BE solvers. The model is utilized to calculate plasma characteristics of a pulsed, parallel plate discharge. A photoionization model provides the simulation with additional electrons after the initial seeded electron density has drifted towards the anode. Comparison of the performance benefits between the GPU-implementation versus a CPU-implementation is considered, and a speed-up factor of 13 for a 3D relaxation Poisson solver is obtained. Furthermore, a factor 60 speed-up is realized for parallelization of the electron processes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Orbital-motion-limited theory of dust charging and plasma response

The foundational theory for dusty plasmas is the dust charging theory that provides the dust potential and charge arising from the dust interaction with a plasma. The most widely used dust charging theory for negatively charged dust particles is the so-called orbital motion limited (OML) theory, which predicts the dust potential and heat collection accurately for a variety of applications, but was previously found to be incapable of evaluating the dust charge and plasma response in any situation. Here, we report a revised OML formulation that is able to predict the plasma response and hence the dust charge. Numerical solutions of the new OML model show that the widely used Whipple approximation of dust charge-potential relationship agrees with OML theory in the limit of small dust radius compared with plasma Debye length, but incurs large (order-unity) deviation from the OML prediction when the dust size becomes comparable with or larger than plasma Debye length. This latter case is expected for the important application of dust particles in a tokamak plasma.

Tang, Xian-Zhu↗
DOE OSTI2014

Verification of gyrokinetic particle simulation of current-driven instability in fusion plasmas. II. Resistive tearing mode

In this work, global gyrokinetic particle simulation of resistive tearing modes has been developed and verified in the gyrokinetic toroidal code (GTC). GTC linear simulations in the fluid limit of the kink-tearing and resistive tearing modes in the cylindrical geometry agree well with the resistive magnetohydrodynamic eigenvalue and initial value codes. Ion kinetic effects are discovered to reduce the radial width of the tearing modes. GTC simulations of the resistive tearing modes in the toroidal geometry find that the toroidicity reduces the growth rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Plasma response measurements of non-axisymmetric magnetic perturbations on DIII-D via soft x-ray imaging

Recent observations on DIII-D have advanced the understanding of plasma response to applied resonant magnetic perturbations (RMPs) in both H-mode and L-mode plasmas. Three distinct 3D features localized in minor radius are imaged via filtered soft x-ray emission: (i) the formation of lobes extending from the unperturbed separatrix in the X-point region at the plasma boundary, (ii) helical kink-like perturbations in the steep-gradient region inside the separatrix, and (iii) amplified islands in the core of a low-rotation L-mode plasma. In this study, these measurements are used to test and to validate plasma response models, which are crucial for providing predictive capability of edge-localized mode control. In particular, vacuum and two-fluid resistive magnetohydrodynamic(MHD) responses are tested in the regions of these measurements. At the plasma boundary in H-mode discharges with n = 3 RMPs applied, measurements compare well to vacuum-field calculations that predict lobe structures. Yet in the steep-gradient region, measurements agree better with calculations from the linear resistive two-fluid MHD code, M3D-C1. Relative to the vacuum fields, the resistive two-fluid MHD calculations show a reduction in the pitch-resonant components of the normal magnetic field (screening), and amplification of non-resonant components associated with ideal kink modes. However, the calculations still over-predict the amplitude of the measuredperturbation by a factor of 4. In a slowly rotating L-mode plasma with n = 1 RMPs, core islands are observed amplified from vacuum predictions. Finally, these results indicate that while the vacuum approach describes measurements in the edge region well, it is important to include effects of extended MHD in the pedestal and deeper in the plasma core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Compact tunable Compton x-ray source from laser-plasma accelerator and plasma mirror

We present an in-depth experimental-computational study of the parameters necessary to optimize a tunable, quasi-monoenergetic, efficient, low-background Compton backscattering (CBS) x-ray source that is based on the self-aligned combination of a laser-plasma accelerator (LPA) and a plasma mirror (PM). The main findings are (1) an LPA driven in the blowout regime by 30 TW, 30 fs laser pulses produce not only a high-quality, tunable, quasi-monoenergetic electron beam, but also a high-quality, relativistically intense (a{sub 0} ∼ 1) spent drive pulse that remains stable in profile and intensity over the LPA tuning range. (2) A thin plastic film near the gas jet exit retro-reflects the spent drive pulse efficiently into oncoming electrons to produce CBS x-rays without detectable bremsstrahlung background. Meanwhile, anomalous far-field divergence of the retro-reflected light demonstrates relativistic “denting” of the PM. Exploiting these optimized LPA and PM conditions, we demonstrate quasi-monoenergetic (50% FWHM energy spread), tunable (75–200 KeV) CBS x-rays, characteristics previously achieved only on more powerful laser systems by CBS of a split-off, counter-propagating pulse. Moreover, laser-to-x-ray photon conversion efficiency (∼6 × 10{sup −12}) exceeds that of any previous LPA-based quasi-monoenergetic Compton source. Particle-in-cell simulations agree well with the measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Surface hopping with a manifold of electronic states. II. Application to the many-body Anderson-Holstein model

We investigate a simple surface hopping (SH) approach for modeling a single impurity level coupled to a single phonon and an electronic (metal) bath (i.e., the Anderson-Holstein model). The phonon degree of freedom is treated classically with motion along–and hops between–diabatic potential energy surfaces. The hopping rate is determined by the dynamics of the electronic bath (which are treated implicitly). For the case of one electronic bath, in the limit of small coupling to the bath, SH recovers phonon relaxation to thermal equilibrium and yields the correct impurity electron population (as compared with numerical renormalization group). For the case of out of equilibrium dynamics, SH current-voltage (I-V) curve is compared with the quantum master equation (QME) over a range of parameters, spanning the quantum region to the classical region. In the limit of large temperature, SH and QME agree. Furthermore, we can show that, in the limit of low temperature, the QME agrees with real-time path integral calculations. As such, the simple procedure described here should be useful in many other contexts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Properties of toroidal Alfvén eigenmode in DIII-D plasma

Here, linear properties of the toroidal Alfvén eigenmode (TAE) excited by energetic particles (EP) in a DIII-D tokamak experiment have been examined in global gyrokinetic particle simulations treating self-consistently kinetic effects of EP, thermal ions, and electrons. Simulation results of the TAE frequency and mode structure agree very well with the experimental measurements. The non-perturbative EP contribution induces a radial localization of the TAE mode structure, a break-down of mode radial symmetry, as well as a frequency dependence on the toroidal mode number. The simulations further exhibit the dependence of the growth rate and mode structure on EP pressure gradients. The in-out asymmetry of the mode structure and the experimental identification of the poloidal harmonics have also been observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Simulation of injector dynamics during steady inductive helicity injection current drive in the HIT-SI experiment

We present simulations of inductive helicity injection in the Helicity Injected Torus with Steady Inductive helicity injection (HIT-SI) device that treats the entire plasma volume in a single dynamic MHD model. A new fully 3D numerical tool, the PSI-center TETrahedral mesh code, was developed that provides the geometric flexibility required for this investigation. Implementation of a zero-β Hall MHD model using PSI-TET will be presented including formulation of a new self-consistent magnetic boundary condition for the wall of the HIT-SI device. Results from simulations of HIT-SI are presented focusing on injector dynamics that are investigated numerically for the first time. Asymmetries in the plasma loading between the two helicity injectors and progression of field reversal in each injector are observed. Analysis indicates cross-coupling between injectors through confinement volume structures. Injector impedance is found to scale with toroidal current at fixed density, consistent with experimental observation. Comparison to experimental data with an injector drive frequency of 14.5 kHz shows good agreement with magnetic diagnostics. Global mode structures from Bi-Orthogonal decomposition agree well with experimental data for the first four modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Structural stability and defect energetics of ZnO from diffusion quantum Monte Carlo

We have applied the many-body ab initio diffusion quantum Monte Carlo (DMC) method to study Zn and ZnO crystals under pressure and the energetics of the oxygen vacancy, zinc interstitial, and hydrogen impurities in ZnO. We show that DMC is an accurate and practical method that can be used to characterize multiple properties of materials that are challenging for density functional theory (DFT) approximations. DMC agrees with experimental measurements to within 0.3 eV, including the band-gap of ZnO, the ionization potential of O and Zn, and the atomization energy of O{sub 2}, ZnO dimer, and wurtzite ZnO. DMC predicts the oxygen vacancy as a deep donor with a formation energy of 5.0(2) eV under O-rich conditions and thermodynamic transition levels located between 1.8 and 2.5 eV from the valence band maximum. Our DMC results indicate that the concentration of zinc interstitial and hydrogen impurities in ZnO should be low under n-type and Zn- and H-rich conditions because these defects have formation energies above 1.4 eV under these conditions. Comparison of DMC and hybrid functionals shows that these DFT approximations can be parameterized to yield a general correct qualitative description of ZnO. However, the formation energy of defects in ZnO evaluated with DMC and hybrid functionals can differ by more than 0.5 eV.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Numerical studies and metric development for validation of magnetohydrodynamic models on the HIT-SI experiment

We present application of three scalar metrics derived from the Biorthogonal Decomposition (BD) technique to evaluate the level of agreement between macroscopic plasma dynamics in different data sets. BD decomposes large data sets, as produced by distributed diagnostic arrays, into principal mode structures without assumptions on spatial or temporal structure. These metrics have been applied to validation of the Hall-MHD model using experimental data from the Helicity Injected Torus with Steady Inductive helicity injection experiment. Each metric provides a measure of correlation between mode structures extracted from experimental data and simulations for an array of 192 surface-mounted magnetic probes. Numerical validation studies have been performed using the NIMROD code, where the injectors are modeled as boundary conditions on the flux conserver, and the PSI-TET code, where the entire plasma volume is treated. Initial results from a comprehensive validation study of high performance operation with different injector frequencies are presented, illustrating application of the BD method. Using a simplified (constant, uniform density and temperature) Hall-MHD model, simulation results agree with experimental observation for two of the three defined metrics when the injectors are driven with a frequency of 14.5 kHz.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

First principles-based multiparadigm, multiscale strategy for simulating complex materials processes with applications to amorphous SiC films

Progress has recently been made in developing reactive force fields to describe chemical reactions in systems too large for quantum mechanical (QM) methods. In particular, ReaxFF, a force field with parameters that are obtained solely from fitting QM reaction data, has been used to predict structures and properties of many materials. Important applications require, however, determination of the final structures produced by such complex processes as chemical vapor deposition, atomic layer deposition, and formation of ceramic films by pyrolysis of polymers. This requires the force field to properly describe the formation of other products of the process, in addition to yielding the final structure of the material. We describe a strategy for accomplishing this and present an example of its use for forming amorphous SiC films that have a wide variety of applications. Extensive reactive molecular dynamics (MD) simulations have been carried out to simulate the pyrolysis of hydridopolycarbosilane. The reaction products all agree with the experimental data. After removing the reaction products, the system is cooled down to room temperature at which it produces amorphous SiC film, for which the computed radial distribution function, x-ray diffraction pattern, and the equation of state describing the three main SiC polytypes agree with the data and with the QM calculations. Extensive MD simulations have also been carried out to compute other structural properties, as well the effective diffusivities of light gases in the amorphous SiC film.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗