Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SOLUTES”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Dirac solutions for quaternionic potentials

The Dirac equation is solved for quaternionic potentials, i V{sub 0} + j W{sub 0} (V{sub 0}∈R , W{sub 0}∈C). The study shows two different solutions. The first one contains particle and anti-particle solutions and leads to the diffusion, tunneling, and Klein energy zones. The standard solution is recovered taking the complex limit of this solution. The second solution, which does not have a complex counterpart, can be seen as a V{sub 0}-antiparticle or |W{sub 0}|-particle solution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Impact of solutes on the lattice parameters and elastic stiffness coefficients of hcp Fe from first-principles calculations

The hexagonal close-packed (hcp) -martensite phase in steels nucleates from the -austenite parent phase and can undergo further transformation to the -martensite phase or exist as a metastable phase depending on temperature, mechanical loading, and alloy chemistry. The solute-dependent lattice parameters and elastic stiffness coefficients of hcp Fe influence the mechanical properties of steels containing the -martensite phase, as well as the martensitic transformations between the phases. We use density functional theory to calculate the lattice parameters and of single-crystal hcp Fe as functions of solute concentration in the dilute limit for the substitutional solutes Al, B, Cu, Mn, and Si, and the octahedral interstitial solutes C and N. Our computationally efficient methodology separates the solute dependence of the into lattice strain and chemical bonding contributions. The computed data can be used to estimate the effect of solutes on polycrystalline elastic moduli and the strain energy associated with martensitic transformations. The data can also serve as inputs to microstructure-based models of multiphase steels containing the -martensite phase.

Fellinger, Michael R.↗
DOE OSTI2016

Strengthening of the Coordination Shell by Counter Ions in Aqueous Th 4+ Solutions

The presence of counter ions in solutions containing highly charged metal cations can trigger processes such as ion-pair formation, hydrogen bond breakages and subsequent reformation, and ligand exchanges. In this work, it is shown how halide (Cl-, Br-) and perchlorate (ClO4-) anions affect the strength of the primary solvent coordination shells around Th4+ using explicit solvent and finite temperature ab initio molecular dynamics modeling methods. The 9-fold solvent geometry was found to be the most stable hydration structure in each aqueous solution. Relative to the dilute aqueous solution, the presence of the counter ions did not significantly alter the geometry of the primary hydration shell. However, the free energy analyses indicated that the 10-fold hydrated states were thermodynamically accessible in dilute and bromide aqueous solutions within 1 kcal/mol. Analysis of the results showed that the hydrogen bond lifetimes were longer and solvent exchange energy barriers were larger in solutions with counter ions in comparison with the solution with no counter ions. This implies that the presence of the counter ions induces a strengthening of the Th4+ hydration shell.

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

A Perturbative Solution for Nonlinear Stratified Upwelling over a Frictional Slope

Abstract A perturbative solution of simplified primitive equations for nonlinear weakly stratified upwelling over a frictional slope is found that resolves the vertical structure of velocity fields and can satisfy Ertel’s potential vorticity conservation in the stratified inviscid interior. The solution uses assumptions consistent with the model proposed by Lentz and Chapman, including a steady-state, constant cross-shore density gradient, no alongshore gradients, laterally inviscid, and consideration of cross-shore advection of alongshore momentum. The solution resolves the vertical structure of velocity fields (including subsurface maxima of compensational flow, not resolved by Lentz and Chapman) and can satisfy Ertel’s potential vorticity conservation in the stratified inviscid interior. The dynamics are similar to Lentz and Chapman; bottom stress balances alongshore wind stress in a homogeneous density ocean and is replaced by nonlinear cross-shore transport of alongshore momentum as the Burger number (S=αN/f, whereα,N, andfare the bottom slope, buoyancy frequency, Coriolis frequency, respectively) increases. When the solution uses the empirical relation between cross-shore and vertical density gradients proposed by Lentz and Chapman, vorticity conservation is not satisfied and the nonlinear momentum transport estimated by the solution linearly increases withS, asymptotically matching Lentz and Chapman forS< 1. When the solution conserves interior potential vorticity, the momentum transport is proportional toS 2 forS< 1 and is in better agreement with numerical simulations.

Oceanography↗
DOE OSTI2016

Intercomparison of 3D pore-scale flow and solute transport simulation methods

Multiple numerical approaches have been developed to simulate porous media fluid flow and solute transport at the pore scale. These include methods that 1) explicitly model the three-dimensional geometry of pore spaces and 2) those that conceptualize the pore space as a topologically consistent set of stylized pore bodies and pore throats. In previous work we validated a model of class 1, based on direct numerical simulation using computational fluid dynamics (CFD) codes, against magnetic resonance velocimetry (MRV) measurements of pore-scale velocities. Here we expand that validation to include additional models of class 1 based on the immersed-boundary method (IMB), lattice Boltzmann method (LBM), smoothed particle hydrodynamics (SPH), as well as a model of class 2 (a pore-network model or PNM). The PNM approach used in the current study was recently improved and demonstrated to accurately simulate solute transport in a two-dimensional experiment. While the PNM approach is computationally much less demanding than direct numerical simulation methods, the effect of conceptualizing complex three-dimensional pore geometries on solute transport in the manner of PNMs has not been fully determined. We apply all four approaches (CFD, LBM, SPH and PNM) to simulate pore-scale velocity distributions and nonreactive solute transport, and intercompare the model results with previously reported experimental observations. Experimental observations are limited to measured pore-scale velocities, so solute transport comparisons are made only among the various models. Comparisons are drawn both in terms of macroscopic variables (e.g., permeability, solute breakthrough curves) and microscopic variables (e.g., local velocities and concentrations).

computational fluid dynamics↗
DOE OSTI2012

Structure and dimerization of translation initiation factor aIF5B in solution

Translation initiation factor 5B (IF5B) is required for initiation of protein synthesis. The solution structure of archaeal IF5B (aIF5B) was analysed by small-angle X-ray scattering (SAXS) and dynamic light scattering (DLS) and was indicated to be in both monomeric and dimeric form. Sedimentation equilibrium (SE) analytical ultracentrifugation (AUC) of aIF5B indicated that aIF5B forms irreversible dimers in solution but only to a maximum of 5.0-6.8% dimer. Sedimentation velocity (SV) AUC at higher speed also indicated the presence of two species, and the sedimentation coefficients s{sub 20,w}{sup 0} were determined to be 3.64 and 5.51 {+-} 0.29 S for monomer and dimer, respectively. The atomic resolution (crystallographic) structure of aIF5B (Roll-Mecak et al. [6]) was used to model monomer and dimer, and theoretical sedimentation coefficients for these models were computed (3.89 and 5.63 S, respectively) in good agreement with the sedimentation coefficients obtained from SV analysis. Thus, the structure of aIF5B in solution must be very similar to the atomic resolution structure of aIF5B. SAXS data were acquired in the same buffer with the addition of 2% glycerol to inhibit dimerization, and the resultant monomeric aIF5B in solution did indeed adopt a structure very similar to the one reported earlier for the protein in crystalline form. The p(r) function indicated an elongated conformation supported by a radius of gyration of 37.5 {+-} 0.2 {angstrom} and a maximum dimension of {approx}130 {angstrom}. The effects of glycerol on the formation of dimers are discussed. This new model of aIF5B in solution shows that there are universal structural differences between aIF5B and the homologous protein IF2 from Escherichia coli.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

The importance of Ion Interactions on Electrolyte Solution Viscosities Determined by Comparing Concentrated Sodium Carbonate and Nitrate Solutions

Sodium carbonate is a common electrolyte in high pH aqueous solution, and its contribution to the viscosity of those electrolytes solutions is of practical importance. Nitrite is isostructural with carbonate, and has nearly identical ionic radius. The relative importance of ionic size versus interactions with water on electrolyte solution viscosity at moderate and high electrolyte concentrations has been ambiguous in the literature for a long time. Comparing the viscosity of sodium carbonate and sodium nitrate solution showed that carbonate solutions are considerably more viscous than nitrate solutions at constant sodium molality even though there is half as many anions in the carbonate solution as nitrate solution. This clearly shows that interactions with water are more important than ion size for the viscosity of these solutions because these two electrolytes with nearly identical ion size have very large differences in viscosity. This comparison was aided by new data on the viscosity of sodium carbonate solutions, which was also used to parametrize the Laliberté model of electrolyte solution viscosities between 20 and 98 °C.

Viscosity, Hanford, Sodium Carbonate↗
DOE OSTI2020

Potential-pH diagrams considering complex oxide solution phases for understanding aqueous corrosion of multi-principal element alloys

The potential-pH diagram, a graphical representation of the thermodynamically predominant reaction products in aqueous corrosion, is originally proposed for the corrosion of pure metals. The original approach only leads to stoichiometric oxides and hydroxides as the oxidation products. However, numerous experiments show that non-stoichiometric oxide scales are prevalent in the aqueous corrosion of alloys. In the present study, a room temperature potential-pH diagram considering oxide solid solutions, as a generalization of the traditional potential-pH diagram with stoichiometric oxides, is constructed for an FCC single-phase multi-principal element alloy (MPEA) based on the CALculation of PHAse Diagram method. The predominant reaction products, the ions in aqueous solution, and the cation distribution in oxides are predicted. The oxide solid solution is stabilized by the mixing free energy (or mixing entropy) and the stabilizing effect becomes more significant as the temperature increases. Consequently, solid solution oxides are stable in large regions of the potential-pH diagram and the mixing free energy mostly affects the equilibrium composition of the stable oxides, while the shape of stable regions for oxides is mostly determined by the structure of the stable oxides. Agreements are found for Ni 2+ , Fe 2+ , and Mn 2+ between the atomic emission spectroelectrochemistry measurements and thermodynamic calculations, while deviations exist for Cr 3+ and Co 2+ possibly due to surface complexation with species such as Cl - and the oxide dissolution. By incorporating the solution models of oxides, the current work presents a general and more accurate way to analyze the reaction products during aqueous corrosion of MPEAs.

Materials Science↗
DOE OSTI2017

Solution to the sign problem in a frustrated quantum impurity model

In this work we solve the sign problem of a frustrated quantum impurity model consisting of three quantum spin-half chains interacting through an anti-ferromagnetic Heisenberg interaction at one end. We first map the model into a repulsive Hubbard model of spin-half fermions hopping on three independent one dimensional chains that interact through a triangular hopping at one end. We then convert the fermion model into an inhomogeneous one dimensional model and express the partition function as a weighted sum over fermion worldline configurations. By imposing a pairing of fermion worldlines in half the space we show that all negative weight configurations can be eliminated. This pairing naturally leads to the original frustrated quantum spin model at half filling and thus solves its sign problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Nonadiabatic rate constants for proton transfer and proton-coupled electron transfer reactions in solution: Effects of quadratic term in the vibronic coupling expansion

Rate constant expressions for vibronically nonadiabatic proton transfer and proton-coupled electron transfer reactions are presented and analyzed. The regimes covered include electronically adiabatic and nonadiabatic reactions, as well as high-frequency and low-frequency proton donor-acceptor vibrational modes. These rate constants differ from previous rate constants derived with the cumulant expansion approach in that the logarithmic expansion of the vibronic coupling in terms of the proton donor-acceptor distance includes a quadratic as well as a linear term. The analysis illustrates that inclusion of this quadratic term in the framework of the cumulant expansion framework may significantly impact the rate constants at high temperatures for proton transfer interfaces with soft proton donor-acceptor modes that are associated with small force constants and weak hydrogen bonds. The effects of the quadratic term may also become significant in these regimes when using the vibronic coupling expansion in conjunction with a thermal averaging procedure for calculating the rate constant. In this case, however, the expansion of the coupling can be avoided entirely by calculating the couplings explicitly for the range of proton donor-acceptor distances sampled. The effects of the quadratic term for weak hydrogen-bonding systems are less significant for more physically realistic models that prevent the sampling of unphysical short proton donor-acceptor distances. Additionally, the rigorous relation between the cumulant expansion and thermal averaging approaches is clarified. In particular, the cumulant expansion rate constant includes effects from dynamical interference between the proton donor-acceptor and solvent motions and becomes equivalent to the thermally averaged rate constant when these dynamical effects are neglected. This analysis identifies the regimes in which each rate constant expression is valid and thus will be important for future applications to proton transfer and proton-coupled electron transfer in chemical and biological processes.

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

Countercations Control Local Specific Bonding Interactions and Nucleation Mechanisms in Concentrated Water-in-Salt Solutions

One of the continuing challenges presented in aqueous salt solutions is to understand ion association reactions driving dynamic demixing from solvation, complexation and solute clustering. The problems understanding this phenomenon are exacerbated in the highly-concentrated water-in-salt electrolytes, where the deficiency of water leads to a dramatic retardation of solvating water and formation of extended solvent-solute clustering networks. By probing microscopic dynamics of water and prenucleation clusters using quasi-elastic neutron scattering and proton nuclear magnetic resonance spectroscopy, we observed contrasting mechanistic specifics of ion-water mobilities in highly-concentrated Na+- vs. K+-based aluminate solutions (diffusion coefficients of 0.2 vs. 2.6?10-10m2s-1 at 293 K, respectively). The magnitude of the differences is far beyond the kosmotropic (Na+) and chaotropic (K+) effects typically understood at low ionic concentrations. The distinct crystallization mechanisms observed in analogous deuterated systems further imply that different prenucleation cluster dynamics and network packing morphologies can either frustrate or promote crystallization, as described by non-classical nucleation theory.

Wang, Hsiu-Wen↗
DOE OSTI2015

Size-dependent error of the density functional theory ionization potential in vacuum and solution

Density functional theory is often the method of choice for modeling the energetics of large molecules and including explicit solvation effects. It is preferable to use a method that treats systems of different sizes and with different amounts of explicit solvent on equal footing. However, recent work suggests that approximate density functional theory has a size-dependent error in the computation of the ionization potential. We here investigate the lack of size-intensivity of the ionization potential computed with approximate density functionals in vacuum and solution. We show that local and semi-local approximations to exchange do not yield a constant ionization potential for an increasing number of identical isolated molecules in vacuum. Instead, as the number of molecules increases, the total energy required to ionize the system decreases. Rather surprisingly, we find that this is still the case in solution, whether using a polarizable continuum model or with explicit solvent that breaks the degeneracy of each solute, and we find that explicit solvent in the calculation can exacerbate the size-dependent delocalization error. We demonstrate that increasing the amount of exact exchange changes the character of the polarization of the solvent molecules; for small amounts of exact exchange the solvent molecules contribute a fraction of their electron density to the ionized electron, but for larger amounts of exact exchange they properly polarize in response to the cationic solute. In vacuum and explicit solvent, the ionization potential can be made size-intensive by optimally tuning a long-range corrected hybrid functional.

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

Local atomic environment analysis of short and long-range solute-solute interactions in a symmetric tilt grain boundary

We develop an efficient molecular-statics algorithm for predicting local solute-solute interaction effects on grain-boundary segregation. This algorithm is applied to the segregation and interaction of Ni and Cu solutes in a Σ9 (221) symmetric tilt grain boundary in Ag. It is found that the strength of interaction between two solutes directly correlates with the change of average segregation energy in the local atomic environment of a fixed grain-boundary solute. It is shown that local interaction effects at grain boundaries can be either attractive or repulsive wherein long-range attractive effects increase the propensity for solute cluster nucleation. This study underscores the atomic nature of grain-boundary solute segregation and provides a tool that could yield meaningful insight into heterogeneous solute segregation phenomena in polycrystalline alloys beyond their dilute limit.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Ion aggregation in high salt solutions. III. Computational vibrational spectroscopy of HDO in aqueous salt solutions

The vibrational frequency, frequency fluctuation dynamics, and transition dipole moment of the O—D stretch mode of HDO molecule in aqueous solutions are strongly dependent on its local electrostatic environment and hydrogen-bond network structure. Therefore, the time-resolved vibrational spectroscopy the O—D stretch mode has been particularly used to investigate specific ion effects on water structure. Despite prolonged efforts to understand the interplay of O—D vibrational dynamics with local water hydrogen-bond network and ion aggregate structures in high salt solutions, still there exists a gap between theory and experiment due to a lack of quantitative model for accurately describing O—D stretch frequency in high salt solutions. To fill this gap, we have performed numerical simulations of Raman scattering and IR absorption spectra of the O—D stretch mode of HDO in highly concentrated NaCl and KSCN solutions and compared them with experimental results. Carrying out extensive quantum chemistry calculations on not only water clusters but also ion-water clusters, we first developed a distributed vibrational solvatochromic charge model for the O—D stretch mode in aqueous salt solutions. Furthermore, the non-Condon effect on the vibrational transition dipole moment of the O—D stretch mode was fully taken into consideration with the charge response kernel that is non-local polarizability density. From the fluctuating O—D stretch mode frequencies and transition dipole vectors obtained from the molecular dynamics simulations, the O—D stretch Raman scattering and IR absorption spectra of HDO in salt solutions could be calculated. The polarization effect on the transition dipole vector of the O—D stretch mode is shown to be important and the asymmetric line shapes of the O—D stretch Raman scattering and IR absorption spectra of HDO especially in highly concentrated NaCl and KSCN solutions are in quantitative agreement with experimental results. We anticipate that this computational approach will be of critical use in interpreting linear and nonlinear vibrational spectroscopies of HDO molecule that is considered as an excellent local probe for monitoring local electrostatic and hydrogen-bonding environment in not just salt but also other confined and crowded solutions.

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

Solution Irregularity Remediation for Spatial Discretization Error Estimation for S N Transport Solutions

The discrete ordinates linear Boltzmann transport equation is typically solved in its spatially discretized form, incurring spatial discretization error. Quantification of this error for purposes such as adaptive mesh refinement or error analysis requires an a posteriori estimator, which utilizes the numerical solution to the spatially discretized equation to compute an estimate. Because the quality of the numerical solution informs the error estimate, irregularities, present in the true solution for any realistic problem configuration, tend to cause the largest deviation in the error estimate vis-a-vis the true error. In this paper, an analytical partial singular characteristic tracking (pSCT) procedure for reducing the estimator’s error is implemented within our novel residual source estimator for a zeroth-order discontinuous Galerkin scheme, at the additional cost of a single inner iteration. Here, a metric-based evaluation of the pSCT scheme versus the standard residual source estimator is performed over the parameter range of a Method of Manufactured Solutions test suite. The pSCT scheme generates near-ideal accuracy in the estimate in problems where the dominant source of the estimator’s error is the solution irregularity, namely, problems where the true solution is discontinuous and problems where the true solution’s first derivative is discontinuous and the scattering ratio is low. In problems where the scattering ratio is high and the true solution is discontinuous in the first derivative, the error in the scattering source, which is not converged by the pSCT scheme, is greater than the error incurred due to the irregularity. Ultimately, a pSCT scheme is judged to be useful for error estimation in problems where the computational cost of the scheme is justified. In the presence of many irregularities, such a scheme may be intractable for general use, but in benchmarks, as an analytical tool, or in problems that have nondissipative discontinuities, the scheme may prove invaluable.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Electronic structure of aqueous solutions: Bridging the gap between theory and experiments

Predicting the electronic properties of aqueous liquids has been a long-standing challenge for quantum mechanical methods. However, it is a crucial step in understanding and predicting the key role played by aqueous solutions and electrolytes in a wide variety of emerging energy and environmental technologies, including battery and photoelectrochemical cell design. Here, we propose an efficient and accurate approach to predict the electronic properties of aqueous solutions, on the basis of the combination of first-principles methods and experimental validation using state-of-the-art spectroscopic measurements. We present results of the photoelectron spectra of a broad range of solvated ions, showing that first-principles molecular dynamics simulations and electronic structure calculations using dielectric hybrid functionals provide a quantitative description of the electronic properties of the solvent and solutes, including excitation energies. The proposed computational framework is general and applicable to other liquids, thereby offering great promise in understanding and engineering solutions and liquid electrolytes for a variety of important energy technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗