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 1,099 records · Page 61

Adaptive Neuron Model: An architecture for the rapid learning of nonlinear topological transformations

A method for the rapid learning of nonlinear mappings and topological transformations using a dynamically reconfigurable artificial neural network is presented. This fully-recurrent Adaptive Neuron Model (ANM) network was applied to the highly degenerate inverse kinematics problem in robotics, and its performance evaluation is bench-marked. Once trained, the resulting neuromorphic architecture was implemented in custom analog neural network hardware and the parameters capturing the functional transformation downloaded onto the system. This neuroprocessor, capable of 10(exp 9) ops/sec, was interfaced directly to a three degree of freedom Heathkit robotic manipulator. Calculation of the hardware feed-forward pass for this mapping was benchmarked at approximately 10 microsec.

Tawel, Raoul↗

Harmonic oscillators and resonance series generated by a periodic unstable classical orbit

The presence of an unstable periodic classical orbit allows one to introduce the decay time as a purely classical magnitude: inverse of the Lyapunov index which characterizes the orbit instability. The Uncertainty Relation gives the corresponding resonance width which is proportional to the Planck constant. The more elaborate analysis is based on the parabolic equation method where the problem is effectively reduced to the multidimensional harmonic oscillator with the time-dependent frequency. The resonances form series in the complex energy plane which is equidistant in the direction perpendicular to the real axis. The applications of the general approach to various problems in atomic physics are briefly exposed.

Kazansky, A. K.↗

Inverse Compton Scattering in Mildly Relativistic Plasma

We investigated the effect of inverse Compton scattering in mildly relativistic static and moving plasmas with low optical depth using Monte Carlo simulations, and calculated the Sunyaev-Zel'dovich effect in the cosmic background radiation. Our semi-analytic method is based on a separation of photon diffusion in frequency and real space. We use Monte Carlo simulation to derive the intensity and frequency of the scattered photons for a monochromatic incoming radiation. The outgoing spectrum is determined by integrating over the spectrum of the incoming radiation using the intensity to determine the correct weight. This method makes it possible to study the emerging radiation as a function of frequency and direction. As a first application we have studied the effects of finite optical depth and gas infall on the Sunyaev-Zel'dovich effect (not possible with the extended Kompaneets equation) and discuss the parameter range in which the Boltzmann equation and its expansions can be used. For high temperature clusters (k(sub B)T(sub e) greater than or approximately equal to 15 keV) relativistic corrections based on a fifth order expansion of the extended Kompaneets equation seriously underestimate the Sunyaev-Zel'dovich effect at high frequencies. The contribution from plasma infall is less important for reasonable velocities. We give a convenient analytical expression for the dependence of the cross-over frequency on temperature, optical depth, and gas infall speed. Optical depth effects are often more important than relativistic corrections, and should be taken into account for high-precision work, but are smaller than the typical kinematic effect from cluster radial velocities.

Molnar, S. M.↗

On Improving Efficiency of Differential Evolution for Aerodynamic Shape Optimization Applications

Differential Evolution (DE) is a simple and robust evolutionary strategy that has been provEn effective in determining the global optimum for several difficult optimization problems. Although DE offers several advantages over traditional optimization approaches, its use in applications such as aerodynamic shape optimization where the objective function evaluations are computationally expensive is limited by the large number of function evaluations often required. In this paper various approaches for improving the efficiency of DE are reviewed and discussed. Several approaches that have proven effective for other evolutionary algorithms are modified and implemented in a DE-based aerodynamic shape optimization method that uses a Navier-Stokes solver for the objective function evaluations. Parallelization techniques on distributed computers are used to reduce turnaround times. Results are presented for standard test optimization problems and for the inverse design of a turbine airfoil. The efficiency improvements achieved by the different approaches are evaluated and compared.

Madavan, Nateri K.↗

ISSM: Ice Sheet System Model

In order to have the capability to use satellite data from its own missions to inform future sea-level rise projections, JPL needed a full-fledged ice-sheet/iceshelf flow model, capable of modeling the mass balance of Antarctica and Greenland into the near future. ISSM was developed with such a goal in mind, as a massively parallelized, multi-purpose finite-element framework dedicated to ice-sheet modeling. ISSM features unstructured meshes (Tria in 2D, and Penta in 3D) along with corresponding finite elements for both types of meshes. Each finite element can carry out diagnostic, prognostic, transient, thermal 3D, surface, and bed slope simulations. Anisotropic meshing enables adaptation of meshes to a certain metric, and the 2D Shelfy-Stream, 3D Blatter/Pattyn, and 3D Full-Stokes formulations capture the bulk of the ice-flow physics. These elements can be coupled together, based on the Arlequin method, so that on a large scale model such as Antarctica, each type of finite element is used in the most efficient manner. For each finite element referenced above, ISSM implements an adjoint. This adjoint can be used to carry out model inversions of unknown model parameters, typically ice rheology and basal drag at the ice/bedrock interface, using a metric such as the observed InSAR surface velocity. This data assimilation capability is crucial to allow spinning up of ice flow models using available satellite data. ISSM relies on the PETSc library for its vectors, matrices, and solvers. This allows ISSM to run efficiently on any parallel platform, whether shared or distrib- ISSM: Ice Sheet System Model NASA's Jet Propulsion Laboratory, Pasadena, California uted. It can run on the largest clusters, and is fully scalable. This allows ISSM to tackle models the size of continents. ISSM is embedded into MATLAB and Python, both open scientific platforms. This improves its outreach within the science community. It is entirely written in C/C++, which gives it flexibility in its design, and the power/speed that C/C++ allows. ISSM is svn (subversion) hosted, on a JPL repository, to facilitate its development and maintenance. ISSM can also model propagation of rifts using contact mechanics and mesh splitting, and can interface to the Dakota software. To carry out sensitivity analysis, mesh partitioning algorithms are available, based on the Scotch, Chaco, and Metis partitioners that ensure equal area mesh partitions can be done, which are then usable for sampling and local reliability methods.

Larour, Eric↗

Modeling arson - An exercise in qualitative model building

A detailed example is given of the role of von Neumann and Morgenstern's 1944 'expected utility theorem' (in the theory of games and economic behavior) in qualitative model building. Specifically, an arsonist's decision as to the amount of time to allocate to arson and related activities is modeled, and the responsiveness of this time allocation to changes in various policy parameters is examined. Both the activity modeled and the method of presentation are intended to provide an introduction to the scope and power of the expected utility theorem in modeling situations of 'choice under uncertainty'. The robustness of such a model is shown to vary inversely with the number of preference restrictions used in the analysis. The fewer the restrictions, the wider is the class of agents to which the model is applicable, and accordingly more confidence is put in the derived results. A methodological discussion on modeling human behavior is included.

Heineke, J. M.↗

Calculation of the flow field in supersonic inlets using a bicharacteristics method with shock wave fitting

The results of weakly viscous flow analysis are presented. The flow field, including molecular transport, is computed with the aid of a bicharacteristic method. The bow shock wave and the internal shock wave are computed with the aid of a three-dimensional shock wave fitting procedure. Characteristic equations are presented, and numerical integration procedure is discussed. Here, an inverse marching scheme is employed in which the solution is obtained on space-like planes of constant x and on space curves defined by the intersections of the internal shock wave with the solid boundaries. The distance between solution planes is arrived at by the Courant-Friedrichs-Lewy stability criterion.

Vadyak, J.↗

Evaluation and development of new VAS remote sensing algorithms

An algorithm, developed to permit the simultaneous retrieval of temperature/moisture profiles and surface skin temperature by direct analytical solution of the radiative transfer equation, alleviates the problem associated with the interdependencies of water vapor retrieval on temperature retrieval and their dependencies on surface emissions. Simultaneous solution of all these quantities is achieved in one calculation using the available radiance observations. Since only a single matrix inversion is required for the specification of all parameters, the solution is computationally efficient. Ancillary observations of temperature and/or moisture from surface sensors on aircraft can be readily incorporated into the solution. Simulation tests of the method indicate improved performance over the previous iterative technique, particularly for the lower troposphere and for water vapor.

Smith, W. L.↗

Noniterative three-dimensional grid generation using parabolic partial differential equations

A new algorithm for generating three-dimensional grids has been developed and implemented which numerically solves a parabolic partial differential equation (PDE). The solution procedure marches outward in two coordinate directions, and requires inversion of a scalar tridiagonal system in the third. Source terms have been introduced to control the spacing and angle of grid lines near the grid boundaries, and to control the outer boundary point distribution. The method has been found to generate grids about 100 times faster than comparable grids generated via solution of elliptic PDEs, and produces smooth grids for finite-difference flow calculations.

Edwards, T. A.↗

Seismic analysis of the solar interior. I - Can opacity changes improve the theoretical frequencies?

The paper describes the application of seismic inverse theory to the deduction of properties of the solar interior using presently available measured frequencies of the solar oscillations. Only the solar opacity is included in this application. This study has used the spectral expansion method of Lanczos and Jackson to derive changes to the opacity which improve agreement between the theoretical and observed frequencies of oscillation. It is found that a family of opacity changes exists which yields models that improve the frequency agreement by amounts that are indistinguishable among the family members.

Korzennik, Sylvain G.↗

Flight Test of an Adaptive Controller and Simulated Failure/Damage on the NASA NF-15B

The method of flight-testing the Intelligent Flight Control System (IFCS) Second Generation (Gen-2) project on the NASA NF-15B is herein described. The Gen-2 project objective includes flight-testing a dynamic inversion controller augmented by a direct adaptive neural network to demonstrate performance improvements in the presence of simulated failure/damage. The Gen-2 objectives as implemented on the NASA NF-15B created challenges for software design, structural loading limitations, and flight test operations. Simulated failure/damage is introduced by modifying control surface commands, therefore requiring structural loads measurements. Flight-testing began with the validation of a structural loads model. Flight-testing of the Gen-2 controller continued, using test maneuvers designed in a sequenced approach. Success would clear the new controller with respect to dynamic response, simulated failure/damage, and with adaptation on and off. A handling qualities evaluation was conducted on the capability of the Gen-2 controller to restore aircraft response in the presence of a simulated failure/damage. Control room monitoring of loads sensors, flight dynamics, and controller adaptation, in addition to postflight data comparison to the simulation, ensured a safe methodology of buildup testing. Flight-testing continued without major incident to accomplish the project objectives, successfully uncovering strengths and weaknesses of the Gen-2 control approach in flight.

Buschbacher, Mark↗

Tracking Vector Magnetograms with the Magnetic Induction Equation

The differential affine velocity estimator (DAVE) that we developed in 2006 for estimating velocities from line-of-sight magnetograms is modified to directly incorporate horizontal magnetic fields to produce a differential affine velocity estimator for vector magnetograms (DAVE4VM). The DAVE4VM's performance is demonstrated on the synthetic data from the anelastic pseudospectral ANMHD simulations that were used in the recent comparison of velocity inversion techniques by Welsch and coworkers. The DAVE4VM predicts roughly 95% of the helicity rate and 75% of the power transmitted through the simulation slice. Intercomparison between DAVE4VM and DAVE and further analysis of the DAVE method demonstrates that line-of-sight tracking methods capture the shearing motion of magnetic footpoints but are insensitive to flux emergence - the velocities determined from line-of-sight methods are more consistent with horizontal plasma velocities than with flux transport velocities. These results suggest that previous studies that rely on velocities determined from line-of-sight methods such as the DAVE or local correlation tracking may substantially misrepresent the total helicity rates and power through the photosphere.

Schuck, P.↗

Overview of the CHarring Ablator Response (CHAR) Code

An overview of the capabilities of the CHarring Ablator Response (CHAR) code is presented. CHAR is a one-, two-, and three-dimensional unstructured continuous Galerkin finite-element heat conduction and ablation solver with both direct and inverse modes. Additionally, CHAR includes a coupled linear thermoelastic solver for determination of internal stresses induced from the temperature field and surface loading. Background on the development process, governing equations, material models, discretization techniques, and numerical methods is provided. Special focus is put on the available boundary conditions including thermochemical ablation, surface-to-surface radiation exchange, and flowfield coupling. Finally, a discussion of ongoing development efforts is presented.

Amar, Adam J.↗

Overview of the CHarring Ablator Response (CHAR) Code

An overview of the capabilities of the CHarring Ablator Response (CHAR) code is presented. CHAR is a one-, two-, and three-dimensional unstructured continuous Galerkin finite-element heat conduction and ablation solver with both direct and inverse modes. Additionally, CHAR includes a coupled linear thermoelastic solver for determination of internal stresses induced from the temperature field and surface loading. Background on the development process, governing equations, material models, discretization techniques, and numerical methods is provided. Special focus is put on the available boundary conditions including thermochemical ablation and contact interfaces, and example simulations are included. Finally, a discussion of ongoing development efforts is presented.

Amar, Adam J.↗

Method development for compensating temperature effects in pressure sensitive paint measurements

Pressure sensitive luminescent paints (PSP) have recently emerged as a viable technique for aerodynamic pressure measurements. The technique uses a surface coating which contains probe molecules that luminesce when excited by light of an appropriate wavelength. The photoluminescence of these materials is known to be quenched by the presence of molecular oxygen. Since oxygen is a fixed mole fraction of the air, the coating's luminescence intensity varies inversely with air pressure. Digital imaging of the luminescence varying across a coated surface produces a pressure distribution map over that surface. One difficulty encountered with this technique is the temperature effect on the luminescence intensity. Present PSP formulations have significant sensitivity to temperature. At the moment, the most practical way of correcting for temperature effects is to calibrate the paint in place at the operating temperatures by using a few well-placed pressure taps. This study is looking at development of temperature indicating coatings that can be applied and measured concurrently with PSP, and use the temperature measurement to compute the correct pressure. Two methods for this dual paint formulation are proposed. One method will use a coating that consists of temperature sensitive phosphors in a polymer matrix. This is similar in construction to PSP, except that the probe molecules used are selected primarily for their temperature sensitivity. Both organic phosphors (e.g., europium thenoyltrifluoroacetonate, bioprobes) and inorganic phosphors (e.g., Mg4(F)GeO6:Mn, La2O2S:Eu, Radelin Type phosphors, Sylvania Type phosphors) will be evaluated for their temperature sensing potential. The next method will involve a novel coating composing of five membered heterocyclic conducting polymers which are known to show temperature dependent luminescence (e.g., poly(3-alkylthiopene), poly(3-alkylselenophene), poly(3-alkylfuran)). Both methods will involve applying a bottom layer of temperature sensitive coating followed by a top coating of PSP. An oxygen-impermeable polymer can be used as the temperature sensitive coating matrix, or it can be layered in between the coatings to prevent oxygen quenching of the bottom coating's luminescence. The probe molecules of the coatings will be excited by a broad band of light, with the different emissions detected and measured at their distinct wavelengths. Developmental research of these coatings is still in progress; however, preliminary results look very promising.

Demandante, Carlo Greg N.↗

Orion-arm magnetic monopoles and gamma rays

It is shown that galactic-disk monopoles could yield air showers which are stochastically aligned with the local Orion-arm magnetic field and which are above the energies of similarly aligned galactic heavy nuclei. In the absence of a cosmic-ray anisotropy, this method sets limits on the flux of galactic monopoles and their production cross sections. Limits based on the cosmic-ray anisotropy at 10 to the 19th power eV are presented. The possibility of stochastic alignment of inverse Compton photons from monopoles and 2.7 K black-body photons is considered. The resulting monopole-flux limits, based on the same data as those calculated by Osborne (1970), differ significantly from Osborne's limits, where the mechanisms of spiral-arm alignment and galactic containment are ignored.

Tompkins, D. R., Jr.↗

The analysis and design of transonic two-element airfoil systems

The multiphase effort in the development of tools for the analysis and design of two-element airfoil systems, that is, airfoils with a slat or a flap at transonic speeds is described. The first phase involved the development of a method to compute the inviscid flow over such configurations. In the second phase the inviscid code was coupled to a boundary layer calculation program in order to compute the loss in performance due to viscous effects. An inverse code that constructs the airfoil system corresponding to a desired pressure distribution is described.

Volpe, G.↗

A general theory of two- and three-dimensional rotational flow in subsonic and transonic turbomachines

This report represents a general theory applicable to axial, radial, and mixed flow turbomachines operating at subsonic and supersonic speeds with a finite number of blades of finite thickness. References reflect the evolution of computational methods used, from the inception of the theory in the 50's to the high-speed computer era of the 90's. Two kinds of relative stream surfaces, S(sub 1) and S(sub 2), are introduced for the purpose of obtaining a three-dimensional flow solution through the combination of two-dimensional flow solutions. Nonorthogonal curvilinear coordinates are used for the governing equations. Methods of computing transonic flow along S(sub 1) and S(sub 2) stream surfaces are given for special cases as well as for fully three-dimensional transonic flows. Procedures pertaining to the direct solutions and inverse solutions are presented. Information on shock wave locations and shapes needed for computations are discussed. Experimental data from a Deutsche Forschungs- und Versuchsanstalt fur Luft- und Raumfahrt e.V. (DFVLR) rotor and from a Chinese Academy of Sciences (CAS) transonic compressor rotor are compared with the computed flow properties.

Wu, Chung-Hua↗