Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “abstract 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

LaRIS: Targeting Portability and Productivity for LAPACK Codes on Extreme Heterogeneous Systems by Using IRIS

In keeping with the trend of heterogeneity in high-performance computing, hardware manufacturers and vendors are developing new architectures and associated software stacks (e.g., libraries) to harness the best possible performance from commonly used kernels (e.g., linear algebra kernels). However, kernels tuned for one architecture are not portable to others. Moreover, the coexistence of different architectures in a single node makes orchestration difficult. To address these challenges, we introduce LaRIS, a portable framework for LAPACK functionalities. LaRIS ensures a separation between linear algebra algorithms and vendor-library kernels by using the IRIS run time and IRIS-BLAS library. Such abstraction at the algorithm level makes the implementation completely agnostic to the vendor library and architecture. LaRIS uses the IRIS run time to dynamically select the vendor-library kernel and suitable processor architecture at run time. Through LU factorization, we demonstrate that LaRIS can fully utilize different heterogeneous systems by launching and orchestrating different vendor-library kernels without any change in the source code.

Monil, M. A. H.↗
DOE OSTI2024

Effective rationality for local unitary invariants of mixed states of two qubits

Abstract We calculate the field of rational local unitary invariants for mixed states of two qubits, by employing methods from algebraic geometry. We prove that this field is rational (i.e. purely transcendental), and that it is generated by nine algebraically independent polynomial invariants. We do so by constructing a relative section, in the sense of invariant theory, whose Weyl group is a finite abelian group. From this construction, we are able to give explicit expressions for the generating invariants in terms of the Bloch matrix representation of mixed states of two qubits. We also prove similar rationality results for the local unitary invariants of symmetrically mixed states of two qubits. We also provide a sketch of how to generalize our results to the case of an arbitrary number of qubits. Our results apply to both complex-valued and real-valued invariants.

Physics↗
DOE OSTI2024

3-Manifolds and VOA Characters

Abstract By studying the properties ofq-series$$\widehat{Z}$$ Z ^ -invariants, we develop a dictionary between 3-manifolds and vertex algebras. In particular, we generalize previously known entries in this dictionary to Lie groups of higher rank, to 3-manifolds with toral boundaries, and to BPS partition functions with line operators. This provides a new physical realization of logarithmic vertex algebras in the framework of the 3d-3d correspondence and opens new avenues for their future study. For example, we illustrate how invoking a knot-quiver correspondence for$$\widehat{Z}$$ Z ^ -invariants leads to many infinite families of new fermionic formulae for VOA characters.

Physics↗
DOE OSTI2022

Deep Learning for Subgrid‐Scale Turbulence Modeling in Large‐Eddy Simulations of the Convective Atmospheric Boundary Layer

Abstract In large‐eddy simulations, subgrid‐scale (SGS) processes are parameterized as a function of filtered grid‐scale variables. First‐order, algebraic SGS models are based on the eddy‐viscosity assumption, which does not always hold for turbulence. Here we apply supervised deep neural networks (DNNs) to learn SGS stresses from a set of neighboring coarse‐grained velocity from direct numerical simulations of the convective boundary layer at friction Reynolds numbers Re τ up to 1243 without invoking the eddy‐viscosity assumption. The DNN model was found to produce higher correlation between SGS stresses compared to the Smagorinsky model and the Smagorinsky‐Bardina mixed model in the surface and mixed layers and can be applied to different grid resolutions and various stability conditions ranging from near neutral to very unstable. The DNN model can capture key statistics of turbulence in a posteriori (online) tests when applied to large‐eddy simulations of the atmospheric boundary layer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A non-commutative Bayes' theorem

Using a diagrammatic reformulation of Bayes' theorem, we provide a necessary and sufficient condition for the existence of Bayesian inference in the setting of finite-dimensional C* -algebras. In other words, we prove an analogue of Bayes' theorem in the joint classical and quantum context. Our analogue is justified by recent advances in categorical probability theory, which have provided an abstract formulation of the classical Bayes' theorem. In the process, we further develop non-commutative almost everywhere equivalence and illustrate its important role in non-commutative Bayesian inversion. The construction of such Bayesian inverses, when they exist, involves solving a positive semidefinite matrix completion problem for the Choi matrix. This gives a solution to the open problem of constructing Bayesian inversion for completely positive unital maps acting on density matrices that do not have full support. In conclusion, we illustrate how the procedure works for several examples relevant to quantum information theory.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Pyomo.DOE: An open-source package for model-based design of experiments in Python

Predictive mathematical models are a cornerstone of science and engineering. Yet selecting, calibrating, and validating said science-based models often remains an art in practice. Model-based design of experiments (MBDoE) provides a systematic framework to maximize information gain from experiments while minimizing time and resource costs. But MBDoE remains limited to niche application areas, in part because practitioners must integrate expertise in statistics, computational optimization, and modeling. To help reduce this barrier, we introduce Pyomo.DOE, an open-source package for MBDoE. Pyomo.DOE uses a nonlinear sensitivity analysis code k_aug to quickly approximate the Fisher information matrix and leverages a new stochastic programming abstraction. We demonstrate Pyomo.DOE with the first application of MBDoE to fixed-bed breakthrough experiments, which highlights the power of Pyomo.DOE to quantify the value of experimental modifications a priori for large-scale partial differential-algebraic equation (PDAE) models. Here we also provide a mathematical primer on MBDoE targeted at general chemical engineers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Modular Finite Element Methods

MFEM is a free, lightweight, scalable C++ library for finite element methods. The goal of MFEM is to enable research and development of scalable finite element discretization and solver algorithms through general finite element abstractions, accurate and flexible visualization, and tight integration with the hypre library. Conceptually, MFEM can be viewed as a finite element toolbox that provides the building blocks for developing finite element algorithms in a manner similar to that of MATLAB for linear algebra methods.

ECP↗
DOE OSTI2022

Causal diamonds, cluster polytopes and scattering amplitudes

The “amplituhedron” for tree-level scattering amplitudes in the bi-adjoint φ 3 theory is given by the ABHY associahedron in kinematic space, which has been generalized to give a realization for all finite-type cluster algebra polytopes, labelled by Dynkin diagrams. In this letter we identify a simple physical origin for these polytopes, associated with an interesting (1 + 1)-dimensional causal structure in kinematic space, along with solutions to the wave equation in this kinematic “spacetime” with a natural positivity property. The notion of time evolution in this kinematic spacetime can be abstracted away to a certain “walk”, associated with any acyclic quiver, remarkably yielding a finite cluster polytope for the case of Dynkin quivers. The A n–3 , B n–1 /C n–1 and D n polytopes are the amplituhedra for n-point tree amplitudes, one-loop tadpole diagrams, and full integrand of one-loop amplitudes. We also introduce a polytope D¯ n , which chops the D n polytope in half along a symmetry plane, capturing one-loop amplitudes in a more efficient way.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

A T-duality of non-supersymmetric heterotic strings and an implication for Topological Modular Forms

Abstract Motivated by recent developments connecting non-supersymmetric heterotic string theory to the theory of Topological Modular Forms (TMF), we show that the worldsheet theory with central charge (17,$$ \frac{3}{2} $$ 3 2 ) obtained by fibering the (E 8 ) 1 × (E 8 ) 1 current algebra over the two$$ \mathcal{N} $$ N = (0, 1) sigma model onS 1 with antiperiodic spin structure (such that theE 8 factors are exchanged as we go around the circle), is continuously connected to the (E 8 ) 2 theory in the Gaiotto Johnson-Freyd Witten sense of going “up and down the RG trajectories”. Combined with the work of Tachikawa and Yamashita, this furnishes a physical derivation of the fact that the (E 8 ) 2 theory corresponds to the unique nontrivial torsion element [(E 8 ) 2 ] of TMF 31 with zero mod-2 elliptic genus.

Physics↗
DOE OSTI2024

A cluster of results on amplituhedron tiles

Abstract The amplituhedron is a mathematical object which was introduced to provide a geometric origin of scattering amplitudes in $$\mathcal {N}=4$$ N = 4 super Yang–Mills theory. It generalizes cyclic polytopes and the positive Grassmannian and has a very rich combinatorics with connections to cluster algebras. In this article, we provide a series of results about tiles and tilings of the $$m=4$$ m = 4 amplituhedron. Firstly, we provide a full characterization of facets of BCFW tiles in terms of cluster variables for $$\text{ Gr}_{4,n}$$ Gr 4 , n . Secondly, we exhibit a tiling of the $$m=4$$ m = 4 amplituhedron which involves a tile which does not come from the BCFW recurrence—the spurion tile, which also satisfies all cluster properties. Finally, strengthening the connection with cluster algebras, we show that each standard BCFW tile is the positive part of a cluster variety, which allows us to compute the canonical form of each such tile explicitly in terms of cluster variables for $$\text{ Gr}_{4,n}$$ Gr 4 , n . This paper is a companion to our previous paper “Cluster algebras and tilings for the $$m=4$$ m = 4 amplituhedron.”

Physics↗
DOE OSTI2023

Structure of pressure-gradient-driven current singularity in ideal magnetohydrodynamic equilibrium

Abstract Singular currents typically appear on rational surfaces in non-axisymmetric ideal magnetohydrodynamic (MHD) equilibria with a continuum of nested flux surfaces and a continuous rotational transition. These currents have two components: a surface current (Dirac δ -function in flux surface labeling) that prevents the formation of magnetic islands, and an algebraically divergent Pfirsch–Schlüter current density when a pressure gradient is present across the rational surface. On flux surfaces adjacent to the rational surface, the traditional treatment gives the Pfirsch–Schlüter current density scaling as J ∼ 1 / Δ ι , where Δ ι is the difference of the rotational transform relative to the rational surface. If the distance s between flux surfaces is proportional to Δ ι , the scaling relation J ∼ 1 / Δ ι ∼ 1 / s will lead to a paradox that the Pfirsch–Schlüter current is not integrable. In this work, we investigate this issue by considering the pressure-gradient-driven singular current in the Hahm–Kulsrud–Taylor problem, which is a prototype for singular currents arising from resonant magnetic perturbations. We show that not only the Pfirsch–Schlüter current density but also the diamagnetic current density are divergent as ∼ 1 / Δ ι . However, due to the formation of a Dirac δ -function current sheet at the rational surface, the neighboring flux surfaces are strongly packed with s ∼ ( Δ ι ) 2 . Consequently, the singular current density J ∼ 1 / s , making the total current finite, thus resolving the paradox. Furthermore, the strong packing of flux surfaces causes a steepening of the pressure gradient near the rational surface, with ∇ p ∼ d p / d s ∼ 1 / s . In general non-axisymmetric MHD equilibrium, contrary to Grad’s conjecture that the pressure profile is flat around densely distributed rational surfaces, our result suggests a pressure profile that densely steepens around them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Light-front puzzles

Abstract Light-front formulations of quantum field theories have many advantages for computing electroweak matrix elements of strongly interacting systems and other quantities that are used to study hadronic structure. The theory can be formulated in Hamiltonian form so non-perturbative calculations of the strongly interacting initial and final states are in principle reduced to linear algebra. These states are needed for calculating parton distribution functions and other types of distribution amplitudes that are used to understand the structure of hadrons. Light-front boosts are kinematic transformations so the strongly interacting states can be computed in any frame. This is useful for computing current matrix elements involving electroweak probes where the initial and final hadronic states are in different frames related by the momentum transferred by the probe. Finally in many calculations the vacuum is trivial so the calculations can be formulated in Fock space. The advantages of light front-field theory would not be interesting if the light-front formulation was not equivalent to the covariant or canonical formulations of quantum field theory. Many of the distinguishing properties of light-front quantum field theory are difficult to reconcile with canonical or covariant formulations of quantum field theory. This paper discusses the resolution of some of the apparent inconsistencies in canonical, covariant and light-front formulations of quantum field theory. The puzzles that will be discussed are (1) the problem of inequivalent representations (2) the problem of the trivial vacuum (3) the problem of ill-posed initial value problems (4) the problem of rotational covariance (5) the problem of zero modes and (6) the problem of spontaneously broken symmetries.

Physics↗
DOE OSTI2023

Monte Carlo goodness-of-fit tests for degree corrected and related stochastic blockmodels

Abstract We construct Bayesian and frequentist finite-sample goodness-of-fit tests for three different variants of the stochastic blockmodel for network data. Since all of the stochastic blockmodel variants are log-linear in form when block assignments are known, the tests for the latent block model versions combine a block membership estimator with the algebraic statistics machinery for testing goodness-of-fit in log-linear models. We describe Markov bases and marginal polytopes of the variants of the stochastic blockmodel and discuss how both facilitate the development of goodness-of-fit tests and understanding of model behaviour. The general testing methodology developed here extends to any finite mixture of log-linear models on discrete data, and as such is the first application of the algebraic statistics machinery for latent-variable models.

Karwa, Vishesh↗
DOE OSTI2021

QLiG: Query Like a Graph For Subgraph Matching

A graph is a natural and flexible modeling approach to represent entities and relationships between them in real-world. A Knowledge Graphs (KG) is a specialized graph with formal and structured representation of facts, relationships, annotated with semantic descriptions. Subgraph matching is one of the fundamental graph problems to identify relationships, interactions and activities of interest within a large graph. A query specification is a collection of abstract components, operations, and constraints to express a pattern. The specification can be implemented in different ways based on underlying data model. Various graph query specifications have been developed over the years and have led to the development of different open-sourced and vendor-specific query languages. Such specification are modeled as an extension of relational algebra used to develop relational query languages such as SQL. Such relational concepts do not inherently support graph queries. There is a need to represent graph queries in terms on graph-based components to expedite query construction by non-database experts. We present a graph-based query approach QLiG (pronounced cleeg), to perform subgraph matching in Labeled Property Graph. We present the query specifications, salient features, and a use case to show functional examples.

Purohit, Sumit↗
DOE OSTI2026

Phenomena-based graph representations and applications to chemical process simulation

Rapid and robust simulation of chemical processes is critical to conduct process design, optimization, techno-economic analysis, and sustainability analysis. Yet, efficiently solving simulation models remains a challenge due to the highly coupled and nonlinear nature of the underlying algebraic equations that capture the physical phenomena taking place in the process (e.g., material and energy conservation, phase equilibrium, reactions). In this work, we show that graph-theoretic representations of the physical phenomena within unit operations can help navigate and decompose equations to systematically identify alternative approaches for fast and robust numerical solutions. Specifically, we present a graph-theoretic abstraction that captures the connectivity between the model variables/equations and use this abstraction to group variables/equations into fundamental phenomena. We show that phenomena-based decomposition of the underlying equations can help decouple nonlinearities and enforce material/energy conservation at the process level to accelerate convergence. The proposed decomposition approach differs from the more traditional sequential modular simulation approach, in which equations are grouped and decomposed by unit operations. We implemented the phenomena-based decomposition in BioSTEAM—an open-source process simulation platform in Python—and demonstrated that this approach can converge a variety of separation process models. Compared to sequential modular simulation, the phenomena-based approach can converge idealized systems faster, but it can be slower for (or even fail to converge) highly coupled and nonideal process systems.

Convergence↗
DOE OSTI2024

From large to small $$ \mathcal{N} $$ = (4, 4) superconformal surface defects in holographic 6d SCFTs

Abstract Two-dimensional (2d)$$ \mathcal{N} $$ N = (4, 4) Lie superalgebras can be either “small” or “large”, meaning their R-symmetry is either$$ \mathfrak{so} $$ so (4) or$$ \mathfrak{so} $$ so (4) ⊕$$ \mathfrak{so} $$ so (4), respectively. Both cases admit a superconformal extension and fit into the one-parameter family$$ \mathfrak{d} $$ d (2,1;γ) ⊕$$ \mathfrak{d} $$ d (2,1;γ), with parameterγ∈ (−∞,∞). The large algebra corresponds to generic values ofγ, while the small case corresponds to a degeneration limit withγ→ −∞. In 11d supergravity, we study known solutions with superisometry algebra$$ \mathfrak{d} $$ d (2,1;γ) ⊕$$ \mathfrak{d} $$ d (2,1;γ) that are asymptotically locally AdS 7 ×𝕊 4 . These solutions are holographically dual to the 6d maximally superconformal field theory with 2d superconformal defects invariant under$$ \mathfrak{d} $$ d (2,1;γ) ⊕$$ \mathfrak{d} $$ d (2,1;γ). We show that a limit of these solutions, in whichγ→ −∞, reproduces another known class of solutions, holographically dual tosmall$$ \mathcal{N} $$ N = (4, 4) superconformal defects. We then use this limit to generate new small$$ \mathcal{N} $$ N = (4, 4) solutions with finite Ricci scalar, in contrast to the known small$$ \mathcal{N} $$ N = (4, 4) solutions. We then use holography to compute the entanglement entropy of a spherical region centered on these small$$ \mathcal{N} $$ N = (4, 4) defects, which provides a linear combination of defect Weyl anomaly coefficients that characterizes the number of defect-localized degrees of freedom. We also comment on the generalization of our results to include$$ \mathcal{N} $$ N = (0,4) surface defects through orbifolding.

Physics↗
DOE OSTI2024

Null Raychaudhuri: canonical structure and the dressing time

Abstract We initiate a study of gravity focusing on generic null hypersurfaces, non-perturbatively in the Newton coupling. We present an off-shell account of the extended phase space of the theory, which includes the expected spin-2 data as well as spin-0, spin-1 and arbitrary matter degrees of freedom. We construct the charges and the corresponding kinematic Poisson brackets, employing a Beltrami parameterization of the spin-2 modes. We explicitly show that the constraint algebra closes, the details of which depend on the non-perturbative mixing between spin-0 and spin-2 modes. Finally we show that the spin zero sector encodes a notion of a clock, called dressing time, which is dynamical and conjugate to the constraint. It is well-known that the null Raychaudhuri equation describes how the geometric data of a null hypersurface evolve in null time in response to gravitational radiation and external matter. Our analysis leads to three complementary viewpoints on this equation. First, it can be understood as a Carrollian stress tensor conservation equation. Second, we construct spin-0, spin-2 and matter stress tensors that act as generators of null time reparametrizations for each sector. This leads to the perspective that the null Raychaudhuri equation can be understood as imposing that the sum of CFT-like stress tensors vanishes. Third, we solve the Raychaudhuri constraint non-perturbatively. The solution relates the dressing time to the spin-2 and matter boost charge operators. Finally we establish that the corner charge corresponding to the boost operator in the dressing time frame is monotonic. These results show that the notion of an observer can be thought of as emerging from the gravitational degrees of freedom themselves. We briefly mention that the construction offers new insights into focusing conjectures.

Physics↗
DOE OSTI2021

The IDAES process modeling framework and model library—Flexibility for process simulation and optimization

Abstract Energy systems and manufacturing processes of the 21st century are becoming increasingly dynamic and interconnected, which require new capabilities to effectively model and optimize their design and operations. Such next generation computational tools must leverage state‐of‐the‐art techniques in optimization and be able to rapidly incorporate new advances. To address these requirements, we have developed the Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform, which builds on the strengths of both process simulators (model libraries) and algebraic modeling languages (advanced solvers). This paper specifically presents the IDAES Core Modeling Framework (IDAES‐CMF), along with a case study demonstrating the application of the framework to solve process optimization problems. Capabilities provided by this framework include a flexible, modifiable, open‐source platform for optimization of process flowsheets utilizing state‐of‐the‐art solvers and solution techniques, fully open and extensible libraries of dynamic unit operations models and thermophysical property models, and integrated support for superstructure‐based conceptual design and optimization under uncertainty.

42 ENGINEERING↗