Search NASA⌕ Search

SEARCH · Search NASA

Results for “HEAT EQUATION”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Coupled model for liquid lithium plasma facing components

Numerical analysis provides the design choice and operating window of liquid metal Plasma Facing Components (PFC) concepts. Coupled analysis of boundary plasma together with the surrounding boundary structures is required. Here, to achieve this goal, PPPL is developing a comprehensive multi-physics model for modeling of PFCs in fusion devices. The model includes the fluid-kinetic code SOLPS-ITER and the flow and heat transfer code CFX from ANSYS. SOLPS-ITER was augmented with a liquid metal boundary condition algorithm, allowing direct two-way coupling of the plasma analysis with the two-dimensional analytical slab flow model which includes heat convection in the liquid metal PFC. The target heat flux resulting from this coupled analysis is used as a boundary condition for detailed 3D Computational Fluid Dynamics (CFD) Magneto Hydro Dynamics (MHD) and heat transfer analysis. A new formulation of MHD equations is introduced in the numerical procedure ensuring current conservation of the discretized equations. Results of the 3D analysis are used for final validation of the coupled model. A PFC design where a porous wall is used to stabilize the liquid metal surface, while MHD drive is used to push the liquid metal flow inside the PFC, will be investigated in the regimes where vapor shielding is created for enhanced volumetric plasma heat dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Simple Heat Transfer Model for Film Cooling Applications

This report describes the development of a simple engineering model for film cooling. This model is used to derive a relationship between local wall temperature variations and key cooling performance parameters like local heat transfer coefficients and film effectiveness. This relation and method new and different from previously published models. The scope of this report includes the derivation of regression model equations for a flat plate with and without film cooling. The model equation for a flat plate without film cooling can be used to estimate local heat transfer coefficients using surface temperatures measured from infrared thermography. The model equation for the flat plate with film cooling can be used to estimate film cooling effectiveness, $η_f$, and heat transfer augmentation from the film cooling jet(s).

42 ENGINEERING↗
DOE OSTI2024

High-Fidelity Arc-Discharge Model for Hydrogen-Plasma-Smelting-Reduction of Iron Ore

Electrification and use of renewable hydrogen is currently a necessity for decarbonizing the iron-and-steel industry. In this regard, hydrogen plasma smelting reduction (HPSR) is a novel pathway that is being explored for reduction of iron ore. HPSR provides several decarbonization merits compared to conventional blast furnaces. Firstly, the use of renewable hydrogen drastically reduces the CO2 emissions compared to the use of coke. Secondly, renewable electricity in the form of a thermal plasma for making reactive hydrogen species (radicals, ions) are more efficient at reducing iron ore compared to neutral H2. Thirdly, a molten product compatible with downstream processes is obtained from the intense heat transfer from the plasma. However, the scale-up of this technology requires fundamental exploration of hydrogen plasma dynamics and its interaction with complex solid material that include phase changing iron-ore and slag. In this work, we present a first principles continuum scale model for thermal plasmas in Ar/H2 gas mixtures typically used for HPSR. The thermal plasma governing equations for mass, momentum and energy with Lorentz force and Joule heating source terms are solved along with electromagnetic equations for electrostatic and magnetic vector potential. Our solver will be based on Pele, a suite of reacting flow solvers designed for advanced scientific computing architectures (Henry De Frahan et al., Proceedings of SIAM Parallel Processing, 13-25, 2024), and will utilize adaptive mesh generation for enhanced resolutions at locations of intense physicochemical interactions. This study will present the impact of Ar to H2 ratios on excited/dissociated hydrogen species concentrations, plasma temperature and conductivity along with the impact of outgassed species (water, metal vapor, O, OH radicals) from ore surface on gas phase chemistry. Furthermore, the heat and species flux to the surface will be quantified as a function of applied voltages in a transferred arc configuration.

hydrogen plasma↗
DOE OSTI2023

A First Principles Approach to Spectral Phonon Transport in Heterostructures

Understanding thermal transport across interfaces which give rise to a thermal resistance (also known as Kapitza resistance) is a critical issue affecting the development of nanotechnologies. Much modern and emergent nanotechnology consist of adjacent materials, and phonon mediated heat transfer governs thermal behavior across internal interfaces in these devices. The physics of thermal transport in solids are governed both by phenomena occurring at the atomic scale and interactions with the material's microstructure. The forecasting of fundamental quantities such as temperature, heat flux and thermal conductivity typically employs the semi-classical Boltzmann transport equation to predict the macroscopic behavior of materials in terms of the microscopic dynamics of its heat carriers. Kapitza resistance was first discovered in liquid helium experiments and has led to a fundamental research thrust in micro and nano-scale heat transport, the behavior of thermal carriers across internal interfaces. Thermal interfacial resistance (TIR) is a widely studied phenomenon, first engaged by Swartz and Pohl through their development of the acoustic and diffuse mismatch methods, then continued through myriad efforts with varying methods and approaches in an attempt to resolve carrier behavior at thermal interfaces. Many of the fundamental approaches to TIR have been at the nanoscale, and research is conducted with molecular dynamics (MD) and density functional theory (DFT) methods. The limitations of these methods is system size; atomistic methods tend to be limited to system sizes of 100,000 atoms or less. Larger length-scale methods have also been pursued, based on the principles of acoustic or diffuse mismatch, but not all include simulation of TIR using a full phonon band spectrum, or temperature dependent methods. Our approach to enabling phonon transport in layered materials draws upon our previous work of demonstrating spectrally coupled phonon transport in homogeneous and heterogeneous materials. We use a semi-analytical approach in which the Bose-Einstein (B-E) statistics set the strength of the phonon radiance in a frequency group, but the B-E statistics are informed with information from the transport system. The B-E statistics in a single frequency group feels the influence of all the groups through the spatial temperature. We also include a new field term which is an indicator of the amount of non-equilibrium behavior of the phonon spectrum---this is added to the phonon source term in all groups to ensure closure and conservation of energy, as the phonon groups in the transport system and the analytical systems are coupled. This work builds upon our previous approach by adding a phonon coupling term at an internal interface, using the principles of the DMM through transmission and reflection coefficients. In this work, the coefficients are determined through computing a common temperature at the interface, influenced by the phonon band structure of both materials, in effect, providing mixing between the two material systems and using the common temperature to set the strength of the phonon radiance at the boundaries on either side of the interface. Our approach uses material properties computed along various crystallographic orientations, and while some isotropy is built into the interface condition, the material properties weight the phonon distributions in the proper crystalline direction. Greater resolution of phonon behavior in proximity to an interface, and more accurate predictions of TIR are obtained. While it is true the assumption of diffuse mismatch can yield inconsistent results compared to experiment especially at low temperatures, this work focuses on room temperature and beyond effects, for future applications in nuclear fuel, or thermoelectric devices; a modified mismatch approach may be feasible if applied properly. Additionally, our methods focus on bridging mesoscale to engineering scale

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fueling limits in a cylindrical viscosity-limited reactor

Recently, a method to achieve a “natural hot-ion mode” was suggested by utilizing ion viscous heating in a rotating plasma with a fixed boundary. We explore the steady-state solution to the Braginskii equations and find the parameter regime in which a significant temperature difference between ions and electrons can be sustained in a driven steady state. Here, the threshold for this effect occurs at ρ i ≳ 0.1R. An analytic, leading order low flow solution is obtained, and a numerical, moderate Mach number M ≲ 2 is investigated. The limitation is found to be at moderate Mach numbers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A Package of Momentum and Heat Transfer Coefficients for the Stable Surface Layer Extended by New Coefficients over Sea Ice

The bulk transfer coefficients of momentum, heat, and humidity belong to the main ingredients of numerical weather prediction and climate models. They are needed for the calculation of turbulent fluxes in the surface layer and often rely on the Monin–Obukhov similarity theory requiring universal stability functions. The problem of a derivation of transfer coefficients based on different stability functions has been considered by many researchers over the years but it remains to this day. In this work, dedicated to the memory of S.S. Zilitinkevich, we also address this task, and obtain transfer coefficients from three pairs of theoretically derived stability functions suggested by Zilitinkevich and co-authors for stable conditions. Additionally, we construct non-iterative parametrizations of these transfer coefficients based on earlier work. Results are compared with state-of-the-art coefficients for land, ocean, and sea ice. The combined parametrizations form a package in a universal framework relying on a semi-analytical solution of the Monin-Obukhov similarity theory equations. A comparison with data of the Surface Heat Budget of the Arctic Ocean campaign (SHEBA) over sea ice reveals large differences between the coefficients for land conditions and the measurements over sea ice. However, two schemes of Zilitinkevich and co-authors show, after slight modification, good agreement with SHEBA although they had not been especially developed for sea ice. One pair of the modified transfer coefficients is superior and is compatible to earlier SHEBA-based parametrizations. Finally, an algorithm for practical use of all transfer coefficients in climate models is given.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Theory for transport in magnetized plasmas

We present a two fluid model for anomalous thermal, particle and momentum transport which has been derived from kinetic theory using a small parameter ε~10-2. This small parameter can come from either the ratio of wave frequency to cyclotron frequency or the ratio of potential to thermal energies. This fluid model is strongly nonlinear and has an exact closure within the accuracy ε. It can be used to simulate the full radius of a tokamak including internal and edge transport barriers at the same time without introducing any local assumptions for physics model or grid size. It involves a full transport matrix, thus giving possibility for both density and heat pinches. The model can be derived from the Vlasov equation, adding also close collisions and includes particle pinches of the type found in the magnetosphere and in the Massachusetts Institute of Technology elevated dipole machine. Pinches are also in agreement for the heat pinch experiments on Doublet III D including simultaneous particle and heat pinches. For tokamaks it gives the right scaling of confinement time versus heating power. Finally, we have also recovered the Turbulent Equipartition result that the density scales as inverse safety factor q in typical simulations. Thus this model is both fundamental and of high practical use.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Design of a Combined Proportional Integral Derivative Controller to Regulate the Temperature Inside a High-Temperature Tubular Solar Reactor

Solar fuels are proven to be promising candidates for thermochemical energy storage. However, the transient nature of solar radiation is an obstacle to maintaining a stable operational temperature inside a solar reactor. To overcome this challenge, the temperature of a solar reactor can be regulated by controlling the incoming solar radiation or the feedstock flowrate inside the reactor. In this work, a combined proportional integral derivative (PID) controller is implemented to regulate the temperature inside a high-temperature tubular solar reactor with counter-current flowing gas/particles. The control model incorporates two control systems to regulate incoming solar radiation and gas flow simultaneously. The design of the controller is based on a reduced-order numerical model of a high-temperature tubular solar reactor that is vertically oriented with an upward gas flow and downward particle flow. The reactor receives heat circumferentially through its wall over a finite segment of its length. Formulation of the heat transfer model is presented by applying the energy balance for the reactor tube and considering heat and mass transfer inside. A set of governing differential equations are solved numerically by using the finite volume method to obtain reactor wall, particles, and gas temperatures along the reactor length with various boundary conditions. Simulation results are used to tune the PID controller parameters by utilizing the Ziegler–Nichols tuning method. Both the simulation results and the controller performance are visualized on the LABVIEW platform. The controller is challenged to track different temperature setpoints with different scenarios of transient solar radiation. Here, the performance of the PID controller was compared to experimental results obtained from an industrial PID controller embedded in a 7 kW electric furnace. Results show that the combined PID controller is successful in maintaining a stable temperature inside the reactor by regulating the incoming solar radiation and the flowrate via small steady-state error and reasonable settling time and overshoot.

14 SOLAR ENERGY↗
DOE OSTI2024

Site-Projected Thermal Conductivity: Application to Defects, Interfaces, and Homogeneously Disordered Materials

With the rapid advance of high-performance computing and electronic technologies, understanding thermal conductivity in materials has become increasingly important. This study presents a novel method: the site-projected thermal conductivity that quantitatively estimates the local (atomic) contribution to heat transport, leveraging the Green–Kubo thermal transport equations. The effectiveness of this approach on disordered and amorphous graphene, amorphous silicon, and grain boundaries in silicon–germanium alloys is demonstrated. Amorphous graphene reveals a percolation behavior for thermal transport. The results highlight the potential of the method to provide new insights into the thermal behavior of materials, offering a promising avenue for materials design and performance optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Conservative discontinuous Galerkin schemes for nonlinear Dougherty–Fokker–Planck collision operators

Here, we present a novel discontinuous Galerkin algorithm for the solution of a class of Fokker–Planck collision operators. These operators arise in many fields of physics, and our particular application is for kinetic plasma simulations. In particular, we focus on an operator often known as the ‘Lenard–Bernstein’ or ‘Dougherty’ operator. Several novel algorithmic innovations, based on the concept of weak equality, are reported. These weak equalities are used to define weak operators that compute primitive moments, and are also used to determine a reconstruction procedure that allows an efficient and accurate discretization of the diffusion term. We show that when two integrations by parts are used to construct the discrete weak form, and finite velocity-space extents are accounted for, a scheme that conserves density, momentum and energy exactly is obtained. One novel feature is that the requirements of momentum and energy conservation lead to unique formulas to compute primitive moments. Careful definition of discretized moments also ensure that energy is conserved in the piecewise linear case, even though the kinetic-energy term, $v^2$ is not included in the basis set used in the discretization. A series of benchmark problems is presented and shows that the scheme conserves momentum and energy to machine precision. Empirical evidence also indicates that entropy is a non-decreasing function. The collision terms are combined with the Vlasov equation to study collisional Landau damping and plasma heating via magnetic pumping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Determining the electrical conductivity of metals using the 2 MA Thor pulsed power driver

We present the development of a pulsed power experimental technique to infer the electrical conductivity of metals from ambient to high energy density conditions. The method is implemented on Thor, a moderate scale (1–2 MA) pulsed power driver. The electrical conductivity of copper at elevated temperature (>4000 K) and pressure (>10 GPa) is determined, and a new tabular material model is developed, guided by density functional theory, which preserves agreement with existing experimental data. Minor modifications (<10%) are found to be necessary to the previous Lee–More–Desjarlais model isotherms in the vicinity of the melt transition in order to account for observed discrepancies with the new experimental data. An analytical model for magnetic direct drive flyer acceleration and Joule heating induced vaporization based on the Tsiolkovsky “rocket equation” is presented to assess sensitivity of the method to minor changes in electrical conductivity.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Transport physics dependence of Bohm speed in presheath–sheath transition

We report the ion exit flow speed at the sheath entrance is constrained by the Bohm criterion, which is used as a boundary condition for simulations that do not resolve the sheath region. Traditional Bohm criterion analysis invokes the equation of state and, thus, ignores transport physics in the sheath transition problem. An expression for the Bohm speed away from the asymptotic limit is derived from a set of anisotropic plasma transport equations. The thermal force, collisional temperature isotropization, and heat flux enter into the evaluation of the Bohm speed. By comparison with kinetic simulation results, this expression is shown to be accurate in the presheath–sheath transition region rather than a single point at the sheath entrance over a broad range of collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A Tailored Convolutional Neural Network for Nonlinear Manifold Learning of Computational Physics Data Using Unstructured Spatial Discretizations

In this work, we propose a nonlinear manifold learning technique based on deep convolutional autoencoders that is appropriate for model order reduction of physical systems in complex geometries. Convolutional neural networks have proven to be highly advantageous for compressing data arising from systems demonstrating a slow-decaying Kolmogorov n-width. However, these networks are restricted to data on structured meshes. Unstructured meshes are often required for performing analyses of real systems with complex geometry. Our custom graph convolution operators based on the available differential operators for a given spatial discretization effectively extend the application space of deep convolutional autoencoders to systems with arbitrarily complex geometry that are typically discretized using unstructured meshes. We propose sets of convolution operators based on the spatial derivative operators for the underlying spatial discretization, making the method particularly well suited to data arising from the solution of partial differential equations. We demonstrate the method using examples from heat transfer and fluid mechanics and show better than an order of magnitude improvement in accuracy over linear methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Progress report on analytic and numerical studies of x-ray-induced impulse

In a variety of high-energy-density (HED) systems, x-rays of a given energy are used to generate shockwaves, bulk motion, and impulse in materials. As radiation energy is deposited within the material, the material is heated. The number of photons of a certain wavelength absorbed is determined by the spectral intensity of the radiation and the material’s opacity evaluated at that wavelength. This heating results in a pressure increase dictated by the material’s equation of state. Depending on the intensity of the absorbed radiation, it may also cause the material to change phase into a liquid, gas, or plasma. The increased pressure drives the heated surface layer to blow off, imparting impulse and sending a compressive wave into the bulk of the material. Additionally, the compression wave interacts with the solid boundary of the material, resulting in a tensile wave that may cause the material to spall. The impulse generated by the deposition of x-ray energy within the sample can be modeled using purely analytical methods, e.g. the Bethe, Bade, Averell, and Yost (BBAY) model. However, the blow-off process is rather complicated, and proper modeling efforts must account for material ejected by spallation, vaporization, jetting, and plasma ablation. For this reason, analytical models have an unclosed term describing the final energy of the blown-off material Ef(z). Prior modeling efforts have arbitrarily fixed this at some value or modeled it using a limiting set of thermodynamic assumptions. The work we are currently performing uses validated simulations using sophisticated photon transport, equation of state, and strength models/data to provide a fit for Ef(z) that is useful for predictive calculation of impulse. We will apply our methodology and show results for different materials and x-ray sources. This work is particularly useful for the design of experiments studying x-ray impulse generation. The present report is outlined as follows. Section 1.2 describes a series of HED experiments investigating x-ray-generated impulse in materials. Section 1.3 describes the computational approach we employ in this study, and presents validation results against the aforementioned experiments. Section 1.4 introduces analytical models for impulse generation, as well as our method for utilizing impulse from simulations to close the models. We discuss the concept of impulse-spectrum sensitivity, it’s application to uncertainty quantification, and derive a very useful analytical expression for it in 1.5. We then summarize recent progress in this project and discuss future work in section 1.6.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The discrete Green's function paradigm for two-way coupled Euler–Lagrange simulation

We outline a methodology for the simulation of two-way coupled particle-laden flows. The drag force that couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. In addition to hydrodynamics, the method can be extended to other physics present in particle-laden flows such as heat transfer and electromagnetism. The discrete Green's functions for the Navier–Stokes equations are obtained at low particle Reynolds number in a two-plane channel geometry. We perform verification at different Reynolds numbers for a particle settling under gravity parallel to a plane wall, for different wall-normal separations. Compared with other point-particle schemes, the Stokesian discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations. To account for degradation in accuracy away from the wall at finite Reynolds number, we extend the present methodology to an Oseen-like discrete Green's function. The extended discrete Green's function method is found to be accurate within 6% at all wall-normal separations for particle Reynolds numbers up to 24. Furthermore, the discrete Green's function approach is well suited to dilute systems with significant mass loading and this is highlighted by comparison against other Euler–Lagrange as well as particle-resolved simulations of gas–solid turbulent channel flow. Strong particle–turbulence coupling is observed in the form of turbulence modification and turbophoresis suppression, and these observations are placed in context of the different methods.

42 ENGINEERING↗
DOE OSTI2021

Transient heat transfer of superfluid He 4 in nonhomogeneous geometries: Second sound, rarefaction, and thermal layer

Transient heat transfer in superfluid 4 He (He II) is a complex process that involves the interplay of the unique counterflow heat-transfer mode, the emission of second-sound waves, and the creation of quantized vortices. Many past researches focused on homogeneous heat transfer of He II in a uniform channel driven by a planar heater. Here, we report our systematic study of He II transient heat transfer in nonhomogeneous geometries that are pertinent to emergent applications. By solving the He II two-fluid equations of motion coupled with Vinen's equation for vortex-line density, we examine and compare the characteristics of transient heat transfer from planar, cylindrical, and spherical heaters in He II. Our results show that as the heater turns on, an outgoing second-sound pulse emerges, within which the vortex-line density grows rapidly. These vortices attenuate the second sound and result in a heated He II layer in front of the heater, i.e., the thermal layer. In the planar case where the vortices are created throughout the space, the second-sound pulse is continuously attenuated, leading to a thick thermal layer that diffusely spreads following the heat pulse. On the contrary, in the cylindrical and the spherical heater cases, vortices are created mainly in a thin thermal layer near the heater surface. As the heat pulse ends, a rarefaction tail develops following the second-sound pulse, in which the temperature drops. This rarefaction tail can promptly suppress the thermal layer and take away the deposited thermal energy. The effects of the heater size, heat flux, pulse duration, and temperature on the thermal-layer dynamics are discussed. We also show how the peak heat flux for the onset of boiling in He II can be studied in our model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Unification of nonequilibrium molecular dynamics and the mode-resolved phonon Boltzmann equation for thermal transport simulations

Nano-size confinement induces many intriguing non-Fourier heat conduction phenomena, such as nonlinear temperature gradients, temperature jumps near the contacts, and size-dependent thermal conductivity. Over the past decades, these effects have been studied and interpreted by nonequilibrium molecular dynamics (NEMD) and phonon Boltzmann transport equation (BTE) simulations separately, but no theory that unifies these two methods has ever been established. In this work, we unify these methods using a quantitative mode-level comparison and demonstrate that they are equivalent for various thermostats. We show that different thermostats result in different non-Fourier thermal transport characteristics due to the different mode-level phonon excitations inside the thermostats, which explains the different size-dependent thermal conductivities calculated using different reservoirs, even though they give the same bulk thermal conductivity. Specifically, the Langevin thermostat behaves like a thermalizing boundary in phonon BTE and provides mode-level thermal-equilibrium phonon outlets, while the Nose-Hoover chain thermostat and velocity rescaling method behave like biased reservoirs, which provide a spatially uniform heat generation and mode-level nonequilibrium phonon outlets. These findings explain why different experimental measurement methods can yield different size-dependent thermal conductivity. They also indicate that the thermal conductivity of materials can be tuned for various applications by specifically designing thermostats. Finally, the unification of NEMD and phonon BTE will largely facilitate the study of thermal transport in complex systems in the future by, e.g., replacing computationally unaffordable first-principles NEMD simulations with computationally less expensive spectral BTE simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Implicit and coupled fluid plasma solver with adaptive Cartesian mesh and its applications to non-equilibrium gas discharges

In this work, we present a new fluid plasma solver with adaptive Cartesian mesh (ACM) based on a full-Newton (nonlinear, implicit) scheme for non-equilibrium gas discharge plasma. The electrons and ions are described using drift-diffusion approximation coupled to Poisson equation for the electric field. The electron-energy transport equation is solved to account for electron thermal conductivity, Joule heating, and energy loss of electrons in collisions with neutral species. The rate of electron-induced ionization is a function of electron temperature and could also depend on electron density (important for plasma stratification). The ion and gas temperature are kept constant. The transport equations are discretized using a non-isothermal Scharfetter-Gummel scheme to resolve possible large temperature gradients in the sheaths. We demonstrate the new solver for simulations of direct current (DC) and radiofrequency (RF) discharges. The implicit treatment of the coupled equations allows using large time steps. The full-Newton method (FNM) enables fast nonlinear convergence at each time step, offering significantly improved simulation efficiency. We discuss the selection of time steps for solving different plasma problems. The new solver enables solving several problems we could not solve before with existing software: two- and three-dimensional structures of the entire DC discharges including cathode and anode regions, electric field reversals and double-layer formation, the normal cathode spot and an anode ring, moving striations in diffuse and constricted DC discharges, and standing striations in RF discharges. The developed FNM-ACM technique offers many benefits for tackling the disparity of gas discharge plasma systems' time scales and nonlinearity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗