Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inversion methods”

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 775 records · Page 43

Stress analysis method for clearance-fit joints with bearing-bypass loads

Within a multi-fastener joint, fastener holes may be subjected to the combined effects of bearing loads and loads that bypass the hole to be reacted elsewhere in the joint. The analysis of a joint subjected to search combined bearing and bypass loads is complicated by the usual clearance between the hole and the fastener. A simple analysis method for such clearance-fit joints subjected to bearing-bypass loading has been developed in the present study. It uses an inverse formulation with a linear elastic finite-element analysis. Conditions along the bolt-hole contact arc are specified by displacement constraint equations. The present method is simple to apply and can be implemented with most general purpose finite-element programs since it does not use complicated iterative-incremental procedures. The method was used to study the effects of bearing-bypass loading on bolt-hole contact angles and local stresses. In this study, a rigid, frictionless bolt was used with a plate having the properties of a quasi-isotropic graphite/epoxy laminate. Results showed that the contact angle as well as the peak stresses around the hole and their locations were strongly influenced by the ratio of bearing and bypass loads. For single contact, tension and compression bearing-bypass loading had opposite effects on the contact angle. For some compressive bearing-bypass loads, the hole tended to close on the fastener leading to dual contact. It was shown that dual contact reduces the stress concentration at the fastener and would, therefore, increase joint strength in compression. The results illustrate the general importance of accounting for bolt-hole clearance and contact to accurately compute local bolt-hole stresses for combined bearings and bypass loading.

Naik, R. A.↗

A new explicit method for the numerical solution of parabolic differential equations

A new method is derived for solving parabolic partial differential equations arising in transient heat conduction or in boundary-layer flows. The method is based on a combination of the modified differential quadrature (MDQ) method with the rational Runge-Kutta time-integration scheme. It is fully explicit, requires no matrix inversion, and is stable for any time-step for the heat equations. Burgers equation and the one- and two-dimensional heat equations are solved to demonstrate the accuracy and efficiency of the proposed algorithm. The present method is found to be very accurate and efficient when results are compared with analytic solutions.

Satofuka, N.↗

A spectral domain method for remotely probing stratified media

The problem of remotely probing a stratified, lossless, dielectric medium is formulated using the spectral domain method of probing. The response of the medium to a spectrum of plane waves incident at various angles is used to invert the unknown profile. For TE polarization, the electric field satisfies a Helmholtz equation. The inverse problem is solved by means of a new representation for the wave function. The principal step in this inversion is solving a second kind Fredholm equation which is very amenable to numerical computations. Several examples are presented including some which indicate that the method can be used with experimentally obtained data. When the fields exhibit a surface wave behavior, a unique inversion can be obtained only if information about the magnetic field is also available. In this case, the inversion is accomplished by a two-step procedure which employs a formula of Jost and Kohn. Some examples are presented, and an approach which greatly shortens the computations without greatly deteriorating the results is discussed.

Schaubert, D. H.↗

Development of replicated optics for AXAF-1 XDA testing

Advanced optical systems for applications such as grazing incidence Wolter I x-ray mirror assemblies require extraordinary mirror surfaces in terms of fine finish and surface figure. The impeccable mirror surface is on the inside of the rotational mirror form. One practical method of producing devices with these requirements is to first fabricate an exterior surface for the optical device then replicate that surface to have the inverse component with lightweight characteristics. The replicated optic is not better than the master or mandrel from which it is made. This task identifies methods and materials for forming these extremely low roughness optical components. The objectives of this contract were to (1) prepare replication samples of electroless nickel coated aluminum, and determine process requirements for plating XDA test optic; (2) prepare and assemble plating equipment required to process a demonstration optic; (3) characterize mandrels, replicas and test samples for residual stress, surface contamination and surface roughness and figure using equipment at MSFC and; (4) provide technical expertise in establishing the processes, procedures, supplies and equipment needed to process the XDA test optics.

Engelhaupt, Darell↗

Aerodynamic Shape Optimization Using Hybridized Differential Evolution

An aerodynamic shape optimization method that uses an evolutionary algorithm known at Differential Evolution (DE) in conjunction with various hybridization strategies is described. DE is a simple and robust evolutionary strategy that has been proven effective in determining the global optimum for several difficult optimization problems. Various hybridization strategies for DE are explored, including the use of neural networks as well as traditional local search methods. A Navier-Stokes solver is used to evaluate the various intermediate designs and provide inputs to the hybrid DE optimizer. The method is implemented on distributed parallel computers so that new designs can be obtained within reasonable turnaround times. Results are presented for the inverse design of a turbine airfoil from a modern jet engine. (The final paper will include at least one other aerodynamic design application). The capability of the method to search large design spaces and obtain the optimal airfoils in an automatic fashion is demonstrated.

Madavan, Nateri K.↗

Quantitative Rainbow Schlieren Deflectometry as a Temperature Diagnostic for Spherical Flames

Numerical analysis and experimental results are presented to define a method for quantitatively measuring the temperature distribution of a spherical diffusion flame using Rainbow Schlieren Deflectometry in microgravity. First, a numerical analysis is completed to show the method can suitably determine temperature in the presence of spatially varying species composition. Also, a numerical forward-backward inversion calculation is presented to illustrate the types of calculations and deflections to be encountered. Lastly, a normal gravity demonstration of temperature measurement in an axisymmetric laminar, diffusion flame using Rainbow Schlieren deflectometry is presented. The method employed in this paper illustrates the necessary steps for the preliminary design of a Schlieren system. The largest deflections for the normal gravity flame considered in this paper are 7.4 x 10(-4) radians which can be accurately measured with 2 meter focal length collimating and decollimating optics. The experimental uncertainty of deflection is less than 5 x 10(-5) radians.

Feikema, Douglas A.↗

Shaping Diffraction-Grating Grooves to Optimize Efficiency

A method of shaping diffraction-grating grooves to optimize the spectral efficiency, spectral range, and image quality of a spectral imaging instrument is under development. The method is based on the use of an advanced design algorithm to determine the possibly complex shape of grooves needed to obtain a desired efficiency-versus-wavelength response (see figure). Then electron- beam fabrication techniques are used to realize the required groove shape. The method could be used, for example, to make the spectral efficiency of the grating in a given wavelength range proportional to the inverse of the spectral efficiency of a photodetector array so that the overall spectral efficiency of the combination of the grating and the photodetector array would be flat. The method has thus far been applied to one-dimensional gratings only, but in principle, it is also applicable to two-dimensional gratings. The algorithm involves calculations in the spatial-frequency domain. The spatial-frequency spectrum of a grating is represented as a diffraction-order spectral-peak-width function multiplied by an efficiency function for a single grating groove. This representation affords computational efficiency and accuracy by making it possible to consider only the response from one grating groove (one period of the grating), instead of from the whole grating area, in determining the response from the entire grating. This combination of efficiency and accuracy is crucial for future extensions of the algorithm to two-dimensional designs and to designs in which polarization must also be taken into account. The algorithm begins with the definition of target values of relative efficiency that represent the desired spectral response of the grating in certain spectral frequencies calculated from the diffraction order and wavelength. The grating period is divided into a number of cells - typically, 100. The phase contribution from each cell is determined from the phase of the incident electromagnetic wave and the height of the grating surface in the cell. The total contribution from all cells to each target value is then calculated. Then a method known to specialists as the optimum-rotation-angle method is used to adjust the height of each cell so that the total response from all cells is optimized. The computation is iterative and continues until the desired response is obtained. In the event that the desired response is unphysical, the algorithm nevertheless strives to generate a grating-grove profile for which the response approximates the desired one as closely as possible.

Backlund, John↗

A Model-Based Probabilistic Inversion Framework for Wire Fault Detection Using TDR

Time-domain reflectometry (TDR) is one of the standard methods for diagnosing faults in electrical wiring and interconnect systems, with a long-standing history focused mainly on hardware development of both high-fidelity systems for laboratory use and portable hand-held devices for field deployment. While these devices can easily assess distance to hard faults such as sustained opens or shorts, their ability to assess subtle but important degradation such as chafing remains an open question. This paper presents a unified framework for TDR-based chafing fault detection in lossy coaxial cables by combining an S-parameter based forward modeling approach with a probabilistic (Bayesian) inference algorithm. Results are presented for the estimation of nominal and faulty cable parameters from laboratory data.

Schuet, Stefan R.↗

A Comparison of Control Allocation Methods in the Presence of Parametric Model Uncertainty

When allocating redundant effectors to virtual control commands, linear (generalized inverse) allocators have historically been used on aircraft and spacecraft. While simple to implement, generalized inverses are unable to realize a significant portion of the attainable moments. To address this drawback, the control allocation problem can also be formulated as a linear programming or quadratic programming problem and solved using convex optimization based solvers. These approaches have been shown to access a larger set of attainable moments, however, little work has been done to understand the performance of convex optimization-based control allocation in the presence of parametric model uncertainty. This paper seeks to compare the performance of several control allocation approaches, including two forms of generalized inverse allocators, the pseudo inverse and minimum-variance, and the linear programming and quadratic programming approach in the presence of parametric model uncertainty. The performance of these four allocators were tested on an aircraft model in the presence of realistic parametric model uncertainty and the convex optimization approaches were shown to outperform the generalized inverses.

Luke J Miller↗

Recursive inverse kinematics for robot arms via Kalman filtering and Bryson-Frazier smoothing

This paper applies linear filtering and smoothing theory to solve recursively the inverse kinematics problem for serial multilink manipulators. This problem is to find a set of joint angles that achieve a prescribed tip position and/or orientation. A widely applicable numerical search solution is presented. The approach finds the minimum of a generalized distance between the desired and the actual manipulator tip position and/or orientation. Both a first-order steepest-descent gradient search and a second-order Newton-Raphson search are developed. The optimal relaxation factor required for the steepest descent method is computed recursively using an outward/inward procedure similar to those used typically for recursive inverse dynamics calculations. The second-order search requires evaluation of a gradient and an approximate Hessian. A Gauss-Markov approach is used to approximate the Hessian matrix in terms of products of first-order derivatives. This matrix is inverted recursively using a two-stage process of inward Kalman filtering followed by outward smoothing. This two-stage process is analogous to that recently developed by the author to solve by means of spatial filtering and smoothing the forward dynamics problem for serial manipulators.

Rodriguez, G.↗

Impact of left ventricular outflow tract area on systolic outflow velocity in hypertrophic cardiomyopathy: a real-time three-dimensional echocardiographic study

OBJECTIVES: The aim of this study was to use real-time three-dimensional echocardiography (3DE) to investigate the quantitative relation between minimal left ventricular (LV) outflow tract area (A(LVOT)) and maximal LV outflow tract (LVOT) velocity in patients with hypertrophic obstructive cardiomyopathy (HCM). BACKGROUND: In patients with HCM, LVOT velocity should change inversely with minimal A(LVOT) unless LVOT obstruction reduces the pumping capacity of the ventricle. METHODS: A total of 25 patients with HCM with systolic anterior motion (SAM) of the mitral valve leaflets underwent real-time 3DE. The smallest A(LVOT) during systole was measured using anatomically oriented two-dimensional "C-planes" within the pyramidal 3DE volume. Maximal velocity across LVOT was evaluated by two-dimensional Doppler echocardiography (2DE). For comparison with 3DE A(LVOT), the SAM-septal distance was determined by 2DE. RESULTS: Real-time 3DE provided unique information about the dynamic SAM-septal relation during systole, with A(LVOT) ranging from 0.6 to 5.2 cm(2) (mean: 2.2 +/- 1.4 cm(2)). Maximal velocity (v) correlated inversely with A(LVOT) (v = 496 A(LVOT)(-0.80), r = -0.95, p < 0.001), but the exponent (-0.80) was significantly different from -1.0 (95% confidence interval: -0.67 to -0.92), indicating a significant impact of small A(LVOT) on the peak LVOT flow rate. By comparison, the best correlation between velocity and 2DE SAM-septal distance was significantly (p < 0.01) poorer at -0.83, indicating the superiority of 3DE for assessing A(LVOT). CONCLUSIONS: Three-dimensional echocardiography-measured A(LVOT) provides an assessment of HCM geometry that is superior to 2DE methods. These data indicate that the peak LVOT flow rate appears to be significantly decreased by reduced A(LVOT). Real-time 3DE is a potentially valuable clinical tool for assessing patients with HCM.

NASA Discipline Cardiopulmonary↗

Model Evaluation of Methods for Estimating Surface Emissions and Chemical Lifetimes from Satellite Data

Column densities from satellite retrievals can provide valuable information for estimating emissions and chemical lifetimes objectively across the globe. To better understand the uncertainties associated with these estimates, we test four methods using simulated column densities from a point source: a box model approach, a 2D Gaussian fit, an Inverse Radius fit and an Exponentially-Modified Gaussian fit. The model results were simulated using the WRF and CAMx models for the year 2005, for a single point source outside Atlanta in Georgia, USA with specified emissions and three chemical scenarios: no chemical reactions, 12 h chemical lifetime and 1 h chemical lifetime. No other sources were included in the simulations. We find that the box model provides reliable estimates irrespective of plume speed and plume direction, if the plume speed and the chemical lifetime are known accurately. The 2D Gaussian fit was found to be sensitive to plume speed and direction, and requires omnidirectional dispersion in order to have a decent fit. However, the 2D Gaussian fit is only an approximate fit to the data, and the discrepancies mean that the results are dependent on the geographical domain used for the optimization. An Inverse Radius fit is introduced to correct this issue, which is found to provide improved emissions and lifetime estimates. The Exponentially-Modified Gaussian fit also gave improved estimates. It is however dependent on accurate plume rotation such that reported chemical lifetimes with this method could be significantly underestimated.

Emission inventory↗

A spectral element approach to wave motion in layered solids

A matrix methodology similar to that of the finite element method is developed for the analysis of stress waves in layered solids. Because the mass distribution is modeled exactly, the approach gives the exact frequency response of each layer. The fast Fourier transform and Fourier series are used for inversion to the time/space domain. The impact of a structured medium with multiple layers is used to demonstrate the method. Comparison with existing propagator and direct global matrix methods show the present approach to be computationally more efficient.

Rizzi, S. A.↗

A new method for analyzing IRAS data to determine the dust temperature distribution

In attempting to analyze the four-band IRAS images of interstellar dust emission, it is found that an inversion theorem recently developed by Chen (1990) enables distribution of the dust to be determined as a function of temperature and thus the total dust column density, for each line of sight. The method and its application to a hypothetical IRAS data set created by assuming a power-law dust temperature distribution, which is characteristic of the actual IRAS data for the Monoceros R2 cloud, are reported. To use the method, the wavelength dependence of the dust emissivity is assumed and a simple function is fitted to the four intensity-wavelength data points. The method is shown to be very successful at retrieving the dust temperature distribution in this case and is expected to have wide applicability to astronomical problems of this type.

Xie, Taoling↗

Inversion of parameters for semiarid regions by a neural network

Microwave brightness temperatures obtained from a passive radiative transfer model are inverted through use of a neural network. The model is applicable to semiarid regions and produces dual-polarized brightness temperatures for 6.6-, 10.7-, and 37-GHz frequencies. A range of temperatures is generated by varying three geophysical parameters over acceptable ranges: soil moisture, vegetation moisture, and soil temperature. A multilayered perceptron (MLP) neural network is trained with a subset of the generated temperatures, and the remaining temperatures are inverted using a backpropagation method. Several synthetic terrains are devised and inverted by the network under local constraints. All the inversions show good agreement with the original geophysical parameters, falling within 5 percent of the actual value of the parameter range.

Zurk, Lisa M.↗

Sky and Elemental Planetary Mapping Via Gamma Ray Emissions

Low-energy gamma ray emissions ((is) approximately 30keV to (is) approximately 30MeV) are significant to astrophysics because many interesting objects emit their primary energy in this regime. As such, there has been increasing demand for a complete map of the gamma ray sky, but many experiments to do so have encountered obstacles. Using an innovative method of applying the Radon Transform to data from BATSE (the Burst And Transient Source Experiment) on NASA's CGRO (Compton Gamma-Ray Observatory) mission, we have circumvented many of these issues and successfully localized many known sources to 0.5 - 1 deg accuracy. Our method, which is based on a simple 2-dimensional planar back-projection approximation of the inverse Radon transform (familiar from medical CAT-scan technology), can thus be used to image the entire sky and locate new gamma ray sources, specifically in energy bands between 200keV and 2MeV which have not been well surveyed to date. Samples of these results will be presented. This same technique can also be applied to elemental planetary surface mapping via gamma ray spectroscopy. Due to our method's simplicity and power, it could potentially improve a current map's resolution by a significant factor.

BATSE (the Burst and Transient Source Experiment)↗

Free vibration analysis of spinning flexible space structures

Efficient computation of natural frequencies and associated free vibration modes of spinning flexible structures is required to accurately determine the nature of interaction between the flexible structure and the attitude control system, which is vital in relating control torques to attitude angles. While structural discretization is effected by the finite element method, the resulting eigenvalue problem is solved by a combined Sturm sequence and inverse iteration procedure that yields a few specified roots and associated vectors. An eigenvalue procedure, based on a simultaneous iteration technique, provides efficient computation of the first few roots and vectors; a modal synthesis procedure proves to be useful for eigenproblem solutions of unusual structures such as spacecraft.

Gupta, K. K.↗

Free vibration analysis of spinning flexible space structures

Efficient computation of natural frequencies and associated free vibration modes of spinning flexible structures is required to accurately determine the nature of interaction between the flexible structure and the attitude control system, which is vital in relating control torques to attitude angles. While structural discretization is effected by the finite element method, the resulting eigenvalue problem is solved by a combined Sturm sequence and inverse iteration procedure that yields a few specified roots and associated vectors. An eigenvalue procedure, based on a simultaneous iteration technique, provides efficient computation of the first few roots and vectors; a modal synthesis procedure proves to be useful for eigenproblem solutions of unusual structures such as spacecraft.

Gupta, K. K.↗