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 1,567 records · Page 87

Boundary formulations for sensitivities of three-dimensional stress invariants

The direct, singular, boundary element analysis (BEA) formulation has been shown to provide a basis for a computationally efficient and accurate shape structural design sensitivity analysis (DSA) approach for three-dimensional solid objects. Within the boundary element analysis context, the theoretical formulation for sensitivities of important stress-related quantities including principal and deviatoric stresses, von Mises, maximum shear, and other stress invariants are presented, both for the surface as well as the interior of a continuum structure. Numerical results are given to demonstrate the accuracy of this approach.

Guru Prasad, K.↗

Extreme weather events

Models that try to forecast the detailed geographical distribution of oceanic and atmospheric variables affected by the ENSO cycle are briefly discussed. Combinations of singular-spectrum analysis and the maximum entropy method that hold promise for predicting the ENSO cycle 2-3 yrs in advance are addressed.

Keppenne, Christian L.↗

Robust dynamic inversion control laws for aircraft control

Dynamic inversion is a technique for control law design in which feedback is used to simultaneously cancel system dynamics and achieve desired dynamic response characteristics. However, dynamic inversion control laws lack robustness to modeling errors if improperly designed. This paper examines a simple linear example, control of roll rate about the body axis of high performance aircraft, to illustrate some robustness problems which may occur with a simple dynamic inversion control law. The paper demonstrates how structured singular value synthesis techniques can be used to enhance the robustness properties of the dynamic inversion controller.

Balas, Gary J.↗

Validity of the field line resonance expansion

The field line resonance expansion is studied with a magnetohydrodynamic (MHD) model applicable to earth's magnetotail. Using a Frobenius series solution, this expansion is found to be inconsistent. The inconsistency is seen by evaluating the terms in the model equations that the expansion neglects as small. Their effect is to provide a solution more singular than that determined by the expansion procedure. A tentative alternative expansion procedure, retaining these neglected terms approximately, yields a nonsingular solution, suggesting that resonant mode conversion rather than field line resonance is the appropriate phenomenon.

Hansen, P. J.↗

Viscous/Inviscid Interaction On An Oscillating Airfoil

Viscous/inviscid-interaction procedure developed for computation of two-dimensional flow of air about airfoil, angle of attack of which oscillates at high amplitude. Purpose of procedure to study flow effects in dynamic stall. Bulk of flow treated as inviscid, while effects of viscosity treated only on surface of airfoil and in boundary layer. Inviscid part of flow simulated via time-marching, unsteady-subsonic-flow panel method incorporating appropriate distributions of unsteady surface and wake singularities.

Carr, L. W.↗

Where's the Beaverhead beef?

Only rare quartz grains with single-set planar (1013) deformation features (PDF's) are present in breccia dikes found in association with uniformly oriented shatter cones that occur over an area 8 x 25 km. This suggests that the Beaverhead shocked rocks come from only the outer part of the central uplift of what must have been a large (greater than 100 km diameter) complex impact structure. An impact event of this magnitude on continental crust (thought to have occurred in late Precambrian or ealy Paleozoic time) could be expected to punctuate local geologic history. Furthermore, although it may now be covered, its scar should remain despite all the considerable subsequent erosion/deposition and tectonism since the impact. The following are three large-scale singularities or anomalies that may reflect the event and mark its source. (1) The Lemhi Arch is a major structural uplift that occurred in late Proterozoic-early Paleozoic time in East Central Idaho and caused the erosion of at least 4 km of sedimentary cover. This may be directly related to the impact. (2) Of the many thrust sheets comprising the Cordilleran belt, the Cabin plate that carries the shocked rocks is unique in that it alone intersected the crystalline basement. It also now marks the apex of the Southwest Montana Recess in the Sevier front. The basement uplift remaining from the impact may have constituted a mechanical obstacle to the advancing thrust sheets in Cretaceous time, causing the recess. (3) What could be interpreted as a roughly circular aeromagnetic anomaly approx. 70 km in diameter can be discerned in the state aeromagnetic map centered about 20 km southeast of Challis, Idaho, in the Lost River range. It is in approximately the right place, and ignoring the possibility that the anomalies have diverse causes and the circular pattern is coincidental, it may mark what remains of the buried central uplift structure.

Hargraves, R. B.↗

Calculation of stress intensity factors in an isotropic multicracked plate. Part 1: Theoretical development

An essential part of describing the damage state and predicting the damage growth in a multicracked plate is the accurate calculation of stress intensity factors (SIF's). Here, a methodology and rigorous solution formulation for SIF's of a multicracked plate, with fully interacting cracks, subjected to a far-field arbitrary stress state is presented. The fundamental perturbation problem is derived, and the steps needed to formulate the system of singular integral equations whose solution gives rise to the evaluation of the SIF's are identified. This analytical derivation and numerical solution are obtained by using intelligent application of symbolic computations and automatic FORTRAN generation capabilities (described in the second part of this paper). As a result, a symbolic/FORTRAN package, named SYMFRAC, that is capable of providing accurate SIF's at each crack tip was developed and validated.

Binienda, W. K.↗

Numerical study of the small scale structures in Boussinesq convection

Two-dimensional Boussinesq convection is studied numerically using two different methods: a filtered pseudospectral method and a high order accurate Essentially Nonoscillatory (ENO) scheme. The issue whether finite time singularity occurs for initially smooth flows is investigated. The numerical results suggest that the collapse of the bubble cap is unlikely to occur in resolved calculations. The strain rate corresponding to the intensification of the density gradient across the front saturates at the bubble cap. We also found that the cascade of energy to small scales is dominated by the formulation of thin and sharp fronts across which density jumps.

Weinan, E.↗

The program FANS-3D (finite analytic numerical simulation 3-dimensional) and its applications

In this study, the program named FANS-3D (Finite Analytic Numerical Simulation-3 Dimensional) is presented. FANS-3D was designed to solve problems of incompressible fluid flow and combined modes of heat transfer. It solves problems with conduction and convection modes of heat transfer in laminar flow, with provisions for radiation and turbulent flows. It can solve singular or conjugate modes of heat transfer. It also solves problems in natural convection, using the Boussinesq approximation. FANS-3D was designed to solve heat transfer problems inside one, two and three dimensional geometries that can be represented by orthogonal planes in a Cartesian coordinate system. It can solve internal and external flows using appropriate boundary conditions such as symmetric, periodic and user specified.

Bravo, Ramiro H.↗

A study of the generation of linear energy transfer spectra for space radiations

The conversion of particle-energy spectra into a linear energy transfer (LET) distribution is a guide in assessing biologically significant components. The mapping of LET to energy is triple valued and can be defined only on open subintervals. A well-defined numerical procedure is found to allow generation of LET spectra on the open subintervals that are integrable in spite of their singular nature.

Wilson, John W.↗

Reliable and Efficient Parallel Processing Algorithms and Architectures for Modern Signal Processing

Least-squares (LS) estimations and spectral decomposition algorithms constitute the heart of modern signal processing and communication problems. Implementations of recursive LS and spectral decomposition algorithms onto parallel processing architectures such as systolic arrays with efficient fault-tolerant schemes are the major concerns of this dissertation. There are four major results in this dissertation. First, we propose the systolic block Householder transformation with application to the recursive least-squares minimization. It is successfully implemented on a systolic array with a two-level pipelined implementation at the vector level as well as at the word level. Second, a real-time algorithm-based concurrent error detection scheme based on the residual method is proposed for the QRD RLS systolic array. The fault diagnosis, order degraded reconfiguration, and performance analysis are also considered. Third, the dynamic range, stability, error detection capability under finite-precision implementation, order degraded performance, and residual estimation under faulty situations for the QRD RLS systolic array are studied in details. Finally, we propose the use of multi-phase systolic algorithms for spectral decomposition based on the QR algorithm. Two systolic architectures, one based on triangular array and another based on rectangular array, are presented for the multiphase operations with fault-tolerant considerations. Eigenvectors and singular vectors can be easily obtained by using the multi-pase operations. Performance issues are also considered.

Liu, Kuojuey Ray↗

Inertial-space disturbance rejection for robotic manipulators

The disturbance rejection control problem for a 6-DOF (degree of freedom) PUMA manipulator mounted on a 3-DOF platform is investigated. A control algorithm is designed to track the desired position and attitude of the end-effector in inertial space, subject to unknown disturbances in the platform axes. Conditions for the stability of the closed-loop system are derived. The performance of the controller is compared for step, sinusoidal, and random disturbances in the platform rotational axis and in the neighborhood of kinematic singularities.

Holt, Kevin↗

Computational reacting gas dynamics

In the study of high speed flows at high altitudes, such as that encountered by re-entry spacecrafts, the interaction of chemical reactions and other non-equilibrium processes in the flow field with the gas dynamics is crucial. Generally speaking, problems of this level of complexity must resort to numerical methods for solutions, using sophisticated computational fluid dynamics (CFD) codes. The difficulties introduced by reacting gas dynamics can be classified into three distinct headings: (1) the usually inadequate knowledge of the reaction rate coefficients in the non-equilibrium reaction system; (2) the vastly larger number of unknowns involved in the computation and the expected stiffness of the equations; and (3) the interpretation of the detailed reacting CFD numerical results. The research performed accepts the premise that reacting flows of practical interest in the future will in general be too complex or 'untractable' for traditional analytical developments. The power of modern computers must be exploited. However, instead of focusing solely on the construction of numerical solutions of full-model equations, attention is also directed to the 'derivation' of the simplified model from the given full-model. In other words, the present research aims to utilize computations to do tasks which have traditionally been done by skilled theoreticians: to reduce an originally complex full-model system into an approximate but otherwise equivalent simplified model system. The tacit assumption is that once the appropriate simplified model is derived, the interpretation of the detailed numerical reacting CFD numerical results will become much easier. The approach of the research is called computational singular perturbation (CSP).

Lam, S. H.↗

Remote sensing of earth terrain

In remote sensing, the encountered geophysical media such as agricultural canopy, forest, snow, or ice are inhomogeneous and contain scatters in a random manner. Furthermore, weather conditions such as fog, mist, or snow cover can intervene the electromagnetic observation of the remotely sensed media. In the modelling of such media accounting for the weather effects, a multi-layer random medium model has been developed. The scattering effects of the random media are described by three-dimensional correlation functions with variances and correlation lengths corresponding to the fluctuation strengths and the physical geometry of the inhomogeneities, respectively. With proper consideration of the dyadic Green's function and its singularities, the strong fluctuation theory is used to calculate the effective permittivities which account for the modification of the wave speed and attenuation in the presence of the scatters. The distorted Born approximation is then applied to obtain the correlations of the scattered fields. From the correlation of the scattered field, calculated is the complete set of scattering coefficients for polarimetric radar observation or brightness temperature in passive radiometer applications. In the remote sensing of terrestrial ecosystems, the development of microwave remote sensing technology and the potential of SAR to measure vegetation structure and biomass have increased effort to conduct experimental and theoretical researches on the interactions between microwave and vegetation canopies. The overall objective is to develop inversion algorithms to retrieve biophysical parameters from radar data. In this perspective, theoretical models and experimental data are methodically interconnected in the following manner: Due to the complexity of the interactions involved, all theoretical models have limited domains of validity; the proposed solution is to use theoretical models, which is validated by experiments, to establish the region in which the radar response is most sensitive to the parameters of interest; theoretically simulated data will be used to generate simple invertible models over the region. For applications to the remote sensing of sea ice, the developed theoretical models need to be tested with experimental measurements. With measured ground truth such as ice thickness, temperature, salinity, and structure, input parameters to the theoretical models can be obtained to calculate the polarimetric scattering coefficients for radars or brightness temperature for radiometers and then compare theoretical results with experimental data. Validated models will play an important role in the interpretation and classification of ice in monitoring global ice cover from space borne remote sensors in the future. We present an inversion algorithm based on a recently developed inversion method referred to as the Renormalized Source-Type Integral Equation approach. The objective of this method is to overcome some of the limitations and difficulties of the iterative Born technique. It recasts the inversion, which is nonlinear in nature, in terms of the solution of a set of linear equations; however, the final inversion equation is still nonlinear. The derived inversion equation is an exact equation which sums up the iterative Neuman (or Born) series in a closed form and, thus, is a valid representation even in the case when the Born series diverges; hence, the name Renormalized Source-Type Integral Equation Approach.

Yueh, Herng-Aung↗

Hierarchical Poly Tree Configurations for the Solution of Dynamically Refined Finte Element Models

This paper demonstrates how a multilevel substructuring technique, called the Hierarchical Poly Tree (HPT), can be used to integrate a localized mesh refinement into the original finite element model more efficiently. The optimal HPT configurations for solving isoparametrically square h-, p-, and hp-extensions on single and multiprocessor computers is derived. In addition, the reduced number of stiffness matrix elements that must be stored when employing this type of solution strategy is quantified. Moreover, the HPT inherently provides localize 'error-trapping' and a logical, efficient means with which to isolate physically anomalous and analytically singular behavior.

Gute, G. D.↗

Optimal guidance law development for an advanced launch system

The proposed investigation on a Matched Asymptotic Expansion (MAE) method was carried out. It was concluded that the method of MAE is not applicable to launch vehicle ascent trajectory optimization due to a lack of a suitable stretched variable. More work was done on the earlier regular perturbation approach using a piecewise analytic zeroth order solution to generate a more accurate approximation. In the meantime, a singular perturbation approach using manifold theory is also under current investigation. Work on a general computational environment based on the use of MACSYMA and the weak Hamiltonian finite element method continued during this period. This methodology is capable of the solution of a large class of optimal control problems.

Calise, Anthony J.↗

Use of system identification techniques for improving airframe finite element models using test data

A method for using system identification techniques to improve airframe finite element models using test data was developed and demonstrated. The method uses linear sensitivity matrices to relate changes in selected physical parameters to changes in the total system matrices. The values for these physical parameters were determined using constrained optimization with singular value decomposition. The method was confirmed using both simple and complex finite element models for which pseudo-experimental data was synthesized directly from the finite element model. The method was then applied to a real airframe model which incorporated all of the complexities and details of a large finite element model and for which extensive test data was available. The method was shown to work, and the differences between the identified model and the measured results were considered satisfactory.

Hanagud, Sathya V.↗

A global path planning approach for redundant manipulators

A new approach for global path planning of redundant manipulators is proposed. It poses the path planning problem as a finite time nonlinear control problem. The solution is found by a Newton-Raphson type algorithm. This technique is capable of handling various goal task descriptions as well as incorporating both joint and task space constraints. The algorithm has shown promising preliminary results in planning joint path sequences for 3R and 4R planar robots to meet Cartesian tip tracking and goal endpoint planning. It is robust with respect to local path planning problems such as singularity considerations and local minimum problems. Repetitive joint path solutions for cyclic end-effector tasks are also generated. Eventual goals of this work include implementation on full spatial robots, as well as provision of an interface for supervisory input to aid in path planning for more complex problems.

Seereeram, Sanjeev↗