Engineering Papers⌕ Search

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.

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Random circuit block-encoded matrix and a proposal of quantum LINPACK benchmark

The LINPACK benchmark reports the performance of a computer for solving a system of linear equations with dense random matrices. Although this task was not designed with a real application directly in mind, the LINPACK benchmark has been used to define the list of TOP500 supercomputers since the debut of the list in 1993. We propose that a similar benchmark, called the quantum LINPACK benchmark, could be used to measure the whole machine performance of quantum computers. The success of the quantum LINPACK benchmark should be viewed as the minimal requirement for a quantum computer to perform a useful task of solving linear algebra problems, such as linear systems of equations. We propose an input model called the Random Circuit Block-Encoded Matrix (RACBEM), which is a proper generalization of a dense random matrix in the quantum setting. The RACBEM model is efficient to be implemented on a quantum computer and can be designed to optimally adapt to any given quantum architecture, with relying on a black-box quantum compiler. Besides solving linear systems, the RACBEM model can be used to perform a variety of linear algebra tasks relevant to many physical applications, such as computing spectral measures, time series generated by a Hamiltonian simulation, and thermal averages of the energy. We implement these linear algebra operations on IBM Q quantum devices as well as quantum virtual machines, and demonstrate their performance in solving scientific computing problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1992

Spacecraft attitude and velocity control system

A spacecraft attitude and/or velocity control system includes a controller which responds to at least attitude errors to produce command signals representing a force vector F and a torque vector T, each having three orthogonal components, which represent the forces and torques which are to be generated by the thrusters. The thrusters may include magnetic torquer or reaction wheels. Six difference equations are generated, three having the form ##EQU1## where a.sub.j is the maximum torque which the j.sup.th thruster can produce, b.sub.j is the maximum force which the j.sup.th thruster can produce, and .alpha..sub.j is a variable representing the throttling factor of the j.sup.th thruster, which may range from zero to unity. The six equations are summed to produce a single scalar equation relating variables .alpha..sub.j to a performance index Z: ##EQU2## Those values of .alpha. which maximize the value of Z are determined by a method for solving linear equations, such as a linear programming method. The Simplex method may be used. The values of .alpha..sub.j are applied to control the corresponding thrusters.

Paluszek, Michael A.↗
NASA NTRS1987

The Jeffcott equations in nonlinear rotordynamics

The Jeffcott equations are a system of coupled differential equations representing the behavior of a rotating shaft. This is a simple model which allows investigation of the basic dynamic behavior of rotating machinery. Nolinearities can be introduced by taking into consideration deadband, side force, and rubbing, among others. The properties of the solutions of the Jeffcott equations with deadband are studied. In particular, it is shown how bounds for the solution of these equations can be obtained from bounds for the solutions of the linearized equations. By studying the behavior of the Fourier transforms of the solutions, we are also able to predict the onset of destructive vibrations. These conclusions are verified by means of numerical solutions of the equations, and of power spectrum density (PSD) plots. This study offers insight into a possible detection method to determine pump stability margins during flight and hot fire tests, and was motivated by the need to explain a phenomenon observed in the development phase of the cryogenic pumps of the Space Shuttle, during hot fire ground testing; namely, the appearance of vibrations at frequencies that could not be accounted for by means of linear models.

Zalik, R. A.↗
NASA NTRS1989

Performance optimization for rotors in hover and axial flight

Performance optimization for rotors in hover and axial flight is a topic of continuing importance to rotorcraft designers. The aim of this Phase 1 effort has been to demonstrate that a linear optimization algorithm could be coupled to an existing influence coefficient hover performance code. This code, dubbed EHPIC (Evaluation of Hover Performance using Influence Coefficients), uses a quasi-linear wake relaxation to solve for the rotor performance. The coupling was accomplished by expanding of the matrix of linearized influence coefficients in EHPIC to accommodate design variables and deriving new coefficients for linearized equations governing perturbations in power and thrust. These coefficients formed the input to a linear optimization analysis, which used the flow tangency conditions on the blade and in the wake to impose equality constraints on the expanded system of equations; user-specified inequality contraints were also employed to bound the changes in the design. It was found that this locally linearized analysis could be invoked to predict a design change that would produce a reduction in the power required by the rotor at constant thrust. Thus, an efficient search for improved versions of the baseline design can be carried out while retaining the accuracy inherent in a free wake/lifting surface performance analysis.

Quackenbush, T. R.↗
NASA NTRS1989

Optimization of rotor performance in hover and axial flight using a free wake analysis

Performance optimization for rotors in hover and axial flight is a topic of continuing importance to rotorcraft designers. The aim of this effort was to demonstrate that a numerical performance optimization algorithm could be coupled to an existing free wake hover code. This code, dubbed EHPIC (Evaluation of Hover Performacne using Influence Coefficients), uses a quasi-linear wake relaxation to solve for the rotor performance. The coupling was accomplished by expanding of the matrix of linearized influence coefficients in EHPIC to accommodate design variables and by deriving new coefficients for linearized equations governing perturbations in power and thrust. These coefficients formed the input to a linear optimization analysis, which used the flow tangency conditions on the blade and in the wake to impose equality constraints on the expanded system of equations; user-specified inequality constraints were also employed to bound the changes in the design. It was found that this locally linearized analysis could be invoked to predict a design change that would produce a reduction in the power required by the rotor at constant thrust. Thus, an efficient search for improved versions of the baseline design can be carried out while retaining the accuracy inherent in a free wake/lifting surface performance analysis. A variety of sample problems were undertaken to demonstrate the success of this approach in reducing the power required at a specified thrust for several representative rotor configurations in hover and axial flight.

Quackenbush, T. R.↗
NASA NTRS1975

Flapping response characteristics of hingeless rotor blades by a gereralized harmonic balance method

Linearized equations of motion for the flapping response of flexible rotor blades in forward flight are derived in terms of generalized coordinates. The equations are solved using a matrix form of the method of linear harmonic balance, yielding response derivatives for each harmonic of the blade deformations and of the hub forces and moments. Numerical results and approximate closed-form expressions for rotor derivatives are used to illustrate the relationships between rotor parameters, modeling assumptions, and rotor response characteristics. Finally, basic hingeless rotor response derivatives are presented in tabular and graphical form for a wide range of configuration parameters and operating conditions.

Peters, D. A.↗
NASA NTRS1977

Application of the TRPLT1 element to large amplitude free vibrations of plates

A finite element formulation is developed for analyzing large amplitude free flexural vibrations of thin plates in NASTRAN. Stress distributions in the plate, in addition to deflection shapes and nonlinear frequencies are determined. Linearized equations of motion governing large amplitude oscillations of plates and a linearized geometrical stiffness matrix are presented. The solution procedure and convergence characteristics are discussed. The quasi-linear geometrical stiffness matrix for an eighteen degree-of-freedom higher order triangular plate element is evaluated by using a seven-point numerical integration. Nonlinear frequencies for square, rectangular, rhombic, and isosceles triangular plates, with edges simply supported or clamped, are compared with earlier solutions. The present formulation is found to give results entirely adequate for engineering purposes.

Mei, C.↗
NASA NTRS1988

Rotordynamic analysis of a bearing tester

The properties of the solutions of a system of four coupled nonlinear differential equations that model the behavior of the rotating shaft of a bearing tester are studied. In particular, it is shown how the bounds for the rotations of these equations can be obtained from bounds for the solutions of the linearized equations. By studying the behavior of the Fourier transforms of the solution, the approach to the stability boundary can also be predicted. These conclusions are verified by means of numerical solutions of the equations, and of power spectrum density (PSD) plots.

Zalik, Richard A.↗
NASA NTRS2002

Marangoni Effects of a Drop in an Extensional Flow: The Role of Surfactant Physical Chemistry

While the changes in stresses caused by surfactant adsorption on non-deforming interfaces have been fairly well established, prior to this work, there were few studies addressing how surfactants alter stresses on strongly deforming interfaces. We chose the model problem of a drop in a uniaxial extensional flow to study these stress conditions To model surfactant effects at fluid interfaces, a proper description of the dependence of the surface tension on surface concentration, the surface equation of state, is required. We have adopted a surface equation of state that accounts for the maximum coverage limit; that is, because surfactants have a finite cross sectional area, there is an upper bound to the amount of surfactant that can adsorb in a monolayer. The surface tension reduces strongly only when this maximum coverage is approached. Since the Marangoni stresses go as the derivative of the surface equation of state times the surface concentration gradient, the non-linear equation of state determines both the effect of surfactants in the normal stress jump, (which is balanced by the product of the mean curvature of the interface times the surface tension), and the tangential stress jump, which is balanced by Marangoni stresses. First, the effects of surface coverage and intermolecular interactions among surfactants which drive aggregation of surfactants in the interface were studied. (see Pawar and Stebe, Physics of Fluids).

Stebe, Kathleen J.↗
NASA NTRS1990

Attitude control/momentum management of the Space Station Freedom for large angle torque-equilibrium-attitude configurations

An attitude-control and momentum-management (ACMM) system for the Space Station in a large-angle torque-equilibrium-attitude (TEA) configuration is developed analytically and demonstrated by means of numerical simulations. The equations of motion for a rigid-body Space Station model are outlined; linearized equations for an arbitrary TEA (resulting from misalignment of control and body axes) are derived; the general requirements for an ACMM are summarized; and a pole-placement linear-quadratic regulator solution based on scheduled gains is proposed. Results are presented in graphs for (1) simulations based on configuration MB3 (showing the importance of accounting for the cross-inertia terms in the TEA estimate) and (2) simulations of a stepwise change from configuration MB3 to the 'assembly complete' stage over 130 orbits (indicating that the present ACCM scheme maintains sufficient control over slowly varying Space Station dynamics).

Parlos, Alexander G.↗
NASA NTRS1975

A formalism for cosmic ray propagation studies

The continuity equation for cosmic ray propagation is used to derive a set of linear equations interrelating the fluxes of multiply charged nuclei as observed at any particular part of the galaxy. The derivation leads to model independent definitions for cosmic ray storage time, mean density of target nuclei and effective mass traversed. The set of equations form a common framework for comparisons of theories and observations. As an illustration, it is shown that there exists a large class of propagation models which give the same result as the exponential path length model. The formalism is shown to accommodate dynamic as well as equilibrium models of production and propagation.

Golden, R. L.↗
NASA NTRS1982

Iterative algorithms for large sparse linear systems on parallel computers

Algorithms for assembling in parallel the sparse system of linear equations that result from finite difference or finite element discretizations of elliptic partial differential equations, such as those that arise in structural engineering are developed. Parallel linear stationary iterative algorithms and parallel preconditioned conjugate gradient algorithms are developed for solving these systems. In addition, a model for comparing parallel algorithms on array architectures is developed and results of this model for the algorithms are given.

Adams, L. M.↗
NASA NTRS1983

Character and stability of axisymmetric thermal convection in spheres and spherical shells

The influence of shell size and mode of heating on the behavior and stability of axisymmetric, infinite Prandtl number convection in a spherical geometry is studied. Heating from within and below features convection onset governed by a self-adjoint system of equations and boundary conditions. For heating only from within or from below, linearized equations and boundary conditions are non-self-adjoint. Identification of the parameter which initiates the departure from self-adjointness, together with the properties of the self-adjoint solution, provide a basis for calculating the heat transfer characteristics of the non-self-adjoint situations. The investigations are an effort to develop a model for heat transfer in planetary interiors. Further development of the technique by modifying the Galerkin method by the introduction of diagonal mode truncation is suggested to permit the consideration of higher values of the Rayleigh numbers, i.e., those more commensurate with terrestrial planet mantles.

Zebib, A.↗
NASA NTRS1959

Loads and Deformations of Buckled Rectangular Plates

The nonlinear large-deflection equations of von Karman for plates are converted into a set of linear equations by expanding the displacements Into a power series in terms of an arbitrary parameter. The postbuckling behavior of simply supported rectangular plates subjected to longitudinal compression and subject to a uniform temperature rise is investigated in detail by solving the first few of the equations. Experimental data are presented for the compression problem. Comparisons are made for total shortening and local strains and deflections which indicate good agreement between experimental and theoretical results.

Stein, Manuel↗
NASA NTRS1975

Almost periodic solutions to difference equations

The theory of Massera and Schaeffer relating the existence of unique almost periodic solutions of an inhomogeneous linear equation to an exponential dichotomy for the homogeneous equation was completely extended to discretizations by a strongly stable difference scheme. In addition it is shown that the almost periodic sequence solution will converge to the differential equation solution. The preceding theory was applied to a class of exponentially stable partial differential equations to which one can apply the Hille-Yoshida theorem. It is possible to prove the existence of unique almost periodic solutions of the inhomogeneous equation (which can be approximated by almost periodic sequences) which are the solutions to appropriate discretizations. Two methods of discretizations are discussed: the strongly stable scheme and the Lax-Wendroff scheme.

Bayliss, A.↗
NASA NTRS1976

Oscillations and stability of numerical solutions of the heat conduction equation

The mathematical model and results of numerical solutions are given for the one dimensional problem when the linear equations are written in a rectangular coordinate system. All the computations are easily realizable for two and three dimensional problems when the equations are written in any coordinate system. Explicit and implicit schemes are shown in tabular form for stability and oscillations criteria; the initial temperature distribution is considered uniform.

Kozdoba, L. A.↗