SEARCH · Engineering Papers
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.
A new formulation of the one-dimensional Vlasov equation is derived which is analogous to the Kalman-transformed cold-plasma equations. The equations are shown to yield nonsecular, nonlinear approximations to a source or boundary-value problem. It is suggested that the formulation may have other applications in nonlinear plasma theory.
The H(sub z)-polarized inverse problem has rarely been studied before due to the complicated way in which the unknown permittivity appears in the wave equation. This problem is equivalent to the acoustic inverse problem with variable density. We have recently reported the solution to the nonlinear variable-permittivity H(sub z)-polarized inverse problem using the Born iterative method. Here, the linear inverse problem is solved for permittivity (epsilon) and permeability (mu) using a different approach which is an extension of the basic ideas of diffraction tomography (DT). The key to solving this problem is to utilize frequency diversity to obtain the required independent measurements. The receivers are assumed to be in the far field of the object, and plane wave incidence is also assumed. It is assumed that the scatterer is weak, so that the Born approximation can be used to arrive at a relationship between the measured pressure field and two terms related to the spatial Fourier transform of the two unknowns, epsilon and mu. The term involving permeability corresponds to monopole scattering and that for permittivity to dipole scattering. Measurements at several frequencies are used and a least squares problem is solved to reconstruct epsilon and mu. It is observed that the low spatial frequencies in the spectra of epsilon and mu produce inaccuracies in the results. Hence, a regularization method is devised to remove this problem. Several results are shown. Low contrast objects for which the above analysis holds are used to show that good reconstructions are obtained for both permittivity and permeability after regularization is applied.
In this paper transient fields for antennas with more general geometries are calculated directly using Finite Difference Time Domain methods. In each FDTD cell which contains a nonlinear load, a nonlinear equation is solved at each time step. As a test case the transient current in a long dipole antenna with a nonlinear load excited by a pulsed plane wave is computed using this approach. The results agree well with both calculated and measured results previously published. The approach given here extends the applicability of the FDTD method to problems involving scattering from targets including nonlinear loads and materials, and to coupling between antennas containing nonlinear loads. It may also be extended to propagation through nonlinear materials.
The proposed method for the treatment of two dimensional solidification problems is based on quasilinearization of the transformed heat conduction equation and latent heat condition at the interface and an iterative sequence in which these are solved simultaneously. Modern algorithms for solving such sparse systems mean that most of the storage advantage of other methods are reduced and the speed of solution can be improved.
The determination of ignition or thermal explosion in an oxidizing porous body of material, as described by a dimensionless reaction-diffusion equation of the form .tu = .2u + .e-1/u over the bounded region O, is critically reexamined from a modern perspective using numerical methodologies. First, the classic stationary model is revisited to establish the proper reference frame for the steady-state solution space, and it is demonstrated how the resulting nonlinear two-point boundary value problem can be reexpressed as an initial value problem for a system of first-order differential equations, which may be readily solved using standard algorithms. Then, the numerical procedure is implemented and thoroughly validated against previous computational results based on sophisticated path-following techniques. Next, the transient nonstationary model is attacked, and the full nonlinear form of the reaction-diffusion equation, including a generalized convective boundary condition, is discretized and expressed as a system of linear algebraic equations. The numerical methodology is implemented as a computer algorithm, and validation computations are carried out as a prelude to a broad-ranging evaluation of the assembly problem and identification of the watershed critical initial temperature conditions for thermal ignition. This numerical methodology is then used as the basis for studying the relationship between the shape of the critical initial temperature distribution and the corresponding spatial moments of its energy content integral and an attempt to forge a fundamental conjecture governing this relation. Finally, the effects of dynamic boundary conditions on the classic storage problem are investigated and the groundwork is laid for the development of an approximate solution methodology based on adaptation of the standard stationary model.
A simple, two dimensional, incompressible and inviscid model for the problem posed by a two dimensional wing with a jet issuing from its lower surface is considered and a parametric analysis is carried out to observe how the aerodynamic characteristics depend on the different parameters. The mathematical problem constitutes a boundary value problem where the position of part of the boundary is not known a priori. A nonlinear optimization approach was used to solve the problem, and the analysis reveals interesting characteristics that may help to better understand the physics involved in more complex situations in connection with high lift systems.
This paper develops numerical strategies which enable the constrained incremental Newton-Raphson scheme to handle the static response of structure to loading fields with completely generalized histories. This is made possible through the use of specially warped hyperelliptic constraint surfaces which control successive or clustered load steps in the vicinity of loading events with specific timing schedules. Such an approach enables improved convergence and stability characteristics. Due to the generality of the methodology, pre- and postbuckling behavior caused by both kinematic and material nonlinearity can be handled. To demonstrate the scheme, the results of several bench-mark problems are also presented. These include situations involving nonlinear kinematics as well as highly history-dependent elastic-plastic and thermoelastic-plastic material behavior.
Hertzian contact theory is applied to a butt joint with specially mismatched bearing surfaces to devise a simple mathematical model of nonlinear axial force-displacement behavior in jointed members. Normalized tangent stiffness-force plots, for several values of a joint imperfection parameter, are presented for the sample case of solid structural members of circular section. The results illustrate the potential problem of high joint compliance at low axial-force levels, as well as the generally desirable stiffening and linearizing effects of preload. A nonlinear oscillator problem based on the static model is also formulated and solved to illustrate the effect of amplitude on natural frequency. As expected, natural frequency is low when amplitude is small. The results call attention to the important roles that tight tolerances and preload are expected to play in the design and fabrication of deployable and erectable truss-type space structures.
In aerospace applications, magnetic suspension systems may be required to operate over large variations in air-gap. Thus the nonlinearities inherent in most types of suspensions have a significant effect. Specifically, large variations in operating point may make it difficult to design a linear controller which gives satisfactory stability and performance over a large range of operating points. One way to address this problem is through the use of nonlinear compensation techniques such as feedback linearization. Nonlinear compensators have received limited attention in the magnetic suspension literature. In recent years, progress has been made in the theory of nonlinear control systems, and in the sub-area of feedback linearization. The idea is demonstrated of feedback linearization using a second order suspension system. In the context of the second order suspension, sampling rate issues in the implementation of feedback linearization are examined through simulation.
Development was extended of a simplified inelastic analysis computer program (ANSYMP) for predicting the stress-strain history at the critical location of a thermomechanically cycled structure from an elastic solution. The program uses an iterative and incremental procedure to estimate the plastic strains from the material stress-strain properties and a plasticity hardening model. Creep effects can be calculated on the basis of stress relaxation at constant strain, creep at constant stress, or a combination of stress relaxation and creep accumulation. The simplified method was exercised on a number of problems involving uniaxial and multiaxial loading, isothermal and nonisothermal conditions, dwell times at various points in the cycles, different materials, and kinematic hardening. Good agreement was found between these analytical results and nonlinear finite-element solutions for these problems. The simplified analysis program used less than 1 percent of the CPU time required for a nonlinear finite-element analysis.
The paper presents a two-stage approach for designing optimal reconfiguration maneuvers for multiple spacecraft. These maneuvers involve well-coordinated and highly-coupled motions of the entire fleet of spacecraft while satisfying an arbitrary number of constraints. This problem is particularly difficult because of the nonlinearity of the attitude dynamics, the non-convexity of some of the constraints, and the coupling between the positions and attitudes of all spacecraft. As a result, the trajectory design must be solved as a single 6N DOF problem instead of N separate 6 DOF problems. The first stage of the solution approach quickly provides a feasible initial solution by solving a simplified version without differential constraints using a bi-directional Rapidly-exploring Random Tree (RRT) planner. A transition algorithm then augments this guess with feasible dynamics that are propagated from the beginning to the end of the trajectory. The resulting output is a feasible initial guess to the complete optimal control problem that is discretized in the second stage using a Gauss pseudospectral method (GPM) and solved using an off-the-shelf nonlinear solver. This paper also places emphasis on the importance of the initialization step in pseudospectral methods in order to decrease their computation times and enable the solution of a more complex class of problems. Several examples are presented and discussed.
The use of two Class 2/fixed mesh/fixed order/multistep integration packages of the PECE type for the numerical integration of the second order, nonlinear, ordinary differential equation of the satellite orbit problem. These two methods are referred to as the general and the second sum formulations. The derivation of the basic equations which characterize each formulation and the role of the basic equations in the PECE algorithm are discussed. Possible starting procedures are examined which may be used to supply the initial set of values required by the fixed mesh/multistep integrators. The results of the general and second sum integrators are compared to the results of various fixed step and variable step integrators.
Three techniques for the solution of bilinear estimation problems are presented. First, finite dimensional optimal nonlinear estimators are presented for certain bilinear systems evolving on solvable and nilpotent lie groups. Then the use of harmonic analysis for estimation problems evolving on spheres and other compact manifolds is investigated. Finally, an approximate estimation technique utilizing cumulants is discussed.
The boundary value problem for vortex separation at zero sideslip on cones and tangent ogives is set up by means of a discrete vortex model. The nonlinear algebraic equations for the boundary value problem admit multiple, physically feasible solutions, including the symmetric and asymmetric vortex solutions. Multiple solutions are proposed as an alternative explanation of the existence of asymmetric vortex separation at zero sideslip.
This paper discusses the development of large scale numerical optimal control algorithms for nonlinear systems and their application to finite element models of structures. This work is based on our expansion of the optimal control algorithm (DDP) in the following steps: improvement of convergence for initial policies in non-convex regions, development of a numerically accurate penalty function method approach for constrained DDP problems, and parallel processing on supercomputers. The expanded constrained DDP algorithm was applied to the control of a four-bay, two dimensional truss with 12 soft members, which generates geometric nonlinearities. Using an explicit finite element model to describe the structural system requires 32 state variables and 10,000 time steps. Our numerical results indicate that for constrained or unconstrained structural problems with nonlinear dynamics, the results obtained by our expanded constrained DDP are significantly better than those obtained using linear-quadratic feedback control.
At the present time, there are three general analytical techniques available to study problems of unsteady motions in rocket motors: linear stability analysis; approximate nonlinear analysis, founded on examining the behavior of coupled normal modes; and numerical calculations based on the conservation equations for one-dimensional flows. The last two yield the linear results as a limit. It is the main purpose of this paper to check the accuracy of the approximate analysis against the numerical analysis for some special cases. The results provide some justification for using the approximate analysis to study three-dimensional problems.
Failure mechanisms of materials under very high strains experienced at and ahead of the crack tip such as formation, growth, and interaction of microvoids in ductile materials, microcracks in brittle solids or crazes in polymers and adhesives are represented by one-dimensional, nonlinear stress-strain relations possessing different ways by which the material loses capacity to carry load up to fracture or total separation. A double cantilever beam (DCB) type specimen is considered. The nonlinear material is confined to a thin strip between the two elastic beams loaded by a wedge. The problem is first modeled as a beam on a nonlinear foundation. The pertinent equation is solved numerically as a two-point boundary value problem for both the stationary and the quasi-stationay propagating crack. A finite element model is then used to model the problem in more detail in order to assess the adequacy of the beam model for the reduction of experimental data to determine in-situ properties of the thin interlayer.
A method of using an electronic nose to detect an airborne mixture of known chemical compounds and measure the temporally varying concentrations of the individual compounds is undergoing development. In a typical intended application, the method would be used to monitor the air in an inhabited space (e.g., the interior of a building) for the release of solvents, toxic fumes, and other compounds that are regarded as contaminants. At the present state of development, the method affords a capability for identifying and quantitating one or two compounds that are members of a set of some number (typically of the order of a dozen) known compounds. In principle, the method could be extended to enable monitoring of more than two compounds. An electronic nose consists of an array of sensors, typically made from polymer carbon composites, the electrical resistances of which change upon exposure to a variety of chemicals. By design, each sensor is unique in its responses to these chemicals: some or all of the sensitivities of a given sensor to the various vapors differ from the corresponding sensitivities of other sensors. In general, the responses of the sensors are nonlinear functions of the concentrations of the chemicals. Hence, mathematically, the monitoring problem is to solve the set of time-dependent nonlinear equations for the sensor responses to obtain the time dependent concentrations of individual compounds. In the present developmental method, successive approximations of the solution are generated by a learning algorithm based on independent-component analysis (ICA) an established information theoretic approach for transforming a vector of observed interdependent signals into a set of signals that are as nearly statistically independent as possible.