Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear equation systems”

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 1,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Simulation of the electron acoustic instability for a finite-size electron beam system

Simulations of the electron acoustic instability for a finite-size electron beam system are performed with a particle-in-cell code to investigate the heating phenomena associated with the instability and the width of the heating region. The results show that the beam radiates electrostatic electron acoustic waves whose decay time outside the beam agrees with the spatial decay length derived from the linear dispersion equation. The ambient cold electrons in a diffusion region surrounding the beam are heated to a higher temperature by absorbing the radiated electron acoustic waves, with the heating occurring mainly in the parallel direction. In the heat diffusion region, the temperature of the cold electrons decreases with distance from the beam with a temperature gradient length smaller than the decay length of the wave energy. The results are discussed with respect to the DE 1 plasma and wave observations in the polar cusp region.

Lin, C. S.↗
NASA NTRS1991

Closed-form solutions for linear regulator design of mechanical systems including optimal weighting matrix selection

Vibration in modern structural and mechanical systems can be reduced in amplitude by increasing stiffness, redistributing stiffness and mass, and/or adding damping if design techniques are available to do so. Linear Quadratic Regulator (LQR) theory in modern multivariable control design, attacks the general dissipative elastic system design problem in a global formulation. The optimal design, however, allows electronic connections and phase relations which are not physically practical or possible in passive structural-mechanical devices. The restriction of LQR solutions (to the Algebraic Riccati Equation) to design spaces which can be implemented as passive structural members and/or dampers is addressed. A general closed-form solution to the optimal free-decay control problem is presented which is tailored for structural-mechanical system. The solution includes, as subsets, special cases such as the Rayleigh Dissipation Function and total energy. Weighting matrix selection is a constrained choice among several parameters to obtain desired physical relationships. The closed-form solution is also applicable to active control design for systems where perfect, collocated actuator-sensor pairs exist.

Hanks, Brantley R.↗
NASA NTRS1992

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.↗
NASA NTRS1993

An algorithm for solving the perturbed gas dynamic equations

The present application of a compact, higher-order central-difference approximation to the linearized Euler equations illustrates the multimodal character of these equations by means of computations for acoustic, vortical, and entropy waves. Such dissipationless central-difference methods are shown to propagate waves exhibiting excellent phase and amplitude resolution on the basis of relatively large time-steps; they can be applied to wave problems governed by systems of first-order partial differential equations.

Davis, Sanford↗
NASA NTRS1993

Mechanically fastened composite laminates subjected to combined bearing-bypass and shear loading

Bolts and rivets provide a means of load transfer in the construction of aircraft. However, they give rise to stress concentrations and are often the source and location of static and fatigue failures. Furthermore, fastener holes are prone to cracks during take-off and landing. These cracks present the most common origin of structural failures in aircraft. Therefore, accurate determination of the contact stresses associated with such loaded holes in mechanically fastened joints is essential to reliable strength evaluation and failure prediction. As the laminate is subjected to loading, the contact region, whose extent is not known, develops between the fastener and the hole boundary through this contact region, which consists of slip and no-slip zones due to friction. The presence of the unknown contact stress distribution over the contact region between the pin and the composite laminate, material anisotropy, friction between the pin and the laminate, pin-hole clearance, combined bearing-bypass and shear loading, and finite geometry of the laminate result in a complex non-linear problem. In the case of bearing-bypass loading in compression, this non-linear problem is further complicated by the presence of dual contact regions. Previous research concerning the analysis of mechanical joints subjected to combined bearing-bypass and shear loading is non-existent. In the case of bearing-bypass loading only, except for the study conducted by Naik and Crews (1991), others employed the concept of superposition which is not valid for this non-linear problem. Naik and Crews applied a linear finite element analysis with conditions along the pin-hole contact region specified as displacement constraint equations. The major shortcoming of this method is that the variation of the contract region as a function of the applied load should be known a priori. Also, their analysis is limited to symmetric geometry and material systems, and frictionless boundary conditions. Since the contact stress distribution and the contact region are not known a priori, they did not directly impose the boundary conditions appropriate for modelling the contact and on-contact regions between the fastener and the hole. Furthermore, finite element analysis is not suitable for iterative design calculations for optimizing laminate construction in the presence of fasteners under complex loading conditions. In this study, the solution method developed by Madenci and Ileri (1992a,b) has been extended to determine the contact stresses in mechanical joints under combined bearing-bypass and shear loading, and bearing-bypass loading in compression resulting in dual contact regions.

Madenci, Erdogan↗
NASA NTRS2001

High Order Numerical Simulation of Sound Generated by the Kirchhoff Vortex

An improved high order finite difference method for low Mach number computational aeroacoustics (CAA) is described. The improvements involve the conditioning of the Euler equations in perturbation form to minimize numerical cancellation error, and the use of a stable non-dissipative sixth-order central spatial differencing for the interior points and third-order at the boundary points. The spatial difference operator satisfies the summation-by-parts property to guarantee strict stability for linear hyperbolic systems. Spurious high frequency oscillations are damped by a third-order characteristic-based filter. The objective of this paper is to apply these improvements in the simulation of sound generated by the Kirchhoff vortex.

Mueller, Bernhard↗
NASA NTRS2024

Core Body Temperature Predictions Using Metabolic Energy Expenditure and Heart Rate During Simulated Extravehicular Activity

Long duration spaceflight missions will require crew to become more autonomous in conducting extravehicular activities (EVA) without direct communication with Mission Control for biomedical support. To enable such autonomy, we are developing a Crew State and Risk Model (CSRM) as a collection of key physiology domains that drive EVA crew capabilities and workloads. One model component of CSRM is human thermal regulation. In this paper, customized development of a baseline model to predict core body temperature is presented using physiology inputs of heart rate and metabolic rate. The model development dataset included a baseline study where participants (n=6, equal male and female) performed a 5-hour EVA in a two-part session while wearing a hybrid space suit simulator (HS3). The first session included an end-to-end EVA traversing 1500 meters to a geology site and traversing back to a habitat conducting geology, payload relocation, and maintenance operations every 500 meters. The second session included standalone tasks of a 2000-meter traverse followed by geology tasks. Thermal measures of core body temperature, local skin temperature, liquid cooling garment temperature, heart rate, and metabolic rate were collected through the test duration. Multiple regression was used to build a linear equation to predict core body temperature from inputs of heart rate and metabolic rate. Heat storage was calculated via predicted core body temperature plus LCG and skin temperatures. The model was tested against a dataset from a pressurized suited test (n=6, equal male and female) in the NASA Active Response Gravity Offload System (ARGOS) conducting similar end-to-end EVA and standalone tasks. Predicted error of the model against the raw test cases was 0.2±0.15 °C. The baseline prediction of core body temperature using heart rates and metabolic rates allows for simple real-time tracking from data collected in-flight to monitor crew consumables and thermal flight limits during EVA.

Bradley Hoffmann↗
DOE OSTI2025

Examining astrophysical gas cloud collapse using an optical depth-scaled, x-ray-irradiated, carbon-foam sphere

When stellar radiation interacts with a molecular cloud, the cloud's fate depends on the strength of the incident radiation and the radiation's mean-free-path within the cloud [F. Bertoldi, Astrophys. J. 346, 735–755 (1989)]. Under the right conditions, the radiation compresses the cloud and a star formation may occur. Where and when the stellar formation occurs in the cloud's collapse are open questions. Direct observation of the complete star–cloud lifecycle is nearly impossible due to the immense timescales and distances over which the interaction occurs. Laboratory astrophysics offers a way to investigate such a system by scaling the important astrophysical parameters to the laboratory. This work describes laboratory experiments to study the radiation-driven implosion of clouds, using x rays from a laser-irradiated, thin, gold foil as a surrogate star and a carbon-foam sphere as a surrogate cloud. An optically thick system, theoretically corresponding to a star-forming regime, was selected by choice of the foam density. Gold foil and sphere motions were imaged by x-ray radiography. Radiographic images show the formation of an interface between rarefied gold and carbon plasmas, a shock moving into the sphere, and a blunting of the initial sphere's shape. Measurements show that the shock moved linearly around 64 μm/ns into the sphere, and the gold–carbon interface formed by 2 ns at the sphere edge remained stationary. The deformation of the sphere was driven by the incident radiation and not by mechanical pressures applied by gold plasma. The blunting of the sphere was likely due to the geometric reduction of flux near the sphere's poles. Higher x-ray flux near the sphere's equator caused high compression and a faster shock, which flattened the sphere. We will discuss the results and implications of our observations.

VanDervort, R. W. [University of Michigan 1 , Ann ↗
DOE OSTI2026

Large-scale white noise and cosmology

The generation of white noise on large scales is a generic property of the dynamics of physical systems described by local nonlinear partial differential equations. Nonlinearities prevent the small-scale dynamics from being erased by smoothing. Unresolved small-scale dynamics act as an uncorrelated (white or Poissonian) noise (seemingly stochastic but actually deterministic) contribution to large-scale dynamics. This white noise exists even when the dynamics is very nearly linear. In cases where the power spectrum is sub-Poissonian on large scales, this noise will dominate on the largest scale power no matter the amplitude of the inhomogeneities. Such is the case in the standard model of cosmology, where the primordial density power spectrum is expected to have an almost Harrison-Zel’dovich, P [ k ] ∼ k , spectrum on a much broader range of scales than can be observed. Even though linear gravitational evolution dominates nonlinear corrections by a factor of ∼ 10 5 , the nonobservation of white noise on the Hubble scale precludes the extrapolation of this power law below the comoving 1 pc scale. More generally, observation or nonobservation of large-scale white noise provides a powerful probe of the Universe on very small scales in the early Universe. Gravitational radiation, phase transitions, vorticity, and running of the spectral index are all phenomena that can be probed with large-scale white noise. Large-scale white noise is a nonoptional feature of all cosmological models but one which has not heretofore been appreciated.

Barenboim, Gabriela [Valencia U., IFIC; Valencia U↗
NASA NTRS1981

Parameter estimation of large flexible aerospace structures with application to the control of the Maypole Deployable Reflector

Systems such as the Maypole deployable reflector have a distributed parameter nature. The flexible column and hoop structure and the circular antenna of 30-100 meter diameter which it supports are described by partial, rather than ordinary, differential equations. Progress completed in reduced order modelling andd controller design and digital parameter estimation and control is summarized. Topics covered include depolyment and on-orbit operation; quasi-static (steady state) operation; dynamic distributed parameter system; autoregressive moving average identification; frequency domain procedures; direct or implicit active control; adaptive observers; parameter estimation using a linear reinforcement learning factor; feedback control; and reduced order modeling for nonlinear systems.

Balas, M. J.↗
NASA NTRS1979

Study of fatigue durability of advanced composite materials under conditions of accelerated loading

The effect of temperature on the tension-tension fatigue life of the T300/5208 graphite/epoxy angle-ply laminate system was investigated in an effort to develop an acceptable and reliable method of accelerated loading. Typical S log sub 10 N curves were determined experimentally at 25 C, 75 C, and 115 C. The time-temperature superposition principle was employed to find the shift factors of uniaxial fatigue strength, and a general linear equation of S log sub 10 N for shifting purpose was established. The combined techniques of cyclic creep measurements and optical microscopy upon fatigue failure allow some assessment of the possible physical basis of S log 10 N curve shifting. Before fatigue, the laminates at all test temperatures and stress levels undergo a unique damage mechanism during fatigue loading. It is concluded that an accelerated loading method is feasible.

Shih, H. M.↗
NASA NTRS1983

Simulating sunflower canopy temperatures to infer root-zone soil water potential

A soil-plant-atmosphere model for sunflower (Helianthus annuus L.), together with clear sky weather data for several days, is used to study the relationship between canopy temperature and root-zone soil water potential. Considering the empirical dependence of stomatal resistance on insolation, air temperature and leaf water potential, a continuity equation for water flux in the soil-plant-atmosphere system is solved for the leaf water potential. The transpirational flux is calculated using Monteith's combination equation, while the canopy temperature is calculated from the energy balance equation. The simulation shows that, at high soil water potentials, canopy temperature is determined primarily by air and dew point temperatures. These results agree with an empirically derived linear regression equation relating canopy-air temperature differential to air vapor pressure deficit. The model predictions of leaf water potential are also in agreement with observations, indicating that measurements of canopy temperature together with a knowledge of air and dew point temperatures can provide a reliable estimate of the root-zone soil water potential.

Choudhury, B. J.↗
NASA NTRS1986

Legendre-tau approximations for functional differential equations

The numerical approximation of solutions to linear retarded functional differential equations are considered using the so-called Legendre-tau method. The functional differential equation is first reformulated as a partial differential equation with a nonlocal boundary condition involving time-differentiation. The approximate solution is then represented as a truncated Legendre series with time-varying coefficients which satisfy a certain system of ordinary differential equations. The method is very easy to code and yields very accurate approximations. Convergence is established, various numerical examples are presented, and comparison between the latter and cubic spline approximation is made.

Ito, K.↗
NASA NTRS1988

Rapid slewing of the orbiting Spacecraft Control Laboratory Experiment (SCOLE) using LQR techniques

The rotational equations of motion, describing the dynamics of the (rigidized) proposed orbiting Spacecraft Control Laboratory Experiment (SCOLE) during the station keeping phase, are derived using the Eulerian formulation. When the attitude angles (roll, pitch, and yaw) are assumed small, a stability analysis is conducted for the system. It is seen that the pitch equation decouples from the roll and yaw equations when the interface between the mast on the reflector is not offset or the offset is only along the Shuttle roll axis. When a second offset is introduced along the pitch axis of the system and when the gravity-gradient torques are present in the dynamics, the system assumes a new equilibrium position. The linear regulator theory is used to derive a control law for the linear model of the rigidized SCOLE. This law is applied to the nonlinear model of the same configuration of the system and preliminary single axis slewing maneuvers (20 deg amplitude) are simulated.

Diarra, Cheick M.↗
NASA NTRS1992

Nonlinear stability and control study of highly maneuverable high performance aircraft, phase 2

This research should lead to the development of new nonlinear methodologies for the adaptive control and stability analysis of high angle-of-attack aircraft such as the F18 (HARV). The emphasis has been on nonlinear adaptive control, but associated model development, system identification, stability analysis and simulation is performed in some detail as well. Various models under investigation for different purposes are summarized in tabular form. Models and simulation for the longitudinal dynamics have been developed for all types except the nonlinear ordinary differential equation model. Briefly, studies completed indicate that nonlinear adaptive control can outperform linear adaptive control for rapid maneuvers with large changes in alpha. The transient responses are compared where the desired alpha varies from 5 degrees to 60 degrees to 30 degrees and back to 5 degrees in all about 16 sec. Here, the horizontal stabilator is the only control used with an assumed first-order linear actuator with a 1/30 sec time constant.

Mohler, R. R.↗
NASA NTRS1993

ISAC: A tool for aeroservoelastic modeling and analysis

The capabilities of the Interaction of Structures, Aerodynamics, and Controls (ISAC) system of program modules is discussed. The major modeling, analysis, and data management components of ISAC are identified. Equations of motion are displayed for a Laplace-domain representation of the unsteady aerodynamic forces. Options for approximating a frequency-domain representation of unsteady aerodynamic forces with rational functions of the Laplace variable are shown. Linear time invariant state-space equations of motion that result are discussed. Model generation and analyses of stability and dynamic response characteristics are shown for an aeroelastic vehicle which illustrates some of the capabilities of ISAC as a modeling and analysis tool for aeroelastic applications.

Adams, William M., Jr.↗
NASA NTRS1999

Eigenmodes of Ducted Flows With Radially-Dependent Axial and Swirl Velocity Components

This report characterizes the sets of small disturbances possible in cylindrical and annular ducts with mean flow whose axial and tangential components vary arbitrarily with radius. The linearized equations of motion are presented and discussed, and then exponential forms for the axial, circumferential, and time dependencies of any unsteady disturbances are assumed. The resultant equations form a generalized eigenvalue problem, the solution of which yields the axial wavenumbers and radial mode shapes of the unsteady disturbances. Two numerical discretizations are applied to the system of equations: (1) a spectral collocation technique based on Chebyshev polynomial expansions on the Gauss-Lobatto points, and (2) second and fourth order finite differences on uniform grids. The discretized equations are solved using a standard eigensystem package employing the QR algorithm. The eigenvalues fall into two primary categories: a discrete set (analogous to the acoustic modes found in uniform mean flows) and a continuous band (analogous to convected disturbances in uniform mean flows) where the phase velocities of the disturbances correspond to the local mean flow velocities. Sample mode shapes and eigensystem distributions are presented for both sheared axial and swirling flows. The physics of swirling flows is examined with reference to hydrodynamic stability and completeness of the eigensystem expansions. The effect of assuming exponential dependence in the axial direction is discussed.

Kousen, Kenneth A.↗
NASA NTRS1987

A study of some approximations in the modeling of a coupled rotor/fuselage aeromechanical system

Equations of motion for a three-bladed hingeless rotor coupled to a gimbal-supported fuselage are developed. A brief derivation of the equations is presented, including the assumptions made, the ordering scheme, and a discussion of both the source of the aerodynamic axial load term and the reason for differences due to hinge sequence. The solution of the equations is discussed, including linearization, the use of multiblade coordinates to eliminate periodic coefficients, and specialization to a root-sprung model. The treatment of the blade root offset by using an ordering scheme to neglect higher-order terms is found to affect results in pitch and roll. A small aerodynamic loading term in the blade axial direction has a significant effect at large collective pitch settings, as does the choice of a flap-lag or a lag-flap transformation sequence for blade motion. The results for a hingeless rotor also differ significantly from those for the root-sprung model when the nonrotating natural frequencies are matched.

Blelloch, P. A.↗