Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solution theory”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Effect of surface properties and polymer chain length on polymer adsorption in solution

In polymer nanoparticle composites (PNCs) with attractive interactions between nanoparticles (NPs) and polymers, a bound layer of the polymer forms on the NP surface, with significant effects on the macroscopic properties of the PNCs. The adsorption and wetting behaviors of polymer solutions in the presence of a solid surface are critical to the fabrication process of PNCs. In this study, we use both classical density functional theory (cDFT) and molecular dynamics (MD) simulations to study dilute and semi-dilute solutions of short polymer chains near a solid surface. Using cDFT, we calculate the equilibrium properties of polymer solutions near a flat surface while varying the solvent quality, surface–fluid interactions, and the polymer chain lengths to investigate their effects on the polymer adsorption and wetting transitions. Using MD simulations, we simulate polymer solutions near solid surfaces with three different curvatures (a flat surface and NPs with two radii) to study the static conformation of the polymer bound layer near the surface and the dynamic chain adsorption process. Additionally, we find that the bulk polymer concentration at which the wetting transition in the poor solvent system occurs is not affected by the difference in surface–fluid interactions; however, a threshold value of surface–fluid interaction is needed to observe the wetting transition. We also find that with good solvent, increasing the chain length or the difference in the surface–polymer interaction relative to the surface–solvent interaction increases the surface coverage of polymer segments and independent chains for all surface curvatures. Finally, we demonstrate that the polymer segmental adsorption times are heavily influenced only by the surface–fluid interactions, although polymers desorb more quickly from highly curved surfaces.

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

Exact two-body expansion of the many-particle wave function

Progress toward the solution of the strongly correlated electron problem has been stymied by the exponential complexity of the wave function. Previous work established an exact two-body exponential product expansion for the ground-state wave function. By developing a reduced density-matrix analog of Dalgarno-Lewis perturbation theory, we prove here that (i) the two-body exponential product expansion is rapidly and globally convergent with each operator representing an order of a renormalized perturbation theory, (ii) the energy of the expansion converges quadratically near the solution, and (iii) the expansion is exact for both ground and excited states. The two-body expansion offers a reduced parametrization of the many-particle wave function as well as the two-particle reduced density matrix with potential applications on both conventional and quantum computers for the study of strongly correlated quantum systems. In this work, we demonstrate the result with the exact solution of the contracted Schrödinger equation for the molecular chains H 4 and H 5 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Molecular thermodynamics for scaling prediction: Case of membrane distillation

Membrane distillation (MD) for water treatment is significantly impaired by the scaling of dissolved minerals. The type and content of minerals generally measured as total dissolved solids (TDS) in hypersaline brines not only reduce the MD flux but also control the scaling behavior on the membrane surface. The scaling-induced pore blockage further reduces water flux and eventually leads to membrane wetting. The scaling problem is even more pronounced in the treatment of produced water (PW) as it contains 3–7 times higher TDS concentrations, compared to seawater. Theoretically, the necessary conditions for a salt to precipitate can be traced from its solubility product constant and activity of the constituents within the solution. Therefore, a comprehensive thermodynamic model is necessary to represent the electrolyte behavior and to predict the precipitation of different salts in a complex solution like PW. We pursued electrolyte Nonrandom Two Liquid Theory (eNRTL), a state-of-the-art electrolyte model, to address the PW fluid phase equilibria. With a fully parameterized eNRTL model, we predicted salt precipitations in two different PW samples and compared the results against the experimental findings. Furthermore, here we forecast the precipitation of salts in response to the change in PW concentration and temperature in the MD operation.

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

Dynamic Barriers to Crystallization of Calcium Barium Carbonates

Metastable carbonates play important roles in geochemistry, biomineralization and serve as model systems for nonclassical theories of nucleation and growth. Balcite (Ca 0.5 Ba 0.5 CO 3 ) is a remarkable carbonate phase that is isostructural with a high-temperature modification of calcite (CaCO 3 ), yet can be synthesized at ambient conditions. Here, we investigate crystallization pathways in the Ba-Ca-CO 3 -H 2 O system, with a focus on the transformation of amorphous calcium barium carbonate (ACBC) to balcite and subsequent decomposition into the equilibrium calcite (CaCO3) and witherite (BaCO 3 ) phases. Density functional theory calculations show that balcite is an unstable solid solution (Ca 1-x BaxCO 3 , R 3¯ m) in the range 0.17 < x < 0.5, but is accessible through the amorphous ACBC precursor for x ≲ 0.5, and predict its decomposition into calcite and witherite. We confirm this pathway experimentally but found demixing to proceed slowly and remain incomplete even after 9 months. Nucleation kinetics of balcite from ACBC was assessed using a microfluidic assay, where increasing barium content led to a surprising increase in the balcite nucleation rate, despite decreasing thermodynamic driving force. We attribute crystallization rates that dramatically accelerate with time to changes in interfacial structure and composition during coarsening of the amorphous precipitate. By carefully quantifying the thermodynamic and kinetic contributions in the multistep crystallization of a metastable carbonate, we produce insights that allow us to better interpret the formation and persistence of metastable minerals in natural and synthetic environments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

First-Principles Modeling of the Repassivation of Corrosion Resistant Alloys: Part I. O and Cl Adsorption Energy

To understand the mechanism of repassivation, a key process that controls the performance of corrosion resistant alloys (CRAs), it is useful to obtain some fundamental parameters that describe the atomistic processes. Oxygen and chlorine adsorption energies were systematically calculated with density functional theory for Ni-22Cr alloys with seven different solute elements: Co, Cu, Fe, Mn, Mo, Ru and W. All elements other than Cu enhance both O and Cl adsorption to alloy surfaces. Cr has the strongest effect among these elements. Surface configurations that give stronger O adsorption also tend to give stronger Cl adsorption. While this makes it difficult to differentiate the effect of these solute elements on the oxidation/repassivation tendency of alloys, these fundamental properties provide the foundation for higher scale models that could be used to study the oxidation process during corrosion and potentially inform the design of CRAs. Bader charge analysis reveals a strong relationship between the adatom adsorption energy and the Bader charges of both O and Cl adatom. In conclusion, the work function of alloy surfaces is increased by O while decreased by Cl, suggesting that Cl adsorption can have an intrinsic negative impact on corrosion resistance by decreasing the barrier to electron transfer.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Improved measurements of nonlinear integrable optics at IOTA

Nonlinear integrable optics (NIO) are a promising novel approach at improving the stability of high intensity beams. Implementations of NIO based on specialized magnetic elements are being tested at the Integrable Optics Test Accelerator (IOTA) at Fermilab. One method of verifying proper implementation of these solutions is by measuring the analytic invariants predicted by theory. The initial measurements of nonlinear invariants were performed during IOTA run in 2019/20, however the covid-19 pandemic prevented the full-scale experimental program from being completed. Several important improvements were implemented in IOTA for the 2022/23 run, including the operation at higher beam energy of 150 MeV, improved optics control, and chromaticity correction. This report presents on the improved calibrations of the NIO for nonlinear invariant measurements.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Improved Measurements of Nonlinear Integrable Optics Invariants at IOTA

Nonlinear integrable optics (NIO) are a promising novel approach at improving the stability of high intensity beams. Implementations of NIO based on specialized magnetic elements are being tested at the Integrable Optics Test Accelerator (IOTA) at Fermilab. One method of verifying proper implementation of these solutions is by measuring the analytic invariants predicted by theory. The initial measurements of nonlinear invariants were performed during IOTA run in 2019/20, however the covid-19 pandemic prevented the full-scale experimental program from being completed. Several important improvements were implemented in IOTA for the 2022/23 run, including the operation at higher beam energy of 150 MeV, improved optics control, and chromaticity correction. This report presents the results of improved measurements of nonlinear invariants.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

On the choice of reference orbitals for linear-response calculations of solution-phase K-edge X-ray absorption spectra

The simplest response theory methods for computing vertical excitation spectra in condensed-phase are configuration interaction with single excitations (CIS) and linear-response time-dependent density functional theory (TDDFT) within the Tamm–Dancoff approximation. In applications to X-ray absorption spectroscopy (XAS), methods like CIS and TDDFT that codify only single excitations into the wave function are prone to catastrophic errors in main-edge and post-edge features whose shapes act as a crucial fingerprint in structural analyses of liquids. Here, we show that these errors manifest primarily due to a lack of orbital relaxation in conventional linear-response theories and that core-ionized (n - 1-electron) references, like those of electron-affinity TDDFT, can eliminate the errors in the spectral profile, even in the highest-energy parts of the post-edge. Crucially, we find that single excitations atop core-ionized references are sufficient to elucidate liquid-phase XAS spectra with semi-quantitative accuracy, opening the door for methods like electron-affinity CIS/TDDFT to be used as efficient alternatives to higher-order wave function approaches.

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

Numerical solution of large scale Hartree–Fock–Bogoliubov equations

The Hartree–Fock–Bogoliubov (HFB) theory is the starting point for treating superconducting systems. However, the computational cost for solving large scale HFB equations can be much larger than that of the Hartree–Fock equations, particularly when the Hamiltonian matrix is sparse, and the number of electrons N is relatively small compared to the matrix size N b . We first provide a concise and relatively self-contained review of the HFB theory for general finite sized quantum systems, with special focus on the treatment of spin symmetries from a linear algebra perspective. We then demonstrate that the pole expansion and selected inversion (PEXSI) method can be particularly well suited for solving large scale HFB equations. For a Hubbard-type Hamiltonian, the cost of PEXSI is at most $\mathcal{O}$(N b 2 ) for both gapped and gapless systems, which can be significantly faster than the standard cubic scaling diagonalization methods. We show that PEXSI can solve a two-dimensional Hubbard-Hofstadter model with N b up to 2.88 × 10 6 , and the wall clock time is less than 100 s using 17 280 CPU cores. Finally, this enables the simulation of physical systems under experimentally realizable magnetic fields, which cannot be otherwise simulated with smaller systems.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2020

Progress in Understanding the Origins of Excellent Corrosion Resistance in Metallic Alloys: From Binary Polycrystalline Alloys to Metallic Glasses and High Entropy Alloys

Some of the factors responsible for good corrosion resistance of select polycrystalline and emerging alloys in chloride solutions are discussed with a goal of providing some perspectives on the current status and future directions. Traditional metallic glass alloys, single phase high entropy alloys (HEAs), early metallic glasses, and high entropy metallic glasses are all emerging corrosion-resistant alloys (CRAs) that utilize traditional strategies for improved corrosion resistance as well as take advantage of some other novel beneficial attributes. These materials enjoy many degrees of freedom as far as choice of both composition and structure, providing great flexibility in the pursuit of superior corrosion resistance. The new materials depart from classical solvent-solute type polycrystalline binary or ternary alloys. Thus, such emerging materials provide significant opportunities to achieve even greater improvements in corrosion resistance in harsh environments. Several examples of the unique corrosion properties of selected materials in the context of modern theories of corrosion are discussed herein. Discussion is restricted to solid-solution binary or ternary polycrystalline alloys, several metallic glass alloys, and single phase HEAs. A common feature of many CRAs is that composition and microstructure often affect both passivity and resistance to localized corrosion that can be divided into initiation, stabilization, and propagation stages. Enormous complexities in protective oxide structures and chemistries and the large number of combinatorial possibilities in newer materials such as HEAs preclude trial-and-error approaches and perhaps even combinatorial experimental design. Computational materials methodologies will be required in the search for new corrosion-resistant alloys in these material classes. The search must consider the best scientific insights available regarding how major and minor alloy additions, as well as various microstructural attributes, contribute to corrosion mitigation. Additional scientific insights, as they emerge, will enable choices beyond the reliance on high concentrations of alloying elements that are known to affect passivity breakdown and pit stabilization. A challenge is to connect the “basic attributes” of an alloy with its properties. The strength of this connection will likely require new scientific principles enabling deep multiphysics insights in order to link feature(s) such as composition and metallurgical phases to the desired corrosion properties. Application of data informatics will likely also play a role given the plethora of variables that are important in corrosion and the difficulty in assessing all relationships. Here, the opportunity exists to accelerate the design of emerging materials for high corrosion resistance.

36 MATERIALS SCIENCE↗
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 OSTI2020

C–H Bond Activation via U(II) in the Reduction of Heteroleptic Bis(trimethylsilyl)amide U(III) Complexes

Reduction of (C 5 Me 5 ) 2 U III (NR 2 ) and (C 5 Me 5 )U III (NR 2 ) 2 (R = SiMe 3 ) with potassium graphite in the presence of 2.2.2-cryptand (crypt) generates dark solutions that have UV-visible spectra consistent with time-dependent density functional theory (TDDFT) calculations on the U(II) products, [(C 5 Me 5 ) 2 U II (NR 2 )] - and [(C 5 Me 5 )U II (NR 2 ) 2 ] - . However, the solutions quickly change color and form the U(III) C–H bond activation products [K(crypt)][(C 5 Me 5 ) 2 U III (CH 2 SiMe 2 NSiMe 3 –κC,κN)], 1, and [K(crypt)][(C 5 Me 5 )U III (NR 2 )(CH 2 SiMe 2 NSiMe 3 –κC,κN)], 2, that were identified by X-ray crystallography. DFT calculations on the putative [(C 5 Me 5 ) 2 U II (NR 2 )] - and [(C 5 Me 5 )U II (NR 2 ) 2 ] - complexes revealed 5f 3 6d 1 ground state electron configurations as previously found in isolable [(C 5 H 4 SiMe 3 ) 3 U II ] - , which indicated that these low symmetry heteroleptic complexes are reasonable precursors for new U(II) complexes.

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

Benchmarking magnetised three-wave coupling for laser backscattering: analytic solutions and kinetic simulations

Understanding magnetised laser–plasma interactions is important for controlling magneto-inertial fusion experiments and developing magnetically assisted radiation and particle sources. For nanosecond pulses at non-relativistic intensities, interactions are dominated by coherent three-wave interactions, whose nonlinear coupling coefficients became known only recently when waves propagate at oblique angles with the magnetic field. In this paper, backscattering coupling coefficients predicted by warm-fluid theory are benchmarked using particle-in-cell simulations in one spatial dimension, and excellent agreements are found for a wide range of plasma temperatures, magnetic field strengths and laser propagation angles, when the interactions are mediated by electron-dominant hybrid waves. Systematic comparisons between theory and simulations are made possible by a rigorous protocol. On the theory side, the initial boundary value problem of linearised three-wave equations is solved, and the transient-time solutions allow the effects of growth and damping to be distinguished. On the simulation side, parameters are carefully chosen and calibration runs are performed to ensure that comparisons are well controlled. Fitting simulation data to analytical solutions yields numerical growth rates that match theory predictions within error bars. Although warm-fluid theory is found to be valid for a wide parameter range, genuine kinetic effects have also been observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Alternative formulation of weak magnetohydrodynamic turbulence theory

In a recent paper, the weak turbulence theory for incompressible magnetohydrodynamics is formulated by employing the method customarily applied in the context of kinetic weak plasma turbulence theory. Such an approach simplified certain mathematical procedures including achieving the closure relationship. The formulation in the above-cited paper starts from the equations of incompressible magnetohydrodynamic (MHD) theory expressed via Elsasser variables. The derivation of nonlinear wave kinetic equation therein is obtained via a truncated solution at the second-order of iteration following the standard practice. In the present paper, the weak MHD turbulence theory is alternatively formulated by employing the pristine form of incompressible MHD equation rather than that expressed in terms of Elsasser fields. The perturbative expansion of the nonlinear momentum equation is carried out up to the third-order iteration rather than imposing the truncation at the second order. Furthermore, it is found that while the resulting wave kinetic equation is identical to that obtained in the previous paper cited above, the third-order nonlinear correction plays an essential role for properly calculating derived quantities such as the total and residual energies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Density-functional theory of material design: fundamentals and applications-I

This article is part-I of a review of density-functional theory (DFT) that is the most widely used method for calculating electronic structure of materials. The accuracy and ease of numerical implementation of DFT methods has resulted in its extensive use for materials design and discovery and has thus ushered in the new field of computational material science. In this article, we start with an introduction to Schrödinger equation and methods of its solutions. After presenting exact results for some well-known systems, difficulties encountered in solving the equation for interacting electrons are described. How these difficulties are handled using the variational principle for the energy to obtain approximate solutions of the Schrödinger equation is discussed. The resulting Hartree and Hartree–Fock theories are presented along with results they give for atomic and solid-state systems. We then describe Thomas–Fermi theory and its extensions which were the initial attempts to formulate many-electron problem in terms of electronic density of a system. Having described these theories, we introduce modern DFT by discussing Hohenberg–Kohn theorems that form its foundations. We then go on to discuss Kohn–Sham (KS) formulation of DFT in its exact form. Next, local density approximation (LDA) is introduced and solutions of KS equation for some representative systems, obtained using the LDA, are presented. We end part-I of the review describing the contents of part-II.

36 MATERIALS SCIENCE↗
DOE OSTI2026

KBKit: A Python Toolkit for Kirkwood–Buff Theory from Molecular Dynamics

Thermodynamic properties of liquid mixtures govern processes that range from drug delivery to energy storage, yet extracting these properties from molecular simulations remains challenging. Kirkwood–Buff (KB) theory offers a rigorous route by linking microscopic pair distribution functions to macroscopic free energies, but practical use of the theory has been hindered by two obstacles: (i) the long simulations needed to obtain well-converged Kirkwood-Buff integrals (KBIs) and (ii) the specialized corrections required to translate finite-size data to the thermodynamic limit. $\texttt{KBKit}$ is an open-source Python package that removes these barriers. It automatically computes KBIs and derived thermodynamic quantities from GROMACS input files, applies state-of-the-art finite-size corrections, and provides built-in diagnostic tools to quantify statistical uncertainty. Written with modern software-engineering practices—continuous integration, extensive unit testing, and thorough documentation—$\texttt{KBKit}$ is both reliable and easy to extend. By condensing complex KBI analysis into a few intuitive commands, $\texttt{KBKit}$ enables researchers to incorporate KB theory into routine simulation workflows and accelerate the discovery of solution-phase thermodynamics.

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

Wilson loops in 5d long quiver gauge theories

Quiver gauge theories with a large number of nodes host a wealth of Wilson loop operators. Expectation values are obtained, using supersymmetric localization, for Wilson loops in the antisymmetric representations associated with each individual gauge node, for a sample of 5d long quiver gauge theories whose UV fixed points have holographic duals in Type IIB. The sample includes the $T_N$ theories and the results are uniformly given in terms of Bloch-Wigner functions. The holographic representation of the Wilson loops is identified. It comprises, for each supergravity solution, a two-parameter family of D3-branes which exactly reproduce the field theory results and identify points in the internal space with the faces of the associated 5-brane web. The expectation values of (anti)fundamental Wilson loops exhibit an enhanced scaling for many operators, which matches between field theory and supergravity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗