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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Solving the Linear Balance Equation on the Globe as a Generalized Inverse Problem

A generalized (pseudo) inverse technique was developed to facilitate a better understanding of the numerical effects of tropical singularities inherent in the spectral linear balance equation (LBE). Depending upon the truncation, various levels of determinancy are manifest. The traditional fully-determined (FD) systems give rise to a strong response, while the under-determined (UD) systems yield a weak response to the tropical singularities. The over-determined (OD) systems result in a modest response and a large residual in the tropics. The FD and OD systems can be alternatively solved by the iterative method. Differences in the solutions of an UD system exist between the inverse technique and the iterative method owing to the non- uniqueness of the problem. A realistic balanced wind was obtained by solving the principal components of the spectral LBE in terms of vorticity in an intermediate resolution. Improved solutions were achieved by including the singular-component solutions which best fit the observed wind data.

Lu, Huei-Iin↗
NASA NTRS2000

Cooperative Robotics and the Search for Extraterrestrial Life

If we think tenuous abodes of life may be hiding in remote extraterrestrial environmental niches, and if we want to assess the biological status of a given locale or entire planet before sending humans (perhaps because of contamination concerns or other motivations) then we face the challenge of robotically exploring a large space efficiently and in enough detail to have confidence in our assessment of the biological status of the environment in question. On our present schedule of perhaps two or so missions per opportunity, we will likely need a different exploratory approach than singular stationary landers or singular rover missions or sample return, because there appear to be fundamental limitations in these mission profiles to-obtain the many samples we will likely need if we want to have confidence in assessing the biological status of an environment in which life could be hiding in remote environmental niches. Singular rover missions can potentially accommodate sampling over a fairly large area, but are still limited by range and can be a single point of failure. More importantly, such mission profiles have limited payload capabilities which are unlikely to meet the demanding requirements of life-detection. Sample return has the advantage of allowing sophisticated analysis of the sample, but also has the severe limitations associated with only being able to bring back a few samples. This presentation will suggest two cooperative robotic approaches for exploration that have the potential to overcome these difficulties and facilitate efficient and thorough life-detecting exploration of a large space. Given the two premises state above, it appears at least two fundamental challenges have to be met simultaneously: coverage of a large space and bringing to bear a sophisticated suite of detection and experimental payloads on any specific location in order to address a major challenge in looking for extraterrestrial life: namely, executing a wide variety of detection scenarios and in situ experiments in order to gather the required data for a confident assessment that life has been detected and to, more generally, cover a wide range of extraterrestrial life possibilities. Cooperative robotics ]ends itself to this kind of problem because cooperation among the combined capabilities of a variety of simple single function agents can give rise to fairly complex task execution such as the search for and detection of extraterrestrial life.

Lupisella, Mark L.↗
NASA NTRS2000

Higher Order Modeling In the BEM/FEM Hybrid Formulation

Hybrid formulations using low order curl-conforming bases to represent the total electric field within a finite element region and low order divergence-conforming bases to represent equivalent electric and magnetic currents on the boundary are well known. However, higher-order divergence and curl-conforming bases have been shown to provide significant benefits in convergence rates and accuracy when employed in strictly integral equation and strictly finite element formulations. In this paper, a hybrid electric field formulation employing higher order bases is presented. The paper addresses benefits and issues associated with using higher order divergence-and curl-conforming bases in the hybrid finite element/boundary element electric field formulation. The method of singularity subtraction may be used to compute the self terms of the boundary integral when the bases are of low order. But this method becomes laborious and requires great care when the divergence conforming bases are of higher order. In order to handle these singularities simply and accurately, a generalized Gaussian quadrature method is employed in which the expansion functions account for the singularity. In preliminary tests of the higher order hybrid formulation, the equivalent electric current induced by scattering of a plane wave from a square dielectric cylinder is examined. Accurate results are obtained using only a two-triangle mesh when the current basis is of order 4 or 5. Additional results are presented comparing the error obtained using higher order bases to that obtained using lower order bases when the number of unknowns is approximately equal. Also, convergence rates obtained with higher order bases are compared to those obtained with lower order bases for selected sample problems.

Fink, Patrick W.↗
NASA NTRS2000

Cooperative Robotics and the Search for Extraterrestrial Life

If we think tenuous abodes of life may be hiding in remote extraterrestrial environmental niches, and if we want to assess the biological status of a given locale or entire planet before sending humans (perhaps because of contamination concerns or other motivations) then we face the challenge of robotically exploring a large space efficiently and in enough detail to have confidence in our assessment of the biological status of the environment in question. On our present schedule of perhaps two or so missions per opportunity, we will likely need a different exploratory approach than singular stationary landers or singular rover missions or sample return, because there appear to be fundamental limitations in these mission profiles to obtain the many samples we will likely need if we want to have confidence in assessing the biological status of an environment in which life could be hiding in remote environmental niches. Singular rover missions can potentially accommodate sampling over a fairly large area, but are still limited by range and can be a single point of failure. More importantly, such mission profiles have limited payload capabilities which are unlikely to meet the demanding requirements of life-detection. Sample return has the advantage of allowing sophisticated analysis of the sample, but also has the severe limitations associated with only being able to bring back a few samples. This presentation will suggest two cooperative robotic approaches for exploration that have the potential to overcome these difficulties and facilitate efficient and thorough life-detecting exploration of a large space. Given the two premises stated above, it appears at least two fundamental challenges have to be met simultaneously: (1) coverage of a large space and (2) bringing to bear a sophisticated suite of detection and experimental payloads on any specific location in order to address a major challenge in looking for extraterrestrial life: namely, executing a wide variety of detection scenarios and in situ experiments in order to gather the required data for a confident assessment that life has been detected and to, more generally, cover a wide range of extraterrestrial life possibilities. Cooperative robotics lends itself to this kind of problem because cooperation among the combined capabilities of a variety of simple single function agents can give rise to fairly complex task execution such as the search for and detection of extraterrestrial life.

Lupisella, M. L.↗
NASA NTRS2003

Attitude Estimation or Quaternion Estimation?

The attitude of spacecraft is represented by a 3x3 orthogonal matrix with unity determinant, which belongs to the three-dimensional special orthogonal group SO(3). The fact that all three-parameter representations of SO(3) are singular or discontinuous for certain attitudes has led to the use of higher-dimensional nonsingular parameterizations, especially the four-component quaternion. In attitude estimation, we are faced with the alternatives of using an attitude representation that is either singular or redundant. Estimation procedures fall into three broad classes. The first estimates a three-dimensional representation of attitude deviations from a reference attitude parameterized by a higher-dimensional nonsingular parameterization. The deviations from the reference are assumed to be small enough to avoid any singularity or discontinuity of the three-dimensional parameterization. The second class, which estimates a higher-dimensional representation subject to enough constraints to leave only three degrees of freedom, is difficult to formulate and apply consistently. The third class estimates a representation of SO(3) with more than three dimensions, treating the parameters as independent. We refer to the most common member of this class as quaternion estimation, to contrast it with attitude estimation. We analyze the first and third of these approaches in the context of an extended Kalman filter with simplified kinematics and measurement models.

F Landis Markley↗
NASA NTRS2003

Flutter and Divergence Analysis using the Generalized Aeroelastic Analysis Method

The Generalized Aeroelastic Analysis Method (GAAM) is applied to the analysis of three well-studied checkcases: restrained and unrestrained airfoil models, and a wing model. An eigenvalue iteration procedure is used for converging upon roots of the complex stability matrix. For the airfoil models, exact root loci are given which clearly illustrate the nature of the flutter and divergence instabilities. The singularities involved are enumerated, including an additional pole at the origin for the unrestrained airfoil case and the emergence of an additional pole on the positive real axis at the divergence speed for the restrained airfoil case. Inconsistencies and differences among published aeroelastic root loci and the new, exact results are discussed and resolved. The generalization of a Doublet Lattice Method computer code is described and the code is applied to the calculation of root loci for the wing model for incompressible and for subsonic flow conditions. The error introduced in the reduction of the singular integral equation underlying the unsteady lifting surface theory to a linear algebraic equation is discussed. Acknowledging this inherent error, the solutions of the algebraic equation by GAAM are termed 'exact.' The singularities of the problem are discussed and exponential series approximations used in the evaluation of the kernel function shown to introduce a dense collection of poles and zeroes on the negative real axis. Again, inconsistencies and differences among published aeroelastic root loci and the new 'exact' results are discussed and resolved. In all cases, aeroelastic flutter and divergence speeds and frequencies are in good agreement with published results. The GAAM solution procedure allows complete control over Mach number, velocity, density, and complex frequency. Thus all points on the computed root loci can be matched-point, consistent solutions without recourse to complex mode tracking logic or dataset interpolation, as in the k and p-k solution methods.

Edwards, John W.↗
NASA NTRS2005

Micromechanics, Fracture Mechanics and Gas Permeability of Composite Laminates for Cryogenic Storage Systems

A micromechanics method is developed to investigate microcrack propagation in a liquid hydrogen composite tank at cryogenic temperature. The unit cell is modeled using square and hexagonal shapes depends on fiber and matrix layout from microscopic images of composite laminates. Periodic boundary conditions are applied to the unit cell. The temperature dependent properties are taken into account in the analysis. The laminate properties estimated by the micromechanics method are compared with empirical solutions using constituent properties. The micro stresses in the fiber and matrix phases based on boundary conditions in laminate level are calculated to predict the formation of microcracks in the matrix. The method is applied to an actual liquid hydrogen storage system. The analysis predicts micro stresses in the matrix phase are large enough to cause microcracks in the composite. Stress singularity of a transverse crack normal to a ply-interface is investigated to predict the fracture behavior at cryogenic conditions using analytical and finite element analysis. When a transverse crack touches a ply-interface of a composite layer with same fiber orientation, the stress singularity is equal to 1/2. When the transverse crack propagates to a stiffer layer normal to the ply-direction, the singularity becomes less than 1/2 and vice versa. Finite element analysis is performed to predict the fracture toughness of a laminated beam subjected to fracture loads measured by four-point bending tests at room and cryogenic temperatures. As results, the fracture load at cryogenic temperature is significantly lower than that at room temperature. However, when thermal stresses are taken into consideration, for both cases of room and cryogenic temperatures, the difference of the fracture toughness becomes insignificant. The result indicates fracture toughness is a characteristic property, which is independent to temperature changes. The experimental analysis is performed to investigate the effect of cryogenic cycling on permeability for various composite material systems. Textile composites have lower permeability than laminated composites even with increasing number of cryogenic cycle. Nano-particles dispersed in laminated composites do not show improvement on permeability. The optical inspection is performed to investigate the microcrack propagation and void content in laminated composites and compared the microscopic results before and after cryogenic cycling.

Choi, Sukjoo↗
NASA NTRS2006

Refinement of Methods for Evaluation of Near-Hypersingular Integrals in BEM Formulations

In this paper, we present advances in singularity cancellation techniques applied to integrals in BEM formulations that are nearly hypersingular. Significant advances have been made recently in singularity cancellation techniques applied to 1 R type kernels [M. Khayat, D. Wilton, IEEE Trans. Antennas and Prop., 53, pp. 3180-3190, 2005], as well as to the gradients of these kernels [P. Fink, D. Wilton, and M. Khayat, Proc. ICEAA, pp. 861-864, Torino, Italy, 2005] on curved subdomains. In these approaches, the source triangle is divided into three tangent subtriangles with a common vertex at the normal projection of the observation point onto the source element or the extended surface containing it. The geometry of a typical tangent subtriangle and its local rectangular coordinate system with origin at the projected observation point is shown in Fig. 1. Whereas singularity cancellation techniques for 1 R type kernels are now nearing maturity, the efficient handling of near-hypersingular kernels still needs attention. For example, in the gradient reference above, techniques are presented for computing the normal component of the gradient relative to the plane containing the tangent subtriangle. These techniques, summarized in the transformations in Table 1, are applied at the sub-triangle level and correspond particularly to the case in which the normal projection of the observation point lies within the boundary of the source element. They are found to be highly efficient as z approaches zero. Here, we extend the approach to cover two instances not previously addressed. First, we consider the case in which the normal projection of the observation point lies external to the source element. For such cases, we find that simple modifications to the transformations of Table 1 permit significant savings in computational cost. Second, we present techniques that permit accurate computation of the tangential components of the gradient; i.e., tangent to the plane containing the source element.

Fink, Patricia W.↗
NASA NTRS2006

Fibonacci Grids

Recent years have seen a resurgence of interest in a variety of non-standard computational grids for global numerical prediction. The motivation has been to reduce problems associated with the converging meridians and the polar singularities of conventional regular latitude-longitude grids. A further impetus has come from the adoption of massively parallel computers, for which it is necessary to distribute work equitably across the processors; this is more practicable for some non-standard grids. Desirable attributes of a grid for high-order spatial finite differencing are: (i) geometrical regularity; (ii) a homogeneous and approximately isotropic spatial resolution; (iii) a low proportion of the grid points where the numerical procedures require special customization (such as near coordinate singularities or grid edges). One family of grid arrangements which, to our knowledge, has never before been applied to numerical weather prediction, but which appears to offer several technical advantages, are what we shall refer to as "Fibonacci grids". They can be thought of as mathematically ideal generalizations of the patterns occurring naturally in the spiral arrangements of seeds and fruit found in sunflower heads and pineapples (to give two of the many botanical examples). These grids possess virtually uniform and highly isotropic resolution, with an equal area for each grid point. There are only two compact singular regions on a sphere that require customized numerics. We demonstrate the practicality of these grids in shallow water simulations, and discuss the prospects for efficiently using these frameworks in three-dimensional semi-implicit and semi-Lagrangian weather prediction or climate models.

Swinbank, Richard↗
NASA NTRS2008

Integrating the Gradient of the Thin Wire Kernel

A formulation for integrating the gradient of the thin wire kernel is presented. This approach employs a new expression for the gradient of the thin wire kernel derived from a recent technique for numerically evaluating the exact thin wire kernel. This approach should provide essentially arbitrary accuracy and may be used with higher-order elements and basis functions using the procedure described in [4].When the source and observation points are close, the potential integrals over wire segments involving the wire kernel are split into parts to handle the singular behavior of the integrand [1]. The singularity characteristics of the gradient of the wire kernel are different than those of the wire kernel, and the axial and radial components have different singularities. The characteristics of the gradient of the wire kernel are discussed in [2]. To evaluate the near electric and magnetic fields of a wire, the integration of the gradient of the wire kernel needs to be calculated over the source wire. Since the vector bases for current have constant direction on linear wire segments, these integrals reduce to integrals of the form

Champagne, Nathan J.↗
NASA NTRS2009

Magnetopause Reconnection Impact Parameters from Multiple Spacecraft Magnetic Field Measurements

We present a novel technique that exploits multiple spacecraft data to determine the impact parameters of the most general form of magnetic reconnection at the magnetopause. The method consists of a superposed epoch of multiple spacecraft magnetometer measurements that yields the instantaneous magnetic spatial gradients near a magnetopause reconnection site. The gradients establish the instantaneous positions of the spacecraft relative to the reconnection site. The analysis is well suited to evaluating the spatial scales of singular field line reconnection, which is characterized by a two-dimensional x-type topology adjacent and perpendicular to a reconnecting singular field line. Application of the method to Cluster data known to lie in the vicinity of a northward IMF reconnection site establishes a field topology consistent with singular field line reconnection and a normal magnetic field component of 20 nT. The corresponding current structure consists of a 130 km sheet possibly embedding a thinner. bifurcated sheet.

Wendel, Deirdre E.↗
NASA NTRS1972

Lifting-surface theory for calculating the loading induced on a wing by a flap

A method is described for using lifting-surface theory to obtain the pressure distribution on a wing with a trailing-edge flap or control surface. The loading has a logarithmic singularity at the flap edges, which may be determined directly by the method of matched asymptotic expansions. Expressions are given for the singular flap loading for various flap hinge line and side edge geometries, both for steady and unsteady flap deflection. The regular part of the flap loading must be obtained by inverting the lifting-surface-theory integral equation relating the pressure and the downwash on the wing: procedures are described to accomplish this for a general wing and flap geometry. The method is applied to several example wings, and the results are compared with experimental data. Theory and test correlate well.

Johnson, W. A.↗
NASA NTRS1972

Prediction of unsteady aerodynamic loadings caused by trailing edge control surface motions in subsonic compressible flow: Analysis and results

A theoretical analysis and a computer program have been developed for the prediction of unsteady lifting surface loadings caused by motions of trailing edge control surfaces having sealed gaps. The final form of the downwash integral equation has been formulated by isolating the singularities from the non-singular terms and establishing a preferred solution process to remove and evaluate the downwash discontinuities in a systematic manner. Comparisons of theoretical and experimental pressure data are made for several control surface configurations. The comparisons indicate that reasonably accurate theoretical pressure distributions and generalized forces may be obtained for a wide variety of control surface configurations. Spanwise symmetry or antisymmetry of motion, and up to four control surfaces on each half span can be accommodated.

Rowe, W. S.↗
NASA NTRS1972

Mode coupling and wave particle interactions for unstable ion acoustic waves

A theory for the spatial development of linearly unstable, coupled waves is presented in which both quasi-linear and mode coupling effects are treated in a self-consistent manner. Steady state excitation of two waves is assumed at the boundary x = 0, the plasma being homogeneous in the y and z directions. Coupled equations are derived for the x dependence of the amplitudes of the primary waves and the secondary waves, correct through second order terms in the wave amplitude, but without usual approximation of small growth rates. This general formalism is then applied to the case of coupled ion acoustic waves driven unstable by an ion beam streaming in the direction of the x axis. If the modifications of the ion beam by the waves (quasi-linear effects) are ignored, explosive instabilities (singularities in all of the amplitudes at finite x) are found, even when all of the waves have positive energy. If these wave-particle interactions are included, the solutions are no longer singular, and all of the amplitudes have finite maxima.

Martin, P.↗
NASA NTRS1972

Variational solutions by the use of stepwise constant functions. I - Linear case.

Demonstration that accurate solutions to the integral equation or to the extremal function can be easily obtained by the variational method in many cases when a stepwise constant function is used for the trial function. The evaluation procedure is simple and straightforward; integrals of the kernel function can be evaluated analytically; the method provides, when the solution involves singularities, the best mean value across the singularity; the results are accurate in both the detailed physical quantities and their averages; the resultant solution can be further integrated analytically over the parameters associated with the problem; and the method can be readily applied to nonlinear integral equations.

Yoshikawa, K. K.↗
NASA NTRS1972

A circumferential crack in a cylindrical shell under tension.

A closed cylindrical shell under uniform internal pressure has a slit around a portion of its circumference. Linear shallow shell theory predicts inverse square-root-type singularities in certain of the stresses at the crack tips. This paper reports the computed strength of these singularities for different values of a dimensionless parameter based on crack length, shell radius and shell thickness.

Duncan-Fama, M. E.↗
NASA NTRS1972

A method for integrating the boundary-layer equations through a region of reverse flow.

Development of a procedure for numerically integrating the boundary-layer equations through a region of reverse flow which takes downstream influence into account. This method is applied to the problem of uniform flow past a parallel flat plate of finite length whose surface has a constant velocity directed opposite to that of the main stream. Although singularities occur at both the point of detachment (x sub s) and reattachment (x sub r) of the psi = 0 streamline, this integration technique provides a solution which ceases to apply only in the close proximity of these singular points. From this solution it is evident that, throughout a large portion of the separated region, the flow is strongly affected by conditions near x sub r, thereby demonstrating the importance of allowing information to be transmitted upstream in a region of backflow. Near x sub s, however, it is found that, in spite of the presence of reverse flow, the solution has a self-similar form in this particular example.

Klemp, J. B.↗
NASA NTRS1972

General theory of spherically symmetric boundary-value problems of the linear transport theory.

A general theory of spherically symmetric boundary-value problems of the one-speed neutron transport theory is presented. The formulation is also applicable to the 'gray' problems of radiative transfer. The Green's function for the purely absorbing medium is utilized in obtaining the normal mode expansion of the angular densities for both interior and exterior problems. As the integral equations for unknown coefficients are regular, a general class of reduction operators is introduced to reduce such regular integral equations to singular ones with a Cauchy-type kernel. Such operators then permit one to solve the singular integral equations by the standard techniques due to Muskhelishvili. We discuss several spherically symmetric problems. However, the treatment is kept sufficiently general to deal with problems lacking azimuthal symmetry. In particular the procedure seems to work for regions whose boundary coincides with one of the coordinate surfaces for which the Helmholtz equation is separable.

Kanal, M.↗