Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solutions”

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

Elucidating Actinide–Pertechnetate and Actinide–Perrhenate Bonding via a Family of Th–TcO 4 and Th–ReO 4 Frameworks and Solutions

Technetium-99, a β-emitter produced from 235 U fission, poses a challenge for the nuclear industry due to co-extraction of pertechnetate (TcO 4 – ) with the actinides (An) during nuclear fuel reprocessing. Previous studies suggested that direct coordination of pertechnetate with An plays an important role in the coextraction process. However, few studies have provided direct evidence for An–TcO 4 – bonding in the solid state, and even fewer in solution. The present study describes synthesis and structural elucidation of a family of thorium(IV)-pertechnetate/perrhenate (ReO 4 – , nonradioactive surrogate) compounds, which is obtained by dissolution of thorium oxyhydroxide in perrhenic/pertechnic acid followed by crystallization, with or without heating. For reaction ratios of 3:1, 4:1, and 6:1 MO 4 – /Th(IV) (M = Tc, Re), the crystallized compounds reflect the same ratio, suggesting facile and flexible coordination. Furthermore, nine structures reveal 1-dimensional and 2-dimensional frameworks with varying topologies. While a multitude of compounds isolated from 4:1 (and 6:1) reaction solutions feature Th monomers linked by MO 4 – , the 3:1 reaction solution yielded the well-known dihydroxide-bridged thorium dimer, linked, and capped by MO 4 – . Density functional theory calculations on ReO 4 – /TcO 4 – isomorphs suggest similar bonding characteristics in the solid state, but experimental solution characterization noted differences. Specifically, small-angle X-ray scattering studies suggest the bonding of Th–TcO 4 – persists in solution, while Th–ReO 4 – bonding is less apparent.

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

Dynamics of Aqueous Electrolyte Solutions: Challenges for Simulations

This Perspective article focuses on recent simulation work on the dynamics of aqueous electrolytes. It is well-established that full-charge, nonpolarizable models for water and ions generally predict solution dynamics that are too slow in comparison to experiments. Models with reduced (scaled) charges do better for solution diffusivities and viscosities but encounter issues describing other dynamic phenomena such as nucleation rates of crystals from solution. Polarizable models show promise, especially when appropriately parametrized, but may still miss important physical effects such as charge transfer. First-principles calculations are starting to emerge for these properties that are in principle able to capture polarization, charge transfer, and chemical transformations in solution. Finally, while direct ab initio simulations are still too slow for simulations of large systems over long time scales, machine-learning models trained on appropriate first-principles data show significant promise for accurate and transferable modeling of electrolyte solution dynamics.

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

Beyond Simple Dilution: Superior Conductivities from Cosolvation of Acetonitrile/LiTFSI Concentrated Solution with Acetone

Concentrated solutions of Li salts in acetonitrile are promising alternative electrolytes for the next generation of Li batteries as they may exhibit superior electrochemical properties. However, the reduced mobility of the chemical species is a barrier yet to be overcome, and for this, we explore the utilization of acetone as a cosolvent. Although acetone is a polar compound, we find that its addition to the LiTFSI/acetonitrile solution does not follow the trends expected for a simple dilution process. At a low concentration, acetone subtly shifts acetonitrile from the first to extended solvation sheaths of the ions. Still, most of the original structure of the solution is preserved, and mobile high-concentration clusters are formed in the solution. At higher concentrations, the cosolvation promotes cation–anion interactions but with a different nature from those in the original solution and still allows for a further increase in conductivity. Additionally, the non-coordinating fraction of acetonitrile acquires features resembling the pure solvent, which is a possible additional facilitating factor for ionic diffusion.

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

NMR Studies of Block Copolymer-Based Supramolecules in Solution

Hierarchical assemblies from block copolymer (BCP)-based supramolecules have shown immense potential as programmable materials owing to their versatility for incorporating functional molecules and provide access to arrays of hierarchical structures. However, there remains a knowledge gap on the formation of the supramolecule in solution. Here, we applied NMR techniques to investigate the solution-phase behavior of the most studied supramolecular systems, polystyrene-block-poly(4-vinylpyridine)(3-pentadecylphenol) (PS-b-P4VP(PDP) r ). The results show that the supramolecule likely adopts a coil-comb conformation, despite the small molecule's (PDP) rapid exchange between the bonded and free states. The exchange rate (>10 4 s -1 ) exceeds the NMR time scale at the frequency of interest. The supramolecules form under dilute conditions (~2 vol %) and are attributed to the enthalpic gain of the hydrogen bonding between the PDP and 4VP. As the solute concentration increases (>10 vol %), the supramolecule forms micelle-like aggregates with PDP accumulated within the comb-block's pervaded volume based on analysis of the apparent molecular weight, viscosity, and chain dynamics. Overall, this work sheds light on the long-standing question regarding the evolution of the constituents in the BCP-based supramolecule in solution and provides valuable guidance toward their solution-based processing and morphological control.

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

Solution epitaxy of polarization-gradient ferroelectric oxide films with colossal photovoltaic current

Solution growth of single-crystal ferroelectric oxide films has long been pursued for the low-cost development of high-performance electronic and optoelectronic devices. However, the established principles of vapor-phase epitaxy cannot be directly applied to solution epitaxy, as the interactions between the substrates and the grown materials in solution are quite different. Here, we report the successful epitaxy of single-domain ferroelectric oxide films on Nb-doped SrTiO 3 single-crystal substrates by solution reaction at a low temperature of ~200 o C. The epitaxy is mainly driven by an electronic polarization screening effect at the interface between the substrates and the as-grown ferroelectric oxide films, which is realized by the electrons from the doped substrates. Atomic-level characterization reveals a nontrivial polarization gradient throughout the films in a long range up to ~500 nm because of a possible structural transition from the monoclinic phase to the tetragonal phase. This polarization gradient generates an extremely high photovoltaic short-circuit current density of ~2.153 mA/cm 2 and open-circuit voltage of ~1.15 V under 375 nm light illumination with power intensity of 500 mW/cm 2 , corresponding to the highest photoresponsivity of ~4.306×10 -3 A/W among all known ferroelectrics. Our results establish a general low-temperature solution route to produce single-crystal gradient films of ferroelectric oxides and thus open the avenue for their broad applications in self-powered photo-detectors, photovoltaic and optoelectronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

One-shot learning for solution operators of partial differential equations

Learning and solving governing equations of a physical system, represented by partial differential equations (PDEs), from data is a central challenge in many areas of science and engineering. Traditional numerical methods can be computationally expensive for complex systems and require complete governing equations. Existing data-driven machine learning methods require large datasets to learn a surrogate solution operator, which could be impractical. Here, we propose a solution operator learning method that requires only one PDE solution, i.e., one-shot learning, along with suitable initial and boundary conditions. Leveraging the locality of derivatives, we define a local solution operator in small local domains, train it using a neural network, and use it to predict solutions of new input functions via mesh-based fixed-point iteration or meshfree neural-network based approaches. We test our method on various PDEs, complex geometries, and a practical spatial infection spread application, demonstrating its effectiveness and generalization capabilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Machine-learning informed prediction of high-entropy solid solution formation: Beyond the Hume-Rothery rules

The empirical rules for the prediction of solid solution formation proposed so far in the literature usually have very compromised predictability. Some rules with seemingly good predictability were, however, tested using small data sets. Based on an unprecedented large dataset containing 1252 multicomponent alloys, machine-learning methods showed that the formation of solid solutions can be very accurately predicted (93%). The machine-learning results help identify the most important features, such as molar volume, bulk modulus, and melting temperature. As such a new thermodynamics-based rule was developed to predict solid–solution alloys. The new rule is nonetheless slightly less accurate (73%) but has roots in the physical nature of the problem. The new rule is employed to predict solid solutions existing in the three blocks, each of which consists of 9 elements. The predictions encompass face-centered cubic (FCC), body-centered cubic (BCC), and hexagonal closest packed (HCP) structures in a high throughput manner. The validity of the prediction is further confirmed by CALculations of PHAse Diagram (CALPHAD) calculations with high consistency (94%). Since the new thermodynamics-based rule employs only elemental properties, applicability in screening for solid solution high-entropy alloys is straightforward and efficient.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Surface interfacial analysis of simulant high level nuclear waste glass dissolved in synthetic cement solutions

Abstract The corrosion mechanisms and kinetics of a Mg-rich alkali aluminoborosilicate glass simulating UK high-level waste (CaZn28) were investigated upon dissolution in synthetic cement solutions. Dissolution varied as a function the different pH and alkali/alkaline earth content of each cement solution. High resolution microscopy and spectroscopy techniques ascertained the nature of the interface between the glass and the cement solutions. TEM-EDS revealed alkali- and alkaline earth-rich silica gels, into which K, Ca and Mg were incorporated. TEM-SAED, combined with synchrotron micro-focus XRD, identified the ubiquitous precipitation of the Mg-aluminate layered double hydroxide phase, meixnerite (Mg 6 Al 2 (OH) 18 ·4H 2 O), in addition to goethite (FeOOH) and crystalline silica. The C-S-H phase, tobermorite (Ca 5 Si 6 O 16 (OH) 2 ·4H 2 O), was identified in the most Ca-rich solution only. These data give insight to the role of alkali/alkaline earth-rich solutions in the dissolution or radioactive waste glasses, of importance to the final disposition in a geological disposal facility.

Corkhill, Claire L. (ORCID:0000000274883219)↗
DOE OSTI2026

Aqueous electrolyte solutions with anion-bridged secondary solvation sheaths for highly efficient zinc metal batteries

Aqueous zinc metal batteries are low-cost electrochemical devices suitable for safe grid energy storage. However, water decomposition and Zn dendrite formation detrimentally affect their coulombic efficiency. Conventional aqueous electrolyte solutions, with a concentration around 1 M, are cost-effective and exhibit high bulk ionic conductivity but cannot form a stable solid electrolyte interphase. Water-in-salt and aqueous-organic hybrid electrolyte solutions can form robust solid electrolyte interphases, but they are not kinetically efficient and cost-effective. Here, to circumvent these issues, we design variously concentrated aqueous electrolyte solutions using several salts with different donor numbers to extend anion coordination into the secondary solvation sheath. We show that salt-derived anions with donor number > 18 enter the Zn2+ first solvation sheath, and ensure a strong binding energy between the Zn2+(H2O)5-anion nanometric clusters and water molecules in the secondary solvation sheath. In particular, 2 M aqueous electrolyte solutions containing fluorinated anions exhibit bulk ionic conductivities of 26-35 mS cm−1 at 25 °C and form a ZnF2-rich solid electrolyte interphase. When tested in Zn||NaV3O8·1.5H2O Swagelok cells, the best-performing electrolyte solution enables an average coulombic efficiency of 99.99% for 1,000 cycles at 1.5 mA cm−2, corresponding to an initial specific energy of 130 Wh kg−1 (based on the combined weight of the positive and negative electrodes).

25 ENERGY STORAGE↗
DOE OSTI2022

Revealing temperature-dependent polymer aggregation in solution with small-angle X-ray scattering

We report that improving the morphology of bulk heterojunction active layers remains a primary challenge for organic photovoltaics (OPVs), and much research has been devoted to achieving this through modifying OPV casting solutions to control film formation and crystallinity. Yet, the solution conformation of conjugated polymers used in OPVs is largely unknown. Here, we report observations of temperature dependent aggregation (TDA) through small-angle X-ray scattering (SAXS) investigations of polymer conformation in chlorobenzene:dichlorobenzene casting solvent as a function of temperature for PffBT4T-2OD, a polymer known to display TDA, and its derivative PffBT3T-2OD which displays significantly reduced TDA. We find that, upon cooling below 80 °C, PffBT4T-2OD forms large crystalline aggregates in solution, while its derivative PffBT3T-2OD forms mostly amorphous aggregates of similar size with some evidence of short-range order. This change in solution aggregation behavior is reflected in the lack of gelation by PffBT3T-2OD upon film deposition by spin coating. Grazing-incidence wide-angle X-ray scattering (GIWAXS) revealed a preferred face-on π–π stacking orientation for PffBT3T-2OD films while PffBT4T-2OD's π–π stacking peak was isotropic. We combine these findings to suggest that the presence of crystalline seed aggregates in PffBT4T-2OD solution quickly form an isotropic crystallite network upon cooling while PffBT3T-2OD's amorphous aggregates more slowly crystallize resulting in improved processability of PffBT3T-2OD.

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

Non-thermal atmospheric pressure plasma–liquid synthesis of organic acids in aqueous solution from carbon monoxide

This work aims at understanding the conversion of CO to organic acids, namely oxalic acid and formic acid, using non-thermal atmospheric pressure plasma over aqueous solutions. CO exhibited significantly higher conversion to organic acids (more than 15×) compared to CO 2 under the same reaction conditions. The result bolsters a proposed two-step process for CO 2 fixation, whereby CO 2 is first converted to CO, and then CO is converted to organic acids. The organic acids produced from CO are intermediates in the water–gas shift (WGS) reaction of CO in the presence of an aqueous solution to dissolved CO 2 and hydrogen gas. Based on a simple thermodynamic analysis, the organic acid yield was increased by lowering the plasma–liquid reaction temperature using an ice bath to cool the reaction flask. The composition of the organic acids could be varied by changing the pH of the solution. Oxalate was formed in higher concentrations with increasing solution pH above the pK a of the radical species (CO 2 )˙ − . Below the pK a value, formate was the exclusive organic acid formed. The production of formate has a rather weak pH dependence but is enhanced slightly at a basic pH above 10. Furthermore, at basic pH, the effect of electrolyte concentration comes into play. Higher electrolyte concentrations, leading to shorter electrolyte Debye lengths, resulted in lowered organic acid yields. The highest yields of organic acids obtained in our system were 122 mg L −1 for oxalate and 77 mg L −1 for formate at an optimum 1 mM NaOH concentration in the starting solution. This work is a successful pioneering example of CO to organic acids conversion using non-thermal plasmas, which opens the pathway for a promising two-step conversion process of CO 2 to organic acids.

Sudagar, Alcina Johnson [Washington University in ↗
DOE OSTI2026

Solution structures in alkali nitrates and nitrites at high concentrations

In highly concentrated electrolyte solutions, where classical models often fail, specific ion–solvent interactions dictate bulk properties. Here, in this study, we combine small-angle X-ray scattering (SAXS) and Raman spectroscopy to link molecular coordination to mesoscale structure across alkali nitrate and nitrite solutions. Our results reveal a structural hierarchy driven by cation identity in that smaller cations (Li + , Na + ) form discrete contact ion pairs, while larger cations (K + , Rb + ) promote increasingly disordered local coordination and extended solute–solvent networks. Using nitrite salts as a structural control, we confirm this organizing principle across different anion geometries. Cs + undergoes a concentration-induced transition to a uniquely ordered state, resulting in a highly structured solution. This work provides a direct, mechanistic explanation for how cation choice dictates the ‘rules’ of solution architecture, offering a predictive basis for understanding phase behavior in complex industrial and environmental systems.

Mergelsberg, Sebastian T. [Pacific Northwest Natio↗
DOE OSTI2023

Multiscale modeling of solute diffusion in triblock copolymer membranes

We develop a multiscale simulation model for diffusion of solutes through porous triblock copolymer membranes. The approach combines two techniques: self-consistent field theory (SCFT) to predict the structure of the self-assembled, solvated membrane and on-lattice kinetic Monte Carlo (kMC) simulations to model diffusion of solutes. Solvation is simulated in SCFT by constraining the glassy membrane matrix while relaxing the brush-like membrane pore coating against the solvent. The kMC simulations capture the resulting solute spatial distribution and concentration-dependent local diffusivity in the polymer-coated pores; we parameterize the latter using particle-based simulations. We apply our approach to simulate solute diffusion through nonequilibrium morphologies of a model triblock copolymer, and we correlate diffusivity with structural descriptors of the morphologies. We also compare the model’s predictions to alternative approaches based on simple lattice random walks and find our multiscale model to be more robust and systematic to parameterize. Furthermore, our multiscale modeling approach is general and can be readily extended in the future to other chemistries, morphologies, and models for the local solute diffusivity and interactions with the membrane.

36 MATERIALS SCIENCE↗
DOE OSTI2023

On the number of stable solutions in the Kuramoto model

We consider a system of n coupled oscillators described by the Kuramoto model with the dynamics given by θ˙=ω+Kf(θ). In this system, an equilibrium solution θ∗ is considered stable when ω+Kf(θ∗)=0, and the Jacobian matrix Df(θ∗) has a simple eigenvalue of zero, indicating the presence of a direction in which the oscillators can adjust their phases. Additionally, the remaining eigenvalues of Df(θ∗) are negative, indicating stability in orthogonal directions. A crucial constraint imposed on the equilibrium solution is that |Γ(θ∗)|≤π, where |Γ(θ∗)| represents the length of the shortest arc on the unit circle that contains the equilibrium solution θ∗. We provide a proof that there exists a unique solution satisfying the aforementioned stability criteria. This analysis enhances our understanding of the stability and uniqueness of these solutions, offering valuable insights into the dynamics of coupled oscillators in this system.

Arenas, Alex (ORCID:0000000309370334)↗
DOE OSTI2024

An implicit solution for Asay foil trajectories generated by separable, sustained-production ejecta source models

We present a simple implicit solution for the time-dependent trajectory of a thin Asay foil ejecta diagnostic for the general case where the impinging ejecta cloud is generated by a source function characterized by an arbitrary (sustained) time dependence and a time-independent (stationary) particle velocity distribution. In the limit that the source function time dependence becomes a delta function, this solution—which is amenable to rapid numerical calculations of arbitrary accuracy—exactly recovers a previously published solution for the special case of instantaneous ejecta production. We also derive simple expressions for the free-surface arrival (catch-up) time as well as the true ejecta areal mass accumulation on the accelerating foil and place bounds on the level of error incurred when applying instant-production mass solutions to a sustained-production trajectory. We demonstrate these solutions with example calculations for hypothetical source functions spanning a wide range of ejecta production durations, velocity distributions, and temporal behaviors. These calculations demonstrate how the foil trajectory is often insensitive to the temporal dependence of the source function, instead being dominated by the velocity distribution. We quantify this insensitivity using a “compatibility score” metric. Under certain conditions, one may capitalize upon this insensitivity to obtain a good approximation of the second integral of the velocity distribution from the observed foil trajectory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Stability of exact solutions of the (2 + 1)-dimensional nonlinear Schrödinger equation with arbitrary nonlinearity parameter κ

In this work, we consider the nonlinear Schrödinger equation (NLSE) in 2+1 dimensions with arbitrary nonlinearity exponent κ in the presence of an external confining potential. Exact solutions to the system are constructed, and their stability as we increase the 'mass' (i.e., the L 2 norm) and the nonlinearity parameter κ is explored. Here we observe both theoretically and numerically that the presence of the confining potential leads to wider domains of stability over the parameter space compared to the unconfined case. Our analysis suggests the existence of a stable regime of solutions for all κ as long as their mass is less than a critical value M*(κ). Furthermore, we find that there are two different critical masses, one corresponding to width perturbations and the other one to translational perturbations. The results of Derrick's theorem are also obtained by studying the small amplitude regime of a four-parameter collective coordinate (4CC) approximation. A numerical stability analysis of the NLSE shows that the instability curve M*(κ) versus κ lies below the two curves found by Derrick's theorem and the 4CC approximation. In the absence of the external potential, κ = 1 demarcates the separation between the blowup regime and the stable regime. In this 4CC approximation, for κ < 1, when the mass is above the critical mass for the translational instability, quite complicated motions of the collective coordinates are possible. Energy conservation prevents the blowup of the solution as well as confines the center of the solution to a finite spatial domain. We call this regime the 'frustrated' blowup regime and give some illustrations. In an appendix, we show how to extend these results to arbitrary initial ground state solution data and arbitrary spatial dimension d.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Global centroid moment tensor solutions in a heterogeneous earth: the CMT3D catalogue

SUMMARY For over 40 yr, the global centroid-moment tensor (GCMT) project has determined location and source parameters for globally recorded earthquakes larger than magnitude 5.0. The GCMT database remains a trusted staple for the geophysical community. Its point-source moment-tensor solutions are the result of inversions that model long-period observed seismic waveforms via normal-mode summation for a 1-D reference earth model, augmented by path corrections to capture 3-D variations in surface wave phase speeds, and to account for crustal structure. While this methodology remains essentially unchanged for the ongoing GCMT catalogue, source inversions based on waveform modelling in low-resolution 3-D earth models have revealed small but persistent biases in the standard modelling approach. Keeping pace with the increased capacity and demands of global tomography requires a revised catalogue of centroid-moment tensors (CMT), automatically and reproducibly computed using Green's functions from a state-of-the-art 3-D earth model. In this paper, we modify the current procedure for the full-waveform inversion of seismic traces for the six moment-tensor parameters, centroid latitude, longitude, depth and centroid time of global earthquakes. We take the GCMT solutions as a point of departure but update them to account for the effects of a heterogeneous earth, using the global 3-D wave speed model GLAD-M25. We generate synthetic seismograms from Green's functions computed by the spectral-element method in the 3-D model, select observed seismic data and remove their instrument response, process synthetic and observed data, select segments of observed and synthetic data based on similarity, and invert for new model parameters of the earthquake’s centroid location, time and moment tensor. The events in our new, preliminary database containing 9382 global event solutions, called CMT3D for ‘3-D centroid-moment tensors’, are on average 4 km shallower, about 1 s earlier, about 5 per cent larger in scalar moment, and more double-couple in nature than in the GCMT catalogue. We discuss in detail the geographical and statistical distributions of the updated solutions, and place them in the context of earlier work. We plan to disseminate our CMT3D solutions via the online ShakeMovie platform.

58 GEOSCIENCES↗
DOE OSTI2022

Exact general solutions for cosmological scalar field evolution in a background-dominated expansion

We derive exact general solutions (as opposed to attractor particular solutions) and corresponding first integrals for the evolution of a scalar field $\phi$ in a universe dominated by a background fluid with equation of state parameter wB. In addition to the previously examined linear [V($\phi$)=V 0 $\phi$] and quadratic [V($\phi$)=V 0 $\phi$ 2 ] potentials, we show that exact solutions exist for the power-law potential V($\phi$)=V 0 $\phi$ n with n=4(1+w B )/(1-w B )+2 and n=2(1+w B )/(1-w B ). These correspond to the potentials V($\phi$)=V 0 $\phi$ 6 and V($\phi$)=V 0 $\phi$ 2 for matter domination and V($\phi$)=V 0 $\phi$ 10 and V($\phi$)=V 0 $\phi$ 4 for radiation domination. The $\phi$ 6 and $\phi$ 10 potentials can yield either oscillatory or nonoscillatory evolution, and we use the first integrals to determine how the initial conditions map onto each form of evolution. The exponential potential yields an exact solution for a stiff/kination (w B =1) background. We use this exact solution to derive an analytic expression for the evolution of the equation of state parameter, $w_\phi$, for this case.

79 ASTRONOMY AND ASTROPHYSICS↗