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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 811 records · Page 45

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Linear, multivariable robust control with a mu perspective

The structured singular value is a linear algebra tool developed to study a particular class of matrix perturbation problems arising in robust feedback control of multivariable systems. These perturbations are called linear fractional, and are a natural way to model many types of uncertainty in linear systems, including state-space parameter uncertainty, multiplicative and additive unmodeled dynamics uncertainty, and coprime factor and gap metric uncertainty. The structured singular value theory provides a natural extension of classical SISO robustness measures and concepts to MIMO systems. The structured singular value analysis, coupled with approximate synthesis methods, make it possible to study the tradeoff between performance and uncertainty that occurs in all feedback systems. In MIMO systems, the complexity of the spatial interactions in the loop gains make it difficult to heuristically quantify the tradeoffs that must occur. This paper examines the role played by the structured singular value (and its computable bounds) in answering these questions, as well as its role in the general robust, multivariable control analysis and design problem.

Packard, Andy↗

NASA NTRS1993

The Generalized Legendre-Clebsch Condition on state/control constrained arcs

An extension of the Generalized Legendre-Clebsch Condition is obtained for problems with singular control along arcs with active state or control constraints. This is achieved by first transforming the Accessory Minimum Problem associated with constrained singular arcs into an unconstrained singular, linear quadratic problem. In a second step whose theorems and proofs are largely based on Goh's work (1966), necessary conditions are derived for such singular linear quadratic problems to yield nonnegative cost.

Seywald, Hans↗

NASA NTRS1993

Transverse cracking and stiffness reduction in composite laminates

A study of transverse cracking mechanism in composite laminates is presented using a singular hybrid finite element model. The model provides the global structural response as well as the precise local crack-tip stress fields. An elasticity basis for the problem is established by employing Lekhnitskii's complex variable potentials and method of eigenfunction expansion. Stress singularities associated with the transverse crack are obtained by decomposing the deformation into the symmetric and antisymmetric modes and proper boundary conditions. A singular hybrid element is thereby formulated based on the variational principle of a modified hybrid functional to incorporate local crack singularities. Axial stiffness reduction due to transverse cracking is studied. The results are shown to be in very good agreement with the existing experimental data. Comparison with simple shear lag analysis is also given. The effects of stress intensity factors and strain energy density on the increase of crack density are analyzed. The results reveal that the parameters approach definite limits when crack densities are saturated, an evidence of the existence of characteristic damage state.

Yuan, F. G.↗

NASA NTRS1993

Robust, nonlinear, high angle-of-attack control design for a supermaneuverable vehicle

High angle-of-attack flight control laws are developed for a supermaneuverable fighter aircraft. The methods of dynamic inversion and structured singular value synthesis are combined into an approach which addresses both the nonlinearity and robustness problems of flight at extreme operating conditions. The primary purpose of the dynamic inversion control elements is to linearize the vehicle response across the flight envelope. Structured singular value synthesis is used to design a dynamic controller which provides robust tracking to pilot commands. The resulting control system achieves desired flying qualities and guarantees a large margin of robustness to uncertainties for high angle-of-attack flight conditions. The results of linear simulation and structured singular value stability analysis are presented to demonstrate satisfaction of the design criteria. High fidelity nonlinear simulation results show that the combined dynamics inversion/structured singular value synthesis control law achieves a high level of performance in a realistic environment.

Adams, Richard J.↗

NASA NTRS1995

Study of nonlinear MHD equations governing the wave propagation in twisted coronal loops

The solar corona, modelled by a low beta, resistive plasma slab, sustains MHD wave propagations due to shearing footpoint motions in the photosphere. By using a numerical algorithm the excitation and nonlinear development of MHD waves in twisted coronal loops are studied. The plasma responds to the footpoint motion by sausage waves if there is no twist. The twist in the magnetic field of the loop destroys initially developed sausage-like wave modes and they become kinks. The transition from sausage to kink modes is analyzed. The twist brings about mode degradation producing high harmonics and this generates more complex fine structures. This can be attributed to several local extrema in the perturbed velocity profiles. The Alfven wave produces remnants of the ideal 1/x singularity both for zero and non-zero twist and this pseudo-singularity becomes less pronounced for larger twist. The effect of nonlinearity is clearly observed by changing the amplitude of the driver by one order of magnitude. The magnetosonic waves also exhibit smoothed remnants of ideal logarithmic singularities when the frequency of the driver is correctly chosen. This pseudo-singularity for fast waves is absent when the coronal loop does not undergo any twist but becomes pronounced when twist is included. On the contrary, it is observed for slow waves even if there is no twist. Increasing the twist leads to a higher heating rate of the loop. The larger twist shifts somewhat uniformly distributed heating to layers inside the slab corresponding to peaks in the magnetic field strength.

Parhi, S.↗

NASA NTRS1997

Exponential Approximations Using Fourier Series Partial Sums

The problem of accurately reconstructing a piece-wise smooth, 2(pi)-periodic function f and its first few derivatives, given only a truncated Fourier series representation of f, is studied and solved. The reconstruction process is divided into two steps. In the first step, the first 2N + 1 Fourier coefficients of f are used to approximate the locations and magnitudes of the discontinuities in f and its first M derivatives. This is accomplished by first finding initial estimates of these quantities based on certain properties of Gibbs phenomenon, and then refining these estimates by fitting the asymptotic form of the Fourier coefficients to the given coefficients using a least-squares approach. It is conjectured that the locations of the singularities are approximated to within O(N(sup -M-2), and the associated jump of the k(sup th) derivative of f is approximated to within O(N(sup -M-l+k), as N approaches infinity, and the method is robust. These estimates are then used with a class of singular basis functions, which have certain 'built-in' singularities, to construct a new sequence of approximations to f. Each of these new approximations is the sum of a piecewise smooth function and a new Fourier series partial sum. When N is proportional to M, it is shown that these new approximations, and their derivatives, converge exponentially in the maximum norm to f, and its corresponding derivatives, except in the union of a finite number of small open intervals containing the points of singularity of f. The total measure of these intervals decreases exponentially to zero as M approaches infinity. The technique is illustrated with several examples.

Banerjee, Nana S.↗

NASA NTRS2001

A Highly Accurate Technique for the Treatment of Flow Equations at the Polar Axis in Cylindrical Coordinates using Series Expansions

Numerical methods for solving the flow equations in cylindrical or spherical coordinates should be able to capture the behavior of the exact solution near the regions where the particular form of the governing equations is singular. In this work we focus on the treatment of these numerical singularities for finite-differences methods by reinterpreting the regularity conditions developed in the context of pseudo-spectral methods. A generally applicable numerical method for treating the singularities present at the polar axis, when nonaxisymmetric flows are solved in cylindrical, coordinates using highly accurate finite differences schemes (e.g., Pade schemes) on non-staggered grids, is presented. Governing equations for the flow at the polar axis are derived using series expansions near r=0. The only information needed to calculate the coefficients in these equations are the values of the flow variables and their radial derivatives at the previous iteration (or time) level. These derivatives, which are multi-valued at the polar axis, are calculated without dropping the accuracy of the numerical method using a mapping of the flow domain from (0,R)*(0,2pi) to (-R,R)*(0,pi), where R is the radius of the computational domain. This allows the radial derivatives to be evaluated using high-order differencing schemes (e.g., compact schemes) at points located on the polar axis. The proposed technique is illustrated by results from simulations of laminar-forced jets and turbulent compressible jets using large eddy simulation (LES) methods. In term of the general robustness of the numerical method and smoothness of the solution close to the polar axis, the present results compare very favorably to similar calculations in which the equations are solved in Cartesian coordinates at the polar axis, or in which the singularity is removed by employing a staggered mesh in the radial direction without a mesh point at r=0, following the method proposed recently by Mohseni and Colonius (1). Extension of the method described here for incompressible flows or for any other set of equations that are solved on a non-staggered mesh in cylindrical or spherical coordinates with finite-differences schemes of various level of accuracy is immediate.

Constantinescu, George S.↗

NASA NTRS1998

Scaling Properties of Particle Density Fields Formed in Simulated Turbulent Flows

Direct numerical simulations (DNS) of particle concentrations in fully developed 3D turbulence were carried out in order to study the nonuniform structure of the particle density field. Three steady-state turbulent fluid fields with Taylor microscale Reynolds numbers (Re(sub lambda)) of 40, 80 and 140 were generated by solving the Navier-Stokes equations with pseudospectral methods. Large scale forcing was used to drive the turbulence and maintain temporal stationarity. The response of the particles to the fluid was parameterized by the particle Stokes number St, defined as the ratio of the particle's stopping time to the mean period of eddies on the Kolmogorov scale (eta). In this paper, we consider only passive particles optimally coupled to these eddies (St approx. = 1) because of their tendency to concentrate more than particles with lesser or greater St values. The trajectories of up to 70 million particles were tracked in the equilibrated turbulent flows until the particle concentration field reached a statistically stationary state. The nonuniform structure of the concentration fields was characterized by the multifractal singularity spectrum, f(alpha), derived from measures obtained after binning particles into cells ranging from 2(eta) to 15(eta) in size. We observed strong systematic variations of f(alpha) across this scale range in all three simulations and conclude that the particle concentration field is not statistically self similar across the scale range explored. However, spectra obtained at the 2(eta), 4(eta), and 8(eta) scales of each flow case were found to be qualitatively similar. This result suggests that the local structure of the particle concentration field may be flow-Independent. The singularity spectra found for 2n-sized cells were used to predict concentration distributions in good agreement with those obtained directly from the particle data. This Singularity spectrum has a shape similar to the analogous spectrum derived for the inertial-range energy dissipation fields of experimental turbulent flows at Re(sub lambda) = 110 and 1100. Based on this agreement, and the expectation that both dissipation and particle concentration are controlled by the same cascade process, we hypothesize that singularity spectra similar to the ones found in this work provide a good characterization of the spatially averaged statistical properties of preferentially concentrated particles in higher Re(sub lambda) turbulent flows.

Hogan, Robert C.↗

NASA NTRS1957

Wall Interference in a Perforated Wind Tunnel

The theory of cascades, made up of a series of flat plates placed one behind the other, is extended to the case where the impinging stream is not uniform, and the deduced properties of this cascade-flow are then applied to the study of the wall interference between such as cascade-like boundary and a vortex-source type of singularity. It is shown that the induced velocities, produced by the presence of such a wall, are equal to what is obtained by action of a suitably chosen 'reflected' singularity situated on one side of the wall, together with the action of another suitably chosen 'transmitted' singularity placed on the other side. The concepts of a reflection factor and a transmission factor are introduced to characterize various types of tunnel boundary, whether this is composed of solid wall, open and closed sequences, or just a free fluid surface. These ideas are then extended to cover the situation in which a pair of such walls are allowed to coalesce, especially in the event that one of the walls consists of open and closed portions and the other is a fluid surface. This latter particular combination of boundaries is called a perforated wall. Finally, the interference arising from an arbitrary general singularity placed symmetrically in between two such perforated walls is analyzed, and these results are then applied to the determination of the particular kinds of wall geometry that will produce no interference effects in the case of a slender lifting wing and also in the case of a symmetric profile having a finite thickness.

Brescia, Riccardo↗

NASA NTRS2000

3D Higher Order Modeling in the BEM/FEM Hybrid Formulation

Higher order divergence- and curl-conforming bases have been shown to provide significant benefits, in both convergence rate and accuracy, in the 2D hybrid finite element/boundary element formulation (P. Fink and D. Wilton, National Radio Science Meeting, Boulder, CO, Jan. 2000). A critical issue in achieving the potential for accuracy of the approach is the accurate evaluation of all matrix elements. These involve products of high order polynomials and, in some instances, singular Green's functions. In the 2D formulation, the use of a generalized Gaussian quadrature method was found to greatly facilitate the computation and to improve the accuracy of the boundary integral equation self-terms. In this paper, a 3D, hybrid electric field formulation employing higher order bases and higher order elements is presented. The improvements in convergence rate and accuracy, compared to those resulting from lower order modeling, are established. Techniques developed to facilitate the computation of the boundary integral self-terms are also shown to improve the accuracy of these terms. Finally, simple preconditioning techniques are used in conjunction with iterative solution procedures to solve the resulting linear system efficiently. In order to handle the boundary integral singularities in the 3D formulation, the parent element- either a triangle or rectangle-is subdivided into a set of sub-triangles with a common vertex at the singularity. The contribution to the integral from each of the sub-triangles is computed using the Duffy transformation to remove the singularity. This method is shown to greatly facilitate t'pe self-term computation when the bases are of higher order. In addition, the sub-triangles can be further divided to achieve near arbitrary accuracy in the self-term computation. An efficient method for subdividing the parent element is presented. The accuracy obtained using higher order bases is compared to that obtained using lower order bases when the number of unknowns is approximately equal. Also, convergence rates obtained using higher order bases are compared to those obtained with lower order bases for selected sample

Fink, P. W.↗

NASA NTRS2011

Normalization of Gravitational Acceleration Models

Unlike the uniform density spherical shell approximations of Newton, the con- sequence of spaceflight in the real universe is that gravitational fields are sensitive to the nonsphericity of their generating central bodies. The gravitational potential of a nonspherical central body is typically resolved using spherical harmonic approximations. However, attempting to directly calculate the spherical harmonic approximations results in at least two singularities which must be removed in order to generalize the method and solve for any possible orbit, including polar orbits. Three unique algorithms have been developed to eliminate these singularities by Samuel Pines [1], Bill Lear [2], and Robert Gottlieb [3]. This paper documents the methodical normalization of two1 of the three known formulations for singularity-free gravitational acceleration (namely, the Lear [2] and Gottlieb [3] algorithms) and formulates a general method for defining normalization parameters used to generate normalized Legendre Polynomials and ALFs for any algorithm. A treatment of the conventional formulation of the gravitational potential and acceleration is also provided, in addition to a brief overview of the philosophical differences between the three known singularity-free algorithms.

Eckman, Randy A.↗

NASA NTRS2016

Normalization and Implementation of Three Gravitational Acceleration Models

Unlike the uniform density spherical shell approximations of Newton, the consequence of spaceflight in the real universe is that gravitational fields are sensitive to the asphericity of their generating central bodies. The gravitational potential of an aspherical central body is typically resolved using spherical harmonic approximations. However, attempting to directly calculate the spherical harmonic approximations results in at least two singularities that must be removed to generalize the method and solve for any possible orbit, including polar orbits. Samuel Pines, Bill Lear, and Robert Gottlieb developed three unique algorithms to eliminate these singularities. This paper documents the methodical normalization of two of the three known formulations for singularity-free gravitational acceleration (namely, the Lear and Gottlieb algorithms) and formulates a general method for defining normalization parameters used to generate normalized Legendre polynomials and Associated Legendre Functions (ALFs) for any algorithm. A treatment of the conventional formulation of the gravitational potential and acceleration is also provided, in addition to a brief overview of the philosophical differences between the three known singularity-free algorithms.

Eckman, Randy A.↗

DOE OSTI2019

Interaction of Harmonic Longitudinal Shear Waves with V-Shaped Inclusions

We solve the problem of determination of the stressed state in the vicinity of a tunnel rigid inclusion whose cross section consists of two segments originating out from the same point. The inclusion is placed in an infinite elastic medium with propagating plane harmonic longitudinal shear waves. The problem is reduced to the solution of a system of two singular integral equations with fixed singularities. For the approximate solution of this system, we apply a numerical method that takes into account the true asymptotics of the unknown functions and is based on the use of special quadrature formulas for the evaluation of singular integrals.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2014

Excited-state quantum phase transitions in systems with two degrees of freedom: Level density, level dynamics, thermal properties

Quantum systems with a finite number of freedom degrees f develop robust singularities in the energy spectrum of excited states as the system’s size increases to infinity. We analyze the general form of these singularities for low f, particularly f=2, clarifying the relation to classical stationary points of the corresponding potential. Signatures in the smoothed energy dependence of the quantum state density and in the flow of energy levels with an arbitrary control parameter are described along with the relevant thermodynamical consequences. The general analysis is illustrated with specific examples of excited-state singularities accompanying the first-order quantum phase transition. -- Highlights: •ESQPTs found in infinite-size limit of systems with low numbers of freedom degrees f. •ESQPTs related to non-analytical evolutions of classical phase–space properties. •ESQPT signatures analyzed for general f, particularly f=2, extending known case f=1. •ESQPT signatures identified in smoothened density and flow of energy spectrum. •ESQPTs shown to induce a new type of thermodynamic anomalies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2014

Electromagnetic fields and potentials generated by massless charged particles

We provide for the first time the exact solution of Maxwell’s equations for a massless charged particle moving on a generic trajectory at the speed of light. In particular we furnish explicit expressions for the vector potential and the electromagnetic field, which were both previously unknown, finding that they entail different physical features for bounded and unbounded trajectories. With respect to the standard Liénard–Wiechert field the electromagnetic field acquires singular δ-like contributions whose support and dimensionality depend crucially on whether the motion is (a) linear, (b) accelerated unbounded, (c) accelerated bounded. In the first two cases the particle generates a planar shock-wave-like electromagnetic field traveling along a straight line. In the second and third cases the field acquires, in addition, a δ-like contribution supported on a physical singularity-string attached to the particle. For generic accelerated motions a genuine radiation field is also present, represented by a regular principal-part type distribution diverging on the same singularity-string. - Highlights: • First exact solution of Maxwell’s equations for massless charges in arbitrary motion. • Explicit expressions of electromagnetic fields and potentials. • Derivations are rigorous and based on distribution theory. • The form of the field depends heavily on whether the motion is bounded or unbounded. • The electromagnetic field contains unexpected Dirac-delta-function contributions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2017

Relativistic model for anisotropic strange stars

Highlights: We find exact solutions for the Einstein field equations for anisotropic strange stars. Maximum anisotropy at the surface is an inherent property for the singularity free compact stars. Maximum value for the mass and radius is presented for the corresponding values of the bag constants. Using observed value of mass different physical parameters for the strange stars are predicted. - Abstract: In this article, we attempt to find a singularity free solution of Einstein’s field equations for compact stellar objects, precisely strange (quark) stars, considering Schwarzschild metric as the exterior spacetime. To this end, we consider that the stellar object is spherically symmetric, static and anisotropic in nature and follows the density profile given by Mak and Harko (2002) , which satisfies all the physical conditions. To investigate different properties of the ultra-dense strange stars we have employed the MIT bag model for the quark matter. Our investigation displays an interesting feature that the anisotropy of compact stars increases with the radial coordinate and attains its maximum value at the surface which seems an inherent property for the singularity free anisotropic compact stellar objects. In this connection we also perform several tests for physical features of the proposed model and show that these are reasonably acceptable within certain range. Further, we find that the model is consistent with the energy conditions and the compact stellar structure is stable with the validity of the TOV equation and Herrera cracking concept. For the masses below the maximum mass point in mass vs radius curve the typical behavior achieved within the framework of general relativity. We have calculated the maximum mass and radius of the strange stars for the three finite values of bag constant B{sub g}.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

New scattering features in non-Hermitian space fractional quantum mechanics

Highlights: • Space fraction quantum mechanics (SFQM) has been investigated with complex delta and barrier potentials. • New scattering features due to these potentials are obtained in SFQM. • Spectral Singular (SS) points are blue/red shifted when Levy index is decreased. • The oscillations in scattering coefficients around SS have interesting behaviors in SFQM. • Deeper relations between the scattering properties in standard QM and those of in SFQM are observed. - Abstract: The spectral singularity (SS) and coherent perfect absorption (CPA) have been extensively studied over the last one and half decade for different non-Hermitian potentials in non-Hermitian standard quantum mechanics (SQM) governed by Schrödinger equation. In the present work we explore these scattering features in the domain of non-Hermitian space fractional quantum mechanics (SFQM) governed by fractional Schrödinger equation which is characterized by Levy index α (1<α≤2). We observe that non-Hermitian SFQM systems have more flexibility for SS and CPA and display some new features of scattering. For the delta potential V(x)=iρδ(x−x{sub 0}), ρ>0, the SS energy, E{sub ss}, is blue or red shifted with decreasing α depending on the strength of the potential. For complex rectangular barrier in non-Hermitian SQM, it is known that the reflection and transmission amplitudes are oscillatory near the spectral singular point. It is found that these oscillations eventually develop SS in non-Hermitian SFQM. The similar features are also reported for the case of CPA phenomena from complex rectangular barrier in non-Hermitian SFQM. These observations suggest a deeper relation between scattering features of non-Hermitian SQM and non-Hermitian SFQM.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

Applications of the worldline Monte Carlo formalism in quantum mechanics

Highlights: Universal approach to the numerical computation of the Euclidean path integral. Inspired by recent work in relativistic quantum field theory. Here adapted to non-relativistic quantum mechanics. Worked out for the computation of propagators and ground-state energies. Special smoothing procedure for singular potentials. -- Abstract: In recent years efficient algorithms have been developed for the numerical computation of relativistic single-particle path integrals in quantum field theory. Here, we adapt this “worldline Monte Carlo” approach to the standard problem of the numerical approximation of the non-relativistic path integral, resulting in a formalism whose characteristic feature is the fast, non-recursive generation of an ensemble of trajectories that is independent of the potential, and is thus universally applicable. The numerical implementation discretises the trajectories with respect to their time parametrisation but maintains a continuous spatial domain. In the case of singular potentials, the discretised action gets adapted to the singularity through a “smoothing” procedure. We show for a variety of examples (the harmonic oscillator in various dimensions, the modified Pöschl–Teller potential, delta-function potentials, the Coulomb and Yukawa potentials) that the method allows one to obtain fast and reliable estimates for the Euclidean propagator and use them in a certain time window suitable for extracting the ground state energy. Usually Monte Carlo calculations in quantum mechanics are performed using algorithms that involve a more specific adaptation to the potential, and we briefly compare our method with the most widely used such algorithms, which are of Metropolis type. As an aside, we apply it for studying the classical limit where nearly classical trajectories are expected to dominate in the path integral. We expect the advances made here to be useful also in the relativistic case.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗