Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear equations”

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,207 records · Page 67

The RAPT control method for nonlinear systems

The Rapid Advancement Preview Tracking (RAPT) control method uses both a nonlinear system model and a linearized system model to provide a form of feedback/preview control for nonlinear systems. It uses measurements of the states to calculate the current preview control signal from a local linearization of the state equations, and then iteratively calculates a perturbation correction for that control using the nonlinear equations. The corrected control signal is then applied to the system, driving the end-effector to the next position along a specified trajectory. Computer simulations of tracking control for two-link planar manipulators indicate that this algorithm may be fast enough for real-time control.

Long, T. W.↗

Nonlinear interaction between a pair of oblique modes in a supersonic mixing layer: Long-wave limit

The nonlinear interaction between a pair of symmetric, oblique, and spatial instability modes is studied in the long-wave limit using asymptotic methods. The base flow is taken to be a supersonic mixing layer whose Mach number is such that the corresponding vortex sheet is marginally stable according to Miles' criterion. It is shown that the amplitude of the mode obeys a nonlinear integro-differential equation. Numerical solutions of this equation show that, when the obliqueness angle is less than pi/4, the effect of the nonlinearity is to enhance the growth rate of the instability. The solution terminates in a singularity at a finite streamwise location. This result is reminiscent of that obtained in the vicinity of the neutral point by other authors in several different types of flows. On the other hand, when the obliqueness angle is more than pi/4, the streamwise development of the amplitude is characterized by a series of modulations. This arises from the fact that the nonlinear term in the amplitude equation may be either stabilizing or destabilizing, depending on the value of the streamwise coordinate. However, even in this case the amplitude of the disturbance increases, though not as rapidly as in the case for which the angle is less than pi/4. Quite generally then, the nonlinear interaction between two oblique modes in a supersonic mixing layer enhances the growth of the disturbance.

Balsa, Thomas F.↗

An Exponential Finite Difference Technique for Solving Partial Differential Equations

An exponential finite difference algorithm, as first presented by Bhattacharya for one-dimensianal steady-state, heat conduction in Cartesian coordinates, has been extended. The finite difference algorithm developed was used to solve the diffusion equation in one-dimensional cylindrical coordinates and applied to two- and three-dimensional problems in Cartesian coordinates. The method was also used to solve nonlinear partial differential equations in one (Burger's equation) and two (Boundary Layer equations) dimensional Cartesian coordinates. Predicted results were compared to exact solutions where available, or to results obtained by other numerical methods. It was found that the exponential finite difference method produced results that were more accurate than those obtained by other numerical methods, especially during the initial transient portion of the solution. Other applications made using the exponential finite difference technique included unsteady one-dimensional heat transfer with temperature varying thermal conductivity and the development of the temperature field in a laminar Couette flow.

Handschuh, Robert F.↗

Annular suspension and pointing system with controlled DC electromagnets

The Annular Suspension and Pointing System (ASPS) developed by the Flight System division of Sperry Corporation is a six-degree of freedom payload pointing system designed for use with the space shuttle. This magnetic suspension and pointing system provides precise controlled pointing in six-degrees of freedom, isolation of payload-carrier disturbances, and end mount controlled pointing. Those are great advantages over the traditional mechanical joints for space applications. In this design, we first analyzed the assumed model of the single degree ASPS bearing actuator and obtained the plant dynamics equations. By linearizing the plant dynamics equations, we designed the cascade and feedback compensators such that a stable and satisfied result was obtained. The specified feedback compensator was computer simulated with the nonlinearized plant dynamics equations. The results indicated that an unstable output occurred. In other words, the designed feedback compensator failed. The failure of the design is due to the Taylor's series expansion not converging.

Vu, Josephine Lynn↗

Magnetic and thermodynamic structure of the solar corona during sunspot minimum

A model of the solar corona during sunspot minimum was constructed. It was suggested that pole-equator asymmetry is caused by interaction of global azimuthal electric currents in the corona and potential magnetic fields originating below the corona. The problem was reduced to a nonlinear ordinary differential equation for a structural function. Taking the difference in electron density between the equatorial plane and the polar direction from observations, the basic equation was solved and the global electric current distribution around the Sun along with the magnetic and thermodynamic structures of the solar corona was derived. The asymmetry between the two hemispheres, which can be due to the magnetic quadrupole, is considered in detail. It is shown that an additional quadrupole term in the magnetic field representation affects the position of polar coronal hole boundaries differently in the two hemispheres. As a result, the outflow from the two polar regions is found to be different.

Osherovich, V. A.↗

Gust response analysis for cascades operating in nonuniform mean flows

The linearized unsteady aerodynamic response of a cascade of airfoils subjected to entropic, vortical, and acoustic gusts is analyzed. Field equations for the first-order unsteady perturbation flow are obtained by linearizing the full time-dependent mass, momentum, and energy conservation equations about a nonlinear, isentropic, and irrotational mean or steady flow. A splitting technique is then used to decompose the unsteady velocity field into irrotational and rotational parts leading to field equations for the unsteady entropy, rotational velocity, and irrotational velocity fluctuations that are coupled only sequentially. The entropic and rotational velocity fluctuations can be described in terms of the mean-flow drift and stream functions which can be computed numerically. The irrotational unsteady velocity is described by an inhomogeneous linearized potential equation which contains a source term that depends on the rotational velocity field. This equation is solved via a finite difference technique. Results are presented to indicate the status of the numerical solution procedure and to demonstrate the impact of blade geometry and mean blade loading on the aerodynamic response of cascades to vortical gust excitations. The analysis described leads to very efficient predictions of cascade unsteady aerodynamic phenomena making it useful for turbomachinery aeroelastic and aeroacoustic design applications.

Hall, Kenneth C.↗

Gust response analysis for cascades operating in nonuniform mean flows

The linearized unsteady aerodynamic response of a cascade of airfoils subjected to entropic, vortical, and acoustic gusts is analyzed. Field equations for the first-order unsteady perturbation flow are obtained by linearizing the full time-dependent mass, momentum, and energy conservation equations about a nonlinear, isentropic, and irrotational mean or steady flow. A splitting technique is then used to decompose the unsteady velocity field into irrotational and rotational parts leading to field equations for the unsteady entropy, rotational velocity, and irrotational velocity fluctuations that are coupled only sequentially. The entropic and rotational velocity fluctuations can be described in terms of the mean-flow drift and stream functions which can be computed numerically. The irrotational unsteady velocity is described by an inhomogeneous linearized potential equation which contains a source term that depends on the rotational velocity field. This equation is solved via a finite difference technique. Results are presented to indicate the status of the numerical solution procedure and to demonstrate the impact of blade geometry and mean blade loading on the aerodynamic response of cascades to vortical gust excitations. The analysis described leads to very efficient predictions of cascade unsteady aerodynamics phenomena making it useful for turbomachinery aeroelastic and aeroacoustic design applications.

Hall, Kenneth C.↗

Nonlinear flap-lag-extensional vibrations of rotating, pretwisted, preconed beams including Coriolis effects

The effects of pretwist, precone, setting angle, Coriolis forces and second degree geometric nonlinearities on the natural frequencies, steady state deflections and mode shapes of rotating, torsionally rigid, cantilevered beams were studied. The governing coupled equations of flap lag extensional motion are derived including the effects of large precone and retaining geometric nonlinearities up to second degree. The Galerkin method, with nonrotating normal modes, is used for the solution of both steady state nonlinear equations and linear perturbation equations. Parametric indicating the individual and collective effects of pretwist, precone, Coriolis forces and second degree geometric nonlinearities on the steady state deflection, natural frequencies and mode shapes of rotating blades are presented. It is indicated that the second degree geometric nonlinear terms, which vanish for zero precone, can produce frequency changes of engineering significance. Further confirmation of the validity of including those generated by MSC NASTRAN. It is indicated that the linear and nonlinear Coriolis effects must be included in analyzing thick blades. The Coriolis effects are significant on the first flatwise and the first edgewise modes.

Subrahmanyam, K. B.↗

Second-order non-iterative ADI solution of non-linear partial differential equations

A new method for the solution of non-linear partial differential equations by an ADI procedure is described. Although the method is second order accurate in time, it does not require either iterations or predictor corrector methods to overcome the nonlinearity of the equations. Thus the computational effort required for the solution of the non-linear problem becomes similar to that required for the linear case. The method is applied to a two-dimensional 'extended Burgers equation'. Linear stability is studied, and some numerical solutions obtained. The improved accuracy obtained by the 2nd order truncation error is clearly manifested.

Wolfshtein, M.↗

Analysis of turbulent free-jet hydrogen-air diffusion flames with finite chemical reaction rates

A numerical analysis is presented of the nonequilibrium flow field resulting from the turbulent mixing and combustion of an axisymmetric hydrogen jet in a supersonic parallel ambient air stream. The effective turbulent transport properties are determined by means of a two-equation model of turbulence. The finite-rate chemistry model considers eight elementary reactions among six chemical species: H, O, H2O, OH, O2 and H2. The governing set of nonlinear partial differential equations was solved by using an implicit finite-difference procedure. Radial distributions were obtained at two downstream locations for some important variables affecting the flow development, such as the turbulent kinetic energy and its dissipation rate. The results show that these variables attain their peak values on the axis of symmetry. The computed distribution of velocity, temperature, and mass fractions of the chemical species gives a complete description of the flow field. The numerical predictions were compared with two sets of experimental data. Good qualitative agreement was obtained.

Sislian, J. P.↗

Control of nonlinear systems with applications to constrained robots and spacecraft attitude stabilization

This thesis is organized in two parts. In Part 1, control systems described by a class of nonlinear differential and algebraic equations are introduced. A procedure for local stabilization based on a local state realization is developed. An alternative approach to local stabilization is developed based on a classical linearization of the nonlinear differential-algebraic equations. A theoretical framework is established for solving a tracking problem associated with the differential-algebraic system. First, a simple procedure is developed for the design of a feedback control law which ensures, at least locally, that the tracking error in the closed loop system lies within any given bound if the reference inputs are sufficiently slowly varying. Next, by imposing additional assumptions, a procedure is developed for the design of a feedback control law which ensures that the tracking error in the closed loop system approaches zero exponentially for reference inputs which are not necessarily slowly varying. The control design methodologies are used for simultaneous force and position control in constrained robot systems. The differential-algebraic equations are shown to characterize the slow dynamics of a certain nonlinear control system in nonstandard singularly perturbed form. In Part 2, the attitude stabilization (reorientation) of a rigid spacecraft using only two control torques is considered. First, the case of momentum wheel actuators is considered. The complete spacecraft dynamics are not controllable. However, the spacecraft dynamics are small time locally controllable in a reduced sense. The reduced spacecraft dynamics cannot be asymptotically stabilized using continuous feedback, but a discontinuous feedback control strategy is constructed. Next, the case of gas jet actuators is considered. If the uncontrolled principal axis is not an axis of symmetry, the complete spacecraft dynamics are small time locally controllable. However, the spacecraft attitude cannot be asymptotically stabilized using continuous feedback, but a discontinuous stabilizing feedback control strategy is constructed. If the uncontrolled principal axis is an axis of symmetry, the complete spacecraft dynamics cannot be stabilized. However, the spacecraft dynamics are small time locally controllable in a reduced sense. The reduced spacecraft dynamics cannot be asymptotically stabilized using continuous feedback, but again a discontinuous feedback control strategy is constructed.

Krishnan, Hariharan↗

The role of shock induced trailing-edge separation in limit cycle oscillations

The potential role of shock induced trailing edge separation (SITES) in limit cycle oscillations (LCO) was established. It was shown that the flip-flop characteristics of transition to and from SITES as well as its hysteresis could couple with wing modes with torsional motion and low damping. This connection led to the formulation of a very simple nonlinear math model using the linear equations of motion with a nonlinear step forcing function with hysteresis. A finite difference solution with time was developed and calculations were made for the F-111 TACT were used to determine the step forcing function due to SITES transition. Since no data were available for the hysteresis, a parameter study was conducted allowing the hysteresis effect to vary. Very small hysteresis effects, which were within expected bounds, were required to obtain reasonable response levels that essentially agreed with flight test results. Also in agreement with wind tunnel tests, LCO calculations for the 1/6 scale F-111 model showed that the model should have not experienced LCO.

Cunningham, Atlee M., Jr.↗

The effects of molecular weight on the single lap shear creep and constant strain rate behavior of thermoplastic polyimidesulfone adhesive

The bonded shear creep and constant strain rate behaviors of zero, one, and three percent endcapped thermoplastic polyimidesulfone adhesive were examined at room and elevated temperatures. Endcapping was accomplished by the addition of phthalic anhydrides. The primary objective was to determine the effects of molecular weight on the mechanical properties of the adhesive. Viscoelastic and nonlinear elastic constitutive equations were utilized to model the adhesive. Ludwik's and Crochet's relations were used to describe the experimental failure data. The effects of molecular weight changes on the above mentioned mechanical behavior were assessed. The viscoelastic Chase-Goldsmith and elastic nonlinear relations gave a good fit to the experimental stress strain behavior. Crochet's relations based on Maxwell and Chase-Goldsmith models were fit to delayed failure data. Ludwik's equations revealed negligible rate dependence. Ultimate stress levels and the safe levels for creep stresses were found to decrease as molecular weight was reduced.

Dembosky, Stanley K.↗

Theory of wing rock

Wing rock is one type of lateral-directional instabilities at high angles of attack. To predict wing rock characteristics and to design airplanes to avoid wing rock, parameters affecting wing rock characteristics must be known. A new nonlinear aerodynamic model is developed to investigate the main aerodynamic nonlinearities causing wing rock. In the present theory, the Beecham-Titchener asymptotic method is used to derive expressions for the limit-cycle amplitude and frequency of wing rock from nonlinear flight dynamics equations. The resulting expressions are capable of explaining the existence of wing rock for all types of aircraft. Wing rock is developed by negative or weakly positive roll damping, and sustained by nonlinear aerodynamic roll damping. Good agreement between theoretical and experimental results is obtained.

Hsu, C.-H.↗

Vectorizable implicit algorithms for the flux-difference split, three-dimensional Navier-Stokes equations

The computational efficiency of four vectorizable implicit algorithms is assessed when applied to calculate steady-state solutions to the three-dimensional, incompressible Navier-Stokes equations in general coordinates. Two of these algorithms are characterized as hybrid schemes; that is, they combine some approximate factorization in two coordinate directions with relaxation in the remaining spatial direction. The other two algorithms utilize an approximate factorization approach which yields two-factor algorithms for three-dimensional systems. All four algorithms are implemented in identical high-resolution upwind schemes for the flux-difference split Navier-Stokes equations. These highly nonlinear schemes are obtained by extending an implicit Total Variation Diminishing (TVD) scheme recently developed for linear one-dimensional systems of hyperbolic conservation laws to the three-dimensional Navier-Stokes equations. The computation of vortical flow over a sharp-edged, thin delta wing has been chosen as a common numerical test case. The convergence of the algorithms is discussed and the accuracy of the computed flow-field results is assessed. The validity of the present results are demonstrated by a comparison with experimental data.

Hartwich, P. M.↗

Nonlinear Oscillators in Space Physics

We discuss dynamical systems that produce an oscillation without an external time dependent source. Numerical results are presented for nonlinear oscillators in the Em1h's atmosphere, foremost the quasi-biennial oscillation (QBOl. These fluid dynamical oscillators, like the solar dynamo, have in common that one of the variables in a governing equation is strongly nonlinear and that the nonlinearity, to first order, has particular form. of 3rd or odd power. It is shown that this form of nonlinearity can produce the fundamental li'equency of the internal oscillation. which has a period that is favored by the dynamical condition of the fluid. The fundamental frequency maintains the oscillation, with no energy input to the system at that particular frequency. Nonlinearities of 2nd or even power could not maintain the oscillation.

Lester,Daniel↗

Secondary Instability of Second Modes in Hypersonic Boundary Layers

Second mode disturbances dominate the primary instability stage of transition in a number of hypersonic flow configurations. The highest amplification rates of second mode disturbances are usually associated with 2D (or axisymmetric) perturbations and, therefore, a likely scenario for the onset of the three-dimensionality required for laminar-turbulent transition corresponds to the parametric amplification of 3D secondary instabilities in the presence of 2D, finite amplitude second mode disturbances. The secondary instability of second mode disturbances is studied for selected canonical flow configurations. The basic state for the secondary instability analysis is obtained by tracking the linear and nonlinear evolution of 2D, second mode disturbances using nonlinear parabolized stability equations. Unlike in previous studies, the selection of primary disturbances used for the secondary instability analysis was based on their potential relevance to transition in a low disturbance environment and the effects of nonlinearity on the evolution of primary disturbances was accounted for. Strongly nonlinear effects related to the self-interaction of second mode disturbances lead to an upstream shift in the upper branch neutral location. Secondary instability computations confirm the previously known dominance of subharmonic modes at relatively small primary amplitudes. However, for the Purdue Mach 6 compression cone configuration, it was shown that a strong fundamental secondary instability can exist for a range of initial amplitudes of the most amplified second mode disturbance, indicating that the exclusive focus on subharmonic modes in the previous applications of secondary instability theory to second mode primary instability may not have been fully justified.

Li, Fei↗