Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Singularities in general relativity”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Distributed Embedded Energy Converter Technologies for Marine Renewable Energy (A Technical Report)

The domain of distributed embedded energy converter technologies (DEEC-Tec) is a nascent and underexplored paradigm for harvesting and converting marine renewable energy. The paradigm distinguishes itself through its use of many small distributed embedded energy converters (DEECs) that, ultimately, are assembled through the creation of "DEEC-Tec metamaterials" to create an overall larger marine renewable energy harvesting and converting structure. As an example, such a structure could be an ocean wave energy converter - a converter whose structure is made from various types of DEEC-Tec metamaterials that harvests ocean wave energy and converts that energy into something more useful such as electricity. To that end, DEEC-Tec can be viewed at three different technology levels: (1) individual distributed embedded energy converters, also known as DEECs; (2) DEEC-Tec metamaterials-essentially, pseudo-materials made from the interconnection of many DEECs; and (3) overall larger complete marine renewable energy harvesting-converting structures-these structures being made from DEEC-Tec metamaterials. Arising directly from the application of DEEC-Tec to harvest and convert ocean wave energy are several noteworthy benefits, some of which include: (1) the lack of load concentrations into singular components or subsystems, (2) broad-banded ocean wave energy frequency harvesting and conversion, and (3) inherent redundancy-failure of some individual DEECs does not represent a failure of an entire DEEC-Tec-based WEC. This report describes DEEC-Tec by way of descriptions of those three technology levels: individual DEECs, DEEC-Tec metamaterials, and DEEC-Tec-based WECs. Moreover, the report describes corresponding research approaches and methodologies for related concepts such as DEEC-Tec-based WEC topologies and morphologies in addition to manufacturing and fabrication techniques found suitable for the application of DEEC-Tec within the general domain of marine renewable energy-moving beyond only ocean wave energy conversion.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

A quantum hamiltonian simulation benchmark

Hamiltonian simulation is one of the most important problems in quantum computation, and quantum singular value transformation (QSVT) is an efficient way to simulate a general class of Hamiltonians. However, the QSVT circuit typically involves multiple ancilla qubits and multi-qubit control gates. In order to simulate a certain class of n-qubit random Hamiltonians, we propose a drastically simplified quantum circuit that we refer to as the minimal QSVT circuit, which uses only one ancilla qubit and no multi-qubit controlled gates. We formulate a simple metric called the quantum unitary evolution score (QUES), which is a scalable quantum benchmark and can be verified without any need for classical computation. Under the globally depolarized noise model, we demonstrate that QUES is directly related to the circuit fidelity, and the potential classical hardness of an associated quantum circuit sampling problem. Under the same assumption, theoretical analysis suggests there exists an ‘optimal’ simulation time t opt ≈ 4.81, at which even a noisy quantum device may be sufficient to demonstrate the potential classical hardness.

97 MATHEMATICS AND COMPUTING↗
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↗
NASA NTRS1975

Steady subsonic flow around finite-thickness wings

The general method for analyzing steady subsonic potential aerodynamic flow around a lifting body having arbitrary shape is presented. By using the Green function method, an integral representation for the potential is obtained. Under small perturbation assumption, the potential at any point, P, in the field depends only upon the values of the potential and its normal derivative on the surface of the body. Hence if the point P approaches the surface of the body, the representation reduces to an integral equation relating the potential and its normal derivative (which is known from the boundary conditions) on the surface. The question of uniqueness is examined and it is shown that, for thin wings, the operator becomes singular as the thickness approaches zero. This fact may yield numerical problems for very thin wings. However, numerical results obtained for a rectangular wing in subsonic flow show that these problems do not appear even for thickness ratio tau = .001. Comparison with existing results shows that the proposed method is at least as fast and accurate as the lifting surface theories.

Kuo, C. C.↗
DOE OSTI2022

Linking the singularities of cosmological correlators

Much of the structure of cosmological correlators is controlled by their singularities, which in turn are fixed in terms of flat-space scattering amplitudes. An important challenge is to interpolate between the singular limits to determine the full correlators at arbitrary kinematics. This is particularly relevant because the singularities of correlators are not directly observable, but can only be accessed by analytic continuation. In this paper, we study rational correlators — including those of gauge fields, gravitons, and the inflaton — whose only singularities at tree level are poles and whose behavior away from these poles is strongly constrained by unitarity and locality. We describe how unitarity translates into a set of cutting rules that consistent correlators must satisfy, and explain how this can be used to bootstrap correlators given information about their singularities. We also derive recursion relations that allow the iterative construction of more complicated correlators from simpler building blocks. In flat space, all energy singularities are simple poles, so that the combination of unitarity constraints and recursion relations provides an efficient way to bootstrap the full correlators. In many cases, these flat-space correlators can then be transformed into their more complex de Sitter counterparts. As an example of this procedure, we derive the correlator associated to graviton Compton scattering in de Sitter space, though the methods are much more widely applicable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1979

An examination of a group-velocity criterion for the breakdown of an idealized vortex flow

The phenomenon of vortex breakdown is believed to be associated with a finite amplitude wave that has become trapped at the critical or breakdown location. The conditions at which the propagating waves become trapped at a certain axial location were examined by use of a group-velocity criterion implied by Landahl's general theory of wave trapping. An ideal vortex having constant vorticity and uniform axial velocity at the inlet of a slowly diverging duct was studied. The linear wave propagation analysis is applied to the base flow at several axial stations for several values of the ratio of swirl velocity to axial velocity at the inlet of the divergent duct, assuming a locally parallel flow. The dipsersion relations and hence the group velocities of both the symmetric (n = 0) and asymmetric modes (n = + or - 1) were investigated. The existence of a critical state in the flow (at which the group velocity vanishes), and its relationship to the stagnation point on the axis of the duct and to the occurrence of an irregular singularity in the equations governing wave propagation in the flow field are discussed.

Tsai, C. Y.↗
NASA NTRS2002

A Proposed Model for Protein Crystal Nucleation and Growth

How does one take a molecule, strongly asymmetric in both shape and charge distribution, and assemble it into a crystal? We propose a model for the nucleation and crystal growth process for tetragonal lysozyme, based upon fluorescence, light, neutron, and X-ray scattering data, size exclusion chromatography experiments, dialysis kinetics, AFM, and modeling of growth rate data, from this and other laboratories. The first species formed is postulated to be a 'head to side' dimer. Through repeating associations involving the same intermolecular interactions this grows to a 4(sub 3) helix structure, that in turn serves as the basic unit for nucleation and subsequent crystal growth. High salt attenuates surface charges while promoting hydrophobic interactions. Symmetry facilitates subsequent helix-helix self-association. Assembly stability is enhanced when a four helix structure is obtained, with each bound to two neighbors. Only two unique interactions are required. The first are those for helix formation, where the dominant interaction is the intermolecular bridging anion. The second is the anti-parallel side-by-side helix-helix interaction, guided by alternating pairs of symmetry related salt bridges along each side. At this stage all eight unique positions of the P4(sub3)2(sub 1),2(sub 1) unit cell are filled. The process is one of a) attenuating the most strongly interacting groups, such that b) the molecules begin to self-associate in defined patterns, so that c) symmetry is obtained, which d) propagates as a growing crystal. Simple and conceptually obvious in hindsight, this tells much about what we are empirically doing when we crystallize macromolecules. By adjusting the growth parameters we are empirically balancing the intermolecular interactions, preferentially attenuating the dominant strong (for lysozyme the charged groups) while strengthening the lesser strong (hydrophobic) interactions. In the general case for proteins the lack of a singularly defined association pathway may lead to formation of multiple species, i.e., amorphous precipitation. Weak interactions, such as hydrogen bonds, are promiscuous, serving to strengthen rather than define specific interactions. Participation in an interaction sequesters that surface from subsequent interactions, and we expect the strongest bonds to form first. This model, its basis, how it fits into the currently understood osmotic second virial coefficient approach to crystallization, and what it suggests will be discussed.

Pusey, Marc↗
DOE OSTI2023

Hidden singularities in spontaneously polarized molecular solids

Films of dipolar molecules formed by physical vapor deposition are, in general, spontaneously polarized, with implications ranging from electron transfer in molecular optoelectronic devices to the properties of astrochemical ices in the interstellar medium. Polarization arises from dipole orientation, which should intuitively decrease with increasing deposition temperature, T. However, it is experimentally found that minimum or maximum values in polarization vs T may be observed for cis-methyl formate, 1-propanol, and ammonia. A continuous analytic form of polarization vs T is developed, which has the property that it is not differentiable at all T. The minima and maxima in polarization vs T are marked by singularities in the differential of this analytic form. In conclusion, this exotic behavior is presently unique to films of dipolar species and has not been reported, for example, in the related magnetic phases of spin glasses.

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

The Fate of Nucleated Black Holes in de Sitter Quantum Gravity

The Euclidean Nariai geometry has long been proposed as the instanton describing the nucleation of maximal-mass black holes in de Sitter space. We place this interpretation on firmer footing by showing that, once an observer is included, the gravitational path integral produces the imaginary phase required for a transition rate. As a warmup, we revisit the Hawking-Moss instanton and, as a byproduct, find that scalar fields can enhance black-hole nucleation, suggesting a quantum-gravity bound on scalar potentials with de Sitter solutions. We then study the subsequent semiclassical evolution of the nucleated black hole. We show that the previously claimed "anti-evaporation" channel is unphysical, arising from a quantum state with singular horizons. In a smooth state, the black hole instead undergoes standard thermal Hawking evaporation. We verify explicit agreement with the no-boundary state and argue that this evaporation is not subject to large quantum-gravity corrections. The nucleated black hole thus evaporates completely back to the maximally-entropic empty de Sitter vacuum, making the full process a Boltzmann fluctuation.

FOS: Physical sciences↗
DOE OSTI2025

A case study in contrastive learning information combination: Application to technical forensics of additive manufacturing filament source identification

Combination of information from disparate data sources into a single decision is a core challenge in many fields, including the field of technical forensics. Technical forensics (TF) utilizes technical characterization of questioned samples to determine properties of that sample; these properties are then used to infer information of forensic interest, such as provenance, age, or attribution. TF is utilized in traditional forensic applications, such as the attribution of material fragments from an explosive, and in nuclear forensic applications, such as the attribution of actinides which have been interdicted out of regulatory control. The challenge of combining information from disparate sources, described alternately by many terms including “Data Fusion” and “Data Integration”, is exacerbated in the technical forensics domain due to at least two factors: the challenge of interpreting each information source singularly, and the relatively small data set sizes available. Extensive literature exists attempting to combine technical forensics information sources, both in manual and automated processes. These attempts are often bespoke to the specific information sources (such as the bi-, tri-, or quad-isotope chart (Moody, Grant, and Hutcheon 2005)), with some emerging examples of simple early- and late- fusion (, respectively). Simultaneous to the information combination efforts described in the previous paragraph, the field of natural language processing attempted (and largely succeeded) in combining information from multiple non-technical information sources. The ecosystem of “multi-modal” language models, which can take text and images as input, and generate text and images as output, became large and diverse by 2025 (Khan et al. 2025). In a generalized sense, many of these methods are trained by learning neural networks which can convert raw text or images into a vector of numbers describing the text or image, hereafter called “embeddings” and the neural networks performing the conversion are called “embedders”. By using a separate embedder for text and images, finding coincident text and images (such as images with their captions), and optimizing the parameters of the embedders such that the embeddings for the text and the image are similar, the field has found a bridge between text and images (Girdhar et al. 2023). It is the contention of the authors of this report that this insight is not limited to text and images but instead can be extended to any modality which can be found coincidently. The subject of the rest of this report is the application of this method to example multi-modal technical forensic data. Some details about the data used in this report are not appropriate for this report, and are included in a companion report (PNNL-38669).

36 MATERIALS SCIENCE↗
NASA NTRS1982

Creep crack-growth: A new path-independent integral (T sub c), and computational studies

The development of valid creep fracture criteria is considered. Two path-independent integral parameters which show some degree of promise are the C* and (Delta T)sub c integrals. The mathematical aspects of these parameters are reviewed by deriving generalized vector forms of the parameters using conservation laws which are valid for arbitrary, three dimensional, cracked bodies with crack surface tractions (or applied displacements), body forces, inertial effects, and large deformations. Two principal conclusions are that (Delta T)sub c has an energy rate interpretation whereas C* does not. The development and application of fracture criteria often involves the solution of boundary/initial value problems associated with deformation and stresses. The finite element method is used for this purpose. An efficient, small displacement, infinitesimal strain, displacement based finite element model is specialized to two dimensional plane stress and plane strain and to power law creep constitutive relations. A mesh shifting/remeshing procedure is used for simulating crack growth. The model is implemented with the quartz-point node technique and also with specially developed, conforming, crack-tip singularity elements which provide for the r to the n-(1+n) power strain singularity associated with the HRR crack-tip field. Comparisons are made with a variety of analytical solutions and alternate numerical solutions for a number of problems.

Stonesifer, R. B.↗
NASA NTRS1989

Development of a finite element based delamination analysis for laminates subject to extension, bending, and torsion

Delamination is a common failure mode of laminated composite materials. This type of failure frequently occurs at the free edges of laminates where singular interlaminar stresses are developed due to the difference in Poisson's ratios between adjacent plies. Typically the delaminations develop between 90 degree plies and adjacent angle plies. Edge delamination has been studied by several investigators using a variety of techniques. Recently, Chan and Ochoa applied the quasi-three-dimensional finite element model to the analysis of a laminate subject to bending, extension, and torsion. This problem is of particular significance relative to the structural integrity of composite helicopter rotors. The task undertaken was to incorporate Chan and Ochoa's formulation into a Raju Q3DG program. The resulting program is capable of modeling extension, bending, and torsional mechanical loadings as well as thermal and hygroscopic loadings. The addition of the torsional and bending loading capability will provide the capability to perform a delamination analysis of a general unsymmetric laminate containing four cracks, each of a different length. The solutions obtained using this program are evaluated by comparing them with solutions from a full three-dimensional finite element solution. This comparison facilitates the assessment of three dimensional affects such as the warping constraint imposed by the load frame grips. It wlso facilitates the evaluation of the external load representation employed in the Q3D formulation. Finally, strain energy release rates computed from the three-dimensional results are compared with those predicted using the quasi-three-dimensional formulation.

Hooper, Steven J.↗
NASA NTRS1992

Robust adaptive kinematic control of redundant robots

The paper presents a general method for the resolution of redundancy that combines the Jacobian pseudoinverse and augmentation approaches. A direct adaptive control scheme is developed to generate joint angle trajectories for achieving desired end-effector motion as well as additional user defined tasks. The scheme ensures arbitrarily small errors between the desired and the actual motion of the manipulator. Explicit bounds on the errors are established that are directly related to the mismatch between actual and estimated pseudoinverse Jacobian matrix, motion velocity and the controller gain. It is shown that the scheme is tolerant of the mismatch and consequently only infrequent pseudoinverse computations are needed during a typical robot motion. As a result, the scheme is computationally fast, and can be implemented for real-time control of redundant robots. A method is incorporated to cope with the robot singularities allowing the manipulator to get very close or even pass through a singularity while maintaining a good tracking performance and acceptable joint velocities. Computer simulations and experimental results are provided in support of the theoretical developments.

Tarokh, M.↗
NASA NTRS2015

Robust Stability Analysis of the Space Launch System Control Design: A Singular Value Approach

Classical stability analysis consists of breaking the feedback loops one at a time and determining separately how much gain or phase variations would destabilize the stable nominal feedback system. For typical launch vehicle control design, classical control techniques are generally employed. In addition to stability margins, frequency domain Monte Carlo methods are used to evaluate the robustness of the design. However, such techniques were developed for Single-Input-Single-Output (SISO) systems and do not take into consideration the off-diagonal terms in the transfer function matrix of Multi-Input-Multi-Output (MIMO) systems. Robust stability analysis techniques such as H(sub infinity) and mu are applicable to MIMO systems but have not been adopted as standard practices within the launch vehicle controls community. This paper took advantage of a simple singular-value-based MIMO stability margin evaluation method based on work done by Mukhopadhyay and Newsom and applied it to the SLS high-fidelity dynamics model. The method computes a simultaneous multi-loop gain and phase margin that could be related back to classical margins. The results presented in this paper suggest that for the SLS system, traditional SISO stability margins are similar to the MIMO margins. This additional level of verification provides confidence in the robustness of the control design.

Pei, Jing↗
NASA NTRS1980

On Rosen's theory of gravity and cosmology

Formal similarities between general relativity and Rosen's bimetric theory of gravity were used to analyze various bimetric cosmologies. The following results were found: (1) physically plausible model universes which have a flat static background metric, have a Robertson-Walker fundamental metric, and which allow co-moving coordinates do not exist in bimetric cosmology. (2) it is difficult to use the Robertson-Walker metric for both the background metric (gamma mu nu) and the fundamental metric tensor of Riemannian geometry( g mu nu) and require that g mu nu and gamma mu nu have different time dependences. (3) A consistency relation for using co-moving coordinates in bimetric cosmology was derived. (4) Certain spatially flat bimetric cosmologies of Babala were tested for the presence of particle horizons. (5) An analytic solution for Rosen's k = +1 model was found. (6) Rosen's singularity free k = +1 model arises from what appears to be an arbitary choice for the time dependent part of gamma mu nu.

Barnes, R. C.↗
NASA NTRS1991

Robust control of systems with real parameter uncertainty and unmodelled dynamics

During this research period we have made significant progress in the four proposed areas: (1) design of robust controllers via H infinity optimization; (2) design of robust controllers via mixed H2/H infinity optimization; (3) M-delta structure and robust stability analysis for structured uncertainties; and (4) a study on controllability and observability of perturbed plant. It is well known now that the two-Riccati-equation solution to the H infinity control problem can be used to characterize all possible stabilizing optimal or suboptimal H infinity controllers if the optimal H infinity norm or gamma, an upper bound of a suboptimal H infinity norm, is given. In this research, we discovered some useful properties of these H infinity Riccati solutions. Among them, the most prominent one is that the spectral radius of the product of these two Riccati solutions is a continuous, nonincreasing, convex function of gamma in the domain of interest. Based on these properties, quadratically convergent algorithms are developed to compute the optimal H infinity norm. We also set up a detailed procedure for applying the H infinity theory to robust control systems design. The desire to design controllers with H infinity robustness but H(exp 2) performance has recently resulted in mixed H(exp 2) and H infinity control problem formulation. The mixed H(exp 2)/H infinity problem have drawn the attention of many investigators. However, solution is only available for special cases of this problem. We formulated a relatively realistic control problem with H(exp 2) performance index and H infinity robustness constraint into a more general mixed H(exp 2)/H infinity problem. No optimal solution yet is available for this more general mixed H(exp 2)/H infinity problem. Although the optimal solution for this mixed H(exp 2)/H infinity control has not yet been found, we proposed a design approach which can be used through proper choice of the available design parameters to influence both robustness and performance. For a large class of linear time-invariant systems with real parametric perturbations, the coefficient vector of the characteristic polynomial is a multilinear function of the real parameter vector. Based on this multilinear mapping relationship together with the recent developments for polytopic polynomials and parameter domain partition technique, we proposed an iterative algorithm for coupling the real structured singular value.

Chang, Bor-Chin↗
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↗
NASA NTRS1993

Tension fracture of laminates for transport fuselage. Part 2: Large notches

Tests were conducted on over 200 center-crack specimens to evaluate: (a) the tension-fracture performance of candidate materials and laminates for commercial fuselage applications; and (b) the accuracy of several failure criteria in predicting response. Crack lengths of up to 12 inches were considered. Other variables included fiber/matrix combination, layup, lamination manufacturing process, and intraply hybridization. Laminates fabricated using the automated tow-placement process provided significantly higher tension-fracture strengths than nominally identical tape laminates. This confirmed earlier findings for other layups, and possibly relates to a reduced stress concentration resulting from a larger scale of repeatable material inhomogeneity in the tow-placed laminates. Changes in material and layup result in a trade-off between small-notch and large-notch strengths. Toughened resins and 0 deg-dominate layups result in higher small-notch strengths but lower large-notch strengths than brittle resins, 90 deg and 45 deg dominated layups, and intraply S2-glass hybrid material forms. Test results indicate that strength-prediction methods that allow for a reduced order singularity of the crack-tip stress field are more successful at predicting failure over a range of notch sizes than those relying on the classical square-root singularity. The order of singularity required to accurately predict large-notch strength from small-notch data was affected by both material and layup. Measured crack-tip strain distributions were generally higher than those predicted using classical methods. Traditional methods of correcting for finite specimen width were found to be lacking, confirming earlier findings with other specimen geometries. Fracture tests of two stiffened panels, identical except for differing materials, with severed central stiffeners resulted in nearly identical damage progression and failure sequences. Strain-softening laws implemented within finite element models appear attractive to account for load redistribution in configured structure due to damage-induced crack tip softening

Walker, Tom H.↗