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

Sparse matrix‐vector and matrix‐multivector products for the truncated SVD on graphics processors

Summary Many practical algorithms for numerical rank computations implement an iterative procedure that involves repeated multiplications of a vector, or a collection of vectors, with both a sparse matrix and its transpose. Unfortunately, the realization of these sparse products on current high performance libraries often deliver much lower arithmetic throughput when the matrix involved in the product is transposed. In this work, we propose a hybrid sparse matrix layout, named CSRC, that combines the flexibility of some well‐known sparse formats to offer a number of appealing properties: (1) CSRC can be obtained at low cost from the popular CSR (compressed sparse row) format; (2) CSRC has similar storage requirements as CSR; and especially, (3) the implementation of the sparse product kernels delivers high performance for both the direct product and its transposed variant on modern graphics accelerators thanks to a significant reduction of atomic operations compared to a conventional implementation based on CSR. This solution thus renders considerably higher performance when integrated into an iterative algorithm for the truncated singular value decomposition (SVD), such as the randomized SVD or, as demonstrated in the experimental results, the block Golub–Kahan–Lanczos algorithm.

Aliaga, José I.↗
DOE OSTI2025

BCARS Simulated Phantom Dataset for Evaluation of Processing Pipelines

Broadband coherent anti-Stokes Raman scattering (BCARS) microscopy is a powerful label-free biological imaging technique, but the raw signal requires careful processing. The vibrationally resonant (Raman) fingerprint signal is usually small compared with instrumental noise sources and the nonresonant background (NRB) inherent in the BCARS signal. Fortunately, the NRB exhibits a systematic phase relationship with the coherent Raman response, acting as a heterodyne amplifier for the weak fingerprint signal. Due to this heterodyne effect, the Raman response can be recovered quantitatively and invariantly across different instruments, provided the NRB shape is known. Even with heterodyne amplification, the amplitudes of fingerprint signal components are often comparable to system noise. Singular value decomposition (SVD), which utilizes spatial information, is often employed for additional noise filtering. Consequently, finding optimal processing parameters to properly distinguish the NRB and Raman responses and suppress noise in the complex BCARS signal requires a reference system that realistically represents the spectral and spatial properties of BCARS signals obtained from biological samples. We present a digital tissue phantom that meets these criteria as a tool for testing candidate signal processing pipelines. The digital phantom is generated with simulated hyperspectral Raman images having system-specific noise and background characteristics. Here, we analyze phantom datasets with differing background and signal-to-noise conditions to evaluate their impact on the performance of multiple signal processing pipelines. Specifically, we investigate the application of a Butterworth filter-based routine to directly estimate the NRB from the BCARS signal. Additionally, we evaluate a Lorentzian wavelet transform as an alternative to the Hilbert transform for extracting the Raman spectrum from the BCARS signal. While we demonstrate this phantom for BCARS, it can be used for any spectroscopic Raman imaging approach.

Dixon, Jessica Z. [Georgia Institute of Technology↗
DOE OSTI2023

CP decomposition for tensors via alternating least squares with QR decomposition

The CP tensor decomposition is used in applications such as machine learning and signal processing to discover latent low-rank structure in multidimensional data. Computing a CP decomposition via an alternating least squares (ALS) method reduces the problem to several linear least squares problems. The standard way to solve these linear least squares subproblems is to use the normal equations, which inherit special tensor structure that can be exploited for computational efficiency. However, the normal equations are sensitive to numerical ill-conditioning, which can compromise the results of the decomposition. In this paper, we develop versions of the CP-ALS algorithm using the QR decomposition and the singular value decomposition, which are more numerically stable than the normal equations, to solve the linear least squares problems. Our algorithms utilize the tensor structure of the CP-ALS subproblems efficiently, have the same complexity as the standard CP-ALS algorithm when the input is dense and the rank is small, and are shown via examples to produce more stable results when ill-conditioning is present. Our MATLAB implementation achieves the same running time as the standard algorithm for small ranks, and we show that the new methods can obtain lower approximation error.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

High‐order cubature rules for tetrahedra

Summary In this paper, we construct new high‐order numerical integration schemes for tetrahedra, with positive weights and integration points that are in the interior of the domain. The construction of cubature rules is a challenging problem, which requires the solution of strongly nonlinear algebraic (moment) equations with side conditions given by affine inequality constraints. We present a robust algorithm based on a sequence of three modified Newton procedures to solve the constrained minimization problem. In the literature, numerical integration rules for the tetrahedron are available up to order p =15. We obtain integration rules for the tetrahedron from p =2 to p =20, which are computed using multiprecision arithmetic. For p ≤15, our approach provides integration rules that have the same or fewer number of integration points than existing rules; for p =16 to p =20, our rules are new. Numerical tests are presented that verify the polynomial‐precision of the cubature rules. Convergence studies are performed for the integration of exponential, rational, weakly singular and trigonometric test functions over tetrahedra with flat and curved faces. In all tests, improvements in accuracy is realized as p is increased, though in some cases nonmonotonic convergence is observed.

Jaśkowiec, Jan↗
DOE OSTI2020

High‐order symmetric cubature rules for tetrahedra and pyramids

Abstract In this article, we present an algorithm to construct high‐order fully symmetric cubature rules for tetrahedral and pyramidal elements, with positive weights and integration points that are in the interior of the domain. Cubature rules are fully symmetric if they are invariant to affine transformations of the domain. We divide the integration points into symmetry orbits where each orbit contains all the points generated by the permutation stars. These relations are represented by equality constraints. The construction of symmetric cubature rules require the solution of nonlinear polynomial equations with both inequality and equality constraints. For higher orders, we use an algorithm that consists of five sequential phases to produce the cubature rules. In the literature, symmetric numerical integration rules are available for the tetrahedron for orders p = 1 – 10, 14 , and for the pyramid up to p = 10 . We have obtained fully symmetric cubature rules for both of these elements up to order p = 20 . Numerical tests are presented that verify the polynomial‐precision of the cubature rules. Convergence studies are performed for the integration of exponential, weakly singular, and trigonometric test functions over both elements with flat and curved faces. With increase in p , improvements in accuracy is realized, though nonmonotonic convergence is observed.

Jaśkowiec, Jan↗
DOE OSTI2024

Numerical integration in the virtual element method with the scaled boundary cubature scheme

Abstract The virtual element method (VEM) is a stabilized Galerkin method on meshes that consist of arbitrary (convex and nonconvex) polygonal and polyhedral elements. A crucial ingredient in the implementation of low‐ and high‐order VEM is the numerical integration of monomials and nonpolynomial functions over such elements. In this article, we apply the recently proposed scaled boundary cubature (SBC) scheme to compute the weak form integrals in various virtual element formulations over polygonal and polyhedral meshes. In doing so, we demonstrate the flexibility of the approach and the accuracy that it delivers on a broad suite of boundary‐value problems in 2D and 3D over polytopes with affine faces as well as on elements with curved boundaries. In addition, the use of the SBC scheme is exemplified in an enriched Poisson formulation of the VEM in which weakly singular functions are required to be integrated. This study establishes the SBC method as a simple, accurate and efficient integration scheme for use in the VEM.

Chin, Eric B.↗
DOE OSTI2024

Design of a symmetry‐broken tetrahedral protein cage by a method of internal steric occlusion

Abstract Methods in protein design have made it possible to create large and complex, self‐assembling protein cages with diverse applications. These have largely been based on highly symmetric forms exemplified by the Platonic solids. Prospective applications of protein cages would be expanded by strategies for breaking the designed symmetry, for example, so that only one or a few (instead of many) copies of an exterior domain or motif might be displayed on their surfaces. Here we demonstrate a straightforward design approach for creating symmetry‐broken protein cages able to display singular copies of outward‐facing domains. We modify the subunit of an otherwise symmetric protein cage through fusion to a small inward‐facing domain, only one copy of which can be accommodated in the cage interior. Using biochemical methods and native mass spectrometry, we show that co‐expression of the original subunit and the modified subunit, which is further fused to an outward‐facing anti‐GFP DARPin domain, leads to self‐assembly of a protein cage presenting just one copy of the DARPin protein on its exterior. This strategy of designed occlusion provides a facile route for creating new types of protein cages with unique properties.

Gladkov, Nika↗
DOE OSTI2022

The Branes Behind Generalized Symmetry Operators

The modern approach to m-form global symmetries in a d-dimensional quantum field theory (QFT) entails specifying dimension d–m–1 topological generalized symmetry operators which non-trivially link with m-dimensional defect operators. In QFTs engineered via string constructions on a non-compact geometry X, these defects descend from branes wrapped on non-compact cycles which extend from a localized source / singularity to the boundary ∂X. The generalized symmetry operators which link with these defects arise from magnetic dual branes wrapped on cycles in ∂X. This provides a systematic way to read off various properties of such topological operators, including their worldvolume topological field theories, and the resulting fusion rules. We illustrate these general features in the context of 6D superconformal field theories, where we use the F-theory realization of these theories to read off the worldvolume theory on the generalized symmetry operators. Defects of dimension 3 which are charged under a suitable 3-form symmetry detect a non-invertible fusion rule for these operators. We also sketch how similar considerations hold for related systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Mapping resonant inelastic x-ray scattering onto electronic structure of iridium compounds

Resonant inelastic X-ray scattering (RIXS) is an indispensable tool that can selectively probe the electronic structure of active sites in energy conversion systems, for example, iridium oxides. However, decoding the RIXS spectra remains challenging due to its inherently complex many-body interactions. Here, we analyze Ir L 3 edge RIXS spectra of Ir, IrCl 3 , and IrO 2 , employing the joint density of states (JDOS) calculated directly from electronic band-structure calculations. The overall energy-loss features in the RIXS spectra are well reproduced by the JDOS, reaffirming the close correspondence between electron–hole excitations and RIXS spectra observed in prior studies. Intriguingly, the RIXS spectra above ∼4 eV in metallic Ir and IrO 2 follow a power-law behavior, 𝐼 RIXS ~Δ −𝑝 , where Δ is the energy loss and p is the associated power-law exponent. This is consistent with edge-singularity behavior commonly found in resonant X-ray scattering from metallic samples. Furthermore, orbital-projected JDOS enables a decomposition of the IrO 2 spectrum into specific dd transitions, providing a clear interpretation of orbital excitations and an efficient strategy for decoding RIXS spectra in iridium-based energy conversion systems.

5d↗
DOE OSTI2020

Learning compact physics‐aware delayed photocurrent models using dynamic mode decomposition

Abstract Radiation‐induced photocurrent in semiconductor devices can be simulated using complex physics‐based models, which are accurate, but computationally expensive. This presents a challenge for implementing device characteristics in high‐level circuit simulations where it is computationally infeasible to evaluate detailed models for multiple individual circuit elements. In this work we demonstrate a procedure for learning compact delayed photocurrent models that are efficient enough to implement in large‐scale circuit simulations, but remain faithful to the underlying physics. Our approach utilizes dynamic mode decomposition (DMD), a system identification technique for learning reduced‐order discrete‐time dynamical systems from time series data based on singular value decomposition. To obtain physics‐aware device models, we simulate the excess carrier density induced by radiation pulses by solving numerically the ambipolar diffusion equation, then use the simulated internal state as training data for the DMD algorithm. Our results show that the significantly reduced‐order delayed photocurrent models obtained via this method accurately approximate the dynamics of the internal excess carrier density—which can be used to calculate the induced current at the device boundaries—while remaining compact enough to incorporate into larger circuit simulations.

Hanson, Joshua↗
DOE OSTI2025

Randomized Algorithms for Low-Rank Matrix and Tensor Decompositions

This paper surveys randomized algorithms in numerical linear algebra for low-rank decompositions of matrices and tensors. The survey begins with a review of classical matrix algorithms that can be accelerated by randomized dimensionality reduction, such as the singular value decomposition (SVD) or interpolative (ID) and CUR decompositions. Recent advances in randomized dimensionality reduction are discussed, including new methods of fast matrix sketching and sampling techniques, which are incorporated into classical matrix algorithms for fast low-rank matrix approximations. The extension of randomized matrix algorithms to tensors is then explored for several low-rank tensor decompositions in the CP and Tucker formats, including the higher-order SVD, ID, and CUR decomposition.

Pearce, Katherine J. [The University of Texas at A↗
DOE OSTI2022

Celestial operator product expansions and w 1+∞ symmetry for all spins

The operator product expansion of massless celestial primary operators of arbitrary spin is investigated. Poincaré symmetry is found to imply a set of recursion relations on the operator product expansion coefficients of the leading singular terms at tree-level in a holomorphic limit. The symmetry constraints are solved by an Euler beta function with arguments that depend simply on the right-moving conformal weights of the operators in the product. These symmetry-derived coefficients are shown not only to match precisely those arising from momentum-space tree-level collinear limits, but also to obey an infinite number of additional symmetry transformations that respect the algebra of w 1+∞ . In tree-level minimally-coupled gravitational theories, celestial currents are constructed from light transforms of conformally soft gravitons and found to generate the action of w 1+∞ on arbitrary massless celestial primaries. Results include operator product expansion coefficients for fermions as well as those arising from higher-derivative non-minimal couplings of gluons and gravitons.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Renormalization group improved photon impact factors and the high energy virtual photon scattering

We perform the renormalization group improved collinear resummation of the photon-gluon impact factors. We construct the resummed cross section for virtual photon-photon (γ ∗ γ ∗ ) scattering which incorporates the impact factors and BFKL gluon Green’s function up to the next-to-leading logarithmic accuracy in energy. The impact factors include important kinematical effects which are responsible for the most singular poles in Mellin space at next-to-leading order. Further conditions on the resummed cross section are obtained by requiring the consistency with the collinear limits. Our analysis is consistent with previous impact factor calculations at NLO, apart from a new term proportional to C F that we find for the longitudinal polarization. Finally, we use the resummed cross section to compare with the LEP data on the γ ∗ γ ∗ cross section and with previous calculations. The resummed result is lower than the leading logarithmic approximation but higher than the pure next-to-leading one, and is consistent with the experimental data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Gluon transverse-momentum-dependent distributions from large-momentum effective theory

We demonstrate that gluon transverse-momentum-dependent parton distribution functions (TMDPDFs) can be extracted from lattice calculations of appropriate Euclidean correlations in large-momentum effective theory (LaMET). Based on perturbative calculations of gluon unpolarized and helicity TMDPDFs, we present a matching formula connecting them and their LaMET counterparts, where the latter are renormalized in a scheme facilitating lattice calculations and converted to the $\overline{MS}$ scheme. The hard matching kernel is given up to one-loop level. We also show that the perturbative result is independent of the prescription used for the pinch-pole singularity in the relevant correlations. Our results offer a guidance for the extraction of gluon TMDPDFs from lattice simulations, and have the potential to greatly facilitate perturbative calculations of the hard matching kernel.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Integrable field theories and their CCFT duals

We compute the Mellin transforms of various two-dimensional integrable S-matrices, providing the first explicit, non-perturbative realizations of celestial CFT. In two dimensions, the Mellin transform is simply the Fourier transform in rapidity space, and the “celestial correlator” has no position dependence. The simplified setting allows us to study the analytic properties of CCFT correlators exactly as a function of the conformal dimensions. We find that the correlators exist as real distributions of the conformal weights, with asymptotics controlled by the mass spectrum and three-point couplings of the model. Coupling these models to a flat space limit of JT gravity preserves integrability and dresses the amplitudes by a rapidly varying gravitational phase. We find that the coupling to gravity smooths out certain singular aspects of the Mellin-transformed correlators.

2D Gravity↗
DOE OSTI2023

The reflected entanglement spectrum for free fermions

We consider the reflected entropy and the associated entanglement spectrum for free fermions reduced to two intervals in 1 + 1 dimensions. Working directly in the continuum theory the reflected entropy can be extracted from the spectrum of a singular integral equation whose kernel is determined by the known free fermion modular evolved correlation function. We find the spectrum numerically and analytically in certain limits. For intervals that almost touch the reflected entanglement spectrum approaches the spectrum of the thermal density matrix. This suggests that the reflected entanglement spectrum is well suited to the task of extracting physical data of the theory directly from the ground state wave function.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Thermal energy of a charm-meson molecule in a pion gas

The thermal corrections to the propagator of a loosely bound charm-meson molecule in a pion gas are calculated to next-to-leading order in the heavy-meson expansion using a zero-range effective field theory. Ultraviolet divergences in the charm-meson-pair self energy are canceled by corrections to the charm-meson-pair contact vertex. Terms that are singular at the charm-meson-pair threshold can be absorbed into thermal corrections to the rest energies and kinetic masses of the charm-meson constituents. The remaining terms reduce to a thermal correction to the binding momentum that is proportional to the pion number density and suppressed by the pion/charm-meson mass ratio. The correction gives a tiny decrease in the binding energy of the charm-meson molecule relative to the charm-meson-pair threshold in the pion gas and a change in its thermal width that is small compared to the thermal widths of the charm-meson constituents. These results are encouraging for the prospects of observing X(3872) and $T$$^{+}_{cc}$ (3875) in the expanding hadron gas produced by heavy-ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Fibers add flavor. Part II. 5d SCFTs, gauge theories, and dualities

In [1, 2] we proposed an approach based on graphs to characterize 5d superconformal field theories (SCFTs), which arise as compactifications of 6d $$ \mathcal{N} $$ N = (1 , 0) SCFTs. The graphs, so-called combined fiber diagrams (CFDs), are derived using the realization of 5d SCFTs via M-theory on a non-compact Calabi-Yau threefold with a canonical singularity. In this paper we complement this geometric approach by connecting the CFD of an SCFT to its weakly coupled gauge theory or quiver descriptions and demonstrate that the CFD as recovered from the gauge theory approach is consistent with that as determined by geometry. To each quiver description we also associate a graph, and the embedding of this graph into the CFD that is associated to an SCFT provides a systematic way to enumerate all possible consistent weakly coupled gauge theory descriptions of this SCFT. Furthermore, different embeddings of gauge theory graphs into a fixed CFD can give rise to new UV-dualities for which we provide evidence through an analysis of the prepotential, and which, for some examples, we substantiate by constructing the M-theory geometry in which the dual quiver descriptions are manifest.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗