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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

What happens to Petrov classification, on horizons of axisymmetric dirty black holes

We consider axisymmetric stationary dirty black holes with regular non-extremal or extremal horizons, and compute their on-horizon Petrov types. The Petrov type (PT) in the frame of the observer crossing the horizon can be different from that formally obtained in the usual (but singular in the horizon limit) frame of an observer on a circular orbit. We call this entity the boosted Petrov type (BPT), as the corresponding frame is obtained by a singular boost from the regular one. The PT off-horizon can be more general than PT on-horizon and that can be more general than the BPT on horizon. This is valid for all regular metrics, irrespective of the extremality of the horizon. We analyze and classify the possible relations between the three characteristics and discuss the nature and features of the underlying singular boost. The three Petrov types can be the same only for space-times of PT D and O off-horizon. The mutual alignment of principal null directions and the generator in the vicinity of the horizon is studied in detail. As an example, we also analyze a special class of metrics with utra-extremal horizons (for which the regularity conditions look different from the general case) and compare their off-horizon and on-horizon algebraic structure in both frames.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

The direct criterion of Newcomb for the ideal MHD stability of an axisymmetric toroidal plasma

A method is presented for determining the ideal magnetohydrodynamic stability of an axisymmetric toroidal plasma, based on a toroidal generalization of the method developed by Newcomb for fixed-boundary modes in a cylindrical plasma. For toroidal mode number n≠0, the stability problem is reduced to the numerical integration of a high-order complex system of ordinary differential equations, the Euler-Lagrange equation for extremizing the potential energy, for the coupled amplitudes of poloidal harmonics m as a function of the radial coordinate ψ in a straight-fieldline flux coordinate system. Unlike the cylindrical case, different poloidal harmonics couple to each other, which introduces coupling between adjacent singular intervals. A boundary condition is used at each singular surface, where m = nq and q(ψ) is the safety factor, to cross the singular surface and continue the solutions beyond it. Fixed-boundary instability is indicated by the vanishing of a real determinant of a Hermitian complex matrix constructed from the fundamental matrix of solutions, the generalization of Newcomb's crossing criterion. In the absence of fixed-boundary instabilities, an M × M plasma response matrix W P with M the number of poloidal harmonics used, is constructed from the Euler-Lagrange solutions at the plasma-vacuum boundary. This is added to a vacuum response matrix W V to form a total response matrix W T . Finally, the existence of negative eigenvalues of W T indicates the presence of free-boundary instabilities. The method is implemented in the fast and accurate DCON code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Deformation of Lie–Poisson algebras and chirality

Linearization of a Hamiltonian system around an equilibrium point yields a set of Hamiltonian symmetric spectra: If λ is an eigenvalue of the linearized generator, -λ and $\barλ$ (hence, -$\barλ$) are also eigenvalues—the former implies a time-reversal symmetry, while the latter guarantees the reality of the solution. However, linearization around a singular equilibrium point (which commonly exists in noncanonical Hamiltonian systems) works out differently, resulting in breaking of the Hamiltonian symmetry of spectra; time-reversal asymmetry causes chirality. This interesting phenomenon was first found in analyzing the chiral motion of the rattleback, a boat-shaped top having misaligned axes of inertia and geometry [Z. Yoshida et al., Phys. Lett. A 381, 2772–2777 (2017)]. To elucidate how chiral spectra are generated, we study the three-dimensional Lie–Poisson systems and classify the prototypes of singularities that cause symmetry breaking. The central idea is the deformation of the underlying Lie algebra; invoking Bianchi’s list of all three-dimensional Lie algebras, we show that the so-called class-B algebras, which are produced by asymmetric deformations of the simple algebra s o ( 3 ) , yield chiral spectra when linearized around their singularities. The theory of deformation is generalized to higher dimensions, including the infinite-dimensional Poisson manifolds relevant to fluid mechanics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Generalized Mercier stability criterion for stellarators

The Mercier criterion is a well-known stability criterion for tokamaks. It is derived from a 2 × 2 matrix eigenvalue problem arising from the expansion of resonant solutions about a singular surface where m−nq=0, with m and n being the poloidal and toroidal mode numbers, respectively, and q being the safety factor. The stability criterion is that the eigenvalues must be real, otherwise, the solution oscillates, violating the Newcomb crossing criterion. Because of the non-axisymmetry of stellarators, different toroidal as well as poloidal harmonics couple to each other. It follows that each singular surface can have multiple resonant harmonics, with multiplicity M≥1. The corresponding matrix eigenvalue problem involves a 2M×2M matrix, resulting in M pairs of positive and negative eigenvalues. The generalized stability criterion is that all eigenvalues must be real. While the original Mercier criterion can be expressed in terms of quadratures of equilibrium quantities over the singular surface, which can be evaluated anywhere, the generalized Mercier criterion can only be evaluated on rational q surfaces with a given set of resonant harmonics.

Physics↗
DOE OSTI2015

Plethora of transitions during breakup of liquid filaments

Significance Fluid flows, governed by nonlinear equations, permit formation of singularities. Often, singularities are artifacts of neglecting physical effects. However, free-surface flows exhibit observable singularities including filament pinch-off. As filaments thin, slightly (highly) viscous filaments are expected from theory to transition from an inertial (viscous) regime where viscosity (density) is negligible to an inertial–viscous regime where viscous and inertial effects are important. Previous works show this transition either does not occur or occurs for filament radii well below theoretical predictions. We demonstrate that thinning filaments unexpectedly pass through a number of intermediate transient regimes, thereby delaying onset of the final regime. The findings raise the question if similar dynamical transitions arise in problems that are not necessarily hydrodynamic in nature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Optimal renormalization of multiscale systems

While model order reduction is a promising approach in dealing with multi-scale time-dependent systems that are too large or too expensive to simulate for long times, the resulting reduced order models can suffer from instabilities. We have recently developed a time-dependent renormalization approach to stabilize such reduced models. In the current work, we extend this framework by introducing a parameter that controls the time-decay of the memory of such models and optimally selecting this parameter based on limited fully resolved simulations. First, we demonstrate our framework on the inviscid Burgers equation whose solution develops a finite-time singularity. Our renormalized reduced order models are stable and accurate for long times while using for their calibration only data from a full order simulation before the occurrence of the singularity. Furthermore, we apply this framework to the 3D Euler equations of incompressible fluid flow, where the problem of finite-time singularity formation is still open and where brute force simulation is only feasible for short times. Our approach allows us to obtain for the first time a perturbatively renormalizable model which is stable for long times and includes all the complex effects present in the 3D Euler dynamics. We find that, in each application, the renormalization coefficients display algebraic decay with increasing resolution, and that the parameter which controls the time-decay of the memory is problem-dependent.

Price, Jacob↗
DOE OSTI2020

Review: knots and other new topological effects in liquid crystals and colloids

Abstract Humankind has been obsessed with knots in religion, culture and daily life for millennia, while physicists like Gauss, Kelvin and Maxwell already involved them in models centuries ago. Nowadays, colloidal particles can be fabricated to have shapes of knots and links with arbitrary complexity. In liquid crystals, closed loops of singular vortex lines can be knotted by using colloidal particles and laser tweezers, as well as by confining nematic fluids into micrometer-sized droplets with complex topology. Knotted and linked colloidal particles induce knots and links of singular defects, which can be interlinked (or not) with colloidal particle knots, revealing the diversity of interactions between topologies of knotted fields and topologically nontrivial surfaces of colloidal objects. Even more diverse knotted structures emerge in nonsingular molecular alignment and magnetization fields in liquid crystals and colloidal ferromagnets. The topological solitons include hopfions, skyrmions, heliknotons, torons and other spatially localized continuous structures, which are classified based on homotopy theory, characterized by integer-valued topological invariants and often contain knotted or linked preimages, nonsingular regions of space corresponding to single points of the order parameter space. A zoo of topological solitons in liquid crystals, colloids and ferromagnets promises new breeds of information displays and a plethora of data storage, electro-optic and photonic applications. Their particle-like collective dynamics echoes coherent motions in active matter, ranging from crowds of people to schools of fish. This review discusses the state of the art in the field, as well as highlights recent developments and open questions in physics of knotted soft matter. We systematically overview knotted field configurations, the allowed transformations between them, their physical stability and how one can use one form of knotted fields to model, create and imprint other forms. The large variety of symmetries accessible to liquid crystals and colloids offer insights into stability, transformation and emergent dynamics of fully nonsingular and singular knotted fields of fundamental and applied importance. The common thread of this review is the ability to experimentally visualize these knots in real space. The review concludes with a discussion of how the studies of knots in liquid crystals and colloids can offer insights into topologically related structures in other branches of physics, with answers to many open questions, as well as how these experimentally observable knots hold a strong potential for providing new inspirations to the mathematical knot theory.

Physics↗
DOE OSTI2020

Annihilation-to-nothing: a quantum gravitational boundary condition for the Schwarzschild black hole

The interior of a static Schwarzschild metric can be written in terms of two functions, similar to some models of anisotropic cosmology. With a suitable choice of canonical variables, we solve the Wheeler-DeWitt equation (WDW) inside the horizon of a Schwarzschild black hole. By imposing classicality near the horizon, and requiring boundedness of the wave function, we get a rather generic solution of the WDW equation, whose steepest-descent solution, i.e., the ridge of the wave function, coincides nicely with the classical trajectory. However, there is an ambiguity in defining the arrow of time which leads to two possible interpretations—(i) if there is only one arrow of time, one can infer that the steepest-descent of the wave function follows the classical trajectory throughout: coming from the event horizon and going all the way down to the singularity, while (ii) if there are two different arrows of time in two separate regimes, it can be inferred that the steepest-descent of the wave function comes inwards from the event horizon in one region while it moves outwards from the singularity in the other region, and there exists an annihilation process of these two parts of the wave function inside the horizon. Additionally, adopting the second interpretation could shed light on the information loss paradox: as time goes on, probabilities for histories that include black holes and singularities decay to zero and eventually only trivial geometries dominate.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Magneto-thermal reconnection processes, related mode momentum and formation of high energy particle populations

In the context of a two-fluid theory of magnetic reconnection, when the longitudinal electron thermal conductivity is relatively large, the perturbed electron temperature tends to become singular in the presence of a reconnected field component and an electron temperature gradient. A finite transverse thermal diffusivity removes this singularity while a finite ‘inductivity’ can remove the singularity of the relevant plasma displacement. Then (i) a new ‘magneto-thermal’ reconnection producing mode, is found with characteristic widths of the reconnection layer remaining significant even when the macroscopic distances involved are very large; (ii) the mode phase velocities can be both in the direction of the electron diamagnetic velocity as well in the opposite (ion) direction. A numerical solution of the complete set of equations has been carried out with a simplified analytical reformulation of the problem. A sequence of processes is analyzed to point out that high-energy particle populations can be produced as a result of reconnection events. These processes involve mode-particle resonances transferring energy of the reconnecting mode to a superthermal ion population and the excitation of lower hybrid waves that can lead to a significant superthermal electron population. The same modes excited in axisymmetric (e.g. toroidal) confinement configurations can extract angular momentum from the main body of the plasma column and thereby sustain a local ‘spontaneous rotation’ of it.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Discovery of two bright high-redshift gravitationally lensed quasars revealed by Gaia

We present the discovery and preliminary characterisation of two high-redshift gravitationally lensed quasar systems in Gaia Data Release 2 (DR2). Candidates with multiple close-separation Gaia detections and quasar-like colours in WISE, Pan-STARRS, and DES are selected for follow-up spectroscopy with the New Technology Telescope. We confirm DES J215028.71-465251.3 as a $z$ = 4.130 ± 0.006 asymmetric, doubly imaged lensed quasar system and model the lensing mass distribution as a singular isothermal sphere. The system has an Einstein radius of 1.202 ± 0.005 arcsec and a predicted time delay of ~122.0 d between the quasar images, assuming a lensing galaxy redshift of $z$ = 0.5, making this a priority system for future optical monitoring. We confirm PS J042913.17+142840.9 as a $z$ = 3.866 ± 0.003 four-image quasar system in a cusp configuration, lensed by two foreground galaxies. The system is well modelled using a singular isothermal ellipsoid for the primary lens and a singular isothermal sphere for the secondary lens with Einstein radii 0.704 ± 0.006 and 0.241 ± 0.030 arcsec, respectively. A maximum predicted time delay of 9.6 d is calculated, assuming lensing galaxy redshifts of $z$ = 1.0. Furthermore, PS J042913.17+142840.9 exhibits a large flux ratio anomaly, up to a factor of 2.66 ± 0.37 in i band, that varies across optical and near-infrared wavelengths. We discuss LSST and its implications for future high-redshift lens searches and outline an extension to the search using supervised machine learning techniques.

Astronomy & Astrophysics↗
DOE OSTI2025

Interacting Dirac magnons in the van der Waals ferromagnet CrBr 3

We study the effects of magnon-magnon interactions in the two-dimensional van der Waals ferromagnet CrBr 3 focusing on its honeycomb lattice structure. Motivated by earlier theoretical predictions of temperature-induced spectral shifts and van Hove singularities in the magnon dispersion [S. S. Pershoguba et al., Phys. Rev. X 8, 011010 (2018)], we go beyond the commonly used thermal magnon approximation by applying second-order perturbation theory in a fully numerical framework. Our analysis uncovers significant deviations from previous analysis: in particular, the predicted singularities are absent, consistent with recent inelastic neutron scattering measurements [S. E. Nikitin et al ., Phys. Rev. Lett. 129, 127201 (2022)]. Moreover, we find that the temperature dependence of the renormalized magnon spectrum exhibits a distinct 𝑇 3 behavior for the optical magnon branch, while retaining 𝑇 2 behavior for the acoustic or down magnon band. This feature sheds light on the collective dynamics of Dirac magnons and their interactions. We further compare the honeycomb case with a triangular Bravais lattice, relevant for ferromagnetic monolayer MnBi 2 ⁢Te 4 , and show that both systems lack singular features while displaying quite distinct thermal trends.

Holstein-Primakoff method↗
DOE OSTI2024

Pairing and Pair Breaking by Gauge Fluctuations in Bilayer Composite Fermion Metals

We study interlayer pairing of composite fermions in the total $\nu=1/2+1/2$ quantum Hall bilayer as a possible framework for understanding the experimentally observed transition from a compressible state at large layer spacing to a bilayer quantum Hall state at small layer spacing. We consider a model in which the effective interlayer composite fermion pairing interaction mediated by the Chern-Simons gauge fields in the two layers is singular with both attractive (out-of-phase) and repulsive (in-phase) components diverging at low frequency. If only the more singular attractive interaction is included the pairing gap obtained by solving the gap equation is proportional to the inverse of the layer spacing squared. In the so-called local approximation, we find that when the less singular repulsive interactions are also included the pairing gap still falls off as inverse layer spacing squared, consistent with recent analyses, but is strongly suppressed to a degree that may account for the fact that this predicted inverse square dependence is not observed experimentally. The analytically obtained local approximation solutions are then used as a starting point to numerically iterate the full gap equation to assess the validity of the approximation in this limit.

Deng, Haoyun↗
DOE OSTI2024

Random projection using random quantum circuits

The random sampling task performed by Google's Sycamore processor gave us a glimpse of the “quantum supremacy era.” This has definitely shed some light on the power of random quantum circuits in this abstract task of sampling outputs from the (pseudo)random circuits. In this paper, we explore a practical near-term use of local random quantum circuits in dimensional reduction of large low-rank data sets. We make use of the well-studied dimensionality reduction technique called the random projection method. This method has been extensively used in various applications such as image processing, logistic regression, entropy computation of low-rank matrices, etc. We prove that the matrix representations of local random quantum circuits with sufficiently shorter depths [ ∼ O ( n ) ] serve as good candidates for random projection. We demonstrate numerically that their projection abilities are not far off from the computationally expensive classical principal components analysis on MNIST and CIFAR-100 image datasets. We also benchmark the performance of quantum random projection against the commonly used classical random projection in the tasks of dimensionality reduction of image data sets and computing von Neumann entropies of large low-rank density matrices. And finally, using variational quantum singular value decomposition, we demonstrate a near-term implementation of extracting the singular vectors with dominant singular values after quantum random projecting a large low-rank matrix to lower dimensions. All such numerical experiments unequivocally demonstrate the ability of local random circuits to randomize a large Hilbert space at sufficiently shorter depths with robust retention of properties of large data sets in reduced dimensions. Published by the American Physical Society 2024

Kumaran, Keerthi (ORCID:0009000949125721)↗
DOE OSTI2016

WE-AB-207A-04: Random Undersampled Cone Beam CT: Theoretical Analysis and a Novel Reconstruction Method

Purpose: Reducing x-ray exposure and speeding up data acquisition motived studies on projection data undersampling. It is an important question that for a given undersampling ratio, what the optimal undersampling approach is. In this study, we propose a new undersampling scheme: random-ray undersampling. We will mathematically analyze its projection matrix properties and demonstrate its advantages. We will also propose a new reconstruction method that simultaneously performs CT image reconstruction and projection domain data restoration. Methods: By representing projection operator under the basis of singular vectors of full projection operator, matrix representations for an undersampling case can be generated and numerical singular value decomposition can be performed. We compared properties of matrices among three undersampling approaches: regular-view undersampling, regular-ray undersampling, and the proposed random-ray undersampling. To accomplish CT reconstruction for random undersampling, we developed a novel method that iteratively performs CT reconstruction and missing projection data restoration via regularization approaches. Results: For a given undersampling ratio, random-ray undersampling preserved mathematical properties of full projection operator better than the other two approaches. This translates to advantages of reconstructing CT images at lower errors. Different types of image artifacts were observed depending on undersampling strategies, which were ascribed to the unique singular vectors of the sampling operators in the image domain. We tested the proposed reconstruction algorithm on a Forbid phantom with only 30% of the projection data randomly acquired. Reconstructed image error was reduced from 9.4% in a TV method to 7.6% in the proposed method. Conclusion: The proposed random-ray undersampling is mathematically advantageous over other typical undersampling approaches. It may permit better image reconstruction at the same undersampling ratio. The novel algorithm suitable for this random-ray undersampling was able to reconstruct high-quality images.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Solving the Electrodynamic Problem for a Microstrip Radiating Structure with a Chiral Substrate

We present a solution of the electrodynamic problem for a microstrip radiating structure with a substrate of a chiral metamaterial using the singular integral representation of the field, which in turn is reduced to a singular integral equation with the Cauchy-type singularity relative to the longitudinal component of the surface current density. Graphs of the current distribution for different types of substrates and the chirality parameters of a substrate are given.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Anomalous dips in reflection spectra of optical polymers deposited on plasmonic metals

Abstract We have studied reflection spectra of dye-doped and undoped polymers deposited onto Ag and Au substrates and found anomalous dips in the UV spectral range. On top of Ag substrates, the λ ∼ 375 nm dips were observed in undoped PMMA, PVP, and PS polymers as well as PMMA doped with Rh590 and HITC laser dyes. In silver-based samples, the spectral positions of the observed reflection dips were close to singularities in the refractive indexes of surface plasmon polaritons (SPPs) propagating at the interface between Ag and polymer. The latter singularities can tentatively explain the λ ∼ 375 nm reflection dips, if the scattering of Ag and polymeric films is large enough to launch SPP without any prism or grating. The dips observed in reflection of Rh590:PMMA and HITC:PMMA on top of Au, were more pronounced than those on Ag, broader, shifted to shorter wavelengths, and their spectral positions had large standard deviations. Furthermore, no anomalous dips in gold-based samples were observed in the reflection spectra of undoped PMMA, PVP, and PS polymers, and a modest singularity in the SPP refractive index, predicted theoretically at λ ∼ 500 nm, cannot explain the dips in the UV reflection spectra observed experimentally. It appears likely that the origin of the reflection dips on top of Au substrates is different from that on top of Ag substrates.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Diving into a holographic superconductor

Charged black holes in anti-de Sitter space become unstable to forming charged scalar hair at low temperatures T < T_\text{c} T < T c . This phenomenon is a holographic realization of superconductivity. We look inside the horizon of these holographic superconductors and find intricate dynamical behavior. The spacetime ends at a spacelike Kasner singularity, and there is no Cauchy horizon. Before reaching the singularity, there are several intermediate regimes which we study both analytically and numerically. These include strong Josephson oscillations in the condensate and possible `Kasner inversions’ in which after many e-folds of expansion, the Einstein-Rosen bridge contracts towards the singularity. Due to the Josephson oscillations, the number of Kasner inversions depends very sensitively on T T , and diverges at a discrete set of temperatures \{T_n\} { T n } that accumulate at T_c T c . Near these T_n T n , the final Kasner exponent exhibits fractal-like behavior.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Natural boundaries for scattering amplitudes

Singularities, such as poles and branch points, play a crucial role in investigating the analytic properties of scattering amplitudes that inform new computational techniques. In this note, we point out that scattering amplitudes can also have another class of singularities called natural boundaries of analyticity. They create a barrier beyond which analytic continuation cannot be performed. More concretely, we use unitarity to show that 2 \to 2 2 → 2 scattering amplitudes in theories with a mass gap can have a natural boundary on the second sheet of the lightest threshold cut. There, an infinite number of ladder-type Landau singularities densely accumulates on the real axis in the center-of-mass energy plane. We argue that natural boundaries are generic features of higher-multiplicity scattering amplitudes in gapped theories.

97 MATHEMATICS AND COMPUTING↗