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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Unveiling horizons in quantum critical collapse

Critical gravitational collapse offers a unique window into regimes of arbitrarily high curvature, culminating in a naked singularity arising from smooth initial data — thus providing a dynamical counterexample to weak cosmic censorship. Near the critical regime, quantum effects from the collapsing matter are expected to intervene before full quantum gravity resolves the singularity. Despite its fundamental significance, a self-consistent treatment has so far remained elusive. In this work, we perform a one-loop semiclassical analysis using the robust anomaly-based method in the canonical setup of Einstein gravity minimally coupled to a free, massless scalar field. Focusing on explicitly solvable near-critical solutions in both 2 + 1 and 3 + 1 dimensions, we analytically solve the semiclassical Einstein equations and obtain controlled, quantitative results for several long-standing questions within the dominant s-wave sector. We find that regularity uniquely selects a Boulware-like quantum state, encoding genuine vacuum polarization effects from the collapsing matter. Remarkably, the resulting quantum corrections manifest as a growing mode. Horizon-tracing analyses, incorporating both classical and quantum modes, reveal the emergence of a finite mass gap, signaling a phase transition from classical Type II to quantum-modified Type I behavior, thereby providing a quantum enforcement of the weak cosmic censorship. The most non-trivial aspect of our analysis involves dealing with non-conformal matter fields in explicitly time-dependent critical spacetimes. Along the way, we uncover intriguing and previously underexplored features of quantum field theory in curved spacetime.

2D Gravity↗
DOE OSTI2020

Gravitational duals to the grand canonical ensemble abhor Cauchy horizons

The gravitational dual to the grand canonical ensemble of a large N holographic theory is a charged black hole. These spacetimes - for example Reissner- Nordström-AdS - can have Cauchy horizons that render the classical gravitational dynamics of the black hole interior incomplete. We show that a (spatially uniform) deformation of the CFT by a neutral scalar operator generically leads to a black hole with no inner horizon. There is instead a spacelike Kasner singularity in the interior. For relevant deformations, Cauchy horizons never form. For certain irrelevant deformations, Cauchy horizons can exist at one specific temperature. We show that the scalar field triggers a rapid collapse of the Einstein-Rosen bridge at the would-be Cauchy horizon. Finally, we make some observations on the interior of charged dilatonic black holes where the Kasner exponent at the singularity exhibits an attractor mechanism in the low temperature limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Observation of Flat Bands in Gated Semiconductor Artificial Graphene

Flat bands near M points in the Brillouin zone are key features of honeycomb symmetry in artificial graphene (AG) where electrons may condense into novel correlated phases. Here we report the observation of van Hove singularity doublet of AG in GaAs quantum well transistors, which presents the evidence of flat bands in semiconductor AG. Two emerging peaks in photoluminescence spectra tuned by backgate voltages probe the singularity doublet of AG flat bands and demonstrate their accessibility to the Fermi level. As the Fermi level crosses the doublet, the spectra display dramatic stability against electron density, indicating interplays between electron-electron interactions and honeycomb symmetry. Our results provide a new flexible platform to explore intriguing flat band physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Three-dimensional higher-order saddle-point-induced flatbands in Co-based kagome metals

The saddle point (Van Hove singularity) exhibits a divergent density of states in two-dimensional systems, leading to fascinating phenomena such as strong correlations and unconventional superconductivity, yet it is seldom observed in three-dimensional (3D) systems. In this work we find two types of 3D higher-order saddle points (HOSPs) in emerging 3D kagome metals YbCo 6⁢ Ge 6 and MgCo 6 ⁢Ge 6 . Both HOSPs exhibit a singularity in their density of states, which is significantly enhanced compared to the ordinary saddle point. The HOSP near the Fermi energy generates a flatband extending a large area in the Brillouin zone, potentially amplifying the correlation effect and fostering electronic instabilities. Two types of HOSPs exhibit distinct robustness upon element substitution and lattice distortions in these kagome compounds. Our work paves the way for engineering exotic band structures, such as saddle points and flatbands, and exploring interesting phenomena in Co-based kagome materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Investigation of Point-Contact Strategies for CFD Simulations of Pebble-Bed Reactor Cores

This study numerically investigated the effects of various contact strategies on the thermal hydraulic behavior within a structured bed of 100 explicitly modeled pebbles. Four contact strategies and two thermal hydraulic conditions were considered. The strategies to avoid contact singularities include decreasing the pebble diameter, increasing the pebble diameter, bridging the pebble surfaces near the contact region, and capping the pebble surfaces near the contact region. One strategy, Strategy 3a, which involves bridging with a cylinder equal to 10% of the pebble diameter, was selected as the baseline strategy because it addressed the contact singularity while minimizing the geometric changes that affect the bed porosity. The two thermal hydraulic conditions were full-power operation (Case 1) and pressurized loss of forced cooling or PLOFC (Case 2). Simulations of the conjugate heat transfer within the structured bed were performed using the Reynolds-averaged Navier–Stokes approach with the realizable k-ϵ turbulence model and two-layer all y + wall treatment. The thermal-fluid quantities of interest were compared between the contact strategies for each case. In Case 1, the hydraulic behavior was sensitive to the contact strategy, with large differences in the pressure drop (30%) and volume-average velocity (4%). The thermal behavior was not sensitive, with less than a 0.5% difference across the strategies. To better understand the separate effects of each heat transfer mode, Case 2 was divided into the following subcases: conduction (Case 2a); conduction/radiation (Case 2b); and conduction/radiation/convection (Case 2c). Case 2a represents an early phase of the PLOFC transient. Case 2b represents an intermediate phase of the PLOFC transient, with the pebble temperatures sufficiently high for the radiative heat transfer to be non-negligible. Case 2c represents a late phase of the PLOFC transient after the establishment of the natural circulation of the heat transfer fluid. For Case 2, large differences in the contact strategy were observed only in Case 2a with only conduction. The difference in the maximum pebble temperature was 23% in Case 2a, 2% in Case 2b, and 0.3% in Case 2c.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

On the BPS Sector in AdS 3 /CFT 2 Holography

Abstract The BPS sector in duality has been fertile ground for the exploration of gauge/gravity duality, from the match between black hole entropy and the CFT elliptic genus to the construction of large families of geometrical microstates and the identification of the corresponding states in the CFT. Worldsheet methods provide a tool to further explore the relation between string theory in the bulk and corresponding CFT quantities. We show how to match individual BPS strings to their counterparts in the symmetric product orbifold CFT. In the process, we find an exact match between known constructions of microstate geometries and condensates of BPS supergraviton strings, and discuss their role in the broader collection of BPS states. In particular, we explore how microstate geometries develop singularities; and how string theory resolves these singularities through the appearance of “tensionless” string dynamics, which is the continuation of structures found in the weak‐coupling CFT into the strongly coupled regime described by string theory in the bulk. We argue that such “tensionless” strings are responsible for black hole microstructure in the bulk description.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Factorization at subleading power and endpoint divergences in h → γγ decay. Part II. Renormalization and scale evolution

Building on the recent derivation of a bare factorization theorem for the b-quark induced contribution to the h → γγ decay amplitude based on soft-collinear effective theory, we derive the first renormalized factorization theorem for a process described at subleading power in scale ratios, where λ = m b /M h « 1 in our case. We prove two refactorization conditions for a matching coefficient and an operator matrix element in the endpoint region, where they exhibit singularities giving rise to divergent convolution integrals. The refactorization conditions ensure that the dependence of the decay amplitude on the rapidity regulator, which regularizes the endpoint singularities, cancels out to all orders of perturbation theory. We establish the renormalized form of the factorization formula, proving that extra contributions arising from the fact that “endpoint regularization” does not commute with renormalization can be absorbed, to all orders, by a redefinition of one of the matching coefficients. We derive the renormalization-group evolution equation satisfied by all quantities in the factorization formula and use them to predict the large logarithms of order $αα^2_sL^k$ in the three-loop decay amplitude, where $L=ln(-M^2_h/m^2_b)$ and k = 6, 5, 4, 3. We find perfect agreement with existing numerical results for the amplitude and analytical results for the three-loop contributions involving a massless quark loop. On the other hand, we disagree with the results of previous attempts to predict the series of subleading logarithms $~αα^n_sL^{2n+1}$

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Islands in the fluid: islands are common in cosmology

We discuss the possibility of entanglement islands in cosmological spacetimes with a general perfect fluid with an equation of state w. We find that flat universes with time-symmetric slices where the Hubble parameter vanishes always have islands on that slice. We then move away from such slices, considering still universes with a general perfect fluid. Under the local thermal equilibrium assumption, the comoving entropy density s c is constant. As a result, the conditions for an island become an inequality between the energy density (or Hubble parameter) and the temperature at some time of normalization. The consequences are that islands can exist for practically all fluids that are not radiation, i.e. w ≠ 1/3. We also discuss the ramifications of our results for universes with spatial curvature. Finally, we show that islands occur in the Simple Harmonic Universe model which has no classical singularity at the background level, in contrast to all previous examples where islands occurred only in space-times with singularities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Bootstrapping M-theory orbifolds

Abstract We analyze correlation functions of SU(k) × SU(2) F flavor currents in a family of three-dimensional$$ \mathcal{N} $$ N = 4 superconformal field theories, combining analytic bootstrap methods with input from supersymmetric localization. Via holographic duality, we extract gluon and graviton scattering amplitudes of M-theory on AdS 4 ×S 7 /ℤ k which contains a ℂ 2 /ℤ k orbifold singularity. From these results, we derive aspects of the effective description of M-theory on the orbifold singularity beyond its leading low energy limit. We also determine a threshold correction to the holographic correlator from the combined contribution of two-loop gluon and tree-level bulk graviton exchange.

Physics↗
DOE OSTI2026

Toward the analytic bootstrap of energy correlators

In this paper, we present a framework for the analytic bootstrap of three-point energy correlators, a crucial observable in $\mathcal{N}$ = 4 super Yang-Mills theory and quantum chromodynamics (QCD). Our approach combines spherical contour techniques, general physical constraints such as pole cancellations, and power correction data in the singular limits to determine its analytic expression. In contrast to previous bootstrap studies restricted to scattering amplitudes for supersymmetric theories, our framework makes use of the properties of Feynman integrals, marking a significant step toward bootstrapping realistic QCD observables. Using this method, we derive analytic expressions for leading-order three-point energy correlators in the multi-collinear limit with equal and unequal energy weights, where the latter are crucial ingredients for projected N -point energy correlators. We also apply the recently developed technique of analytic regression with lattice reduction as a way to bypass needing explicit expressions for the singular limits. Bridging theoretical advances in scattering amplitudes with the renewed interest in weighted cross-sections, our work opens the door to precision tests of QCD dynamics through analytic event-shape predictions.

jet substructure↗
DOE OSTI2022

Landau discriminants

Scattering amplitudes in quantum field theories have intricate analytic properties as functions of the energies and momenta of the scattered particles. In perturbation theory, their singularities are governed by a set of nonlinear polynomial equations, known as Landau equations, for each individual Feynman diagram. The singularity locus of the associated Feynman integral is made precise with the notion of the Landau discriminant, which characterizes when the Landau equations admit a solution. In order to compute this discriminant, we present approaches from classical elimination theory, as well as a numerical algorithm based on homotopy continuation. These methods allow us to compute Landau discriminants of various Feynman diagrams up to 3 loops, which were previously out of reach. For instance, the Landau discriminant of the envelope diagram is a reducible surface of degree 45 in the three-dimensional space of kinematic invariants. We investigate geometric properties of the Landau discriminant, such as irreducibility, dimension and degree. In particular, we find simple examples in which the Landau discriminant has codimension greater than one. Furthermore, we describe a numerical procedure for determining which parts of the Landau discriminant lie in the physical regions. In order to study degenerate limits of Landau equations and bounds on the degree of the Landau discriminant, we introduce Landau polytopes and study their facet structure. Finally, we provide an efficient numerical algorithm for the computation of the number of master integrals based on the connection to algebraic statistics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Gravitational index of the heterotic string

The fundamental heterotic string has a tower of BPS states whose supersymmetric index has an exponential growth in the charges. We construct the saddle-point of the gravitational path integral corresponding to this index. The saddle-point configuration is a supersymmetric rotating non-extremal Euclidean black hole. This configuration is singular in the two-derivative theory. We show that the addition of higher-derivative terms in four-dimensional N = 2 supergravity resolves the singularity. In doing so, we extend the recently-developed “new attractor mechanism” to include the effect of higher-derivative terms. Remarkably, the one-loop, four-derivative F-term contribution to the prepotential leads to a precise match of the gravitational and microscopic index. We also comment, using the effective theory near the horizon, on the possibility of a string-size near-extremal black hole. Our results clarify the meaning of different descriptions of this system in the literature. The thermal state transitions to a winding condensate and a gas of strings without ever reaching a small black hole, while the index is captured by the rotating Euclidean black hole solution and is constant and thus smoothly connected to the microscopic ensemble.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Strings in irrelevant deformations of AdS 3 /CFT 2

We generalize our recent analysis [1] of probe string dynamics to the case of general single-trace TT¯,JT¯ and TJ¯ deformations. We show that in regions in coupling space where the bulk geometry is smooth, the classical trajectories of such strings are smooth and approach the linear dilaton boundary at either the far past or the far future. These trajectories give rise to quantum scattering states with arbitrarily high energies. When the bulk geometry has closed timelike curves (CTC’s), the trajectories are singular for energies above a critical value E c . This singularity occurs in the region with CTC’s, and the value of E c agrees with that read off from the dual boundary theory for all values of the couplings and charges.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Towards a classification of rank r $\mathscr{N}$ = 2 SCFTs. Part II. Special Kahler stratification of the Coulomb branch

We study the stratification of the singular locus of four dimensional $\mathscr{N}$ = 2 Coulomb branches. We present a set of self-consistency conditions on this stratification which can be used to extend the classification of scale-invariant rank 1 Coulomb branch geometries to two complex dimensions, and beyond. The calculational simplicity of the arguments presented here stems from the fact that the main ingredients needed - the rank 1 deformation patterns and the pattern of inclusions of rank 2 strata - are discrete topological data which satisfy strong self-consistency conditions through their relationship to the central charges of the SCFT. This relationship of the stratification data to the central charges is used here, but is derived and explained in a companion paper by one of the authors. We illustrate the use of these conditions by re-analyzing many previously-known examples of rank 2 SCFTs, and also by finding examples of new theories. The power of these conditions stems from the fact that for Coulomb branch stratifications a conjecturally complete list of physically allowed “elementary slices” is known. By contrast, constraining the possible elementary slices of symplectic singularities relevant for Higgs branch stratifications remains an open problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Improved Statistics for F-theory Standard Models

Much of the analysis of F-theory-based Standard Models boils down to computing cohomologies of line bundles on matter curves. By varying parameters one can degenerate such matter curves to singular ones, typically with many nodes, where the computation is combinatorial and straightforward. The question remains to relate the (a priori possibly smaller) value on the original curve to the singular one. In this work, we introduce some elementary techniques (pruning trees and removing interior edges) for simplifying the resulting nodal curves to a small collection of terminal ones that can be handled directly. When applied to the QSMs, these techniques yield optimal results in the sense that obtaining more precise answers would require currently unavailable information about the QSM geometries. This provides us with an opportunity to enhance the statistical bounds established in earlier research regarding the absence of vector-like exotics on the quark-doublet curve.

Bies, Martin↗
DOE OSTI2024

An unstructured body-of-revolution electromagnetic particle-in-cell algorithm with radial perfectly matched layers and dual polarizations

A novel electromagnetic particle-in-cell algorithm has been developed for fully kinetic plasma simulations on unstructured (irregular) meshes in complex body-of-revolution geometries. The algorithm, implemented in the BORPIC++ code, utilizes a set of field scalings and a coordinate mapping, reducing the Maxwell field problem in a cylindrical system to a Cartesian finite element Maxwell solver in the meridian plane. The latter obviates the cylindrical coordinate singularity in the symmetry axis. The choice of an unstructured finite element discretization enhances the geometrical flexibility of the BORPIC++ solver compared to the more traditional finite difference solvers. Symmetries in Maxwell’s equations are explored to decompose the problem into two dual polarization states with isomorphic representations that enable code reuse. The particle-in-cell scatter and gather steps preserve charge conservation at the discrete level. Our previous algorithm (BORPIC+) discretized the E and B field components of TE Φ and TM Φ polarizations on the finite element (primal) mesh. Here, we employ a new field-update scheme. Using the same finite element (primal) mesh, this scheme advances two sets of field components independently: (1) E and B of TE Φ polarized fields, (E z , E ρ , B Φ ) and (2) D and H of TM Φ polarized fields, (D Φ , H z , H ρ ). Since these field updates are not explicitly coupled, the new field solver obviates the coordinate singularity, which otherwise arises at the cylindrical symmetric axis, ρ = 0 when defining the discrete Hodge matrices (generalized finite element mass matrices). Here, a cylindrical perfectly matched layer is implemented as a boundary condition in the radial direction to simulate open space problems, with periodic boundary conditions in the axial direction. We investigate effects of charged particles moving next to the cylindrical perfectly matched layer. We model azimuthal currents arising from rotational motion of charged rings, which produce TMΦ polarized fields. Several numerical examples are provided to illustrate the first application of the algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Accelerating multigrid with streaming chiral SVD for Wilson fermions in lattice QCD

A modification to the setup algorithm for the multigrid preconditioner of Wilson fermions in lattice QCD is presented. A larger basis of test vectors than that used in regular multigrid is calculated by the smoother and truncated by singular value decomposition on the chiral components of the test vectors. The truncated basis is used to form the prolongation and restriction matrices of the multigrid hierarchy. This modification of the setup method is demonstrated to increase the convergence of linear solvers on an anisotropic lattice with m π ≈ 239 MeV from the Hadron Spectrum Collaboration and an isotropic lattice with m π ≈ 220 MeV from the MILC Collaboration. The lattice volume dependence of the method is also examined. Increasing the number of test vectors improves speedup up to a point, but storing these vectors becomes impossible in limited memory resources such as GPUs. To address storage cost, we implement a streaming singular value decomposition of the basis of test vectors on the chiral components and demonstrate a decrease in the number of fine level iterations by a factor of 1.7 for m q ≈ m crit

Iterative methods↗
DOE OSTI2023

The physics of desiccation cracks 2: Modeling and prediction of the crack patterns

Here, this paper extends the experimental results of our companion paper by modeling and predicting the onset and pattern formation of desiccation cracks in geomaterials (Ruoyu et al., 2023). Thin-layer silt samples in controlled atmospheric conditions were tested to obtain the surface strain maps with the digital image correlation (DIC) method during dehydration. Support experiments, including consolidation and displacement-controlled triaxial tests, were conducted for the properties of geomaterials. These experimental results were used to validate a viscoplastic theoretical model by comparing cracks and singularities locations that distribute following the Cnoidal wave pattern. A critical value exists in the viscoplastic model that determines the number of singularities, which accurately predicts the number of cracks in the experiment. The critical values contain two crucial parameters: rate sensitivity and λ. Rate sensitivity describes the rate-dependent stress–strain relation, while λ shows the ability of pore pressure redistribution under the external mechanical loading rate. These results provide a new view to analyze the desiccation cracks considering the rate-dependent viscoplasticity.

42 ENGINEERING↗