Numerical simulations of axisymmetric nanosecond pulsed discharges at atmospheric pressure and temperature are performed with a novel fully implicit time integration approach. The plasma fluid equations with a drift-diffusion model and local-field approximation are made dimensionless and solved using a preconditioned Jacobian Free Newton-Krylov method. A simplified kinetics model is employed, including electrons, one positive ion, and one negative ion. The chemical processes of ionization, attachment, detachment, and recombination are considered along with photoionization. The newly developed fully-implicit integration scheme with physics-based preconditioning allows for the efficient simulations capable of describing the cathode sheath over time-scales of O(10 us). The implicit solver overcomes the limiting time scales related to electron drift, diffusion, dielectric relaxation, and ionization.
more »
« less
Sheath formation around a dielectric droplet in a He atmospheric pressure plasma
Interactions at the interface between atmospheric pressure plasmas and liquids are being investigated to address applications ranging from nanoparticle synthesis to decontamination and fertilizer production. Many of these applications involve activation of droplets wherein the droplet is fully immersed in the plasma and synergistically interacts with the plasma. To better understand these interactions, two-dimensional modeling of radio frequency (RF) glow discharges at atmospheric pressure operated in He with an embedded lossy dielectric droplet (tens of microns in size) was performed. The properties of the sheath that forms around the droplet were investigated over the RF cycle. The electric field in the bulk plasma polarizes the dielectric droplet while the electron drift in the external electric field is shadowed by the droplet. The interaction between the bulk and sheath electric fields produces a maximum in E/N (electric field/gas number density) at the equator on one side of the droplet where the bulk and sheath fields are aligned in the same direction and a minimum along the opposite equator. Due to resistive heating, the electron temperature T e is maximum 45° above and below the equator of the droplet where power deposition per electron is the highest. Although the droplet is, on the average, negatively charged, the charge density on the droplet is positive on the poles and negative on the equator, as the electron motion is primarily due to diffusion at the poles but due to drift at the equator.
more »
« less
- PAR ID:
- 10382510
- Date Published:
- Journal Name:
- Journal of Applied Physics
- Volume:
- 132
- Issue:
- 8
- ISSN:
- 0021-8979
- Page Range / eLocation ID:
- 083303
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Since the 1950s, high frequency and very high frequency radars near the magnetic equator have frequently detected strong echoes caused ultimately by the Farley‐Buneman instability (FBI) and the gradient drift instability (GDI). In the 1980s, coordinated rocket and radar campaigns made the astonishing observation of flat‐topped electric fields coincident with both meter‐scale irregularities and the passage of kilometer‐scale waves. The GDI in the daytimeEregion produces kilometer‐scale primary waves with polarization electric fields large enough to drive meter‐scale secondary FBI waves. The meter‐scale waves propagate nearly vertically along the large‐scale troughs and crests and act as VHF tracers for the large‐scale dynamics. This work presents a set of hybrid numerical simulations of secondary FBIs, driven by a primary kilometer‐scale GDI‐like wave. Meter‐scale density irregularities develop in the crest and trough of the kilometer‐scale wave, where the total electric field exceeds the FBI threshold, and propagate at an angle near the direction of total Hall drift determined by the combined electric fields. The meter‐scale irregularities transport plasma across the magnetic field, producing flat‐topped electric fields similar to those observed in rocket data and reducing the large‐scale wave electric field to just above the FBI threshold value. The self‐consistent reduction in driving electric field helps explain why echoes from the FBI propagate near the plasma acoustic speed.more » « less
-
Abstract Non-equilibrium plasmas derive their low temperature reactivity from producing and driving energetic electrons and active species under large electric fields. Therefore, the impact of reactants on the plasma properties including electron number density, electric field, and electron temperature is critical for applications such as plasma methane (CH 4 ) reforming. Due to experimental complexity, electron properties and the electric field are rarely measured together in the same discharge. In this work, we combine time-resolved Thomson scattering and electric field induced second harmonic generation to probe electron temperature, electron density, and electric field strength in a 60 Torr CH 4 /Ar nanosecond-pulsed dielectric barrier discharge while varying the CH 4 mole fraction from 0% to 8%. These measurements are compared to a 1D numerical model to benchmark its predictions and identify areas of uncertainty. Nonlinear coupling between CH 4 addition, electron temperature, electron density, and the electric field was directly observed. Contrary to previous measurements in He, the electron temperature increased with CH 4 mole fraction. This rise in electron temperature is identified as electron heating by residual electric fields that increased with larger CH 4 mole fraction. Moreover, the electron number density has been found to decrease rapidly with the increase of methane mole fraction. Comparison of these measurements with the model yielded better agreement at higher CH 4 mole fractions and with the usage of ab initio calculated Ar electron-impact cross-sections from the B-spline R-matrix database. Furthermore, the calculated plasma properties are shown to be sensitive to the residual surface charge implanted on the quartz dielectric surfaces. Without considering surface charge in the simulations, the calculated electric field profiles agreed well with the measurements, but the electron properties were underpredicted by more than a factor of three. Therefore, measurements of either the electric field or electron properties measurements alone are insufficient to fully validate modeling predictions.more » « less
-
Abstract Atmospheric pressure plasma jets (APPJs) are increasingly being used to functionalize polymers and dielectric materials for biomedical and biotechnology applications. Once such application is microfluidic labs-on-a-chip consisting of dielectric slabs with microchannel grooves hundreds of microns in width and depth. The periodic channels, an example of a complex surface, present challenges in terms of directly and uniformly exposing the surface to the plasma. In this paper, we discuss results from computational and experimental investigations of negative APPJs sustained in Ar/N2mixtures flowing into ambient air and incident onto a series of microchannels. Results from two-dimensional plasma hydrodynamics modeling are compared to experimental measurements of electric field and fast-camera imaging. The propagation of the plasma across dry microchannels largely consists of a sequence of surface ionization waves (SIWs) on the top ridges of the channels and bulk ionization waves (IWs) crossing over the channels. The IWs are directed into electric field enhanced vertices of the next ridge. The charging of these ridges produce reverse IWs responsible for the majority of the ionization. The propagation of the plasma across water filled microchannels evolve into hopping SIWs between the leading edges of the water channels, regions of electric enhancement due to polarization of the water. Positive, reverse IWs follow the pre-ionized path of the initial negative waves.more » « less
-
Abstract We present a number of unique observations of ionospheric anomalies following the Hunga‐Tonga Hunga‐Ha'apai (HTHH) volcanic eruption on 15 January 2022. All are based on non‐dedicated geodetic satellite systems: Global Positioning System tracking of Low Earth Orbit (LEO) CubeSats, intersatellite tracking between two GRACE Follow‐On satellites, satellite radar altimeters to the ocean surface, and Doppler radio beacons from ground stations to LEO geodetic satellites. Their observations revealed the development of anomalously large trough‐like plasma depletions, along with plasma bubbles, in the equatorial regions of the Pacific and East Asian sectors. Trough‐like plasma depletions appeared to be confined within approximately ±20° magnetic latitude, accompanied by density enhancements just outside this latitude range. These plasma depletions and enhancements were aligned with the magnetic equator and occurred across broad longitudes. They were detected in regions where atmospheric waves from the HTHH eruption passed through around the time of the sunset terminator. We interpret these phenomena in terms of theEdynamo electric fields driven by atmospheric waves from the eruption. The uplift of the ionosphere beyond satellite altitudes, followed by subsequent plasma diffusion to higher latitudes along magnetic field lines, results in the formation of trough‐like plasma depletions around the magnetic equator and density enhancement at higher latitudes. The detection of plasma bubbles in the Asian sector during the non‐bubble season (January) is likely associated with the uplift of the ionosphere at the sunset terminator.more » « less
An official website of the United States government

