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

Topological symplectic Kondo effect

Multiple conduction channels interacting with a quantum impurity—a spin in the conventional “multichannel Kondo effect” or a topological mesoscopic device (“topological Kondo effect”)—has been proposed as a platform to realize anyonic quasiparticles. However, the above implementations require either perfect channel symmetry or the use of Majorana fermions. Here, in this work, we propose a Majorana-free mesoscopic setup which implements the Kondo effect of the symplectic Lie group and can harbor emergent anyons (including Majorana fermions, Fibonacci anyons, and parafermions) even in the absence of perfect channel symmetry. In addition to the detailed prescription of the implementation, we present the strong coupling solution by mapping the model to the multichannel Kondo effect associated to an internal symmetry and exploit conformal field theory to predict the nontrivial scaling of a variety of observables, including conductance, as a function of temperature. This work does not only open the door for robust Kondo-based anyon platforms, but also sheds light on the physics of strongly correlated materials with competing order parameters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Dissolution enables dolomite crystal growth near ambient conditions

Crystals grow in supersaturated solutions. A mysterious counterexample is dolomite CaMg(CO 3 ) 2 , a geologically abundant sedimentary mineral that does not readily grow at ambient conditions, not even under highly supersaturated solutions. Using atomistic simulations, we show that dolomite initially precipitates a cation-disordered surface, where high surface strains inhibit further crystal growth. However, mild undersaturation will preferentially dissolve these disordered regions, enabling increased order upon reprecipitation. Furthermore, our simulations predict that frequent cycling of a solution between supersaturation and undersaturation can accelerate dolomite growth by up to seven orders of magnitude. We validated our theory with in situ liquid cell transmission electron microscopy, directly observing bulk dolomite growth after pulses of dissolution. This mechanism explains why modern dolomite is primarily found in natural environments with pH or salinity fluctuations. More generally, it reveals that the growth and ripening of defect-free crystals can be facilitated by deliberate periods of mild dissolution.

Kim, Joonsoo↗
DOE OSTI2020

From Salt in Solution to Solely Ions: Solvation of Methyl Viologen in Deep Eutectic Solvents and Ionic Liquids

Solvation and transport properties of methlyviologen dichloride in mixtures of choline chloride (ChCl) and ethylene glycol (EG) including the deep eutectic solvent (DES) ethaline were studied through physical property measurements, cyclic voltammetry, and Raman spectroscopy. The transport properties of the solvate structures in ethaline and ChCl:EG mixtures were compared to a well-studied ionic liquid (IL), 1-methyl-1-propylpyrrolidinium bis(trifluoromethylsulfonyl)imide, [PYR 13 ][TFSI] through the lens of hole theory. Here by exchanging the chloride anion of ethaline with bis(trifuoromethylsufonyl)imide) anion, [TFSI], is found to increase the fluidity of the solvent and promoted stronger solute-solvent interactions. Raman spectroscopy shows that methylviologen is solvated by EG in solutions with ChTFSI but interstitially accommodated in void spaces in both ChCl solutions and [PYR 13 ][TFSI]. Here, complex solvents such as DESs are regarded as “designer solvents” and we demonstrate that the physical properties and solvation characteristics of these fluids can be varied by anion exchange of hydrogen bond acceptors.

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

Effective theories for a nonrelativistic field in an expanding universe: Induced self-interaction, pressure, sound speed, and viscosity

A massive, nonrelativistic scalar field in an expanding spacetime is usually approximated by a pressureless perfect fluid, which leads to the standard conclusion that such a field can play the role of cold dark matter. In this paper, we systematically study these approximations, incorporating subleading corrections. We provide two equivalent effective descriptions of the system, each of which offers its own advantages and insights: (i) A nonrelativistic effective field theory (EFT) with which we show that the relativistic corrections induce an effective self-interaction for the nonrelativistic field. As a byproduct, our EFT also allows one to construct the exact solution, including oscillatory behavior, which is often difficult to achieve from the exact equations. (ii) An effective (imperfect) fluid description, with which we demonstrate that, for a perturbed Friedmann-Lemaître- Robertson-Walker (FLRW) universe: (a) The pressure is small but nonzero (and positive), even for a free theory with no tree-level self-interactions. (b) The sound speed of small fluctuations is also nonzero (and positive), reproducing already known leading-order results, correcting a subdominant term, and identifying a new contribution that had been omitted in previous analyses. (c) The fluctuations experience a negative effective bulk viscosity. The positive sound speed and the negative bulk viscosity act in favor of and against the growth of overdensities, respectively. The net effect may be considered a smoking gun for ultra-light dark matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Numerical Study of Mineral Nucleation and Growth on a Substrate

We report nucleation and growth processes of minerals and other crystals can significantly affect one another due to the transport limitations and local depletion of reactive ions in the solution. Most numerical models and experimental measurements have typically focused on either growth or nucleation, but not both. In this work, we incorporate a heterogeneous nucleation process based on classical nucleation theory into a microcontinuum model that implements the Darcy–Brinkman–Stokes approach to study the interplay between nucleation and crystal growth on a substrate in diffusive systems. We demonstrate how the Damköhler number (reaction rate) and nucleation rate prefactor change the effective nucleation rate on a substrate. Higher surface growth rates deplete the solute concentration around the nuclei that appear initially on the substrate, creating islands that screen against further nucleation. The model predicts that measured nucleation rates may be affected by the history of crystal nucleation on the substrate. In the extreme case of high growth rates relative to diffusion, it predicts that the rate of subsequent nucleation is limited by reactant depletion. We introduce a nondimensional number α to represent the relation between surface propagation rate during growth and the heterogeneous nucleation rate. We show that it is important to control Damköhler number and α to achieve similar precipitation regimes at different reaction and nucleation rates. We suggest that the observed universality can guide the interpretation of experimental results on nucleation rates, since matching experiment can be achieved by tuning transport, reaction, and nucleation parameters simultaneously. In addition, we show how the bulk solution concentration affects the structure and topology of precipitation on a substrate.

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

When Do Extended Physics-Informed Neural Networks (XPINNs) Improve Generalization?

Physics-informed neural networks (PINNs) have become a popular choice for solving high-dimensional partial differential equations (PDEs) due to their excellent approximation power and generalization ability. Recently, extended PINNs (XPINNs) based on domain decomposition methods have attracted considerable attention due to their effectiveness in modeling multiscale and multiphysics problems and their parallelization. However, theoretical understanding of their convergence and generalization properties remains unexplored. In this study, we take an initial step towards understanding how and when XPINNs outperform PINNs. Specifically, for general multilayer PINNs and XPINNs, we first provide a prior generalization bound via the complexity of the target functions in the PDE problem and a posterior generalization bound via the posterior matrix norms of the networks after optimization. Moreover, based on our bounds, we analyze the conditions under which XPINNs improve generalization. Concretely, our theory shows that the key building block of XPINN, namely, the domain decomposition, introduces a tradeoff for generalization. On the one hand, XPINNs decompose the complex PDE solution into several simple parts, which decreases the complexity needed to learn each part and boosts generalization. On the other hand, decomposition leads to less training data being available in each subdomain, and hence such a model is typically prone to overfitting and may become less generalizable. Empirically, we choose five PDEs to show when XPINNs perform better than, similar to, or worse than PINNs, hence demonstrating and justifying our new theory.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Theory meets experiment: Insights into structure and magnetic properties of Fe 1-x Ni x B alloy

Here, we studied the structural and magnetic properties of the solid solution Fe 1-x Ni x B through theoretical and experimental approaches. Powder X-ray diffraction, X-ray Pair Distribution Function analysis, and energy dispersive X-ray spectroscopy reveal that the Fe 1-x Ni x B solid solution crystallizes in the β-FeB structure type up to x = 0.6–0.7 and exhibits anisotropic unit cell volume contraction with increasing Ni concentration. Magnetic measurements showed a transition from ferromagnetism to paramagnetism around x = 0.7. For x = 0.5, the low (< 0.3 μ B ) magnetic moments suggest itinerant magnetism despite the relatively high Curie temperature (up to 225 K). Theoretical calculations indicated different types of magnetic orderings depending on the Fe/Ni atomic order, with the antiferromagnetic state being stable for ordered FeNiB 2 , whereas the ground state is ferromagnetic for the disordered alloy. Calculations also predicted the coexistence of low- and high-spin states in Fe atoms around the composition with x = 0.5, in line with the experimental evidence from 57 Fe Mössbauer spectroscopy. The two magnetically distinct Fe sites for x = 0.3, 0.4, and 0.5 observed by 57 Fe Mössbauer spectroscopy can also be interpreted as two magnetically different regions or clusters. The formation of these clusters could affect the critical behavior near a quantum magnetic transition based on a potential ferromagnetic quantum critical point identified computationally and experimentally near x = 0.64. This work highlights the complex interplay between structure and magnetism in Fe 1-x Ni x B alloys, suggesting areas for future research on quantum critical behavior.

Bhaskar, Gourab [Iowa State Univ., Ames, IA (Unite↗
DOE OSTI2023

Exceptional enhancement of mechanical properties in high-entropy alloys via thermodynamically guided local chemical ordering

Understanding the local chemical ordering propensity in random solid solutions, and tailoring its strength, can guide the design and discovery of complex, paradigm-shifting multicomponent alloys. First, we present a simple thermodynamic framework, based solely on binary enthalpies of mixing, to select optimal alloying elements to control the nature and extent of chemical ordering in high-entropy alloys (HEAs). Next, we couple high-resolution electron microscopy, atom probe tomography, hybrid Monte-Carlo, special quasirandom structures, and density functional theory calculations to demonstrate how controlled additions of Al and Ti and subsequent annealing drive chemical ordering in nearly random equiatomic face-centered cubic CoFeNi solid solution. We establish that short-range ordered domains, the precursors of long-range ordered precipitates, inform mechanical properties. Specifically, a progressively increasing local order boosts the tensile yield strengths of the parent CoFeNi alloy by a factor of four while also substantially improving ductility, which breaks the so-called strength–ductility paradox. Finally, we validate the generality of our approach by predicting and demonstrating that controlled additions of Al, which has large negative enthalpies of mixing with the constituent elements of another nearly random body-centered cubic refractory NbTaTi HEA, also introduces chemical ordering and enhances mechanical properties.

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

Symmetry Determining Equations of the Rankine-Hugoniot Equations for Variable Velocity Shock Waves

The “constant velocity piston” problem (Fig. 1), also known as the “piston problem,” is a standard model for a one dimensional, in our case linear, symmetric shock wave moving through an inviscid, perfect gas. The model can be divided into two regions - a perturbed section on the left and an unperturbed section on the right - by a moving shock wave moving left to right. Both the perturbed and unperturbed sections, i.e. the shocked and unshocked regions, respectively, obey the Eulerian conservation equations; however, at the exact location of the shock, there is a mathematical discontinuity not satisfied by the Euler equations. To ensure continuity and conservation of certain quantities when crossing between the unshocked and shocked regions, we evoke a series of equations derived from the Eulerian conservation equations, called the Rankine-Hugoniot equations, or “jump” equations as it is often referred to in the literature on the topic. The classical constant-velocity piston problem assumes the piston features a constant driving velocity (among many other willing suspensions of belief required in the pursuit of a first principles equation model); consequent to this assumption is a constant-velocity shock and a constant-velocity shocked flow state. However, using Lie Group Theory (LGT), also known as symmetry analysis, we can attempt to reinterpret the model with a shock wave of variable velocity in time and space. An extension of the model in this way opens up the possibility for obtaining new analytical solutions to the piston problem for certain shock velocity models. In this report, we use LGT to derive the symmetry determining equations (SDEs), whose solutions are Lie groups, which permit analytical solutions. In the future, we can then use the SDEs to define constraint equations on the shock velocity model and what the successive solutions to the Euler equations might be based off such constraints. This report is structured as follows: Section 2 provides a brief derivation of the Rankine-Hugoniot (“jump”) equations; Section 3 gives an overview of Lie group theory; Section 4 derives the SDEs of the jump equations; Section 5 derives the Euler conservation equations for fluids; and Section 6 presents concluding remarks and opportunities for future studies.

42 ENGINEERING↗
DOE OSTI2024

Bespoke Liquid/Liquid Interfaces (Final Technical Report)

The goal of DE-SC0001815 was to advance the basic science of liquid:liquid interface formation, to develop a deeper understanding of the mechanisms of phase separation and the essential relationships between solution composition, organization and dynamics that underlie the kinetic regime of solvent extraction. This included learning how interfacial organization and dynamics alters the properties of the primary coordination sphere of ions and the free energy of transport of ions complexes across a phase boundary. We relied primarily upon classical molecular dynamics studies to determine the equilibrium ensembles of these complex systems, but also utilized ab-initio MD and cluster-based density functional theory (DFT) calculations when more detailed investigation of the electronic structure was needed. We continued development of graph-theory based analyses to elucidate hierarchical correlations and expanded into geometric topology methods to quantify the collectively organized structures that can organize at a liquid/liquid interface during solute transport. One of the main conclusions was from the observation of two distinct mechanisms for solute transport - those that derive from amplifications of interfacial heterogeneity and surface roughness, and those wherein surface roughness has been dampened and instead collectively organized macrostructures work to bring solutes into the organic phase. It was our aim to create a concrete chemical model of the underlying driving forces behind interfacial primary and secondary structure formation and to map out the energetic features of solute transport so that tailored liquid/liquid can be developed that have characteristic kinetic features associated with mass transport.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Osmotic control of the spacing of parallel shear cracks in shale growing subcritically in geologic past

The geological genesis of natural cracks in sedimentary rocks such as shale is a problem that needs to be understood to improve the technology of hydraulic fracturing as well as deep sequestration of harmful fluids. Why are the vertical natural cracks roughly parallel and equidistant, and why is the spacing roughly 10 cm rather than 1 cm or 100 cm? Fracture mechanics of critical cracks cannot answer this question. Neither can the material heterogeneity. The growth of critical parallel cracks is impossible because the relative crack face displacements would immediately localize into one crack, leading to an earthquake. The cracks must have formed, on the tectonic time scale, by a slow growth of subcritical shear cracks governed by the Charles-Evans law. The idea advanced here is that what controls the crack spacing is the balance between the reduction, due to shear dilatancy, of the concentration of ions such as Na + and Cl - in each fracture process zone (PFZ), which decelerates the cracks, and the restoration of ion concentration by diffusion of ions from the space between the cracks into the FPZ. This diffusion of water is driven mainly by the osmotic pressure gradient, which offsets the deceleration and depends strongly on the crack spacing. A simple analytical solution of the steady state is rendered possible by approximating the ion concentration profiles between adjacent cracks by parabolic arcs. Applying this theory to Woodford shale yields the approximate crack spacing of 10 cm, which is realistic. Furthermore, the stability of unlimited parallel mode II frictional crack growth is proven by examining the second variation of the free energy. Water concentration drop in the FPZ due to shear dilatancy and its restoration by water diffusion from the inter-crack space have similar effect, although probably much weaker.

42 ENGINEERING↗
DOE OSTI2022

Vibrational Entropy Contribution to Mixing Free Energy of Ni-Rich LiNi 1- y Co y O 2

First principles density-functional-theory calculations are performed of the mixing free energy of the layered oxide LiNi 1-y Co y O 2 , a model Li-ion-battery cathode material, at 1000 K. Analysis of internal energy, LN exchange, configurational entropy and vibrational entropy are combined to obtain the net mixing free energy. The vibrational entropy requires a more indirect treatment than the other terms; analysis based on heat capacity is presented here, to show that the vibrational entropy is of comparable magnitude to the configurational entropy, and it tends to promote solid-solution behavior in LiNi 1-y Co y O 2 .

25 ENERGY STORAGE↗
DOE OSTI2024

In situ Formation of Platinum-Carbon Catalysts in Propane Dehydrogenation

The catalytic production of propylene via propane dehydrogenation (PDH) is a key reaction in the chemical industry. By combining operando transmission electron microscopy with density functional theory analysis, we show that the intercalation and ordering of carbon on Pt interstitials to form Pt-C solid solutions is relevant for increasing propylene production. More specifically, we found that at the point of enhanced propylene formation, the structure of platinum nanoparticles is transformed into a transient caesium chloride-type Pt-C polymorph. At more elevated temperatures, the zincblende and rock salt polymorphs seemingly coexist. When propylene production was highest, multiple crystal structures consisting of Pt and carbon were occasionally found to coexist in one individual nanoparticle, distorting the Pt lattice. Catalyst coking was detected at all stages of the reaction, but did initially not affect all particles. These findings could lead to the development of novel synthesis strategies towards tailoring highly efficient PDH catalysts.

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

General solution to the Kohn–Luttinger nonconvergence problem

A simple, but general solution is proposed for the Kohn–Luttinger problem, i.e., the nonconvergence of the finite-temperature many-body perturbation theory with its zero-temperature counterpart as temperature is lowered to zero under some circumstances. How this nonconvergence can be avoided by altering the reference wave function is illustrated numerically by using up to the fifth order of the perturbation theory.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Opposing Effects of Impurity Ion Sr 2+ on the Heterogeneous Nucleation and Growth of Barite (BaSO 4 )

Nucleation of barite (BaSO 4 ) has broad implications in geological, environmental, and materials sciences. While impurity metals are common, our understanding of how they impact nucleation remains dim. In this work, we used classical optical microscopy compared to fast X-ray nanotomography (XnT) to investigate heterogeneous nucleation of barite on silica in situ with Sr 2+ as an impurity ion. The observed barite nucleation rates were consistent with classical nucleation theory (CNT), where barite crystals displayed a nonuniform size distribution, exhibiting distinct morphologies and incubation periods in Sr-free solutions. While undetectable with optical microscopy, nanotomography revealed that addition of Sr2+ enhances nucleation rates driven by the pre-factor in CNT, likely because both adsorbed Ba 2+ and Sr 2+ act as precursor sites on which nucleation occurs. Sr 2+ simultaneously inhibits growth, however, leading to a homogeneous distribution of smaller crystals. This finding will enable an improved predictive understanding of nucleation in natural and synthetic environments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

N 5 -Scaling Excited-State-Specific Perturbation Theory

In this paper, we show that by working in a basis similar to that of the natural transition orbitals and using a modified zeroth-order Hamiltonian, the cost of a recently introduced perturbative correction to excited-state mean field theory can be reduced from seventh to fifth order in the system size. The (occupied) 2 (virtual) 3 asymptotic scaling matches that of ground-state second-order Møller-Plesset theory but with a significantly higher prefactor because the bottleneck is iterative: it appears in the Krylov-subspace-based solution of the linear equation that yields the first-order wave function. Here, we discuss the details of the modified zeroth-order Hamiltonian we use to reduce the cost and the automatic code generation process we used to derive and verify the cost scaling of the different terms. Overall, we find that our modifications have little impact on the method's accuracy, which remains competitive with singles and doubles equation-of-motion coupled cluster.

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

Hydrogen Stabilized RhPdH 2D Bimetallene Nanosheets for Efficient Alkaline Hydrogen Evolution

The design of catalysts with high activity and robust stability for alkaline hydrogen evolution reaction (HER) remains a great challenge. Here, we report an efficient catalyst of two-dimensional bimetallene hydrides, in which H atoms stabilize the rhodium palladium bimetallene. The system exists because of the introduction of H that is in-situ chemically released from the formaldehyde solution during the synthesis. This provides a highly stable catalyst based on an unstable combination of metal elements. Density functional theory calculations show the H is confined by electronic interactions and the Miedema rule of reverse stability of the RhPd alloy.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Achieving high tensile strength and ductility in refractory alloys by tuning electronic structure

The energy efficiency of heat engines (gas and steam turbines) for electricity production and propulsion is determined by the Carnot cycle and scales with operating temperature. Commercial nickel- and cobalt-based superalloys melt near 1,500 °C and rapidly lose mechanical strength beyond 1,000 °C. Refractory metals melt well above 2,000 °C but have inherent manufacturability challenges that are barriers to adoption, such as high ductile-to-brittle transition temperatures. Using density functional theory-guided design, we demonstrate tailored local lattice distortions that promote phase-stable, non-equiatomic refractory concentrated solid solutions with both high ductility and strength. Here, we exemplify this for single-phase, body-centred cubic Nb 4 Ta 4 V 3 Ti that exhibits castability, excellent room-temperature tensile yield strength (∼1 GPa) and ductility (approaching 20% uniform strain), and exceptional high-temperature tensile strength (500 MPa at 1,000 °C). These findings illustrate a path for designing materials that hold great potential for advancing next-generation technologies such as Generation IV fission reactors, first-generation fusion-plasma reactors, and more efficient gas turbines for electricity generation and propulsion.

DFT↗