A linearization tool for use with matrix formalisms of rotational dynamics.
Rotational dynamics equations in linearized matrix form applied to satellite librations in circular orbit and gyroscope motion
Likins, P. W.↗
SEARCH · Engineering Papers
Results for “linear 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.
Chaotic /strange/ and periodic behavior in instability saturation by the oscillating two-stream instability
The nonlinear Schroedinger equation with linear growth and damping is truncated to three waves. The resulting system of nonlinear ordinary differential equations describes the excitation of linearly damped waves by the oscillating two-stream instability driven by a linearly unstable pump wave. This system represents a simple model for the nonlinear saturation of a linearly unstable wave. The model is examined analytically and numerically as a function of the dimensionless parameters of the system. It is found that the model can exhibit a wealth of characteristic dynamical behavior including stationary equilibria, Hopf bifurcations to periodic orbits, period doubling bifurcations, chaotic solutions characteristic of a strange attractor, tangent bifurcations from chaotic to periodic solutions, transient chaos, and hysteresis. Many of these features are shown to be explainable on the basis of one-dimensional maps. In the case of chaotic solutions, evidence for the presence of a strange attractor is provided by demonstrating Cantor set-like structure (i.e., scale invariance) in the surface of section.
Russell, D. A.↗
Parallels between control PDE's (Partial Differential Equations) and systems of ODE's (Ordinary Differential Equations)
System theorists understand that the same mathematical objects which determine controllability for nonlinear control systems of ordinary differential equations (ODEs) also determine hypoellipticity for linear partial differentail equations (PDEs). Moreover, almost any study of ODE systems begins with linear systems. It is remarkable that Hormander's paper on hypoellipticity of second order linear p.d.e.'s starts with equations due to Kolmogorov, which are shown to be analogous to the linear PDEs. Eigenvalue placement by state feedback for a controllable linear system can be paralleled for a Kolmogorov equation if an appropriate type of feedback is introduced. Results concerning transformations of nonlinear systems to linear systems are similar to results for transforming a linear PDE to a Kolmogorov equation.
Hunt, L. R.↗
Parallels between control PDE's and systems of ODE's
System theorists understand that the same mathematical objects which determine controllability for nonlinear control systems of ordinary differential equations (ODEs) also determine hypoellipticity for linear partial differential equations (PDEs). Moreover, almost any study of ODE systems begins with linear systems. It is remarkable that Hormander's paper on hypoellipticity of second order linear p.d.e.'s starts with equations due to Kolmogorov, which are shown to be analogous to the linear PDEs. Eigenvalue placement by state feedback for a controllable linear system can be paralled for a Kolmogorov equation if an appropriate type of feedback is introduced. Results concerning transformations of nonlinear systems to linear systems are similar to results for transforming a linear PDE to a Kolmogorov equation.
Hunt, L. R.↗
Qualitative properties of large buckled states of spherical shells
A system of 6th-order quasi-linear Ordinary Differential Equations is analyzed to show the global existence of axisymmetrically buckled states. A surprising nodal property is obtained which shows that everywhere along a branch of solutions that bifurcates from a simple eigenvalue of the linearized equation, the number of simultaneously vanishing points of both shear resultant and circumferential bending moment resultant remains invariant, provided that a certain auxiliary condition is satisfied.
Shih, K. G.↗
Conformal BK equation at QCD Wilson-Fisher point
High-energy scattering in pQCD in the Regge limit is described by the evolution of Wilson lines governed by the BK equation. In the leading order, the BK equation is conformally invariant and the eigenfunctions of the linearized BFKL equation are powers. It is a common belief that at d ≠ 4 the BFKL equation is useless since unlike d = 4 case it cannot be solved by usual methods. However, we demonstrate that at critical Wilson-Fisher point of QCD the relevant part of NLO BK restores the conformal invariance so the solutions are again powers. As a check of our approach to high-energy amplitudes at the Wilson-Fisher point, we calculate the anomalous dimensions of twist-2 light-ray operators in the Regge limit j → 1.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
A high order accurate finite element algorithm for high Reynolds number flow prediction
A Galerkin-weighted residuals formulation is employed to establish an implicit finite element solution algorithm for generally nonlinear initial-boundary value problems. Solution accuracy, and convergence rate with discretization refinement, are quantized in several error norms, by a systematic study of numerical solutions to several nonlinear parabolic and a hyperbolic partial differential equation characteristic of the equations governing fluid flows. Solutions are generated using selective linear, quadratic and cubic basis functions. Richardson extrapolation is employed to generate a higher-order accurate solution to facilitate isolation of truncation error in all norms. Extension of the mathematical theory underlying accuracy and convergence concepts for linear elliptic equations is predicted for equations characteristic of laminar and turbulent fluid flows at nonmodest Reynolds number. The nondiagonal initial-value matrix structure introduced by the finite element theory is determined intrinsic to improved solution accuracy and convergence. A factored Jacobian iteration algorithm is derived and evaluated to yield a consequential reduction in both computer storage and execution CPU requirements while retaining solution accuracy.
Baker, A. J.↗
Instability and sound emission from a flow over a curved surface
The growth and decay of a wavepacket convecting in a boundary layer over a concave-convex surface is studied numerically using direct computations of the Navier-Stokes equations. The resulting sound radiation is computed using the linearized Euler equations with the pressure from the Navier-Stokes solution as a time-dependent boundary condition. It is shown that on the concave portion the amplitude of the wavepacket increases and its bandwidth broadens while on the convex portion some of the components in the packet are stabilized. The pressure field decays exponentially away from the surface and then algebraically exhibits a decay characteristic of acoustic waves in two dimensions. The far-field acoustic pressure exhibits a peak at a frequency corresponding to the inflow instability frequency.
Maestrello, L.↗
Preferred latitudes of the intertropical convergence zone
The latitude preference of the intertropical convergence zone (ITCZ) is examined on the basis of observations, theory, and a modeling analysis. Observations show that convection is enhanced at latitudes of about 4 deg to 10 deg relative to the equator, even in regions where the sea surface temperature (SST) is maximum on the equator. Both linear shallow-water theory and a moist primitive equation model suggest a new explanation for the off-equatorial latitude preference of the ITCZ that requires neither the existence of zonally propagating disturbances nor an off-equatorial maximum in SST. The shallow-water theory indicates that a finite-width, zonally oriented, midtropospheric heat source (i.e., an ITCZ) produces the greatest local low-level convergence when placed a finite distance away from the equator. This result suggests that an ITCZ is most likely to be supported via low-level convergence of moist energy when located at these "preferred" latitudes away from the equator. For a plausible range of heating widths and damping parameters, the theoretically predicted latitude is approximately equal to the observed position (s) of the ITCZ (s). Analysis with an axially symmetric, moist, primitive equation model indicates that when the latent heating field is allowed to be determined internally, a positive feedback develops between the midtropospheric latent heating and the low-level convergence, with the effect of enhancing the organization of convection at latitudes of about 4 deg to 12 deg. Numerical experiments show that (1) two peaks in convective precipitation develop straddling the equator when the SST maximum is located on the equator; (2) steady ITCZ-like structures form only when the SST maximum is located away from the equator; and (3) peaks in convection can develop away from the maximum in SST, with a particular preference for latitudes of about 4 deg to 12 deg, even in the ('cold') hemisphere without the SST maximum. The relationship between this mechanism and earlier theories is discussed, as are implications for the coupled ocean-atmosphere system and the roles played by midlevel latent heating and SST gradients in forcing the low-level atmospheric circulation in the tropics.
Waliser, Duane E.↗
Flutter and Divergence Analysis using the Generalized Aeroelastic Analysis Method
The Generalized Aeroelastic Analysis Method (GAAM) is applied to the analysis of three well-studied checkcases: restrained and unrestrained airfoil models, and a wing model. An eigenvalue iteration procedure is used for converging upon roots of the complex stability matrix. For the airfoil models, exact root loci are given which clearly illustrate the nature of the flutter and divergence instabilities. The singularities involved are enumerated, including an additional pole at the origin for the unrestrained airfoil case and the emergence of an additional pole on the positive real axis at the divergence speed for the restrained airfoil case. Inconsistencies and differences among published aeroelastic root loci and the new, exact results are discussed and resolved. The generalization of a Doublet Lattice Method computer code is described and the code is applied to the calculation of root loci for the wing model for incompressible and for subsonic flow conditions. The error introduced in the reduction of the singular integral equation underlying the unsteady lifting surface theory to a linear algebraic equation is discussed. Acknowledging this inherent error, the solutions of the algebraic equation by GAAM are termed 'exact.' The singularities of the problem are discussed and exponential series approximations used in the evaluation of the kernel function shown to introduce a dense collection of poles and zeroes on the negative real axis. Again, inconsistencies and differences among published aeroelastic root loci and the new 'exact' results are discussed and resolved. In all cases, aeroelastic flutter and divergence speeds and frequencies are in good agreement with published results. The GAAM solution procedure allows complete control over Mach number, velocity, density, and complex frequency. Thus all points on the computed root loci can be matched-point, consistent solutions without recourse to complex mode tracking logic or dataset interpolation, as in the k and p-k solution methods.
Edwards, John W.↗
Direct finite element equation solving algorithms
This paper presents and examines direct solution algorithms for the linear simultaneous equations that arise when finite element models represent an engineering system. It identifies the mathematical processing of four solution methods and assesses their data processing implications using concurrent processing.
Melosh, R. J.↗
Linearized frequency domain Landau-Lifshitz-Gilbert equation formulation
We present a general finite element linearized Landau-Lifshitz-Gilbert equation (LLGE) solver for magnetic systems under weak time-harmonic excitation field. The linearized LLGE is obtained by assuming a small deviation around the equilibrium state of the magnetic system. Inserting such expansion into LLGE and keeping only first order terms gives the linearized LLGE, which gives a frequency domain solution for the complex magnetization amplitudes under an external time-harmonic applied field of a given frequency. We solve the linear system with an iterative solver using generalized minimal residual method. We construct a preconditioner matrix to effectively solve the linear system. The validity, effectiveness, speed, and scalability of the linear solver are demonstrated via numerical examples.
36 MATERIALS SCIENCE↗
Fourth order difference methods for hyperbolic IBVP's
Fourth order difference approximations of initial-boundary value problems for hyperbolic partial differential equations are considered. We use the method of lines approach with both explicit and compact implicit difference operators in space. The explicit operator satisfies an energy estimate leading to strict stability. For the implicit operator we develop boundary conditions and give a complete proof of strong stability using the Laplace transform technique. We also present numerical experiments for the linear advection equation and Burgers' equation with discontinuities in the solution or in its derivative. The first equation is used for modeling contact discontinuities in fluid dynamics, the second one for modeling shocks and rarefaction waves. The time discretization is done with a third order Runge-Kutta TVD method. For solutions with discontinuities in the solution itself we add a filter based on second order viscosity. In case of the non-linear Burger's equation we use a flux splitting technique that results in an energy estimate for certain different approximations, in which case also an entropy condition is fulfilled. In particular we shall demonstrate that the unsplit conservative form produces a non-physical shock instead of the physically correct rarefaction wave. In the numerical experiments we compare our fourth order methods with a standard second order one and with a third order TVD-method. The results show that the fourth order methods are the only ones that give good results for all the considered test problems.
Gustafsson, Bertil↗
Absorbing boundary conditions for second-order hyperbolic equations
A uniform approach to construct absorbing artificial boundary conditions for second-order linear hyperbolic equations is proposed. The nonlocal boundary condition is given by a pseudodifferential operator that annihilates travelling waves. It is obtained through the dispersion relation of the differential equation by requiring that the initial-boundary value problem admits the wave solutions travelling in one direction only. Local approximation of this global boundary condition yields an nth-order differential operator. It is shown that the best approximations must be in the canonical forms which can be factorized into first-order operators. These boundary conditions are perfectly absorbing for wave packets propagating at certain group velocities. A hierarchy of absorbing boundary conditions is derived for transonic small perturbation equations of unsteady flows. These examples illustrate that the absorbing boundary conditions are easy to derive, and the effectiveness is demonstrated by the numerical experiments.
Jiang, Hong↗
Absorbing boundary conditions for second-order hyperbolic equations
A uniform approach to construct absorbing artificial boundary conditions for second-order linear hyperbolic equations is proposed. The nonlocal boundary condition is given by a pseudodifferential operator that annihilates travelling waves. It is obtained through the dispersion relation of the differential equation by requiring that the initial-boundary value problem admits the wave solutions travelling in one direction only. Local approximation of this global boundary condition yields an nth-order differential operator. It is shown that the best approximations must be in the canonical forms which can be factorized into first-order operators. These boundary conditions are perfectly absorbing for wave packets propagating at certain group velocities. A hierarchy of absorbing boundary conditions is derived for transonic small perturbation equations of unsteady flows. These examples illustrate that the absorbing boundary conditions are easy to derive, and the effectiveness is demonstrated by the numerical experiments.
Jiang, Hong↗
Numerical computation of exponential matrices using the Cayley-Hamilton theorem
A method for computing exponential matrices, which often arise naturally in the solution of systems of linear differential equations, is developed. An exponential matrix is generated as a linear combination of a finite number (equal to the matrix order) of matrices, the coefficients of which are scalar infinite sums. The method can be generalized to apply to any formal power series of matrices. Attention is focused upon the exponential function, and the matrix exponent is assumed tri-diagonal in form. In such cases, the terms in the coefficient infinite sums can be extracted, as recursion relations, from the characteristic polynomial of the matrix exponent. Two numerical examples are presented in some detail: (1) the three dimensional infinitesimal rotation rate matrix, which is skew symmetric, and (2) an N-dimensional tri-diagonal and symmetric finite difference matrix which arises in the numerical solution of the heat conduction partial differential equation. In the second example, the known eigenvalues and eigenvectors of the finite difference matrix permit an analytical solution for the exponential matrix, through the theory of diagonalization and similarity transformations, which is used for independent verification. The convergence properties of the scalar infinite summations are investigated for finite difference matrices of various orders up to ten, and it is found that the number of terms required for convergence increases slowly with the order of the matrix.
Walden, H.↗
Reduced-order modeling on a near-term quantum computer
Quantum computing is an advancing area of research in which computer hardware and algorithms are developed to take advantage of quantum mechanical phenomena. In recent studies, quantum algorithms have shown promise in solving linear systems of equations as well as systems of linear ordinary differential equations (ODEs) and partial differential equations (PDEs). Reducedorder modeling (ROM) algorithms for studying fluid dynamics have shown success in identifying linear operators that can describe flowfields, where dynamic mode decomposition (DMD) is a particularly useful method in which a linear operator is identified from data. In this work, DMD is reformulated as an optimization problem to propagate the state of the linearized dynamical system on a quantum computer. This reformulation was chosen as a means of facilitating implementation on a near-term quantum computer. Quadratic unconstrained binary optimization (QUBO), a technique for optimizing quadratic polynomials in binary variables, allows for quantum annealing algorithms to be applied. A quantum circuit model (quantum approximation optimization algorithm, QAOA) is utilized to obtain predictions of the state trajectories. Results are shown for the quantum-ROM predictions for flow over a 2D cylinder at Re = 220 and flow over a NACA0009 airfoil at Re = 500 and α = 15°. The quantum-ROM predictions are found to depend on the number of bits utilized for a fixed point representation and the truncation level of the DMD model. Comparisons with DMD predictions from a classical computer algorithm are made, as well as an analysis of the computational complexity and prospects for future, more fault-tolerant quantum computers.
97 MATHEMATICS AND COMPUTING↗
Linear feedback guidance
Determination of closed loop guidance function with linear, time-variable feedback studied by behavior of motion equations
FEEDBACK↗