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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Realistic gravitational focusing of meteoroid streams

A meteor shower can be modeled, to first order, as a parallel stream of particles encountering the Earth. Gravity bends the trajectories of these particles inward, producing an increase in flux both near the planet and along the anti-radiant line. This effect is known as gravitational focusing, and the interception of a portion of meteoroid trajectories by the Earth or another massive body is known as planetary shielding. For a perfectly parallel meteoroid stream, gravitational focusing produces a flux singularity along the anti-radiant line, near which intense enhancements in flux occur. In reality, meteoroids will exhibit some dispersion in their motion that removes this singularity. We present a modified analytic treatment of gravitational focusing and planetary shielding that takes this dispersion into account and discuss its ramifications for several problems in meteor astronomy.

Althea Valkyrie Moorhead↗
NASA NTRS2023

A 1D Model for Nucleation of Ice From Aerosol Particles: An Application to a Mixed-Phase Arctic Stratus Cloud Layer

Mixed-phase clouds (MPCs) have been identified as significant contributors to uncertainties in climate projections, attributable to model representation of processes controlling the formation and loss of supercooled water droplets and ice particles from the atmosphere. Arctic MPCs are commonly widespread and long-lived, with sustained ice crystal formation processes that challenge current understanding. This study examines the ice-nucleating particle (INP) reservoir dynamics governing immersion-mode heterogeneous freezing in an observed case of Arctic MPCs using a simplified 1D aerosol-cloud model. The model setup includes prescribed dynamical forcings and thermodynamic profiles, and represents INPs as multicomponent and polydisperse particle size distributions. Diagnostic and prognostic approaches to immersion freezing parameterization are compared, including time-independent (singular) number- and surface area-based descriptions and a time-dependent description following classical nucleation theory (CNT). The choice of freezing parameterization defines the size of the INP reservoir. The CNT-based description yields an orders of magnitude larger INP reservoir than the singular parameterizations, which is the dominant factor for sustained ice crystal formation. The efficiency of the freezing process and cloud cooling are of secondary importance. A diagnostic treatment neglecting INP loss is only accurate when the INP reservoir size is large and INP depletion weak. Since a larger INP reservoir sustains ice crystal formation substantially longer, and ice water path scales with ice crystal concentrations for the conditions considered, resolving the source of differences in INP reservoir dynamics due to model implementation is a high priority for advancing climate model physics.

aerosol↗
NASA NTRS2024

Spacecraft Modeling, Attitude Determination, and Control: Quaternion-Based Approach

Spacecraft attitude determination and control is an important part of a spacecraft to achieve its designed mission. As of today, many spacecrafts have been successfully launched, and most of them have performed well as they were designed. Many research papers have been published to address the attitude determination and control design problems. Several text books are available for students to learn the technology and for engineers to use as references. The most popular spacecraft models for attitude determination algorithms and control design methods are the Euler angle models and the quaternion models . The Euler angle models have been proved very efficient as the linearized models are controllable, and all standard linear control system design methods are directly applicable. The drawbacks related to the Euler angle methods are (a) the designs based on linearized models may not globally stabilize the original nonlinear spacecraft, i.e., the design may not work when the attitude of the spacecraft is far away from the point where the linearization is performed; (b) the models depend on the rotational sequences, this can be error prone if several teams work on the same project and they use different rotational sequences; (c) for any rotational sequence, there is a singular point where the model is not applicable; and (d) since most attitude determination methods use quaternion to represent the spacecraft attitude, there is a need to transform quaternion into Euler angles. On the other hand, for quaternion models, people have found controllers that can globally stabilize nonlinear spacecraft systems; the models do not depend on rotational sequences and they have no singular point; and the quaternion is provided by attitude determination system and ready to use. The main problem with the quaternion model based control system design is that the linearized quaternion model is not controllable. Therefore, most published design methods heavily rely on Lyapunov functions for the nonlinear spacecraft system. But there is no systematic way to obtain a desired Lyapunov functions. Moreover, the Lyapunov function based designs focus on the closed-loop system stability but pay little attention to the closed-loop system performance. In a series of papers, the author proposed some reduced quaternion models which lead to some controllable linearized spacecraft models. Therefore, all standard linear system theory can be directly applied to analyze and design the spacecraft control systems. We showed that, in some cases, the designed control system is not only optimal for the linearized system, but also globally stabilize the original nonlinear system . Clearly, the reduced quaternion models do not depend on rotational sequences. Due to the special structure of the linearized spacecraft model, some most important design methods, such as LQR design and robust pole assignment design are very simple, enjoy the analytical solutions for some problems, have direct connection to the performance measures, such as settling time, rising time , and percentage of overshoot . All these features are attractive for high quality control system designs. The idea mentioned above is then extended to more spacecraft control problems using specific actuators such as magnetic torque bars and control momentum gyroscopes. These types of actuators may not provide exactly desired torques. Most existing methods use different conversions to get approximate solutions, meaning that these actuators may generate a torque close to but not equal to the desired one. Using the reduced quaternion models that incorporate the actuators into the system model, the control inputs are not torques but the operational parameters. The main benefit of this idea is that the control actions are not approximate but accurate. As all actuators have their operational limit, design with input constraints are also considered in this book by using recently developed interior-point optimization techniques. This book grows up from my research on the spacecraft attitude determination and control design methods in more than a decade which is focused on using reduced quaternion models because of their merits stated above. The book provides all necessary background materials on orbital dynamics, rotations and quaternion, frequently used reference frames, transformations between reference frames, space environment and disturbance torques, ephemeris astronomical vector calculations and measurement instruments, spacecraft control actuators and their models, so that the readers will get a global picture and can apply all these information into the spacecraft system modeling, attitude determination, and spacecraft control system designs, which is the main purpose of this book. This book is different from existing books in that we focus on quaternion based spacecraft control system designs and we consider only attitude control system design related problems, from spacecraft modeling, to attitude determination and estimation, to control system design method selection, to control algorithm development, and to the simulation of the control system designs. Moreover, this book addresses different attitude control tasks in the spacecraft life cycle, including spacecraft maneuver, orbit raising, attitude control, and rendezvous. Finally, this book emphasizes the state space design methods rather than the classical frequency design methods.

Attitude determination and control↗
DOE OSTI2020

Computationally Efficient Decompositions of Oblique Projection Matrices

Oblique projection matrices arise in problems in weighted least squares, signal processing, and optimization. While these matrices can be potentially very large, their low-rank structure can be exploited for efficient computation. Here, we propose fast and scalable algorithms for computing their eigendecomposition and singular value decomposition (SVD). Numerical experiments that compare our proposed approaches to existing methods, including randomized SVD, are presented. In addition, we test their accuracy on linear systems from equality constrained optimization problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Nanoscale Mapping and Defect-Assisted Manipulation of Surface Plasmon Resonances in 2D Bi 2 Te 3 /Sb 2 Te 3 In-Plane Heterostructures

Here, the Bi 2 Te 3 /Sb 2 Te 3 in-plane heterostructure is reported as a low-dimensional tunable chalcogenide well suited as plasmonic building block for the visible-UV spectral range. Electron-driven plasmon excitations of low-dimensional Bi 2 Te 3 /Sb 2 Te 3 are investigated by monochromated electron energy loss spectroscopy spectrum imaging. To resolve the nanoscale spatial distribution of various local plasmonic resonances, singular value decomposition is used to disentangle the spectral data and identify the individual spectral contributions of various corner, edge, and face modes. Furthermore, defect-plasmon interactions are investigated both for nanoscale intrinsic and thermally induced extrinsic polygonal defects (in situ sublimation). Signature of defect-induced red shift ranging from a several hundreds of millielectronvolts to a few electronvolts, broadening of various plasmon response, together with selective enhancement and significant variations in their intensity are detected. This study highlights the presence of a heterointerface and identifies defects as physical tuning pathways to modulate the plasmonic response over a broad spectral range. Finally, the experimental observations are compared qualitatively and validated with numerical simulations using the electron-driven discrete dipole approximation. Low-dimensional Bi 2 Te 3 /Sb 2 Te 3 as a less explored plasmonic system holds great promises as emerging platform for integrated plasmonics. Furthermore, introducing controlled structural defects can open the door for nanoengineering of plasmonic properties in such systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Obstacles to constructing de Sitter space in string theory

There have been many attempts to construct de Sitter space-times in string theory. While arguably there have been some successes, this has proven challenging, leading to the de Sitter swampland conjecture: quantum theories of gravity do not admit stable or metastable de Sitter space. Here we explain that, within controlled approximations, one lacks the tools to construct de Sitter space in string theory. Such approximations would require the existence of a set of (arbitrarily) small parameters, subject to severe constraints. But beyond this one also needs an understanding of big-bang and big-crunch singularities that is not currently accessible to standard approximations in string theory. The existence or non-existence of metastable de Sitter space in string theory remains a matter of conjecture.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Rescuing a black hole in the large- q coupled SYK model

In this paper, we develop a general effective theory for two copies of the Sachdev-Ye-Kitaev (SYK) model with a time-dependent bilinear coupling. For a quantum quench problem with an initial state of the thermofield double state, we show how the evolution of the system is described by a complex reparametrization field with a classical Hamiltonian. We study correlation functions in this system and compare the large-q theory with the bulk low energy effective theory. In particular, we study the special case of a “rescued black hole”, which describes how a time-evolved thermofield double state can evolve to the ground state of a coupled SYK model by a carefully tuned time-dependent coupling. In the low energy region, there is a holographic dual interpretation, which is a geometry that crosses over from an eternal black hole to a global AdS2 vacuum. This family of geometries allow us to access the bulk region that would be the black hole interior without the rescue process. By comparing the large-q and low energy theory, we find that even in the low energy region the deviation from the low energy theory cannot be neglected if the rescue process starts late. This provides evidence that the low energy effective theory of the bulk fails near the inner horizon of the black hole. We note the possibility of a connection to a two-dimensional analog of the higher-dimensional black hole singularity.

2D Gravity↗
DOE OSTI2021

Holography abhors visible trapped surfaces

We prove that consistency of the holographic dictionary implies a hallmark prediction of the weak cosmic censorship conjecture: that in classical gravity, trapped surfaces lie behind event horizons. In particular, the existence of a trapped surface implies the existence of an event horizon, and that furthermore this event horizon must be outside of the trapped surface. More precisely, we show that the formation of event horizons outside of a strong gravity region is a direct consequence of causal wedge inclusion, which is required by entanglement wedge reconstruction. We make few assumptions beyond the absence of evaporating singularities in strictly classical gravity. We comment on the implication that spacetimes with naked trapped surfaces do not admit a holographic dual, note a possible application to holographic complexity, and speculate on the dual CFT interpretation of a trapped surface.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Shape-constrained Input Estimation for Efficient Multi-shaker Vibration Testing

Multi-shaker vibration testing is gaining interest from structural dynamics test engineers as it can provide a much more accurate match to complicated field vibration responses than traditional single-axis shaker tests. However, the force capabilities of the small modal shakers typically used in multi-shaker vibration tests has limited the achievable response levels. To date, most multi-shaker vibration tests have been performed using a variety of standard, commercially-available control systems. While these control systems are adequate for a wide range of multiple-input/multiple-output tests, their control algorithms have not been tailored for the specific problem of multi-shaker vibration tests: efficiently coordinating the various shakers to work together to achieve a desired response. Here, a new input estimation algorithm is developed and demonstrated using simulations and actual test data. This algorithm, dubbed shape-constrained input estimation, is shown to effectively coordinate multiple shakers using a set of constraint vectors based on the deflection shapes of the test structure. This is accomplished by using the singular vector shapes of the system frequency response matrix, which allows the constraint vectors to automatically change as a function of frequency. Simulation and test results indicate a significant reduction in the input forces required to achieve a desired response. Finally, the results indicate that shape-constrained input estimation is an effective method to achieve higher response levels from limited shaker forces which will enable higher level multi-shaker vibration tests to be performed.

42 ENGINEERING↗
DOE OSTI2021

Some aspects of the five-dimensional Lovelock black hole spacetime: Strong homotopy retract, perihelion precession and quasistationary levels

Highlights: • Some topological aspects are explored. • The perihelion precession is calculated. • The general Heun functions are applied to obtain analytical solutions. • The quasibound states and the radial wave eigenfunctions are obtained. In this work we explore some mathematical physics aspects of the spherically symmetric Lovelock black hole in high dimensions. Intended for this aim, we thoroughly consider the metric corresponding to the five-dimensional Lovelock black hole spacetime. We construct the strong retractions by the geodesic equations on the background under consideration. As a result, from the topological point of view, we construct the theory of strong homotopy retract, which will allow us, in principle, to better understand some of its suitable applications on astrophysics and cosmology, in particular, in the analysis of the spacetime singularities. We find the solutions of the equation of motion for both radial and angular coordinates, and then we describe the outer (“exterior”) and lower (“interior”) apparent horizons. Indeed, the outer apparent horizon is the last surface from which the light waves could still escape from the black hole. Thus, it is meaningful to analyze some physical phenomena related to quantum particles propagating outside the exterior apparent horizon, in particular, we discuss the quasistationary levels of scalar fields and their radial wave functions, which are given in terms of the general Heun functions. We also calculate the perihelion precession in this background.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Advantages of the Samarskii-type schemes on the Shishkin mesh

The schemes of the Samarskii type are simple modifications of the upwind scheme. We use them on the Shishkin mesh and discuss their advantages over the upwind scheme when applied to the linear one-dimensional singularly perturbed convection–diffusion problem. One of the advantages is that the Samarskii-type schemes have exact first-order accuracy uniform in the perturbation parameter, as opposed to the upwind scheme which is almost first-order uniformly accurate because its accuracy is diminished by logarithmic factors. Although this is not a new result, we re-emphasize it in the paper. We also demonstrate another advantage, that the Samarskii-type schemes are almost second-order uniformly accurate on the layer component of the solution. Motivated by this fact, we present a further improvement of the numerical method.

Convection–diffusion↗
DOE OSTI2022

Towards sharp error analysis of extended Lagrangian molecular dynamics

The extended Lagrangian molecular dynamics (XLMD) method provides a useful framework for reducing the computational cost of a class of molecular dynamics simulations with constrained latent variables. The XLMD method relaxes the constraints by introducing a fictitious mass ε for the latent variables and solving a set of singularly perturbed ordinary differential equations. While favorable numerical performance of XLMD has been demonstrated in several different contexts in the past decade, mathematical analysis of the method remains scarce. Here, we propose the first error analysis of the XLMD method in the context of a classical polarizable force field model. While the dynamics with respect to the atomic degrees of freedom are general and nonlinear, the key mathematical simplification of the polarizable force field model is that the constraints on the latent variables are given by a linear system of equations. We prove that when the initial value of the latent variables is compatible in a sense that we define, XLMD converges as the fictitious mass ε is made small with $\mathscr{O}$(ε) error for the atomic degrees of freedom and with $\mathscr{O}$($\sqrt{ε}$) error for the latent variables, when the dimension of the latent variable d' is 1. Furthermore, when the initial value of the latent variables is improved to be optimally compatible in a certain sense, we prove that the convergence rate can be improved to $\mathscr{O}$(ε) for the latent variables as well. Numerical results verify that both estimates are sharp not only for d'=1, but also for arbitrary d'. In the setting of general d', we do obtain convergence, but with the non-sharp rate of $\mathscr{O}$($\sqrt{ε}$) for both the atomic and latent variables.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

An efficient quantum circuit for block encoding a pairing Hamiltonian

We present an efficient quantum circuit for block encoding a pairing Hamiltonian often studied in nuclear physics. Our block encoding scheme does not require mapping the creation and annihilation operators to the Pauli operators and representing the Hamiltonian as a linear combination of unitaries. Instead, we show how to encode the Hamiltonian directly using controlled swap operations. We analyze the gate complexity of the block encoding circuit and show that it scales polynomially with respect to the number of qubits required to represent a quantum state associated with the pairing Hamiltonian. We also show how the block encoding circuit can be combined with the quantum singular value transformation to construct an efficient quantum circuit for approximating the density of states of a pairing Hamiltonian. The techniques presented can be extended to encode more general second-quantized Hamiltonians.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Causal horizons, geodesic completeness and stability in slow contraction cosmology

We show that cosmological models with a semi-infinite phase of slow contraction (ekpyrosis) possess a combination of properties that can address several fundamental problems in cosmology, otherwise faced in contracting de Sitter phases or standard big bang expansion. In particular, flat or open slow contraction admits a stable past attractor that asymptotes to Minkowski space and is past geodesically complete, as well as a stable, flat, homogeneous, and isotropic future attractor with negligible Weyl curvature (and, therefore, negligible gravitational entropy). In bouncing cosmologies, this contracting attractor is terminated by a smooth, non-singular bounce that transforms the attractor properties at the end of contraction into the initial conditions for the subsequent expanding phase. Cosmologies incorporating a slow contraction phase have no particle horizon and therefore avoid the causal horizon problem. The past Minkowski attractor also generates an initial spectrum of vacuum-like quantum fluctuations on all wavelengths. Moreover, because the averaged expansion rate along past-directed geodesics is non-positive, models incorporating a semi-infinite phase of slow contraction also evade the Borde–Guth–Vilenkin theorem. By contrast, contracting de Sitter space possesses a finite particle horizon and becomes unstable in the presence of scalar fields, matter, or radiation.

Cosmology and Nongalactic Astrophysics (astro-ph.C↗
DOE OSTI2020

A holey cavity for single-transducer 3D ultrasound imaging with physical optimization

Within the compressive sensing (CS) framework, one effective way to increase the likelihood of successful signal reconstruction is to employ random processes in the construction of the sensing matrix. This work presents a 3D holey cavity, with diverse frequency modes, to spectrally code, that is, randomize, the ultrasound wave fields. The simulated results show that the use of such a cavity enables imaging simple or complex targets, such as spheres or the letter E, by only a single transceiver—something that is not possible without the use of a coding structure like the cavity. The effect of noise on imaging results and the size of the targets on the first-order Born approximation (BA) are also investigated. Moreover, this study attempts to optimize the cavity, based on a single numerical metric, such as the sum of singular values (SSV) or mutual coherence (MC). Yet, it will be shown that neither of these metrics can consistently compare the norm-one imaging performance between two cavities of different materials or hole sizes. This leaves finding a quantitative metric for these kinds of optimizations an open problem.

42 ENGINEERING↗
DOE OSTI2024

A flexible class of priors for orthonormal matrices with basis function-specific structure

Statistical modeling of high-dimensional matrix-valued data motivates the use of a low-rank representation that simultaneously summarizes key characteristics of the data and enables dimension reduction. Low-rank representations commonly factor the original data into the product of orthonormal basis functions and weights, where each basis function represents an independent feature of the data. However, the basis functions in these factorizations are typically computed using algorithmic methods that cannot quantify uncertainty or account for basis function correlation structure a priori. While there exist Bayesian methods that allow for a common correlation structure across basis functions, empirical examples motivate the need for basis function-specific dependence structure. We propose a prior distribution for orthonormal matrices that can explicitly model basis function-specific structure. The prior is used within a general probabilistic model for singular value decomposition to conduct posterior inference on the basis functions while accounting for measurement error and fixed effects. We discuss how the prior specification can be used for various scenarios and demonstrate favorable model properties through synthetic data examples. Finally, we apply our method to two-meter air temperature data from the Pacific Northwest, enhancing our understanding of the Earth system’s internal variability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

High-Fidelity Large-Signal Order Reduction Approach for Composite Load Model

With the increasing penetration of electronic loads and distributed energy resources, conventional load models cannot capture their dynamics. Therefore, a new comprehensive composite load model is developed by Western Electricity Coordinating Council (WECC). However, this model is a complex high-order non-linear system with multi-time-scale property, which poses challenges on power system studies with heavy computational burden. In order to reduce the model complexity, the authors firstly develop a large-signal order reduction (LSOR) method using singular perturbation theory. In this method, the fast dynamics are integrated into the slow ones to preserve transient characteristics of the former. Then, accuracy assessment conditions are proposed and embedded into the LSOR to improve and guarantee the accuracy of reduced-order model. Finally, the reduced-order WECC composite load model is derived by using the proposed algorithm. Overall, simulation results show that the reduced-order large-signal model significantly alleviates the computational burden while maintaining similar dynamic responses as the original composite load model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Operator-level quantum acceleration of non-logconcave sampling

Sampling from probability distributions of the form 𝝈 ∝ e −𝜷V , where V is a continuous potential, is a fundamental task across physics, chemistry, biology, computer science, and statistics. However, when V is nonconvex, the resulting distribution becomes non-logconcave, and classical methods such as Langevin dynamics often exhibit poor performance. We introduce a quantum algorithm that provably accelerates a broad class of continuous-time sampling dynamics. For Langevin dynamics, our method encodes the target Gibbs measure into the amplitudes of aquantum state, identified as the kernel of a block matrix derived from a factorization of the Witten Laplacian operator. This connection enables Gibbs sampling via singular value thresholding and yields up to a quartic quantum speedup over best-knownclassical Langevin-based methods in the non-logconcave setting. Building on this framework, we further develop the first quantum algorithm that accelerates replica exchange Langevin diffusion, a widely used method for sampling from complex, rugged energy landscapes.

97 MATHEMATICS AND COMPUTING↗