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 577 records · Page 32

Use of the maximum entropy method to retrieve the vertical atmospheric ozone profile and predict atmospheric ozone content

A method is developed for prediction of ozone levels in planetary atmospheres. This method is formulated in terms of error covariance matrices, and is associated with both direct measurements, a priori first guess profiles, and a weighting function matrix. This is described by the following linearized equation: y = A(matrix) x X + eta, where A is the weighting matrix and eta is noise. The problems to this approach are: (1) the A matrix is near singularity; (2) the number of unknowns in the profile exceeds the number of data points, therefore, the solution may not be unique; and (3) even if a unique solution exists, eta may cause the solution to be ill conditioned.

Turner, B. Curtis↗

A theoretical approach for analyzing the restabilization of wakes

Recently reported experimental results demonstrate that restabilization of the low-Reynolds-number flow past a circular cylinder can be achieved by the placement of a smaller cylinder in the wake of the first at particular locations. Traditional numerical procedures for modeling such phenomena are computationally expensive. An approach is presented here in which the properties of the adjoint solutions to the linearized equations of motion are exploited to map quickly the best positions for the small cylinder's placement. Comparisons with experiment and previous computations are favorable. The approach is shown to be applicable to general flows, illustrating how strongly control mechanisms that involve sources of momentum couple to unstable (or stable) modes of the system.

Hill, D. C.↗

Boundary-layer receptivity due to a wall suction and control of Tollmien-Schlichting waves

A numerical study of the generation of Tollmien-Schlichting (T-S) waves due to the interaction between a small free-stream disturbance and a small localized suction slot on an otherwise flat surface was carried out using finite difference methods. The nonlinear steady flow is of the viscous-inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier-Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of T-S waves generated by the interaction between the free-stream disturbance and the suction slot, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the T-S waves and the demonstration of the possible active control of the growth of T-S waves.

Bodonyi, R. J.↗

Adaptive attitude stability and control for Space Station/Orbiter berthing operations

The fully coupled, systematically linearized equations of motion are used for attitude control system design of an early Space Station Freedom flight configuration. These equations are valid around any estimated spacecraft equilibrium point. A full state feedback controller, incorporating gain scheduled adaptation of the attitude gains, is developed for use during spacecraft on-orbit operations characterized by significant mass properties variations. Feasibility of the gain adaptation is demonstrated via a Space Station Freedom/Orbiter berthing operation scenario. The attitude controller transient performance during gain adaptation appear satisfactory.

Parlos, Alexander G.↗

Modelling issues of the SCOLE control problem

Nonlinear and linear equations of motions were derived. The preliminary investigation consisted of model beam as truss structure, effects of truss structure on control design, and effects of reflector offset on control design. It was concluded that the offset of the reflector c.g. from the beam reflector attach point is dynamically significant. Also, truss effects may also significantly effect the performance of the controller if ignored. If the truss is included in the modeling of the NASA/SCOLE configuration, a practically implementable scheme is available to reduce the model order.

Hotz, A. F.↗

A Combined XRD/XRF Instrument for Lunar Resource Assessment

Robotic surface missions to the Moon should be capable of measuring mineral as well as chemical abundances in regolith samples. Although much is already known about the lunar regolith, our data are far from comprehensive. Most of the regolith samples returned to Earth for analysis had lost the upper surface, or it was intermixed with deeper regolith. This upper surface is the part of the regolith most recently exposed to the solar wind; as such it will be important to resource assessment. In addition, it may be far easier to mine and process the uppermost few centimeters of regolith over a broad area than to engage in deep excavation of a smaller area. The most direct means of analyzing the regolith surface will be by studies in situ. In addition, the analysis of the impact-origin regolith surfaces, the Fe-rich glasses of mare pyroclastic deposits, are of resource interest, but are inadequately known; none of the extensive surface-exposed pyroclastic deposits of the Moon have been systematically sampled, although we know something about such deposits from the Apollo 17 site. Because of the potential importance of pyroclastic deposits, methods to quantify glass as well as mineral abundances will be important to resource evaluation. Combined x ray diffraction (XRD) and x ray fluorescence (XRF) analysis will address many resource characterization problems on the Moon. XRF methods are valuable for obtaining full major-element abundances with high precision. Such data, collected in parallel with quantitative mineralogy, permit unambiguous determination of both mineral and chemical abundances where concentrations are high enough to be of resource grade. Collection of both XRD and XRF data from a single sample provides simultaneous chemical and mineralogic information. These data can be used to correlate quantitative chemistry and mineralogy as a set of simultaneous linear equations, the solution of which can lead to full characterization of the sample. The use of Rietveld methods for XRD data analysis can provide a powerful tool for quantitative mineralogy and for obtaining crystallographic data on complex minerals.

Vaniman, D. T.↗

Analysis, modelling and simulation of the large-angle magnetic suspension test fixture

The large-angle magnetic suspension test fixture (LAMSTF) research project is described which is based on a sophisticated finite element computer program, VF/GFUN, for calculating magnetic fields. The LAMSTF design includes eddy current paths the effect of which on system dynamics is being studied to modify the system dynamic model and the digital controller. Linearized equations of motions have been developed for system modelling and analysis of controllers.

Britcher, Colin P.↗

Aerodynamic shape optimization via sensitivity analysis on decomposed computational domains

Direct and iterative method considered to be most applicable to large systems of linear equations arising in discrete sensitivity analysis are assessed. Based on a single-domain grid, computations are performed using a banded matrix solver and an iterative solver, the generalized minimum residual (GMRES) method. The banded matrix solver is found to be generally the most economical method for those applications where the number of right-hand sides is large (i.e., a large number of design variables or a large number of adjoint vectors). For systems of equations that are too large to be solved by direct methods, an approach is proposed whereby the computational domain is divided into small subdomains, and each subdomain is solved separately.

Eleshaky, Mohamed E.↗

The evolving intergalactic medium - The uncollapsed baryon fraction in a cold dark matter universe

The time-varying density of the intergalactic medium (IGM) is calculated by coupling detailed numerical calculations of the thermal and ionization balance and radiative transfer in a uniform IGM of H and He to the linearized equations for the growth of density fluctuations in both gases and a dark component in a cold dark matter universe. The IGM density is identified with the collapsed baryon fraction. It is found that even if the IGM is never reheated, a significant fraction of the baryons remain uncollapsed at redshifts of four. If instead the collapsed fraction releases enough ionizing radiation or thermal energy to reionize the IGM by z greater than four as required by the Gunn-Peterson (GP) constraint, the uncollapsed fraction at z of four is even higher. The known quasar distribution is insufficient to supply the ionizing radiation necessary to satisfy the GP constraint in this case and, if stars are instead responsible, a substantial metallicity must have been produced by z of four.

Shapiro, Paul R.↗

Numerical quasi-linear study of the critical ionization velocity phenomenon

The critical ionization velocity (CIV) for a neutral barium (Ba) gas cloud moving across the static magnetic field is studied numerically using quasi-linear equations and a parameter range which is typical for the shaped-charge Ba gas release experiments in space. For consistency the charge exchange between the background oxygen ions and neutral atoms and its reverse process, as well as the excitation of the neutral Ba atoms, are included. The numerical results indicate that when the ionization rate due to CIV becomes comparable to the charge exchange rate the energy lost to the ionization and excitation collisions by the superthermal electrons exceeds the energy gain from the waves that are excited by the ion beam. This results in a CIV yield less than the yield by the charge exchange process.

Moghaddam-Taaheri, E.↗

A new parallel-vector finite element analysis software on distributed-memory computers

A new parallel-vector finite element analysis software package MPFEA (Massively Parallel-vector Finite Element Analysis) is developed for large-scale structural analysis on massively parallel computers with distributed-memory. MPFEA is designed for parallel generation and assembly of the global finite element stiffness matrices as well as parallel solution of the simultaneous linear equations, since these are often the major time-consuming parts of a finite element analysis. Block-skyline storage scheme along with vector-unrolling techniques are used to enhance the vector performance. Communications among processors are carried out concurrently with arithmetic operations to reduce the total execution time. Numerical results on the Intel iPSC/860 computers (such as the Intel Gamma with 128 processors and the Intel Touchstone Delta with 512 processors) are presented, including an aircraft structure and some very large truss structures, to demonstrate the efficiency and accuracy of MPFEA.

Qin, Jiangning↗

Geometrically nonlinear design sensitivity analysis on parallel-vector high-performance computers

Parallel-vector solution strategies for generation and assembly of element matrices, solution of the resulted system of linear equations, calculations of the unbalanced loads, displacements, stresses, and design sensitivity analysis (DSA) are all incorporated into the Newton Raphson (NR) procedure for nonlinear finite element analysis and DSA. Numerical results are included to show the performance of the proposed method for structural analysis and DSA in a parallel-vector computer environment.

Baddourah, Majdi A.↗

Finite element procedures for coupled linear analysis of heat transfer, fluid and solid mechanics

Coupled finite element formulations for fluid mechanics, heat transfer, and solid mechanics are derived from the conservation laws for energy, mass, and momentum. To model the physics of interactions among the participating disciplines, the linearized equations are coupled by combining domain and boundary coupling procedures. Iterative numerical solution strategy is presented to solve the equations, with the partitioning of temporal discretization implemented.

Sutjahjo, Edhi↗

Application of GPS attitude determination to gravity gradient stabilized spacecraft

Recent advances in the Global Positioning System (GPS) technology have initiated a new era in aerospace navigation and control. GPS receivers have become increasingly compact and affordable, and new developments have made attitude determination using subcentimeter positioning among two or more antennas feasible for real-time applications. GPS-based attitude control systems will become highly portable packages which provide time, navigation, and attitude information of sufficient accuracy for many aerospace needs. A typical spacecraft application of GPS attitude determination is a gravity gradient stabilized satellite in low Earth orbit that employs a GPS receiver and four body mounted patch antennas. The coupled, linearized equations of motion enable complete position and attitude information to be extracted from only two antennas. A discussion of the various error sources for spaceborne GPS attitude measurement systems is included. Attitude determination of better than 0.3 degrees is possible for 1 meter antenna separation. Suggestions are provided to improve the accuracy of the attitude solution.

Lightsey, E. G.↗

A hydrodynamic model of the oscillating screen viscometer

The viscometer consists of an oscillating screen immersed in a fluid and free to rotate about an axis in its plane. The viscosity can be determined from the measured ratio of the periodic driving force to the screen motion when an adequate hydrodynamical model of the immersed oscillator is available. The screen is formed by a square mesh of thin wire whose dimensions invite comparison with asymptotic results for narrow hollow bodies translating in Stokes flow. These indicate that the closed hole structure of the grid plays an important role in determining its motion. It is shown that this role diminishes as the frequency increases. The computed results, obtained from a system of linear equations, are consistent with experimental values over the appropriate range of frequency.

Davis, A. M. J.↗

On subsonic compressible flows by a method of correspondence II : application of methods to studies of flow with circulation about a circular cylinder

A general method for studying the flow of a compressible fluid around a closed body was discussed in Part I of this report. Here, application is made to the specific case in which the linearized equation of state is used. For a given incompressible flow around a specific profile, a corresponding compressible flow is found. The flow at infinity remains unchanged. Detailed studies are made of the flow with circulation around a unit circle, and velocity distributions are found for a wide range of Mach number and angle of attack. Comparisons are made with other methods.

GELBART ABE↗

Theoretical lift and drag of thin triangular wings at supersonic speeds

A method is derived for calculating the lift and the drag due to lift of point-forward triangular wings and a restricted series of sweptback wings at supersonic speeds. The elementary or "supersonic sources" solution of the linearized equation of motion is used to find the potential function of a line of doublets. The flow about the triangular flat plate is then obtained by a surface distribution of these doublet lines. The lift-curve slope of triangular wings is found to be a function of the ratio of the tangent of the apex angle to the tangent of the Mach angle. As the apex angle approaches and becomes greater than the Mach angle, the lift coefficient of the triangular wing becomes equal to that of a two-dimensional supersonic airfoil at the same Mach number. The drag coefficient due to lift of triangular wings with leading edges well behind the Mach cone is shown to be close to that of elliptically loaded wings of the same aspect ratio in subsonic flight. The resultant force on wings with leading edges outside the Mach cone, however, is shown to act normal to the surfaces and thus an induced drag equal to the lift times the angle of attack is obtained.

Brown, Clinton E↗

Theoretical study of air forces on an oscillating or steady thin wing in a supersonic main stream

A theoretical study, based on the linearized equations of motion for small disturbance, is made of the air forces on wings of general plan forms moving forward at a constant supersonic speed. The boundary problem is set up for both the harmonically oscillating and the steady conditions. Two types of boundary conditions are distinguished, which are designated "purely supersonic" and "mixed supersonic." the method is illustrated by applications to a number of examples for both the steady and the oscillating conditions. The purely supersonic case involves independence of action of the upper and lower surfaces of the airfoil and present analysis is mainly concerned with this case. A discussion is first given of the fundamental or elementary solution corresponding to a moving source. The solutions for the velocity potential are then synthesized by means of integration of the fundamental solution for the moving source. The method is illustrated by applications to a number of examples for both the steady and the oscillating cases and for various plan forms, including swept wings and rectangular and triangular plan forms. The special results of a number of authors are shown to be included in the analysis.

Garrick, I E↗