Search NASA⌕ Search

SEARCH · Search NASA

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 55 records · Page 3

An Optimal Orthogonal Decomposition Method for Kalman Filter-Based Turbofan Engine Thrust Estimation

A new linear point design technique is presented for the determination of tuning parameters that enable the optimal estimation of unmeasured engine outputs, such as thrust. The engine s performance is affected by its level of degradation, generally described in terms of unmeasurable health parameters related to each major engine component. Accurate thrust reconstruction depends on knowledge of these health parameters, but there are usually too few sensors to be able to estimate their values. In this new technique, a set of tuning parameters is determined that accounts for degradation by representing the overall effect of the larger set of health parameters as closely as possible in a least-squares sense. The technique takes advantage of the properties of the singular value decomposition of a matrix to generate a tuning parameter vector of low enough dimension that it can be estimated by a Kalman filter. A concise design procedure to generate a tuning vector that specifically takes into account the variables of interest is presented. An example demonstrates the tuning parameters ability to facilitate matching of both measured and unmeasured engine outputs, as well as state variables. Additional properties of the formulation are shown to lend themselves well to diagnostics.

Litt, Jonathan S.↗

Visualization of 3-D tensor fields

Second-order tensor fields have applications in many different areas of physics, such as general relativity and fluid mechanics. The wealth of multivariate information in tensor fields makes them more complex and abstract than scalar and vector fields. Visualization is a good technique for scientists to gain new insights from them. Visualizing a 3-D continuous tensor field is equivalent to simultaneously visualizing its three eigenvector fields. In the past, research has been conducted in the area of two-dimensional tensor fields. It was shown that degenerate points, defined as points where eigenvalues are equal to each other, are the basic singularities underlying the topology of tensor fields. Moreover, it was shown that eigenvectors never cross each other except at degenerate points. Since we live in a three-dimensional world, it is important for us to understand the underlying physics of this world. In this report, we describe a new method for locating degenerate points along with the conditions for classifying them in three-dimensional space. Finally, we discuss some topological features of three-dimensional tensor fields, and interpret topological patterns in terms of physical properties.

Hesselink, L.↗

Multifractal characterizations of nonstationarity and intermittency in geophysical fields: Observed, retrieved, or simulated

Geophysical data rarely show any smoothness at any scale, and this often makes comparison with theoretical model output difficult. However, highly fluctuating signals and fractal structures are typical of open dissipative systems with nonlinear dynamics, the focus of most geophysical research. High levels of variability are excited over a large range of scales by the combined actions of external forcing and internal instability. At very small scales we expect geophysical fields to be smooth, but these are rarely resolved with available instrumentation or simulation tools; nondifferentiable and even discontinuous models are therefore in order. We need methods of statistically analyzing geophysical data, whether measured in situ, remotely sensed or even generated by a computer model, that are adapted to these characteristics. An important preliminary task is to define statistically stationary features in generally nonstationary signals. We first discuss a simple criterion for stationarity in finite data streams that exhibit power law energy spectra and then, guided by developments in turbulence studies, we advocate the use of two ways of analyzing the scale dependence of statistical information: singular measures and qth order structure functions. In nonstationary situations, the approach based on singular measures seeks power law behavior in integrals over all possible scales of a nonnegative stationary field derived from the data, leading to a characterization of the intermittency in this (gradient-related) field. In contrast, the approach based on structure functions uses the signal itself, seeking power laws for the statistical moments of absolute increments over arbitrarily large scales, leading to a characterization of the prevailing nonstationarity in both quantitative and qualitative terms. We explain graphically, step by step, both multifractal statistics which are largely complementary to each other. The geometrical manifestations of nonstationarity and intermittency, 'roughness' and 'sparseness', respectively, are illustrated and the associated analytical (differentiability and continuity) properties are discussed. As an example, the two techniques are applied to a series of recent measurements of liquid water distributions inside marine stratocumulus decks; these are found to be multifractal over scales ranging from approximately 60 m to approximately 60 km. Finally, we define the 'mean multifractal plane' and show it to be a simple yet comprehensive tool with many applications including data intercomparison, (dynamical or stochastic) model and retrieval validations.

Davis, Anthony↗

The topology of fractal universes

It is shown how the genus per unit volume of isodensity surfaces in general nonlinear universes is related to the entire hierarchy of correlation functions. The general relation between the correlation function, the probability distribution of densities at several points, and the probability distributions of density and its derivatives at a point are given. Formulas for the area and genus per unit volume of isodensity surfaces are presented. As an application, after first reviewing the case of Gaussian fields, analytic results are reported for one particular example of a thoroughly nonlinear universe, Mandelbrot's Rayleigh-Levy random-walk fractal. While this fractal bears little resemblance to the real universe of galaxies, it possesses the singular and theoretically interesting property that in it cluster-cluster correlations are identically equal to galaxy-galaxy correlations to all orders.

Hamilton, A. J. S.↗

To perform a gyro test of general relativity in a satellite and develop associated control technology

A satellite configuration having two gyroscopes with axes parallel to the boresight of a telescope and two at right angles to the telescope and approximately parallel and perpendicular to the earth's axis is proposed for measuring geodetic precessions due to the earth's motion about the sun, higher order geodetic terms calculated from the earth's quadrapole mass moment (0.010 arc-sec/year in a 400 nautical mile polar orbit), and deflection by the sun of the starlight signal for the reference telescope. Data from the experiment also contain large periodic signals due to the annual and orbital aberrations of starlight which are useful in providing a built in reference signal of known amplitude for scaling the relativity signals, and should yield a singularly precise measurement of the parallax of the reference star. The development of the gyroscope and its readout system are discussed, as well as signal integration, drag-free control, and attitude control.

Fairbank, W. M.↗

Interlaminar stresses in composite laminates: A perturbation analysis

A general method of solution for an elastic balanced symmetric composite laminate subject to a uniaxial extension was developed based upon a perturbation analysis of a limiting free body containing an interfacial plane. The solution satisfies more physical requirements and boundary conditions than previous investigations, and predicts smooth continuous interlaminar stresses with no instabilities. It determines the finite maximum intensity for the interlaminar normal stress in all laminates, provides mathematical evidences for the singular stresses in angle-ply laminates, suggests the need for the experimental determination of an important problem parameter, and introduces a viable means for solving related problems of practical interest.

Hsu, P. W.↗

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.↗

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.↗

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↗

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.↗

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.↗

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.↗

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↗

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.↗

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↗

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.↗

Contrasting Indian Ocean SST Variability With and Without ENSO Influence: A Coupled Atmosphere-Ocean GCM Study

In this study, we perform experiments with a coupled atmosphere-ocean general circulation model (CGCM) to examine ENSO's influence on the interannual sea surface temperature (SST) variability of the tropical Indian Ocean. The control experiment includes both the Indian and Pacific Oceans in the ocean model component of the CGCM (the Indo-Pacific Run). The anomaly experiment excludes ENSOs influence by including only the Indian Ocean while prescribing monthly-varying climatological SSTs for the Pacific Ocean (the Indian-Ocean Run). In the Indo-Pacific Run, an oscillatory mode of the Indian Ocean SST variability is identified by a multi-channel singular spectral analysis (MSSA). The oscillatory mode comprises two patterns that can be identified with the Indian Ocean Zonal Mode (IOZM) and a basin-wide warming/cooling mode respectively. In the model, the IOZM peaks about 3-5 months after ENSO reaches its maximum intensity. The basin mode peaks 8 months after the IOZM. The timing and associated SST patterns suggests that the IOZM is related to ENSO, and the basin- wide warming/cooling develops as a result of the decay of the IOZM spreading SST anomalies from western Indian Ocean to the eastern Indian Ocean. In contrast, in the Indian-Ocean Run, no oscillatory modes can be identified by the MSSA, even though the Indian Ocean SST variability is characterized by east-west SST contrast patterns similar to the IOZM. In both control and anomaly runs, IOZM-like SST variability appears to be associated with forcings from fluctuations of the Indian monsoon. Our modeling results suggest that the oscillatory feature of the IOZM is primarily forced by ENSO.

Yu, Jin-Yi↗

A General Theory of Unsteady Compressible Potential Aerodynamics

The general theory of potential aerodynamic flow around a lifting body having arbitrary shape and motion is presented. By using the Green function method, an integral representation for the potential is obtained for both supersonic and subsonic flow. 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, sigma, of the body. Hence, if the point P approaches the surface of the body, the representation reduces to an integro-differential equation relating the potential and its normal derivative (which is known from the boundary conditions) on the surface sigma. For the important practical case of small harmonic oscillation around a rest position, the equation reduces to a two-dimensional Fredholm integral equation of second-type. It is shown that this equation reduces properly to the lifting surface theories as well as other classical mathematical formulas. 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.

Morino, L.↗