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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Elucidating and contrasting the mechanisms for Mg and Ca sulfate ion-pair formation with multi-level embedded quantum mechanics/molecular dynamics simulations

Here, solutions and minerals containing sulfate (SO 4 2- ), and Ca 2+ and Mg 2+ cations, are ubiquitous throughout the lithosphere and are significant components of seawater, thus presenting a prototypical system for the study of strong electrolytes and crystal nucleation mechanisms. However, despite their relative abundance, key questions remain unanswered about the most fundamental atomic-level steps of their mineralization pathways and aqueous dynamics. Here, we carry out enhanced sampling multi-level molecular dynamics (MD) embedded correlated wavefunction theory simulations to elucidate ion-pairing mechanisms for Mg–SO 4 and Ca–SO 4 in concentrated aqueous solution, accurately capturing effects arising from both structural dynamics and electron exchange–correlation. We predict contact-ion-pair formation to be barrierless and highly exoergic for Ca–SO 4 , in agreement with its minimal solubility, whereas for Mg–SO 4 , solvent-shared and contact ion pairs have similar free energies, qualitatively consistent with its higher solubility. Finally, we demonstrate that brief high-temperature pre-equilibration may be utilized to accelerate convergence of free energies in blue-moon-ensemble enhanced-sampling MD.

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

Exact solutions of burst pressure for thick-walled cylinders in power-law strain hardening steels

Exact solutions of burst pressure for defect-free, thick-walled cylindrical pressure vessels with capped ends are developed using the flow theory of plasticity in terms of the Tresca, von Mises, and Zhu-Leis yield criteria. The material is assumed to obey a power-law strain hardening rule, and the finite strain theory is adopted to describe large plastic deformation. On this basis, internal pressure is obtained as a complex polynomial function of effective strains on the inside and outside surfaces of the pressure vessel with the use of the Bernoulli numbers. At burst failure, the effective strains and stresses on the inside and outside surfaces are first determined, and then three flow solutions of burst pressure are obtained as a power series function of the diameter ratio (D o /D i ), strain hardening exponent (n), and ultimate tensile strength (UTS). The power series solution is confirmed to agree well with the exact flow solution of burst pressure in a closed-form for each yield criterion, and the results show that the von Mises flow solution is an upper bound prediction, the Tresca flow solution is a lower bound prediction, and the Zhu-Leis flow solution is an intermediate prediction of burst pressure. In conclusion, two sets of full-scale burst test data are then utilized to evaluate and validate the proposed flow solutions of burst pressure for both thin-walled pipes and thick-walled cylindrical pressure vessels.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Symplectic effective field theory for nuclear structure studies

Here, a Symplectic Effective Field Theory that unveils the observed emergence of symplectic symmetry in atomic nuclei is advanced. Specifically, starting from a simple extension of the harmonic-oscillator Lagrangian, an effective field theory applied against symplectic basis states is shown to yield a Hamiltonian system with one fitted parameter. The scale of the system can be determined self consistently as the ratio of the average volume of a nucleus assumed to be spherical to its volume as determined by the average number of oscillator quanta, which is stretched by the fact that the plane-wave solution satisfies the equations of motion at every order without the need for perturbative corrections. As an application of the theory, results for 20 Ne, 22 Ne and 22 Mg are presented that yield energy spectra, B(E2) values, and matter radii in good agreement with experimentally measured results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

J eff = $\frac{3}{2}$ metallic phase and unconventional superconductivity in GaTa 4 Se 8

By means of density functional theory plus dynamical mean-field theory (DFT+DMFT) calculations and resonant inelastic x-ray scattering (RIXS) experiments, we investigate the high-pressure phases of the spin-orbit-coupled Jeff = 3/2 insulator GaTa 4 Se 8 . Its metallic phase, derived from the Mott state by applying pressure, is found to carry J eff = 3/2 moments. The characteristic excitation peak in the RIXS spectrum maintains its destructive quantum interference of J eff at the Ta L2-edge up to 10.4 GPa. Our exact diagonalization based DFT+DMFT calculations including spin-orbit coupling also reveal that the J eff = 3/2 character can be clearly identified under high pressure. These results establish the intriguing nature of the correlated metallic magnetic phase, which represents the first confirmed example of J eff =3/2 moments residing in a metal. They also indicate that the pressure-induced superconductivity is likely unconventional and influenced by these J eff = 3/2 moments. Based on a self-energy analysis, we furthermore propose the possibility of doping-induced superconductivity related to a spin-freezing crossover.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Control-Theory-Informed Feature Selection for Detecting Malicious Tampering in Additive Layer Manufacturing Processes

Additive layer manufacturing (ALM) is rapidly becoming an appealing solution to the low-volume manufacturing of metal, polymer, or composite parts. However, ALM’s reliance on digital part specifications, microcontrollers, and modern networking makes these devices vulnerable to malicious tampering by cyber attackers, which can negatively affect part performance and even result in catastrophic failure. We present a hybrid analytic approach to feature discovery using control theoretic techniques and linear modelling on input-output data collected from a representative controller system. Employing this approach, we design, train, and test an anomaly detection system. The preliminary results show that the proposed approach effectively discovers useful input-output relationships for anomaly detection in a simulated ALM process. Application to larger and more complex systems are discussed.

Dawson, Joel↗
DOE OSTI2025

Virtual element approximations of the time-fractional nonlinear convection-diffusion equation on polygonal meshes

We extend the Virtual Element Method to a two-dimensional unsteady nonlinear convection-diffusion equation characterized by a fractional-order derivative with respect to the time variable. Our methodology is based on three fundamental technical components: a fractional version of the Grunwald-Letnikov approximation, discrete maximal regularity, and the regularity theory associated with non-linearity. We prove the method's well-posedness, i.e., the approximate solution's existence and uniqueness to the time-fractional convection-diffusion equation with a Lipschitz nonlinear source term. The fully discrete scheme inherently maintains stability and consistency by leveraging the discrete maximal regularity and the energy projection operator. The convergence in the L 2 -norm and H 1 -norm to various mesh configurations is validated by numerical results, underlining the practical effectiveness of the proposed method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Density Functional Theory Guided Investigation of Ligand‐Induced Neptunyl‐Neptunyl Interactions

Abstract Actinyl‐actinyl interactions are particularly prevalent for the pentavalent neptunyl cation (Np(V)O 2 ) + where these interactions appear either as a T ‐ or D ‐shape (diamond‐shape). T ‐shaped interactions have been previously identified in high concentration Np(V) solutions containing simple anions (NO 3 − , ClO 4 − , Cl − ) whereas D ‐shaped have only been isolated in the solid‐state in the presence of carboxylate ligands. In this study, Density Functional Theory (DFT) calculations were paired with Raman spectroscopy to evaluate the formation of D ‐shaped interactions in the presence of aliphatic (R=H (formate), CH 3 (acetate), CH 2 CH 3 (propionate)) and aromatic (R=C 6 H 5 (benzoate), C 6 H 4 OH (4‐hydroxybenzoate), C 5 H 4 N (isonicotinate)) carboxylate ligands. DFT studies indicate that the Δ G to form hydrated T ‐ and D ‐shaped forms are not spontaneous but become so with the addition of the carboxylate ligands. Raman spectra of the Np(V) carboxylate solutions contained vibrational modes associated with the D ‐shaped interactions, but spectral changes observed over time indicate a dynamic system. Crystallization experiments from the Np(V) carboxylate systems confirmed the presence of D‐ shaped dimers for the aromatic carboxylates, suggesting that the choice of the anion in solution favors actinyl‐actinyl interactions even at low concentrations (≤20 mM) of Np(V).

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

Global optimization of harmonic oscillator basis in covariant density functional theory

The present investigation focuses on the improvement of the accuracy of the description of binding energies within moderately sized fermionic basis. Using the solutions corresponding to infinite fermionic basis it was shown that in the case of meson exchange (ME) covariant energy density functionals (CEDFs) the global accuracy of the description of binding energies in the finite $N_F$ = 16 - 20 bases can be drastically (by a factor ranging from ~3 up to ~9 dependent on the functional and $N_F$) improved by a global optimization of oscillator frequency of the basis. This is a consequence of the unique feature of the ME functionals in which with increasing fermionic basis size fermionic and mesonic energies approach the exact (infinite basis) solution from above and below, respectively. As a consequence, an optimal oscillator frequency $\hbar\omega_0$ of the basis can be defined which provides an accurate reproduction of exact total binding energies by the ones calculated in truncated basis. This leads to a very high accuracy of the calculations in moderately sized $N_F=20$ basis when mass dependent oscillator frequency is used: global rms differences $\delta B_{rms}$ between the binding energies calculated in infinite and truncated bases are only 0.025 MeV and 0.031 MeV for the NL5(Z) and DD-MEZ functionals, respectively. Optimized values of the oscillator frequency $\hbar\omega_0$ are provided for three major classes of CEDFs, i.e. for density dependent meson exchange functionals, nonlinear meson exchange ones and point coupling functionals.

Binding energy & masses↗
DOE OSTI2023

A map of single-phase high-entropy alloys

High-entropy alloys have exhibited unusual materials properties. The stability of equimolar single-phase solid solution of five or more elements is supposedly rare and identifying the existence of such alloys has been challenging because of the vast chemical space of possible combinations. Herein, based on high-throughput density-functional theory calculations, we construct a chemical map of single-phase equimolar high-entropy alloys by investigating over 658,000 equimolar quinary alloys through a binary regular solid-solution model. We identify 30,201 potential single-phase equimolar alloys (5% of the possible combinations) forming mainly in body-centered cubic structures. We unveil the chemistries that are likely to form high-entropy alloys, and identify the complex interplay among mixing enthalpy, intermetallics formation, and melting point that drives the formation of these solid solutions. We demonstrate the power of our method by predicting the existence of two new high-entropy alloys, i.e. the body-centered cubic AlCoMnNiV and the face-centered cubic CoFeMnNiZn, which are successfully synthesized.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A chemical kinetics theory for interpreting the non-monotonic temperature dependence of enzymatic reactions

Abstract. One notable observation of enzymatic chemical reactions is that, for a given abundance of enzymes and substrates, temperature increases cause reaction rates to first increase consistent with the Arrhenius relationship, then plateau, and finally fall off quickly to zero at high temperatures. While many mathematical functions have been used to describe this pattern, we here propose a chemical kinetics theory which successfully replicates this observation and provides insights into the processes responsible for these dynamics. The chemical kinetics theory combines the law of mass action, von Smoluchowski's diffusion-limited chemical reaction theory, and Eyring's transition state theory. This new theory reveals that the thermally reversible enzyme denaturation ensured by the ceaseless thermal motion of molecules and ions in an enzyme solution explains the plateau and subsequent decrease in chemical reaction rates with increasing temperature. The temperature-dependent affinity parameter (K) that relates enzymes and substrates through their binding also affects the shape of the emergent temperature response. We demonstrate that with an increase in substrate availability, K shifts the optimal temperature, where reaction rates plateau, towards higher values. Further, we show that the chemical kinetics theory accurately represents 12 sets of published enzyme assay data and includes the popular mechanistic model by Ratkowsky et al. (2005) as a special case. Given its good performance and solid theoretical underpinning, we believe this new theory will facilitate the construction of more mechanistic-based environmental biogeochemical models.

Tang, Jinyun↗
DOE OSTI2021

Generalising G 2 geometry: involutivity, moment maps and moduli

We analyse the geometry of generic Minkowski N = 1, D = 4 flux compactifications in string theory, the default backgrounds for string model building. In M-theory they are the natural string theoretic extensions of G2 holonomy manifolds. In type II theories, they extend the notion of Calabi-Yau geometry and include the class of flux backgrounds based on generalised complex structures first considered by Graña et al. (GMPT). Using E 7(7) × $\mathbb{R}$ + generalised geometry we show that these compactifications are characterised by an SU(7) ⊂ E 7(7) structure defining an involutive subbundle of the generalised tangent space, and with a vanishing moment map, corresponding to the action of the diffeomorphism and gauge symmetries of the theory. The Kähler potential on the space of structures defines a natural extension of Hitchin’s G 2 functional. Using this framework we are able to count, for the first time, the massless scalar moduli of GMPT solutions in terms of generalised geometry cohomology groups. It also provides an intriguing new perspective on the existence of G 2 manifolds, suggesting possible connections to Geometrical Invariant Theory and stability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Computing the Central Charge of the 3D Ising CFT Using Quantum Finite Elements

The 3D Ising conformal field theory (CFT) describes different physical systems, such as uniaxial magnets or fluids, at their critical points. In absence of an analytical solution for the 3D Ising model, the scaling dimensions and operator product expansion (OPE) coefficients characterizing this CFT must be determined numerically. The currently most-cited values for these quantities have been obtained from the conformal bootstrap, while lattice calculations have so far only produced reliable results for the scaling dimensions involved in calculating the critical exponents. Using Quantum Finite Elements to investigate critical \(\phi^4\)-theory on \(\mathbb{R}\times\mathbb{S}^2\), we show in this work that it is possible to extract scaling dimensions and OPE coefficients of the 3D Ising CFT by fitting the lattice four-point function with expectations from the operator product expansion for the radially quantized CFT and extrapolating to the continuum limit. This way, we have for the first time been able to use Monte Carlo simulations to compute the central charge of the theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Parity solution to the strong CP problem and a unified framework for inflation, baryogenesis, and dark matter

It has been known for some time that asymptotic parity invariance of weak interactions can provide a solution to the strong CP problem without the need for the axion. Left-right symmetric theories which employ a minimal Higgs sector consisting of a left-handed and a right-handed doublet is an example of such a theory wherein all fermion masses arise through a generalized seesaw mechanism. In this paper we present a way to understand the origin of matter-antimatter asymmetry as well as the dark matter content of the universe in these theories using the Affleck-Dine (AD) leptogenesis mechanism and inflaton decay, respectively. Three gauge singlet fermions are needed for this purpose, two of which help to implement the Dirac seesaw for neutrino masses while the third one becomes the non-thermal dark matter candidate. A soft lepton number breaking term involving the AD scalar field is used to generate lepton asymmetry which suffers no wash-out effects and maintains the Dirac nature of neutrinos. This framework thus provides a unified description of many of the unresolved puzzles of the standard model that require new physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Universal entropy and hawking radiation of near-extremal AdS4 black holes

We compute the Bekenstein-Hawking entropy of near-extremal asymptotically AdS 4 electrically charged rotating black holes using three different methods: (i) from the gravity solution, (ii) from the near-horizon Kerr/CFT correspondence and (iii) from the boundary conformal field theory. The results from these three different approaches match exactly, giving us a unique and universal expression for the entropy and the microstate counting of near-extremal AdS black holes via the AdS/CFT correspondence. In the second method, we extend the Kerr/CFT correspondence to the near-extremal case to compute the left and right central charges. We also use hidden conformal symmetry of the near-horizon geometry to compute the Frolov-Thorne temperatures. From the results of the near-extremal AdS 4 black hole entropy, we provide a microscopic foundation for Hawking radiation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Soft de Sitter Effective Theory

Calculating the quantum evolution of a de Sitter universe on superhorizon scales is notoriously difficult. To address this challenge, we introduce the Soft de Sitter Effective Theory (SdSET). This framework holds for superhorizon modes whose comoving momentum is far below the UV scale, which is set by the inverse comoving horizon. The SdSET is formulated using the same approach that yields the Heavy Quark Effective Theory. The degrees of freedom that capture the long wavelength dynamics are identified with the growing and decaying solutions to the equations of motion. The operator expansion is organized using a power counting scheme, and loops can be regulated while respecting the low energy symmetries. For massive quantum fields in a fixed de Sitter background, power counting implies that all interactions beyond the horizon are irrelevant. Alternatively, if the fields are very light, the leading interactions are at most marginal, and resumming the associated logarithms using (dynamical) renormalization group techniques yields the evolution equation for canonical stochastic inflation. The SdSET is also applicable to models where gravity is dynamical, including inflation. In this case, diffeomorphism invariance ensures that all interactions are irrelevant, trivially implying the all-orders conservation of adiabatic density fluctuations and gravitational waves. We briefly touch on the application to slow-roll eternal inflation by identifying novel relevant operators. This work serves to demystify many aspects of perturbation theory outside the horizon, and has a variety of applications to problems of cosmological interest.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Modeling gas release from a Bingham plastic slurry and deconvoluting measured data

Physics-based models for predicting the off-gassing characteristics of agitated slurries are presented. This approach decomposes the system into two separate but coupled slurry and headspace models. The physics driving bubble transport through the slurry and gas mixing within the headspace are discussed. An analytical expression for predicting the time evolution of the headspace concentration using first principles theory is also presented. Predictions from the numerical models, as well as expectations from the analytical solution, both agree with measured off gassing data for two different experimental operating conditions. After benchmarking the numerical model predictions against experimental data, the model was used to make predictions for gas release from a full-scale vessel and to perform sensitivity analyses to examine the sensitivity of gas release to parameters such as yield stress, consistency index, slurry density, bubble size, bubble concentration, and impeller speed.

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

A method for characterization of multiple dynamic constitutive parameters of FRCs

We propose a method to measure multiple dynamic material constitutive parameters of unidirectional fiber reinforced composites (FRCs) in a single experiment. Dynamic short-beam shear (DSBS) experiments were performed on a modified Kolsky compression bar, with integration of high-speed imaging and digital image correlation (DIC). The unidirectional FRCs investigated were S-2 glass fiber reinforced matrix of TGDDM-Jeffamine® D230 with monoamine functionalized partially reacted substructures (mPRS) and commercially available SC-15. Analytical solutions of normal and shear strains of a composite beam were derived based on Timoshenko beam theory, assuming material to be transversely isotropic and have different moduli in tension and compression in each principle material orientation. Tensile and compressive moduli were inversely computed through monitoring normal strain slope when specimen was constantly loaded at a speed of ~7.3 m/s within a small deflection. Non-linear shear stress-strain behavior of the composite was described via Ramberg–Osgood equation. Finite element (FE) analysis was conducted in ABAQUS, simultaneously defining via user subroutine UMAT the transverse isotropy of material, bi-modulus constitutive model, and non-linear shear stress-strain relation. The method proposed in this work was validated by comparing strain distributions computed by FE model and DIC measurements. Comparing with traditional dynamic tensile, compressive, and shear experiments on FRCs, this method significantly simplifies the specimen preparation and design of complicated gripping fixtures for multiple experiments. Here, systematic errors resulting from variations of specimen geometry and dimension, loading direction, and instrumentation are reduced, thereby providing compatible data for numerical studies on impact behavior of composites.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Room temperature sodium diffusion in sodium iron phosphate investigated with neutron scattering

Sodium iron phosphate, Na x FePO 4 , is currently one of the most promising cathode materials for sodium ion batteries. Here we employ a systematic approach using three neutron scattering techniques to investigate temperature dependent properties with a focus on the charge transfer processes. Two distinct Na + diffusive motions are observed with quasielastic neutron scattering, including room temperature motion. These processes correspond to rapid double jumps through extrinsic vacancy sites of ∼6 Å along with the more fundamental nearest neighbor jumps of 3.1 Å. Experimental observations of room temperature ionic diffusion are highly desirable, as the diffusion pathways are naturally well suited for efficient transport. Inelastic neutron scattering was used to probe the vibrational density of states that assist in the ionic and electronic (polaronic) transport processes, which are further supported with density functional theory. Neutron diffraction examines the disordering of the sodium sublattice upon transition to the solid-solution phase, which is believed to enhance the overall kinetics during sodium insertion/extraction. Finally, 57 Fe Mössbauer spectroscopy measurements probes polaron hopping and the local structural evolution that modifies the electron density surrounding the redox active iron sites. This study sheds light on the temperature induced atomic scale dynamics that occur in olivine Na x FePO 4 .

Novak, Eric [Swarthmore College, Swarthmore, PA (U↗