Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “discretization”

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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Holography from lattice $\mathcal{N}$ = 4 super Yang-Mills

In this paper we use lattice simulation to study four dimensional $\mathcal{N}$ = 4 super Yang-Mills (SYM) theory. We have focused on the three color theory on lattices of size 12 4 and for ’t Hooft couplings up to λ = 40.0. Our lattice action is based on a discretization of the Marcus or GL twist of $\mathcal{N}$ = 4 SYM and retains one exact supersymmetry for non-zero lattice spacing. We show that lattice theory exists in a single non-Abelian Coulomb phase for all ’t Hooft couplings. Furthermore the static potential we obtain from correlators of Polyakov lines is in good agreement with that obtained from holography — specifically the potential has a Coulombic form with a coefficent that varies as the square root of the ’t Hooft coupling.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Jet rotational metrics

Abstract Embedding symmetries in the architectures of deep neural networks can improve classification and network convergence in the context of jet substructure. These results hint at the existence of symmetries in jet energy depositions, such as rotational symmetry, arising from the physical features of the underlying processes. We introduce new jet observables, Jet Rotational Metrics (JRMs), which provide insights into the substructure of jets by comparing them to jets with perfect discrete rotational symmetry. We show that JRMs are formidable jet features, achieving good classification scores when used as inputs to deep neural networks. We also show that when used in combination with other jet observables, like N-subjettiness and EFPs, our features increase classification performance. The results suggest that JRMs may capture information not efficiently captured by the other observables, motivating the design of future jet observables for learning the underlying symmetries in the physical processes.

Physics↗
DOE OSTI2025

Entanglement maximization and mirror symmetry in two-Higgs-doublet models

We consider 2-to-2 scatterings of Higgs bosons in a CP-conserving two-Higgs-doublet model (2HDM) and study the implication of maximizing the entanglement in the flavor space, where the two doublets Φ a , a = 1, 2, can be viewed as a qubit: Φ 1 = |0⟩ and Φ 2 = |1⟩. More specifically, we compute the scattering amplitudes for Φ a Φ b → Φ c Φ d and require the outgoing flavor entanglement to be maximal for a full product basis such as the computational basis, which consists of {|00⟩, |01⟩, |10⟩, |11⟩}. In the unbroken phase and turning off the gauge interactions, entanglement maximization results in the appearance of an U(2) × U(2) global symmetry among the quartic couplings, which in general is broken softly by the mass terms. Interestingly, once the Higgs bosons acquire vacuum expectation values, maximal entanglement enforces an exact U(2) × U(2) symmetry, which is spontaneously broken to U(1) × U(1). As a byproduct, this gives rise to Higgs alignment as well as to the existence of 6 massless Nambu-Goldstone bosons. The U(2) × U(2) symmetry can be gauged to lift the massless Goldstones, while maintaining maximal entanglement demands the presence of a discrete Z 2 symmetry interchanging the two gauge sectors. The model is custodially invariant in the scalar sector, and the inclusion of fermions requires a mirror dark sector, related to the standard one by the Z 2 symmetry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Non-invertible global symmetries and completeness of the spectrum

It is widely believed that consistent theories of quantum gravity satisfy two basic kinematic constraints: they are free from any global symmetry, and they contain a complete spectrum of gauge charges. For compact, abelian gauge groups, completeness follows from the absence of a 1-form global symmetry. However, this correspondence breaks down for more general gauge groups, where the breaking of the 1-form symmetry is insufficient to guarantee a complete spectrum. We show that the correspondence may be restored by broadening our notion of symmetry to include non-invertible topological operators, and prove that their absence is sufficient to guarantee a complete spectrum for any compact, possibly disconnected gauge group. In addition, we prove an analogous statement regarding the completeness of twist vortices: codimension-2 objects defined by a discrete holonomy around their worldvolume, such as cosmic strings in four dimensions. We discuss how this correspondence is modified in various, more general contexts, including non-compact gauge groups, Higgsing of gauge theories, and the addition of Chern-Simons terms. Finally, we discuss the implications of our results for the Swampland program, as well as the phenomenological implications of the existence of twist strings.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Symmetry transmutation and anomaly matching

We explore a situation where a global symmetry of the ultraviolet (UV) theory does not act faithfully on the local infrared (IR) degrees of freedom, but instead acts effectively as a higher-form symmetry. We refer to this phenomenon as symmetry transmutation, where the UV symmetry is “transmuted” into a higher-form symmetry in the IR. Notably, unlike emergent (accidental) symmetries, which are approximate, these symmetries are exact. We illustrate the ubiquity of this phenomenon in various continuum and lattice systems and provide examples where the ’t Hooft anomalies of the UV symmetry are matched by those of the new higher-form symmetry in the IR. We also show that in certain phases and for certain energies, the UV baryon-number symmetry of one-flavor QCD is transmuted into a discrete one-form global symmetry. Finally, we compare our symmetry transmutation to the well-known phenomenon of symmetry fractionalization.

Anomalies in Field and String Theories↗
DOE OSTI2023

Double Kerr-Schild spacetimes and the Newman-Penrose map

The Newman-Penrose map, which is closely related to the classical double copy, associates certain exact solutions of Einstein’s equations with self-dual solutions of the vacuum Maxwell equations. Here we initiate an extension of the Newman-Penrose map to a broader class of spacetimes. As an example, we apply principles from the Newman-Penrose map to associate a self-dual gauge field to the Kerr-Taub-NUT-(A)dS spacetime and we show that the result agrees with previously studied examples of classical double copies. The corresponding field strength exhibits a discrete electric-magnetic duality that is distinct from its (Hodge star) self-dual property.

scattering amplitudes↗
DOE OSTI2023

AdS 3 orbifolds, BTZ black holes, and holography

Conical defects of the form (AdS 3 × S 3 )/Z k have an exact orbifold description in worldsheet string theory, which we derive from their known presentation as gauged Wess-Zumino-Witten models. The configuration of strings and fivebranes sourcing this geometry is well-understood, as is the correspondence to states/operators in the dual CFT 2 . One can analytically continue the construction to Euclidean AdS 3 (i.e. the hyperbolic ball $H$$^{+}_{3}$) and consider the orbifold by any infinite discrete (Kleinian) group generated by a set of elliptic elements γi ϵ SL(2, C), $γ$$^{ki}_{i}$ = 1, i = 1, . . . , K. The resulting geometry consists of multiple conical defects traveling along geodesics in $H$$^{+}_{3}$, and provides a semiclassical bulk description of correlation functions in the dual CFT involving the corresponding defect operators, which is nonperturbatively exact in α'. The Lorentzian continuation of these geometries describes a collection of defects colliding to make a BTZ black hole. We comment on a recent proposal to use such correlators to prepare a basis of black hole microstates, and elaborate on a picture of black hole formation and evaporation in terms of the underlying brane dynamics in the bulk.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Towards high-precision parton distributions from lattice QCD via distillation

We apply the Distillation spatial smearing program to the extraction of the unpolarized isovector valence PDF of the nucleon. The improved volume sampling and control of excited-states afforded by distillation leads to a dramatically improved determination of the requisite Ioffe-time Pseudo-distribution (pITD). The impact of higher-twist effects is subsequently explored by extending the Wilson line length present in our non-local operators to one half the spatial extent of the lattice ensemble considered. The valence PDF is extracted by analyzing both the matched Ioffe-time Distribution (ITD), as well as a direct matching of the pITD to the PDF. Through development of a novel prescription to obtain the PDF from the pITD, we establish a concerning deviation of the pITD from the expected DGLAP evolution of the pseudo-PDF. The presence of DGLAP evolution is observed once more following introduction of a discretization term into the PDF extractions. Observance and correction of this discrepancy further highlights the utility of distillation in such structure studies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Graviton topology

Over the past three decades, it has been shown that discrete and continuous media can support topologically nontrivial waves. Recently, it was shown that the same is true of the vacuum, in particular, right (R) and left (L) circularly polarized photons are topologically nontrivial. Here, we study the topology of another class of massless particles, namely gravitons. We show that the collection of all gravitons forms a topologically trivial vector bundle over the lightcone, allowing us to construct a globally smooth basis for gravitons. The graviton bundle also has a natural geometric splitting into two topologically nontrivial subbundles, consisting of the R and L gravitons. The R and L gravitons are unitary irreducible bundle representations of the Poincaré group, and are thus elementary particles; their topology is characterized by the Chern numbers ∓4. This nontrivial topology obstructs the splitting of graviton angular momentum into spin and orbital angular momentum.

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 OSTI2022

Matching Complexes of Trees and Applications of the Matching Tree Algorithm

A matching complex of a simple graph G is a simplicial complex with faces given by the matchings of G. The topology of matching complexes is mysterious; there are few graphs for which the homotopy type is known. Marietti and Testa showed that matching complexes of forests are contractible or homotopy equivalent to a wedge of spheres. We study two specific families of trees. For caterpillar graphs, we give explicit formulas for the number of spheres in each dimension and for perfect binary trees we find a strict connectivity bound. We also use a tool from discrete Morse theory called the Matching Tree Algorithm to study the connectivity of honeycomb graphs, partially answering a question raised by Jonsson.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Level-set topology optimization with PDE generated conformal meshes

This paper presents a level-set topology optimization approach that uses conformal meshes for the analysis of the displacement field. The structure’s boundary is represented by the iso-contour of a level-set field discretized on a fixed background design mesh. The conformal mesh is updated for each design iteration via a PDE based mesh morphing process that identifies the set of facets in the background mesh that are homeomorphic to the boundary and relaxes the homeomorphic mesh to conform to the structure’s boundary and ensure high element quality. The conformal mesh allows for a more accurate computation of the response versus density and some level-set based methods which interpolate material properties using the volume fraction. Numerical examples illustrate the proposed approach by optimizing linear-elastic two- and three-dimensional structures, wherein insight into the performance of the mesh morphing process is provided. The examples also highlight the scalability of the approach.

42 ENGINEERING↗
DOE OSTI2022

Modeling rarefied gas chemistry with QuiPS, a novel quasi-particle method

The goal of this work is to build up the capability of quasi-particle simulation (QuiPS), a novel flow solver, such that it can adequately model the rarefied portion of an atmospheric reentry trajectory. Direct simulation Monte Carlo (DSMC) is the conventional solver for such conditions, but struggles to resolve transient flows, trace species, and high-level internal energy states due to stochastic noise. Quasi-particle simulation (QuiPS) is a novel Boltzmann solver that describes a system with a discretized, truncated velocity distribution function. The resulting fixed-velocity, variable weight quasi-particles enable smooth variation of macroscopic properties. The distribution function description enables the use of a variance-reduced collision model, greatly minimizing expense near equilibrium. This work presents the addition of a neutral air chemistry model to QuiPS and some demonstrative 0D simulations. The explicit representation of internal distributions in QuiPS reveals some of the flaws in existing physics models. Overall, variance reduction, a key feature of QuiPS, can greatly reduce expense of multi-dimensional calculations, but is only cheaper when the gas composition is near chemical equilibrium.

42 ENGINEERING↗
DOE OSTI2021

Conservative numerical schemes with optimal dispersive wave relations: Part I. Derivation and analysis

An energy-conserving and an energy-and-enstrophy conserving numerical schemes are derived by approximating the Hamiltonian formulation of the inviscid shallow water flows based on the vorticity-divergence variables. These schemes also conserve the first-order moments such as mass and vorticity, as usual. The conservative properties of the schemes stem from the skew-symmetry and singularities of the Poisson brackets, which are carefully retained in the discrete approximations. Here, the schemes operate on unstructured orthogonal dual meshes, over bounded or unbounded domains, and they are also shown to possess the same optimal dispersive wave relations as those of the Z-grid scheme, which is a consequence of the use of the vorticity and divergence variables.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Narain to Narnia

We generalize the holographic correspondence between topological gravity coupled to an abelian Chern–Simons theory in three dimensions and an ensemble average of Narain’s family of massless free bosons in two dimensions, discovered by Afkhami-Jeddi et al. and by Maloney and Witten. We find that the correspondence also works for toroidal orbifolds but not for K3 or Calabi–Yau sigma-models and not always for the minimal models. Here, we conjecture that the correspondence requires that the central charge is equal to the critical central charge defined by the asymptotic density of states of the chiral algebra. For toroidal orbifolds, we extend the holographic correspondence to correlation functions of twist operators by using topological properties of rational tangles in the three-dimensional ball, which represent configurations of vortices associated to a discrete gauge symmetry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

On the Trotter Error in Many-body Quantum Dynamics with Coulomb Potentials

Efficient simulation of many-body quantum systems is central to advances in physics, chemistry, and quantum computing, with a key question being whether the simulation cost scales polynomially with the system size. Here, in this work, we analyze many-body quantum systems with Coulomb interactions, which are fundamental to electronic and molecular systems. We prove that Trotterization for such unbounded Hamiltonians achieves a 1/4-order convergence rate, with explicit polynomial dependence on the number of particles. The result holds for all initial wavefunctions in the domain of the Hamiltonian, and the 1/4-order convergence rate is optimal, as previous work has numerically demonstrated that it can be saturated by a specific initial ground state. The main challenges arise from the many-body structure and the singular nature of the Coulomb potential. Our proof strategy differs from prior state-of-the-art Trotter analyses, addressing both difficulties in a unified framework. Our analysis treats the Coulomb potential as an unbounded operator without modification or regularization, and does not rely on spatial discretization, making it compatible with both first- and second-quantized circuit constructions.

Fang, Di [Duke Univ., Durham, NC (United States)]↗
DOE OSTI2021

Spatially Quasi-Periodic Water Waves of Infinite Depth

Abstract We formulate the two-dimensional gravity-capillary water wave equations in a spatially quasi-periodic setting and present a numerical study of solutions of the initial value problem. We propose a Fourier pseudo-spectral discretization of the equations of motion in which one-dimensional quasi-periodic functions are represented by two-dimensional periodic functions on a torus. We adopt a conformal mapping formulation and employ a quasi-periodic version of the Hilbert transform to determine the normal velocity of the free surface. Two methods of time-stepping the initial value problem are proposed, an explicit Runge–Kutta (ERK) method and an exponential time-differencing (ETD) scheme. The ETD approach makes use of the small-scale decomposition to eliminate stiffness due to surface tension. We perform a convergence study to compare the accuracy and efficiency of the methods on a traveling wave test problem. We also present an example of a periodic wave profile containing vertical tangent lines that is set in motion with a quasi-periodic velocity potential. As time evolves, each wave peak evolves differently, and only some of them overturn. Beyond water waves, we argue that spatial quasi-periodicity is a natural setting to study the dynamics of linear and nonlinear waves, offering a third option to the usual modeling assumption that solutions either evolve on a periodic domain or decay at infinity.

Wilkening, Jon (ORCID:0000000327827596)↗
DOE OSTI2022

Local bilinear computation of Jacobi sets

Abstract We propose a novel method for the computation of Jacobi sets in 2D domains. The Jacobi set is a topological descriptor based on Morse theory that captures gradient alignments among multiple scalar fields, which is useful for multi-field visualization. Previous Jacobi set computations use piecewise linear approximations on triangulations that result in discretization artifacts like zig-zag patterns. In this paper, we utilize a local bilinear method to obtain a more precise approximation of Jacobi sets by preserving the topology and improving the geometry. Consequently, zig-zag patterns on edges are avoided, resulting in a smoother Jacobi set representation. Our experiments show a better convergence with increasing resolution compared to the piecewise linear method. We utilize this advantage with an efficient local subdivision scheme. Finally, our approach is evaluated qualitatively and quantitatively in comparison with previous methods for different mesh resolutions and across a number of synthetic and real-world examples.

97 MATHEMATICS AND COMPUTING↗