Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “differential geometry”

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 289 records · Page 16

Adhesive-Bonded Composite Joint Analysis with Delaminated Surface Ply Using Strain-Energy Release Rate

This paper presents an analytical model to determine the strain energy release rate due to an interlaminar crack of the surface ply in adhesively bonded composite joints subjected to axial tension. Single-lap shear-joint standard test specimen geometry with thick bondline is followed for model development. The field equations are formulated by using the first-order shear-deformation theory in laminated plates together with kinematics relations and force equilibrium conditions. The stress distributions for the adherends and adhesive are determined after the appropriate boundary and loading conditions are applied and the equations for the field displacements are solved. The system of second-order differential equations is solved to using the symbolic computation tool Maple 9.52 to provide displacements fields. The equivalent forces at the tip of the prescribed interlaminar crack are obtained based on interlaminar stress distributions. The strain energy release rate of the crack is then determined by using the crack closure method. Finite element analyses using the J integral as well as the crack closure method are performed to verify the developed analytical model. It has been shown that the results using the analytical method correlate well with the results from the finite element analyses. An attempt is made to predict the failure loads of the joints based on limited test data from the literature. The effectiveness of the inclusion of bondline thickness is justified when compared with the results obtained from the previous model in which a thin bondline and uniform adhesive stresses through the bondline thickness are assumed.

Chadegani, Alireza↗

Nature of convection-stabilized dc arcs in dual-flow nozzle geometry. I - The cold flow field and dc arc characteristics. II - Optical diagnostics and theory

Steady-state low-current air arcs in a dual-flow nozzle system are studied experimentally. The cold flow field with no arc is investigated using a 12.7-mm diameter dual-flow nozzle in a steady-flow facility. Mach number and mass flux distributions are determined for various nozzle-pressure ratios and nozzle-gap spacing. It is found that the shock waves in the converging-diverging nozzles result in a decrease in overal resistance by about 15 percent. Also, Schlieren and differential interferometry techniques are used to visualize the density gradients within the arc plasma and thermal mantle. Both optical techniques reveal a laminar arc structure for a reservoir pressure of 1 atm at various current levels. Experimentally determined axial static pressure and cold-flow mass flux rate distributions and a channel-flow model with constant arc temperatre are used to solve the energy integral for the arc radius as a function of axial distance. The arc electric field strength, voltage, resistance, and power are determined with Ohm's law and the total heat transfer is related to arc power.

Serbetci, Ilter↗

Current-State Constrained Filter Bank for Wald Testing of Spacecraft Conjunctions

We propose a filter bank consisting of an ordinary current-state extended Kalman filter, and two similar but constrained filters: one is constrained by a null hypothesis that the miss distance between two conjuncting spacecraft is inside their combined hard body radius at the predicted time of closest approach, and one is constrained by an alternative complementary hypothesis. The unconstrained filter is the basis of an initial screening for close approaches of interest. Once the initial screening detects a possibly risky conjunction, the unconstrained filter also governs measurement editing for all three filters, and predicts the time of closest approach. The constrained filters operate only when conjunctions of interest occur. The computed likelihoods of the innovations of the two constrained filters form a ratio for a Wald sequential probability ratio test. The Wald test guides risk mitigation maneuver decisions based on explicit false alarm and missed detection criteria. Since only current-state Kalman filtering is required to compute the innovations for the likelihood ratio, the present approach does not require the mapping of probability density forward to the time of closest approach. Instead, the hard-body constraint manifold is mapped to the filter update time by applying a sigma-point transformation to a projection function. Although many projectors are available, we choose one based on Lambert-style differential correction of the current-state velocity. We have tested our method using a scenario based on the Magnetospheric Multi-Scale mission, scheduled for launch in late 2014. This mission involves formation flight in highly elliptical orbits of four spinning spacecraft equipped with antennas extending 120 meters tip-to-tip. Eccentricities range from 0.82 to 0.91, and close approaches generally occur in the vicinity of perigee, where rapid changes in geometry may occur. Testing the method using two 12,000-case Monte Carlo simulations, we found the method achieved a missed detection rate of 0.1%, and a false alarm rate of 2%.

Carpenter, J. Russell↗

Analysis of Triangular Rollable and Collapsible Composite Booms Under the Effects of Gravity and Twist

Finite element analysis is used to investigate Triangular Rollable and Collapsible (TRAC) composite booms. The TRAC booms are subjected to a suite of loading conditions, including in-plane bending, out-of-plane bending, and axial compression. Simulations containing as-designed boom geometry are performed, as are simulations with geometric deviations associated with the as-built conditions of a sample TRAC boom. Namely, a cross-sectional twist along the length of the boom is incorporated into the analysis to study the effect of the twist on the deformation characteristics of the boom. Twisting in the sample boom was observed after short-term storage in a rolled configuration and may be attributed to differential creep deformation between the inner and outer flanges of the TRAC boom. Simulations are conducted with and without the effects of gravity to understand the significance of gravity on testing conditions. Booms of three lengths are evaluated: 3 meters, 7.5 meters and 30 meters. The mechanical responses of the TRAC booms are discussed with reference to the effects of gravity, boom length and cross-sectional twisting. Gravity is found to have a more significant effect on boom mechanical response as boom length is increased. Cross-sectional twisting is found to have a more significant effect on boom mechanical response in the shorter booms.

High-strain composites↗

Analysis of Triangular Rollable and Collapsible (TRAC) Composite Booms under the Effects of Gravity and Twist

Finite element analysis is used to investigate Triangular Rollable and Collapsible (TRAC) composite booms. The TRAC booms are subjected to a suite of loading conditions, including in-plane bending, out-of-plane bending, and axial compression. Simulations containing as-designed boom geometry are performed, as are simulations with geometric deviations associated with the as-built conditions of a sample TRAC boom. Namely, a cross-sectional twist along the length of the boom is incorporated into the analysis to study the effect of the twist on the deformation characteristics of the boom. Twisting in the sample boom was observed after short-term storage in a rolled configuration and may be attributed to differential creep deformation between the inner and outer flanges of the TRAC boom. Simulations are conducted with and without the effects of gravity to understand the significance of gravity on testing conditions. Booms of three lengths are evaluated: 3 meters, 7.5 meters and 30 meters. The mechanical responses of the TRAC booms are discussed with reference to the effects of gravity, boom length and cross-sectional twisting. Gravity is found to have a more significant effect on boom mechanical response as boom length is increased. Cross-sectional twisting is found to have a more significant effect on boom mechanical response in the shorter booms.

High-strain composites↗

Recognition of 3-D symmetric objects from range images in automated assembly tasks

A new technique is presented for the three dimensional recognition of symmetric objects from range images. Beginning from the implicit representation of quadrics, a set of ten coefficients is determined for symmetric objects like spheres, cones, cylinders, ellipsoids, and parallelepipeds. Instead of using these ten coefficients trying to fit them to smooth surfaces (patches) based on the traditional way of determining curvatures, a new approach based on two dimensional geometry is used. For each symmetric object, a unique set of two dimensional curves is obtained from the various angles at which the object is intersected with a plane. Using the same ten coefficients obtained earlier and based on the discriminant method, each of these curves is classified as a parabola, circle, ellipse, or hyperbola. Each symmetric object is found to possess a unique set of these two dimensional curves whereby it can be differentiated from the others. It is shown that instead of using the three dimensional discriminant which involves evaluation of the rank of its matrix, it is sufficient to use the two dimensional discriminant which only requires three arithmetic operations.

Alvertos, Nicolas↗

Three-dimensional aspects of interplanetary shock waves

The three-dimensional geometries of two classes of interplanetary shock waves are analyzed: (1) the more or less spherical shock wave originating from some short-lived solar event, such as a solar flare, and (2) the corotating shock pair formed by the interaction of long-lived solar wind streams. The systematic distortion of a spherical wave due to persistent latitudinal solar wind structure should amount to 25 deg for an equator-to-pole solar wind speed differential of 200 km/sec. A geometrical argument is given for determining the heliocentric distance of the onset of corotating shocks formed by contiguous long-lived solar streams. Estimates are obtained for the heliocentric distance to the inner edges of shock waves as a function of speed difference between the streams and as a function of heliographic latitude. Means by which the three-dimensional structures studied can be observed by out-of-the-ecliptic space missions are discussed.

Siscoe, G. L.↗

Energy transfer in volume-reflecting heat shields

An approximate analysis of radiative transfer in highly scattering materials was developed based on the Kubelka-Munk differential equations--a set of two differential equations representing the spatial rate of change of radiative half-fluxes within the scattering media. These approximate solutions of the Kubelka-Munk equations together with analytic solutions for the steady state temperature distribution for two types of boundary conditions are given. These solutions show the influence of back surface reflectance, scattering power, incident radiative flux parameter and boundary conductive flux parameter on overall reflectance and temperature distributions. This radiation field analysis, adapted to spherical geometry, was applied to the evaluation of the thermal performance of teflon and fritted quartz as heat protection materials for entry into the atmosphere of Jupiter.

Weston, K. C.↗

Dynamics of early planetary gear trains

A method to analyze the static and dynamic loads in a planetary gear train was developed. A variable-variable mesh stiffness (VVMS) model was used to simulate the external and internal spur gear mesh behavior, and an equivalent conventional gear train concept was adapted for the dynamic studies. The analysis can be applied either involute or noninvolute spur gearing. By utilizing the equivalent gear train concept, the developed method may be extended for use for all types of epicyclic gearing. The method is incorporated into a computer program so that the static and dynamic behavior of individual components can be examined. Items considered in the analysis are: (1) static and dynamic load sharing among the planets; (2) floating or fixed Sun gear; (3) actual tooth geometry, including errors and modifications; (4) positioning errors of the planet gears; (5) torque variations due to noninvolute gear action. A mathematical model comprised of power source, load, and planetary transmission is used to determine the instantaneous loads to which the components are subjected. It considers fluctuating output torque, elastic behavior in the system, and loss of contact between gear teeth. The dynamic model has nine degrees of freedom resulting in a set of simultaneous second order differential equations with time varying coefficients, which are solved numerically. The computer program was used to determine the effect of manufacturing errors, damping and component stiffness, and transmitted load on dynamic behavior. It is indicated that this methodology offers the designer/analyst a comprehensive tool with which planetary drives may be quickly and effectively evaluated.

August, R.↗

Initial On-Orbit Spatial Resolution Characterization of OrbView-3 Panchromatic Images

Characterization was conducted under the Memorandum of Understanding among Orbital Sciences Corp., ORBIMAGE, Inc., and NASA Applied Sciences Directorate. Acquired five OrbView-3 panchromatic images of the permanent Stennis Space Center edge targets painted on a concrete surface. Each image is available at two processing levels: Georaw and Basic. Georaw is an intermediate image in which individual pixels are aligned by a nominal shift in the along-scan direction to adjust for the staggered layout of the panchromatic detectors along the focal plane array. Georaw images are engineering data and are not delivered to customers. The Basic product includes a cubic interpolation to align the pixels better along the focal plane and to correct for sensor artifacts, such as smile and attitude smoothing. This product retains satellite geometry - no rectification is performed. Processing of the characterized images did not include image sharpening, which is applied by default to OrbView-3 image products delivered by ORBIMAGE to customers. Edge responses were extracted from images of tilted edges in two directions: along-scan and cross-scan. Each edge response was approximated with a superposition of three sigmoidal functions through a nonlinear least-squares curve-fitting. Line Spread Functions (LSF) were derived by differentiation of the analytical approximation. Modulation Transfer Functions (MTF) were obtained after applying the discrete Fourier transform to the LSF.

Blonski, Slawomir↗

Streamlined Convergence Acceleration for CFD Codes

Enigma, a simplified interface to the PETSc library, is shown to enable the rapid solution of discrete partial differential equations. Two CFD codes, LAURA and HyperSolve, use Enigma to compute steady solutions of the Navier-Stokes equations. Using PETSc, Enigma is shown to provide a Jacobian-Free Newton-Krylov method (JFNK), globalized with pseudotransient continuation, that improves efficiency over the point-implicit relaxation method traditionally used by LAURA. It is shown that iterative error has a large impact on surface heat transfer predicted by LAURA on an axisymmetric sphere-cone geometry. Also, the convergence rate of HyperSolve simulating subsonic flow over a delta wing geometry with the JFNK method is shown to be more efficient than employing a defect correction method as the nonlinear solver.

Thompson, Kyle B.↗

Stresses in a three-dimensional unidirectional composite containing broken fibers

An approximate solution is developed for the determination of the interlaminar normal and shear stresses in the vicinity of a crack in a three dimensional composite containing unidirectional linearly elastic fibers in an infinite linearly elastic matrix. In order to reduce the complexity of the formulation, certain assumptions are made as to the physically significant stresses to be retained. These simplifications reduce the partial differential equations of elasticity to differential-difference equations which are tractable using Fourier transform techniques. This 'material modeling' approach is in contrast with solutions developed by considering each lamina as a homogeneous, orthotropic layer. The resulting solution does not contain the classical singular stress field for the fibers and the influence of broken fibers on unbroken fibers is felt by a change in stress concentration factors. The matrix stresses however, are unbounded as the fiber spacing vanishes and an equivalent fiber-matrix geometry is proposed which gives the correct singular behavior. The numerical solution is considered in detail and several specific examples are presented. The potential for damaged or debonded zones to be generated by an embedded crack is discussed, and stress concentration factors for fibers near the crack are given. Detailed comparisons are made between the present solution, the analogous two-dimensional problem, and corresponding shear-lag models.

Goree, J. G.↗

Preferred orbit planes in the gravitational field of a tumbling spheroidal galaxy

An analytic method is described by means of which the preferred planes into which gas will settle by dissipative differential precession in a slightly nonspherical, nonstatic potential can be determined. Attention is given to the analysis of circular orbit precession in the external potential of a slightly prolate spheroidal mass that is tumbling about a short axis, providing a crude model for barred spiral galaxies as well as for prolate elliptical ones. The results obtained agree with previous work on the existence, orientation, and stability of periodic orbits in triaxial potentials. Depending on the source and injection geometry of the gas and on the nature of dissipation in the gas, the orbital dynamics presented could lead to enhanced, or even catastrophic, radial inflow and to steady state warped structures.

Durisen, R. H.↗

Ray trajectories in a torus: An application of MACSYMA to a complex numerical computation

The study of ray trajectories of plasma waves in a torodial geometry using MACSYMA is an example of how symbolic, numerical, and graphical facilities can be used in concert to accomplish a complex computational goal. Computational features of this study which are of particular significance include: the derivation of code (i.e. writing functions to generate program fragments), the use of array functions to simplify the specification of a numerical iteration scheme, and the graphical presentation of the results. Mathematically, this study originates in the solution of a linear inhomogeneous partial differential equation in 3 dimensions by the method of characteristics. While it is possible to describe this equation compactly by using vector notation, and by specifying the spatial variation of the coefficients in terms of intermediate parameters, the transformation of the equation into a form amenable to solution is very tedious. A MACSYMA program is presented for obtaining description of the rf field structure excited by a waveguide located at the edge of a toroidal plasma confinement device.

Kulp, J. L.↗

Tracking Landsat-5 by a differential GPS technique

As part of an international campaign to develop precise geodetic applications of the Global Positioning System (GPS), the Jet Propulsion Laboratory is conducting a demonstration of differential GPS tracking using Landsat-5. Two strategies have been investigated: one in which only the Landsat-5 orbit is estimated, and one in which both the Landsat-5 and GPS orbits are estimated together. Error studies show that under the limited conditions of the experiment, three-dimensional Landsat-5 position accuracies of about 5 m with the first strategy and 2 m with the second strategy can be achieved over a 20-min period of good observing geometry. Orbit determination results using a version of the first strategy appear to achieve the 5 m goal. This is supported by various formal error measures and independent comparisons. The more powerful strategy has yet to be carried out.

Yunck, T. P.↗

Study of a pursuit-evasion guidance law for high performance aircraft

The study of a one-on-one aircraft pursuit-evasion guidance scheme for high-performance aircraft is discussed. The research objective is to implement a guidance law derived earlier using differential game theory in conjunction with the theory of feedback linearization. Unlike earlier research in this area, the present formulation explicitly recognizes the two-sided nature of the pursuit-evasion scenario. The present research implements the guidance law in a realistic model of a modern high-performance fighter aircraft. Also discussed are the details of the guidance law, implementation in a highly detailed simulation of a high-performance fighter, and numerical results for two engagement geometries. Modifications of the guidance law for onboard implementation is also discussed.

Williams, Peggy S.↗

Helicopter approach capability using the differential global positioning system

The results of flight tests to determine the feasibility of using the Global Positioning System (GPS) in the differential mode (DGPS) to provide high accuracy, precision navigation, and guidance for helicopter approaches to landing are presented. The airborne DGPS receiver and associated equipment is installed in a NASA UH-60 Black Hawk helicopter. The ground-based DGPS reference receiver is located at a surveyed test site and is equipped with a real-time VHF data link to transmit correction information to the airborne DGPS receiver. The corrected airborne DGPS information, together with the preset approach geometry, is used to calculate guidance commands which are sent to the aircraft's approach guidance instruments. The use of DGPS derived guidance for helicopter approaches to landing is evaluated by comparing the DGPS data with the laser tracker truth data. Both standard (3 deg) and steep (6 deg and 9 deg) glideslope straight-in approaches were flown. DGPS positioning accuracy based on a time history analysis of the entire approach was 0.2 m (mean) +/- 1.8 m (2 sigma) laterally and -2.0 m (mean) +/- 3.5 m (2 sigma) vertically for 3 deg glideslope approaches, -0.1 m (mean) +/- 1.5 m (2 sigma) laterally and -1.1 m (mean) +/- 3.5 m (2 sigma) vertically for 6 deg glideslope approaches and 0.2 m (mean) +/- 1.3 m (2 sigma) laterally and -1.0 m (mean) +/- 2.8 m (2 sigma) vertically for 9 deg glideslope approaches. DGPS positioning accuracy at the 200 ft decision height (DH) on a standard 3 deg slideslope approach was 0.3 m (mean) +/- 1.5 m (2 sigma) laterally and -2.3 m (mean) +/- 1.6 m (2 sigma) vertically. These errors indicate that the helicopter position based on DGPS guidance satisfies the International Civil Aviation Organization (ICAO) Category 1 (CAT 1) lateral and vertical navigational accuracy requirements.

Kaufmann, David N.↗

Helicopter approach capability using the differential Global Positioning System

The results of flight tests to determine the feasibility of using the Global Positioning System (GPS) in the differential mode (DGPS) to provide high accuracy, precision navigation and guidance for helicopter approaches to landing are presented. The airborne DGPS receiver and associated equipment is installed in a NASA UH-60 Black Hawk helicopter. The ground-based DGPS reference receiver is located at a surveyed test site and is equipped with a real-time VHF data link to transmit correction information to the airborne DGPS receiver. The corrected airborne DGPS information, together with the preset approach geometry, is used to calculate guidance commands which are sent to the aircraft's approach guidance instruments. The use of DGPS derived guidance for helicopter approaches to landing is evaluated by comparing the DGPS data with the laser tracker truth data. Both standard (3 degrees) and steep (6 degrees and 9 degrees) glidescope straight-in approaches were flown. DGPS positioning accuracy based on a time history analysis of the entire approach was 0.2 m (mean) +/- 1.8 m (2 sigma) laterally and -2.0 m (mean) +/- 3.5 m (2 sigma) vertically for 3 degree glidescope approaches, -0.1 m (mean) +/- 1.5 m (2 sigma) laterally and -1.1 m (mean) +/- 3.5 m (2 sigma) vertically for 6 degree glidescope approaches, and 0.2 m (mean) +/- 1.3 m (2 sigma) laterally and -1.0 m (mean) +/- 2.8 (2 sigma) vertically for 9 degree glidescope approaches. DGPS positioning accuracy at the 200 ft decision height on a standard 3 degree glidescope approach was 0.3 m (mean) +/- 1.5 m (2 sigma) laterally and -2.3 m (mean) +/- 1.6 m (2 sigma) vertically. These errors indicate that the helicopter position based on DGPS guidance satisfies the International Civil Aviation Organization Category 1 lateral and vertical accuracy requirements.

Kaufmann, David N.↗