Search NASA⌕ Search

SEARCH · Search NASA

Results for “Solution thermodynamics”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Colloidal Behavior of Plutonium Oxide in Concentrated Electrolyte Solutions

The Hanford Site in Washington State manages legacy high-level radioactive waste streams that display major chemistry and engineering challenges, including the high salt levels and pH values that correspond to conditions under which many classical concepts describing chemical reactivity cannot be applied. One particular challenge that needs to be tackled at Hanford is that Pu concentrations, [Pu], in the soluble phases of the tank wastes are higher than expected based on the solubility of crystalline PuO2, which is widely accepted to be caused by the formation of PuO2 colloid, consisting of nano- to submicron-sized particles (PuO2 NPs). Fundamental research underpinning the behavior of PuO2 NPs under conditions not only relevant to the Hanford tank waste but at high ionic strength in general is needed to reliably predict the chemical reactivity of PuO2 NPs and develop engineering solutions to safely and efficiently process high-level radioactive waste into forms suitable for long-term storage. In this work, we study the behavior of PuO2 NPs (particle size ~100 nm) under high ionic strength conditions by reacting it with highly concentrated (up to 5 M) salt solutions. We explore different electrolyte compositions to elucidate the impact of different anions (NO3-, Cl-, ClO4-, SO42-, C2O42-, CO32-) on the stability of PuO2 NP in the acidic and alkaline pH regime. PuO2 NP aggregation and precipitation as function electrolyte concentration is tracked by a combination of liquid scintillation counting, dynamic light scattering for determination of particle size distributions, and zeta potentials as a proxy for particle charge. At acidic pH, electrolytes containing non-coordinating anions, such as NaNO3, NaCl, and NaClO4 mostly stabilize PuO2 NPs over a large electrolyte concentration range, showing only subtle differences in their reactivity. Other electrolyte anions show a more pronounced effect on the PuO2 NP stability: SO42-, binds directly to the particles’ surface, reverses the particle charge, and precipitates the PuO2 NPs efficiently even at intermediate sulfate concentrations (>0.1 M). In contrast, C2O42- is found to lead to high [Pu] in solution, in the milli-molar range, even at mildly acidic pH (~4). Thermodynamic modeling of the dissolved Pu concentrations using PHREEQC is unable to predict the observed [Pu] in the acidic pH regime, supporting the influence of colloids in maintaining elevated [Pu]. It is noteworthy that the current thermodynamic databases do not include constants for colloidal Pu phases and cannot accurately predict many of the high ionic strength solutions relevant to this work. The mechanisms and models responsible for these observations will need further investigation in the future. At high pH values (~12), PuO2 NPs exhibits classical sol-gel chemistry, meaning that upon destabilization of the colloidal sol, for example by addition of concentrated NaOH, highly porous and viscid PuO2 coagulates are formed that consist of a three-dimensional network likely held together by physical interactions. The PuO2 NP coagulate shows no significant reversibility of the aggregation when contacted with concentrated brines; however, PuO2 NPs can be efficiently resuspended in solution by addition of diluted electrolytes, alkaline solutions containing high amounts of carbonate, or simple addition of water. Especially carbonate is shown to stabilize PuO2 NPs in solution at high pH, characterized by stable colloidal suspensions that are resistant against sedimentation during centrifugation. Thermodynamic modeling of the carbonate system was able to predict an increasing dissolved Pu concentration with increasing carbonate concentration. However, the model was profoundly sensitive to the fixed redox potential and does not include any thermodynamic constants for colloidal Pu species.

Neumann, Julia↗
NASA NTRS2015

Thermochemistry of Silicates

The thermodynamic properties of vapor and condensed phases of silicates are crucial in many fields of science. These quantities address fundamental questions on the formation, stability, transformation, and physical properties of silicate minerals and silicate coating compositions. Here the thermodynamic activities of silica and other species in solid solution have been measured by the analysis of the corresponding high temperature vapors using Knudsen Effusion Mass Spectrometry (KEMS). In first set of experiments KEMS has been used to examine the volatility sequence of species (Fe, SiO, Mg, O2 and O) present in the vapor phase during heating of fosterite-rich olivine (Fo93Fa7) up to 2400 C and to measure the Fe, SiO and Mg activities in its solid solution. The data of fosterite-rich olivine are essential for thermochemical equilibrium models to predict the atmospheric and surface composition of hot, rocky exoplanets (Lava Planets). In the second set of experiments the measured thermodynamic activities of the silica in Y2O3-SiO2 and Yb2O3-SiO2 systems are used to assess their reactivity and degradation recession as environmental barrier coatings (EBCs) in combustion environments (e.g. non-moveable parts of gas turbine engine).

Silicates↗
DOE OSTI2020

Equilibrium solute segregation to matrix- θ' precipitate interfaces in Al-Cu alloys from first principles

Particular combinations of solute atoms segregated to the interface of the Al and θ ' - Al 2 Cu in Al-Cu alloys can help stabilize θ ' precipitates at high temperatures. Stabilization of such precipitates is determined by a combination of thermodynamics (including driving forces for coarsening and transformation and solute segregation tendencies) and kinetic effects (including solute diffusion and interfacial mobility in the presence of interfacial solute segregation). For some alloys such as recent Al-Cu-based alloys, multiple solutes segregate in significant quantities to interfaces, and solute-solute interactions at the interface are important, with multiple types of solutes competing for similar interfacial sites. To treat this situation, we develop and apply a statistical mechanics approach to calculate the temperature-dependent equilibrium solute atoms distribution near the coherent and semicoherent interfaces between the Al matrix and the θ ' - Al 2 Cu precipitates. The developed approach is applied to the investigation of Si, Mn, and Zr segregation at the interface, as particular combinations of these elements affect the thermal stability of the θ ' precipitates. We demonstrate that because Si and Mn atoms segregate on the same semicoherent interface, the presence of Si reduces the concentration of Mn solutes at the interface. Si atoms preferably occupy the first layer of the interface and compete with Mn atoms for one type of particular sites in the layer. Mn atoms preferably occupy the second layer of the semicoherent interface, and the Mn-Mn interaction plays an important role in their distribution. Zr atoms mostly segregate on one of the two nonequivalent sites of the second layer of the coherent interface. Due to symmetry properties of the coherent interface, the calculations show that the segregation Zr of atoms to this interface will likely lead to the formation of L 1 2 ordered Al 3 Zr layer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Synthesis and thermodynamics of uranium-incorporated α-Fe 2 O 3 nanoparticles

Hematite nanoparticles were synthesized with U(VI) in circumneutral water through a coprecipitation and hydrothermal treatment process. XRD, TEM, and EXAFS analyses reveal that uranium may aggregate along grain boundaries and occupy Fe sites within hematite. The described synthesis method produces crystalline, single-phase iron oxide nanoparticles absent of surface-bound uranyl complexes. EXAFS data were comparable to spectra from existing studies whose syntheses were more representative of naturally occurring, extended aging processes. Herein this work provides and validates an accelerated method of synthesizing uranium-immobilized iron oxide nanoparticles for further mechanistic studies. High temperature oxide melt solution calorimetry measurements were performed to calculate the thermodynamic stability of uranium-incorporated iron oxide nanoparticles. Increasing uranium content within hematite resulted in more positive formation enthalpies. Standard formation enthalpies of U x Fe 2–2x O 3 were as high as 76.88 ± 2.83 kJ/mol relative to their binary oxides, or -764.04 ± 3.74 kJ/mol relative to their constituent elements, at x = 0.037. Data on the thermodynamic stability of uranium retention pathways may assist in predicting waste uranyl remobilization, as well as in developing more effective methods to retain uranium captured from aqueous environments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Energetics of hydroxylbastnäsite solid solutions, La 1-$\chi$ Nd $\chi$ CO 3 OH

Bastnäsites (LnCO 3 (F,OH)) are a group of common rare earth elements (REE)-bearing minerals and are one of the primary global sources of REE. Due to the chemical similarities among REE, bastnӓsites tend to occur as solid solutions instead of end members in REE containing ores. To better understand the processes and the mechanisms of formation of such deposits, it is essential to determine the thermodynamic properties of bastnӓsites, including hydroxylbastnӓsite (LnCO 3 OH) solid solutions. In this work, we performed detailed structural and calorimetric investigations on synthetic hexagonal La–Nd hydroxylbastnӓsite (La 1-$\chi$ NdCO 3 OH, $\chi$ = 0, 0.25, 0.5, 0.75, 1) solid solutions. X-ray diffraction confirms the crystal structure of the solid solution series in the $\ P$6 space group, and pair distribution function (PDF) analysis reveals local bonding environments characterized by three different types of 9-coordinated metal-oxygen polyhedra. Unit cell parameters of La1–xNdxCO3OH exhibit a nearly linear relation with the Nd content x, suggesting a random distribution of La and Nd in the structure. Their standard enthalpies of formation (Δ$\ H$° f ) were determined by high temperature oxide melt drop solution calorimetry, from which the enthalpies of mixing (Δ$\ H$ mix ) were derived. The Δ$\ H$ mix can be fitted by a regular solution model with an interaction parameter of 12.58 ± 0.16 kJ/mol, suggesting enthalpic metastability of La1–xNdxCO3OH relative to the two endmembers. Combining entropy and enthalpy, we further estimated the Gibbs free energies of mixing (Δ$\ G$ mix ) at relevant temperatures, revealing favorable temperatures under which the intermediate La 1-$\chi$ NdCO 3 OH phases can be stabilized. Such entropy-driven stabilization, as is consistent with our geochemical modeling results, may explain the enhancement of thermal stability of the solid solutions in nature. Additionally, the temperature range constrained from this study may be used to estimate the thermal history of REE bastnӓsite deposit.

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

Mapping Interfacial Solution Forces that Drive Fluorescent Nanoparticle Incorporation into Crystals

Interfacial solution structure governs processes ranging from catalytic and electrochemical reactions to particle aggregation and crystallization. Here, we design bio-inspired nanocomposites by mapping interfacial forces and controlling nanoparticle incorporation during crystal growth. Our analytical model contains measured kinetic barriers for surface approaches and equilibrium binding constants with the crystal surface. We validate this model using fluorescent silica nanoparticles and calcite. Our results show that surface chemistry dictates incorporation pathways: methoxy groups are least favorable, hydroxyls reduce kinetic barriers by interacting favorably with hydration layers, carboxylates bind strongly but must overcome a kinetic barrier, and amines outperform all kinetically and thermodynamically. This framework resolves the interplay between interfacial solution structures, particle forces and dynamics, and growth kinetics, enabling predicative control of nanocomposite formation and multiplexed mixed-particle systems.

Zhang, Mingyi [Pacific Northwest National Laborato↗
NASA NTRS1993

Experimental partitioning of Zr, Ti, and Nb between silicate liquid and a complex noble metal alloy and the partitioning of Ti between perovskite and platinum metal

El Goresy et al.'s observation of Nb, Zr, and Ta in refractory platinum metal nuggets (RPMN's) from Ca-Al-rich inclusions (CAI's) in the Allende meteorite led them to propose that these lithophile elements alloyed in the metallic state with noble metals in the early solar nebula. However, Grossman pointed out that the thermodynamic stability of Zr in the oxide phase is vastly greater than metallic Zr at estimated solar nebula conditions. Jones and Burnett suggested this discrepancy may be explained by the very non-ideal behavior of some lithophile transition elements in noble metal solutions and/or intermetallic compounds. Subsequently, Fegley and Kornacki used thermodynamic data taken from the literature to predict the stability of several of these intermetallic compounds at estimated solar nebula conditions. Palme and Schmitt and Treiman et al. conducted experiments to quantify the partitioning behavior of certain lithophile elements between silicate liquid and Pt-metal. Although their results were somewhat variable, they did suggest that Zr partition coefficients were too small to explain the observed 'percent' levels in some RPMN's. Palme and Schmitt also observed large partition coefficients for Nb and Ta. No intermetallic phases were identified. Following the work of Treiman et al., Jurewicz and Jones performed experiments to examine Zr, Nb, and Ti partitioning near solar nebula conditions. Their results showed that Zr, Nb, and Ti all have an affinity for the platinum metal, with Nb and Ti having a very strong preference for the metal. The intermetallic phases (Zr,Fe)Pt3, (Nb,Fe)Pt3, and (Ti,Fe)Pt3 were identified. Curiously, although both experiments and calculations indicate that Ti should partition strongly into Pt-metal (possibly as TiPt3), no Ti has ever been observed in any RPMN's. Fegley and Kornacki also noticed this discrepancy and hypothesized that the Ti was stabilized in perovskite which is a common phase in Allende CAI's.

Jurewicz, Stephen R.↗
NASA NTRS1984

A thermodynamic model for subsolidus equilibria in the system CaO-MgO-Al2O3-SiO2

It is shown that the high temperature enthalpy of solution data for pure phases and solid solutions in the CMAS system are, for 'gabbroic' and 'peridotitic' compositions, consistent with available phase equilibrium data for the MAS, CAS and CMAS systems. A refined set of values of thermodynamic properties for these phases and solid solutions is tabulated. The small differences between the new data set and that of Helgeson et al. (1978) arise from new data on heat capacity and enthalpy being incorporated. The important constraints applied and the major difficulties in fitting which arose are summarized.

Wood, B. J.↗
DOE OSTI2020

Actinide N-Donor Thermodynamics: Expanding the f-element Covalency Dialogue. Final report

Resolving the chemistry and physics of f-electrons is one of the grand challenges of science for energy technology. A central component of this challenge is describing covalency in actinide-ligand interactions. Crystallography, x-ray spectroscopy and computational techniques have provided reasonable, and in some cases remarkable, support for covalency in actinide-ligand interactions. Courtesy the quality of the results obtained, a significant portion of the covalency dialogue has been forwarded through solid-state and computational studies of the more stable actinides (Th, U, Pu, Np). However; recent reports have questioned whether covalency in actinide interactions increases or decreases across the series. Examination of less-stable trans-actinides (Am, Cm, Bk, Cf and Es) by the aforementioned approaches can be rapidly limited due to material availability, radiological hazards and experimental data to validate computational models. Furthermore, limited thermodynamic data exists to confirm covalency in actinide solution phase interactions that are highly relevant to the remediation and reprocessing of used nuclear fuel. The best-defined interactions of actinides with soft donors in aqueous solutions involve (poly)aminopolycarboxylate (APC) ligands. The chelate effect and binding affinity with the acetic acid APCs subgroups encourages amine interactions with the actinide metal center. In the absence of these factors, amine interactions with actinides are too weak to overcome the protective hydration shell of the dissolved ion. The ability for APCs to force actinide interactions with soft nitrogen donors (as defined by Pearson’s Hard Soft Acid Base theory) encourages the application of these ligands in a variety of processes for actinide recovery from nearly chemically identical lanthanides. The ability to functionalize the amine center of the ligand in a variety of capacities (adding additional amine groups, exchanging the conventional acetate group for an acetate group, etcetera), allows for the APC ligand to serve as a thermodynamic probe for actinide-nitrogen interactions. Perhaps the single most significant breakthrough during the previous funding cycle was the observation that covalency for the transplutonium part of the actinide series can be increasingly influenced with energy degeneracy driven covalency as the actinides become heavier. This was observed most predominantly with dipicolinic acid, but extensions of f-orbital degeneracy were found to affect aliphatic aminopolycarboxylate-actinide complexes through einsteinium.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1981

An evaluation of computer codes for simulating the Galileo Probe aerothermal entry environment

The approaches of three computer flow field codes (HYVIS, COLTS, and RASLE), used to determine the Galileo Probe aerothermal environment and its effect on the design of the thermal protection system, are analyzed in order to resolve differences in their predicted results. All three codes account for the hypersonic, massively blown, radiation shock layers, characteristic of Jupiter entry. Significant differences, however, are evident in their solution procedures: the governing conservation equations, the numerical differencing methods, the governing physics (chemical, radiation, diffusion, and turbulence models), and the basic physical data (thermodynamic, transport, chemical, and spectral properties for atomic and molecular species). Solutions are compared for two near peak heating entry conditions for a Galileo Probe baseline configuration, having an initial mass of 242 kg and simulating entry into the Orton nominal atmosphere. The modern numerical methodology of COLTS and RASLE appear to provide an improved capability for coupled flow-field solutions.

Menees, G. P.↗
DOE OSTI2024

BioRT‐HBV 1.0: A Biogeochemical Reactive Transport Model at the Watershed Scale

Abstract Reactive Transport Models (RTMs) are essential tools for understanding and predicting intertwined ecohydrological and biogeochemical processes on land and in rivers. While traditional RTMs have focused primarily on subsurface processes, recent watershed‐scale RTMs have integrated ecohydrological and biogeochemical interactions between surface and subsurface. These emergent, watershed‐scale RTMs are often spatially explicit and require extensive data, computational power, and computational expertise. There is however a pressing need to create parsimonious models that require minimal data and are accessible to scientists with limited computational background. To that end, we have developed BioRT‐HBV 1.0, a watershed‐scale, hydro‐biogeochemical RTM that builds upon the widely used, bucket‐type HBV model known for its simplicity and minimal data requirements. BioRT‐HBV uses the conceptual structure and hydrology output of HBV to simulate processes including advective solute transport and biogeochemical reactions that depend on reaction thermodynamics and kinetics. These reactions include, for example, chemical weathering, soil respiration, and nutrient transformation. The model uses time series of weather (air temperature, precipitation, and potential evapotranspiration) and initial biogeochemical conditions of subsurface water, soils, and rocks as input, and output times series of reaction rates and solute concentrations in subsurface waters and rivers. This paper presents the model structure and governing equations and demonstrates its utility with examples simulating carbon and nitrogen processes in a headwater catchment. As shown in the examples, BioRT‐HBV can be used to illuminate the dynamics of biogeochemical reactions in the invisible, arduous‐to‐measure subsurface, and their influence on the observed stream or river chemistry and solute export. With its parsimonious structure and easy‐to‐use graphical user interface, BioRT‐HBV can be a useful research tool for users without in‐depth computational training. It can additionally serve as an educational tool that promotes pollination of ideas across disciplines and foster a diverse, equal, and inclusive user community.

Sadayappan, Kayalvizhi↗
DOE OSTI2021

Relating phase transitions to pore size distributions and mechanical mortar properties in CSA-OPC-C$ based systems – The potential impact of delayed straetlingite formation

The scope of this research is to investigate the performance evolution of five ternary CSA-OPC-Cs systems. Time resolved thermodynamic modeling was employed to gain insight into changes in the volume fractions of the hydrating pastes and link these changes to trends in the pore size distributions of the mortars. Pore solution analyses were performed to investigate which hydrates have the potential to exert high crystallization pressures. In three out of five blends the resulting tensile stress exceeded the tensile strength of the materials, resulting in the occurrence of coarse pores, which could be linked to the measured compressive strengths. Two damaging crystallization processes were identified in the investigated blends. One is the continuous precipitation of ettringite at high supersaturation in a system with high anhydrite addition. The other is the delayed crystallization of straetlingite at high supersaturation and crystallization pressure in the confined environment of mature pastes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Performance enhancement of aqueous ionic liquids with lower critical solution temperature (LCST) behavior through ternary mixtures

Thermally responsive ionic liquids (ILs) exhibit liquid–liquid phase separation when mixed with water and heated above a lower critical solution temperature (LCST), resulting in a water-rich (WR) and an IL-rich (ILR) phase. These binary IL–water mixtures can be employed in a variety of thermodynamic processes such as forward osmosis (FO) desalination, for which two solution properties are desirable: low phase separation temperature and high osmotic strength (osmolality). However, these two properties are interlinked, with ILs that exhibit higher osmotic strengths typically requiring higher phase separation temperatures. This behavior tends to arise from the hydrophilicity of the IL cations, which enhances osmotic strength while also elevating the phase separation temperature. In this work, we highlight a pathway to overcome this tradeoff by developing ternary IL mixtures (two ILs with varying cation hydrophilicity mixed with water), which lowers the phase separation temperature while maintaining and even enhancing the osmotic strength of the solution. We characterize the mixing behavior (osmolality, phase separation temperature, WR phase purity, and WR to ILR phase mass ratio) of four ILs as a function of their concentration in solution. We find that an enhancement of up to 81.6% in the osmolality with a concomitant reduction of up to 15.4% in the phase separation temperature can be achieved using this approach. The ternary mixture is also shown to improve the phase separation kinetics by nearly 95% compared to the binary mixture. Overall, this work highlights a new pathway to improve the performance of LCST ILs for water and energy applications.

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

Effect of particle size on thermodynamics and lithium ion transport in electrodes made of Ti 2 Nb 2 O 9 microparticles or nanoparticles

This study compares the charging mechanisms, thermodynamics, lithium ion transport, and operando isothermal calorimetry in lithium-ion battery electrodes made of Ti 2 Nb 2 O 9 microparticles or nanoparticles synthesized by solid-state or sol-gel methods, respectively. First, electrochemical testing showed that electrodes made of Ti 2 Nb 2 O 9 nanoparticles exhibited larger specific capacity, smaller polarization, and better capacity retention at large currents than those made of Ti 2 Nb 2 O 9 microparticles. Furthermore, potentiometric entropy measurements revealed that electrodes made of either Ti 2 Nb 2 O 9 microparticles or nanoparticles showed similar thermodynamics behavior governed by lithium intercalation in solid solution, as confirmed by in situ XRD measurements. However, electrodes made of Ti 2 Nb 2 O 9 nanoparticles featured smaller overpotential and faster lithium ion transport than those made of Ti 2 Nb 2 O 9 microparticles. In fact, operando isothermal calorimetry revealed smaller instantaneous and time-averaged irreversible heat generation rates at electrodes made of Ti 2 Nb 2 O 9 nanoparticles, highlighting their smaller resistive losses and larger electrical conductivity. Finally, the measured total heat generation over a charging/discharging cycle matched the measured net electrical energy loss. Overall, Ti 2 Nb 2 O 9 nanoparticles synthesized by the novel sol-gel method displayed excellent cycling performance and reduced heat generation as a fast-charging lithium-ion battery anode material. These features present major advantages for actual battery systems including larger energy and power densities, simpler thermal management, and enhanced safety.

25 ENERGY STORAGE↗
DOE OSTI2023

Ab initio Geochemistry of Hydrous Phases (Final Report)

Hydrous phases are among the most important Earth components. They are important for a broad suite of Earth processes, including the origin of life. From both thermodynamic and structural perspectives, however, they represent some of the most complex naturally occurring materials: their bonding often includes a combination of covalent solid, weak ionic, van der Waals, and hydrogen bonding, all within large unit cells. Most are solid solutions, and many are prone to variations in layer packing. Thermodynamic modeling of these materials is fundamental for understanding present and past natural processes, including those that shaped—and continue to shape—the structure and evolution of our planet. Yet, the thermodynamic properties of these materials at appropriate conditions are challenging to measure. To make significant progress and attain a deep understanding of these materials requires an atomistic theoretical approach.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Impact of Radiolysis on Iodine Speciation in a Variety of Matrices

This report details the work conducted in fiscal year 2023 with the aim of providing a better understanding of the speciation of iodine in support of domestic 99 Mo production Shine who uses a sulfate matrix in their process. Laboratory experiments to study the behavior of iodine speciation were conducted using spectroscopic methods including UV-vis and Raman as well as electrochemical studies. Radiolysis occurring in solution presents a unique challenge that is largely unaddressed with thermodynamic modeling, therefore a kinetic model was developed to address these issues. The focus of the studies conducted this year used surrogate radiolytic conditions induced by hydrogen peroxide for comparison to radiolytic conditions from direct gamma irradiation and neat acid solutions. For each case a model incorporating the kinetics was used to predict iodine speciation in the defined conditions where the output was directly compared to the experimental conditions. Further modeling was conducted incorporating the process conditions that are specific to Shine with incorporation of radiolysis and gas sparging.

07 ISOTOPE AND RADIATION SOURCES↗
NASA NTRS1970

Structure of the blast wave with ionization

Plane, cylindrical and spherical blast waves structure with ionization at local thermodynamic equilibrium analyzed using successive approximation to nonsimilar solution

Fu, T. S.↗