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 289 records · Page 16

Dynamic model for an instrument pointing system

The instrument pointing system concept discussed in the present paper was developed for Spacelab missions. The system is mounted on soft shockmounts to minimize Shuttle disturbances (such as man motion, thruster firings, etc.) and to reduce the effect of large center of mass offsets between the instrument and the pointing system gimbal axes. The instrument pointing system is soft mounted to the Spacelab pallet located in the Orbiter payload bay. The nonlinear model, incorporating the dynamic equations of motion, the controller equations and the definition of sensor dynamics and gimbal bearing friction, is given, along with a system of linearized equations of motion of the Space Shuttle Orbiter and instrument pointing system (linearized model).

Howell, J. T.↗

Linear systems of equations solved using mathematical algorithms

New mathematical algorithm solves linear systems of equations, AX equals B, and preserves the integer properties of the coefficients. The algorithms presented can also be used for the efficient evaluation of determinates and their leading minors.

Bareiss, E. H.↗

Theoretical Determination of Lifetime of Compressed Plates at Elevated Temperatures

A method for the theoretical determination of the lifetime of com- pressed plates at elevated temperatures is presented. In this approach, linearized equations are used throughout with the assumption that the plate material is a standard linear solid. The critical time (lifetime) is determined by reducing the time-dependent behavior to the time- independent response of purely elastic buckling. Theoretically predicted lifetimes of 2024-T3 (formerly 24S-T3) aluminum-alloy plates at 450 F are compared with experimental values obtained in previous work.

Herrmann, George↗

Verification of a quasi-linear model for gyrokinetic turbulent transport

The verification and calibration of a new quasi-linear transport model with a large database of gyrokinetic turbulence simulations is presented in this paper. In a previous paper (Staebler et al 2020 Plasma Phys. Control. Fusion 63 015013), a model for the saturated spectrum of electric potential fluctuations was developed based on the properties of the non-linear 3D spectrum. In this paper, a modification to the overall multiplicative factor in this model is found to be necessary to improve the fit to scans of the temperature and density gradients. The error in the fit of the quasi-linear fluxes of electron and ion energy fluxes is significantly better than for previous saturation models. The spectral shift model for the impact of equilibrium E × B velocity shear (Staebler et al 2013 Phys. Rev. Lett. 110 055003) and the zonal flow mixing model for electron-scale turbulence (Staebler et al 2016 Phys. Plasmas 23 062518) are both revised to be compatible with this new model. Furthermore, the models for the loss of bounce averaging and electron collisions in the TGLF reduced linear equations (Staebler et al 2005 Phys. Plasmas 12 102508) are also changed to improve the linear eigenmodes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Unsteady thermocapillary migration of bubbles

Upon the introduction of a gas bubble into a liquid possessing a uniform thermal gradient, an unsteady thermo-capillary flow begins. Ultimately, the bubble attains a constant velocity. This theoretical analysis focuses upon the transient period for a bubble in a microgravity environment and is restricted to situations wherein the flow is sufficiently slow such that inertial terms in the Navier-Stokes equation and convective terms in the energy equation may be safely neglected (i.e., both Reynolds and Marangoni numbers are small). The resulting linear equations were solved analytically in the Laplace domain with the Prandtl number of the liquid as a parameter; inversion was accomplished numerically using a standard IMSL routine. In the asymptotic long-time limit, the theory agrees with the steady-state theory of Young, Goldstein, and Block. The theory predicts that more than 90 percent of the terminal steady velocity is achieved when the smallest dimensionless time, i.e., the one based upon the largest time scale-viscous or thermal-equals unity.

Dill, Loren H.↗

An efficient formulation and implementation of the analytic energy gradient method to the single and double excitation coupled-cluster wave function - Application to Cl2O2

The analytic energy gradient for the single and double excitation coupled-cluster (CCSD) wave function has been reformulated and implemented in a new set of programs. The reformulated set of gradient equations have a smaller computational cost than any previously published. The iterative solution of the linear equations and the construction of the effective density matrices are fully vectorized, being based on matrix multiplications. The new method has been used to investigate the Cl2O2 molecule, which has recently been postulated as an important intermediate in the destruction of ozone in the stratosphere. In addition to reporting computational timings, the CCSD equilibrium geometries, harmonic vibrational frequencies, infrared intensities, and relative energetics of three isomers of Cl2O2 are presented.

Rendell, Alistair P.↗

On the prediction of velocity fields from redshift space galaxy samples

We present a new method for recovering the underlying velocity field from an observed distribution of galaxies in redshift space. The method is based on a kinematic Zel'dovich relation between the velocity and density fields in redshift space. This relation is expressed in a differential equation slightly modified from the usual Poisson equation and which depends nontrivially on Beta identically equal to Omega(exp 0.6)/b. The linear equation can be readily solved by standard techniques of separation of variables by means of spherical harmonics. One can also include a term describing the 'rocket effect' discussed by Kaiser (1987). From this redshift space information alone, one can generate a prediction of the peculiar velocity field for each harmonic (l, m) as a function of distance. We note that for the quadrupole and higher order moments, the equation is a boundary value problem with solutions dependent on both the interior and exterior mass distribution. However, for a shell at distance r, the dipole, as well as the monopole, of the velocity field in the Local Group frame is fully determined by the interior mass distribution. This implies that the shear of the measured velocity field, when fitted to a dipole distortion, should be aligned and consistent with the gravity field inferred from the well determined local galaxy distribution. As a preliminary application we compute the velocity dipole of distant shells as predicted from the 1.2 Jy IRAS survey compared to the measured velocity dipole on shells, as inferred from a recent POTENT analysis. The coherence between the two fields is good, yielding a best estimate of Beta = 0.6 +/- 0.2.

Nusser, Adi↗

Interpolation using surface splines.

A surface spline is a mathematical tool for interpolating a function of two variables. It is based upon the small deflection equation of an infinite plate. The surface spline depends upon the solution of a system of linear equations, and thus, will ordinarily require the use of a digital computer. The closed form solution involves no functions more complicated than logarithms, and is easily coded. Several modifications which can be incorporated are discussed.

Harder, R. L.↗

A study of the control problem of the shoot side environment delivery system of a closed crop growth research chamber

The details of our initial study of the control problem of the crop shoot environment of a hypothetical closed crop growth research chamber (CGRC) are presented in this report. The configuration of the CGRC is hypothetical because neither a physical subject nor a design existed at the time the study began, a circumstance which is typical of large scale systems control studies. The basis of the control study is a mathematical model which was judged to adequately mimic the relevant dynamics of the system components considered necessary to provide acceptable realism in the representation. Control of pressure, temperature, and flow rate of the crop shoot environment, along with its oxygen, carbon dioxide, and water concentration is addressed. To account for mass exchange, the group of plants is represented in the model by a source of oxygen, a source of water vapor, and a sink for carbon dioxide. In terms of the thermal energy exchange, the group of plants is represented by a surface with an appropriate temperature. Most of the primitive equations about an experimental operating condition and a state variable representation which was extracted from the linearized equations are presented. Next, we present the results of a real Jordan decomposition and the repositioning of an undesirable eigenvalue via full state feedback. The state variable representation of the modeling system is of the nineteenth order and reflects the eleven control variables and eight system disturbances. Five real eigenvalues are very near zero, with one at zero, three having small magnitude positive values, and one having a small magnitude negative value. A Singular Value Decomposition analysis indicates that these non-zero eigenvalues are not results of numerical error.

Blackwell, C. C.↗

Thermohydrodynamic Analysis of Cryogenic Liquid Turbulent Flow Fluid Film Bearings

Computational programs developed for the thermal analysis of tilting and flexure-pad hybrid bearings, and the unsteady flow and transient response of a point mass rotor supported on fluid film bearings are described. The motion of a cryogenic liquid on the thin film annular region of a fluid film bearing is described by a set of mass and momentum conservation, and energy transport equations for the turbulent bulk-flow velocities and pressure, and accompanied by thermophysical state equations for evaluation of the fluid material properties. Zeroth-order equations describe the fluid flow field for a journal static equilibrium position, while first-order (linear) equations govern the fluid flow for small amplitude-journal center translational motions. Solution to the zeroth-order flow field equations provides the bearing flow rate, load capacity, drag torque and temperature rise. Solution to the first-order equations determines the rotordynamic force coefficients due to journal radial motions.

Luis San Andres↗

Preconditioned conjugate gradient methods for the compressible Navier-Stokes equations

The compressible Navier-Stokes equations are solved for a variety of two-dimensional inviscid and viscous problems by preconditioned conjugate gradient-like algorithms. Roe's flux difference splitting technique is used to discretize the inviscid fluxes. The viscous terms are discretized by using central differences. An algebraic turbulence model is also incorporated. The system of linear equations which arises out of the linearization of a fully implicit scheme is solved iteratively by the well known methods of GMRES (Generalized Minimum Residual technique) and Chebyschev iteration. Incomplete LU factorization and block diagonal factorization are used as preconditioners. The resulting algorithm is competitive with the best current schemes, but has wide applications in parallel computing and unstructured mesh computations.

Venkatakrishnan, V.↗

Dynamics of a spinning space station with a counterweight connected by multiple cables.

A unique and relatively simple approach is presented for obtaining the linearized equations of motion. A conservative system consisting of two rigid bodies connected by any number of massless cables with linear axial stiffness is assumed. The cable forces and torques are expanded in a Taylor series about the equilibrium values of the system coordinates which results in a cable stiffness matrix. The method of obtaining the equilibrium values of the coordinates is discussed and results are presented. The range of validity of the linear model is determined by comparing results with a digital simulation of the nonlinear system.

Nixon, D. D.↗

Bulk-Flow Analysis of Hybrid Thrust Bearings for Advanced Cryogenic Turbopumps

A bulk-flow analysis and computer program for prediction of the static load performance and dynamic force coefficients of angled injection, orifice-compensated hydrostatic/hydrodynamic thrust bearings have been completed. The product of the research is an efficient computational tool for the design of high-speed thrust bearings for cryogenic fluid turbopumps. The study addresses the needs of a growing technology that requires of reliable fluid film bearings to provide the maximum operating life with optimum controllable rotordynamic characteristics at the lowest cost. The motion of a cryogenic fluid on the thin film lands of a thrust bearing is governed by a set of bulk-flow mass and momentum conservation and energy transport equations. Mass flow conservation and a simple model for momentum transport within the hydrostatic bearing recesses are also accounted for. The bulk-flow model includes flow turbulence with fluid inertia advection, Coriolis and centrifugal acceleration effects on the bearing recesses and film lands. The cryogenic fluid properties are obtained from realistic thermophysical equations of state. Turbulent bulk-flow shear parameters are based on Hirs' model with Moody's friction factor equations allowing a simple simulation for machined bearing surface roughness. A perturbation analysis leads to zeroth-order nonlinear equations governing the fluid flow for the thrust bearing operating at a static equilibrium position, and first-order linear equations describing the perturbed fluid flow for small amplitude shaft motions in the axial direction. Numerical solution to the zeroth-order flow field equations renders the bearing flow rate, thrust load, drag torque and power dissipation. Solution to the first-order equations determines the axial stiffness, damping and inertia force coefficients. The computational method uses well established algorithms and generic subprograms available from prior developments. The Fortran9O computer program hydrothrust runs on a Windows 95/NT personal computer. The program, help files and examples are licensed by Texas A&M University Technology License Office. The study of the static and dynamic performance of two hydrostatic/hydrodynamic bearings demonstrates the importance of centrifugal and advection fluid inertia effects for operation at high rotational speeds. The first example considers a conceptual hydrostatic thrust bearing for an advanced liquid hydrogen turbopump operating at 170,000 rpm. The large axial stiffness and damping coefficients of the bearing should provide accurate control and axial positioning of the turbopump and also allow for unshrouded impellers, therefore increasing the overall pump efficiency. The second bearing uses a refrigerant R134a, and its application in oil-free air conditioning compressors is of great technological importance and commercial value. The computed predictions reveal that the LH2 bearing load capacity and flow rate increase with the recess pressure (i.e. increasing orifice diameters). The bearing axial stiffness has a maximum for a recess pressure rati of approx. 0.55. while the axial damping coefficient decreases as the recess pressure ratio increases. The computer results from three flow models are compared. These models are a) inertialess, b) fluid inertia at recess edges only, and c) full fluid inertia at both recess edges and film lands. The full inertia model shows the lowest flow rates, axial load capacity and stiffness coefficient but on the other hand renders the largest damping coefficients and inertia coefficients. The most important findings are related to the reduction of the outflow through the inner radius and the appearance of subambient pressures. The performance of the refrigerant hybrid thrust bearing is evaluated at two operating speeds and pressure drops. The computed results are presented in dimensionless form to evidence consistent trends in the bearing performance characteristics. As the applied axial load increases, the bearing film thickness and flow rate decrease while the recess pressure increases. The axial stiffness coefficient shows a maximum for a certain intermediate load while the damping coefficient steadily increases. The computed results evidence the paramount of centrifugal fluid inertia at low recess pressures (i.e. low loads), and where there is actually an inflow through the bearing inner diameter, accompanied by subambient pressures just downstream of the bearing recess edge. These results are solely due to centrifugal fluid inertia and advection transport effects. Recommendations include the extension of the computer program to handle flexure pivot tilting pad hybrid bearings and the ability to calculate moment coefficients for shaft angular misalignments.

SanAndres, Luis↗

The Elasto-Plastic Stability of Plates

This article explains results developed from the following research: 'The Stability of Plates and Shells beyond the Elastic Limit.' A significant improvement is found in the derivation of the relations between the stress factors and the strains resulting from the instability of plates and shells. In a strict analysis, the problem reduces to the solution of two simultaneous nonlinear partial differential equations of the fourth order in the deflection and stress function, and in the approximate analysis to a single linear equation of the Bryan type. Solutions are given for the special cases of a rectangular plate buckling into a cylindrical form, and of an arbitrarily shaped plate under uniform compression. These solutions indicate that the accuracy obtained by the approximate method is satisfactory.

Ilyushin, A. A.↗

Least-squares finite element solution of 3D incompressible Navier-Stokes problems

Although significant progress has been made in the finite element solution of incompressible viscous flow problems. Development of more efficient methods is still needed before large-scale computation of 3D problems becomes feasible. This paper presents such a development. The most popular finite element method for the solution of incompressible Navier-Stokes equations is the classic Galerkin mixed method based on the velocity-pressure formulation. The mixed method requires the use of different elements to interpolate the velocity and the pressure in order to satisfy the Ladyzhenskaya-Babuska-Brezzi (LBB) condition for the existence of the solution. On the other hand, due to the lack of symmetry and positive definiteness of the linear equations arising from the mixed method, iterative methods for the solution of linear systems have been hard to come by. Therefore, direct Gaussian elimination has been considered the only viable method for solving the systems. But, for three-dimensional problems, the computer resources required by a direct method become prohibitively large. In order to overcome these difficulties, a least-squares finite element method (LSFEM) has been developed. This method is based on the first-order velocity-pressure-vorticity formulation. In this paper the LSFEM is extended for the solution of three-dimensional incompressible Navier-Stokes equations written in the following first-order quasi-linear velocity-pressure-vorticity formulation.

Jiang, Bo-Nan↗