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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

An Overview of Gemma FY20 Verification Activities

Gemma verification activities for FY20 can be divided into three categories: the development of specialized quadrature rules, initial progress towards the development of manufactured solutions for code verification, and automated code-verification testing. In the method-of-moments implementation of the electric-field integral equation, the presence of a Green’s function in the four-dimensional integrals yields singularities in the integrand when two elements are nearby. To address these challenges, we have developed quadrature rules to integrate the functions through which the singularities can be characterized. Code verification is necessary to develop confidence in the implementation of the numerical methods in Gemma. Therefore, we have begun investigating the use of manufactured solutions to more thoroughly verify Gemma. Manufactured solutions provide greater flexibility for testing aspects of the code; however, the aforementioned singularities provide challenges, and existing work is limited in rigor and quantity. Finally, we have implemented automated code-verification testing using the VVTest framework to automate the mesh refinement and execution of a Gemma simulation to generate mesh convergence data. This infrastructure computes the observed order of accuracy from these data and compares it with the theoretical order of accuracy to either develop confidence in the implementation of the numerical methods or detect coding errors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Using Pilot Jobs and CernVM File System for Simplified Use of Containers and Software Distribution

High Energy Physics (HEP) experiments entail an abundance of computing resources, i.e. sites, to run simulations and analyses by processing data. This requirement is fulfilled by local batch farms, grid sites, private/commercial clouds, and supercomputing centers via High Throughput Computing (HTC). The growing needs of such experiments and resources being prone to trends of heterogeneity make it difficult for physicists to handle these resources directly. Additionally, HEP collaborations heavily rely on data and software releases, typically in the order of tens of gigabytes, while conducting simulations and analyses. Hence, aspects of scalability, reliability, and maintenance become crucial with regards to the distribution of the necessary data and software stack. The GlideinWMS [4] framework helps with the resource management problem by using pilot jobs, aka Glideins, to provision reliable elastic virtual clusters. Glideins are submitted to unreliable heterogeneous resources which are validated and customized by the Glideins to make the worker nodes available for end-user job execution. On the other hand, the CernVM File System (CernVM-FS or CVMFS) [1] helps with data distribution. It is a write-once, read-everywhere filesystem used to deploy scientific software to thousands of nodes on a worldwide distributed computing infrastructure. CVMFS is based on the Hyper Text Transfer Protocol and has been widely used within the particle physics community for (1) distributing experiment software and data such as calibrations, and (2) facilitating containerization by efficiently hosting container images along with providing containerization software, especially Singularity [3] GlideinWMS relies on CVMFS installed locally on the computing resources to satisfy the experiments' software needs. This requires system administrators' effort to install and maintain CVMFS at the sites and limits the use of sites, especially HPC resources, that do not have CVMFS installed. This poster presents a solution, taking advantage of Glideins to provide CVMFS at most sites without the need for a local installation. Doing so expands the pool of resources available for HEP experiments and reduces the effort of system administrators for current resources. Additionally, the proposed solution allows GlideinWMS to also start Singularity [3], a containerization software that can run unprivileged, on sites where neither CVMFS nor Singularity are available, including HPC sites. The benefits provided by this solution are: (1) lower overhead for site administrators in that they have less software to install, (2) an expanded pool of resources that run user jobs with easy access to software and data provided by CVMFS, thus making life easier for the scientists, and (3) improved flexibility to use HPC resources by enabling GlideinWMS pilot jobs to support HPC sites.

Urs, Namratha↗
DOE OSTI2024

Black Box Equations of State: Creating Semi-analytic Solutions to the Noh Problem and Verifying Equation of State Interfaces

The objective of this report is threefold. First, it details a method for deriving a semi-analytic solution to the Noh Problem when using a “black-box” equation of state. Such capability allows us to perform verification on complicated, more realistic equations of state. Examples include Steinberg equations of state for materials and tabulated equations of state. The second objective is to apply the methodology to verify the singularity-eos equation of state library. We do so by solving the Rankine-Hugoinot jump conditions for the Noh Problem, ensuring singularity derives the correct solution and comparing the error to an exact implementation of the equation of state. The third objective is to perform verification of the xRAGE Eulerian hydrodynamics code when interfaced with singularity. We provide the theory, analysis, documentation for a python implementation of the proposed solver, and verification results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Guidelines for predicting stress in cemented doublets undergoing temperature change

This work explores quick predictive methods for calculating potentially risky stresses and deflections in cemented doublets experiencing temperature change that agree well with finite element analysis. There are three failure modes of interest: cohesive failure of the adhesive, delamination (surface bond failure or debonding), and glass fracture. Adhesive theory, confirmed by finite element analysis, predicts stress singularities that complicate interpretation of the stress calculations. The presence of a stress singularity indicates the breakdown of linear elastic assumptions, but damage initiation and stress singularities are related. The authors find that geometry details near a bond edge can exacerbate or minimize damage initiation and stress concentrations. Because the interpretation of the stress results is complicated, the authors investigated predicted stresses in doublets that have been successfully tested between -40 and 85 °C. This study found that the product (ΔT ∙ Δα) should be less than 189 ppm, where ΔT is the temperature excursion and Δα is the difference in glass coefficient of thermal expansion. If the product (ΔT ∙ Δα) is equal to or greater than 189 ppm, further analysis and testing is warranted. But the authors also show that the fabrication process can significantly influence stress failure, particularly with large diameter doublets.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Guidelines for predicting stress in cemented doublets undergoing temperature change

This work explores quick predictive methods for calculating potentially risky stresses and deflections in cemented doublets experiencing temperature change that agree well with finite element analysis. There are three failure modes of interest: cohesive failure of the adhesive, delamination (surface bond failure or debonding), and glass fracture. Adhesive theory, confirmed by finite element analysis, predicts stress singularities that complicate interpretation of the stress calculations. The presence of a stress singularity indicates the breakdown of linear elastic assumptions, but damage initiation and stress singularities are related. The authors find that geometry details near a bond edge can exacerbate or minimize damage initiation and stress concentrations. Because the interpretation of the stress results is complicated, the authors investigated predicted stresses in doublets that have been successfully tested between -40 and 85 °C. This study found that the product (∆T ∙ ∆α) should be less than 189 ppm, where ∆T is the temperature excursion and ∆α is the difference in glass coefficient of thermal expansion. If the product (∆T ∙ ∆α) is equal to or greater than 189 ppm, further analysis and testing is warranted. But the authors also show that the fabrication process can significantly influence stress failure, particularly with large diameter doublets

42 ENGINEERING↗
DOE OSTI2025

Verification of a Fluid-Based Plasma-Edge Model Within the Multiphysics Object-Oriented Simulation Environment (MOOSE) Framework

As the goal of achieving fusion power on the grid comes closer to fruition, fully coupled multiphysics models of fusion devices will be crucial. These models must incorporate the interconnected phenomena of these devices, including plasma physics, neutronics, first wall interactions, and tritium transport. Currently, there are two main approaches to developing these platforms: (1) loosely coupled, where one couples existing codes and solvers together through input and output parameters and data, and (2) tightly coupled, where one develops the necessary models within a singular, integrated framework. This work focuses on the latter approach for magnetically confined fusion devices by developing a fluid-based plasma-edge model within the Multiphysics Object Oriented Simulation Environment (MOOSE) Framework. This effort is coordinated with other efforts to develop, test, demonstrate, and deploy fusion relevant multiphysics capabilities including electromagnetics, particle-in-cell plasma, tritium transport, and fusion blanket design. This new model is an expansion of the MOOSE-based plasma application, Zapdos, which was originally formulated to model low-temperature, non-magnetized plasma processes. Verification studies have been conducted using newly developed magnetic plasma capabilities. These involved convergence analyses utilizing the method of manufactured solutions to verify new operators and case studies. A modular approach was taken here to demonstrate increasingly complicated simulation scenarios, which included a singular fluid with uniform magnetic field case, a singular fluid with spatially varying magnetic field case, and a coupled multifluid case.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗