Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conformal algebra”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Perturbatively exact $w_{1+∞}$ asymptotic symmetry of quantum self-dual gravity

The infinite tower of positive-helicity soft gravitons in any minimally coupled, tree-level, asymptotically flat four-dimensional (4D) gravity was recently shown to generate a w 1+∞ asymptotic symmetry algebra. It is natural to ask whether this classical algebra acquires quantum corrections at loop level. We explore this in quantum self-dual gravity, whose amplitudes acquire known one-loop exact all-plus helicity quantum corrections. We show using collinear splitting formulae that, remarkably, the w 1+∞ algebra persists in quantum self-dual gravity without corrections.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Bootstrapping conformal QED 3 and deconfined quantum critical point

We bootstrap the deconfined quantum critical point (DQCP) and 3D Quantum Electrodynamics (QED 3 ) coupled to N f flavors of two-component Dirac fermions. We show the lattice and perturbative results on the SO(5) symmetric DQCP are excluded by the bootstrap bounds with an assumption that the lowest singlet scalar is irrelevant. Remarkably, we discover a new family of kinks in the 3D SO(N) vector bootstrap bounds with N ⩾ 6. We demonstrate coincidences between SU(N f ) adjoint and $SO (N^2_f-1$) vector bootstrap bounds due to a novel algebraic relation between the crossing equations. By introducing gap assumptions breaking the $SO (N^2_f-1$) symmetry, the SU(N f ) adjoint bootstrap bounds with large N f converge to the 1/N f perturbative results of QED 3 . Our results provide strong evidence that the SO(5) DQCP is not continuous and the critical flavor number of QED 3 is slightly above 2: $N^*_f\in$ (2,4). Bootstrap results near $N^*_f$ are well consistent with the merger and annihilation mechanism for the loss of conformality in QED 3 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Milestone 49 Report: Batched Sparse LA Phase 5 Implementation

Batched sparse linear algebra operations in general, and solvers in particular, have become the major algorithmic development activity and foremost performance engineering effort in the numerical software libraries work on modern hardware with accelerators such as GPUs. Many applications, ECP and non-ECP alike, require simultaneous solutions of many small linear systems of equations that are structurally sparse in one form or another. In order to move towards high hardware utilization levels, it is important to provide these applications with appropriate interface designs to be both functionally efficient and performance portable and give full access to the appropriate batched sparse solvers running on modern hardware accelerators prevalent across DOE supercomputing sites since the inception of ECP. To this end, we present here a summary of recent advances on the interface designs in use by HPC software libraries supporting batched sparse linear algebra and the development of sparse batched kernel codes for solvers and preconditioners. We also address the potential interoperability opportunities to keep the corresponding software portable between the major hardware accelerators from AMD, Intel, and NVIDIA, while maintaining the appropriate disclosure levels conforming to the active NDA agreements. The presented interface specifications include a mix of batched band, sparse iterative, and sparse direct solvers with their accompanying functionality that is already required by the application codes or we anticipated to be needed in the near future. This report summarizes progress in Kokkos Kernels and the xSDK libraries MAGMA, Ginkgo, hypre, PETSc, and SuperLU.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

Classical geometries defined by exterior differential systems on higher frame bundles

Exterior differential ideals are discussed, and sets of invariant generators presented, for Reimannian, conformal and projective geometries, and for specializations such as Ricci-flat, self-dual and Einstein-Maxwell theories. The Cartan characteristic integers are explicitly calculated, and involutory basis forms found, for each of these (specialized to four dimensions), exposing their algebraic structure and showing how they generate well-posed sets of partial differential equations.

Estabrook, Frank B.↗
NASA NTRS1983

A finite element method for diffusion dominated unsteady viscous flows

A general conforming finite element scheme for computing viscous flows is presented which is of second-order accuracy in space and time. Viscous terms are treated implicitly and advection terms are treated explicitly in the time marching segment of the algorithm. A method for solving the algebraic equations at each time step is given. The method is demonstrated on two test problems, one of them being a plane vortex flow for which asymptotic methods are used to obtain suitable numerical boundary conditions at each time step.

Gunzburger, M. D.↗
DOE OSTI2024

Celestial gluon amplitudes from the outside in

We show that, given a two-dimensional realization of the celestial OPE in self-dual Yang-Mills, we can find a scalar source around which scattering amplitudes replicate correlation functions computed from the 2D ‘gluon’ operators in a limit where a dynamic massless scalar decouples. We derive conditions on the two-dimensional three-point correlation function so that such a source exists and give two particular examples of this construction, one in which gluons are constructed from vertex operators in the semiclassical limit of Liouville theory and another in which the soft gluons arise from generalized free fields. Finally, we identify a bulk dual to the level of the boundary Kac-Moody algebra and discuss moving beyond the decoupling limit.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1988

Interactive grid generation for turbomachinery flow field simulations

The control point form of algebraic grid generation presented provides the means that are needed to generate well structured grids for turbomachinery flow simulations. It uses a sparse collection of control points distributed over the flow domain. The shape and position of coordinate curves can be adjusted from these control points while the grid conforms precisely to all boundaries. An interactive program called TURBO, which uses the control point form, is being developed. Basic features of the code are discussed and sample grids are presented. A finite volume LU implicit scheme is used to simulate flow in a turbine cascade on the grid generated by the program.

Choo, Yung K.↗
NASA NTRS1988

Interactive grid generation for turbomachinery flow field simulations

The control point form of algebraic grid generation presented provides the means that are needed to generate well structured grids of turbomachinery flow simulations. It uses a sparse collection of control points distributed over the flow domain. The shape and position of coordinate curves can be adjusted from these control points while the grid conforms precisely to all boundaries. An interactive program called TURBO, which uses the control point form, is being developed. Basic features of the code are discussed and sample grids are presented. A finite volume LU implicit scheme is used to simulate flow in a turbine cascade on the grid generated by the program.

Choo, Yung K.↗
DOE OSTI2026

Quantum simulation of massive Thirring and Gross--Neveu models for arbitrary number of flavors

The study of fermionic quantum field theories is an important problem for realizing the standard model of particle physics on a quantum computer. As a step towards this goal, we consider the massive Thirring and Gross--Neveu models with arbitrary number of fermion flavors, $N_f$, discretized on a spatial one-dimensional lattice of size $L$ in the Hamiltonian formulation. We compute the gate complexity using the higher-order product formula and using block-encoding/qubitization and quantum singular value transformations in the limit of large $N_f$ and $L$. We also prepare the ground states of both models with excellent fidelity for system sizes up to 20 qubits with $N_f = 1,2,3,4$ using the adaptive-variational quantum imaginary time algorithm. In addition, we also classify the dynamical Lie algebras of these relativistic fermionic models and show that they belong to the same isomorphism class. Our work is a concrete step towards the quantum simulation of real-time dynamics of large $N_f$ fermionic quantum field theories models relevant for chiral symmetry breaking, understanding dimensional transmutation, and exploring the conformal window of field theories on near-term and early fault-tolerant quantum computers.

FOS: Physical sciences↗
NASA NTRS1990

Algebraic grid generation with boundary orthogonality control

This paper presents a new approach in applying the blending function method to the generation of structured finite-difference meshes. It is shown that the surface equation developed by Coons (1967) can be easily adapted to generate two-dimensional meshes, and its extension to hyperspace allows the construction of three-dimensional meshes. The robust method in its basic form blends among the boundary curves/surfaces conforming to their intrinsic slopes. When parameter-controlled boundary-slope corrections are introduced, it is possible to establish varying degrees of orthogonality along the boundaries. Additional controls are made possible by moving points in the parameter space, which affects the movement of points in the physical space. The combination of these controls provides a simple and yet powerful means for exploiting the blending-function method to satisfy the needs of practical CFD applications. Results from a validation program for two-dimensions are shown for a variety of geometries.

Luh, Raymond Ching-Chung↗
DOE OSTI2022

An adaptive scalable fully implicit algorithm based on stabilized finite element for reduced visco-resistive MHD

The magnetohydrodynamics (MHD) equations are continuum models used in the study of a wide range of plasma physics systems, including the evolution of complex plasma dynamics in tokamak disruptions. However, efficient numerical solution methods for MHD are extremely challenging due to disparate time and length scales, strong hyperbolic phenomena, and nonlinearity. Additionally, therefore the development of scalable, implicit MHD algorithms and high-resolution adaptive mesh refinement strategies is of considerable importance. In this work, we develop a high-order stabilized finite-element algorithm for the reduced visco-resistive MHD equations based on the MFEM finite element library (mfem.org). The scheme is fully implicit, solved with the Jacobian-free Newton-Krylov (JFNK) method with a physics-based preconditioning strategy. Our preconditioning strategy is a generalization of the physics-based preconditioning methods in Chacón et al. (2002) to adaptive, stabilized finite elements. Algebraic multigrid methods are used to invert sub-block operators to achieve scalability. A parallel adaptive mesh refinement scheme with dynamic load-balancing is implemented to efficiently resolve the multi-scale spatial features of the system. Our implementation uses the MFEM framework, which provides arbitrary-order polynomials and flexible adaptive conforming and non-conforming meshes capabilities. Results demonstrate the accuracy, efficiency, and scalability of the implicit scheme in the presence of large scale disparity. The potential of the AMR approach is demonstrated on an island coalescence problem in the high Lundquist-number regime (≥ 10 7 ) with the successful resolution of plasmoid instabilities and thin current sheets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1983

A fully vectorized numerical solution of the incompressible Navier-Stokes equations

A vectorizable algorithm is presented for the implicit finite difference solution of the incompressible Navier-Stokes equations in general curvilinear coordinates. The unsteady Reynolds averaged Navier-Stokes equations solved are in two dimension and non-conservative primitive variable form. A two-layer algebraic eddy viscosity turbulence model is used to incorporate the effects of turbulence. Two momentum equations and a Poisson pressure equation, which is obtained by taking the divergence of the momentum equations and satisfying the continuity equation, are solved simultaneously at each time step. An elliptic grid generation approach is used to generate a boundary conforming coordinate system about an airfoil. The governing equations are expressed in terms of the curvilinear coordinates and are solved on a uniform rectangular computational domain. A checkerboard SOR, which can effectively utilize the computer architectural concept of vector processing, is used for iterative solution of the governing equations.

Patel, N.↗
DOE OSTI2020

Unitarity methods in AdS/CFT

We develop a systematic unitarity method for loop-level AdS scattering amplitudes, dual to non-planar CFT correlators, from both bulk and boundary perspectives. We identify cut operators acting on bulk amplitudes that put virtual lines on shell, and show how the conformal partial wave decomposition of the amplitudes may be efficiently computed by gluing lower-loop amplitudes. A central role is played by the double discontinuity of the amplitude, which has a direct relation to these cuts. We then exhibit a precise, intuitive map between the diagrammatic approach in the bulk using cutting and gluing, and the algebraic, holographic unitarity method of [1] that constructs the non-planar correlator from planar CFT data. Our analysis focuses mostly on four-point, one-loop diagrams — we compute cuts of the scalar bubble, triangle and box, as well as some one-particle reducible diagrams — in addition to the five-point tree and four-point double-ladder. Analogies with S-matrix unitarity methods are drawn throughout.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A Condensed Constrained Nonconforming Mortar-Based Approach for Preconditioning Finite Element Discretization Problems

This paper presents and studies an approach for constructing auxiliary space preconditioners for finite element problems using a constrained nonconforming reformulation that is based on a proposed modified version of the mortar method. The well-known mortar finite element discretization method is modified to admit a local structure, providing an element-by-element or subdomain-by-subdomain assembly property. This is achieved via the introduction of additional trace finite element spaces and degrees of freedom (unknowns) associated with the interfaces between adjacent elements or subdomains. The resulting nonconforming formulation and a reduced-via-static-condensation Schur complement form on the interfaces are used in the construction of auxiliary space preconditioners for a given conforming finite element discretization problem. Overall, the properties of these preconditioners are studied and their performance is illustrated on model second order scalar elliptic problems utilizing high order elements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Auxiliary Space Preconditioning of Finite Element Equations Using a Nonconforming Interior Penalty Reformulation and Static Condensation

We modify the well-known interior penalty finite element discretization method so that it allows for element-by-element assembly. This is possible due to the introduction of additional unknowns associated with the interfaces between neighboring elements. The resulting bilinear form, and a Schur complement (reduced) version of it, are utilized in a number of auxiliary space preconditioners for the original conforming finite element discretization problem. Furthermore, these preconditioners are analyzed on the fine scale, and their performance is illustrated on model second order scalar elliptic problems discretized with high order elements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Supersymmetry and the celestial Jacobi identity

In this paper we study the simplifying effects of supersymmetry on celestial OPEs at both tree and loop level. We find at tree level that theories with unbroken supersymmetry around a stable vacuum have celestial soft current algebras satisfying the Jacobi identity, and we show at one loop that celestial OPEs in these theories have no double poles.

Supersymmetric Gauge Theory↗
DOE OSTI2020

Scattering equations in AdS: scalar correlators in arbitrary dimensions

We introduce a bosonic ambitwistor string theory in AdS space. Even though the theory is anomalous at the quantum level, one can nevertheless use it in the classical limit to derive a novel formula for correlation functions of boundary CFT operators in arbitrary space-time dimensions. The resulting construction can be treated as a natural extension of the CHY formalism for the flat-space S-matrix, as it similarly expresses tree-level amplitudes in AdS as integrals over the moduli space of Riemann spheres with punctures. These integrals localize on an operator-valued version of scattering equations, which we derive directly from the ambitwistor string action on a coset manifold. As a testing ground for this formalism we focus on the simplest case of ambitwistor string coupled to two cur- rent algebras, which gives bi-adjoint scalar correlators in AdS. In order to evaluate them directly, we make use of a series of contour deformations on the moduli space of punctured Riemann spheres and check that the result agrees with tree level Witten diagram computations to all multiplicity. We also initiate the study of eigenfunctions of scattering equations in AdS, which interpolate between conformal partial waves in different OPE channels, and point out a connection to an elliptic deformation of the Calogero-Sutherland model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

CFT correlators, $\mathcal{W}$-algebras and generalized Catalan numbers

In two spacetime dimensions the Virasoro heavy-heavy-light-light (HHLL) vacuum block in a certain limit is governed by the Catalan numbers. The equation for their generating function can be generalized to a differential equation which the logarithm of the block satisfies. We show that a similar story holds for the HHLL $\mathcal{W}_N$ vacuum blocks, where a suitable generalization of the Catalan numbers plays the main role. Moreover, the $\mathcal{W}_N$ blocks have the same form as the stress tensor sector of HHLL near lightcone conformal correlators in 2(N – 1) spacetime dimensions. In the latter case the Catalan numbers are generalized to the numbers of linear extensions of certain partially ordered sets.

1/N expansion↗