Kinetic theory of plasma in a magnetic field
Guiding center equations for charged particle motion and statistical fluctuations of plasma in magnetic field, using Krylov-Bogoliubov transformation
Price, J. C.↗
SEARCH · Engineering Papers
Results for “Krylov”
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.
The influence of the great inequality on the secular disturbing function of the planetary system.
This paper derives the contribution by the great inequality to the secular disturbing function of the principal planets. Andoyer's expansion of the planetary disturbing function and von Zeipel's method of eliminating the periodic terms is employed; thereby, the corrected secular disturbing function for the planetary system is derived. The conclusion is drawn that the canonicity of the equations for the secular variation of the heliocentric elements can be preserved if there be retained, in the secular disturbing function, terms only of the second and fourth order relative to the eccentricity and inclinations. The Krylov-Bogoliubov method is suggested for eliminating periodic terms, if it is desired to include the secular perturbations of the fifth and higher order in the heliocentric elements. The additional part of the secular disturbing function derived in this paper can be included in existing theories of the secular effects of principal planets.
Musen, P.↗
Flutter of articulated pipes at finite amplitude
The plane motion of an articulated pipe made of two segments is examined and the flow velocity at which flutter manifests itself is sought. The pressure in the reservoir feeding the pipe is kept constant. In contrast to previous works, the flow velocity is not taken as a prescribed parameter of the system but is left to follow the laws of motion. This approach requires a nonlinear formulation of the problem and the equations of motion are solved using Krylov-Bogoliubov's method. A graph of the amplitude of the limit cycles, as a function of the fluid-system mass ratio, is presented and conclusions are drawn as to the necessity of considering nonlinearities in the analysis.
Rousselet, J.↗
A note on obtaining an approximate value for the period of a periodic solution of x-double prime = -V prime /x/
The first integral of x-double prime = - V prime (x) yields an integral for the period of a periodic solution, if such exists. In general, this integral cannot be evaluated. By means of an approximate solution along with the minimization of a mean-square error, one can obtain an approximate value for the period in terms of the amplitude of the motion. The calculated period agrees very well with the period obtained by means of numerical integration for the case of orbit-orbit resonance involving the motion of two satellites of a planet. The same method is applied to obtain an alternative derivation of the first Krylov-Bogoliuboff averaging method in nonlinear mechanics.
Lass, H.↗
Dynamic behavior of continuous cantilevered pipes conveying fluid near critical velocities
The plane motion of a cantilevered pipe conveying fluid is examined when the flow velocity is in the neighborhood of that generating flutter. In contrast to previous studies, the flow velocity is not prescribed as a constant, but is determined from the laws of motion. A system of two nonlinear partial differential equations which are coupled through the nonlinear terms is thereby obtained. The solution is found by the use of the Krylov-Bogoliubov averaging method and the results are discussed indicating the effect of nonlinearities.
Rousselet, J.↗
The algebraic criteria for the stability of control systems
This paper critically examines the standard algebraic criteria for the stability of linear control systems and their proofs, reveals important previously unnoticed connections, and presents new representations. Algebraic stability criteria have also acquired significance for stability studies of non-linear differential equation systems by the Krylov-Bogoljubov-Magnus Method, and allow realization conditions to be determined for classes of broken rational functions as frequency characteristics of electrical network.
Cremer, H.↗
TRANAIR - A computer code for transonic analyses of arbitrary configurations
Attention is given to a new approach to solving full potential equations about arbitrary configurations. Numerical algorithms from such fields as finite elements, preconditioned Krylov subspace methods, discrete Fourier analysis, and integral equations are combined to take advantage of the size and speed of current and emerging supercomputers. On the basis of this appraoch, a robust, efficient and easy to use computer code referred to as TRANAIR has been developed for transonic analysis of complex geometries.
Samant, S. S.↗
A new approach to the solution of boundary value problems involving complex configurations
A new approach for solving certain types of boundary value problems about complex configurations is presented. Numerical algorithms from such diverse fields as finite elements, preconditioned Krylov subspace methods, discrete Fourier analysis, and integral equations are combined to take advantage of the memory, speed and architecture of current and emerging supercomputers. Although the approach has application to many branches of computational physics, the present effort is concentrated in areas of Computational Fluid Dynamics (CFD) such as steady nonlinear aerodynamics, time harmonic unsteady aerodynamics, and aeroacoustics. The most significant attribute of the approach is that it can handle truly arbitrary boundary geometries and eliminates the difficult task of generating surface fitted grids.
Rubbert, P. E.↗
Effects of nonlinear damping in flexible space structures
The classical Krylov-Bogoliubov "averaging" technique is used to study a class of nonlinear damping models, for which the damping force is proportional to the product of positive integer or fractional power of absolute values of displacement and that of velocity. The results are compared with linear viscous damping models. The amplitude decrement of free vibration for a single mode system with nonlinear models depends not only on damping ratio, but also on the initial amplitude, the time to measure the response, frequency of the system, and the powers of displacement and velocity. For the distributed system, the action of nonlinear damping is found to reduce energy of the system as well as to pass energy to higher modes. Experimental evidence such as in Spacecraft Control Laboratory Experiment seems to support the need for nonlinear models.
Hu, Anren↗
Some nonlinear damping models in flexible structures
A class of nonlinear damping models is introduced with application to flexible flight structures characterized by low damping. Approximate solutions of engineering interest are obtained for the model using the classical averaging technique of Krylov and Bogoliubov. The results should be considered preliminary pending further investigation.
Balakrishnan, A. V.↗
Application of the TRANAIR rectangular grid approach to the aerodynamic analysis of complex configurations
A numerical method is described which uses a rectangular grid to solve the nonlinear full potential equation about complex configurations. The grid is locally refined to resolve high velocity gradients arising from leading edge expansions or shock waves. The grid penetrates the boundary (described by networks of quadrilateral panels) and is generated automatically. Discrete operators are constructed using the finite element method. The system of nonlinear discrete equations is solved iteratively using a Krylov subspace method preconditioned by an exterior Poisson solver and a direct sparse solver. The primary emphasis is to provide design engineers with an aerodynamic analysis tool (the TRANAIR code) which is accurate, reliable, economical, and flexible to use. Computational results for many interesting configurations are presented.
Johnson, Forrester T.↗
Parallel iterative methods for sparse linear and nonlinear equations
As three-dimensional models are gaining importance, iterative methods will become almost mandatory. Among these, preconditioned Krylov subspace methods have been viewed as the most efficient and reliable, when solving linear as well as nonlinear systems of equations. There has been several different approaches taken to adapt iterative methods for supercomputers. Some of these approaches are discussed and the methods that deal more specifically with general unstructured sparse matrices, such as those arising from finite element methods, are emphasized.
Saad, Youcef↗
A decentralized linear quadratic control design method for flexible structures
A decentralized suboptimal linear quadratic control design procedure which combines substructural synthesis, model reduction, decentralized control design, subcontroller synthesis, and controller reduction is proposed for the design of reduced-order controllers for flexible structures. The procedure starts with a definition of the continuum structure to be controlled. An evaluation model of finite dimension is obtained by the finite element method. Then, the finite element model is decomposed into several substructures by using a natural decomposition called substructuring decomposition. Each substructure, at this point, still has too large a dimension and must be reduced to a size that is Riccati-solvable. Model reduction of each substructure can be performed by using any existing model reduction method, e.g., modal truncation, balanced reduction, Krylov model reduction, or mixed-mode method. Then, based on the reduced substructure model, a subcontroller is designed by an LQ optimal control method for each substructure independently. After all subcontrollers are designed, a controller synthesis method called substructural controller synthesis is employed to synthesize all subcontrollers into a global controller. The assembling scheme used is the same as that employed for the structure matrices. Finally, a controller reduction scheme, called the equivalent impulse response energy controller (EIREC) reduction algorithm, is used to reduce the global controller to a reasonable size for implementation. The EIREC reduced controller preserves the impulse response energy of the full-order controller and has the property of matching low-frequency moments and low-frequency power moments. An advantage of the substructural controller synthesis method is that it relieves the computational burden associated with dimensionality. Besides that, the SCS design scheme is also a highly adaptable controller synthesis method for structures with varying configuration, or varying mass and stiffness properties.
Su, Tzu-Jeng↗
Special class of nonlinear damping models in flexible space structures
A special class of nonlinear damping models is investigated in which the damping force is proportional to the product of positive integer or the fractional power of the absolute values of displacement and velocity. For a one-degree-of-freedom system, the classical Krylov-Bogoliubov 'averaging' method is used, whereas for a distributed system, both an ad hoc perturbation technique and the finite difference method are employed to study the effects of nonlinear damping. The results are compared with linear viscous damping models. The amplitude decrement of free vibration for a single mode system with nonlinear models depends not only on the damping ratio but also on the initial amplitude, the time to measure the response, the frequency of the system, and the powers of displacement and velocity. For the distributed system, the action of nonlinear damping is found to reduce the energy of the system and to pass energy to lower modes.
Hu, Anren↗
Numerical methods in Markov chain modeling
Several methods for computing stationary probability distributions of Markov chains are described and compared. The main linear algebra problem consists of computing an eigenvector of a sparse, usually nonsymmetric, matrix associated with a known eigenvalue. It can also be cast as a problem of solving a homogeneous singular linear system. Several methods based on combinations of Krylov subspace techniques are presented. The performance of these methods on some realistic problems are compared.
Philippe, Bernard↗
Projection methods for the numerical solution of Markov chain models
Projection methods for computing stationary probability distributions for Markov chain models are presented. A general projection method is a method which seeks an approximation from a subspace of small dimension to the original problem. Thus, the original matrix problem of size N is approximated by one of dimension m, typically much smaller than N. A particularly successful class of methods based on this principle is that of Krylov subspace methods which utilize subspaces of the form span(v,av,...,A(exp m-1)v). These methods are effective in solving linear systems and eigenvalue problems (Lanczos, Arnoldi,...) as well as nonlinear equations. They can be combined with more traditional iterative methods such as successive overrelaxation, symmetric successive overrelaxation, or with incomplete factorization methods to enhance convergence.
Saad, Youcef↗
Rational approximations from power series of vector-valued meromorphic functions
Let F(z) be a vector-valued function, F: C yields C(sup N), which is analytic at z = 0 and meromorphic in a neighborhood of z = 0, and let its Maclaurin series be given. In this work we developed vector-valued rational approximation procedures for F(z) by applying vector extrapolation methods to the sequence of partial sums of its Maclaurin series. We analyzed some of the algebraic and analytic properties of the rational approximations thus obtained, and showed that they were akin to Pade approximations. In particular, we proved a Koenig type theorem concerning their poles and a de Montessus type theorem concerning their uniform convergence. We showed how optical approximations to multiple poles and to Laurent expansions about these poles can be constructed. Extensions of the procedures above and the accompanying theoretical results to functions defined in arbitrary linear spaces was also considered. One of the most interesting and immediate applications of the results of this work is to the matrix eigenvalue problem. In a forthcoming paper we exploited the developments of the present work to devise bona fide generalizations of the classical power method that are especially suitable for very large and sparse matrices. These generalizations can be used to approximate simultaneously several of the largest distinct eigenvalues and corresponding eigenvectors and invariant subspaces of arbitrary matrices which may or may not be diagonalizable, and are very closely related with known Krylov subspace methods.
Sidi, Avram↗
Domain decomposition: A bridge between nature and parallel computers
Domain decomposition is an intuitive organizing principle for a partial differential equation (PDE) computation, both physically and architecturally. However, its significance extends beyond the readily apparent issues of geometry and discretization, on one hand, and of modular software and distributed hardware, on the other. Engineering and computer science aspects are bridged by an old but recently enriched mathematical theory that offers the subject not only unity, but also tools for analysis and generalization. Domain decomposition induces function-space and operator decompositions with valuable properties. Function-space bases and operator splittings that are not derived from domain decompositions generally lack one or more of these properties. The evolution of domain decomposition methods for elliptically dominated problems has linked two major algorithmic developments of the last 15 years: multilevel and Krylov methods. Domain decomposition methods may be considered descendants of both classes with an inheritance from each: they are nearly optimal and at the same time efficiently parallelizable. Many computationally driven application areas are ripe for these developments. A progression is made from a mathematically informal motivation for domain decomposition methods to a specific focus on fluid dynamics applications. To be introductory rather than comprehensive, simple examples are provided while convergence proofs and algorithmic details are left to the original references; however, an attempt is made to convey their most salient features, especially where this leads to algorithmic insight.
Keyes, David E.↗