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

Dimer piling problems and interacting field theory

The dimer tiling problem asks in how many ways can the edges of a graph be covered by dimers so that each site is covered once. In the special case of a planar graph, this problem has a solution in terms of a free fermionic field theory. We rediscover and explore an expression for the number of coverings of an arbitrary graph by arbitrary objects in terms of an interacting fermionic field theory first proposed by Samuel. Generalizations of the dimer tiling problem, which we call “dimer piling problems,” demand that each site be covered N times by indistinguishable dimers. Our field theory provides a solution of these problems in the large- N limit. We give a similar path integral representation for certain lattice coloring problems. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2021

Generating rotating spacetime in Ricci-based gravity: naked singularity as a black hole mimicker

Motivated by the lack of rotating solutions sourced by matter in General Relativity as well as in modified gravity theories, in this work we extend a recently discovered exact rotating solution of the minimal Einstein-scalar theory to its counterpart in Eddington-inspired Born-Infeld gravity coupled to a Born-Infeld scalar field. This is accomplished with the implementation of a well-developed mapping between solutions of Ricci-Based Palatini theories of gravity and General Relativity. The new solution is parametrized by the scalar charge and the Born-Infeld coupling constant apart from the mass and spin of the compact object. Compared to the spacetime prior to the mapping, we find that the high-energy modifications at the Born-Infeld scale are able to suppress but not remove the curvature divergence of the original naked null singularity. Depending on the sign of the Born-Infeld coupling constant, these modifications may even give rise to an additional timelike singularity exterior to the null one. In spite of that, both of the naked singularities before and after the mapping are capable of casting shadows, and as a consequence of the mapping relation, their shadows turn out to be identical as seen by a distant observer on the equatorial plane. Even though the scalar field induces a peculiar oblateness to the appearance of the shadow with its left and right endpoints held fixed, the closedness condition for the shadow contour sets a small upper bound on the absolute value of the scalar charge, which leads to observational features of the shadow closely resembling those of a Kerr black hole.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Ab initio investigation of elastic properties of dilute Cu alloys for high-gradient accelerating structures

The elastic properties and solid solution strengthening of several Cu alloys that may be promising for high gradient applications were studied using ab initio calculations based on density functional theory. Specifically, the temperature and solute concentration dependence of elastic properties including bulk and shear moduli, thermal expansion coefficient, and heat capacity of these alloys were calculated. We found that the bulk and shear moduli are linearly dependent on the solute concentration in the dilute regime and change significantly for different Cu alloys. The thermal expansion coefficient and heat capacity are comparatively less sensitive to solute concentrations. Solid solution strengthening is observed to be linearly correlated with the change of lattice parameter of Cu. This work provides an in-depth understanding of the elastic properties of Cu alloys, which is key for predicting their performance at high gradient in normal conducting accelerating structures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Logarithmic potential with super-super-exponential kink profiles and tails

In this work, we consider a novel one dimensional model of a logarithmic potential which has super-super-exponential kink profiles as well as kink tails. We provide analytic kink solutions of the model—it has 3 kinks, 3 mirror kinks and the corresponding antikinks. While some of the kink tails are super-super-exponential, some others are super-exponential whereas the remaining ones are exponential. The linear stability analysis reveals that there is a gap between the zero mode and the onset of continuum. Finally, we compare this potential and its kink solutions with those of very high order field theories harboring seven degenerate minima and their attendant kink solutions, specifically φ 14 , φ 16 and φ 18 .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Multiparametric Programming in Process Systems Engineering: Recent Developments and Path Forward

The inevitable presence of uncertain parameters in critical applications of process optimization can lead to undesirable or infeasible solutions. For this reason, optimization under parametric uncertainty was, and continues to be a core area of research within Process Systems Engineering. Multiparametric programming is a strategy that offers a holistic perspective for the solution of this class of mathematical programming problems. Specifically, multiparametric programming theory enables the derivation of the optimal solution as a function of the uncertain parameters, explicitly revealing the impact of uncertainty in optimal decision-making. By taking advantage of such a relationship, new breakthroughs in the solution of challenging formulations with uncertainty have been created. Apart from that, researchers have utilized multiparametric programming techniques to solve deterministic classes of problems, by treating specific elements of the optimization program as uncertain parameters. In the past years, there has been a significant number of publications in the literature involving multiparametric programming. The present review article covers recent theoretical, algorithmic, and application developments in multiparametric programming. Additionally, several areas for potential contributions in this field are discussed, highlighting the benefits of multiparametric programming in future research efforts.

Pappas, Iosif↗
DOE OSTI2021

Electronic Structure Theory Calculations Using Modern Architectures: KNL vs Haswell

The time to solution and parallel efficiency of several commonly used electronic structure methods (Hartree–Fock, density functional theory, second order perturbation theory, resolution of the identity second order perturbation theory, coupled cluster) are evaluated on both the Intel Xeon Haswell and the Intel Xeon Phi Knights Landing (KNL) architectures. The Haswell completes the benchmark calculations with a faster time to solution than the KNL for all molecules and methods tested. While the Haswell exhibits an average speedup of at least 3.5 relative to the KNL for all nonthreaded computations, the KNL has a better parallel efficiency than the Haswell with increasing core counts. The architectures are further tested using a more computationally costly coupled cluster method on a transition state reaction. The Haswell appears to be the best choice to minimize the time to solution, though for very large systems and high levels of theory that require memory intensive processes the superior memory hierarchy and larger on node memory of the KNL can make it a better choice. These results are used to showcase aspects of novel architectures that will increase efficiency for quantum chemistry applications.

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

Impacts of hydrogen bonding interactions with Np(v/vi)O 2 Cl 4 complexes: vibrational spectroscopy, redox behavior, and computational analysis

The neptunyl (Np(V)O 2 + /Np(VI)O 2 2+ ) cation is the dominant form of 237 Np in acidic aqueous solutions and the stability of the Np(V) and Np(VI) species is driven by the specific chemical constituents present in the system. Hydrogen bonding with the oxo group may impact the stability of these species, but there is limited understanding of how these intermolecular interactions influence the behavior of both solution and solid-state species. In the current study, we systematically evaluate the interactions between the neptunyl tetrachloride species and hydrogen donors in coordination complexes and in the related aqueous solutions. Both Np(V) compounds (N 2 C 4 H 12 ) 2 [Np(V)O 2 Cl 4 ]Cl (Np(V)pipz) and (NOC 4 H 10 ) 3 [Np(V)O 2 Cl 4 ] (Np(V)morph) exhibit directional hydrogen bonding to the neptunyl oxo group while Np(VI) compounds (NC 5 H 6 ) 2 [Np(VI)O 2 Cl 4 ] (Np(VI)pyr) and (NOC 4 H 10 ) 4 [Np(VI)O 2 Cl 4 ]·2Cl (Np(VI)morph) assemble via halogen interactions. Here, the Raman spectra of the solid-state phases indicate the activation of vibrational bands when there is asymmetry of the neptunyl bond, while these spectral features are not observed within the related solution phase spectra. Density functional theory calculations of the Np(V)pipz system suggest that activation of the ν 3 asymmetric stretch and other combination modes lead to additional complexity within the solid-state spectra. Electrochemical analyses of complexes in the solution phases are consistent with the results of the crystallization experiments as the voltammetric potentials of Np(V)/Np(VI) complexes in the presence of protonated heterocycles differ from the potentials of pure Np(V) and may correlate with the hydrogen bonding interactions.

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

Computational modeling and neutron imaging to understand interface shape and solute segregation during the vertical gradient freeze growth of BaBrCl:Eu

In this work we apply continuum models to analyze phase change, heat transfer, fluid flow, solute transport, and segregation in order to understand prior neutron imaging observations of the vertical gradient freeze growth of Eu-doped BaBrCl. The models provide a rigorous framework in which to understand the mechanisms that are responsible for the complicated evolution of interface shape and dopant distribution in the growth experiment. We explain how a transition in the solid/liquid interface shape from concave to convex is driven by changes in radial heat transfer caused by furnace design. We also provide a mechanistic explanation of how dynamic growth conditions and changes of the flow structure in the melt result in complicated segregation patterns in this system. A growth pause caused by controller lock-up is shown to result in a band of solute depletion in accordance with classical theory. However, changing flow patterns during growth result in a non-monotonic axial distribution of solute that cannot be explained by simple application of classical segregation models. We assert that the approach presented here, namely the use of rigorous models in conjunction advanced diagnostics, such as neutron imaging, provides an exciting path forward for process optimization and control, accelerating the incremental advances that have, in the past, typically relied on empiricism, experience, and intuition.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Burst Pressure Solutions of Thin and Thick-Walled Cylindrical Vessels

Pressure vessels (PVs) are widely used in the energy industry. Accurate burst pressure is critical to structural design and safe operation for both thin and thick-walled PVs. The traditional strength theories utilized a single-parameter material property, such as the yield stress or the ultimate tensile stress (UTS) to develop failure models for determining the yield or ultimate pressure carrying capacity in the PV design. The UTS-based Barlow formula is a typical burst pressure model developed from the Tresca strength theory that provides the basis for developing regulation rules and failure models for different industry design codes, such as ASME BPVC, ASME B31.3, and ASME B31G, among others. In order to reduce the conservatism of the Tresca strength model, ASME BPVC recently adapted failure models developed from the von Mises strength theory for the PV design and analysis. It has been commonly accepted that the burst pressure of pipelines depends on the UTS and strain hardening exponent, n, of the pipeline steel. An average shear stress yield theory was thus developed, and the Zhu-Leis solution of burst pressure was obtained as a function of UTS and n for thin-walled line pipes. Experiments showed that the Zhu-Leis solution provides an accurate, reliable prediction of burst pressure for defect-free thin-walled pipes. In order to extend the Zhu–Leis solution to thick-walled cylindrical PVs, this paper defined three new flow stresses, modified the traditional strength theories, and obtained three new burst pressure solutions that are valid for both thin and thick-walled cylindrical vessels. The proposed flow stresses are able to describe the tensile strength and the plastic flow response of PVs for a strain hardening steel. The associated strength theories were then developed in terms of the Tresca, von Mises, and Zhu-Leis yield criteria. From these new strength theories, three burst pressure solutions were obtained for thick-walled cylinders, where the von Mises solution is an upper bound prediction, the Tresca solution is a lower bound prediction, and the Zhu-Leis solution is an intermediate prediction of burst pressure for thick-walled cylinders. Lastly, the proposed burst pressure solutions were evaluated and validated by two large datasets of full-scale burst tests for thick-walled tubes and for thin-walled pipes.

42 ENGINEERING↗
DOE OSTI2024

A Unified Description of Salt Effects on the Liquid–Liquid Phase Separation of Proteins

Protein aggregation via liquid-liquid phase separation (LLPS) is ubiquitous in nature and is intimately connected to many human diseases. Although it is widely known that the addition of salt has crucial impacts on the LLPS of proteins, full understanding of the salt effects remains an outstanding challenge. Here, we develop a molecular theory that systematically incorporates the self-consistent field theory for charged macromolecules into the solution thermodynamics. The electrostatic interaction, hydrophobicity, ion solvation, and translational entropy are included in a unified framework. Our theory fully captures the long-standing puzzles of the nonmonotonic salt concentration dependence and the specific ion effect. We find that proteins show salting-out at low salt concentrations due to ionic screening. The solubility follows the inverse Hofmeister series. In the high salt concentration regime, protein continues salting-out for small ions but turns to salting-in for larger ions, accompanied by the reversal of the Hofmeister series. We reveal that the solubility at high salt concentrations is determined by the competition between the solvation energy and translational entropy of the ion. Furthermore, we derive an analytical criterion for determining the boundary between the salting-in and salting-out regimes, which is in good agreement with experimental results for various proteins and salt ions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

NSTX-U theory, modeling and analysis results

Here, the mission of the low aspect ratio spherical tokamak NSTX-U is to advance the physics basis and technical solutions required for optimizing the configuration of next-step steady-state tokamak fusion devices. NSTX-U will ultimately operate at up to 2 MA of plasma current and 1 T toroidal field on axis for 5 s, and has available up to 15 MW of neutral beam injection power at different tangency radii and 6 MW of high harmonic fast wave heating. With these capabilities NSTX-U will develop the physics understanding and control tools to ramp-up and sustain high performance fully non-inductive plasmas with large bootstrap fraction and enhanced confinement enabled via the low aspect ratio, high beta configuration. With its unique capabilities, NSTX-U research also supports ITER and other critical fusion development needs. Super-Alfvénic ions in beam-heated NSTX-U plasmas access energetic particle (EP) parameter space that is relevant for both α-heated conventional and low aspect ratio burning plasmas. NSTX-U can also generate very large target heat fluxes to test conventional and innovative plasma exhaust and plasma facing component solutions. This paper summarizes recent analysis, theory and modelling progress to advance the tokamak physics basis in the areas of macrostability and 3D fields, EP stability and fast ion transport, thermal transport and pedestal structure, boundary and plasma material interaction, RF heating, scenario optimization and real-time control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A self-consistent electrostatic electrified shallow water model for charged liquid surface dynamics

The dynamics of an electrified liquid surface are investigated using a shallow water model that is self-consistently coupled with an electrostatic solver. To account for the spatiotemporal variation of a curved liquid surface, the electric field on curved surfaces is calculated by solving a two-dimensional electrostatic equation using the weighted least squares (WLSQ) interpolation method. First, the WLSQ implementation is verified with analytical theory obtained from a Laplace solution assuming a half-cylinder liquid surface profile. Then, the coupled electrified shallow water model is used to study the instability of liquid surface perturbation as a function of the potential drop between an electrode and conducting liquid, surface tension, and gravity. We present a linear dispersion theory of liquid surface instability for long-wavelength perturbations, including the effects of liquid viscosity. Furthermore, the effects of multiple sinusoidal surface perturbations on the electrified liquid surface instability are investigated.

Capillary waves↗
DOE OSTI2022

Polar liquids at charged interfaces: A dipolar shell theory

The structure of polar liquids and electrolytic solutions, such as water and aqueous electrolytes, at interfaces underlies numerous phenomena in physics, chemistry, biology, and engineering. In this work, we develop a continuum theory that captures the essential features of dielectric screening by polar liquids at charged interfaces, including decaying spatial oscillations in charge and mass, starting from the molecular properties of the solvent. The theory predicts an anisotropic dielectric tensor of interfacial polar liquids previously studied in molecular dynamics simulations. We explore the effect of the interfacial polar liquid properties on the capacitance of the electrode/electrolyte interface and on hydration forces between two plane-parallel polarized surfaces. In the linear response approximation, we obtain simple formulas for the characteristic decay lengths of molecular and ionic profiles at the interface.

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

Generalized models for inflationary preheating: Oscillations and symmetries

The paradigm of the inflationary universe provides a possible explanation for several observed cosmological properties. In order for such solutions to be successful, the universe must convert the energy stored in the inflaton potential into standard model particles through a process known as reheating. In this paper, we reconsider the reheating process for the case where the inflaton potential respects an approximate (but spontaneously broken) conformal symmetry during the reheating epoch. After reviewing the Effective Field Theory of Reheating, we present solutions for the nonlinear oscillations of the inflaton field, derive the corresponding Hill’s equation for the coupled reheating field, and determine the stability diagram for parametric resonance. For this class of models —the simplest realization being a scalar field with a quartic term—the expansion of the universe drives the coupled field toward a more unstable part of parameter space, in contrast to the standard case. We also generalize this class of models to include quadratic breaking terms in the potential during the reheating epoch and address the process of stability in that universality class of models.

Barrowes, Leia [University of Michigan, Ann Arbor]↗
DOE OSTI2024

A Gaussian field approach to the solvation of spherical ions in electrolyte solutions

In this paper, the electrostatic response of an electrolyte solution to a spherical ion is studied with a Gaussian field theory. In order to capture the ionic correlation effect in concentrated solutions, the bulk dielectric response function is described by a two-Yukawa response function. The modified response function of the solution is solved analytically in the spherical geometry, from which the induced charge density and the electrostatic energy are also derived analytically. Comparisons with results for small ions in electrolyte solutions from the hyper-netted chain theory demonstrate the validity of the Gaussian field theory.

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

Dissertation Proposal: Innovations in Pulsed Power and Plasma Science Theory Simulation and Experiment

Pulsed power and plasma physics are topics of great study at both Sandia National Laboratories (SNL or Sandia) and the University of New Mexico (UNM). The goal of this research is to further knowledge and understanding of these fields using the resources of both SNL and UNM in three ways. The first way is through the comprehension, application, and testing of theory. Reading and analytically deriving theoretical solutions of problems both real-world and simplified will allow for a fresh perspective and the furthering of the theory. One such theory is Ottinger's generalized theory for voltage measurement in magnetically insulated transmission lines (MITLs). By working through the math, a deeper understanding of the theory is gained from which one may add more physically accurate and/or more detailed physics into the theory. Additionally, understanding the theory lays a good foundation from which one can analyze, test, and compare results to the theory in the following two ways that will advance the fields of pulsed power and plasma physics. The second way is through the modeling and simulation of real-world and simplified problems that utilize and test the afore mentioned theories. Theory can be applied to a simulation domain by using the unstructured time-domain electromagnetic (UTDEM) codes EMPHASIS and EMPIRE as well as the physical modeling software CUBIT, all of which were developed at SNL. Problems such as the modeling and design of the extended MITL on HERMES III, the understanding of space-charge-limited emission from vacuum cathodes, and the interaction between a relativistic electron beam and an ideal gas can all be modeled, simulated, and analyzed with this set of codes. Here the advantage is three-fold. Firstly, theory that describes our understanding of these problems can be put to the test and advanced through iterative simulation and analysis. Secondly, the understanding of these problems will have a positive impact on national security through the advancement of the technological capability of the United States of America. Thirdly, and not unrelated to the prior advantage, is the validation and verification of EMPIRE and EMPHASIS. This segues into the third way, which is through experiment and the comparison of experiment to simulated and theoretical results. Performing experimental comparisons completes the scientific method and grounds all of the work in reality. Being able to physically test theory and simulation is necessary for any real conclusions to be drawn. Another advantage for carrying out experimental work is to advance the physical testing capabilities of SNL. Several systems will be developed and tested through the course of this work that positively impact technological advancement of Sandia National Labs. Lastly, all of the above work will converge to yield a well-rounded perspective that ties the three categories of research together.

42 ENGINEERING↗
DOE OSTI2023

Traversable wormholes in four dimensions

Here, we present a wormhole solution in four dimensions. It is a solution of an Einstein Maxwell theory plus charged massless fermions. The fermions give rise to a negative Casimir-like energy, which makes the wormhole possible. It is a long wormhole that does not lead to causality violations in the ambient space. It can be viewed as a pair of entangled near extremal black holes with an interaction term generated by the exchange of fermion fields. The solution can be embedded in the Standard Model by making its overall size small compared to the electroweak scale.

79 ASTRONOMY AND ASTROPHYSICS↗