Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex solids”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Major element composition controls rare earth element solubility during leaching of coal fly ash and coal by-products

Coal combustion ash and pre-combustion coal refuse are currently under consideration as potential sources of rare earth elements. One of the early steps in recovering REEs from coal by-products is often acid leaching, which can result in low pH leachates with complex aqueous chemistry. The aim of this work was to understand the connection between REE solubility, pH, and major elemental components of leachates for coal by-products. To accomplish this, we investigated the effects of solids concentration (i.e., pulp density) and pH adjustment on REE solubility in acid leachates of coal fly ashes from the Powder River Basin (PRB) and Appalachian Basin in the United States, and a coal processing refuse from the Southwestern U.S. For PRB ashes, the concentrations of soluble REEs generally increased with increasing pulp density; however, at pulp density values above 80–100 g/L, the soluble REE concentrations in the leachates were markedly lower. Similarly, the soluble concentrations of other major solutes (Fe, Al, Si) that leached from PRB fly ashes were also non-linear with pulp density. These major elements tended to reach maximum concentration values at 60–70 g/L pulp density. In contrast, for the Appalachian fly ashes and the coal by-product, soluble concentrations of REE and major elements in leachates increased linearly with pulp density. Chemical equilibria calculations of mineral saturation indices indicated that trends in soluble REE concentrations could be explained by saturation conditions for Fe and Al-(hydr)oxides and possibly sulfate minerals, but not lanthanide hydroxides. Furthermore, pH adjustment of the acid leachates showed that REEs and many major solutes were removed from solution at pH values above 4.5, also consistent with Fe- and Al-(hydr)oxide precipitation. These results highlight the importance of understanding the chemical composition of leachates when designing REE recovery processes for low-grade geologic feedstocks and that precipitation of hydr(oxide) or sulfate minerals of major elements rather than discreet formation of REE mineral phases could be used for process optimization.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

Strategy for Nonenzymatic Harvesting of Cells via Decoupling of Adhesive and Disjoining Domains of Nanostructured Stimulus-Responsive Polymer Films

The nanostructured polymer film introduces a novel mechanism of nonenzymatic cell harvesting by decoupling solid cell-adhesive and soft stimulus-responsive cell-disjoining areas on the surface. The key characteristics of this architecture are the decoupling of adhesion from detachment and the impermeability to the integrin protein complex of the adhesive domains. This surface design eliminates inherent limitations of thermoresponsive coatings, namely, the necessity for the precise thickness of the coating, grafting or cross-linking density, and material of the basal substrate. The concept is demonstrated with nanostructured thermoresponsive films made of cell-adhesive epoxy photoresist domains and cell-disjoining poly(N-isopropylacrylamide) brush domains.

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

Use of vibrational spectroscopy to identify the formation of neptunyl–neptunyl interactions: a paired density functional theory and Raman spectroscopy study

Actinyl–Actinyl interactions (AAIs) occur in pentavalent actinide systems, particularly for neptunium (Np), and lead to complex vibrational signals that are challenging to analyze and interpret. Previous studies have focused on neptunyl–neptunyl dimeric species, but trimers and tetramers have been identified as the primary motif for extended topologies observed in solid-state materials. Our hypothesis is that trimeric and tetrameric AAIs lead to the additional signals in the vibrational spectra, but this has yet to be explored systematically. Herein, we investigate three different neptunyl–neptunyl subunits (dimeric, trimeric, tetrameric) and determine the vibrational frequencies of the O=Np=O stretches using both computational and experimental approaches. Density Functional Theory (DFT) was used to identify distinct vibrational motions related to specific neptunyl oligomers and compared to previous literature precedent from Np(V) in HClO 4 and HCl systems. The vibrational behavior of Np(V) in HNO 3 was then evaluated via Raman spectroscopy. As the solution evaporated signals were linked to trimeric and tetrameric models. Solid phases produced in the evaporation include (NpO 2 ) 2 (NO 3 ) 2 (H 2 O) 5 and newly identified crystalline phase, Na(NpO 2 )(NO 3 ) 2 ·4H 2 O (NpNa). Furthermore, the combined computational studies and vibrational analysis provide evidence for unique observable vibrational bands for each polymerized subunit, allowing us to assign spectral features to trimeric and tetrameric models within three different simple anionic systems.

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

Simulation of the High Temperature Test Facility (HTTF) Core Using the 2D Ring Model with SAM

A two-dimensional ring model is developed with SAM to model the core of the High Temperature Test Facility (HTTF) at the system level. The ring model simplifies the complex structure of the HTTF core by converting the hexagonal rows of heaters and flow channels into layers of concentric annular rings. The ring model is first compared against a 3D-1D model where the solid structures are fully resolved in three dimensions while the fluid structures are modeled as one-dimensional flows. Comparison between the 3D-1D and the ring models shows that the latter can predict major parameters reasonably well under steady-state normal operating condition, but the heater temperatures are underpredicted. Adjustment is made to the effective thermal conductivity of the ceramic core of the ring model to improve the heater temperature predictions. Furthermore, the ring model is also used to simulate a transient Pressurized Conduction Cooldown (PCC) condition and is benchmarked with the experimental data from the HTTF Test PG-27. Good agreement is obtained between the experimental data and the predictions by the ring model.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Assessment of exchange-correlation functionals on oxygen vacancy formation energies of metal oxides

Oxygen vacancies are a common defect in oxides and play vital roles in many technological applications such as oxygen separation, catalytic reactors, solid oxide fuel cells, and solar thermochemical water splitting. The oxygen vacancy formation energy is directly related to the oxide reduction enthalpy and the vacancy concentration and, therefore, is a key quantity underlying these applications. The complexity of measuring oxygen vacancy formation energies experimentally suggests the utility of calculations based on density functional theory (DFT). However, the calculated results are strongly dependent on the exchange-correlation functionals and other parameters such as the Hubbard U . In this work, we compare the performance of the meta-GGA strongly constrained and appropriately normed (SCAN) functional and the commonly used semilocal generalized gradient approximation (GGA) functionals with experimental values of structural parameters, band gaps, magnetic structures, and oxygen vacancy formation energies for six representative oxides, i.e., CaMnO 3 , SrMnO 3 , LaMnO 3 , YMnO 3 , LaFeO 3 , and CeO 2 . Further, our results show that SCAN usually has better agreement with the experimental lattice constants and band gaps and larger magnetic moments than the commonly used GGA functionals. Although SCAN overestimates the oxygen vacancy formation energies of transition metal oxides and therefore requires unusually large Hubbard U values to reproduce the experimental reduction enthalpies, it does predict the correct oxygen vacancy formation energy of CeO 2 , which challenges the commonly used GGA and hybrid functionals. Our results underscore the challenges that exist in describing these complex oxides by DFT and may shed light on developing more accurate exchange-correlation functionals for oxygen vacancy formation energy calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Equation of State of CO 2 Shock Compressed to 1 TPa

Equation-of-state (pressure, density, temperature, internal energy) and reflectivity measurements on shock-compressed CO 2 at and above the insulating-to-conducting transition reveal new insight into the chemistry of simple molecular systems in the warm-dense-matter regime. CO 2 samples were precompressed in diamond-anvil cells to tune the initial densities from 1.35 g/cm 3 (liquid) to 1.74 g/cm 3 (solid) at room temperature and were then shock compressed up to 1 TPa and 93 000 K. Lastly, variation in initial density was leveraged to infer thermodynamic derivatives including specific heat and Gruneisen coefficient, exposing a complex bonded and moderately ionized state at the most extreme conditions studied.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Editorial: Particle interaction with afterglow plasma and non-quasi-neutral plasma

When immersed in a plasma environment, small nano-to micrometer sized solid particles undergo electrical charging by collecting charged species from the plasma. Once charged, these particles interact in peculiar ways with the surrounding plasma, triggering a widespread of interesting physical phenomena. Such systems, also known as complex or dusty plasmas, have been investigated for a few decades already.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Unnamed Pt(Cu 0.67 Sn 0.33 ) from the Bolshoy Khailyk River, Western Sayans, Russia, and a Review of Related Compounds and Solid Solutions

We describe a potentially new species of a platinum cupride–stannide mineral (PCSM) of composition Pt(Cu 0.67 Sn 0.33 ). It occurs in a placer deposit in the River Bolshoy Khailyk, southern Krasnoyarskiy kray, Russia. A synthetic equivalent of PCSM was obtained and characterized. The PCSM occurs as anhedral or subhedral grains up to 15 μm × 30 μm in association with various platinum-group minerals, Rh–Co-rich pentlandite and magnetite, all hosted by a placer grain of Cu–Au–Pt alloy. Synchrotron micro-Laue diffraction studies indicate that the PCSM mineral is tetragonal and belongs to the inferred space-group P4/mmm (#123). Its unit-cell parameters are a = 2.838 (3) Å, c = 3.650 (4) Å, and V = 29.40 (10) Å 3 , and Z = 1. The c:a ratio calculated from the unit-cell parameters is 1.286. These characteristics are in good agreement with those obtained for specimens of synthetic Pt(Cu 0.67 Sn 0.33 ). A review on related minerals and unnamed phases is provided to outline compositional variations and extents of solid solutions in the relevant systems PtNi–PtFe–PtCu, PdCu–PdHg–PdAu, PdHg–PtHg, and AuCu–PtCu. The PCSM-bearing mineralization appears to be related genetically with an ophiolitic source-rock of the Aktovrakskiy complex of the western Sayans. The unnamed phase likely crystallized from microvolumes of a highly fractionated melt rich in Cu and Sn.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Poison or Promoter? Investigating the Dual-Role of Carbon Monoxide in Pincer-Iridium-Based Alkane Dehydrogenation Systems via Operando Diffuse Reflectance Infrared Fourier Transform Spectroscopy

Pincer-ligated iridium complexes of the form [( R 4 PCP)IrL] ( R 4 PCP = κ 3 -C 6 H 3 -2,6-(XPR 2 ) 2 ; X = CH 2 , O; R = tBu, iPr) have previously been shown competent for acceptorless alkane dehydrogenation when supported on silica. It was observed by postcatalysis solid-state NMR that silica-tethered [(≡SiO- tBu 4 POCOP)Ir(C 2 H 4 )] (3-C 2 H 4 ) was converted fully to [(≡SiO- tBu 4 POCOP)Ir(CO)] (3-CO) at 300 °C. In this work, the characterization of species under dehydrogenation reaction conditions far from equilibrium between butane and butenes (approach to equilibrium Q/K eq = 0.3 at 300 °C) is performed with operando Diffuse Reflectance Infrared Fourier–Transform Spectroscopy (DRIFTS) to show the kinetics of species conversion from 3-C 2 H 4 to 3-CO. It is further found that [(≡SiO- tBu 4 POCOP)IrHCl] (3-HCl), a species considered to be a precatalyst for alkane dehydrogenation, is also fully converted to 3-CO. Here, a mechanism of decomposition is proposed that implicates surface silanol groups, while carbon monoxide acts as a “stabilizer” for the catalyst by promoting their reductive elimination and maintaining the complex in the I oxidation state.

02 PETROLEUM↗
DOE OSTI2025

Local chemical ordering of a neutron-irradiated CrFeMnNi compositionally complex alloy

While ion-irradiation studies are a critical first step in studying compositionally complex alloys (CCAs) for nuclear applications, they do not capture all the microstructural changes occurring under the low irradiation dose rates and different particles’ scattering patterns experienced in a nuclear reactor setting. To explore these phenomena in reactor-relevant conditions for the first time in CCA, the single-phase solid-solution Cr 10 Fe 30 Mn 30 Ni 30 was neutron irradiated up to 6.61 displacements per atom at 395 and 579 °C. Irradiation-enhanced local chemical ordering (LCO) well beyond the range of short range ordering was observed, and is predicted to be the precursor to the precipitation of a coherent Ni-Mn L1 0 phase and a Cr-rich α’ phase, though TEM analysis did not indicate the presence of either in any irradiation condition. The line density of faulted dislocation loops decreased from 6.47 to 1.69 ∙ 10 15 m -2 from 3.43 to 6.61 dpa at 579 °C despite no appreciable faulted loop content in the unirradiated material. LCO is expected to increase the complexity of the energy landscape within this alloy, restricting interstitial point defect mobility and creating local regions of greater stacking fault energy. These contribute to the negative correlation between irradiation dose and faulted dislocation loop density in this study, as well as the lack of void swelling observed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

In situ Raman investigation of Dy complexation in Cl-bearing aqueous solutions at 20–300 °C

Raman spectroscopy provides a versatile tool for in situ characterization of aqueous rare earth elements (REE) speciation at the molecular level. Complexation of REE with ligands such as Cl – and OH – is of particular interest for understanding the mobility of REE in NaCl-bearing hydrothermal fluids responsible for enriching REE to economic levels in nature. Raman spectroscopic studies of REE speciation in Cl-bearing aqueous fluids are primarily conducted at ambient temperature, whereas natural systems indicate temperatures of >100–600 °C. In this study, the speciation of Dy in acidic chloride-bearing hydrothermal solutions was investigated using confocal Raman spectroscopy with a new capillary Raman heating stage at 20–300 °C. Background solutions (pure water, NaCl-solutions) and solutions with 0.14–1.8 mol kg –1 dissolved DyCl 3 were sealed in quartz capillary cells. Comparison of the spectra for Dy chloride solutions with those for background solutions and the spectra for reference Dy-bearing solids was used to identify Raman bands specific to Dy–O and Dy–Cl bonds. The Raman band for the Dy–O stretching mode of hydrated Dy 3+ aqua ions was measured at 365–384 cm –1 and a Raman band for the Dy–Cl stretching modes of Dy chloride complexes was measured near 240 cm –1 . The Dy–O band decreases systematically with temperature, whereas the Dy–Cl band systematically increases, indicating a systematic increase in the stability of Dy chloride complexes with temperature. Here, this study provides the framework for expanding the use of in situ Raman spectroscopy to investigate the speciation of REE in aqueous solutions to hydrothermal conditions.

58 GEOSCIENCES↗
DOE OSTI2025

Spin‐Orbit Effects in a Thallium Borohydride Stabilized by Coordination to Bis(diisopropylamino)Cyclopropenylidene (BAC)

The reaction of Tl(OTf) with 2 equiv of bis(diisopropylamino)cyclopropenylidene (BAC) in THF results in formation of [Tl(BAC) 2 (OTf)] (1) in moderate yields. Subsequent reaction of 1 with [K][H 2 -9-BBN] ([H 2 -9-BBN] − = dihydrido 9-boratabicyclo[3.3.1]nonane) in THF results in formation of [Tl(BAC)(μ-H 2 -9-BBN)] 2 (3), also in moderate yield. Complex 3 is the first reported thallium borohydride. We attribute its thermal stability to the strong donor ability of the BAC co-ligand. Both 1 and 3 exhibit trigonal pyramidal geometries about Tl + in the solid-state, indicative of the presence of stereochemically active lone pairs. The hydride environment in 3 is calculated to exhibit a 3.9 ppm downfield shift attributed to spin-orbit effects from the adjacent Tl center.

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

Exploring pressure-dependent inelastic deformation and failure in bonded granular composites: An energetic materials perspective

In polymer-filled granular composites, damage may develop in mechanical loading prior to material failure. Damage mechanisms such as microcracking or plastic deformation in the binder phase can substantially alter the material’s mesostructure. For energetic materials, such as solid propellants and plastic bonded explosives, these mesostructural changes can have far reaching effects including degraded mechanical properties, potentially increased sensitivity to further insults, and changes in expected performance. Unfortunately, predicting damage is nontrivial due to the complex nature of these composites and the entangled interactions between inelastic mechanisms. In this work, we assess the current literature of experimental knowledge, focusing on the pressure-dependent shear response, and propose a simple simulation framework of bonded particles to study four limiting-case material formulations at both meso- and macro-scales. Further, to construct the four cases, we systematically vary the relative interfacial strength between the polymer binder and granular filler phase and also vary the polymer’s glass transition temperature relative to operating temperature which determines how much the binder can plastically deform. These simulations identify key trends in global mechanical response, such as the emergence of strain hardening or softening regimes with increasing pressure which qualitatively resemble experimental results. By quantifying the activation of different inelastic mechanisms, such as bonds breaking and plastically straining, we identify when each mechanism becomes relevant and provide insight into potential origins for changes in mechanical responses. The locations of broken bonds are also used to define larger, mesoscopic cracks to test various metrics of damage. We primarily focus on triaxial compression, but also test the opposite case of triaxial extension to highlight the impact of Lode angle on mechanical behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Mesoscale Modeling of Hydrogels Under Frictional Shear Stress

Hydrogels are three-dimensional networks of hydrophilic polymers often used as a simplified model of hydrated biological materials, from cartilaginous joints to the ocular tear film. However, the lubrication mechanisms of hydrogels remain poorly understood, partly due to their complex polymeric structure, which creates blurred interfaces during sliding that are challenging to study experimentally. In this study, we employ dissipative particle dynamics (DPD) to investigate the frictional behavior of a polymeric hydrogel network sliding against a solid wall in an explicit viscous solvent. This computational approach enables us to model hydrodynamic interactions and mesoscale polymer dynamics, capturing key aspects of hydrogel friction. Our simulations reveal that hydrogel friction is governed by the interplay between polymer relaxation and viscous shear, characterized by the Weissenberg number (Wi). At low Wi, friction coefficient remain nearly constant, dominated by polymer relaxation. However, at higher Wi, friction is dominated by viscous drag within a near-wall solvent layer, leading to a linear increase in friction coefficient with Wi. Furthermore, our results demonstrate an inverse relationship between the friction coefficient and the applied normal load, consistent with experimental observations. This work provides new insights into the fundamental tribological properties of hydrogels, shedding light on the micromechanics of hydrogel friction. Improving our understanding of hydrogel structure and dynamics under friction advances our knowledge of the mechanisms regulating biological lubrication in health and disease.

36 MATERIALS SCIENCE↗
DOE OSTI2025

ThO 2 and Th 1– x U x O 2 Nanoscale Materials and Thin Films for Nuclear Science Applications

This study investigates the dynamics and mechanisms of solution combustion synthesis (SCS) for the preparation of nanoscale ThO 2 and Th 1–x U x O 2 materials, utilizing metal nitrates (Th(NO 3 ) 4 and UO 2 (NO 3 ) 2 ) and acetylacetone (C 5 H 8 O 2 ) as reactants dissolved in a 2-methoxyethanol (C 3 H 8 O 2 ) solvent. By combining thermodynamic calculations, dynamic time–temperature profile measurements with differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA), this research reveals how variations in acetylacetone concentration and uranium content influence the structural parameters of the synthesized oxides. The time–temperature measurements show that the heating rate and maximum combustion temperatures are sensitive to acetylacetone concentration. DSC-TGA results indicate shifts in exothermic peak temperatures as the uranium content changes. The complexation between thorium and acetylacetone emerges as a critical factor, impacting combustion parameters and the structural characteristics of the final products. The uniform distribution of Th and U in the Th 1–x U x O 2 solid solution and the formation of nanoscale particles with strained crystallites are considered essential for the low-temperature densification of these materials for nuclear fuel pellet applications. Additionally, high-quality ThO 2 and Th 1–x U x O 2 thin (100–150 nm) films are successfully synthesized via electrospray deposition of combustible solutions followed by a brief period of heat treatment. Furthermore, these films exhibit excellent structural and morphological uniformity, making them ideal candidates for nuclear measurements, irradiation damage studies, and investigations into the physical properties of both pure and mixed oxides.

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

Impact of surface hydrophilicity on the ordering and transport properties of bicontinuous microemulsions

Microemulsions (MEs) have many industrial applications, where recent developments have shown that MEs can be utilized for electrochemical applications, including potentially in redox flow batteries. However, understanding the structure and dynamics of these systems, including at a surface, is needed to direct and rationally control their electrochemical behavior. While bulk solution measurements have provided insight into their structure, their assembly at an interface also impacts the electron (to the electrode) and ion (across the surfactant) charge transfer processes in the system. To address this shortcoming, neutron reflectivity experiments and molecular simulations have been completed that document the near surface structure of a series of deuterated water (D 2 O)/polysorbate-20/toluene MEs on hydrophilic and amphiphilic surfaces. These results show that the microemulsions form complex layered structures near a hard electrode surface, where most layers are not purely one component. Decreasing the D 2 O concentration in the ME increases the number of and purity of the layers established on the solid surface. These lamellar-type layers transition from the surface to the bulk microemulsion as a series of mixed layers (i.e., containing oil, water, and surfactant) that are consistent with perforated lamellae. Additionally, these mixed lamellae appear to become more perforated with oil and water pathways on an amphiphilic surface. Furthermore, the purity and thickness of these layers will influence the accessibility of an electrode by redox active species, as well as ion transport required to satisfy the electroneutrality condition.

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

Simulation of Reverse Osmosis Membrane Compaction Using Material Point Method

Reverse Osmosis (RO) is a promising technology to address the impending water-crisis in the upcoming decades. RO at high salinities and high pressures is challenging due to membrane compaction that changes its porosity and permeability. In this study, we present a simulation methodology for membrane structural mechanics to understand pore size distribution and permeability variations under high pressure. We use the material point method (MPM), that solves the solid mechanics equations in a Lagrangian framework. MPM provides many features that make it well-suited for simulating mesoporous membranes. The Lagrangian framework allows for large deformations, easy integration of constitutive models and direct import of complex geometries as particles. The spatial discretization in our MPM solver is achieved using linear or cubic-spline shape functions while the time integration is carried out using the explicit Euler method. A series of images containing detailed pore structures obtained from X-ray tomography is first converted to a collection of material points. Compressive loads are applied to the top layer of the membrane to simulate the application of pressure. The membrane deformation and pore size distribution before and after load application are reported and compared with the experimental measurements. The presentation discusses the numerical methods used, the performance of the solver on high-performance computing machines, and the results of membrane compact in detail.

compaction↗
DOE OSTI2023

Prediction of alkaline earth metal ion adsorption on goethite for various background electrolytes with the CD-MUSIC model

As water scarcity drives the use of more saline water sources, contaminant fate and transport models must capture the impact of high concentrations of alkaline earth metal ions (AEMs) and background electrolytes in these more complex waters. By utilizing macroscopic adsorption data from various electrolyte systems, a Charge Distribution – Multisite Complexation (CD-MUSIC) model, capable of incorporating electrolyte adsorption, was able to accurately simulate the adsorption behavior of alkaline earth metal ions onto goethite. The modeling effort was guided by previous spectroscopic and surface complexation modeling of alkaline earth metal adsorption and built on previous CD-MUSIC modeling that accounted for changes in crystal face contributions to the surface site density as a function of specific surface area. Further, the model was constrained to consider only two dominant surface complex species for each metal ion adsorption reaction. These two species were selected from 44 possible species through objective curve fitting of single-solute macroscopic adsorption data. While most of the alkaline earth metal surface complexes formed outer-sphere complexes at the goethite surface, an inner-sphere species was utilized for Mg 2+ . With the surface complex species and equilibrium constants obtained from this study, the calibrated model successfully predicted alkaline earth metal ion adsorption over a wide range of solution and surface conditions; the model predictions encompassed a wide range of pH (5–11), solute/solid ratio (1.37 × 10 -5 – 8.33 × 10 -4 mol -solute /g -solid ), ionic strengths (0.01 M – 0.7 M), and background electrolytes (Na + , Cs + , Rb + , Cl - , and NO 3 - ) using the same crystal face contribution methodology for site density, capacitance values, and surface acidity constants adopted for proton and cadmium adsorption in previous work (Han and Katz, 2019). Model simulations for a range of background water chemistries demonstrated the potential for Mg 2+ to reduce Cd 2+ adsorption to goethite in model seawater and oil- and gas-produced waters.

42 ENGINEERING↗