Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Doubling method”

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

Structural mode coupling in perovskite oxides using hypothesis-driven active learning

Abstract Finding the ground-state structure with minimum energy is paramount to designing any material. In ABO 3 -type perovskite oxides with Pnma symmetry, the lowest energy phase is driven by an inherent trilinear coupling between the two primary order parameters such as rotation and tilt with antiferroelectric displacement of the A-site cations as established via hybrid improper ferroelectric mechanism. Conventionally, finding the relevant mode coupling driving phase transition requires performing first-principles calculations which is computationally time-consuming as well as expensive. It involves following an intuitive iterative hit and trial method of (a) adding two or multiple mode vectors, followed by (b) evaluating which combination would lead to the ground-state energy. In this study, we show how a hypothesis-driven active learning framework can identify suitable mode couplings within the Landau free energy expansion with minimal information on amplitudes of modes for a series of double perovskite oxides with A-site layered, columnar and rocksalt ordering. This scheme is expected to be applicable universally for understanding atomistic mechanisms derived from various structural mode couplings behind functionalities, for e.g. polarization, magnetization and metal–insulator transitions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Electron Scattering from 1-Methyl-5-Nitroimidazole: Cross-Sections for Modeling Electron Transport through Potential Radiosensitizers

In this study, we present a complete set of electron scattering cross-sections from 1-Methyl-5-Nitroimidazole (1M5NI) molecules for impact energies ranging from 0.1 to 1000 eV. This information is relevant to evaluate the potential role of 1M5NI as a molecular radiosensitizers. The total electron scattering cross-sections (TCS) that we previously measured with a magnetically confined electron transmission apparatus were considered as the reference values for the present analysis. Elastic scattering cross-sections were calculated by means of two different schemes: The Schwinger multichannel (SMC) method for the lower energies (below 15 eV) and the independent atom model-based screening-corrected additivity rule with interferences (IAM-SCARI) for higher energies (above 15 eV). The latter was also applied to calculate the total ionization cross-sections, which were complemented with experimental values of the induced cationic fragmentation by electron impact. Double differential ionization cross-sections were measured with a reaction microscope multi-particle coincidence spectrometer. Using a momentum imaging spectrometer, direct measurements of the anion fragment yields and kinetic energies by the dissociative electron attachment are also presented. Cross-sections for the other inelastic channels were derived with a self-consistent procedure by sampling their values at a given energy to ensure that the sum of the cross-sections of all the scattering processes available at that energy coincides with the corresponding TCS. This cross-section data set is ready to be used for modelling electron-induced radiation damage at the molecular level to biologically relevant media containing 1M5NI as a potential radiosensitizer. Nonetheless, a proper evaluation of its radiosensitizing effects would require further radiobiological experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Characterization of a SiPM-based monolithic neutron scatter camera using dark counts

The Single Volume Scatter Camera (SVSC) Collaboration aims to develop portable neutron imaging systems for a variety of applications in nuclear non-proliferation. Conventional double-scatter neutron imagers are composed of several separate detector volumes organized in at least two planes. A neutron must scatter in two of these detector volumes for its initial trajectory to be reconstructed. As such, these systems typically have a large footprint and poor geometric efficiency. We report on the design and characterization of a prototype monolithic neutron scatter camera that is intended to significantly improve upon the geometrical shortcomings of conventional neutron cameras. The detector consists of a 50 mm×56 mm× 60 mm monolithic block of EJ-204 plastic scintillator instrumented on two faces with arrays of 64 Hamamatsu S13360-6075PE silicon photomultipliers (SiPMs). The electronic crosstalk is limited to < 5% between adjacent channels and < 0.1% between all other channel pairs. SiPMs introduce a significantly elevated dark count rate over PMTs, as well as correlated noise from after-pulsing and optical crosstalk. In this article, we characterize the dark count rate and optical crosstalk and present a modified event reconstruction likelihood function that accounts for them. We find that the average dark count rate per SiPM is 4.3 MHz with a standard deviation of 1.5 MHz among devices. The analysis method we employ to measure internal optical crosstalk also naturally yields the mean and width of the single-electron pulse height. Here, we calculate separate contributions to the width of the single-electron pulse-height from electronic noise and avalanche fluctuations. We demonstrate a timing resolution for a single-photon pulse to be (128 ± 4) ps. Finally, coincidence analysis is employed to measure external (pixel-to-pixel) optical crosstalk. We present a map of the average external crosstalk probability between 2×4 groups of SiPMs, as well as the in-situ timing characteristics extracted from the coincidence analysis. Further work is needed to characterize the performance of the camera at reconstructing single- and double-site interactions, as well as image reconstruction.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Suppressing CO formation in low-temperature methanol steam reforming via Ce-modified CuZnGa layered oxide catalysts

Cu-based layered double hydroxides (LDHs) are widely recognized as effective catalysts for low-temperature methanol steam reforming, yet achieving high hydrogen productivity together with near-complete suppression of CO formation remains challenging. Here, we report the synthesis and evaluation of a series of CuZnGa LDH-derived catalysts and Ce-modified analogues prepared via an aqueous miscible organic method, which enables high metal dispersion and precise structural control. The optimized CuZnGa catalyst exhibits a hydrogen production rate of 16.9 µmol H 2 ·g cat −1 ·s −1 at 180 °C with an H 2 /CO ratio exceeding 3500, outperforming many state-of-the-art low-temperature systems. Importantly, the incorporation of small amounts of Ce further suppresses CO formation while maintaining high hydrogen productivity. Combined spectroscopic characterization and density functional theory calculations reveal that Ce is incorporated into the LDH lattice by substituting Ga 3+ sites up to a critical threshold, beyond which highly dispersed CeO x species are formed. These species provide mobile lattice oxygen that participates in a Mars-van Krevelen-type pathway, selectively oxidizing CO and suppressing the reverse water-gas shift reaction. This study establishes a clear relationship between Ce speciation, oxygen mobility, and catalytic selectivity in LDH-derived systems. The resulting catalysts demonstrate the potential of interface-engineered Cu-based materials for efficient low-temperature hydrogen production with minimal CO contamination.

09 BIOMASS FUELS↗
DOE OSTI2022

Catalyst layer formulations for slot-die coating of PEM fuel cell electrodes

The series of electrodes were fabricated by the scalable and manufacturable slot-die coating method for proton exchange membrane fuel cell (PEMFC) application. The inks with different amounts of solids were studied by rheological methods in order to establish a coating window with minimum manufacturing defects. The obtained electrodes were characterized by SEM, AFM, and optical microscopy, which showed that they were uniform and homogeneous with minimum defects. The electrochemical evaluation of the manufactured gas diffusion electrodes (GDE) showed that the main characteristics of the electrodes, like electrochemical surface area, proton resistivity, and double layer capacitance, were found to be close for all samples confirming the reproducibility of the slot-die process. Additionally, we studied the effects of membrane thickness on the performance of the GDE membrane electrode assemblies and determined that a decrease in membrane thickness favored the performance. Furthermore, the obtained results clearly demonstrated the applicability and feasibility of the approach for the Manufacturing of catalyst layers for the fuel cell application with potential for future mass production.

08 HYDROGEN↗
DOE OSTI2020

$t_1$-Noise eliminated dipolar heteronuclear multiple-quantum coherence solid-state NMR spectroscopy

Heteronuclear correlation (HETCOR) spectroscopy is one of the key tools in the arsenal of the solid-state NMR spectroscopist to probe chemical and spatial proximities between two different nuclei and enhance spectral resolution. Dipolar heteronuclear multiple-quantum coherence (D-HMQC) is a powerful technique that can be potentially utilized to obtain 1 H detected 2D HETCOR solid-state NMR spectra of any NMR active nucleus. A long-standing problem in 1 H detected D-HMQC solid-state NMR experiments is the presence of t 1 -noise which reduces sensitivity and impedes spectral interpretation. In this contribution, we describe novel pulse sequences, termed t 1 -noise eliminated (TONE) D-HMQC, that minimize t 1 -noise and can provide higher sensitivity and resolution than conventional D-HMQC. Monte-Carlo and numerical simulations confirm that t 1 -noise in conventional D-HMQC primarily occurs because random MAS frequency fluctuations cause variations in the NMR signal amplitude from scan to scan, leading to imperfect cancellation of uncorrelated signals by phase cycling. The TONE D-HMQC sequence uses 1 H π-pulses to refocus the evolution of 1 H CSA across each SR4 2 1 recoupling block, improving the stability of the pulse sequence to random MAS frequency fluctuations. The 1 H refocusing pulses also restore the orthogonality of in-phase and anti-phase magnetization for all crystallite orientations at the end of each recoupling block, enabling the use of 90° flip-back or LG spin-lock trim pulses to reduce the intensity of uncorrelated signals. We demonstrate the application of these methods to acquire 1 H detected 2D 1 H{ 35 Cl} and 1 H{ 13 C} HETCOR spectra of histidine·HCl·H 2 O with reduced t 1 -noise. To show generality, we also apply these methods to obtain 2D 1 H{ 17 O} spectra of 20%- 17 O fmoc-alanine and for the first time at natural abundance, 2D 1 H{ 25 Mg} HETCOR spectra of magnesium hydroxide. The TONE D-HMQC sequences are also used to probe 1 H– 25 Mg and 1 H– 27 Al proximities in Mg–Al layered double hydroxides and confirm the even mixing of Mg and Al in these materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

ASMT-783246, Assessment Report for Department 635 Analysis of Fiscal Year 2021 ES&H Events

This assessment analyzes Environment, Safety, and Health (ES&H) occurrences and Non-Occurrence Trackable Events (NOTEs) from fiscal year (FY) 2021. For this report, assessors used three primary methods for categorizing occurrence and NOTE data: issue categorization, DOE reporting criteria groups, and DOE cause codes. The FY 2021 Q1 occurrence and NOTE total was the lowest since this type of analysis began in FY 2018 Q4, following a downward trend from the FY 2019 Q3 high point. The FY 2021 Q2 occurrence and NOTE total was nearly double the FY 2021 Q1 total; occurrence totals declined slightly in Q3 and Q4, while NOTE totals remained the same. NOTEs in each of the final three quarters were more than double the amount from Q1. This assessment resulted in one observation. As COVID-19 vaccination rates increase and COVID 19 impacts on operations decrease, the number of workers on-site and the amount of activity-level work will increase. With these changes, focused attention on the following areas related to work planning and controls may reduce the probability of future events, hazard identification and analysis, compliance with standards, formality of operations, and job scoping.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

Ultradiffuse Galaxies (UDGs) with Hyper Suprime-Cam. I. Revised Catalog of Coma Cluster UDGs*

This is the first in a series of papers on the properties of ultradiffuse galaxies (UDGs) in clusters of galaxies. We present an updated catalog of UDGs in the Coma Cluster using g- and r-band images obtained with Hyper Suprime-Cam (HSC) of the Subaru telescope. We develop a method to find UDGs even in the presence of contaminating objects, such as halos and background galaxies. This study expands upon our previous works that covered about half the area of the Coma Cluster. The HSC observations covered the whole Coma Cluster up to the virial radius and beyond (an area twice as large as the previous studies) and doubled the numbers of UDGs (r eff , r ≥ 1.5 kpc) and sub-UDGs (1.0 ≤ r eff , r < 1.5 kpc) to 774 and 729, respectively. The new UDGs show internal properties consistent with those of previous studies (e.g., a Sérsic index of approximately 1), and are distributed across the cluster, with a concentration around the cluster center. The whole cluster coverage clearly revealed an excess of their distribution toward the east to southwest direction along the cluster center, where Coma connects to other large-scale structure, and where a known substructure exists (the NGC 4839 subgroup). The alignment of the UDG distribution along the large-scale structure around Coma supports the interpretation that most of them lie at the distance of the Coma Cluster and the NGC 4839 subgroup.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Numerical modelling of impact seismic sources using the stress glut theory

SUMMARY Meteorite impacts have proved to be a significant source of seismic signal on the Moon, and have now been recorded on Mars by InSight seismometers. Understanding how impacts produce seismic signal is key to the interpretation of this unique data, and to improve their identification in continuous seismic records. Here, we use the seismic Representation Theorem, and particularly the stress glut theory, to model the seismic motion resulting from impact cratering. The source is described by equivalent forces, some resulting from the impactor momentum transfer, and others from the stress glut, which represents the mechanical effect of plasticity and non linear processes in the source region. We condense these equivalent forces into a point-source with a time-varying single force and nine-component moment tensor. This analytical representation bridges the gap between the complex dynamics of crater formation, and the linear point-source representation classically used in seismology. Using the multiphysics modelling software HOSS, we develop a method to compute the stress glut of an impact, and the associated point-source from hypervelocity impact simulations. For a vertical and an oblique impact at 1000 m s−1, we show that the moment tensor presents a significant deviatoric component. Hence, the source is not an ideal isotropic explosion contrary to previous assumptions, and draws closer to a double couple for the oblique impact. The contribution of the point force to the seismic signal appears negligible. We verify this model by comparing two signals: (1) HOSS is coupled to SPECFEM3D to propagate the near-source signal elastically to remote seismic stations; (2) the point-source model derived from the stress-glut theory is used to generate displacements at the same distance. The comparison shows that the point-source model is accurately simulating the low-frequency impact seismic waveform, and its seismic moment is in trend with Lunar and Martian impact data. High-frequencies discrepancies exist, which are partly related to finite-source effects, but might be further explained by the difference in mathematical framework between classical seismology and HOSS’ numerical modelling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Time propagation of electronic wavefunctions using nonorthogonal determinant expansions

The use of truncated configuration interaction in real-time time-dependent simulations of electron dynamics provides a balance of computational cost and accuracy, while avoiding some of the failures associated with real-time time-dependent density functional theory. However, low-order truncated configuration interaction also has limitations, such as overestimation of polarizability in configuration interaction singles, even when perturbative doubles are included. Increasing the size of the determinant expansion may not be computationally feasible, and so, in this work, we investigate the use of nonorthogonality in the determinant expansion to establish the extent to which higher-order substitutions can be recovered, providing an improved description of electron dynamics. Model systems are investigated to quantify the extent to which different methods accurately reproduce the (hyper)polarizability, including the high-harmonic generation spectrum of H 2 , water, and butadiene.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Determining the Solubility Behavior of Kogarkoite in Simulated Nuclear Waste

Kogarkoite (Na 3 FSO 4 ) is a sparingly soluble fluoride–sulfate double salt that has been identified in high level nuclear waste sludge at the Hanford Site and, more recently, in sludge batch compilation samples at the Savannah River Site (SRS). Due to its complex dissolution behavior, which exhibits an inverse dependence on sodium ion activity, the presence of this mineral poses significant challenges to waste retrieval and processing. Incomplete dissolution during sludge washing can lead to the retention of fluoride and sulfate in the high-level waste feed, potentially causing the formation of corrosive, immiscible molten salt layers, known as "glass gall,” in vitrification melters. Current efforts to optimize flowsheet parameters and wash-water volumes are hindered by the absence of a commercially available, certified reference material, which prevents the accurate calibration of analytical methods and the verification of dissolution kinetics. To address this critical gap, this research focuses on the laboratory synthesis of pure Kogarkoite to serve as a standard for comprehensive solubility and washing performance testing. A coupled synthesis and simulant campaign was executed using an evaporative crystallization protocol designed to replicate the dynamic concentration effects observed in tank farm operations. Thirteen simulant matrices were prepared by dissolving systematically varied ratios of sodium fluoride (NaF) and sodium sulfate (Na 2 SO 4 ) in deionized water under three distinct caustic regimes: 0.0 g (control), 4.0 g (~1 M), and 12.0 g (~3 M) sodium hydroxide (NaOH). While thermodynamic equilibrium models suggest that high-caustic environments should favor the stability of the double salt7, results from this evaporative study at 25 0 C revealed a distinct kinetic divergence. Simulants with high hydroxide loading predominantly yielded large, blocky crystals of sodium sulfate decahydrate (Na 2 SO 4 .10H 2 O). Successful synthesis of pure Kogarkoite was achieved exclusively in specific NaOH-free compositional windows, where the precipitate manifested as fine, opaque granular aggregates. Ion chromatography (IC) analysis confirmed phase purity through the simultaneous stoichiometric depletion of both fluoride and sulfate from the supernatant. This successful synthesis establishes a reproducible route to generate bulk Kogarkoite, enabling the subsequent phase of quantitative dissolution testing using inhibited water to optimize sludge-batch assembly.

Sarker, Md Sharif [Florida International Univ. (FI↗
DOE OSTI2023

Hydrogen-induced degradation dynamics in silicon heterojunction solar cells via machine learning

Abstract Among silicon-based solar cells, heterojunction cells hold the world efficiency record. However, their market acceptance is hindered by an initial 0.5% per year degradation of their open circuit voltage which doubles the overall cell degradation rate. Here, we study the performance degradation of crystalline-Si/amorphous-Si:H heterojunction stacks. First, we experimentally measure the interface defect density over a year, the primary driver of the degradation. Second, we develop SolDeg, a multiscale, hierarchical simulator to analyze this degradation by combining Machine Learning, Molecular Dynamics, Density Functional Theory, and Nudged Elastic Band methods with analytical modeling. We discover that the chemical potential for mobile hydrogen develops a gradient, forcing the hydrogen to drift from the interface, leaving behind recombination-active defects. We find quantitative correspondence between the calculated and experimentally determined defect generation dynamics. Finally, we propose a reversed Si-density gradient architecture for the amorphous-Si:H layer that promises to reduce the initial open circuit voltage degradation from 0.5% per year to 0.1% per year.

14 SOLAR ENERGY↗
DOE OSTI2020

Structure of the (110) Ln ScO 3 ( Ln =Gd,Tb,Dy) surfaces

The surface reconstruction on the (110) surface of LnScO 3 (Ln=Gd,Tb,Dy) has been investigated using x-ray photoelectron spectroscopy, atomic force microscopy, and transmission electron diffraction coupled with higher-level density-functional theory methods. The experimental techniques generate constraints on the surface chemistry and structure that are used to inform the theoretical comparison of a number of potential surface structures. The resulting structure is then compared to experimental aberration-corrected transmission electron microscopy results using multislice simulation. The surfaces of both single-crystalline substrates and hydrosauna synthesized nanoparticles exhibit a Sc-rich double layer with a Sc 3 O 4 termination.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Interpreting Ionic Conductivity for Polymer Electrolyte Fuel Cell Catalyst Layers with Electrochemical Impedance Spectroscopy and Transmission Line Modeling

A cathode catalyst layer containing optimally distributed ionomer is critical to reduce the platinum loading and increase its utilization in polymer electrolyte fuel cells. Here, electrochemical impedance spectroscopy (EIS) was used to measure effective ionic conductivity of pseudo catalyst layers (PCLs) at a relative humidity (RH) range of 50%-120%. These results are compared to previous work using the hydrogen pump (HP) method. EIS effective ionic conductivity results reported here are higher than those from the HP because in the HP set-up ionic pathways must be effectively connected through the PCL to be counted, whereas in the EIS measurement, ionomer segments that are in contact with the membrane but are not effectively connected all the way through the PCL can be detected. Double layer capacitances and effective ionic conductivities of Pt/C catalyst layers with various supports and ionomer to carbon (I/C) ratios were studied. High surface area carbon support resulted in a lower effective ionic conductivity compared to the graphitized carbon support due to worse ionomer dispersion. Effective ionic conductivities of Pt/C layers were compared to that of PCLs. On average, effective ionic conductivities of Pt/C layers were higher than PCLs because of possible carbon agglomeratio n within the PCLs.

25 ENERGY STORAGE↗
DOE OSTI2020

Functionalized Polymer Thin Films for Plutonium Capture and Isotopic Screening from Aqueous Sources

The rapid screening of plutonium from aqueous sources remains a critical challenge for nuclear nonproliferation efforts. The determination of trace-level Pu isotopes in water requires offsite sample preparation and analysis; therefore, new methods that combine plutonium purification, concentration, and isotopic screening in a fieldable detection system will provide an invaluable tool for nuclear safeguards. This contribution describes the development and characterization of thin polymer-ligand films for the isolation and concentration of waterborne Pu for direct spectroscopic analyses. Submicron thin films were prepared through spin coating onto Si wafers and consisted of combinations of polystyrene (PS) with dibenzoylmethane, thenoyltrifluoroacetone, and di(2-ethylhexyl)phosphoric acid (HDEHP). Pu uptake studies from solutions at pH from 2.3 to 6.3 indicated that only films containing HDEHP exhibited significant recovery of Pu. High alpha spectroscopy peak energy resolutions were achieved for PS-HDEHP films over a range of film thicknesses from 30 to 250 nm. Furthermore, a separate study was performed to evaluate uptake from a primarily Pu(V) solution where it was observed that doubling the HDEHP loading in the film increased uptake of Pu by an order of magnitude. X-ray photoelectron spectroscopy (XPS) analysis revealed that HDEHP was highly concentrated within the first few nanometers of the film at the higher loading. XPS analysis also revealed that, in the presence of water, HDEHP was stripped from the surface layer of the film at circumneutral pH. While significant losses of ligand were seen in all samples, higher loadings of HDEHP resulted in measurable amounts of ligand retained after a 12-h soak in water. Findings of this study are being used to guide the development of thin-film composite membrane-based detection methods for the rapid, fieldable analysis of Pu in water.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Conformational Panorama and Chirality Controlled Structure–Energy Relationship in a Chiral Carboxylic Acid Dimer

Chirality recognition in dimers of tetrahydro-2-furoic acid (THFA) was studied in a conformer-specific manner using rotational spectroscopy and theoretical approaches. THFA shows a strong preference for the trans- over the cis-COOH configuration. Two drastically different scenarios are possible for the detectable (THFA) 2 : a kinetically preferred dimer bound by feeble interactions between two trans-THFAs or a thermodynamically favored dimer with a double hydrogen-bonded ring structure between two cis subunits. To identify the conformers responsible for the extremely dense rotational spectra observed, it was essential not only to locate several hundred homo/heterochiral (THFA) 2 minima in ab initio calculations but also to evaluate the energetic connectivities among the minima. The study further reveals an interesting chirality dependent structure-energy ordering relationship. A method for enantiomeric excess (ee) determination of THFA is presented using a recently proposed chiral self-tag approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Improving the Accuracy of Atomistic Simulations of the Electrochemical Interface

Atomistic simulation of the electrochemical double layer is an ambitious undertaking, requiring quantum mechanical description of electrons, phase space sampling of liquid electrolytes, and equilibration of electrolytes over nanosecond time scales. All models of electrochemistry make different trade-offs in the approximation of electrons and atomic configurations, from the extremes of classical molecular dynamics of a complete interface with point-charge atoms to correlated electronic structure methods of a single electrode configuration with no dynamics or electrolyte. Here, we review the spectrum of simulation techniques suitable for electrochemistry, focusing on the key approximations and accuracy considerations for each technique. Additionally, we discuss promising approaches, such as enhanced sampling techniques for atomic configurations and computationally efficient beyond density functional theory (DFT) electronic methods, that will push electrochemical simulations beyond the present frontier.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗