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 523 records · Page 29

Issues and Methods Concerning the Evaluation of Hypersingular and Near-Hypersingular Integrals in BEM Formulations

It is known that higher order modeling of the sources and the geometry in Boundary Element Modeling (BEM) formulations is essential to highly efficient computational electromagnetics. However, in order to achieve the benefits of hIgher order basis and geometry modeling, the singular and near-singular terms arising in BEM formulations must be integrated accurately. In particular, the accurate integration of near-singular terms, which occur when observation points are near but not on source regions of the scattering object, has been considered one of the remaining limitations on the computational efficiency of integral equation methods. The method of singularity subtraction has been used extensively for the evaluation of singular and near-singular terms. Piecewise integration of the source terms in this manner, while manageable for bases of constant and linear orders, becomes unwieldy and prone to error for bases of higher order. Furthermore, we find that the singularity subtraction method is not conducive to object-oriented programming practices, particularly in the context of multiple operators. To extend the capabilities, accuracy, and maintainability of general-purpose codes, the subtraction method is being replaced in favor of the purely numerical quadrature schemes. These schemes employ singularity cancellation methods in which a change of variables is chosen such that the Jacobian of the transformation cancels the singularity. An example of the sin,oularity cancellation approach is the Duffy method, which has two major drawbacks: 1) In the resulting integrand, it produces an angular variation about the singular point that becomes nearly-singular for observation points close to an edge of the parent element, and 2) it appears not to work well when applied to nearly-singular integrals. Recently, the authors have introduced the transformation u(x(prime))= sinh (exp -1) x(prime)/Square root of ((y prime (exp 2))+ z(exp 2) for integrating functions of the form I = Integral of (lambda(r(prime))((e(exp -jkR))/(4 pi R) d D where A (r (prime)) is a vector or scalar basis function and R = Square root of( (x(prime)(exp2) + (y(prime)(exp2) + z(exp 2)) is the distance between source and observation points. This scheme has all of the advantages of the Duffy method while avoiding the disadvantages listed above. In this presentation we will survey similar approaches for handling singular and near-singular terms for kernels with 1/R(exp 2) type behavior, addressing potential pitfalls and offering techniques to efficiently handle special cases.

Fink, P. W.↗

Minimum-fuel rocket trajectories involving intermediate-thrust arcs

The optimal trajectories in the neighborhood of an optimal intermediate-thrust arc are investigated for the minimum-fuel orbit rendezvous problem with fixed specific impulse. Since such an arc is singular, the thrust acceleration magnitude being the singular control component, a second-variation analysis leads to the identification of a field of neighboring, singular arcs in a state space of dimension four rather than six, provided that a suitable Jacobi condition is met. A given neighboring initial six-dimensional state vector does not generally lie on a neighboring singular arc, and junction onto the appropriate singular arc must be accomplished by a short period of strong variations in the acceleration. The neighboring singular arc meets the final condition in 4 dimensions, rather than 6 dimensions, and rendezvous must be completed by another, terminal short period of strong variations in the acceleration. Implications for midcourse guidance near a singular arc are discussed.

Breakwell, J. V.↗

Steering Law Design for Redundant Single Gimbal Control Moment Gyro Systems

The correspondence between robotic manipulators and single gimbal Control Moment Gyro (CMG) systems was exploited to aid in the understanding and design of single gimbal CMG Steering laws. A test for null motion near a singular CMG configuration was derived which is able to distinguish between escapable and unescapable singular states. Detailed analysis of the Jacobian matrix null-space was performed and results were used to develop and test a variety of single gimbal CMG steering laws. Computer simulations showed that all existing singularity avoidance methods are unable to avoid Elliptic internal singularities. A new null motion algorithm using the Moore-Penrose pseudoinverse, however, was shown by simulation to avoid Elliptic type singularities under certain conditions. The SR-inverse, with appropriate null motion was proposed as a general approach to singularity avoidance, because of its ability to avoid singularities through limited introduction of torque error. Simulation results confirmed the superior performance of this method compared to the other available and proposed pseudoinverse-based Steering laws.

Bedrossian, Nazareth Sarkis↗

Magnetic flipping - Reconnection in three dimensions without null points

In three dimensions, magnetic reconnection may take place in a sheared magnetic field at any singular field line, where the nearby field has X-type topology in planes perpendicular to the field line and where an electric field is present parallel to the field line. In the ideal region around the singular line there will, in general, be singularities in the plasma flow and electric field, both at the singular line and at 'magnetic flipping layers', which are remnants of local magnetic separatrices. In the absence of a three-dimensional magnetic neutral point or null point, reconnection of field lines can still occur by a process of magnetic flipping, in which the plasma crosses the flipping layers but the field lines rapidly flip along them by magnetic diffusion. Depending on the boundary conditions, there may be two or four flipping layers which converge on the singular line. A boundary layer analysis of a flipping layer is given, in which the magnetic field parallel to the layer decreases as one crosses it while the plasma pressure (or magnetic pressure associated with the field along the singular line) increases. The width of the flipping layer decreases with distance from the singular line.

Priest, E. R.↗

Simple and Efficient Numerical Evaluation of Near-Hypersingular Integrals

Recently, significant progress has been made in the handling of singular and nearly-singular potential integrals that commonly arise in the Boundary Element Method (BEM). To facilitate object-oriented programming and handling of higher order basis functions, cancellation techniques are favored over techniques involving singularity subtraction. However, gradients of the Newton-type potentials, which produce hypersingular kernels, are also frequently required in BEM formulations. As is the case with the potentials, treatment of the near-hypersingular integrals has proven more challenging than treating the limiting case in which the observation point approaches the surface. Historically, numerical evaluation of these near-hypersingularities has often involved a two-step procedure: a singularity subtraction to reduce the order of the singularity, followed by a boundary contour integral evaluation of the extracted part. Since this evaluation necessarily links basis function, Green s function, and the integration domain (element shape), the approach ill fits object-oriented programming concepts. Thus, there is a need for cancellation-type techniques for efficient numerical evaluation of the gradient of the potential. Progress in the development of efficient cancellation-type procedures for the gradient potentials was recently presented. To the extent possible, a change of variables is chosen such that the Jacobian of the transformation cancels the singularity. However, since the gradient kernel involves singularities of different orders, we also require that the transformation leaves remaining terms that are analytic. The terms "normal" and "tangential" are used herein with reference to the source element. Also, since computational formulations often involve the numerical evaluation of both potentials and their gradients, it is highly desirable that a single integration procedure efficiently handles both.

Fink, Patrick W.↗

$T\overline{T}$, black holes and negative strings

String theory on AdS 3 has a solvable single-trace irrelevant deformation that is closely related to $T\overline{T}$. For one sign of the coupling, it leads to an asymptotically linear dilaton spacetime, and a corresponding Hagedorn spectrum. For the other, the resulting spacetime has a curvature singularity at a finite radial location, and an upper bound on the energies of states. Beyond the singularity, the signature of spacetime is flipped and there is an asymptotically linear dilaton boundary at infinity. We study the properties of black holes and fundamental strings in this spacetime, and find a sensible picture. The singularity does not give rise to a hard ultraviolet wall for excitations -one must include the region beyond it to understand the theory. The size of black holes diverges as their energy approaches the upper bound, as does the location of the singularity. Fundamental strings pass smoothly through the singularity, but if their energy is above the upper bound, their trajectories are singular. From the point of view of the boundary at infinity, this background can be thought of as a vacuum of Little String Theory which contains a large number of negative strings.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Redundancy Parameterization and Inverse Kinematics of 7-DOF Revolute Manipulators

Seven degree-of-freedom (DOF) robot arms have one redundant DOF for obstacle and singularity avoidance. This DOF does not change the end effector motion and must be parameterized to fully specify the joint angles for a given end effector pose. For 7-DOF revolute (7R) manipulators, we introduce the generalized shoulder-elbow-wrist (SEW) angle, a generalization of the widely-used conventional SEW angle but with an arbitrary reference direction function, along with Jacobian expressions and singularity analysis. Other redundancy parameterizations including the conventional SEW angle encounter an algorithmic singularity along a line in the workspace. We introduce a reference direction function choice called the stereographic SEW angle which has a singularity only along a half-line and can be placed out of reach. We prove such a singularity is unavoidable for any parameterization. Finally, using the general SEW angle and the subproblem decomposition method, we provide efficient singularity-robust inverse kinematics solutions for 7R manipulators which are often closed-form but may involve a 1D or (in general) 2D search. Search-based solutions may be converted to finding polynomial roots. Examples are available in a publicly accessible repository.

Kinematics↗

Holographic thermal correlators: a tale of Fuchsian ODEs and integration contours

We analyze real-time thermal correlation functions of conserved currents in holographic field theories using the grSK geometry, which provides a contour prescription for their evaluation. We demonstrate its efficacy, arguing that there are situations involving components of conserved currents, or derivative interactions, where such a prescription is, in fact, essential. To this end, we first undertake a careful analysis of the linearized wave equations in AdS black hole backgrounds and identify the branch points of the solutions as a function of (complexified) frequency and momentum. All the equations we study are Fuchsian with only regular singular points that for the most part are associated with the geometric features of the background. Special features, e.g., the appearance of apparent singular points at the horizon, whence outgoing solutions end up being analytic, arise at higher codimension loci in parameter space. Using the grSK geometry, we demonstrate that these apparent singularities do not correspond to any interesting physical features in higher-point functions. We also argue that the Schwinger-Keldysh collapse and KMS conditions, implemented by the grSK geometry, continue to hold even in the presence of such singularities. For charged black holes above a critical charge, we furthermore demonstrate that the energy density operator does not possess an exponentially growing mode, associated with ‘pole-skipping’, from one such apparent singularity. Our analysis suggests that the connection between the scrambling physics of black holes and energy transport has, at best, a limited domain of validity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Applications of the Landau bootstrap

We advocate a strategy of bootstrapping Feynman integrals from just knowledge of their singular behavior. This approach is complementary to other bootstrap programs, which exploit nonperturbative constraints such as unitarity, or amplitude-level constraints such as gauge invariance. We begin by studying where a Feynman integral can become singular, and the behavior it exhibits near these singularities. We then characterize the space of functions that we expect the integral to evaluate to, in order to formulate an appropriate ansatz. Finally, we derive constraints on where each singularity can appear in this ansatz, and use information about the expansion of the integral around singular points in order to determine the value of all remaining free coefficients. Throughout, we highlight how constraints that have previously only been derived for integrals with generic masses can be extended to integrals involving particles of equal or vanishing mass. We illustrate the effectiveness of this approach by bootstrapping a number of examples, including the four-point double box with a massive internal loop. Published by the American Physical Society 2025

Hannesdottir, Holmfridur S. (ORCID:000000025440208↗

An improved method for the aerodynamic analysis of wing-body-tail configurations in subsonic and supersonic flow. Part 1: Theory and application

A new method has been developed for calculating the pressure distribution and aerodynamic characteristics of wing-body-tail combinations in subsonic and supersonic potential flow. A computer program has been developed to perform the numerical calculations. The configuration surface is subdivided into a large number of panels, each of which contains an aerodynamic singularity distribution. A constant source distribution is used on the body panels, and a vortex distribution having a linear variation in the streamwise direction is used on the wing and tail panels. The normal components of velocity induced at specified control points by each singularity distribution are calculated and make up the coefficients of a system of linear equations relating the strengths of the singularities to the magnitude of the normal velocities. The singularity strengths which satisfy the boundary condition of tangential flow at the control points for a given Mach number and angle of attack are determined by solving this system of equations using an interactive procedure. Once the singularity strengths are known, the pressure coefficients are calculated, and the forces and moments acting on the configuration determined by numerical integration.

Woodward, F. A.↗

The Dirichlet problem for the two-dimensional Helmholtz equation for an open boundary

Development of a complete theory of the two-dimensional Dirichlet problem for an open boundary. It is shown that the solution of the Dirichlet problem for an open boundary requires the solution of a Fredholm integral equation of the first kind. Although a Fredholm integral equation of the first kind usually has no solution if the kernel is continuous, owing to the logarithmic singularity of the kernel, the equation in this case is converted to a singular integral equation with a Cauchy kernel. It is proven that the homogeneous adjoint equation of the singular integral equation has no nonzero solution. By virtue of this result, and with the aid of an existence theorem known in the theory of singular integral equations, the existence of solutions of the singular integral equation, and then of the unique solution of the Fredholm integral equation of the first kind is proved.

Hayashi, Y.↗

Elasto-plastic flow in cracked bodies using a new finite element model

Cracked geometries were studied by finite element techniques with the aid of a new special element embedded at the crack tip. This model seeked to accurately represent the singular stresses and strains associated with the elasto-plastic flow process. The present model was not restricted to a material type and did not predetermine a singularity. Rather the singularity was treated as an unknown. For each step of the incremental process the nodal degrees of freedom and the unknown singularity were found through minimization of an energy-like functional. The singularity and nodal degrees of freedom were determined by means of an iterative process.

Karabin, M. E., Jr.↗

On the solution of integral equations with a generalized cauchy kernel

In this paper a certain class of singular integral equations that may arise from the mixed boundary value problems in nonhomogeneous materials is considered. The distinguishing feature of these equations is that in addition to the Cauchy singularity, the kernels contain terms that are singular only at the end points. In the form of the singular integral equations adopted, the density function is a potential or a displacement and consequently the kernel has strong singularities of the form (t-x) sup-2, x sup n-2 (t+x) sup n, (n or = 2, 0x,tb). The complex function theory is used to determine the fundamental function of the problem for the general case and a simple numerical technique is described to solve the integral equation. Two examples from the theory of elasticity are then considered to show the application of the technique.

Kaya, A. C.↗

On the solution of integral equations with a generalized cauchy kernal

A certain class of singular integral equations that may arise from the mixed boundary value problems in nonhonogeneous materials is considered. The distinguishing feature of these equations is that in addition to the Cauchy singularity, the kernels contain terms that are singular only at the end points. In the form of the singular integral equations adopted, the density function is a potential or a displacement and consequently the kernal has strong singularities of the form (t-x)(-2), x(n-2) (t+x)(n), (n is = or 2, 0 x, t b). The complex function theory is used to determine the fundamental function of the problem for the general case and a simple numerical technique is described to solve the integral equation. Two examples from the theory of elasticity are then considered to show the application of the technique.

Kaya, A. C.↗

Controller design for the improvement of feedback properties: A method for tuning the weight in the evaluation function of the LQG theory

In the controller design of a linear time-invariant system, it is important to improve the feedback properties such as robust stability and sensitivity. In the multi-input multi-output case, these properties can be estimated by using the singular values of return difference matrix. The design method to obtain better singular-value-points is desired. How to tune the weight matrix of performance index and/or covariance matrix of noise in the Linear Quadratic Gaussian (LQC) theory to get a desired singular-value-plot was studied. First, the property of singular-value plots of return difference matrix of a system designed by LQG theory is examined from the viewpoint of tuning weight. Second, a distance between real singular-value-plots and a desired plot is defined, and the weight of performance index is numerically determined by quasi-Newton method so that the distance is minimized.

Saeki, Masami↗

Steering law design for redundant single-gimbal control moment gyroscopes

Two steering laws are presented for single-gimbal control moment gyroscopes. An approach using the Moore-Penrose pseudoinverse with a nondirectional null-motion algorithm is shown by example to avoid internal singularities for unidirectional torque commands, for which existing algorithms fail. Because this is still a tangent-based approach, however, singularity avoidance cannot be guaranteed. The singularity robust inverse is introduced as an alternative to the pseudoinverse for computing torque-producing gimbal rates near singular states. This approach, coupled with the nondirectional null algorithm, is shown by example to provide better steering law performance by allowing torque errors to be produced in the vicinity of singular states.

Bedrossian, Nazareth S.↗