Linear functional-differential equations with constant coefficients
Linear differential-difference equations with constant coefficients
LINEAR EQUATION↗
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.
An ascent existence and uniqueness proof for linear ordinary differential equations
Ascent existence and uniqueness proof for linear ordinary differential equations
Pellicciaro, E. J.↗
On the control of a linear functional- differential equation with quadratic cost
Riccati-like linear functional differential equation with quadratic cost, analyzing feedback control solution existence and uniqueness
Barnea, D. I.↗
Tensor-GMRES method for large sparse systems of nonlinear equations
This paper introduces a tensor-Krylov method, the tensor-GMRES method, for large sparse systems of nonlinear equations. This method is a coupling of tensor model formation and solution techniques for nonlinear equations with Krylov subspace projection techniques for unsymmetric systems of linear equations. Traditional tensor methods for nonlinear equations are based on a quadratic model of the nonlinear function, a standard linear model augmented by a simple second order term. These methods are shown to be significantly more efficient than standard methods both on nonsingular problems and on problems where the Jacobian matrix at the solution is singular. A major disadvantage of the traditional tensor methods is that the solution of the tensor model requires the factorization of the Jacobian matrix, which may not be suitable for problems where the Jacobian matrix is large and has a 'bad' sparsity structure for an efficient factorization. We overcome this difficulty by forming and solving the tensor model using an extension of a Newton-GMRES scheme. Like traditional tensor methods, we show that the new tensor method has significant computational advantages over the analogous Newton counterpart. Consistent with Krylov subspace based methods, the new tensor method does not depend on the factorization of the Jacobian matrix. As a matter of fact, the Jacobian matrix is never needed explicitly.
Feng, Dan↗
Runge-Kutta Methods for Linear Ordinary Differential Equations
Three new Runge-Kutta methods are presented for numerical integration of systems of linear inhomogeneous ordinary differential equations (ODES) with constant coefficients. Such ODEs arise in the numerical solution of the partial differential equations governing linear wave phenomena. The restriction to linear ODEs with constant coefficients reduces the number of conditions which the coefficients of the Runge-Kutta method must satisfy. This freedom is used to develop methods which are more efficient than conventional Runge-Kutta methods. A fourth-order method is presented which uses only two memory locations per dependent variable, while the classical fourth-order Runge-Kutta method uses three. This method is an excellent choice for simulations of linear wave phenomena if memory is a primary concern. In addition, fifth- and sixth-order methods are presented which require five and six stages, respectively, one fewer than their conventional counterparts, and are therefore more efficient. These methods are an excellent option for use with high-order spatial discretizations.
Zingg, David W.↗
Development of an unsteady aerodynamic analysis for finite-deflection subsonic cascades
An unsteady potential flow analysis, which accounts for the effects of blade geometry and steady turning, was developed to predict aerodynamic forces and moments associated with free vibration or flutter phenomena in the fan, compressor, or turbine stages of modern jet engines. Based on the assumption of small amplitude blade motions, the unsteady flow is governed by linear equations with variable coefficients which depend on the underlying steady low. These equations were approximated using difference expressions determined from an implicit least squares development and applicable on arbitrary grids. The resulting linear system of algebraic equations is block tridiagonal, which permits an efficient, direct (i.e., noniterative) solution. The solution procedure was extended to treat blades with rounded or blunt edges at incidence relative to the inlet flow.
Verdon, J. M.↗
Method of Conjugate Radii for Solving Linear and Nonlinear Systems
This paper describes a method to solve a system of N linear equations in N steps. A quadratic form is developed involving the sum of the squares of the residuals of the equations. Equating the quadratic form to a constant yields a surface which is an ellipsoid. For different constants, a family of similar ellipsoids can be generated. Starting at an arbitrary point an orthogonal basis is constructed and the center of the family of similar ellipsoids is found in this basis by a sequence of projections. The coordinates of the center in this basis are the solution of linear system of equations. A quadratic form in N variables requires N projections. That is, the current method is an exact method. It is shown that the sequence of projections is equivalent to a special case of the Gram-Schmidt orthogonalization process. The current method enjoys an advantage not shared by the classic Method of Conjugate Gradients. The current method can be extended to nonlinear systems without modification. For nonlinear equations the Method of Conjugate Gradients has to be augmented with a line-search procedure. Results for linear and nonlinear problems are presented.
Nachtsheim, Philip R.↗
Turbulent disruptions from the Strauss equations
Preliminary results are reported from application of a three-dimensional spectral method model to the solution of the Strauss (1976) reduced MHD equations. The investigation was focused on describing MHD turbulence in a current-carrying bounded magnetofluid. A cylindrical geometry with a square cross-section was considered, with the walls being rigid perfect conductors with free-slip boundary conditions. A uniform magnetic field and the electric current density both point in the z-direction. Initial conditions are specified which feature small amounts of random noise expressed as Fourier modes. Linearized equations are defined for tracing the movement to equlibrium conditions or other temporal development. The model is further refined with nonlinear equations to examine the effects of the appearance of disruptions. Comparisons are drawn between solutions obtained with linear and nonlinear equations, with an eye to the associated physical realities.
Dahlburg, J. P.↗
Simulation of an inductively coupled plasma with a two-dimensional Darwin particle-in-cell code
A two-dimensional particle-in-cell code for the simulation of low-frequency electromagnetic processes in laboratory plasmas has been developed. The code uses the Darwin method omitting the electromagnetic wave propagation. The Darwin method separates the electric field into solenoidal and irrotational parts. The solenoidal electric field is calculated with a new algorithm based on the equation for the electric field vorticity. The system of linear equations in the new algorithm is readily solved using a standard iterative method. The irrotational electric field is the electrostatic field calculated with the direct implicit algorithm. The code is verified by reproducing the two-stream instability, electron electromagnetic waves, and shear Alfvén waves. The code is applied to simulate an inductively coupled plasma with the driving current flowing around the plasma region. In this simulation, a ring of dense plasma forms at the initial stage but then the density becomes maximal in the center and decays monotonically toward the walls. The skin effect is in the transitional mode between local and non-local, and the electron velocity distribution function is non-Maxwellian.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Battery-Charge-State Model
Charge-state model for lead/acid batteries proposed as part of effort to make equivalent of fuel gage for battery-powered vehicles. Models based on equations that approximate observable characteristics of battery electrochemistry. Uses linear equations, easier to simulate on computer, and gives smooth transitions between charge, discharge, and recuperation.
Vivian, H. C.↗
Advanced algorithm for orbit computation
Computational and analytical techniques which simplify the solution of complex problems in orbit mechanics, Astrodynamics and Celestial Mechanics were developed. The major tool of the simplification is the substitution of transformations in place of numerical or analytical integrations. In this way the rather complicated equations of orbit mechanics might sometimes be reduced to linear equations representing harmonic oscillators with constant coefficients.
Szenbehely, V.↗
Linear/nonlinear behavior in unsteady transonic aerodynamics
The accurate calculation of the aerodynamic forces in unsteady transonic flow requires the solution of the nonlinear flow equations. The aeroelastician, on the other hand, seeks to treat his problems (flutter, for example) by means of linear equations whenever possible. He may do this, even when the underlying flow is nonlinear, if the perturbation forces are linear over some (perhaps small) range of unsteady amplitude of motion. This paper assesses the range of parameters over which linear behavior occurs. In particular calculations are made for an NACA 64A006 airfoil oscillating in pitch over a range of amplitudes, frequencies, and Mach numbers. The primary aerodynamic method used is the well known LTRAN2 code of Ballhaus and Goorjian that provides a finite-difference solution to the low frequency, small disturbance, two-dimensional potential flow equation. Comparisons are made with linear subsonic theory, local linearization, and, for steady flow, with the full potential equation code of Bauer, Garabedian, and Korn.
Dowell, E. H.↗
The condition of regular degeneration for singularly perturbed linear differential- difference equations.
Initial value problem for singularly perturbed differential-difference equation which is scalar, linear and has constant coefficients
DIFFERENTIAL EQUATION↗
Simplification in the computation of the sensitivity functions for constant coefficient linear systems
Constant coefficient linear systems sensitivity functions computational procedure using nth order differential equations linear transformations
Denery, D. G.↗
Similarity analysis of differential equations by Lie group.
Methods for transforming partial differential equations into forms more suitable for analysis and solution are investigated. The idea of Lie's infinitesimal contact transformation group is introduced to develop a systematic method which involves mostly algebraic manipulations. A thorough presentation of the application of this general method to the problem of similarity analysis in a broader sense - namely, the similarity between partial and ordinary differential equations, boundary value and initial value problems, and nonlinear and linear equations - is given with new and very general methods evolved for deriving the possible groups of transformations.
Na, T. Y.↗
The causation and sensitivity of the northern winter planetary waves
High-resolution calculations of the wintertime stationary atmospheric response to planetary scale topographic and thermal forcings have been made using a numerical model that solves the spherical primitive equations linearized about observed zonal wind and temperature fields. The model equations are presented and their numerical implementation is discussed. The model inputs and parameters are presented and the model's response to relativistic forcing is analyzed. In middle latitudes the response to topographic forcing strongly dominates the response to thermal forcing. In the midlatitude troposphere, the topographic response is insensitive to changes in the zonal wind. In the stratosphere, both the variance and mean are dominated by the topographic response. In general, the sensitivity of the planetary waves to changes in the zonal wind is found to be much smaller than other recent calculations have indicated.
Jacqmin, D.↗
Towards enhancing and delaying disturbances in free shear flows
The family of shear flows comprising the jet, wake, and the mixing layer are subjected to perturbations in an inviscid incompressible fluid. By modeling the basic mean flows as parallel with piecewise linear variations for the velocities, complete and general solutions to the linearized equations of motion can be obtained in closed form as functions of all space variables and time when posed as an initial value problem. The results show that there is a continuous as well as the discrete spectrum that is more familiar in stability theory and therefore there can be both algebraic and exponential growth of disturbances in time. These bases make it feasible to consider control of such flows. To this end, the possibility of enhancing the disturbances in the mixing layer and delaying the onset in the jet and wake is investigated. It is found that growth of perturbations can be delayed to a considerable degree for the jet and the wake but, by comparison, cannot be enhanced in the mixing layer. By using moving coordinates, a method for demonstrating the predominant early and long time behavior of disturbances in these flows is given for continuous velocity profiles. It is shown that the early time transients are always algebraic whereas the asymptotic limit is that of an exponential normal mode. Numerical treatment of the new governing equations confirm the conclusions reached by use of the piecewise linear basic models. Although not pursued here, feedback mechanisms designed for control of the flow could be devised using the results of this work.
Criminale, W. O.↗
LuGo: An enhanced quantum phase estimation implementation
Quantum Phase Estimation (QPE) is a cardinal algorithm in quantum computing that plays a crucial role in various applications, including cryptography, molecular simulation, and solving systems of linear equations. However, the standard implementation of QPE faces challenges related to time complexity and circuit depth, which limit its practicality for large-scale computations. We introduce LuGo, a novel framework designed to enhance the performance of QPE by reducing circuit duplication, as well as using parallelization techniques to achieve faster generation of the QPE circuit and gate reduction. We validate the effectiveness of our framework by generating quantum linear solver circuits, which require both QPE and inverse QPE, to solve linear systems of equations. LuGo achieves significant improvements in both computational efficiency and hardware requirements without compromising on accuracy. Compared to a standard QPE implementation, LuGo reduces time consumption to generate a circuit that solves a 2 6 × 2 6 system matrix by a factor of 50.68 and over 31× reduction of quantum gates and circuit depth, with no fidelity loss on an ideal quantum simulator. Furthermore, we demonstrated the versatility and scalability of LuGo enabled HHL algorithm by simulating a canonical Hele-Shaw fluid problem using a quantum simulator. With these advantages, LuGo paves the way for more efficient implementations of QPE, enabling broader applications across several quantum computing domains.
Quantum algorithm↗