Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusivity”

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 199 records · Page 11

Effect of concentration dependence of the diffusion coefficient on homogenization kinetics in multiphase binary alloy systems

Diffusion calculations were performed to establish the conditions under which concentration dependence of the diffusion coefficient was important in single, two, and three phase binary alloy systems. Finite-difference solutions were obtained for each type of system using diffusion coefficient variations typical of those observed in real alloy systems. Solutions were also obtained using average diffusion coefficients determined by taking a logarithmic average of each diffusion coefficient variation considered. The constant diffusion coefficient solutions were used as reference in assessing diffusion coefficient variation effects. Calculations were performed for planar, cylindrical, and spherical geometries in order to compare the effect of diffusion coefficient variations with the effect of interface geometries. In most of the cases considered, the diffusion coefficient of the major-alloy phase was the key parameter that controlled the kinetics of interdiffusion.

Tenney, D. R.↗

Experiments on tandem diffusers with boundary-layer suction applied in between

Experiments were performed on conical diffusers of various configurations with the same, but rather unusually large, 16:1 area ratio. Because available performance data on diffusers fall short of very large area ratio configurations, an unconventional design, consisting of two diffusers following each other in tandem, was proposed. Both diffusers had the same area ratio of 4:1, but had different taper angles. While for the first diffuser (called leading) the angle remained constant, for the second (called follower), the taper angle was stepped up to higher values. Boundary layer control, by way of suction, was applied between the diffusers, and a single slot suction ring was inserted between them. The leading diffuser had an enclosed nominal divergence angle 2 theta = 5 degrees, while the follower diffusers had either 10, 20, 30, or 40 degrees, respectively, giving 4 combinations. The experiments were performed at four different Reynolds numbers with various suction rates. The rates indicate a general improvement in the performance of all diffusers with boundary layer suction. It appears that the improvement of the pressure recovery depends on both the Reynolds number and the suction rate, and the largest increase, 0.075, was found at the lowest R sub e when the follower divergence was 2 theta = 40 degrees.

Barna, P. S.↗

Characteristics of acoustic-gravity waves in a diffusively separated atmosphere

The effects of wave-induced diffusion on the characteristics of plane acoustic gravity waves in the earth's thermosphere are investigated by means of a two fluid model atmosphere in diffusive equilibrium. At wave periods greater than the shortest characteristic diffusion time for the minor gas, diffusion induced by the wave is considered to act to eliminate amplitude and phase differences between the wave-induced density fluctuations of individual gases and restore the perturbed fluid to diffusive equilibrium. Vertical diffusion is most important at large scale sizes, but horizontal diffusion dominates for horizontal wavelengths of less than several hundred kilometers. As a result of wave-induced diffusion, AE-C satellite measurements of neutral density fluctuations of thermospheric constituents are only compatible with relatively small scale internal gravity waves. Velocity and temperature differences for acoustic-gravity waves are greatest at periods near the mean collision time and diffusion time and lead to the dissipation of wave energy. The predictions of the theory are consistent with a high latitude, low altitude source for the medium scale waves and a more localized random source for the waves observed by AE-C.

Del Genio, A. D.↗

Ambipolar diffusion in the middle atmosphere

In the middle atmosphere above 60 km, the electron concentration increases with altitude, reaching values of 10 to the 10th per cu m in the daytime ionospheric E region near 100 km. The electrons are more mobile than the ions and diffuse more rapidly through the neutral atmosphere. The electron diffusion polarizes the medium, causing an electric field to develop that acts to retard the electron diffusion and enhance the conduction current of ions. A global zonally averaged numerical model of atmosheric electricity from the ground to 100 km is used to examine the effect of ambipolar diffusion and the earth's geomagnetic field on the currents and fields in the middle atmosphere. The results show that above about 65 km, ambipolar diffusion generates local electric fields and conduction currents that balance electron diffusion currents. The electric fields and conduction currents are a few orders of magnitude larger than the vertical fields and currents calculated from the downward mapping of the ionospheric potential without taking electron diffusion into account. Ambipolar diffusion does not alter the total current flowing in the global circuit. It is a local effect where enhanced conduction currens flow to balance the electron diffusion current.

Tzur, I.↗

Diffusion of a multi-species component and its role in oxygen and water transport in silicates

The diffusion of a multispecies component is complicated by the different diffusion coefficient of each species and the interconversion reactions among the species. A diffusion equation is derived that incorporates the diffusive fluxes of all species contributing to the component's concentration. The effect of speciation on diffusion is investigated experimentally by measuring concentration profiles of all species developed during diffusion experiments. Data on water diffusion in rhyolitic glasses indicate that H2O molecules predominate over OH groups as the diffusing species at very low to high water concentrations. A simple theoretical relationship is drawn between the effective total oxygen diffusion coefficient and the total water concentration of silicates at low water content.

Zhang, Youxue↗

Investigation of Diffusion Flame Tip Thermodiffusive and Hydrodynamic Instability under Microgravity Conditions

We employ the opposed flow flame-spread configuration in order to examine flame-front instability of diffusion flames near cold, solid boundaries. The thermo-diffusive and hydrodynamic instabilities can transform an initially planar flame front into an irregularly curved, corrugated, possibly fragmented front. Under ordinary 1-g conditions, the buoyancy-induced flow masks the thermo-diffusive and hydrodynamic instabilities and produces planar flames. Such stable spreading flames have been observed for decades in laboratory experiments. Experiments in zero gravity are necessary to produce unstable flame fronts. The thermo-diffusive/hydrodynamic microgravity instability appears in diffusion flames such as, for example: the candle flame oscillations observed by Dietrich et al.; smolder instabilities on a recent Space Shuttle flight. Drs. T. Kashiwagi and S. Olson have attributed the latter to a lowered oxygen transport rate to the hot, reactive surface. Consider a burning surface near the flame extinction limit. The flow, or stretch, induced by the diffusion flame is weak, hence buoyancy plays a small role, thereby enabling previously secondary mechanisms, such as differential thermo-diffusion, to become the most important mechanisms. The flame leading edge becomes unstable; and diffusion flame breakup, oscillation, and rejoining all occur at a measurable frequency of approximately O(1 Hz). This project has only begun in January of this year, 1999. To date, there have been no flight experiments on flame spread instabilities. However, we have made numerous experiments in the NASA 2.2 and 5 second drop towers on flame spread over very thin cellulosic fuels. We have been very fortunate through a combination of factors, to be explained, to obtain some interesting, perhaps even compelling, results on diffusion flame instability in the presence of heat losses to cold surfaces.

Wichman, I. S.↗

Temperature Dependence of Diffusivities in Liquid Elements (LMD)

This research was to advance the understanding of diffusion mechanisms in liquid metals and alloys through accurate diffusivity measurements over a wide range of temperatures, including the proximity of the materials melting points. Specifically, it was driven towards developing a methodology (and subsequent flight hardware) to enable several diffusion coefficient measurements (i.e., at several different temperatures) to be performed using a single sample. The Liquid Metal Diffusion (LMD) was funded as a Flight Definition Project in February 1993 in response to NRA 91-OSSA-20 (Microgravity Science and Applications Division). The Science Concept Review for LAID was held during April 1994. In January 1995 we were informed that we had failed this review and the project was change to ground-based activities only. A new proposal was submitted for the next NRA addressing the panels concerns. As part of NASA's Risk Mitigation program, a scaled-down version of the hardware was funded in July of 1995 for a flight opportunity utilizing experiment on the Microgravity Isolation Mount. This experiment was to determine the self-diffusivity of indium at 185 C. The LMD was transferred to the Mir Space Station in STS-81 and returned on STS-84 (January - May 1997). Three, out of five, self-diffusion data sets were returned. A description of this experiment/hardware is included below. This summary is only intended to give the reader an overview of the results obtained for the tasks outlined in the original proposal. Research that was not published is explained in more detail. At the end of this report is a list of refereed publications and invited talks that were given as a result of this work. The reader is directed to these for further details. Attachment: Real-time diffusivity measurements in liquids at several temperatures with one sample, On the insensitivity of liquid diffusivity measurements to deviations from 1D transport, and Numerical simulations of the convective contamination of diffusivity measurements in liquids.

Banish, R. Michael↗

Flow/Soot-Formation Interactions in Nonbuoyant Laminar Diffusion Flames

This is the final report of a research program considering interactions between flow and soot properties within laminar diffusion flames. Laminar diffusion flames were considered because they provide model flame systems that are far more tractable for theoretical and experimental studies than more practical turbulent diffusion flames. In particular, understanding the transport and chemical reaction processes of laminar flames is a necessary precursor to understanding these processes in practical turbulent flames and many aspects of laminar diffusion flames have direct relevance to turbulent diffusion flames through application of the widely recognized laminar flamelet concept of turbulent diffusion flames. The investigation was divided into three phases, considering the shapes of nonbuoyant round laminar jet diffusion flames in still air, the shapes of nonbuoyant round laminar jet diffusion flames in coflowing air, and the hydrodynamic suppression of soot formation in laminar diffusion flames.

Dai, Z.↗

Diffusion Flame Stabilization

Diffusion flames are commonly used for industrial burners in furnaces and flares. Oxygen/fuel burners are usually diffusion burners, primarily for safety reasons, to prevent flashback and explosion in a potentially dangerous system. Furthermore, in most fires, condensed materials pyrolyze, vaporize, and burn in air as diffusion flames. As a result of the interaction of a diffusion flame with burner or condensed-fuel surfaces, a quenched space is formed, thus leaving a diffusion flame edge, which plays an important role in flame holding in combustion systems and fire spread through condensed fuels. Despite a long history of jet diffusion flame studies, lifting/blowoff mechanisms have not yet been fully understood, compared to those of premixed flames. In this study, the structure and stability of diffusion flames of gaseous hydrocarbon fuels in coflowing air at normal earth gravity have been investigated experimentally and computationally. Measurements of the critical mean jet velocity (U(sub jc)) of methane, ethane, or propane at lifting or blowoff were made as a function of the coflowing air velocity (U(sub a)) using a tube burner (i.d.: 2.87 mm) (Fig. 1, left). By using a computational fluid dynamics code with 33 species and 112 elementary reaction steps, the internal chemical-kinetic structures of the stabilizing region of methane and propane flames were investigated (Fig. 1, right). A peak reactivity spot, i.e., reaction kernel, is formed in the flame stabilizing region due to back-diffusion of heat and radical species against an oxygen-rich incoming flow, thus holding the trailing diffusion flame. The simulated flame base moved downstream under flow conditions close to the measured stability limit.

Takahashi, Fumiaki↗

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively. The Mach Number distribution shows small yet distinct differences between the two cases such as locations of shocks/shock reflections and a slightly different impingement point on the wall of the diffuser from the expansion at the exit of the nozzle. Similarly the temperature distribution indicates different flow recirculation patterns in the test cell. Both cases capture all the essential flow phenomena such as the shock-boundary layer interaction, plume expansion, expansion of the first stage ejectors, mixing between the engine plume and the first stage ejector flow and pressurization due to the first stage ejectors. The final paper will discuss thermal loads on the walls of the diffuser and cooling mechanisms investigated.

Allgood, Daniel C.↗

Diffusion Flame Stabilization

Diffusion flames are commonly used for industrial burners in furnaces and flares. Oxygen/fuel burners are usually diffusion burners, primarily for safety reasons, to prevent flashback and explosion in a potentially dangerous system. Furthermore, in most fires, condensed materials pyrolyze, vaporize, and burn in air as diffusion flames. As a result of the interaction of a diffusion flame with burner or condensed-fuel surfaces, a quenched space is formed, thus leaving a diffusion flame edge, which plays an important role in flame holding in combustion systems and fire spread through condensed fuels. Despite a long history of jet diffusion flame studies, lifting/blowoff mechanisms have not yet been fully understood, compared to those of premixed flames. In this study, the structure and stability of diffusion flames of gaseous hydrocarbon fuels in coflowing air at normal earth gravity have been investigated experimentally and computationally. Measurements of the critical mean jet velocity (U(sub jc)) of methane, ethane, or propane at lifting or blowoff were made as a function of the coflowing air velocity (U(sub a)) using a tube burner (i.d.: 2.87 mm). By using a computational fluid dynamics code with 33 species and 112 elementary reaction steps, the internal chemical-kinetic structures of the stabilizing region of methane and propane flames were investigated. A peak reactivity spot, i.e., reaction kernel, is formed in the flame stabilizing region due to back-diffusion of heat and radical species against an oxygen-rich incoming flow, thus holding the trailing diffusion flame. The simulated flame base moved downstream under flow conditions close to the measured stability limit.

Takahashi, Fumiaki↗

Probing the Impact of Vacancy Diffusion on Void Dynamics at the Lithium Metal–Solid Electrolyte Interface

Lithium (Li) metal-based solid-state batteries (SSBs) are considered promising candidates for next-generation energy storage due to their superior energy density and enhanced safety compared to conventional Li-ion systems. However, their practical application is limited by challenges such as void formation at the Li-solid electrolyte (SE) interface, which disrupts ion transport and accelerates interfacial degradation. This work investigates how the coupled effects of electro-dissolution kinetics and surface diffusion at the Li metal surface govern the evolution of interfacial morphology during stripping. This work examines the influence of three distinct surface diffusion modes, which are terrace diffusion, step diffusion, and interlayer diffusion, on maintaining interfacial stability. In addition, how the dominant surface diffusion mechanism can overcome the contact loss due to high reaction kinetics is explored. Furthermore, the roughness of the Li metal anode surface is quantified, and the influence of different diffusion mechanisms on the evolution of the dynamic solid–solid interface is examined. The critical role of temperature in enhancing Li surface diffusivity and expanding the regime of stable contact is highlighted. By identifying distinct regimes of interface stability, this study analyzes how non-uniform electrochemical dynamics dictate void morphology evolution and interfacial contact. These insights offer guiding principles for engineering robust Li–SE interfaces in SSBs.

lithium metal anode↗

Cluster dynamics simulations of tritium and helium diffusion in lithium ceramics

Tritium (T) and He diffusion in LiAlO 2 and LiAl 5 O 8 phases influences the performance of tritium producing burnable absorber rods (TPBARs) by affecting the gas release, swelling and thermal conductivity of Li-bearing ceramic pellets. Frenkel pair defects and clusters created by irradiation can attract T and He interstitials and form clusters of the type He i x Li, He i x Al, He i x O, T i x Li, T i x Al, and T i x O, 1 ≤ x ≤ 4 in a Li, Al or O vacancy site (notation denotes x He or T atoms in a 1 Li, 1 Al or 1 O vacant site). The concentration and mobility of each of these clusters collectively contribute to the diffusion of the He and T gases in LiAlO 2 and LiAl 5 O 8 . Here, in this work, free energy cluster dynamics simulations implemented in the Centipede code, are used to obtain the concentration and diffusivities of these clusters which are then used to calculate the total diffusivity of T and He gases in LiAlO 2 and LiAl 5 O 8 . The results show that diffusivity of T is at least one order of magnitude higher in LiAlO 2 as compared to that in LiAl 5 O 8 whereas He diffusion is 2–13 orders of magnitude higher in LiAlO 2 as compared to that in LiAl 5 O 8 . There is a higher concentration of highly diffusive species (T interstitials and T i 03 Li for the case of tritium and He i 01 Li, He i 02 Li, and He i 03 Li for the case of He) in LiAlO 2 than in LiAl 5 O 8 which increase the total diffusion of T and He in LiAlO 2 .

36 MATERIALS SCIENCE↗

Diffusion-limited hydrolysis in polymeric materials

Polymer degradation under aggressive environmental stressors often develops heterogeneities due to diffusion-limited reaction phenomena. This is well established for diffusion-limited oxidation (DLO), which is known to occur for most polymeric materials at elevated temperatures but has been less summarized for the conditions of diffusion-limited hydrolysis (DLH). Here, an overview of hydrolysis for several materials and a computational model, analogous to the underlying equations for DLO, able to define this diffusion-reaction system is presented. A systematic study of the influence of various parameters, such as water diffusivity, reaction rate and order, and a more in-depth focus on residual isocyanate hydrolysis in a polymeric methylene diphenyl diisocyanate (pMDI) based polyurethane (PU) foam is given. For this system, we present experimental data for model ‘input’ parameters and discuss predictions for different conditions. We conceptually compare the behavior of diffusion-limited oxidation to that of diffusion-limited hydrolysis (DLH). With the mathematical framework and key material properties presented herein, any DLH phenomena following Fickian diffusion behavior can be understood, modeled, and predicted.

36 MATERIALS SCIENCE↗

Study of Tritium Diffusivity in Pure and Sn-Defective Zr: A First-Principles Density Functional Theory Approach

Zirconium alloys (e.g., zircaloy-4) are used as tritium ( 3 H) getter materials in tritium-producing burnable absorber rods (TPBARs) owing to their ability to capture 3 H and chemically convert them into metal hydrides. Understanding of 3 H diffusion mechanisms in zircaloy is crucial for the optimal design of material performance in nuclear technology. Here, for this work, we perform first-principles density functional theory calculations to study the 3 H diffusion mechanism in pure and impure Zr with a low concentration of tin (Sn) atoms to determine the impact of the presence of Sn on the movement of 3 H atoms through the material. First, we calculated the diffusion barriers for 3 H in pure Zr by taking different migration pathways. We then introduced a low concentration of Sn impurity and systematically explored the impurity effect on the diffusion barriers for 3 H. Using our calculated diffusion energy barriers, we further obtained the diffusion coefficients and analyzed the results by comparing them with the experimental and previously calculated values. A diffusion coefficient of the order of 10 –8 m 2 /s is predicted. We also found that the presence of a Sn impurity could reduce the diffusivity up to 4 orders of magnitude. In conclusion, our results could serve as guidelines for further experimental investigations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electrode and Microstructure Dependence of Oxygen Diffusion in Ferroelectric Hafnium Zirconium Oxide Thin Films

Hafnia-based ferroelectrics hold promise to reduce energy demand for computing by enabling compute-in-memory and as non-volatile memories. The ferroelectric phase in this material system is, in part, stabilized by oxygen vacancies. While oxygen vacancies may be a necessity for phase stability, they limit device endurance through diffusion and accumulation into conducting channels. Herein, it is shown that oxygen diffusion is spatially variable within individual grains of ferroelectric hafnium zirconium oxide (HZO). Using 18 O tracers and finite difference modeling, it is shown that grain boundaries and regions near electrode interfaces allow for relatively rapid oxygen diffusion, with values as much as 10 4 larger than the grain cores. Further, the selection of electrode material affects the diffusion coefficients across all microstructural regions. HZO films in contact with TiN electrodes result in more oxygen-deficient HZO films and higher oxygen diffusion coefficients. Tungsten electrodes result in fewer vacancies and lower diffusion coefficients. Diffusion activation energy differences between the HZO with the two electrodes is reconciled by differing populations of charged and uncharged oxygen vacancies. This insight into the local vacancy populations and diffusion pathways provides a platform for designing hafnia-based films, deposition processes, and integration strategies to reduce vacancy gradients and improve performance.

36 MATERIALS SCIENCE↗

Quantifying Anionic Diffusion in 2D Halide Perovskite Lateral Heterostructures

Anionic diffusion strongly impacts the stability of halide perovskite materials, but it is still not well understood. Here, a quantitative investigation of in-plane thermally driven anionic inter-diffusion in a series of novel 2D and quasi-2D halide perovskites lateral heterostructures is reported. The calculated diffusion coefficients (D) reveal the inhibition of Br–I inter-diffusion with bulky π-conjugated organic cations compared with short-chain aliphatic organic cations. Furthermore, halide diffusion is found to be faster in quasi-2D (n > 1) than 2D perovskites (n = 1). The increment becomes less apparent as the “n” number increases, akin to the quantum confinement effect observed for band gaps. These trends are rationalized by molecular dynamics simulations of free energy barriers for halide diffusion that reveal mechanisms for suppressing diffusion. We report this work provides important fundamental insights on the anionic migration and diffusion process in halide perovskite materials.

2D materials↗

Mutual Diffusion Coefficients and Refractive Index Increments of K 2 SO 4 (aq) at 298.15 K from Rayleigh Interferometry

Here, isothermal mutual diffusion coefficients (interdiffusion coefficients) were measured for K 2 SO 4 (aq) at 298.15 ± 0.005 K, at numerous concentrations ranging from dilute solutions to near saturation (0.59648 mol∙dm -3 ; 0.61349 mol∙kg -1 ) under free diffusion boundary conditions, using high precision Rayleigh interferometry. Under the experimental conditions these diffusion coefficients are on the volume-fixed reference frame D v . Two series of experiments were performed, the first using the traditional experimental approach with a mercury lamp source and with the interference patterns being recorded on glass photographic plates, and the second with a He-Ne laser light source and computerized data acquisition using a photodiode array. The results from both series of experiments are in excellent agreement, and generally yield diffusion coefficients precise to at least 0.002 x 10 -9 m 2 ∙s 1 (0.15% to 0.19%) and in most cases to 0.001 x 10 -9 m 2 ∙s -1 (0.07% to 0.1%). These experiments also yield accurate values of the refractive index differences for the solution pairs used in the diffusion experiments. The new diffusion coefficients are compared to two sets of published values of diffusion coefficients for K 2 SO 4 (aq) which are somewhat discrepant from each other. This study extends and complements our earlier work on the diffusion coefficients of the most common brine salts: NaCl(aq), KCl(aq), MgCl 2 (aq), CaCl 2 (aq), Na 2 SO 4 (aq), ] MgSO 4 (aq), NaHCO 3 (aq), and KHCO 3 (aq) at 298.15 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗