Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear differential 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 487 records · Page 27

Calculation of sheath and wake structure about a pillbox-shaped spacecraft in a flowing plasma

A computer program was used for studies of the disturbed zones around bodies in flowing plasmas, particularly spacecraft and their associated sheaths and wakes. The program solved a coupled Poisson-Vlasov system of nonlinear partial differential integral equations to obtain distributions of electric potential and ion and electron density about a finite length cylinder in a plasma flow at arbitrary ion Mach numbers. The approach was applicable to a larger range of parameters than other available approaches. In sample calculations, bodies up to 100 Debye lengths in radius were treated, that is, larger than any previously treated realistically. Applications were made to in-situ satellite experiments.

Parker, L. W.↗

A Galerkin type finite element method for rotary-wing aeroelasticity in hover and forward flight

A Galerkin finite element method for the spatial discretization of the nonlinear, nonselfadjoint, partial differential equations governing rotary-wing aeroelasticity is presented. This method reduces algebraic manipulative labor significantly when compared to the global Galerkin method based on assumed modes. Furthermore, the Galerkin finite element method is ideally suited to treat rotor blades with discontinuous mass and stiffness distribution and structurally redundant configurations as they appear in bearingless rotors. Implementation of the method is illustrated for the coupled flap-lag aeroelastic problem of hingeless rotor blades in hover and forward flight. Numerical results for stability and response illustrate the numerical properties and convergence behavior of the method. It is concluded that the Galerkin finite element method is a practical tool for solving rotary-wing aeroelastic stability and response problems.

Straub, F. K.↗

An approximation theory for the identification of nonlinear distributed parameter systems

An abstract approximation framework for the identification of nonlinear distributed parameter systems is developed. Inverse problems for nonlinear systems governed by strongly maximal monotone operators (satisfying a mild continuous dependence condition with respect to the unknown parameters to be identified) are treated. Convergence of Galerkin approximations and the corresponding solutions of finite dimensional approximating identification problems to a solution of the original finite dimensional identification problem is demonstrated using the theory of nonlinear evolution systems and a nonlinear analog of the Trotter-Kato approximation result for semigroups of bounded linear operators. The nonlinear theory developed here is shown to subsume an existing linear theory as a special case. It is also shown to be applicable to a broad class of nonlinear elliptic operators and the corresponding nonlinear parabolic partial differential equations to which they lead. An application of the theory to a quasilinear model for heat conduction or mass transfer is discussed.

Banks, H. T.↗

An approximation theory for the identification of nonlinear distributed parameter systems

An abstract approximation framework for the identification of nonlinear distributed parameter systems is developed. Inverse problems for nonlinear systems governed by strongly maximal monotone operators (satisfying a mild continuous dependence condition with respect to the unknown parameters to be identified) are treated. Convergence of Galerkin approximations and the corresponding solutions of finite dimensional approximating identification problems to a solution of the original finite dimensional identification problem is demonstrated using the theory of nonlinear evolution systems and a nonlinear analog of the Trotter-Kato appproximation result for semigroups of bounded linear operators. The nonlinear theory developed here is shown to subsume an existing linear theory as a special case. It is also shown to be applicable to a broad class of nonlinear elliptic operators and the corresponding nonlinear parabolic partial differential equations to which they lead. An application of the theory to a quasilinear model for heat conduction or mass transfer is discussed.

Banks, H. T.↗

The interaction of an artificial electron beam with the earth's upper atmosphere - Effects on spacecraft charging and the near-plasma evironment

The ionization created in the earth's upper atmosphere by an electron beam emitted from a low altitude spacecraft is described by two coupled nonlinear differential electron transport equations for the up-going and down-going differential energy flux. The equations are solved numerically, yielding estimates of the differential electron energy flux density at the spacecraft location. At altitudes below 200-250 km and for beam energies around 1 keV, it is shown that secondary electrons supply a significant contribution to the return current to the spacecraft and thereby reduce the spacecraft potential. The numerical results are in good agreement with observations from the Charge 2 sounding rocket experiment.

Neubert, T.↗

Spurious Numerical Solutions Of Differential Equations

Paper presents detailed study of spurious steady-state numerical solutions of differential equations that contain nonlinear source terms. Main objectives of this study are (1) to investigate how well numerical steady-state solutions of model nonlinear reaction/convection boundary-value problem mimic true steady-state solutions and (2) to relate findings of this investigation to implications for interpretation of numerical results from computational-fluid-dynamics algorithms and computer codes used to simulate reacting flows.

Lafon, A.↗

Fast Multigrid Techniques in Total Variation-Based Image Reconstruction

Existing multigrid techniques are used to effect an efficient method for reconstructing an image from noisy, blurred data. Total Variation minimization yields a nonlinear integro-differential equation which, when discretized using cell-centered finite differences, yields a full matrix equation. A fixed point iteration is applied with the intermediate matrix equations solved via a preconditioned conjugate gradient method which utilizes multi-level quadrature (due to Brandt and Lubrecht) to apply the integral operator and a multigrid scheme (due to Ewing and Shen) to invert the differential operator. With effective preconditioning, the method presented seems to require Omicron(n) operations. Numerical results are given for a two-dimensional example.

Oman, Mary Ellen↗

Numerical solution of flame sheet problems with and without multigrid methods

Flame sheet problems are on the natural route to the numerical solution of multidimensional flames, which, in turn, are important in many engineering applications. In order to model the structure of flames more accurately, we use the vorticity-velocity formulation of the fluid flow equations, as opposed to the streamfunction-vorticity approach. The numerical solution of the resulting nonlinear coupled elliptic partial differential equations involves a pseudo transient process and a steady state Newton iteration. Rather than working with dimensionless variables, we introduce scale factors that can yield significant savings in the execution time. In this context, we also investigate the applicability and performance of several multigrid methods, focusing on nonlinear damped Newton multigrid, using either one way or correction schemes.

Douglas, Craig C.↗

Burnett solutions along the stagnation line of a cooled cylinder in low-density hypersonic flows

The Burnett equations in the cylindrical coordinate system have been applied to the hypersonic flow over a cylinder in low density environments. Through a local similarity transformation, the Burnett equations are degenerated into nonlinear third-order ordinary differential equations which are only valid along the stagnation line. These equations are solved by the multi-staged Runge-Kutta integration technique. Steady state solutions are asymptotically reached by the time-marching procedure. A newly-developed computer code is used to solve both the Navier-Stokes and the Burnett equations within the shock layer. The computed heat transfer coefficients are consistent with wind-tunnel and STS flight data. Comparisons with the DSMC solutions and the Boltzmann solutions show that the Burnett equations are more appropriate as the governing equations than the Navier-Stokes equations in the transitional flow regime.

Liaw, G. S.↗

Parametric Identification of Nonlinear Dynamical Systems

In this project, we looked at the application of harmonic balancing as a tool for identifying parameters (HBID) in a nonlinear dynamical systems with chaotic responses. The main idea is to balance the harmonics of periodic orbits extracted from measurements of each coordinate during a chaotic response. The periodic orbits are taken to be approximate solutions to the differential equations that model the system, the form of the differential equations being known, but with unknown parameters to be identified. Below we summarize the main points addressed in this work. The details of the work are attached as drafts of papers, and a thesis, in the appendix. Our study involved the following three parts: (1) Application of the harmonic balance to a simulation case in which the differential equation model has known form for its nonlinear terms, in contrast to a differential equation model which has either power series or interpolating functions to represent the nonlinear terms. We chose a pendulum, which has sinusoidal nonlinearities; (2) Application of the harmonic balance to an experimental system with known nonlinear forms. We chose a double pendulum, for which chaotic response were easily generated. Thus we confronted a two-degree-of-freedom system, which brought forth challenging issues; (3) A study of alternative reconstruction methods. The reconstruction of the phase space is necessary for the extraction of periodic orbits from the chaotic responses, which is needed in this work. Also, characterization of a nonlinear system is done in the reconstructed phase space. Such characterizations are needed to compare models with experiments. Finally, some nonlinear prediction methods can be applied in the reconstructed phase space. We developed two reconstruction methods that may be considered if the common method (method of delays) is not applicable.

Feeny, Brian↗

Monolithic Multigrid Methods for Magnetohydrodynamics

The magnetohydrodynamics equations model a wide range of plasma physics applications and are characterized by a nonlinear system of partial differential equations that strongly couples a charged fluid with the evolution of electromagnetic fields. After discretization and linearization, the resulting system of equations is generally difficult to solve due to the coupling between variables and the heterogeneous coefficients induced by the linearization process. In this paper, we investigate multigrid preconditioners for this system based on specialized relaxation schemes that properly address the system structure and coupling. Here, three extensions of Vanka relaxation are proposed and applied to problems with up to 170 million degrees of freedom and fluid and magnetic Reynolds numbers up to 400 for stationary problems and up to 20,000 for time-dependent problems.

97 MATHEMATICS AND COMPUTING↗

Nonlinear evolution of magnetic flux ropes. I - Low-beta limit

We study the nonlinear self-similar evolution of a cylindrical magnetic flux tube with two components of the magnetic field, axial and azimuthal. We restrict ourselves to the case of a plasma of low beta. Introducing a special class of configurations we call 'separable fields', we reduce the problem to an ordinary differential equation. Two cases are to be distinguished: (1) when the total field minimizes on the symmetry axis, the magnetic configuration inexorably collapses, and (2) when, on the other hand, the total field maximizes on the symmetry axis, the magnetic configuration behaves analogously to a nonlinear oscillator. Here we focus on the latter case. The effective potential of the motion contains two terms: a strong repulsive term and a weak restoring term associated with the pinch. We solve the nonlinear differential equation of motion numerically and find that the period of oscillations grows exponentially with the energy of the oscillator. Our treatment emphasizes the role of the force-free configuration as the lowest potential energy state about which the system oscillates.

Osherovich, V. A.↗

Nonstandard finite difference schemes

The major research activities of this proposal center on the construction and analysis of nonstandard finite-difference schemes for ordinary and partial differential equations. In particular, we investigate schemes that either have zero truncation errors (exact schemes) or possess other significant features of importance for numerical integration. Our eventual goal is to bring these methods to bear on problems that arise in the modeling of various physical, engineering, and technological systems. At present, these efforts are extended in the direction of understanding the exact nature of these nonstandard procedures and extending their use to more complicated model equations. Our presentation will give a listing (obtained to date) of the nonstandard rules, their application to a number of linear and nonlinear, ordinary and partial differential equations. In certain cases, numerical results will be presented.

Mickens, Ronald E.↗

Simulation of multicorrelated random processes using the FFT algorithm

A technique for the digital simulation of multicorrelated Gaussian random processes is described. This technique is based upon generating discrete frequency functions which correspond to the Fourier transform of the desired random processes, and then using the fast Fourier transform (FFT) algorithm to obtain the actual random processes. The main advantage of this method of simulation over other methods is computation time; it appears to be more than an order of magnitude faster than present methods of simulation. One of the main uses of multicorrelated simulated random processes is in solving nonlinear random vibration problems by numerical integration of the governing differential equations. The response of a nonlinear string to a distributed noise input is presented as an example.

Wittig, L. E.↗

Krylov subspace methods - Theory, algorithms, and applications

Projection methods based on Krylov subspaces for solving various types of scientific problems are reviewed. The main idea of this class of methods when applied to a linear system Ax = b, is to generate in some manner an approximate solution to the original problem from the so-called Krylov subspace span. Thus, the original problem of size N is approximated by one of dimension m, typically much smaller than N. Krylov subspace methods have been very successful in solving linear systems and eigenvalue problems and are now becoming popular for solving nonlinear equations. The main ideas in Krylov subspace methods are shown and their use in solving linear systems, eigenvalue problems, parabolic partial differential equations, Liapunov matrix equations, and nonlinear system of equations are discussed.

Sad, Youcef↗

Difference methods for problems with different time scales

The use of difference methods for weakly nonlinear systems of ordinary differential equations with rapidly oscillating solutions is considered, and a general approach which depends only on the smoothness of the coefficients and the nonlinearities is developed. In particular, one is led to a strategy which is suitable for the detection and resolution of turning points in such systems. A computational example is presented.

Scheid, R. E., Jr.↗

Adaptive PID Gain Scheduling Control for Hydropower Turbine Using Neural CDE and Stochastic Distribution Shaping

This paper introduces a gain-scheduling PID controller design strategy for hydroturbine frequency control mode. This scheme first uses real data to learn the nonlinear dynamics of the hydroturbine using neural controlled differential equations and then perturbs the obtained nonlinear system at different equilibrium points, based on which a static output feedback adaptive dynamic programming algorithm is then used to optimize the PID gains for each equilibrium point. Moreover, a continuous-time version of stochastic distribution control is proposed to further fine-tune the optimized PID gains. Finally, the controller is obtained by implementing linear interpolation between the optimized PID control gains. The simulation results show that the proposed gain-scheduling PID controller can control a larger range of operation points compared with the given fixed PID controller and the baseline method. Compared with the given fixed PID controller, the proposed gain-scheduling PID controller can regulate hydroturbine frequency against disturbances induced by power-load variation with over 50% less overshoot for some operation points.

13 HYDRO ENERGY↗