Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinetics”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

UniKP: a unified framework for the prediction of enzyme kinetic parameters

Prediction of enzyme kinetic parameters is essential for designing and optimizing enzymes for various biotechnological and industrial applications, but the limited performance of current prediction tools on diverse tasks hinders their practical applications. Here, we introduce UniKP, a unified framework based on pretrained language models for the prediction of enzyme kinetic parameters, including enzyme turnover number (k cat ), Michaelis constant (K m ), and catalytic efficiency (k cat / K m ), from protein sequences and substrate structures. A two-layer framework derived from UniKP (EF-UniKP) has also been proposed to allow robust k cat prediction in considering environmental factors, including pH and temperature. In addition, four representative re-weighting methods are systematically explored to successfully reduce the prediction error in high-value prediction tasks. We have demonstrated the application of UniKP and EF-UniKP in several enzyme discovery and directed evolution tasks, leading to the identification of new enzymes and enzyme mutants with higher activity. UniKP is a valuable tool for deciphering the mechanisms of enzyme kinetics and enables novel insights into enzyme engineering and their industrial applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Design principles for energy transfer in the photosystem II supercomplex from kinetic transition networks

Photosystem II (PSII) has the unique ability to perform water-splitting. With light-harvesting complexes, it forms the PSII supercomplex (PSII-SC) which is a functional unit that can perform efficient energy conversion, as well as photoprotection, allowing photosynthetic organisms to adapt to the naturally fluctuating sunlight intensity. Achieving these functions requires a collaborative energy transfer network between all subunits of the PSII-SC. In this work, we perform kinetic analyses and characterise the energy landscape of the PSII-SC with a structure-based energy transfer model. With first passage time analyses and kinetic Monte Carlo simulations, we are able to map out the overall energy transfer network. We also investigate how energy transfer pathways are affected when individual protein complexes are removed from the network, revealing the functional roles of the subunits of the PSII-SC. In addition, we provide a quantitative description of the flat energy landscape of the PSII-SC. We show that it is a unique landscape that produces multiple kinetically relevant pathways, corresponding to a high pathway entropy. These design principles are crucial for balancing efficient energy conversion and photoprotection.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Revealing the growth kinetics of atmospheric corrosion pitting in aluminum via in situ microtomography

Abstract Understanding the mechanistic relationship between the environment, microstructure, and local kinetics of atmospheric corrosion damage remains a central challenge. To address this challenge, this study used laboratory-based X-ray tomography to directly observe attack in-operando over an extended period, enabling insights into the evolving growth kinetics and morphology of individual pits over months of exposure. Damage progression associated with nine pits in a 99.9% pure aluminum wire exposed to chloride salts in humid air was characterized. Most pits grew at a nominally linear rate up until pit death, which occurred within 12–24 h of nucleation. Exceptions to this were observed, with three pits exhibiting bimodal growth kinetics and growing for 40 or more hours. This was explained by secondary droplets that formed near the pits, increasing the cathode area. A corrosion-driven drying mechanism likely contributed to pit death in both cases. Pits first grew into the material followed by lateral expansion.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Data analytics approach to predict high-temperature cyclic oxidation kinetics of NiCr-based Alloys

Abstract Although of practical importance, there is no established modeling framework to accurately predict high-temperature cyclic oxidation kinetics of multi-component alloys due to the inherent complexity. We present a data analytics approach to predict the oxidation rate constant of NiCr-based alloys as a function of composition and temperature with a highly consistent and well-curated experimental dataset. Two characteristic oxidation models, i.e., a simple parabolic law and a statistical cyclic oxidation model, have been chosen to numerically represent the high-temperature oxidation kinetics of commercial and model NiCr-based alloys. We have successfully trained machine learning (ML) models using highly ranked key input features identified by correlation analysis to accurately predict experimental parabolic rate constants ( k p ). This study demonstrates the potential of ML approaches to predict oxidation kinetics of alloys over wide composition and temperature ranges. This approach can also serve as a basis for introducing more physically meaningful ML input features to predict the comprehensive cyclic oxidation behavior of multi-component high-temperature alloys with proper constraints based on the known underlying mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Pit growth kinetics in aluminum: effects of salt loading and relative humidity

Abstract The growth kinetics of localized corrosion, e.g. pits, in corrosive environments often controls the service life of metallic components. Yet, our understanding of these kinetics is largely based on coupon-level, e.g. mass-loss, studies which provide limited insights into the evolution of individual damage events. It is critical to relate observed cumulative loss trends, such as links between changing humidity and mass loss rates, to the growth kinetics of individual pits. Towards this goal, we leverage in-situ X-ray computed tomography to measure the growth rates of over sixty pits in aluminum in four different humid, chloride environments over ≈3 days of exposure. Pit growth rates and final volumes increased with increasing droplet volume, which was observed to increase with increasing humidity and salt loading. Two factors, droplet spreading and oxide jacking, dramatically increased pit growth rates and final volumes.

Noell, Philip J. (ORCID:0000000346998999)↗
DOE OSTI2023

Learning heterogeneous reaction kinetics from X-ray videos pixel by pixel

Reaction rates at spatially heterogeneous, unstable interfaces are notoriously difficult to quantify, yet are essential in engineering many chemical systems, such as batteries and electrocatalysts. Experimental characterizations of such materials by operando microscopy produce rich image datasets, but data-driven methods to learn physics from these images are still lacking because of the complex coupling of reaction kinetics, surface chemistry and phase separation. Here we show that heterogeneous reaction kinetics can be learned from in situ scanning transmission X-ray microscopy (STXM) images of carbon-coated lithium iron phosphate (LFP) nanoparticles. Combining a large dataset of STXM images with a thermodynamically consistent electrochemical phase-field model, partial differential equation (PDE)-constrained optimization and uncertainty quantification, we extract the free-energy landscape and reaction kinetics and verify their consistency with theoretical models. We also simultaneously learn the spatial heterogeneity of the reaction rate, which closely matches the carbon-coating thickness profiles obtained through Auger electron microscopy (AEM). Across 180,000 image pixels, the mean discrepancy with the learned model is remarkably small (<7%) and comparable with experimental noise. Our results open the possibility of learning nonequilibrium material properties beyond the reach of traditional experimental methods and offer a new non-destructive technique for characterizing and optimizing heterogeneous reactive surfaces.

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

Sulfate adenylyl transferase kinetics and mechanisms of metabolic inhibitors of microbial sulfate respiration

Abstract Sulfate analog oxyanions that function as selective metabolic inhibitors of dissimilatory sulfate reducing microorganisms (SRM) are widely used in ecological studies and industrial applications. As such, it is important to understand the mode of action and mechanisms of tolerance or adaptation to these compounds. Different oxyanions vary widely in their inhibitory potency and mechanism of inhibition, but current evidence suggests that the sulfate adenylyl transferase/ATP sulfurylase (Sat) enzyme is an important target. We heterologously expressed and purified the Sat from the model SRM, Desulfovibrio alaskensis G20. With this enzyme we determined the turnover kinetics (kcat, KM) for alternative substrates (molybdate, selenate, arsenate, monofluorophosphate, and chromate) and inhibition constants (KI) for competitive inhibitors (perchlorate, chlorate, and nitrate). These measurements enable the first quantitative comparisons of these compounds as substrates or inhibitors of a purified Sat from a respiratory sulfate reducer. We compare predicted half-maximal inhibitory concentrations (IC50) based on Sat kinetics with measured IC50 values against D. alaskensis G20 growth and discuss our results in light of known mechanisms of sensitivity or resistance to oxyanions. This analysis helps with the interpretation of recent adaptive laboratory evolution studies and illustrates the value of interpreting gene–microbe–environment interactions through the lens of enzyme kinetics.

Carlson, Hans K. (ORCID:0000000215835313)↗
DOE OSTI2020

Hydration structure and water exchange kinetics at xenotime–water interfaces: implications for rare earth minerals separation

Hydration of surface ions gives rise to structural heterogeneity and variable exchange kinetics of water at complex mineral–water interfaces. In this work, we employ ab initio molecular dynamics (AIMD) simulations and water adsorption calorimetry to examine the aqueous interfaces of xenotime, a phosphate mineral that contains predominantly Y 3+ and heavy rare earth elements. Consistent with natural crystal morphology, xenotime is predicted to have a tetragonal prismatic shape, dominated by the {100} surface. Hydration of this surface induces multilayer interfacial water structures with distinct OH orientations, which agrees with recent crystal truncation rod measurements. The exchange kinetics between two adjacent water layers exhibits a wide range of underlying timescales (5–180 picoseconds), dictated by ion–water electrostatics. Adsorption of a bidentate hydroxamate ligand reveals that {100} xenotime surface can only accommodate monodentate coordination with water exchange kinetics strongly depending on specific ligand orientation, prompting us to reconsider traditional strategies for selective separation of rare-earth minerals.

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

One-electron redox kinetics of aqueous transition metal couples Zn 2+/+ , Co 2+/+ , and Ni 2+/+ using pulse radiolysis

The one-electron redox potentials for aqueous metal couples Co 2+/+ and Ni 2+/+ have been investigated by using pulse radiolysis using their reactions with a series of reference compounds to establish the most positive upper limits of E 0 . Experiments with Zn + were also carried out to confirm the characteristic shape of the expected reduction kinetics. Both formate ions and t-BuOH were employed to scavenge ˙OH radicals and ˙H atoms. Kinetics and fitted first and second order reaction rates have been reported for reactions with methyl viologen, fluorescein, Ru(NH 3 ) 6 3+ , Co(en) 3 2+ , Co(sepulcrate) 3+ , Ru(bpy) 3 2+ , Cr(bpy)33+, and Ni(Me 6 [14]4,11-dieneN4) 2+ . Previous work demonstrated that both Co 2+ and Ni 2+ can be reduced by CO 2 ˙- radicals, giving a negative E 0 limit of -1.9 V vs. SHE. A definite reaction of Ni + with fluorescein di-anions provides a new upper limit of the Ni 2+/+ couple as -0.906 V vs. SHE. The reaction of Co + with Ru(bpy) 3 2+ has been confirmed, giving E 0 = -1.3 V vs. SHE as a rigorous upper limit of the Co 2+/+ couple. In the case of Co 2+/+ , kinetics were complicated by a self-catalyzed metal clustering phenomenon. Initiation rate constants of this process have also been reported.

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

Kinetic origins of the metastable zone width in the manganese oxide Pourbaix diagram

Pourbaix diagrams show solid–aqueous phase stability as a function of pH and redox potential and can be a valuable tool to guide the hydrothermal synthesis of transition metal oxides. However, the Pourbaix diagram is based on thermodynamics, and nucleation kinetics are not readily apparent in this framework. Here, we conduct a combined experimental and theoretical study to measure the onset of MnO 2 precipitation from an MnO 4 - (aq) solution at various pH conditions, which we then compare against Pourbaix diagram phase boundaries. Using a combination of in situ X-ray absorption spectroscopy and X-ray wide-angle scattering, we directly observe the transformation kinetics from the tetrahedral MnO 4 - (aq) ion to octahedrally coordinated Mn pre-nuclei, as well as its initiation into crystalline δ'-MnO 2 at various pH. The kinetics of octahedral-Mn precursor availability is observed to govern induction times of crystalline MnO 2 nucleation, which results in a metastable zone width in the Mn–H 2 O Pourbaix diagram. These results suggest that synthesis conditions often need to be prepared far beyond the phase boundaries in a Pourbaix diagram to initiate crystallization within a reasonable timescale.

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

Stereoisomer-dependent unimolecular kinetics of 2,4-dimethyloxetanyl peroxy radicals

2,4,dimethyloxetane is an important cyclic ether intermediate that is produced from hydroperoxyalkyl (QOOH) radicals in the low-temperature combustion of n-pentane. However, the reaction mechanisms and rates of consumption pathways remain unclear. In the present work, the pressure- and temperature-dependent kinetics of seven cyclic ether peroxy radicals, which stem from 2,4,dimethyloxetane via H-abstraction and O 2 addition, were determined. The automated kinetic workflow code, KinBot, was used to model the complexity of the chemistry in a stereochemically resolved manner and solve the resulting master equations from 300–1000 K and from 0.01–100 atm. The main conclusions from the calculations include (i) diastereomeric cyclic ether peroxy radicals show significantly different reactivities, (ii) the stereochemistry of the peroxy radical determines which QOOH isomerization steps are possible, (iii) conventional QOOH decomposition pathways, such as cyclic ether formation and HO 2 elimination, compete with ring-opening reactions, which primarily produce OH radicals, the outcome of which is sensitive to stereochemistry. Ring-opening reactions lead to unique products, such as unsaturated, acyclic peroxy radicals, that form direct connections with species present in other chemical kinetics mechanisms through "cross-over" reactions that may complicate the interpretation of experimental results from combustion of n-pentane and, by extension, other alkanes. For example, one cross-over reaction involving 1-hydroperoxy-4-pentanone-2-yl produces 2-(hydroperoxymethyl)-3-butanone-1-yl, which is an iso-pentane-derived ketohydroperoxide (KHP). At atmospheric pressure, the rate of chemical reactions of all seven peroxy radicals compete with that of collisional stabilization, resulting in well-skipping reactions. However, at 100 atm, only one out of seven peroxy radicals undergoes significant well-skipping reactions. Here, the rates produced from the master equation calculations provide the first foundation for the development of detailed sub-mechanisms for cyclic ether intermediates. In addition, analysis of the complex reaction mechanisms of 2,4-dimethyloxetane-derived peroxy radicals provides insights into the effects of stereoisomers on reaction pathways and product yields.

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

Kinetic processes of phase separation and aggregation behaviors in slot-die processed high efficiency Y6-based organic solar cells

Morphological optimization has proved to be one crucial factor contributing to the 19% efficiency of Y6-based organic solar cells (OSCs). Although the relationship between component miscibility and film morphology has been established, it has not been clarified how the film formation processes proceed, especially for methods compatible with large-area fabrication, i.e. slot-die coating, which restricts the further optimization of the morphology. Herein, we comprehensively investigate the effect of miscibility and the film-formation kinetic process on film morphology during processing with different solvents. A highest power conversion efficiency (PCE) of 17.38% can be obtained in a D18:Y6 device processed with slot-die coating in an open-air environment. However, due to their relatively poor miscibility, incorporating a trace amount of chlorobenzene (CB) into CF can sufficiently promote Y6 aggregation, leading to increased phase separation and thus the average PCE drops to 15.16%. In contrast, a PM6:Y6 blend shows insensitive changes in the kinetic process as well as the final morphology, and thus in the comparable PCEs when cast with different solvents, which is ascribed to their relatively good miscibility. As a result, this work provides scientific guidelines for device optimization by combining the intermediate gap from miscibility to kinetic process and their impact on final morphology.

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

Quantifying how the cis/trans ratio of N,N -dimethyl-3,5-dimethylpiperidinium hydroxide impacts the growth kinetics, composition and local structure of SSZ-39

This work integrates experiments and computational methods to quantify how the cis/trans ratio of the OSDA used in SSZ-39 synthesis impacts the crystallization kinetics, material properties, and final product composition. The crystallization kinetics increase by 30% when increasing the trans isomer content from 14% to 80%. Per prior work, in all cases based on the synthesis gel composition and product yield aluminum is the limiting reagent, and the absence of any amorphous material detected in the time resolved PXRD studies leads us to conclude that FAU dissolution is the rate limiting step in the formation of SSZ-39 in this synthesis protocol. The TGA and NMR results suggest that the trans isomer of OSDA is selectively incorporated into the product. The NMR binding studies, and corresponding DFT-based results show that the trans isomer binds to FAU more strongly than the cis isomer, providing one possible explanation for this enhancement in kinetics and preferential uptake of the trans isomer. The EDS analysis indicates that the Si/Al ratios are between 7.7 and 8.6 at low and high trans OSDA content, indicating zeolite composition is mildly sensitive to the trans isomer content. EDS results show this decrease in aluminum content leads to a corresponding decrease in sodium uptake. DFT-based calculations confirm OSDA–sodium interactions cannot explain any decrease in sodium uptake, reinforcing lower aluminum content as the cause of lower sodium uptake. Preliminary cobalt titration experiments show a surprisingly low cobalt uptake but also show a clear dependence of the cobalt uptake on the solution pH.

Cui, Zheng [Tulane University, New Orleans, LA (Un↗
DOE OSTI2025

Dynamics of inverse metal oxides on metal catalysts using spectro-kinetics: reversible Brønsted acid site formation and irreversible reduction

Brønsted acid sites (BASs) in inverse catalysts are vital for the selective hydrogenolysis of polyols, specifically cleaving secondary C–O bonds. These BASs form dynamically in situ in an H 2 environment. While H 2 enables rapid BAS generation on short timescales, it reduces the catalyst at prolonged exposures. The active center for BAS generation, the kinetics of BAS formation, its reverse decomposition, and the irreversible oxide reduction have lacked direct experimental evidence. Here, aided by advanced spectro-kinetic studies, we identify trimeric W 3 O x sites on Pt as the active centers for BAS generation, whereas isolated WO x species on SiO 2 act merely as spectator species, demonstrated using an inverse WO x /Pt catalyst as a representative system. A detailed kinetic profile capturing the dynamics of W 3 O x sites on Pt is also established. The rate constant for BAS formation is two orders of magnitude higher than for its decomposition, which is one order of magnitude faster than the irreversible site reduction. Co-fed H 2 O suppresses the site reduction by ∼50%. Furthermore, the H 2 partial pressure plays an important role. While lower gas-phase H 2 partial pressure does not influence the reversible BAS formation, it can significantly (∼3×) suppress catalyst reduction. Finally, these findings offer critical insights into optimizing reaction conditions through periodic H 2 pulsing, enhancing catalyst stability and performance in hydrogenolysis reactions.

Sourav, Sagar [Indian Inst. of Technology (IIT), M↗
DOE OSTI2025

Probing the proton exchange kinetics of BaZr 0.1 Ce 0.7 Y 0.1 Yb 0.1 O 3− δ ceramic electrolyte by operando diffuse reflectance infrared Fourier transform spectroscopy

Proton exchange kinetics plays an important role in governing the performance of intermediate-temperature protonic ceramic electrolysis cells (PCECs) for hydrogen production. Our understanding of the nature of the surface hydration reaction at the single-cell level, however, remains very limited, hampering further efficiency improvements. Here, in this study, we developed a custom operando diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) platform that operates under high temperature and steam conditions with applied bias. Quantitative investigations of surface H 2 O/D 2 O isotope exchange in a BaZr 0.1 Ce 0.7 Y 0.1 Yb 0.1 O 3−δ (BZCYYb1711) protonic electrolyte-based single cell were conducted under different applied voltages using this DRIFTS platform, to gain molecular-level insight into hydration kinetics. The findings show that the application of an external voltage significantly enhances the surface proton exchange rate, decreasing the apparent activation energy from 29.1 kJ mol −1 at open-circuit voltage (OCV) to 6.8 kJ mol −1 at 1.3 V. In addition, distinct voltage-induced spectral shifts in O–D vibrations point to dynamic changes in surface hydration. These findings demonstrate a sensitive spectroscopic platform for probing interfacial proton processes and reveal strong electrochemical control over surface proton kinetics, offering new opportunities for probing electrolyte hydration behavior in PCECs.

36 - MATERIALS SCIENCE↗
DOE OSTI2026

Solubility-limited depolymerization kinetics in the glycolysis of carbonyl-containing polymers

Chemical recycling of condensation polymers is often rationalized on the basis of the intrinsic reactivity of ester and carbonate functional groups. However, under heterogeneous conditions relevant to plastic waste processing and environmental degradation, bulk depolymerization rates often diverge from trends predicted by homogeneous chemistry. Here, we investigate how polymer–solvent compatibility, catalyst strength, and phase behavior govern the heterogeneous glycolysis of carbonyl-containing polymers. Using poly(ethylene terephthalate) (PET), glycol-modified PET (PETG), and bisphenol-A polycarbonate (PC) as model systems, we examine depolymerization kinetics at 180 °C with ethylene glycol and bisphenol A as diols under both amphoteric organosalt (TBD : MSA) and strong base (TBD) catalysis. Despite substantial differences in crystallinity and glycol uptake, PET and PETG depolymerize at comparable rates under organosalt catalysis, while PC depolymerizes significantly more slowly under identical conditions. Time-resolved molecular weight analysis and thermal characterization demonstrate that these rate differences do not arise from crystallinity, swelling, or inherent carbonyl reactivity, but instead reflect solubility-limited kinetics that constrain the transition from heterogeneous to homogeneous reaction regimes. When polymer solubility is low, depolymerization remains heterogeneous and slow; when solubility is enhanced—either through increased polymer–diol compatibility or stronger base catalysis—rapid homogeneous depolymerization is observed, reversing apparent reactivity trends. These results establish solubility and phase behavior as primary determinants of depolymerization kinetics in heterogeneous polymer recycling systems. By demonstrating how catalyst selection and solvent compatibility can expose or overcome solubility limitations, this work provides mechanistic insight to design more energy-efficient and selective chemical recycling processes. More broadly, these findings suggest that polymers with limited solvent or water compatibility may resist chemical degradation in the environment, favoring fragmentation and persistence as micro- and nanoplastics. Understanding solubility-controlled depolymerization offers a pathway toward more sustainable polymer design and end-of-life chemical recovery.

Watson-Sanders, Shelby [Department of Chemistry, U↗
DOE OSTI2020

Shock-driven kinetic and diffusive mix in high-Z pusher ICF designs

Revolver and Double Shell Inertial Confinement Fusion capsule designs hope to achieve a robust volumetric thermonuclear burn via the use of a high-Z pusher shell filled with a cryogenic D–T fuel. Unfortunately, mix of the pusher material into the fuel (gas) may adversely impact the burn performance. Hydrodynamic instability of the metal/gas interface as the mix source is an obvious concern, but 1D effects may also be detrimental. Such effects include plasma diffusion at material interfaces, which has been the subject of numerous theoretical, computational, and experimental investigations. However, other 1D mix mechanisms may exist, which have yet to be thoroughly explored. In particular, plasma kinetic effects may drive the mix when a shock breaks out of the metal/gas interface. Using the state-of-the-art, hybrid (kinetic-ion/fluid electron), multi-ion Vlasov–Fokker–Planck code, iFP, we show herein that shock-driven kinetic effects can reconfigure the interface and the interfacial width subsequently grows diffusively. Finally, we consider any implications for high-Z pusher designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Mean force kinetic theory applied to self-diffusion in supercritical Lennard-Jones fluids

A new kinetic-theory-based calculation of the self-diffusion coefficient for dense supercritical Lennard-Jones fluids is presented. The mean force kinetic theory, which was recently developed for transport in dense plasmas, is applied for the calculation of diffusion in dense neutral fluids. The calculation only requires the pair distribution function, a quantity that is readily calculable from equilibrium statistical mechanics for many systems, including the Lennard-Jones fluid. The self-diffusion coefficients are compared with calculations from molecular dynamics simulations, and good agreement at high density is demonstrated, even in the vicinity of the solid–fluid coexistence line. A comparison of different kinetic models with molecular dynamics simulations demonstrates that the transport coefficients have important contributions due to particle interaction via a potential of mean force and local correlations, which increase the collision rate. The new calculations compare well to those from free-volume theory and overcome a limitation of this theory that prevents its use in systems that interact via long range monotonic potentials. It is expected that this approach will also apply to other systems, including neutral-plasma and neutral-electrolyte mixtures.

79 ASTRONOMY AND ASTROPHYSICS↗