Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma processing”

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,477 records · Page 82

Combined Experimental/Numerical Study of the Soot Formation Process in a Gasoline Direct-Injection Spray in the Presence of Laser-Induced Plasma Ignition

Here, combustion issued from an eight-hole, direct-injection spray was experimentally studied in a constant-volume pre-burn combustion vessel using simultaneous high-speed diffused back-illumination extinction imaging (DBIEI) and OH* chemiluminescence. DBIEI has been employed to observe the liquid-phase of the spray and to quantitatively investigate the soot formation and oxidation taking place during combustion. The fuel-air mixture was ignited with a plasma induced by a single-shot Nd:YAG laser, permitting precise control of the ignition location in space and time. OH* chemiluminescence was used to track the high-temperature ignition and flame. The study showed that increasing the delay between the end of injection and ignition drastically reduces soot formation without necessarily compromising combustion efficiency. For long delays between the end of injection and ignition (1.9 ms) soot formation was eliminated in the main downstream charge of the fuel spray. However, poorly atomized and large droplets formed at the end of injection (dribble) eventually do form soot near the injector even when none is formed in the main charge. The quantitative soot measurements for these spray and ignition scenarios, resolved in time and space, represents a significant new achievement. Reynolds-averaged Navier-Stokes (RANS) simulations were performed to assess spray mixing and combustion. An analysis of the predicted fuel-air mixture in key regions, defined based upon experimental observations, was used to explain different flame propagation speeds and soot production tendencies when varying ignition timing. The mixture analysis indicates that soot production can be avoided if the flame propagates into regions where the equivalence ratio (Φ) is already below 2. Reactive RANS simulations have also been performed, but with a poor match against the experiment, as the flame speed and heat-release rate are largely over estimated. This modeling weakness appears related to a very high level of turbulent viscosity predicted for the high-momentum spray in the RANS simulations, which is an important consideration for modeling ignition and flame propagation in mixtures immediately created by the spray.

42 ENGINEERING↗

A thermal/nonthermal model for solar hard X-ray bursts

A formalism is developed for analyzing high-resolution hard X-ray spectra, incorporating the coexistence of thermal and nonthermal bremsstrahlung. The two processes are physically linked by the presence of electric currents, which both heat the surrounding plasma via Joule dissipation and accelerate electrons via the runaway process. We use this formalism to analyze the flare of 1980 June 27 and find that both the gradual and spike components of the hard X-ray emission are consistent with runaway acceleration. We also find that significant heating is observed only in the gradual component. The electric field is always sub-Dreicer, the maximum total potential drop in the acceleration region is found to be approximately 100 kV in two of the spikes, and the average accelerated electron flux is approximately 10(exp 34) electrons/s. We argue that classical resistivity is a valid assumption for this event and find the density in the current channels (approximately 10(exp 11)/cc) and a lower limit on the volume of the heated plasma. We find that the ratio of the electric field to the Dreicer field (epsilon = E/E(sub D) varies systematically, whereas the value of E alone does not. We also find that the acceleration region fragmented into many current/return current pairs, and that the fragmentation varied systematically. We also discuss further implications of this model.

Benka, Stephen G.↗

Post-inflationary dark matter bremsstrahlung

Dark matter may only interact with the visible sector efficiently at energy scales above the inflaton mass, such as the Planck scale or the grand unification scale. In such a scenario, the dark matter is mainly produced out of equilibrium during the period of reheating, often referred to as UV freeze-in. Here, we evaluate the abundance of the dark matter generated from bremsstrahlung off the inflaton decay products assuming no direct coupling between the inflaton and the dark matter. This process generally dominates the production of dark matter for low reheating temperatures where the production through the annihilations of particle in the thermal plasma becomes inefficient. We find that the bremsstrahlung process dominates for reheating temperatures T RH ≲ 10 10 GeV, and produces the requisite density of dark matter for a UV scale ≃ 10 16 GeV. As examples, we calculate numerically the yield of the dark matter bremsstrahlung through gravitation and dimension-6 vector portal effective interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Overview of physics results from MAST upgrade towards core-pedestal-exhaust integration

Recent results from MAST Upgrade are presented, emphasising understanding the capabilities of this new device and deepening understanding of key physics issues for the operation of ITER and the design of future fusion power plants. The impact of MHD instabilities on fast ion confinement have been studied, including the first observation of fast ion losses correlated with Compressional and Global Alfvén Eigenmodes. High-performance plasma scenarios have been developed by tailoring the early plasma current ramp phase to avoid internal reconnection events, resulting in a more monotonic q profile with low central shear. The impact of m/n = 3/2, 2/1 and 1/1 modes on thermal plasma confinement and rotation profiles has been quantified, and scenarios optimised to avoid them have transiently reached values of normalised beta approaching 4.2. In pedestal and ELM physics, a maximum pedestal top temperature of ~350 eV has been achieved, exceeding the value achieved on MAST at similar heating power. Mitigation of type-I ELMs with n = 1 RMPs has been observed. Studies of plasma exhaust have concentrated on comparing conventional and Super-X divertor configurations, while X-point target, X-divertor and snowflake configurations have been developed and studied in parallel. In L-mode discharges, the separatrix density required to detach the outer divertors is approximately a factor 2 lower in the Super-X than the conventional configuration, in agreement with simulations. Detailed analysis of spectroscopy data from studies of the Super-X configuration reveal the importance of including plasma-molecule interactions and D 2 Fulcher band emission to properly quantify the rates of ionisation, plasma-molecule interactions and volumetric recombination processes governing divertor detachment. In H-mode with conventional and Super-X configurations, the outer divertors are attached in the former and detached in the latter with no impact on core or pedestal confinement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Diffusion–convection model of runaway electrons due to large magnetohydrodynamic perturbations in post-thermal quench plasmas

Systematic test particle tracing simulations for runaway electrons (REs) are performed for six post-thermal quench equilibria from DIII-D and ITER, where large scale, kink-like n = 1 (n is the toroidal mode number) magnetohydrodynamic (MHD) instabilities are found. The modeled particle guiding center orbits allow extraction of the effective diffusion–convection coefficients of REs in the presence of large three-dimensional (3D) perturbations up to 10% of the equilibrium toroidal field. With a fixed spatial distribution of the field perturbation, the RE transport coefficients along the plasma radial coordinate track reasonably well with the surface-averaged perturbation level. A substantial variation in the value of the transport coefficients—by three orders of magnitude in most cases, however, occurs with varying launching location of REs along the plasma radius. Large 3D perturbations almost always lead to comparable diffusion and convection processes, meaning that diffusion alone is insufficient to describe the particle motion. At lower (but still high) level of perturbation, the RE convection is found to be dominant over diffusion. A similar observation is made when the perturbation is too strong. In the presence of large perturbation, the dependence of the RE transport on the particle energy is sensitive to the spatial distribution of the perturbation. Based on numerically obtained RE transport coefficients, an analytic fitting model is proposed to quantify the particle diffusion and convection processes due to large MHD events in post-thermal quench plasmas. The model is shown to reasonably well reproduce the direct test particle tracing results for the RE loss fraction and can, thus, be useful for incorporating into other kinetic RE codes in order to simulate the RE beam evolution in the presence of large 3D perturbations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SPACE: 3D parallel solvers for Vlasov-Maxwell and Vlasov-Poisson equations for relativistic plasmas with atomic transformations

A parallel, relativistic, three-dimensional particle-in-cell code SPACE has been developed for the simulation of electromagnetic fields, relativistic particle beams, and plasmas. In addition to the standard second-order Particle-in-Cell (PIC) algorithm, SPACE includes efficient novel algorithms to resolve atomic physics processes such as multi-level ionization of plasma atoms, recombination, and electron attachment to dopants in dense neutral gases. SPACE also contains a highly adaptive particle-based method, called Adaptive Particle-in-Cloud (AP-Cloud), for solving the Vlasov-Poisson problems. It eliminates the traditional Cartesian mesh of PIC and replaces it with an adaptive octree data structure. The code's algorithms, structure, capabilities, parallelization strategy, and performance have been discussed. Additionally, typical examples of SPACE applications to accelerator science and engineering problems are described.

43 PARTICLE ACCELERATORS↗

Direct Conversion of Ethane to Oxygenates, Ethylene, and Hydrogen in a Noncatalytic Biphasic Plasma Microreactor

We selectively upgrade ethane (C 2 H 6 ) to ethanol (C 2 H 5 OH), methanol (CH 3 OH), and acetic acid (CH 3 COOH) in a catalyst-free, continuous, argon/water biphasic plasma microreactor. The water (H 2 O) evaporates and electron- dissociates into OH· radicals. OH· recombines with alkyl radicals, produced via electron dissociation of ethane, to generate the oxygenates that absorb into H 2 O. A plasma-assisted path, reminiscent of the low-temperature thermocatalytic ethane steam reforming, leads to significant H 2 coproduction. The gaseous stream also comprises CO 2 and C 2 H 4 . Up to 1.3 and 1 μmol min –1 of liquid C 2 H 5 OH and CH 3 OH are attained, respectively. Compared to CO 2 -assisted ethane plasma conversion, which produces many oxygenates with low selectivity, the carbon selectivity can range from >70% C 2 H 5 OH, CH 3 OH, and CH 3 COOH to 60% C 2 H 4 . Furthermore, the low carbon footprint, electrified, modular, intensified process using a reactive evaporation and separation plasma could pave the way for the valorization of underutilized shale gas resources in remote areas.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Chemical and electrochemical pathways to low-carbon iron and steel

Currently, the iron and steel industry is responsible for 7% of global CO 2 emissions. In this review, we summarize the operational principles of current emissions-intensive steelmaking technologies and review emerging low- and zero-carbon technologies that could substantially reduce emissions. Current technologies that are discussed include blast furnaces, electric arc furnaces, and smelting. Promising low-carbon routes include use of alternative reductants for ore processing (hydrogen direct reduction, hydrogen plasma-smelting, hydrogen smelting, and ammonia-based reduction), electrolytic iron production (with aqueous and molten oxide electrolytes) and biocarbon-based electric arc furnace operation. Advantages of each approach are presented, and remaining research hurdles are identified.

36 MATERIALS SCIENCE↗

Investigations of the stability of GaAs for photoelectrochemical H 2 evolution in acidic or alkaline aqueous electrolytes

The long-term stability of p-GaAs photocathodes has been investigated for the hydrogen-evolution reaction (HER) in contact with either 1.0 M H 2 SO 4 (aq) or 1.0 M KOH(aq). Stability for the HER was evaluated using p-GaAs electrodes that were either etched or coated with active HER catalysts (Pt and CoP). Changes in surface characteristics of GaAs after exposure to electrochemical conditions were monitored by X-ray photoelectron spectroscopy (XPS), and electrode dissolution processes were evaluated by inductively coupled plasma mass spectrometry (ICP-MS). Consistent with thermodynamic predictions, after operation of the HER at pH 0 or pH 14, illuminated etched p-GaAs electrodes exhibited minimal dissolution while preserving a nearly stoichiometric surface. Electrodeposition or sputtering of Pt on the p-GaAs surface promoted the formation of excess As 0 via an interfacial reaction during the HER. The resulting non-stoichiometric As 0 -rich surface of p-GaAs/Pt electrodes caused a loss in photoactivity as well as substantial cathodic dark current. In contrast, p-GaAs electrodes coated with thin-film CoP catalysts did not display an increase in surficial As 0 after operation of the HER in acidic electrolytes. Minimization of deleterious interfacial reactions is thus critical to obtain extended stability in conjunction with high performance from p-GaAs photocathodes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Semi-Dynamic Leach Testing of Densified Silicon-based Iodine Waste Forms

Iodine waste forms (IWF) require a conceptual corrosion release model (CCRM) to provide iodine (I) release rates for performance modeling nuclear waste disposal repositories. To develop a CCRM, an understanding of the corrosion mechanisms of the IWF are required along with data from consistent test methods to parameterize the model. The present study has advanced both areas by providing and assessing the corrosion resistance of IWF types based on their processing history in minor variations of semi-dynamic leach tests. The test included a series of semi-dynamic leach tests using monolithic IWFs in deionized water (leachant). Several experiments were conducted under alternate test conditions with changes to temperature, leachant replacement, leachant pH, leachant volume, masking, and surface finish to elucidate if varying these conditions impacted IWF corrosion behavior. Tests were conducted on two classes of IWFs: (1) I-bearing silver-mordenite (AgZ) materials processed by hot isostatic pressing (HIP) at different temperatures, pressures, sizes, and times; and (2) I-bearing silver-functionalized silica aerogels (SFA) processed by either HIP or spark plasma sintering (SPS). The corrosion susceptibility of AgZ samples was influenced by HIP temperature and pressure. The SPS SFAs retained I far better than HIP SFAs. Additional findings in this study include: (1) The iodine dissolution rate decreased with decreasing temperature, (2) a common ion effect may occur and slow dissolution of the host phase if the leachant is not regularly replaced, (3) pH controls the dissolution rate, and (4) the iodine dissolution rate slows with extended test time (up to 224 days). Based on this work, these parameters should thus be represented when developing a CCRM.

corrosion↗

The Three Hundred Project: Modeling baryon and hot-gas fraction evolution in simulated clusters

The baryon fraction of galaxy clusters, expressed as the ratio between the mass in baryons (including both stars and cold or hot gas) and the total mass, is a powerful tool to provide information on the cosmological parameters, while the hot-gas fraction provides indications on the physics of the intracluster plasma and its interplay with the processes that drive galaxy formation. Using cosmological hydrodynamical simulations of about 300 simulated massive galaxy clusters with a median mass M 500 ≈ 7 × 10 14 M ⊙ at z = 0, we model the relations between total mass and either baryon fraction or the hot gas fractions at overdensities Δ = 2500, 500, and 200 with respect to the cosmic critical density, and their evolution from z ∼ 0 to z ∼ 1.3. We utilized the simulated galaxy clusters from the Three Hundred project, which include star formation and feedback from both supernovae and active galactic nuclei. We fit the simulation results for such scaling relations against three analytic forms (linear, quadratic, and logarithmic in a logarithmic plane) and three forms for the redshift dependence, and we considered as a variable both the inverse of the cosmic scale factor, (1 + z), and the Hubble expansion rate, E(z). We show that power-law dependencies on cluster mass poorly describe the investigated relations. A power law fails to simultaneously capture the flattening of the total baryon and gas fractions at high masses, their drop at low masses, and the transition between these two regimes. The other two functional forms provide a more accurate description of the curvature in mass scaling. The fractions measured within smaller radii exhibit a stronger evolution than those measured within larger radii. From the analysis of these simulations, we evince that as long as we include systems in the mass range herein investigated, the baryon or gas fraction can be accurately related to the total mass through either a parabola or a logarithm in the logarithmic plane. The trends are common to all modern hydro simulations, although the amplitude of the drop at low masses might differ. Being able to observationally determine the gas fraction in groups will thus provide constraints on the baryonic physics.

galaxy clusters↗

Gyrokinetic particle simulations of interactions between energetic particles and magnetic islands induced by neoclassical tearing modes

Interactions between energetic particles (EPs) and neoclassical tearing mode (NTM) islands in the DIII-D tokamak are studied using the global gyrokinetic toroidal code (GTC). GTC simulations find that the EP radial profile is partially flattened within the magnetic island regions and that there are stochastic regions in the particle phase space. Radial particle flux is induced mainly around the magnetic island regions and decreases with time to almost zero when the initial EP distribution achieves a new steady-state in the absence of EP sources. Stochastic regions of magnetic field lines induced by the superposition of multiple islands have weak effects on the particle flux when the width of stochastic regions is smaller than the EP drift orbit width. In conclusion, the perturbed parallel EP current induced by the magnetic islands has weak stabilizing effects on the linear growth rate of the NTM instability in this DIII-D experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation of Alfvénic avalanche onset in NSTX

In some tokamak discharges, a number of Alfvén eigenmodes (AEs) have been observed to cause a large-scale collapse of the high energy particle distribution, a phenomenon referred to as an avalanche. We examine the necessary conditions for an avalanche using the available experimental information from NSTX on the equilibrium and mode properties for two cases, one with a measurable but benign AE activity and one with an AE activity leading up to an avalanche. To produce an avalanche, the modes present in the discharge must possess resonances that can overlap with a modest increase in instability magnitude, providing a path to global particle stochastic motion. In this paper, we find that the modes present in the avalanche-free discharge do not provide such a path even at a very large amplitude. During the discharge which subsequently produces an avalanche, the high energy population is growing and the Alfvén frequency is dropping due to increasing density, and we find that both these changes, producing a small increased drive or an increased resonance width for the Alfvén modes, can lead in this case to uncontrolled mode growth and large-scale beam particle loss.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Photonic Doppler velocimetry of ohmically exploded aluminum surfaces

Photonic Doppler Velocimetry (PDV) was used to measure the exploding surface of coated metal driven by lineal current density increasing at 3 × 10 15 A/m/s. Here, aluminum-6061 rods were coated with three thicknesses (5, 17, and 41 μm) of Parylene-N and driven to 850 kA in approximately 100 ns, with the metal thicker than the skin depth. According to PDV, the metal surface melts at 85 ± 13 T and expands with a constant acceleration. Then, at surface magnetic field 140-180 T, the metal acceleration begins to increase in time. Later, plasma forms on, and is correlated with PDV signal loss from, uncoated aluminum rods, 5-μm-coated rods, and 17-μm-coated rods, while no evidence suggests plasma forms on 41-μm-coated metal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Gyrokinetic understanding of the edge pedestal transport driven by resonant magnetic perturbations in a realistic divertor geometry

Self-consistent simulations of neoclassical and electrostatic turbulent transport in a DIII-D H-mode edge plasma under resonant magnetic perturbations (RMPs) have been performed using the global total-f gyrokinetic particle-in-cell code x-point gyrokinetic code (XGC), in order to study density pump-out and electron heat confinement. The RMP field is imported from the extended magneto-hydrodynamics code M3D-C1, taking into account the linear two-fluid plasma response. With both neoclassical and turbulence physics considered together, the XGC simulation reproduces two key features of experimentally observed edge transport under RMPs: increased radial particle transport in the pedestal region that is sufficient to account for the experimental pump-out rate and suppression of the electron heat flux in the steepest part of the edge pedestal. In the simulation, the density fluctuation amplitude of modes moving in the electron diamagnetic direction increases due to interaction with RMPs in the pedestal shoulder and outward, while the electron temperature fluctuation amplitude decreases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Phase-space dynamics of Alfvén mode chirping

In tokamak discharges, toroidal Alfvén eigenmodes often experience complex semi-periodic frequency modulation known as chirping. These events modify the local high energy particle distribution and are expected to occur in many future fusion devices, which include energetic beams or fusion products. In this paper, we present a study of simulations of mode chirping made in order to better understand its phase-space properties in a realistic tokamak configuration. We find a mechanism that permits rapid repeated chirping with strong amplitude variation in each chirp. Each chirp is then associated with an amplitude crash to low magnitude and local manipulation of the density gradients through a shift of mode phase through π. The chirping produces high density clumps, which propagate down the fast ion density gradient and low density holes that propagate up the density gradient away from the resonance. This flow of particles across the resonance provides an energy source and local gradients for repeated chirping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A generalized forward fit for neutron detectors with energy-dependent response functions

To date, most analysis of neutron time-of-flight data from inertial confinement fusion experiments has focused on the relatively small range of energies corresponding to the primary neutrons from DD and DT fusion, and have therefore employed instrument response functions (IRF’s) corresponding to monoenergetic 2.45-MeV or 14.03-MeV neutrons. For analysis of time-of-flight signals corresponding to broader ranges of neutron energies, accurate treatment of the data requires the use of an energy-dependent IRF. Here, this work describes interpolation of the IRF for neutrons of arbitrary energy, construction of an energy-dependent IRF, and application of this IRF in a forward fit via matrix multiplication. As an example of the application of this method, an analysis of synthetic data relevant to TT fusion experiments at the Omega Laser Facility is discussed. This example is used to illustrate the differences between a forward fit that uses an energy-dependent IRF and a forward fit that uses a monoenergetic IRF. Use of the energy-dependent IRF is shown to result in accurate inference of the fit parameters of interest.

47 OTHER INSTRUMENTATION↗