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 883 records · Page 49

Geopotential Field Anomaly Continuation with Multi-Altitude Observations

Conventional gravity and magnetic anomaly continuation invokes the standard Poisson boundary condition of a zero anomaly at an infinite vertical distance from the observation surface. This simple continuation is limited, however, where multiple altitude slices of the anomaly field have been observed. Increasingly, areas are becoming available constrained by multiple boundary conditions from surface, airborne, and satellite surveys. This paper describes the implementation of continuation with multi-altitude boundary conditions in Cartesian and spherical coordinates and investigates the advantages and limitations of these applications. Continuations by EPS (Equivalent Point Source) inversion and the FT (Fourier Transform), as well as by SCHA (Spherical Cap Harmonic Analysis) are considered. These methods were selected because they are especially well suited for analyzing multi-altitude data over finite patches of the earth such as covered by the ADMAP database. In general, continuations constrained by multi-altitude data surfaces are invariably superior to those constrained by a single altitude data surface due to anomaly measurement errors and the non-uniqueness of continuation.

Kim, Jeong Woo↗

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

A simple direct stress analysis method is proposed for a laminate with a clearance-fit fastener subjected to combined bearing and bypass loads in tension or compression, including bearing and bypass loads on bolt-hole contact and local stresses. The approach uses a linear-elastic finite element analysis with an inverse formulation; conditions along the bolt-hole interface are specified by constraint equations that limit nodal displacements to a circular arc corresponding to the bolt diameter. The application of the method to the analysis of the effects of bearing-bypass loading on bolt-hole contact angles and local stresses, with the plate material properties representing a quasi-isotropic T300/5208 graphite/epoxy laminate, is discussed.

Naik, R. A.↗

Profile inversion in presence of ray bending

Inversion of radio occultation data for planetary atmospheres and ionospheres has been performed using the seismological Herglotz-Wiechert method, as adapted by Phinney and Anderson to the radio-occultation case. Profile reconstruction performed in computer simulated experiments with this approach have been compared with the ones obtained with the straight-ray Abel transform. For a thin atmosphere and ionosphere, like the ones encountered on Mars, microwave occultation data can be inverted accurately with both methods. For a dense ionosphere like the sun's corona, ray bending of microwaves is severe, and recovered refractivity by the Herglotz-Wiechert method provides significant improvement over the straight-ray Abel transform: the error reduces from more than 60% to less than 20% at a height of 60,000 km above the base of the corona.

Wallio, H. A.↗

Contrast reduction by the atmosphere and retrieval of nonuniform surface reflectance

A radiative transfer model is developed which gives the upward radiance at nadir for any 1-D Lambertian surface reflectance. This model is used to depict the atmospheric effect on the transmittance of contrast for any 1-D surface reflectance. Here by contrast we mean a general variation of the radiation field across the image. With the aid of this model an inversion algorithm is developed for retrieval of true surface reflectance from high resolution satellite data (e.g., Landsat). This inversion technique can be a useful tool for extraction of surface reflectance from satellite data in the case of a surface reflectance variable in one dimension only (e.g., seashore or near borders of big fields). A sensitivity study of the inversion procedure on the knowledge of atmospheric parameters and sensor calibration was performed. It is shown that this inversion technique is stable even in the presence of errors in the sensor calibration and the atmospheric parameters. The method was applied to Landsat data in two wavelengths. The results show reasonable dependence of the derived surface reflectance on the distance from the seashore.

Mekler, Y.↗

Method of examining microcircuit patterns

Examination of microstructures of LSI and VLSI devices is facilitated by employing a method in which the device is photographed through a darkfield illumination optical microscope and the resulting negative subjected to inverse processing to form a positive on a photographic film. The film is then developed to form photographic prints or transparencies which clearly illustrate the structure of the device. The entire structure of a device may be examined by alternately photographing the device and selectively etching layers of the device in order to expose underlying layers.

Suszko, S. F.↗

TRANDESNF: A computer program for transonic airfoil design and analysis in nonuniform flow

The use of a transonic airfoil code for analysis, inverse design, and direct optimization of an airfoil immersed in propfan slipstream is described. A summary of the theoretical method, program capabilities, input format, output variables, and program execution are described. Input data of sample test cases and the corresponding output are given.

Chang, J. F.↗

Monte Carlo renormalization-group investigation of the two-dimensional O(4) sigma model

An improved Monte Carlo renormalization-group method is used to determine the beta function of the two-dimensional O(4) sigma model. While for (inverse) couplings beta = greater than about 2.2 agreement is obtained with asymptotic scaling according to asymptotic freedom, deviations from it are obtained at smaller couplings. They are, however, consistent with the behavior of the correlation length, indicating 'scaling' according to the full beta function. These results contradict recent claims that the model has a critical point at finite coupling.

Heller, Urs M.↗

Heat Flow and Magnetization in the Oceanic Lithosphere

Two aspects of the processing and interpretation of satellite measurements of the geomagnetic field are described. One deals with the extraction of the part of the geomagnetic field that originates from sources in the earth's atmosphere. The other investigates the possibility of using the thermal state of the oceanic lithosphere to further constrain modelling and interpretation of magnetic anomalies. It is shown that some of the magnetic signal in crustal anomaly maps can be an artifact of the mathematical algorithms that have been used to separate the crustal field from the observed data. Strong magnetic anomalies can be distorted but are probably real, but weak magnetic anomalies can arise from leakage of power from short wavelengths, and will also appear in anomaly maps as repetitions of the strong crustal anomaly. The distortion and the ghost anomalies follow the magnetic dip lines in a way that is similar to actual MAGSAT anomaly fields. This phenomenon will also affect the lower degree spherical harmonic terms in the power spectrum of the crustal field. A model of the magnetic properties of the oceanic crust that has been derived from direct measurements of the rock magnetic properties of oceanic rocks is presented. The average intensity of magnetization in the oceanic crust is not strong enough to explain magnetic anomalies observed over oceanic areas. This is the case for both near surface observations (ship and aeromagnetic data) and satellite altitude observations. It is shown that magnetic sources in the part of the upper mantle that is situated above the Curie isotherm, if sufficiently strong, can produce satellite magnetic anomalies that are comparable to MAGSAT data. The method developed for the study of depth to the Curie isotherm and magnetic anomalies can also be used in inverse modelling of satellite magnetic anomalies when the model is to be adjusted with an annihilator.

Hayling, Kjell Lennart↗

Mesh-matrix analysis method for electromagnetic launchers

The mesh-matrix method is a procedure for calculating the current distribution in the conductors of electromagnetic launchers with coil or flat-plate geometry. Once the current distribution is known the launcher performance can be calculated. The method divides the conductors into parallel current paths, or meshes, and finds the current in each mesh by matrix inversion. The author presents procedures for writing equations for the current and voltage relations for a few meshes to serve as a pattern for writing the computer code. An available subroutine package provides routines for field and flux coefficients and equation solution.

Elliott, David G.↗

Inverse problems in diffraction

A two-dimensional problem of diffraction of a plane electromagnetic wave on a smooth 2 pi-periodic surface is considered. A numerical algorithm solving this problem is developed. An inverse problem of determination of the shape of 2 pi-periodic surface using the performance data of reverse scattering is considered. The inverse problem was solved by means of minimization of the residual functional with the help of the gradient descent method. The initial data were calculated with the help of the numerical method. On each step of the iterative method of minimization, the residual functional was calculated approximately with the help of the small slope method. The examples of the shape determination are considered.

Mikheev, Andrew G.↗

Integrals and integral equations in linearized wing theory

The formulas of subsonic and supersonic wing theory for source, doublet, and vortex distributions are reviewed and a systematic presentation is provided which relates these distributions to the pressure and to the vertical induced velocity in the plane of the wing. It is shown that care must be used in treating the singularities involved in the analysis and that the order of integration is not always reversible. Concepts suggested by the irreversibility of order of integration are shown to be useful in the inversion of singular integral equations when operational techniques are used. A number of examples are given to illustrate the methods presented, attention being directed to supersonic flight speed.

Lomax, Harvard↗

Radiative interactions in chemically reacting compressible nozzle flows using Monte Carlo simulations

The two-dimensional spatially elliptic Navier-Stokes equations have been used to investigate the radiative interactions in chemically reacting compressible flows of premixed hydrogen and air in an expanding nozzle. The radiative heat transfer term in the energy equation is simulated using the Monte Carlo method (MCM). The nongray model employed is based on the statistical narrow band model with an exponential-tailed inverse intensity distribution. The spectral correlation has been considered in the Monte Carlo formulations. Results obtained demonstrate that the effect of radiation on the flow field is minimal but its effect on the wall heat transfer is significant. Extensive parametric studies are conducted to investigate the effects of equivalence ratio, wall temperature, inlet flow temperature, and the nozzle size on the radiative and conductive wall fluxes.

Liu, J.↗

Simultaneous Planck, Swift, and Fermi Observations of X-ray and Gamma-ray Selected Blazars

We present simultaneous Planck, Swift, Fermi, and ground-based data for 105 blazars belonging to three samples with flux limits in the soft X-ray, hard X-ray, and gamma-ray bands, with additional 5 GHz flux-density limits to ensure a good probability of a Planck detection. We compare our results to those of a companion paper presenting simultaneous Planck and multi-frequency observations of 104 radio-loud northern active galactic nuclei selected at radio frequencies. While we confirm several previous results, our unique data set allows us to demonstrate that the selection method strongly influences the results, producing biases that cannot be ignored. Almost all the BL Lac objects have been detected by the Fermi Large Area Telescope (LAT), whereas 30% to 40% of the flat-spectrum radio quasars (FSRQs) in the radio, soft X-ray, and hard X-ray selected samples are still below the gamma-ray detection limit even after integrating 27 months of Fermi-LAT data. The radio to sub-millimetre spectral slope of blazars is quite flat, with (alpha) approx 0 up to about 70GHz, above which it steepens to (alpha) approx -0.65. The BL Lacs have significantly flatter spectra than FSRQs at higher frequencies. The distribution of the rest-frame synchrotron peak frequency (nu(sup s)(sub peak)) in the spectral energy distribution (SED) of FSRQs is the same in all the blazar samples with (nu(sup s)(sub peak)) = 10(exp 13.1 +/- 0.1) Hz, while the mean inverse Compton peak frequency, (nu(sup IC)(sub peak)), ranges from 10(exp 21) to 10(exp 22) Hz. The distributions of nu(sup s)(sub peak) and nu(sup IC)(sub peak) of BL Lacs are much broader and are shifted to higher energies than those of FSRQs; their shapes strongly depend on the selection method. The Compton dominance of blazars. defined as the ratio of the inverse Compton to synchrotron peak luminosities, ranges from less than 0.2 to nearly 100, with only FSRQs reaching values larger than about 3. Its distribution is broad and depends strongly on the selection method, with gamma-ray selected blazars peaking at approx 7 or more, and radio-selected blazars at values close to 1, thus implying that the common assumption that the blazar power budget is largely dominated by high-energy emission is a selection effect. A comparison of our multi-frequency data with theoretical predictions shows that simple homogeneous SSC models cannot explain the simultaneous SEDs of most of the gamma-ray detected blazars in all samples. The SED of the blazars that were not detected by Fermi~LAT may instead be consistent with SSC emission. Our data challenge the correlation between bolometric luminosity and nu(sup s)(sub peak) predicted by the blazar sequence.

Giommi, P.↗

Krylov Subspace Methods for Complex Non-Hermitian Linear Systems

We consider Krylov subspace methods for the solution of large sparse linear systems Ax = b with complex non-Hermitian coefficient matrices. Such linear systems arise in important applications, such as inverse scattering, numerical solution of time-dependent Schrodinger equations, underwater acoustics, eddy current computations, numerical computations in quantum chromodynamics, and numerical conformal mapping. Typically, the resulting coefficient matrices A exhibit special structures, such as complex symmetry, or they are shifted Hermitian matrices. In this paper, we first describe a Krylov subspace approach with iterates defined by a quasi-minimal residual property, the QMR method, for solving general complex non-Hermitian linear systems. Then, we study special Krylov subspace methods designed for the two families of complex symmetric respectively shifted Hermitian linear systems. We also include some results concerning the obvious approach to general complex linear systems by solving equivalent real linear systems for the real and imaginary parts of x. Finally, numerical experiments for linear systems arising from the complex Helmholtz equation are reported.

Freund, Roland W.↗

Implementation and Testing of Inverse Kinematics on Robotic Arm

COSIE (Coronal Spectrographic Imager in the Extreme Ultraviolet) is a proposed solar tracking ISS imaging payload that will help bridge the theoretical gap between the physics of the low corona and the heliosphere. This scientific instrument requires high pointing accuracy, on the order of arc seconds. The instrument is mounted on to a three revolute joint robotic arm in order to track the roll, pitch and yaw motion of the Sun. The goal of this project is to construct a prototype model of the robotic arm and implement the proposed analytical inverse kinematics algorithm. In robotics, the inverse kinematics problem is solving for the set of joint angles that achieve the desired end effect or location and/or orientation. In this case, orientation is the focus. Depending on the configuration, multiple sets of joint angle solutions may exist. Due to the complexity of robotics, typically iterative methods are used to solve for the joint angle solution sets. However, in this case, an analytical solution exists. A small robotic arm representative of the full size hardware was constructed. The inverse kinematics algorithm, originally in MATLAB/Simulink, was converted into C in order to interface with the motors. This C software was implemented on a Windows PC and micro-controller, and serial communication between the two was established, allowing the motors to be directly controlled by the inverse kinematics algorithm. Testing the inverse kinematics on a physical system will allow the validity and accuracy of the analytic solution to be verified.

Franz, Carter↗

Maxwell's second- and third-order equations of transfer for non-Maxwellian gases

Condensed algebraic forms for Maxwell's second- and third-order equations of transfer are developed for the case of molecules described by either elastic hard spheres, inverse-power potentials, or by Bird's variable hard-sphere model. These hardly reduced, yet exact, equations provide a new point of origin, when using the moment method, in seeking approximate solutions in the kinetic theory of gases for molecular models that are physically more realistic than that provided by the Maxwell model. An important by-product of the analysis when using these second- and third-order relations is that a clear mathematical connection develops between Bird's variable hard-sphere model and that for the inverse-power potential.

Baganoff, D.↗

An asymptotic method for estimating the vertical ozone distribution in the Earth's atmosphere from satellite measurements of backscattered solar UV-radiation

An asymptotic approach to solution of the inverse problems of remote sensing is presented. It consists in changing integral operators characteristic of outgoing radiation into their asymptotic analogues. Such approach does not add new principal uncertainties into the problem and significantly reduces computation time that allows to develop the real (or about) time algorithms for interpretation of satellite measurements. The asymptotic approach has been realized for estimating vertical ozone distribution from satellite measurements of backscatter solar UV radiation in the Earth's atmosphere.

Ishov, Alexander G.↗