Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electrical fluids”

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 145 records · Page 8

Electrification of Boilers in U.S. Manufacturing

Decarbonization of the industrial heat demand through electrification where low/no-carbon electricity is used can contribute significantly to global greenhouse gas (GHG) reduction. In U.S. manufacturing, thermal processes account for approximately 75% of the total final energy demand, of which nearly 17% was consumed by conventional industrial boilers for steam generation in 2018 (this does not include boilers for combined heat and power – CHP). Steam is generally used in industry to regulate temperatures and pressures in industrial processes, dry products, strip impurities from process fluids, etc. Although all kinds of energy sources such as fossil fuels, renewables, nuclear, and electricity can generate steam, fossil fuels’-fired boilers are dominant in U.S. manufacturing. Electric boilers, which are a mature technology, have a small market share for steam generation in the global and U.S. industry (approximately 2% in U.S. manufacturing) due to several techno-economic reasons.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Pumped Thermal Electricity Storage with Supercritical CO2 Cycles and Solar Heat Input

Pumped Thermal Electricity Storage (PTES) is an energy storage device that uses grid electricity to drive a heat pump that generates hot and cold storage reservoirs. This thermal potential is later used to power a heat engine and return electricity to the grid. In this article, a PTES variant that uses supercritical carbon dioxide (sCO2) as the working fluid is introduced. sCO2-PTES cycles have higher work ratios and power densities than the systems based on ideal gases that have been investigated to date. Furthermore, sCO2- PTES cycles may achieve higher round-trip efficiencies for a given hot storage temperature (up to 78% at 560°C). The sensitivity of PTES cycles to loss factors such as isentropic efficiencies and temperature differences between the power cycle and storage fluid is investigated. A second concept whereby an sCO2- PTES cycle is integrated with concentrating solar power (CSP) is introduced. This concept ‘time-shifts’ the recompression of an sCO2 recompression cycle to a period of lower electricity prices and stores the heat. When solar heat is dispatched, the recompressor may be avoided as the required heat is obtained from storage, thereby leading to increased heat engine efficiencies. The net work output of this integrated system is 10-18% greater than the conventional recompression cycle. Combining PTES with a CSP power cycle is therefore shown to improve the dispatch of solar heat as well as providing electricity storage services.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Townsend to glow discharge transition for a nanosecond pulse plasma in helium: space charge formation and resulting electric field dynamics

Stark polarization spectroscopy is used to investigate the temporal evolution of the electric field distribution in the cathode region of a nanosecond pulsed discharge in helium at 120 Torr. The measurements are performed on the He I transition at 492.19 nm, during the early stages of the discharge formation. The experimental results are compared with the predictions of a 1D fluid model. Time-resolved ICCD images show that the discharge develops as a diffuse, cathode-directed ionization wave with a Townsend-like feature before transitioning into a glow-like structure. Near anode instabilities characterized by filament formation were observed near the high voltage electrode. Within 30 ns, a reduction of the sheath thickness to about 250 μm is observed, coinciding with a gradual increase of the discharge current and proportional increase in electric field at the cathode. Additionally, the cathode electric field corresponding to this sheath with a thickness of 250 μm is about 40 kV cm -1 . A subsequent steep increase of the discharge current leads to a further reduction of the sheath width. The electric field evolution as obtained by the fluid model is in excellent agreement with the measurements and shows that an enhanced ionization near the cathode is causing the space charge formation responsible for the increase in electric field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electric field driven assembly of ordered nanocrystal superlattices

An electric field drives nanocrystals dispersed in solvents to assemble into ordered three-dimensional superlattices. A first electrode and a second electrode 214 are in the vessel. The electrodes face each other. A fluid containing charged nanocrystals fills the vessel between the electrodes. The electrodes are connected to a voltage supply which produces an electrical field between the electrodes. The nanocrystals will migrate toward one of the electrodes and accumulate on the electrode producing ordered nanocrystal accumulation that will provide a superlattice thin film, isolated superlattice islands, or coalesced superlattice islands.

Yu, Yixuan↗

A compact design for an electric driven hydraulic gear machine capable of multiple quadrant operation

This paper presents an innovative design of an integrated electro-hydraulic (EH) flow generation unit able to operate in multiple quadrants. The proposed unit is composed of a hydraulic machine which can operate in either pumping or motoring mode, and an electric machine, which can operate as electric motor or generator. The adopted morphological integration of the two unit fulfills the goals of maximizing compactness, power to weight ratio, energy efficiency and quietness. Furthermore, the design results from a numerical optimization procedure that includes detailed simulation of both the electric machine and the hydraulic machine, using simulation tools previously developed by the authors, but combined for the first time in this work. After describing the optimization procedure, the paper focuses on the design integration aspects of the novel unit. A 9 kW prototype of the integrated unit was also fabricated and tested to verify its performance. Tests results show a measured overall energy efficiency up to 69%, demonstrating the potentials of the proposed concept, despite some limitations present in the first prototypes which could be addressed in future design refinements.

42 ENGINEERING↗

Numerical Investigation of Enhanced Dehumidification Processes By Using Dielectrophoresis Principles in Moist Airflows

Dispersed particle-laden flows are encountered in many building and industrial applications, such as flow in a fluidized bed, hydrocarbon transportation in pipelines, and the fouling of air-cooled heat exchangers (Kuruneru et al., 2016; Ray et al., 2019; Wang et al., 2019). Computational fluid dynamic (CFD) models have been developed in recent years to depict particle-fluid and particle-particle interactions in laminar or turbulent flows with increasing accuracy and stability. One particular particle-laden system of interest for moisture control is electrically-enhanced condensation in air and water droplet flows. Electrically-enhanced condensation consists of the use of highly charged water droplets injected in the moist air. The droplets become electric seeds that attract polar water vapor molecules to their surfaces and promote condensation. The nucleation and growth of the charged droplets deplete the vapor phase near a droplet, which is compensated for by the dielectrophoresis flow and diffusion. Dielectrophoresis flow involves surrounding vapor at a distance of about 10 to 100 nm for droplets charged by an electrospray compared to ~2 nm for a single electron charge in a droplet. As the vapor molecules collapse on the surface of the droplets, their initial electrical charge decreases with time due to the neutralization of the ions. While the physics of this phenomena is well known, engineering models for predicting the condensation rates are not available. This work computationally investigates dehumidification of moist airflow in a converging rectangular duct. The objective is to develop an engineering model that predicts water vapor condensation by employing dielectrophoresis principles. We construct a Computational Fluid Dynamics (CFD) model of the duct with electrically-enhanced condensation. The model is implemented in the open-source software OpenFOAM. We utilize the Multi-Phase Particle-In-Cell (MP-PIC) method coupled with a Population Balance Equation (PBE) approach to simulate the particle-laden system. This methodology is an Eulerian-Lagrangian approach used to simulate the droplets' behavior in the humid air. The MP-PIC approach (Andrews and O'Rourke, 1996) mitigates the computational cost by parceling several fundamental particles with similar properties (such as types, sizes, and temperature) into one computational particle. Thus, the billions of particles can be substituted by millions of computational particles without significant loss of information. The PBE was considered with the Lagrangian frame to combine the particle distribution function used in MP-PIC (Kim et al., 2020). This approach preserves mass and energy conservation between the phases in the Eulerian and Lagrangian structures. The PBE in this procedure was directly linked to the discrete parcels, making the simulation of the particle distribution computationally efficient and robust. The MP-PIC-PBE approach used in the present work was applied to the dehumidification of air. Water droplets were injected in the air stream and forced to grow according to experimentally derived correlation. The experiments were conducted on a converging duct with the same geometry and boundary conditions used to build the CFD model. This approach enabled us to approximate the effect of dielectrophoresis phenomena on the droplet and air interface. This presentation will discuss the details of the new CFD model built for the duct, the implementation of the model in OpenFOAM CFD programming language, and the experimental validation of the newly developed model. The results revealed a moderate yet measurable increase in droplet diameter due to water vapor condensation at the vapor-liquid interface of the electrically charged droplets' surface. The seed water droplet particles grew in size by capturing the water vapor in the surrounding air. The OpenFOAM model predicted reductions of humidity in the air from 5 to 10 percent.

Yel Mahi, Maliha↗

Unveiling the Significance of Correlations in K-Space and Configuration Space for Drift Wave Turbulence in Tokamaks

Turbulence and transport phenomena play a crucial role in the confinement and stability of tokamak plasmas. Turbulent fluctuations in certain physical quantities, such as density or temperature fluctuations, can have a wide range of spatial scales, and understanding their correlation length is important for predicting and controlling the behavior of the plasma. The correlation length in the radial direction is identified as the critical length in real space. The dynamics in real space are of significant interest because transport in configuration space is primarily focused on them. When investigating transport caused by the E×B drift, the correlation length in real space represents the size of E×B whirls. It was numerically discovered that in drift wave turbulence, this length is inversely proportional to the normalized mode number of the fastest growing mode relative to the drift frequency. Considerable time was required before a proper analytical derivation of this condition was accomplished. Therefore, a connection has been established between phenomena occurring in real space and those occurring in k-space. Although accompanied by a turbulent spectrum in k-space with a substantial width, transport in real space is uniquely determined by the correlation length, allowing for accurate transport calculations through the dynamics of a single mode. Naturally, the dynamics are subject to nonlinear effects, with resonance broadening in frequency being the most significant nonlinear effect. Thus, mode number space is once again involved. Resonance broadening leads to the detuning of waves from particles, permitting a fluid treatment. It should be emphasized that the consideration here involves the total electric field, including the induction part, which becomes particularly important at higher beta plasmas.

Weiland, Jan (ORCID:0000000253496025)↗

Method for producing nanoparticles and the nanoparticles produced therefrom

Disclosed herein is a method comprising disposing a container containing a metal and/or ferromagnetic solid and abrasive particles in a static magnetic field; where the container is surrounded by an induction coil; activating the induction coil with an electrical current, to heat up the metallic or ferromagnetic solid to form a fluid; generating sonic energy to produce acoustic cavitation and abrasion between the abrasive particles and the container; and producing nanoparticles that comprise elements from the container, the metal and/or the ferromagnetic solid and the abrasive particles. Disclosed herein too is a composition comprising first metal or a first ceramic; and particles comprising carbides and/or nitrides dispersed therein. Disclosed herein too is a composition comprising nanoparticles comprising chromium carbide, iron carbide, nickel carbide, γ-Fe and magnesium nickel.

36 MATERIALS SCIENCE↗

Method for producing nanoparticles and the nanoparticles produced therefrom

Disclosed herein is a method comprising disposing a container containing a metal and/or ferromagnetic solid and abrasive particles in a static magnetic field; where the container is surrounded by an induction coil; activating the induction coil with an electrical current, to heat up the metallic or ferromagnetic solid to form a fluid; generating sonic energy to produce acoustic cavitation and abrasion between the abrasive particles and the container; and producing nanoparticles that comprise elements from the container, the metal and/or the ferromagnetic solid and the abrasive particles. Disclosed herein too is a composition comprising first metal or a first ceramic; and particles comprising carbides and/or nitrides dispersed therein. Disclosed herein too is a composition comprising nanoparticles comprising chromium carbide, iron carbide, nickel carbide, γ-Fe and magnesium nitride.

Manuel, Michele Viola↗

Method for producing nanoparticles and the nanoparticles produced therefrom

Disclosed herein is a method comprising disposing a container containing a metal and/or ferromagnetic solid and abrasive particles in a static magnetic field; where the container is surrounded by an induction coil; activating the induction coil with an electrical current, to heat up the metallic or ferromagnetic solid to form a fluid; generating sonic energy to produce acoustic cavitation and abrasion between the abrasive particles and the container; and producing nanoparticles that comprise elements from the container, the metal and/or the ferromagnetic solid and the abrasive particles. Disclosed herein too is a composition comprising first metal or a first ceramic; and particles comprising carbides and/or nitrides dispersed therein. Disclosed herein too is a composition comprising nanoparticles comprising chromium carbide, iron carbide, nickel carbide, γ-Fe and magnesium nitride.

Manuel, Michele Viola↗

Method for producing nanoparticles and the nanoparticles produced therefrom

Disclosed herein is a method comprising disposing a container containing a metal and/or ferromagnetic solid and abrasive particles in a static magnetic field; where the container is surrounded by an induction coil; activating the induction coil with an electrical current, to heat up the metallic or ferromagnetic solid to form a fluid; generating sonic energy to produce acoustic cavitation and abrasion between the abrasive particles and the container; and producing nanoparticles that comprise elements from the container, the metal and/or the ferromagnetic solid and the abrasive particles. Disclosed herein too is a composition comprising first metal or a first ceramic; and particles comprising carbides and/or nitrides dispersed therein. Disclosed herein too is a composition comprising nanoparticles comprising chromium carbide, iron carbide, nickel carbide, γ-Fe and magnesium nitride.

Manuel, Michele Viola↗

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

L-H transition trigger physics in ITER-similar plasmas with applied n = 3 magnetic perturbations

The L–H transition power threshold PLH is observed to increase with applied n = 3 resonant magnetic perturbations (RMP) in ITER-similar-shape plasmas with balanced neutral beam torque injection in DIII-D. The increase is most pronounced with added electron–cyclotron heating: PLH increases with decreasing edge plasma collisionality as PLH/PLH-08 ~ (ν*)−0.5, where PLH-08 is the 2008 ITPA multi-machine power threshold scaling. This result raises concerns for H-mode access at low edge collisionality in ITER, where RMP may have to be applied before the L–H transition to safely suppress the first edge-localized mode. Non-axisymmetric modifications with RMP include a simultaneous reduction of the radial electric field (Er) well depth and E × B shear. This can be attributed to increasing edge toroidal co-current rotation, and is consistent with substantially increased local long-wavelength turbulence (measured via beam emission spectroscopy). At high RMP perturbation strength the edge electric field Er reverses sign locally (becomes positive), with changes in dominant turbulence modes. Edge magnetic stochasticity provides an attractive explanation of the observed modifications, and the observed changes in toroidal rotation and Er are consistent with a simple fluid model describing radial electron current flow along stochastic fieldlines. The observed collisionality dependence of the L-mode edge electric field with applied RMP is also qualitatively consistent with this model. Reflectometry data indicate a significant reduction of the normalized L-mode radial density gradient a/Ln at high RMP field with simultaneous increase in radial particle flux and electron thermal flux from power balance analysis. We conjecture that the increase of PLH with RMP results from the combined effects of reduced E × B flow shear (increasing turbulent transport levels) and toroidal/poloidal flow modulation due to edge stochasticity. Initial experiments indicate that non-resonant n = 3 magnetic perturbations lead only to relatively small changes in Er, E × B shear and fluctuation characteristics, and have less impact on the L–H transition power threshold. This motivates further exploration of the RMP spectrum dependence of PLH for possible mitigation of the observed threshold increase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mechanically-, Thermally, and Chemically-Robust High-Temperature Ceramic Composites

In an effort to lower the levelized cost of electricity generated by Concentrated Solar Power (CSP) plants, the heat-to-electricity conversion efficiency may be increased by operating such plants with higher turbine inlet temperatures using closed Brayton cycles with high-pressure, supercritical carbon dioxide (sCO 2 ) as the working fluid. For example, if the turbine inlet temperature is raised to ≥750°C, in combination with sCO2 closed Brayton cycles, the relative heat-to-electricity efficiency may be raised by more than 20% (compared to subcritical steam Rankine cycles operating with turbine inlet temperatures ≤550°C). The associated reduction in the cost of dispatchable electricity from such CSP plants (utilizing high-temperature thermal energy storage) would be an important breakthrough on the path towards direct competition with fossil-fuel-based plants. A key barrier to achieve such high-temperature efficient operation has been the limited thermomechanical performance of metal alloys used in compact, primary printed circuit- type heat exchangers (PCHEXs) for heat transfer to high-pressure sCO 2 . The maximum allowed stresses for the use of conventional stainless steels and nickel-based superalloys at high sCO 2 pressures (≥20 MPa) decline rapidly at temperatures ≥550°C. This project has been focused on demonstrating the desired high-temperature properties and attractive manufacturing characteristics of mechanically-, thermally-, and chemically-robust ceramic/metal composites (cermets), in order to allow for the use of such composites in primary PCHEXs for heat transfer to sCO 2 at ≥750°C and ≥20 MPa in CSP plants. Prior work with high-temperature, co-continuous carbide/refractory metal cermets, such as ZrC/W composites, has already demonstrated: i) the attractive combination of high-temperature properties exhibited by such materials (e.g., high values of stiffness, failure strength, and thermal conductivity at 800°C), and ii) cost-effective processes for manufacturing PCHEX plates comprised of such cermets with tailorable channel patterns. However, ZrC/W cermets were not found to be inherently corrosion resistant in sCO 2 at 750°C and 20 MPa. While the application of a Cu coating to ZrC/W surfaces, along with a modest addition (50 ppm) of CO to CO 2 , rendered such composites resistant to corrosion for 1000 h in such sCO 2 -based fluids at 750°C, such additional steps would increase the design complexity of such cermet-based components. An alternative approach (the focus of this project) is to develop mechanically-robust, cost- effective cermets that are inherently resistant to high-temperature oxidation in air and CO 2 . The overall goals of this project have been to demonstrate: i) the oxidation resistance of at least one cermet at 750°C in air and in CO 2 (i.e., with a projected annual corrosion ≤30 μm), ii) the mechanically robust nature of at least one cermet at 750°C (i.e., with an average failure strength ≥200 MPa), and iii) that the selected cermet can be manufactured via low-cost forming techniques.

14 SOLAR ENERGY↗

Measured indoor PM2.5, black carbon, and oxidative potential before and after replacing gas with induction cooking in asthmatic households

Cooking is a major source of fine particulate matter (PM) in homes and evidence to date is inconclusive about the impact of cooking fuel on measures of residential PM exposure. The Cooking Energy and Ventilation Impacts on Children's Asthma (CEVICA) study measured cooking frequency, range hood use, indoor air quality (IAQ) and respiratory health indicators of children with asthma living in homes with gas stoves in California's San Joaquin Valley. Intensive measurements occurred over three 2-week intensive periods: at baseline and at the end of two consecutive 3-month study phases. Participants were randomly assigned to have their gas stoves replaced with electric induction at the start of Phase 1 or Phase 2. As part of the IAQ assessment, we collected particulate matter on Teflon filters using ultrasonic personal air samplers (UPAS). The filters were analyzed to quantify time-integrated PM 2.5 mass concentration, black carbon (BC) oxidative potential (OP). OP was measured using the SLF-OH assay, an acellular chemical assay that quantifies production of OH radicals in simulated lung fluid (SLF). Across all comparisons, PM 2.5 mass tended to be higher during gas cooking than induction electric cooking, although only one transition (Baseline to Phase 1) was statistically discernible at p<0.05. BC and mass-normalized OP showed no measurable differences across stove types or across phases. These preliminary findings suggest that stove fuel alone may not be the dominant determinant of indoor PM 2.5 concentrations in this sample, as previous studies have reported that indoor PM 2.5 during cooking can also be influenced by factors such as cooking practices, food and oil type, kitchen ventilation, and outdoor infiltration.

Lin, Zhuoying↗

Mobility of heavy rare earth elements in magmatic-hydrothermal systems: Experimental determination of DyPO 4 solubility in supercritical fluids

Dysprosium (Dy) is part of the heavy rare earth elements (REE) and a critical component in the production of high-performance super magnets used in wind turbines, electric vehicles, and other green technologies. The HREE phosphate, xenotime, is commonly associated with hydrothermal systems, where the mobility, fractionation, and deposition of REE are controlled by the stability of aqueous complexes and solubility of REE minerals which depend on temperature, pressure, and fluid composition (salinity, pH, and ligand availability). Here, we conducted a series of synthetic DyPO 4 solubility experiments in NaCl-HCl-bearing aqueous solutions (pH of 2 at 25 °C) from 500 to 700 °C and 1.5 kbar using externally heated pressure vessels. DyPO 4 displays a strong prograde solubility with temperature across all studied fluid salinities (0–0.5 m NaCl). The DyPO₄ solubility isotherms at 500 °C, 600 °C, and 700 °C show a complex relationship with salinity where DyPO₄ solubility is highest (39.2–1552 ppm at 500–700 °C) in 0 m NaCl solutions, decreases strongly (21–331 ppm at 500–700 °C) in the 0.06 m NaCl solutions and remains relatively constant in the 0.25 m (42.8–248 ppm) and 0.5 m NaCl (40.2–124 ppm) solutions. Thermodynamic modeling suggests that pH increases as a function of temperature and salinity, and controls solubility due to the stability of chloride complexes (i.e., DyCl +2 and DyCl 2 + ) at 500–600 °C, transitioning to hydroxyl complexes (i.e., Dy(OH) 3 0 ) dominating at 700 °C and relatively alkaline pH. Comparison with NdPO 4 solubility data indicates that the light REE in monazite is significantly more soluble than the heavy REE in xenotime in saline fluids at high temperature (500–700 °C), implying efficient light/heavy REE fractionation in supercritical fluids. The lower solubility of DyPO 4 relative to NdPO 4 suggests that at high temperatures, acidic (pH 2–4), low-salinity (∼5 wt% NaCl) fluids will favor the precipitation of xenotime over monazite as it will cool from 700 °C to below 500 °C. These observations align with paragenetic mineral assemblages in NYF pegmatites (e.g., Baveno, Italy), where xenotime commonly immobilizes heavy REE during the high-temperature hydrothermal stages.

58 GEOSCIENCES↗

Acceleration mechanisms of energetic ion debris in laser-driven tin plasma EUV sources

Laser-driven tin plasmas are driving new-generation nanolithography as sources of extreme ultraviolet (EUV) radiation centered at 13.5 nm. A major challenge facing industrial EUV source development is predicting energetic ion debris produced during the plasma expansion that may damage the sensitive EUV channeling multilayer optics. Gaining a detailed understanding of the plasma dynamics and ion acceleration mechanisms in these sources could provide critical insights for designing debris mitigation strategies in future high-power EUV sources. Here, we develop a fully kinetic model of tin-EUV sources using one-dimensional particle-in-cell simulations to study ion debris acceleration, which will be valuable for cross-validation of radiation-hydrodynamic simulations. An inverse-bremsstrahlung heating operator is used to model the interaction of a tin target with an Nd:YAG laser, and thermal conduction is included through a Monte Carlo Coulomb collision operator. While the large-scale evolution is in reasonable agreement with analogous hydrodynamic simulations, the significant timescale for collisional equilibration between electrons and ions allows for the development of prominent two-temperature features. A collimated flow of energetic ions is produced with a spectrum that is significantly enhanced at high energies compared to fluid simulations. The dominant acceleration mechanism is found to be a large-scale electric field supported mainly by the electron pressure gradient, which is enhanced in the kinetic simulations due to the increased electron temperature. We discuss the implications of these results for future modeling of tin-EUV sources and the development of debris mitigation schemes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A numerical solver for investigating the space charge effect on the electric field in liquid argon time projection chambers

This paper reports the development of a numerical solver aimed to simulate the interaction between the space charge (i.e. ions) distribution and the electric field in liquid argon time projection chamber (LArTPC) detectors. The ion transport equation is solved by a time-accurate, cell-centered finite volume method and the electric potential equation by a continuous finite element method. The electric potential equation updates the electric field which provides the drift velocity to the ion transport equation. The ion transport equation updates the space charge density distribution which appears as the source term in the electric potential equation. The interaction between the space charge distribution and the electric field is numerically simulated within each physical time step. The convective velocity in the ion transport equation can include the background flow velocity in addition to the electric drift velocity. The numerical solver has been parallelized using the Message Passing Interface (MPI) library. Numerical tests show and verify the capability and accuracy of the current numerical solver. It is planned that the developed numerical solver, together with a Computational Fluid Dynamics (CFD) package which provides the flow velocity field, can be used to investigate the space charge effect on the electric field in large-scale particle detectors.

42 ENGINEERING↗