Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inverse problem”

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 829 records · Page 46

Fission-induced plasmas

The possibility of creating a plasma from fission fragments, and to utilize the energy of the particles to create population inversion that would lead to laser action is investigated. An investigation was made of various laser materials which could be used for nuclear-pumped lasing. The most likely candidate for a fissioning material in the gaseous form is uranium hexafluoride - UF6, and experiments were performed to investigate materials that would be compatible with it. One of the central problems in understanding a fission-induced plasma is to obtain a model of the electron behavior, and some preliminary calculations are presented. In particular, the rates of various processes are discussed. A simple intuitive model of the electron energy distribution function is also shown. The results were useful for considering a mathematical model of a nuclear-pumped laser. Next a theoretical model of a (3)He-Ar nuclear-pumped laser is presented. The theory showed good qualitative agreement with the experimental results.

Harries, W. L.↗

Computer modeling of inversion layer MOS solar cells and arrays

A two dimensional numerical model of the inversion layer metal insulator semiconductor (IL/MIS) solar cell is proposed by using the finite element method. The two-dimensional current flow in the device is taken into account in this model. The electrostatic potential distribution, the electron concentration distribution, and the hole concentration distribution for different terminal voltages are simulated. The results of simple calculation are presented. The existing problems for this model are addressed. Future work is proposed. The MIS structures are studied and some of the results are reported.

Ho, Fat Duen↗

Calibration of long term satellite ozone data sets using the space shuttle

Trends in atmospheric ozone continue to be an environmental concern. Drifts in satellite observations are the major obstacle in the detection of changes in global ozone over the long term. Careful re-analysis of satellite ozone data along with groundbased observations have more or less corroborated photochemical models which predict ozone depletion. However, there remains the margin of error in the observations that is as large as the trend itself. The National Plan for Stratospheric Monitoring calls for monitoring global ozone for at lease the next ten years employing the NOAA polar orbiting satellites. Ozone observations will be made with the Solar Backscatter Ultraviolet Spectral Radiometer Mod 2 (SBUV/2) which ia a refinement of the SBUV instrument flying on NASA's Nimbus-7 satellite. The first instrument in the operational series began taking data from the NOAA-9 spacecraft in February 1985. A second instrument was launched on NOAA-11 in September 1988. Both continue to operate. Earlier attempts to calibrate satellite data relied on comparisons with ground based observations. However, differences in instrumental techniques severely complicated these efforts. This problem will be over come by regular flights, about once per year, of the Shuttle Solar Backscatter Ultraviolet radiometer (SSBUV). The data from the SSBUV instrument will be compared with nearly coincident data taken by the NOAA satellite instruments. This procedure will permit a direct calibration transfer in space, since the two instruments observe the same qualities thereby bypassing the inversion algorithm which converts the observations to ozone amounts.

Hilsenrath, Ernest↗

Dynamics in Layer Models of Solid Flame Propagation

Self-propagating high-temperature synthesis (SHS) is a process in which combustion waves, e.g., "solid flames", which are considered here, are employed to synthesize desired materials. Like many other systems, SHS is a pattern forming system. The problem of describing experimentally observed patterns and of predicting new, as yet unobserved, patterns continues to attract the attention of scientists and mathematicians due to the fundamental significance of the phenomena in combustion in particular, and in nonlinear science in general. Here, we analyze the dynamics of solid flame propagation in a 2D region by considering the region to be composed of parallel, identical layers aligned along the direction of propagation and having thermal contact. Each layer is then described by wave propagation in 1D, with the transverse Laplacian replaced by a term describing heat exchange between neighboring layers. This configuration is the simplest model of a 2D system because it accounts, in a simple way, for the principal feature of the problem, i.e., heat exchange between neighbors in the transverse direction. For simplicity, we describe the situation for two layers. Because the layers are identical, uniformly propagating waves in each layer must be identical, independent of the heat exchange rate alpha. When the Zeldovich number Z exceeds a critical value Z(sub c), which depends on alpha, uniformly propagating waves become unstable. The stability diagram for the two coupled layers reproduces that for the full 2D problem after appropriate identification of parameters in the two problems. Depending on parameter values, we determine three different steady-state dynamical behaviors (though additional behaviors are also expected to occur). The three behaviors are: (i) waves in each layer which pulsate in phase as they propagate, so that together they form a single pulsating propagating wave; (ii) the waves in each layer are no longer identical, and antiphase pulsations occur, with the two waves alternately advancing and receding as they propagate. This wave is the analog of the spinning wave on the surface of a circular cylinder; (iii) finally, there is a region of bistability between the in phase and antiphase waves. with each having its own domain of attraction, so that which of the two behaviors occur depends on the condition of initiation of the wave. The results of our computations indicate a qualitative similarity in the behavior of combustion waves in the layer model and in the full 2D model. Specifically, due to the similarity between the small alpha wave behavior in the layer model and the large diameter behavior in the model of waves on the surface of a cylinder, we are able to predict the behavior of the mean velocity for the waves on the cylinder, where computations of the full problem can be rather difficult. We also compute solutions for three or more layers. The results of our computations prompt us to predict that, while planar uniformly propagating waves are unstable, the wave will be quasiplanar, i.e., the resulting spinning waves have very low amplitude hot spots, and travel with the velocity close to that of the uniformly propagating wave. Such waves may be difficult to distinguish from uniformly propagating waves in experiments. We also find that for both the layer model and full 2D problem, steady-state time-dependent waves, e.g., pulsating and spinning wave, have a conserved quantity H which characterizes the excess energy in the wave, just as in the case of uniformly propagating waves. The quantity H, which is generated by dissipation, does not vary in time and is proportional to the diffusivity and caloricity of the system, and inversely proportional to the mean wave velocity.

Aldushin, A. P.↗

Infrared radiation and inversion population of CO2 laser levels in Venusian and Martian atmospheres

Formation mechanisms of nonequilibrium 10 micron CO2 molecule radiation and the possible existence of a natural laser effect in the upper atmospheres of Venus and Mars are theoretically studied. An analysis is made of the excitation process of CO2 molecule vibrational-band levels (with natural isotropic content) induced by direct solar radiation in bands 10.6, 9.4, 4.3, 2.7 and 2.0 microns. The model of partial vibrational-band temperatures was used in the case. The problem of IR radiation transfer in vibrational-rotational bands was solved in the radiation escape approximation.

Gordiyets, B. F.↗

A Polarimetric Scattering Model for the 2-Layer Problem

In this paper, polarimetric signatures from two layers are examined. In the two layer problem a low-loss dielectric layer sits on top of a shallow subsurface layer of scatterers. Volume scattering within the top surface layer is not considered. We show that for polarimetric backscatter data acquired from just one incidence angle, no parameter is proportional to the depth of the subsurface layer. We also show that separation of the two components and inversion of the polarimetric backscatter measurements is not possible, in the latter case due to the mismatch between the number of variables that determine the backscatter behavior, versus the degrees of freedom in the measured values. The interesting results obtained at a sand dune site in SW France courtesy of the group led by Dr. Phillippe Paillou are examined through the prism of this model. It is shown that the model agrees well with results obtained for a subsurface layer at 3 m depth, where the subsurface returns are significant, and that surface scattering alone can account for observations made over a subsurface layer at a depth of 5.8 m, provided the surface slope is adjusted by approx. 10 degrees. The implications for exploration of the near-surface of Mars, a significant proportion of which is obscured by a dust mantle, are discussed.

Freeman, A.↗

Dynamic analysis of astronaut motions in microgravity: Applications for Extravehicular Activity (EVA)

Simulations of astronaut motions during extravehicular activity (EVA) tasks were performed using computational multibody dynamics methods. The application of computational dynamic simulation to EVA was prompted by the realization that physical microgravity simulators have inherent limitations: viscosity in neutral buoyancy tanks; friction in air bearing floors; short duration for parabolic aircraft; and inertia and friction in suspension mechanisms. These limitations can mask critical dynamic effects that later cause problems during actual EVA's performed in space. Methods of formulating dynamic equations of motion for multibody systems are discussed with emphasis on Kane's method, which forms the basis of the simulations presented herein. Formulation of the equations of motion for a two degree of freedom arm is presented as an explicit example. The four basic steps in creating the computational simulations were: system description, in which the geometry, mass properties, and interconnection of system bodies are input to the computer; equation formulation based on the system description; inverse kinematics, in which the angles, velocities, and accelerations of joints are calculated for prescribed motion of the endpoint (hand) of the arm; and inverse dynamics, in which joint torques are calculated for a prescribed motion. A graphical animation and data plotting program, EVADS (EVA Dynamics Simulation), was developed and used to analyze the results of the simulations that were performed on a Silicon Graphics Indigo2 computer. EVA tasks involving manipulation of the Spartan 204 free flying astronomy payload, as performed during Space Shuttle mission STS-63 (February 1995), served as the subject for two dynamic simulations. An EVA crewmember was modeled as a seven segment system with an eighth segment representing the massive payload attached to the hand. For both simulations, the initial configuration of the lower body (trunk, upper leg, and lower leg) was a neutral microgravity posture. In the first simulation, the payload was manipulated around a circular trajectory of 0.15 m radius in 10 seconds. It was found that the wrist joint theoretically exceeded its ulnal deviation limit by as much as 49. 8 deg and was required to exert torques as high as 26 N-m to accomplish the task, well in excess of the wrist physiological limit of 12 N-m. The largest torque in the first simulation, 52 N-m, occurred in the ankle joint. To avoid these problems, the second simulation placed the arm in a more comfortable initial position and the radius and speed of the circular trajectory were reduced by half. As a result, the joint angles and torques were reduced to values well within their physiological limits. In particular, the maximum wrist torque for the second simulation was only 3 N-m and the maximum ankle torque was only 6 N-m.

Newman, Dava J.↗

Optimal mistuning for enhanced aeroelastic stability of transonic fans

An inverse design procedure was developed for the design of a mistuned rotor. The design requirements are that the stability margin of the eigenvalues of the aeroelastic system be greater than or equal to some minimum stability margin, and that the mass added to each blade be positive. The objective was to achieve these requirements with a minimal amount of mistuning. Hence, the problem was posed as a constrained optimization problem. The constrained minimization problem was solved by the technique of mathematical programming via augmented Lagrangians. The unconstrained minimization phase of this technique was solved by the variable metric method. The bladed disk was modelled as being composed of a rigid disk mounted on a rigid shaft. Each of the blades were modelled with a single tosional degree of freedom.

Hall, K. C.↗

A numerical method for solving the equations of compressible viscous flow

Although much progress has already been made in solving problems in aerodynamic design, many new developments are still needed before the equations for unsteady compressible viscous flow can be solved routinely. This paper describes one such development. A method for solving these equations has been devised that: (1) is second-order accurate in space and time; (2) is unconditionally stable; (3) preserves conservation form; (4) requires no block or scalar tridiagonal inversions; (5) is simple and straightforward to program (estimated 10% modification for the update of many existing programs); (6) is more efficient than present methods; and (7) should easily adapt to current and future computer architectures. Computational results for laminar and turbulent flows at Reynolds numbers from 300,000 to 3 x 10 to the 7th and at CFL numbers as high as 1000 are compared with theory and experiment.

Maccormack, R. W.↗

Algorithms and sensitivity analyses for Stratospheric Aerosol and Gas Experiment II water vapor retrieval

The operational inversion algorithm used for the retrieval of the water-vapor vertical profiles from the Stratospheric Aerosol and Gas Experiment II (SAGE II) occultation data is presented. Unlike the algorithm used for the retrieval of aerosol, O3, and NO2, the water-vapor retrieval algorithm accounts for the nonlinear relationship between the concentration versus the broad-band absorption characteristics of water vapor. Problems related to the accuracy of the computational scheme, the accuracy of the removal of other interfering species, and the expected uncertainty of the retrieved profile are examined. Results are presented on the error analysis of the SAGE II water vapor retrieval, indicating that the SAGE II instrument produced good quality water vapor data.

Chu, W. P.↗

A PDE Sensitivity Equation Method for Optimal Aerodynamic Design

The use of gradient based optimization algorithms in inverse design is well established as a practical approach to aerodynamic design. A typical procedure uses a simulation scheme to evaluate the objective function (from the approximate states) and its gradient, then passes this information to an optimization algorithm. Once the simulation scheme (CFD flow solver) has been selected and used to provide approximate function evaluations, there are several possible approaches to the problem of computing gradients. One popular method is to differentiate the simulation scheme and compute design sensitivities that are then used to obtain gradients. Although this black-box approach has many advantages in shape optimization problems, one must compute mesh sensitivities in order to compute the design sensitivity. In this paper, we present an alternative approach using the PDE sensitivity equation to develop algorithms for computing gradients. This approach has the advantage that mesh sensitivities need not be computed. Moreover, when it is possible to use the CFD scheme for both the forward problem and the sensitivity equation, then there are computational advantages. An apparent disadvantage of this approach is that it does not always produce consistent derivatives. However, for a proper combination of discretization schemes, one can show asymptotic consistency under mesh refinement, which is often sufficient to guarantee convergence of the optimal design algorithm. In particular, we show that when asymptotically consistent schemes are combined with a trust-region optimization algorithm, the resulting optimal design method converges. We denote this approach as the sensitivity equation method. The sensitivity equation method is presented, convergence results are given and the approach is illustrated on two optimal design problems involving shocks.

Borggaard, Jeff↗

An Efficient Spectral Method for Ordinary Differential Equations with Rational Function Coefficients

We present some relations that allow the efficient approximate inversion of linear differential operators with rational function coefficients. We employ expansions in terms of a large class of orthogonal polynomial families, including all the classical orthogonal polynomials. These families obey a simple three-term recurrence relation for differentiation, which implies that on an appropriately restricted domain the differentiation operator has a unique banded inverse. The inverse is an integration operator for the family, and it is simply the tridiagonal coefficient matrix for the recurrence. Since in these families convolution operators (i.e. matrix representations of multiplication by a function) are banded for polynomials, we are able to obtain a banded representation for linear differential operators with rational coefficients. This leads to a method of solution of initial or boundary value problems that, besides having an operation count that scales linearly with the order of truncation N, is computationally well conditioned. Among the applications considered is the use of rational maps for the resolution of sharp interior layers.

Coutsias, Evangelos A.↗

A Numerical Method for Solving the Equations of Compressible Viscous Flow

Although much progress has already been made In solving problems in aerodynamic design, many new developments are still needed before the equations for unsteady compressible viscous flow can be solved routinely. This paper describes one such development. A new method for solving these equations has been devised that 1) is second-order accurate in space and time, 2) is unconditionally stable, 3) preserves conservation form, 4) requires no block or scalar tridiagonal inversions, 5) is simple and straightforward to program (estimated 10% modification for the update of many existing programs), 6) is more efficient than present methods, and 7) should easily adapt to current and future computer architectures. Computational results for laminar and turbulent flows at Reynolds numbers from 3 x 10(exp 5) to 3 x 10(exp 7) and at CFL numbers as high as 10(exp 3) are compared with theory and experiment.

MacCormack, R. W.↗

6 centimeter observations of solar active regions with 6 sec resolution

We have used the Westerbork Synthesis Radio Telescope for observations of solar active regions at 6 cm during 1974 May 8-10. The particular problems of using a sidereal instrument for solar observations are discussed. Maps of total intensity and circular polarization were obtained for four active regions with sunspots and a plage region; the linear polarization was below the noise limit. The brightest components of the 6 cm emission are associated with sunspots, while weaker components are associated with the plage magnetic fields, both longitudinal and transverse. An inversion of the sense of circular polarization was observed in the preceding part of a region with bipolar magnetic-field structure, about two days after the region passed through the central meridian. Nine bursts were also observed, and their one-dimensional positions with respect to the associated active regions are given.

Kundu, M. R.↗

Inverting a dispersive scene's side-scanned image

Consideration is given to the problem of using a remotely sensed, side-scanned image of a time-variant scene, which changes according to a dispersion relation, to estimate the structure at a given moment. Additive thermal noise is neglected in the models considered in the formal treatment. It is shown that the dispersion relation is normalized by the scanning velocity, as is the group scanning velocity component. An inversion operation is defined for noise-free images generated by SAR. The method is extended to the inversion of noisy imagery, and a formulation is defined for spectral density estimation. Finally, the methods for a radar system are used for the case of sonar.

Harger, R. O.↗

The dynamics of orbital maneuvering: Design and evaluation of a visual display aid for human controllers

An interactive proximity operations planning system, which allows on-site planning of fuel-efficient, multi-burn maneuvers in a potential multi-spacecraft environment was developed and tested. Though this display system most directly assists planning by providing visual feedback to aid visualization of the trajectories and constraints, its most significant features include an inverse dynamics algorithm that removes control nonlinearities facing the operator and a trajectory planning technique that reduces the order of control and creates, through a geometric spread-sheet the illusion of an inertially stable environment. This synthetic environment provides the user with control of relevant static and dynamic properties of way-points during small orbital changes allowing independent solutions to the normally coupled problems of orbital maneuvering. An experiment was carried out in which experienced operators were required to plan a trajectory to retrieve an object accidently separated from a dual-keel space station. The time required to plan these maneuvers was found to be predicted by the direction of the insertion thrust and did not depend on the point of separation from the space station.

Ellis, Stephen R.↗

The isolation limits of stochastic vibration

The vibration isolation problem is formulated as a 1D kinematic problem. The geometry of the stochastic wall trajectories arising from the stroke constraint is defined in terms of their significant extrema. An optimal control solution for the minimum acceleration return path determines a lower bound on platform mean square acceleration. This bound is expressed in terms of the probability density function on the significant maxima and the conditional fourth moment of the first passage time inverse. The first of these is found analytically while the second is found using a Monte Carlo simulation. The rms acceleration lower bound as a function of available space is then determined through numerical quadrature.

Knopse, C. R.↗

Analytical modeling of structure-soil systems for lunar bases

The study of the behavior of granular materials in a reduced gravity environment and under low effective stresses became a subject of great interest in the mid 1960's when NASA's Surveyor missions to the Moon began the first extraterrestrial investigation and it was found that Lunar soils exhibited properties quite unlike those on Earth. This subject gained interest during the years of the Apollo missions and more recently due to NASA's plans for future exploration and colonization of Moon and Mars. It has since been clear that a good understanding of the mechanical properties of granular materials under reduced gravity and at low effective stress levels is of paramount importance for the design and construction of surface and buried structures on these bodies. In order to achieve such an understanding it is desirable to develop a set of constitutive equations that describes the response of such materials as they are subjected to tractions and displacements. This presentation examines issues associated with conducting experiments on highly nonlinear granular materials under high and low effective stresses. The friction and dilatancy properties which affect the behavior of granular soils with low cohesion values are assessed. In order to simulate the highly nonlinear strength and stress-strain behavior of soils at low as well as high effective stresses, a versatile isotropic, pressure sensitive, third stress invariant dependent, cone-cap elasto-plastic constitutive model was proposed. The integration of the constitutive relations is performed via a fully implicit Backward Euler technique known as the Closest Point Projection Method. The model was implemented into a finite element code in order to study nonlinear boundary value problems associated with homogeneous as well as nonhomogeneous deformations at low as well as high effective stresses. The effect of gravity (self-weight) on the stress-strain-strength response of these materials is evaluated. The calibration of the model is performed via three techniques: (1) physical identification, (2) optimized calibration at the constitutive level, and (3) optimized calibration at the finite element level (Inverse Identification). Activities are summarized in graphic and outline form.

Macari-Pasqualino, Jose Emir↗