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 703 records · Page 39

Prediction of unsteady aerodynamic loads in cascades using the linearized Euler equations on deforming grids

A linearized Euler solver for calculating unsteady flows in turbomachinery blade rows due to both incident gusts and blade motion is presented. Using the linearized Euler technique, one decomposes the flow into a mean (or steady) flow plus an unsteady, harmonically varying, small disturbance flow. Linear variable coefficient equations describe the small disturbance behavior of the flow, and are solved using a pseudo-time marching Lax-Wendroff scheme. For the blade motion problem, a harmonically deforming computational rid that conforms to the motion of vibrating blades eliminates large error producing mean flow gradient terms that would otherwise appear in the unsteady flow tangency boundary condition. The paper also presents a new, numerically exact, nonreflecting far-field boundary condition based on an eigenanalysis of the discretized equations. Computed flow solutions demonstrate the computational accuracy and efficiency of the present method. The solution of the linearized Euler equations requires one to two orders of magnitude less computer time than solution of the nonlinear Euler equations using traditional time-accurate time-marching techniques. In addition, the deformable grid significantly improves the accuracy of the solution.

Hall, Kenneth C.↗

Unified nonlinear analysis for nonhomogeneous anisotropic beams with closed cross sections

A unified methodology for geometrically nonlinear analysis of nonhomogeneous, anisotropic beams is presented. A 2D cross-sectional analysis and a nonlinear 1D global deformation analysis are derived from the common framework of a 3D, geometrically nonlinear theory of elasticity. The only restrictions are that the strain and local rotation are small compared to unity and that warping displacements are small relative to the cross-sectional dimensions. It is concluded that the warping solutions can be affected by large deformation and that this could alter the incremental stiffnes of the section. It is shown that sectional constants derived from the published, linear analysis can be used in the present nonlinear, 1D analysis governing the global deformation of the beam, which is based on intrinsic equations for nonlinear beam behavior. Excellent correlation is obtained with published experimental results for both isotropic and anisotropic beams undergoing large deflections.

Atilgan, Ali R.↗

Vibrational Control of a Nonlinear Elastic Panel

The paper is concerned with the stabilization of the nonlinear panel oscillation by an active control. The control is actuated by a combination of additive and parametric vibrational forces. A general method of vibrational control is presented for stabilizing panel vibration satisfying a nonlinear beam equation. To obtain analytical results, a perturbation technique is used in the case of weak nonlinearity. Possible application to other types of problems is briefly discussed.

Chow, P. L.↗

Computation of a controlled store separation from a cavity

Coupling of the Reynolds-averaged Navier-Stokes equations, rigid-body dynamics, and a pitch attitude control law is demonstrated in two- and three-dimensions. The application problem was the separation of a canard-controlled store from an open-flow rectangular cavity bay at a freestream Mach number of 1.2. The transient flowfield was computed using a diagonal scheme in an overset mesh framework, with the resultant aerodynamic loads used as the forcing functions in the nonlinear dynamics equations. The proportional and rate gyro sensitivities were computed a priori using pole placement techniques for the linearized dynamical equations. These fixed gain values were used in the controller for the nonlinear simulation. Reasonable comparison between the full and linearized equations for a perturbed two-dimensional missile was found. Also in two-dimensions, a controlled store was found to possess improved separation characteristics over a canard-fixed store. In three-dimensions, trajectory comparisons with wind-tunnel data for the canard-fixed case will be made. In addition, it will be determined if a canard-controlled store is an effective means of improving cavity store separation characteristics.

Atwood, Christopher A.↗

Evaluation of improved algorithms for nonlinear transient structural heat transfer

The paper presents the evaluation of algorithms having potential for reducing computer time to integrate large systems of equations governing nonlinear transient heat transfer in complex structures. A series of test problems is being developed to test two algorithms. The Shuttle Orbiter frame model having 190 nodes was the subject of an analytical/experimental study. Plans are to perform calculations for the complete wing model having about 4000 nodes as well as a regional (plug) model having about 200 nodes.

Adelman, H. M.↗

GCKP84--General Chemical Kinetics Code for Gas-Phase Flow and Batch Processes Including Heat Transfer Effects

A general chemical kinetics code is described for complex, homogeneous ideal gas reactions in any chemical system. The main features of the GCKP84 code are flexibility, convenience, and speed of computation for many different reaction conditions. The code, which replaces the GCKP code published previously, solves numerically the differential equations for complex reaction in a batch system or one dimensional inviscid flow. It also solves numerically the nonlinear algebraic equations describing the well stirred reactor. A new state of the art numerical integration method is used for greatly increased speed in handling systems of stiff differential equations. The theory and the computer program, including details of input preparation and a guide to using the code are given.

David A Bittker↗

Stationary nonlinear Alfven waves and solitons

Stationary solutions of the derivative nonlinear Schroedinger equation are discussed and classified by using a pseudopotential formulation. The solutions consist of a rich family of nonlinear Alfven waves and solitons with parallel and oblique propagation directions. Expressions for the envelope and the phase of nonlinear waves with periodic envelope modulation, and 'hyperbolic' and 'algebraic' solitons are given. The propagation angle for the slightly modulated elliptic, periodic waves and for oblique solitons is evaluated.

Hada, T.↗

Calculation of shocked one-dimensional flows on abruptly changing grids by mathematical programming

Cell-centered finite differences on cells of an arbitrarily spaced grid are presently used to discretize the steady-state inviscid and nearly inviscid Burgers' equations for steady-state shocked flow in a quasi-one-dimensional nozzle. The everdetermined system of nonlinear algebraic equations thus obtained is solved by a procedure which minimizes a weighted sum of the residuals' absolute values. Numerical solutions for both viscous and inviscid problems are accurate and nonoscillatory, on grids whose abrupt mesh lengths refinements are as great as a factor of 10,000.

Lavery, John E.↗

The relative diffusion of strong magnetic fields and tenuous gases

General equations are given for treating the diffusion of a magnetic field in a fluid with electrical conductivity and a resistive diffusion coefficient. Consequences of nonnegligible gas motion are investigated by considering the one-dimensional problem where a tenuous isothermal gas with pressure expressed as a function of x-coordinate and time is enmeshed in a magnetic field in the y-direction. Approximate equations are obtained for the case of weak fluctuations in gas or field pressure, and the properties of a nonlinear diffusion equation are illustrated with the example of an initial thin slab of gas confined by a uniform field on either side. The advance of a gas front into a nonvanishing gas density is illustrated by using the progressive wave solution. Application of the results to the evolution of streamers and filaments observed in the solar corona is briefly discussed

Parker, E. N.↗

Computation of Static Shapes and Voltages for Micromachined Deformable Mirrors with Nonlinear Electrostatic Actuators

In modeling micromachined deformable mirrors with electrostatic actuators whose gap spacings are of the same order of magnitude as those of the surface deformations, it is necessary to use nonlinear models for the actuators. In this paper, we consider micromachined deformable mirrors modeled by a membrane or plate equation with nonlinear electrostatic actuator characteristics. Numerical methods for computing the mirror deformation due to given actuator voltages and the actuator voltages required for producing the desired deformations at the actuator locations are presented. The application of the proposed methods to circular deformable mirrors whose surfaces are modeled by elastic membranes is discussed in detail. Numerical results are obtained for a typical circular micromachined mirror with electrostatic actuators.

Wang, P. K. C.↗

An O(Nm(sup 2)) Plane Solver for the Compressible Navier-Stokes Equations

A hierarchical multigrid algorithm for efficient steady solutions to the two-dimensional compressible Navier-Stokes equations is developed and demonstrated. The algorithm applies multigrid in two ways: a Full Approximation Scheme (FAS) for a nonlinear residual equation and a Correction Scheme (CS) for a linearized defect correction implicit equation. Multigrid analyses which include the effect of boundary conditions in one direction are used to estimate the convergence rate of the algorithm for a model convection equation. Three alternating-line- implicit algorithms are compared in terms of efficiency. The analyses indicate that full multigrid efficiency is not attained in the general case; the number of cycles to attain convergence is dependent on the mesh density for high-frequency cross-stream variations. However, the dependence is reasonably small and fast convergence is eventually attained for any given frequency with either the FAS or the CS scheme alone. The paper summarizes numerical computations for which convergence has been attained to within truncation error in a few multigrid cycles for both inviscid and viscous ow simulations on highly stretched meshes.

Thomas, J. L.↗

Dynamic Theory of Relativistic Electrons Stochastic Heating by Whistler Mode Waves with Application to the Earth Magnetosphere

In the Hamiltonian approach an electron motion in a coherent packet of the whistler mode waves propagating along the direction of an ambient magnetic field is studied. The physical processes by which these particles are accelerated to high energy are established. Equations governing a particle motion by group symmetries of the problem were transformed in to a closed pair of nonlinear difference equations. The solutions of these equations have shown there exists the energetic threshold below that the electron motion is regular, and when the initial energy is above the threshold an electron moves stochastically. It is proved that the upper boundary of particle stochastic heating is conditioned by intrinsic property of the particle chaotic motion. Particle energy spectra and pitch angle electron scattering are described by the Fokker-Planck-Kolmogorov equations. It is shown that significant pitch angle diffusion occurs for the Earth radiation belt electrons with energies from a few keV up to a few MeV.

Khazanov, G. V.↗

Computation of Static Shapes and Voltages for Micromachined Deformable Mirrors with Nonlinear Electostatic Actuators

In modeling micromachined deformable mirrors with electrostatic actuators whose gap spacings are of the same order of magnitude as that of the surface deformations, it is neccessary to use nonlinear models for the actuators. In this paper, we consider micromachined deformable mirrors modeled by a membrane or plate equation with nonlinear electrostatic actuator characteristics.

deformations optical systems adaptive optical syst↗

A new class of exact, nonlinear solutions to the Grad-Shafranov equation

We have constructed a new class of exact, nonlinear solutions to the Grad-Shafranov equation, representing force-free magnetic fields with translational symmetry. These exact solutions are pertinent to the study of magnetic structures in the solar corona that are subjected to photospheric shearing motions.

Roumeliotis, George↗

Numerical solution of the Navier-Stokes equations for high Reynolds number incompressible turbulent flow

The full Navier-Stokes equations for incompressible turbulent flow must be solved to accurately represent all flow phenomena which occur in a high Reynolds number incompressible flow. A two layer algebraic eddy viscosity turbulence model is used to represent the Reynolds stress in the primitive variable formulation. The development of the boundary-fitted coordinate systems makes the numerical solution of these equations feasible for arbitrarily shaped bodies. The nondimensional time averaged Navier-Stokes equations, including the turbulence mode, are represented by finite difference approximations in the transformed plane. The resulting coupled system of nonlinear algebraic equations is solved using a point successive over relaxation iteration. The test case considered was a NACA 64A010 airfoil section at an angle of attack of two degrees and a Reynolds number of 2,000,000.

Thompson, D. S.↗

Parallel processors and nonlinear structural dynamics algorithms and software

Techniques are discussed for the implementation and improvement of vectorization and concurrency in nonlinear explicit structural finite element codes. In explicit integration methods, the computation of the element internal force vector consumes the bulk of the computer time. The program can be efficiently vectorized by subdividing the elements into blocks and executing all computations in vector mode. The structuring of elements into blocks also provides a convenient way to implement concurrency by creating tasks which can be assigned to available processors for evaluation. The techniques were implemented in a 3-D nonlinear program with one-point quadrature shell elements. Concurrency and vectorization were first implemented in a single time step version of the program. Techniques were developed to minimize processor idle time and to select the optimal vector length. A comparison of run times between the program executed in scalar, serial mode and the fully vectorized code executed concurrently using eight processors shows speed-ups of over 25. Conjugate gradient methods for solving nonlinear algebraic equations are also readily adapted to a parallel environment. A new technique for improving convergence properties of conjugate gradients in nonlinear problems is developed in conjunction with other techniques such as diagonal scaling. A significant reduction in the number of iterations required for convergence is shown for a statically loaded rigid bar suspended by three equally spaced springs.

Belytschko, Ted↗

A linearized Euler analysis of unsteady flows in turbomachinery

A method for calculating unsteady flows in cascades is presented. The model, which is based on the linearized unsteady Euler equations, accounts for blade loading shock motion, wake motion, and blade geometry. The mean flow through the cascade is determined by solving the full nonlinear Euler equations. Assuming the unsteadiness in the flow is small, then the Euler equations are linearized about the mean flow to obtain a set of linear variable coefficient equations which describe the small amplitude, harmonic motion of the flow. These equations are discretized on a computational grid via a finite volume operator and solved directly subject to an appropriate set of linearized boundary conditions. The steady flow, which is calculated prior to the unsteady flow, is found via a Newton iteration procedure. An important feature of the analysis is the use of shock fitting to model steady and unsteady shocks. Use of the Euler equations with the unsteady Rankine-Hugoniot shock jump conditions correctly models the generation of steady and unsteady entropy and vorticity at shocks. In particular, the low frequency shock displacement is correctly predicted. Results of this method are presented for a variety of test cases. Predicted unsteady transonic flows in channels are compared to full nonlinear Euler solutions obtained using time-accurate, time-marching methods. The agreement between the two methods is excellent for small to moderate levels of flow unsteadiness. The method is also used to predict unsteady flows in cascades due to blade motion (flutter problem) and incoming disturbances (gust response problem).

Hall, Kenneth C.↗

Calculation of unsteady flows in turbomachinery using the linearized Euler equations

A method for calculating unsteady flows in cascades is presented. The model, which is based on the linearized unsteady Euler equations, accounts for blade loading shock motion, wake motion, and blade geometry. The mean flow through the cascade is determined by solving the full nonlinear Euler equations. Assuming the unsteadiness in the flow is small, then the Euler equations are linearized about the mean flow to obtain a set of linear variable coefficient equations which describe the small amplitude, harmonic motion of the flow. These equations are discretized on a computational grid via a finite volume operator and solved directly subject to an appropriate set of linearized boundary conditions. The steady flow, which is calculated prior to the unsteady flow, is found via a Newton iteration procedure. An important feature of the analysis is the use of shock fitting to model steady and unsteady shocks. Use of the Euler equations with the unsteady Rankine-Hugoniot shock jump conditions correctly models the generation of steady and unsteady entropy and vorticity at shocks. In particular, the low frequency shock displacement is correctly predicted. Results of this method are presented for a variety of test cases. Predicted unsteady transonic flows in channels are compared to full nonlinear Euler solutions obtained using time-accurate, time-marching methods. The agreement between the two methods is excellent for small to moderate levels of flow unsteadiness. The method is also used to predict unsteady flows in cascades due to blade motion (flutter problem) and incoming disturbances (gust response problem).

Hall, Kenneth C.↗