Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SINGULARITY”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

NONSINGULAR UNIVERSES IN GAUSS–BONNET GRAVITY’S RAINBOW

In this paper, we study the rainbow deformation of Friedmann-Robertson-Walker (FRW) cosmology in both Einstein gravity and Gauss–Bonnet (GB) gravity. We demonstrate that the singularity in FRW cosmology can be removed because of the rainbow deformation of the FRW metric. We obtain the general constraints required for FRW cosmology to be free of singularities. We observe that the inclusion of GB gravity can significantly change the constraints required to obtain nonsingular universes. We use rainbow functions motivated by the hard spectra of gamma-ray bursts to deform FRW cosmology and explicitly demonstrate that such a deformation removes the singularity in FRW cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Automatically adaptive stabilized finite elements and continuation analysis for compaction banding in geomaterials

Under compressive creep, viscoplastic solids experiencing internal mass transfer processes can accommodate singular cnoidal wave solutions as material instabilities at the stationary wave limit. These instabilities appear when the loading rate is significantly faster than the material's capacity to diffusive internal perturbations, leading to localized failure features (e.g., cracks and compaction bands). These cnoidal waves, generally found in fluids, have strong nonlinearities that produce periodic patterns. Due to the singular nature of the solutions, the applicability of the theory is currently limited. Additionally, practical simulation tools require proper regularization to overcome the challenges that singularity induces. We focus on the numerical treatment of the governing equation using a nonlinear approach building on a recent adaptive stabilized finite element method. This automatic refinement method provides an error estimate that drives mesh adaptivity, a crucial feature for the problem at hand. We compare the performance of this adaptive strategy against analytical and standard finite element solutions. We then investigate the sensitivity of the diffusivity ratio, the parameter controlling the process, and identify multiple possible solutions with several stress peaks. We also show the evolution of the spacing between peaks for all solutions as a function of that parameter.

42 ENGINEERING↗
DOE OSTI2019

Balancing Performance and Portability with Containers in HPC: An OpenSHMEM Example

There is a growing interest in using Linux containers to streamline software development and application deployment. A container enables the user to bundle the salient elements of the software stack from an application’s perspective. In this paper, we discuss initial experiences in using the Open MPI implementation of OpenSHMEM with containers on HPC resources. We provide a brief overview of two container runtimes, Docker & Singularity, highlighting elements that are of interest for HPC users. The Docker platform offers a rich set of services that are widely used in enterprise environments, whereas Singularity is an emerging container runtime that is specifically written for use on HPC systems. We describe our procedure for container assembly and deployment that strives to maintain the portability of the container-based application. We show performance results for the Graph500 benchmark running along the typical continuum of development testbed up to full production supercomputer (ORNL’s Titan). The results show consistent performance between the native and Singularity (container) tests. The results also showed an unexplained drop in performance when using the Cray Gemini network with Open MPI’s OpenSHMEM, which was unrelated to the container usage.

Naughton III, Thomas↗
DOE OSTI2021

Constraints on a massive double-copy and applications to massive gravity

We propose and study a BCJ double-copy of massive particles, showing that it is equivalent to a KLT formula with a kernel given by the inverse of a matrix of massive bi-adjoint scalar amplitudes. For models with a uniform non-zero mass spectrum we demonstrate that the resulting double-copy factors on physical poles and that up to at least 5-particle scattering, color-kinematics duality satisfying numerators always exist. For the scattering of 5 or more particles, the procedure generically introduces spurious singularities that must be cancelled by imposing additional constraints. When massive particles are present, color-kinematics duality is not enough to guarantee a physical double-copy. As an example, we apply the formalism to massive Yang-Mills and show that up to 4-particle scattering the double-copy construction generates physical amplitudes of a model of dRGT massive gravity coupled to a dilaton and a two-form with dilaton parity violating couplings. We show that the spurious singularities in the 5-particle double-copy do not cancel in this example, and the construction fails to generate physically sensible amplitudes. We conjecture sufficient constraints on the mass spectrum, which in addition to massive BCJ relations, guarantee the absence of spurious singularities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

On information loss in AdS 3 /CFT 2

We discuss information loss from black hole physics in AdS 3 , focusing on two sharp signatures infecting CFT 2 correlators at large central charge c: ‘forbidden singularities’ arising from Euclidean-time periodicity due to the effective Hawking temperature, and late-time exponential decay in the Lorentzian region. We study an infinite class of examples where forbidden singularities can be resolved by non-perturbative effects at finite c, and we show that the resolution has certain universal features that also apply in the general case. Analytically continuing to the Lorentzian regime, we find that the non-perturbative effects that resolve forbidden singularities qualitatively change the behavior of correlators at times t ~S BH , the black hole entropy. This may resolve the exponential decay of correlators at late times in black hole backgrounds. By Borel resumming the 1/c expansion of exact examples, we explicitly identify ‘information-restoring’ effects from heavy states that should correspond to classical solutions in AdS 3 . Lastly, our results suggest a line of inquiry towards a more precise formulation of the gravitational path integral in AdS 3 .

1/N expansion↗
DOE OSTI2018

Characteristically Closed Domains for First Order Strictly Hyperbolic Systems in the Plane

We consider a first order strictly hyperbolic system of n equations with constant coefficients in a bounded domain. It is assumed that the domain is strictly convex relative to characteristics, so that the projection along each characteristic is an involution having two fixed singular points. The natural statement of boundary value problems for such systems requires that singular points go to singular points under such transformations. We present a necessary and sufficient condition for the existence of such domains, called characteristically closed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

New implicitly solvable potential produced by second order shape invariance

The procedure proposed recently by Bougie et al. (2010) to study the general form of shape invariant potentials in one-dimensional Supersymmetric Quantum Mechanics (SUSY QM) is generalized to the case of Higher Order SUSY QM with supercharges of second order in momentum. A new shape invariant potential is constructed by this method. It is singular at the origin, it grows at infinity, and its spectrum depends on the choice of connection conditions in the singular point. The corresponding Schrödinger equation is solved explicitly: the wave functions are constructed analytically, and the energy spectrum is defined implicitly via the transcendental equation which involves Confluent Hypergeometric functions. - Highlights: • New potential with 2nd order irreducible shape invariance was constructed. • The connection conditions at the singularity of potential were obtained. • The explicit expressions for all wave functions were derived. • The implicit equation for the energy spectrum was obtained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

The exact discontinuity of a partial wave along the left-hand cut and the exact N∕D method in non-relativistic scattering

We first deduce the analytical continuation in the complex planes of the initial and final three-momenta of the Lippmann–Schwinger equation in coupled or uncoupled partial-wave amplitudes. This result allows us to deduce a master equation whose solution is the exact discontinuity of the on-shell partial-wave amplitudes along the left-hand cut. This equation is always a linear non-singular integral equation whose solution is fixed exclusively by the knowledge of the potential, applicable to either regular or singular potentials. The capability of calculating exactly this discontinuity allows one to settle the exact N∕D method in two-body non-relativistic scattering for coupled and uncoupled waves. We exemplify this new advance in scattering theory by explicitly checking the agreement between the Lippmann–Schwinger equation with the corresponding solutions of the exact N∕D method for some examples that involve regular and singular potentials, either attractive or repulsive.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Superintegrable systems from block separation of variables and unified derivation of their quadratic algebras

Highlights: • A new method for constructing D-dimensional minimally superintegrable systems. • Superintegrable systems from block coordinate separation of variables. • Unified derivation of quadratic symmetry algebras via gauge transformations. - Abstract: We present a new method for constructing D-dimensional minimally superintegrable systems based on block coordinate separation of variables. We give two new families of superintegrable systems with N (N≤D) singular terms of the partitioned coordinates and involving arbitrary functions. These Hamiltonians generalize the singular oscillator and Coulomb systems. We derive their exact energy spectra via separation of variables. We also obtain the quadratic algebras satisfied by the integrals of motion of these models. We show how the quadratic symmetry algebras can be constructed by novel application of the gauge transformations from those of the non-partitioned cases. We demonstrate that these quadratic algebraic structures display a universal nature to the extent that their forms are independent of the functions in the singular potentials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

A spatially adaptive high-order meshless method for fluid–structure interactions

We present a scheme implementing an a posteriori refinement strategy in the context of a high-order meshless method for problems involving point singularities and fluid–solid interfaces. The generalized moving least squares (GMLS) discretization used in this work has been previously demonstrated to provide high-order compatible discretization of the Stokes and Darcy problems, offering a high-fidelity simulation tool for problems with moving boundaries. The meshless nature of the discretization is particularly attractive for adaptive h -refinement, especially when resolving the near-field aspects of variables and point singularities governing lubrication effects in fluid–structure interactions. We demonstrate that the resulting spatially adaptive GMLS method is able to achieve optimal convergence in the presence of singularities for both the div-grad and Stokes problems. Furthermore, we present a series of simulations for flows of colloid suspensions, in which the refinement strategy efficiently achieved highly accurate solutions, particularly for colloids with complex geometries.

42 ENGINEERING↗
DOE OSTI2025

Non-intrusive reduced-order modeling for dynamical systems with spatially localized features

This work presents a non-intrusive reduced-order modeling framework for dynamical systems with spatially localized features characterized by slow singular value decay. The proposed approach builds upon two existing methodologies for reduced and full-order non-intrusive modeling, namely Operator Inference (OpInf) and sparse Full-Order Model (sFOM) inference. We decompose the domain into two complementary subdomains that exhibit fast and slow singular value decay. The dynamics of the subdomain exhibiting slow singular value decay are learned with sFOM while the dynamics with intrinsically low dimensionality on the complementary subdomain are learned with OpInf. The resulting, coupled OpInf-sFOM formulation leverages the computational efficiency of OpInf and the high resolution of sFOM, and thus enables fast non-intrusive predictions for conditions beyond those sampled in the training data set. A novel regularization technique with a closed-form solution based on the Gershgorin disk theorem is introduced to promote stable sFOM and OpInf models. We also provide a data-driven indicator for subdomain selection and ensure solution smoothness over the interface via a post-processing interpolation step. We evaluate the efficiency of the approach in terms of offline and online speedup through a quantitative, parametric computational cost analysis. We demonstrate the coupled OpInf-sFOM formulation for two test cases: a one-dimensional Burgers’ model for which accurate predictions beyond the span of the training snapshots are presented, and a two-dimensional parametric model for the Pine Island Glacier ice thickness dynamics, for which the OpInf-sFOM model achieves an average prediction error on the order of 1% with an online speedup factor of approximately 8$\times$ compared to the numerical simulation.

42 ENGINEERING↗
DOE OSTI2024

Principal Landau determinants

We reformulate the Landau analysis of Feynman integrals with the aim of advancing the state of the art in modern particle-physics computations. We contribute new algorithms for computing Landau singularities, using tools from polyhedral geometry and symbolic/numerical elimination. Inspired by the work of Gelfand, Kapranov, and Zelevinsky (GKZ) on generalized Euler integrals, we define the principal Landau determinant of a Feynman diagram. We illustrate with a number of examples that this algebraic formalism allows to compute many components of the Landau singular locus. We adapt the GKZ framework by carefully specializing Euler integrals to Feynman integrals. For instance, ultraviolet and infrared singularities are detected as irreducible components of an incidence variety, which project dominantly to the kinematic space. We compute principal Landau determinants for the infinite families of one-loop and banana diagrams with different mass configurations, and for a range of cutting-edge Standard Model processes. Furthermore, our algorithms build on the Julia package this http URL and are implemented in the new open-source package this http URL available at this https URL.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Quantum cosmology of Eddington-Born–Infeld gravity fed by a scalar field: The big rip case

We study the quantum avoidance of the big rip singularity in the Eddington-inspired-Born–Infeld (EiBI) phantom model. Instead of considering a simple phantom dark energy component, which is described by a perfect fluid, we consider a more fundamental degree of freedom corresponding to a phantom scalar field with its corresponding potential, which would lead the classical universe to a big rip singularity. We apply a quantum geometrodynamical approach by performing an appropriate Hamiltonian study including an analysis of the constraints of the system. We then derive the Wheeler–DeWitt (WDW) equation and see whether the solutions to the WDW equation satisfy the DeWitt boundary condition. We find that by using a suitable Born–Oppenheimer (BO) approximation, whose validity is proven, the DeWitt condition is satisfied. Furthermore, the big rip singularity is expected to be avoided in the quantum realm.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Multigrid deflation for Lattice QCD

Computing the trace of the inverse of large matrices is typically addressed through statistical methods. Deflating out the lowest eigenvectors or singular vectors of the matrix reduces the variance of the trace estimator. This work summarizes our efforts to reduce the computational cost of computing the deflation space while achieving the desired variance reduction for Lattice QCD applications. Previous efforts computed the lower part of the singular spectrum of the Dirac operator by using an eigensolver preconditioned with a multigrid linear system solver. Despite the improvement in performance in those applications, as the problem size grows the runtime and storage demands of this approach will eventually dominate the stochastic estimation part of the computation. In this work, we propose to compute the deflation space in one of the following two ways. First, by using an inexact eigensolver on the Hermitian, but maximally indefinite, operator. Second, by exploiting the fact that the multigrid prolongator for this operator is rich in components toward the lower part of the singular spectrum. We show experimentally that the inexact eigensolver can approximate the lower part of the spectrum even for ill-conditioned operators. Also, the deflation based on the multigrid prolongator is more efficient to compute and apply, and, despite its limited ability to approximate the fine level spectrum, it obtains similar variance reduction on the trace estimator as deflating with approximate eigenvectors from the fine level operator.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Multicycle flash sintering of cubic Y 2 O 3 -stabilized ZrO 2 : An in situ energy dispersive synchrotron x-ray diffraction study with high temporal resolution

The current induced unit cell volume changes, (111) Bragg peak full width at half maximum (FWHM) and its integrated intensity in 8 % Y 2 O 3 stabilized ZrO 2 (8 %YSZ) solid state electrolyte was monitored during a triple-flash sintering experiment by in situ energy dispersive x-ray diffraction using a polychromatic synchrotron probe (max, photon energy 200 keV) with 2 second temporal resolution. The first spontaneous singularity in the unit cell volume (+0.54 %) was observed at 899 °C under 15 V/mm applied field intensity, which was associated with 13 mA/mm 2 current draw and an increase in density to 97 %. Following anelastic relaxation of the unit cell volume under open circuit conditions, the same applied field was applied twice in a row which resulted in additional induced singularities at 925 °C (+0.48 %) and 944 °C (+0.42 %). A floating baseline, which was above the thermal expansion baseline, was observed from 833 to 969 °C and was attributed to Joule heating. In conclusion, the singularity at 899 °C is associated with a sharp change in (111) FWHM and a 34 % decrease in integrated peak area that was attributed to changes in the distribution of oxygen vacancies and the changes in their concentration as induced by the applied field in the spontaneous transient stage of flash sintering.

36 MATERIALS SCIENCE↗
DOE OSTI2019

The fuzzball nature of two-charge black hole microstates

It has been suggested by A. Sen that the entropy of two-charge supersymmetric bound states in string theory should be accounted for by adding the entropy of source-free horizonless supergravity solutions to the entropy associated with the horizons of small black holes. This would imply that the entropy arises differently depending on the duality frame: in the D1-D5 frame one would count source-free horizonless solutions, while in the NS1-P frame one would compute the area of a horizon. This might lead to the belief that the microstates are described by fuzzball solutions in the D1-D5 duality frame but by a black hole with a horizon in the latter. We argue that this is not the case, and that the microstates are fuzzballs in all duality frames. We observe that the scaling argument used by Sen fails to account for the entropy in the D1-P and other duality frames. We also note that the traditional extremal black hole solution is not a complete string background, since finite-action paths connect the exterior near-horizon extremal throat to the region inside the horizon, including the singularity. The singularity of the traditional black hole solution does not give a valid boundary condition for a fundamental string; correcting this condition by resolving the singularity modifies the black hole to a fuzzball with no horizon. We argue that for questions of counting states, the traditional black hole solution should be understood through its Euclidean continuation as a saddle point, and that the Lorentzian states being counted are fuzzballs in all duality frames.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Solving gyrokinetic systems with higher-order time dependence

We discuss theoretical and numerical aspects of gyrokinetics as a Lagrangian field theory when the field perturbation is introduced into the symplectic part. A consequence is that the field equations and particle equations of motion in general depend on the time derivatives of the field. The most well-known example is when the parallel vector potential is introduced as a perturbation, where a time derivative of the field arises only in the equations of motion, so an explicit equation for the fields may still be written. We will consider the conceptually more problematic case where the time-dependent fields appear in both the field equations and equations of motion, but where the additional term in the field equations is formally small. The conceptual issues were described by Burby (J. Plasma Phys., vol. 82 (3), 2016, 905820304): these terms lead to apparent additional degrees of freedom to the problem, so that the electric field now requires an initial condition, which is not required in low-frequency (Darwin) Vlasov–Maxwell equations. Also, the small terms in the Euler–Lagrange equations are a singular perturbation, and these two issues are interlinked. For well-behaved problems the apparent additional degrees of freedom are spurious, and the physically relevant solution may be directly identified. Because we needed to assume that the system is well behaved for small perturbations when deriving gyrokinetic theory, we must continue to assume that when solving it, and the physical solutions are thus the regular ones. The spurious nature of the singular degrees of freedom may also be seen by changing coordinate systems so the varying field appears only in the Hamiltonian. We then describe how methods appropriate for singular perturbation theory may be used to solve these asymptotic equations numerically. We then describe a proof-of-principle implementation of these methods for an electrostatic strong-flow gyrokinetic system; two basic test cases are presented to illustrate code functionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Direct Spectroscopic Evidence for Phase Competition between the Pseudogap and Superconductivity in Bi 2 Sr 2 CaCu 2 O 8+δ

In the high-temperature (T c ) cuprate superconductors, increasing evidence suggests that the pseudogap, existing below the pseudogap temperature T*, has a distinct broken electronic symmetry from that of superconductivity. Particularly, recent scattering experiments on the underdoped cuprates have suggested that a charge ordering competes with superconductivity. However, no direct link of this physics and the important low-energy excitations has been identified. We report an antagonistic singularity at Tc in the spectral weight of Bi 2 Sr 2 CaCu 2 O 8+δ as a compelling evidence for phase competition, which persists up to a high hole concentration p ~ 0.22. Comparison with a theoretical calculation confirms that the singularity is a signature of competition between the order parameters for the pseudogap and superconductivity. Our observation of the spectroscopic singularity at finite temperatures over a wide doping range provides new insights into the nature of the competitive interplay between the two intertwined phases and the complex phase diagram near the pseudogap critical point.

36 MATERIALS SCIENCE↗