SEARCH · Engineering Papers

Results for “Recombinant”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,423 records · Page 79

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Coordinate parameterisation and spectral method optimisation for Beltrami field solver in stellarator geometry

The numerical solution of the stepped pressure equilibrium (Hudson et al 2012 Phys. Plasmas 19 112502) requires a fast and robust solver to obtain the Beltrami field in three-dimensional geometry such as stellarators. The spectral method implemented in the stepped pressure equilibrium code (SPEC) is efficient when the domain is a hollow torus, but ill-conditioning of the discretised linear equations occurs in the (solid) toroid due to the artificially singular coordinate parameterisation near the axis. Here, we propose an improved choice for the reference axis to prevent coordinates surfaces from overlapping. Then, we examine the parity and asymptotics of the magnetic vector potential near the axis and suggest the use of recombined and rescaled Zernike radial basis functions. The maximum relative error in the magnetic field of the Wendelstein 7-X geometry is shown to reach 10 –9 at high resolution in a series of convergence tests and benchmarks against the boundary integral equation solver for Taylor states. The new method is also reported to significantly improve the accuracy of multi-volume SPEC calculations. A comparison between free-boundary SPEC and the analytical Dommaschk potential is presented with higher-than-usual Fourier resolution. It is illustrated that we are able to resolve low amplitude current sheets when an interface is placed where there is no flux surface in the analytic solution. This was previously concealed because of insufficient numerical resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Kinetic equilibrium reconstructions of plasmas in the MAST database and preparation for reconstruction of the first plasmas in MAST upgrade

Reconstructions of plasma equilibria using magnetic sensors were routine during operation of the Mega Ampere Spherical Tokamak (MAST) device, but reconstructions using kinetic profiles were not. These are necessary for stability and disruption analysis of the MAST database, as well as for operation in the upgrade to the device, MAST-U. The three-dimensional (3D) code VALEN is used to determine eddy currents in the 3D vessel structures for vacuum coil test shots, which are then mapped to effective resistances in the two-dimensional vessel groupings in the EFIT equilibrium reconstruction code to be used in conjunction with nearby loop voltage measurements for estimated currents in the structures during reconstruction. Kinetic equilibrium reconstructions with EFIT, using all available magnetic sensors as well as Thomson scattering measurements of electron temperature and density, charge exchange recombination spectroscopy measurements of ion temperature, and internal magnetic field pitch angle measurements from a motional Stark effect (MSE) diagnostic are performed for a large database of MAST discharges. Excellent convergence errors are obtained for the portions of the discharges where the stored energy was not too low, and it is found that reconstructions performed with temperature and density measurements but without MSE data usually already match the pitch angle measurements well. A database of 275 kinetic equilibria is used to test the ideal MHD stability calculation capability for MAST. In conclusion, the necessary changes to conducting structure in VALEN, and diagnostic setup in EFIT have been completed for the upgrade from MAST to MAST-U, enabling kinetic reconstructions to commence from the first plasma discharges of the upgraded device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Impurity transport studies at the HSX stellarator using active and passive CVI spectroscopy

The transport of carbon impurities has been studied in the helically symmetric stellarator experiment (HSX) using active and passive charge exchange recombination spectroscopy (CHERS). For the analysis of the CHERS signals, the STRAHL impurity transport code has been re-written in the python programming language and optimized for the application in stellarators. In addition, neutral hydrogen densities both along the NBI line of sight as well as for the background plasma have been calculated using the FIDASIM code. By using the basinhopping algorithm to minimize the difference between experimental and predicted active and passive signals, significant levels of impurity diffusion are observed. In this work, comparisons with neoclassical calculations from DKES/PENTA show that the inferred levels exceed the neoclassical transport by about a factor of four in the core and more than 100 times towards the plasma edge, thus indicating a high level of anomalous transport. This observation is in agreement with experimental heat diffusivites determined from a power balance analysis which exhibits strong anomalous transport as well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Experimental distinction of the molecularly induced Balmer emission contribution and its application for inferring molecular divertor density with 2D filtered camera measurements during detachment in JET L-mode plasmas

A previously presented model for generating 2D estimates of the divertor plasma conditions at JET from deuterium Balmer line intensity ratios, obtained from tomographic reconstructions of divertor camera images, was amended to consider also the Balmer emission arising from molecular processes. Utilizing the AMJUEL and H2VIBR atomic and molecular databases of EIRENE enabled also inference of the molecular divertor density from the distinguished molecularly induced emission. Analysis of a JET L-mode density scan suggests the molecularly induced emission accounting for up to 60%–70% and 10%–20% of the Balmer Dα and Dγ intensities, respectively, at the onset of detachment, while electron-ion recombination becomes increasingly dominant with deepening detachment. Similar observations were made by post-processing EDGE2D-EIRENE simulations, which indicated significant roles of molecular D$^{+}_{2}$ ions and vibrational excitation of the D 2 molecules as precursors for the molecularly induced emission. The experimentally inferred molecular density at the outer strike point was found to increase monotonously with decreasing strike point temperature, reaching approximately 30%–50% of the local electron density at n mol, osp = 1–2 X 10 20 m -3 at T e,osp ≈ 0.7 eV. A further steep increase by a factor of 3–5 was observed with decrease of T e,osp to 0.5 eV. The observations are in qualitative and reasonable quantitative agreement with EDGE2D-EIRENE predictions of n mol, osp within the uncertainties of the experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Characterizing W sources in the all-W wall, all-RF WEST tokamak environment

In this work, experimental data, together with interpretive modeling tools, are examined to study trends in the tungsten (W) source in the all-W environment of the WEST tokamak, both from the divertor and from the main chamber. In particular, a poloidal limiter protecting an ion cyclotron resonance heating (ICRH) antenna is used as proxy for main chamber sourcing. The key study is carried out by stepping up lower hybrid current drive (LHCD) power, as the only auxiliary power source. Limiter and divertor W sources exhibit a qualitatively similar proportionality to the total power crossing the separatrix, P SEP , although the main chamber source remains substantially lower than the divertor source, for the range of P SEP accessible in the experiments. Intepretive modeling of the limiter source is carried out with a particle-in-cell (PIC) sheath model coupled to a surface sputtering model. Oxygen is used as a proxy for all light impurity species allowing for characterization of the critical W erosion regions. To get a good quantitative match to the data, it is necessary to assume that the oxygen arrives at the surface mostly at high ionization stages (4+ and above). A separate simulation with SOLEDGE-EIRENE, constrained to measured upstream scrape-off-layer plasma profiles, gives oxygen fractional abundances that are compatible with the PIC simulation result. This is understood to arise from transport processes that dominate over recombination. Substituting the LHCD by ICRH, in an equivalent experiment, the local W source exhibits a 3× enhancement. This can be matched by the simulation, by assuming local RF electric field rectification, based on ~100 eV peak-to-peak, near-antennna electric field. This work has highlighted the particular importance of understanding the ion charge state balance of light impurities as these are most likely the dominant sputtering species in fusion devices with high-Z walls.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Investigation of core impurity transport in DIII-D diverted negative triangularity plasmas

Abstract Tokamak operation at negative triangularity has been shown to offer high energy confinement without the typical disadvantages of edge pedestals (Marinoni et al 2021 Nucl. Fusion 61 116010). In this paper, we examine impurity transport in DIII-D diverted negative triangularity experiments. Analysis of charge exchange recombination spectroscopy reveals flat or hollow carbon density profiles in the core, and impurity confinement times consistently shorter than energy confinement times. Bayesian inferences of impurity transport coefficients based on laser blow-off injections and forward modeling via the Aurora package (Sciortino et al 2021 Plasma Phys. Control. Fusion 63 112001) show core cross-field diffusion to be higher in L-mode than in H-mode. Impurity profile shapes remain flat or hollow in all cases. Inferred radial profiles of diffusion and convection are compared to neoclassical, quasilinear gyrofluid, and nonlinear gyrokinetic simulations. Heat transport is observed to be better captured by reduced turbulence models with respect to particle transport. State-of-the-art gyrokinetic modeling compares favorably with measurements across multiple transport channels. Overall, these results suggest that diverted negative triangularity discharges may offer a path to a highly-radiative L-mode scenario with high core performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Divertor plasma detachment: roles of plasma momentum, energy, and particle balances

Abstract Contrary to statements from some recent papers, the bifurcation of the scrape-off layer and divertor plasma parameters, as described within a 1D model, does exist. The analysis of both numerical simulations and experimental data shows that the real physical reasons that can cause a strong reduction of the plasma particle flux to the target j d are impurity radiation loss and volumetric plasma recombination. Fitting the parameter describing the ‘plasma momentum removal’, based on the results of 2D numerical simulations, then using this fit in subsequent analytic expressions for j d , and finally interpreting the evolution of j d as the result of elastic plasma-neutral interactions break the causality principle and lead to a wrong physical picture of divertor plasma detachment.

Physics↗

DOE OSTI2023

Linear instabilities in the hot-ion regime in a high-field spherical tokamak

Abstract The present operating high-field compact spherical tokamak ST40 is an important step toward an ST-based fusion reactor (Gryaznevich and Asunta 2017 Fusion Eng. Des. 123 177; McNamara 2023 Nucl. Fusion 63 054002). Temperature and density profiles and their uncertainties in recent hot ion ST40 plasmas with central ion temperatures exceeding 8.6 keV, have been inferred using an integrated analysis of several diagnostics including line-of-sight integration and volume average measurements, as well as limited profile information from a charge-exchange-recombination spectrometer. A linear gyrokinetic stability analysis has been carried out to identify the most unstable micro-instabilities in these hot ion plasmas. In one of these plasmas, it is found that linear growth rates of both ion- and electron-scale ( k θ ρ s ⩾ 0.2 where k θ is the poloidal wavenumber and ρ s is the ion gyro-radius with sound speed) modes decrease from the edge toward the core of the plasma (i.e. from ρ = 0.8 to 0.3 where ρ is the square root of the normalized magnetic flux), while the change of linear growth rates at k θ ρ s < 0.2 is non-monotonic. In particular, at ρ = 0.3 no unstable mode was found at k θ ρ s > 5 , and about an order of magnitude reduction in the maximum linear growth rate in the ion-scale wavenumber range of 0.2 ⩽ k θ ρ s < 1 is seen at ρ = 0.3 compared with the ρ = 0.8 location. It is found through parametric scans that while the ion temperature gradient (ITG) mode/trapped electron mode (TEM) and kinetic ballooning mode (KBM)/kinetic shear Alfvén (KSA) mode are more important in the plasma core, ubiquitous mode (UM) is found to be the important ion-scale instability at ρ ⩾ 0.6 . At ρ = 0.8, two instabilities are found at the electron scale, with one being the typical electrostatic electron temperature gradient (ETG) mode at ρ e scale and beyond and the other being an unidentified low-real-frequency instability (at the intermediate scale between electron and ion gyroradii) which can propagate in the ion direction and is driven by ETG. The existence of UM at the ion scale and the low-real-frequency instability and ETG mode at the electron scale at ρ = 0.8 is supported by an extensive linear gyrokinetic stability analysis of another ST40 hot ion plasma in a wider parametric range, taking into account the uncertainties in the experimental profiles.

ubiquitous mode↗

DOE OSTI2024

Benign termination of runaway electron beams on ASDEX Upgrade and TCV

Abstract This paper discusses the development of a benign termination scenario for runaway electron (RE) beams on ASDEX Upgrade and TCV. A systematic study revealed that a low electron density (n e ) companion plasma was required to achieve a large MHD instability, which expelled the confined REs over a large wetted area and allowed for the conversion of magnetic energy to radiation. Control of the companion plasma n e was achieved via neutral pressure regulation and was agnostic to material injection method. The neutral pressure required for recombination was found to be dependent on impurity species, quantity and RE current. On TCV,n e increased at neutral pressures above 1 Pa, indicating that higher collisionality between the REs and neutrals may lead to an upper pressure limit. The conversion of magnetic energy to radiated energy was measured on both machines and a decrease in efficiency was observed at high neutral pressure on TCV. The benign termination technique was able to prevent any significant increase in maximum heat flux on AUG from 200 to 600 kA of RE current, highlighting the ability of this approach to handle fully formed RE beams.

Physics↗

DOE OSTI2024

One-dimensional simulation and validation of divertor detachment induced through nitrogen seeding on HL-2A

Divertor detachment is a promising method to solve the power exhaust problem in tokamak devices or even in future magnetic fusion reactors. In this work, a detachment experiment (HL-2A shot #38008) with mixed gas seeding (60% nitrogen and 40% deuterium) is simulated using the SD1D module in BOUT++. In the process from attachment to detachment, the target electron temperature and the target ion saturation current in simulations are found to be consistent with the experimental results measured by Langmuir probes on the target plate. In order to understand the underlying detachment mechanism on HL-2A, this work analyses the role of different particle species in the cases with different seeding rates. It shows that the plasma density varies little and the density of neutrals (D and D 2 ) slightly decreases with the increase of nitrogen seeding rate, such that plasma–neutral interactions cannot effectively reduce plasma energy and plasma momentum in the divertor. The case with a high seeding rate predicts that increasing seeding rate cannot reach a target temperature lower than 2.5 eV, which is the required temperature for strong plasma–molecule interactions. Thus, the plasma–molecule interactions may not be important in the divertor during nitrogen seeding. Here this work also studies the parallel forces including (1) force due to the parallel electric field, (2) friction forces, (3) ion- and electron-thermal forces, and (4) collisional reactions (e.g. charge exchange recombination and ionisation). It is found that the friction force between nitrogen ions and other particle species (primarily D + ) is the dominant force towards the target, while ion- and electron-thermal forces are the dominant force pushing nitrogen ions back to upstream. Parallel forces determine the parallel distribution of nitrogen impurities, and therefore decide the region of nitrogen radiation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Electron ITB formation in EAST high poloidal beta plasmas under dominant electron heating

Abstract Plasma confinement and transport in tokamaks play a crucial role in the development of high poloidal beta steady-state operation scenarios. Therefore, it is very important to study the relevant mechanism of internal transport barriers (ITBs), which can help plasmas to obtain better confinement in order to achieve higher fusion gain. This paper mainly introduces the analysis of the characteristics of electron heat transport of discharges with ITB in high β P operation regime in Experimental Advanced Superconducting Tokamak (EAST). Based on the statistical analysis of stable discharges with β P > 1.5, it is found that there is an obvious bifurcation of the normalised electron temperature gradient (ETG) ( R / L T e ) in the range of β P = 2–2.2. Then the discharges of lower β P ( β P < 2, where the value of β P is below the bifurcation threshold) and of higher β P ( β P > 2.2, where the value of β P is above the threshold) were selected for analysis. The diagnostic data provided by Thomson scattering, x-ray crystal spectrometry and charge exchange recombination spectroscopy are used to provide reliable parameter profiles and then combined with the data of external magnetic probe measurements and the polarisation interferometer diagnosis system to fully reconstruct the balance. A relevant plasma current calculation model is used to calculate and analyse the current density profiles and power deposition, and then the transport analysis is carried out. Interestingly, in the higher β P discharges, it is found that the turbulence intensity provided by the CO 2 laser collective scattering system gradually decreases and the normalised ETG R / L T e gradually increases. It is also found that the electron heat transport coefficient decreases in the discharge with higher β P and the growth rate of the electron-scale turbulence calculated by transport gyro-Landau fluid (TGLF) is significantly reduced. Meanwhile, similar conclusions are also obtained in the discharges when the β P is further increased.

Wang, Z. H. (ORCID:0009000203431036)↗

DOE OSTI2024

Measurement of small island characteristics using high resolution ECE and CER at DIII-D

Abstract The measurements using the high resolution electron cyclotron emission radiometry and the charge exchange and recombination spectroscopy are processed using analytic formulas to allow for the detection of islands as small as 1.9 cm. In contrast to large, saturated magnetic islands which are relatively well understood to be governed by the loss of bootstrap current inside the island, small islands are less well understood due to the difficulty of their accurate measurement in tokamaks. Here, ‘small’ islands are islands comparable in size to the ion banana width, which can be as small as 0.8 cm at DIII-D. The new measurement methods allow for the detection of small island widths when the predicted increase of mode frequency to match the Doppler shifted ion diamagnetic frequency is observed. Therefore, for the first time, the mode frequency increase can be unambiguously associated to the acceleration of the magnetic island propagation. Such association allows for a further development and validation of the much-debated theory of ion polarization currents, which is thought to govern the small island growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Measurements of radial neutral density profiles from Balmer- α emission in Wendelstein 7-X

Radial neutral density profiles are estimated from measurements of passive H α emission in the Wendelstein 7-X stellarator. To parametrize the generally three-dimensional distribution with a low number of degrees of freedom, the neutral density is reduced to a flux surface quantity. Accounting for emission from excitation and recombination processes, neutral density profiles are derived independently for each of the available lines of sight. Density profiles obtained from the different viewing geometries are found to vary within one order of magnitude. Toroidally oriented lines of sight predict systematically lower neutral densities when compared to poloidally oriented ones. This discrepancy is attributed to the simplifications inherent in the imposed model and significant differences in integration volumes across the viewing geometries. In line with expectations, obtained neutral densities are found to decrease with increasing plasma density. Key restrictions of the model include the reduction of the neutral density to a flux surface quantity, uncertainties in the plasma profiles and instrument function, and line integration effects outside the last closed flux surface.

Wendelstein 7-X↗

DOE OSTI2025

Investigation of radiated-power for low aspect ratio fusion plasmas

Investigating radiated power in fusion plasmas is of utmost importance to understand the effect of undesired impurities, such as metals, in present devices, or also desired impurities, such as noble gasses, to purposefully radiate a large fraction of power in future devices. These studies are especially important for high Z impurities which will play a crucial role in future generation fusion pilot plants. In this work, we have developed a power radiation analysis module, which is used to investigate 2D distributions of impurity densities and radiated power asymmetries caused due to both plasma rotation and the cooling rate dependence on temperature profile for the cases of experimental NSTX plasmas and designed scenarios for the spherical tokamak advanced reactor (STAR), both being low aspect ratio tokamaks. Two different atomic databases have been tested during this work to study their impact on the radiated power distribution, especially for high Z impurities. Also in PRAM, self-consistent calculations of two-dimensional electron, main ion, and impurity ion densities are carried out using one-dimensional input density, temperature, and rotation profiles. In the case of NSTX, discharges with high rotation of ∼ 170 km s, measured with charge exchange recombination spectroscopy, have been investigated. Rotation-induced charge separation, leading to an electrostatic potential, is calculated iteratively to a self-consistent solution while testing high Z impurities to observe any 2D asymmetry in the core radiated power due to centrifugal forces. The STAR design, being much larger (R = 4 m), is projected to have a much lower rotation, and is shown to have low rotation-induced asymmetries, on the order of ten percent or less, between the low field and high field sides. However, another effect not due to rotation but to the dependence of impurity cooling rates on temperature can lead to radiation peaking off-axis, near the plasma edge. This effect is noticeable for argon in NSTX, for example, but can also be enhanced for certain impurities at much higher temperatures projected for STAR (T e0 ~ 32 keV), for example for undesired tungsten or possibly desired xenon.

NSTX↗

DOE OSTI2026

Continuum radiated power density $(P^{\textrm{C}}_{\textrm{rad}})$ and effective charge $(Z_{\textrm{eff}})$ estimates from multi-energy photon-counting measurements

Multi-energy soft x-ray pinhole cameras have been designed, built, calibrated, and operated at Madison symmetric Torus, Alcator C-Mod, and more recently at Tokamak a configuration variable and Tungsten Environment in Steady-State Tokamak (WEST), to measure plasma emission across multiple energy ranges. Here we describe a new methodology to estimate the local continuum radiated power density and the plasma effective charge (Z eff ) directly from photon-counting measurements of the line-free continuum emission (Bremsstrahlung and Radiative Recombination) in several energy bands between 11 and 18 keV. This capability is particularly valuable for confinement systems using metal plasma-facing components, where x-ray losses from interactions with the sputtered wall can represent a significant fraction of the total radiated power (P rad ). The approach leverages a well-characterized detector responsivity, modeled by a complementary error function, and interprets the differential multi-energy measurements between adjacent energy levels through the probability density function of a Gaussian distribution. The implementation of this diagnostic technique is currently under development on the WEST tokamak, aiming at the goal of providing real-time P rad and Z eff measurements during long-pulse operation (up to 1000 s) in the 2026 campaign.

WEST↗

DOE OSTI2026

Implementation of Pfirsch–Schlüter parallel flow in x-ray imaging crystal spectrometer inversion analysis

The x-ray imaging crystal spectrometer (XICS) tomographic inversion code for Wendelstein 7-X (W7-X) has been modified to consider the effects of parallel flows and has been applied to analyze measurements taken during recent experimental campaigns. Previous analysis neglected the effects of parallel flows due to the primarily perpendicular geometry of the sightlines and the small magnitude predicted by neoclassical theory. To reconsider these effects, the incompressibility condition for plasma flows is used to calculate the parallel Pfirsch–Schlüter flow component for the equilibrium configuration. By incorporating this condition along with the geometry of the sightlines—i.e. the fractional contributions of perpendicular and parallel flows—, an updated expression for the measured flow is used for the profile inversion. Application of this modified inversion code to data from W7-X shows that the magnitude of the radial electric field and the flux surface averaged perpendicular flow are reduced by approximately a factor of 2 and brought into better agreement with neoclassical predictions and charge exchange recombination spectroscopy measurements.

Pfirsch–Schlüter flows↗

DOE OSTI2022

Measurements of atoms and metastable species in N 2 and H 2 –N 2 nanosecond pulse plasmas

We report time-resolved, absolute number densities of metastable N 2 (A 3 Σ u + , v = 0, 1) molecules, ground state N 2 and H atoms, and rotational–translational temperature have been measured by tunable diode laser absorption spectroscopy and two-photon absorption laser-induced fluorescence in diffuse N 2 and N 2 –H 2 plasmas during and after a nanosecond pulse discharge burst. Comparison of the measurement results with the kinetic modeling predictions, specifically the significant reduction of the N 2 (A 3 Σ u + ) populations and the rate of N atom generation during the burst, suggests that these two trends are related. The slow N atom decay in the afterglow, on a time scale longer than the discharge burst, demonstrates that the latter trend is not affected by N atom recombination, diffusion to the walls, or convection with the flow. This leads to the conclusion that the energy pooling in collisions of N 2 (A 3 Σ u + ) molecules is a major channel of N 2 dissociation in electric discharges where a significant fraction of the input energy goes to electronic excitation of N 2 . Additional measurements in a 1% H 2 –N 2 mixture demonstrate a further significant reduction of N 2 (A 3 Σ u + , v = 0, 1) populations, due to the rapid quenching by H atoms accumulating in the plasma. Comparison with the modeling predictions suggests that the N 2 (A 3 Σ u + ) molecules may be initially formed in the highly vibrationally excited states. The reduction of the N 2 (A 3 Σ u + ) number density also diminishes the contribution of the energy pooling process into N2 dissociation, thus reducing the N atom number density. The rate of N atom generation during the burst also decreases, due to its strong coupling to N 2 (A 3 Σ u + , v) populations. On the other hand, the rate of H atom generation, produced predominantly by the dissociative quenching of the excited electronic states of N 2 by H 2 , remains about the same during the burst, resulting in a nearly linear rise in the H atom number density. Comparison of the kinetic model predictions with the experimental results suggests that the yield of H atoms during the quenching of the excited electronic state of N 2 by molecular H 2 is significantly less than 100%. The present results quantify the yield of N and H atoms in high-pressure H 2 –N 2 plasmas, which have significant potential for ammonia generation using plasma-assisted catalysis

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Nonequilibrium plasma in water vapor

A kinetic model of a non-equilibrium water vapor plasma created by a nanosecond pulsed discharge at elevated temperatures and at low pressure is presented. The model includes equations for electron, vibrational and gas temperatures combined with equations for neutral, charged, electronically and vibrationally excited species. Simulation results are compared with the available experimental data and show good agreement in dynamics of the plasma decay and electron drift velocities. It is shown that the plasma decay in pure water vapor is controlled by the dissociative recombination. Under the conditions simulated, the concentration of electronically excited species is low, but we observe a formation of vibrationally excited species. Finally, a significant decrease in the electrical conductivity of water plasma with time, which depends on the level of initial ionization, is shown.

42 ENGINEERING↗