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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
NASA NTRS2024

MSFC-TVC | TB-03: Derivation of Thrust Vector Control (TVC) Actuator-Force / Gimbal-Torque Transformation Matrix

Thrust vector control (TVC) systems for rocket engine propulsion traditionally use a simple linear relationship to convert between actuator forces and torques about the engine gimbal’s center-of-rotation (COR). As shown in Equation (1), the torque about the gimbal COR is proportional to the applied actuator force and the TVC moment arm (MA)—the perpendicular distance between the TVC actuator’s line-of-action (LOA) and the engine gimbal’s COR. While this fundamental relationship remains valid and accurate in a two-dimensional (2D), one-degree-of-freedom (1-DOF) context—particularly in its non-linear formulation as described in the ER63 Technical Bulletin TB-02 (Derivation of Thrust Vector Control (TVC) Engine-Gimbal / Actuator Moment-Arm Geometry)—it becomes limited when extended to three-dimensional (3D), two-degree-of-freedom (2-DOF) analyses. In 3D, out-of-plane angular displacements can arise, causing the gimbaled engine plane-of-motion to become non-coplanar with the respective actuator plane-of-motion. Such misalignments occur due to the rod-end (RE) and/or tail-stock (TE) mounting geometry of individual TVC actuators. These geometric complexities lead to inaccuracies in calculating engine torques and corresponding actuator forces if only the traditional TVC MA relationship is employed. For higher gimbal angular displacements (e.g., greater than 5 degrees), these inaccuracies become more pronounced, necessitating a more robust mathematical framework. This document presents a comprehensive 3D (2-DOF) derivation of a transformation matrix that accurately converts between gimbaled engine torques and TVC actuator forces. By incorporating the necessary geometric and rotational considerations, this new approach corrects the limitations of the traditional TVC MA method. Subsequent sections compare the newly formulated approach to the traditional method, demonstrating its enhanced accuracy and reliability for a broad range of gimbaled engine conditions.

Thrust Vector Control↗
NASA NTRS1994

A Galerkin method for linear PDE systems in circular geometries with structural acoustic applications

A Galerkin method for systems of PDE's in circular geometries is presented with motivating problems being drawn from structural, acoustic, and structural acoustic applications. Depending upon the application under consideration, piecewise splines or Legendre polynomials are used when approximating the system dynamics with modifications included to incorporate the analytic solution decay near the coordinate singularity. This provides an efficient method which retains its accuracy throughout the circular domain without degradation at singularity. Because the problems under consideration are linear or weakly nonlinear with constant or piecewise constant coefficients, transform methods for the problems are not investigated. While the specific method is developed for the two dimensional wave equations on a circular domain and the equation of transverse motion for a thin circular plate, examples demonstrating the extension of the techniques to a fully coupled structural acoustic system are used to illustrate the flexibility of the method when approximating the dynamics of more complex systems.

Smith, Ralph C.↗
NASA NTRS1988

The dynamics and control of large-flexible space structures, part 10

A mathematical model is developed to predict the dynamics of the proposed orbiting Spacecraft Control Laboratory Experiment (SCOLE) during the station keeping phase. The equations of motion are derived using a Newton-Euler formulation. The model includes the effects of gravity, flexibility, and orbital dynamics. The control is assumed to be provided to the system through the Shuttle's three torquers, and through six actuators located by pairs at two points on the mast and at the mass center of the reflector. The modal shape functions are derived using the fourth order beam equation. The generic mode equations are derived to account for the effects of the control forces on the modal shape and frequencies. The equations are linearized about a nominal equilibrium position. The linear regulator theory is used to derive control laws for both the linear model of the rigidized SCOLE as well as that of the actual SCOLE including the first four flexible modes. The control strategy previously derived for the linear model of the rigidized SCOLE is applied to the nonlinear model of the same configuration of the system and preliminary single axis slewing maneuvers conducted. The results obtained confirm the applicability of the intuitive and appealing two-stage control strategy which would slew the SCOLE system, as if rigid to its desired position and then concentrate on damping out the residual flexible motions.

Bainum, Peter M.↗
NASA NTRS1973

A guidance and navigation system for continuous low thrust vehicles

A midcourse guidance and navigation system for continuous low thrust vehicles is described. A set of orbit elements, known as the equinoctial elements, are selected as the state variables. The uncertainties are modelled statistically by random vector and stochastic processes. The motion of the vehicle and the measurements are described by nonlinear stochastic differential and difference equations respectively. A minimum time nominal trajectory is defined and the equation of motion and the measurement equation are linearized about this nominal trajectory. An exponential cost criterion is constructed and a linear feedback guidance law is derived to control the thrusting direction of the engine. Using this guidance law, the vehicle will fly in a trajectory neighboring the nominal trajectory. The extended Kalman filter is used for state estimation. Finally a short mission using this system is simulated. The results indicate that this system is very efficient for short missions.

Tse, C. J. C.↗
NASA NTRS1994

Use of rotation to suppress thermosolutal convection in directionally solidified binary alloys

Effects of rotation on onset of convection during plane-front directional solidification of Pb-Sn and the pseudobinary system mercury cadmium telluride (Hg(1-x)Cd(x)Te), and on dendritic solidification of Pb-Sn have been studied by means of linear stability analysis. Incorporating Coriolis and centrifugal accelerations into the momentum equation of Coriell et al., we find that under realistic processing conditions, a large degree of stabilization can be achieved using modest rotation rates for both Pb-Sn and mercury cadmium telluride. At a growth velocity of 5 micron/sec and nominal liquid-side temperature gradient of 200 K/cm in Pb-Sn, rotation at 500 rpm results in a hundredfold increase in the critical Sn concentration. Large increases in the maximum allowable growth velocity at fixed melt composition are also attainable with modest rotation rates. The effect is amplified under conditions of reduced gravitational acceleration. For Hg(1-x)Cd(x)Te, we have also studied the nonrotating case. The key differences are due to the existence of a composition range for Hg(1-x)Cd(x)Te in which the melt density has a local maximum as a function of temperature. When the melt solidifies by cooling from below, the liquid density may initially increase with distance above the interface, before ultimately decreasing as the melt temperature increases above the value at which the local density maximum occurs. In contrast to the Pb-Sn case where density depends monotonically on temperature and composition, for Hg(1-x)Cd(x)Te there exists a critical value of the growth velocity above which plane-front solidification is unstable for all bulk CdTe mole fractions. Again, rotation leads to significant inhibition of onset. We identify the predicted stabilization with the Taylor-Proudman mechanism by which rotation inhibits thermal convection in a single-component fluid heated from below. In a binary liquid undergoing solidification, rotation inhibits the onset of buoyancy-driven convection, and has no effect on the short-wavelength morphological instability. At large growth velocities, the plane-front interface between liquid and solid becomes unstable with respect to a morphological instability and solidification occurs dendritically, with a mushy zone of dendrites and interdendritic fluid separating the solid from the melt. For the Pb-Sn system, rotation substantially suppresses the onset of convection in the mushy zone and in the overlying liquid, holding open the promise that rotation can suppress freckling and other macrosegregation defects.

Pearlstein, Arne J.↗
NASA NTRS1976

High angle-of-attack stability-and-control analysis

Methods of linear systems analysis were applied to mathematical models of aircraft flying at high angle of attack and maneuver rate. First order longitudinal and lateral directional coupling is obtained by linearizing the complete nonlinear equations of motion about a generalized (quasi steady) trim point. Open loop stability boundaries are defined using the linear dynamic equations, and pilot in the loop effects are presented. Stability augmentation structures for maneuvering flight conditions are shown to be defined readily using optimal control theory.

Stengel, R. F.↗
NASA NTRS1987

Determination of rotor derivatives and rotor hub force and moment derivatives from flight measurements with the RSRA compound helicopter

A case study of the application of an integrated rotorcraft identification method to the linear modeling of rotor system dynamics and rotor hub loads is presented. Applying the method to flight data obtained from the RSRA compound helicopter, the rotor derivatives in the rotor state dynamic equation are identified along with blade equivalent damping and spring periodic coefficients. A rigid blade flapping equation of motion can be derived from the identified rotor state equation. It is shown that the concept of rotor hub load derivatives is useful for small maneuvering loads. The hub load derivatives can be used to relate the hub dynamic load to rotor/fuselage motion and applied inputs.

Wang, Ji C.↗
NASA NTRS1989

Control of flexible beams using a free-free active truss

An analytical and experimental study involving controlling flexible beams using a free-free active truss is presented. This work extends previous work in controlling flexible continua with active trusses which were configured with fixed-free boundary conditions. The following describes the Lagrangian approach used to derive the equations of motion for the active truss and the beams attached to it. A partial-state feedback control law is derived for this system based on a full-state feedback Linear Quadratic Regulator method. The analytical model is examined via numerical simulations and the results are compared to a similar experimental apparatus described herein. The results show that control of a flexible continua is possible with a free-free active truss.

Clark, W. W.↗
NASA NTRS1993

An upwind multigrid method for solving viscous flows on unstructured triangular meshes

A multigrid algorithm is combined with an upwind scheme for solving the two dimensional Reynolds averaged Navier-Stokes equations on triangular meshes resulting in an efficient, accurate code for solving complex flows around multiple bodies. The relaxation scheme uses a backward-Euler time difference and relaxes the resulting linear system using a red-black procedure. Roe's flux-splitting scheme is used to discretize convective and pressure terms, while a central difference is used for the diffusive terms. The multigrid scheme is demonstrated for several flows around single and multi-element airfoils, including inviscid, laminar, and turbulent flows. The results show an appreciable speed up of the scheme for inviscid and laminar flows, and dramatic increases in efficiency for turbulent cases, especially those on increasingly refined grids.

Bonhaus, Daryl Lawrence↗
NASA NTRS2003

Mapping Soil pH Buffering Capacity of Selected Fields

Soil pH buffering capacity, since it varies spatially within crop production fields, may be used to define sampling zones to assess lime requirement, or for modeling changes in soil pH when acid forming fertilizers or manures are added to a field. Our objective was to develop a procedure to map this soil property. One hundred thirty six soil samples (0 to 15 cm depth) from three Georgia Coastal Plain fields were titrated with calcium hydroxide to characterize differences in pH buffering capacity of the soils. Since the relationship between soil pH and added calcium hydroxide was approximately linear for all samples up to pH 6.5, the slope values of these linear relationships for all soils were regressed on the organic C and clay contents of the 136 soil samples using multiple linear regression. The equation that fit the data best was b (slope of pH vs. lime added) = 0.00029 - 0.00003 * % clay + 0.00135 * % O/C, r(exp 2) = 0.68. This equation was applied within geographic information system (GIS) software to create maps of soil pH buffering capacity for the three fields. When the mapped values of the pH buffering capacity were compared with measured values for a total of 18 locations in the three fields, there was good general agreement. A regression of directly measured pH buffering capacities on mapped pH buffering capacities at the field locations for these samples gave an r(exp 2) of 0.88 with a slope of 1.04 for a group of soils that varied approximately tenfold in their pH buffering capacities.

Weaver, A. R.↗
NASA NTRSDate not supplied

Transition Modeling Based on the Dual N-factor Method for the CRM-NLF Wind Tunnel Configuration

The dual N-factor method is used to model the boundary-layer transition over the common research model with natural laminar flow (CRM-NLF) aircraft configuration. The flow conditions match selected test conditions from a wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. The paper presents a systematic methodology for transition prediction in the presence of a dual shock system and extends the prior capability for iteratively coupled computational fluid dynamics (CFD) predictions to incorporate three-dimensional, transonic wings. The method employs stability computations based on the linear parabolized stability equations (PSE), along with a dual N-factor criterion. The iterative process begins with the fully turbulent Reynolds-averaged-Navier-Stokes (RANS) mean flow solution. For the first iteration, a mean flow solution is calculated with an imposed transition front that aligns with the shock front from the fully turbulent solution. Subsequently, stability computations are performed along a set of streamlines across the wing to calculate the amplification of planar Tollmien-Schlichting (TS) and stationary crossflow (CF) modes. The transition criterion based on the dual N-factor method is used to infer the updated transition front and the process is successively repeated until convergence of the solution. Within three iterations, the predicted fronts for angles of attack of 1.45, 1.98, 2.46 and 2.94 degrees and a mean-aerodynamic-chord Reynolds number equal to 15 million, approach visual convergence in most regions of the studied cases, and the resulting predictions are in good agreement with the transition fronts deduced from measurements of temperature-sensitive paint. Even though surface pressure measurements based on fully-turbulent flow agree well with the measured pressure coefficient distributions, strong viscous-inviscid interaction effects cause significant shifts in the shock locations based on the imposed transition front, underscoring the intrusive nature of static pressure measurements using surface mounted ports on the CRM-NLF configuration.

Boundary Layer Transition↗
NASA NTRS1983

Fractional representation theory - Robustness results with applications to finite dimensional control of a class of linear distributed systems

This paper reviews and extends the fractional representation theory. In particular, new and powerful robustness results are presented. This new theory is utilized to develop a preliminary design methodology for finite dimensional control of a class of linear evolution equations on a Banach space. The design is for stability in an input-output sense, but particular attention is paid to internal stability as well.

Nett, C. N.↗
NASA NTRS1972

Normal freezing of ideal ternary systems of the pseudobinary type

Perfect liquid mixing but no solid diffusion is assumed in normal freezing. In addition, the molar compositions of the freezing solid and remaining liquid, respectively, follow the solidus and liquidus curves of the constitutional diagram. For the linear case, in which both the liquidus and solidus are perfectly straight lines, the normal freezing equation giving the fraction solidified at each melt temperature and the solute concentration profile in the frozen solid was determined as early as 1902, and has since been repeatedly published. Corresponding equations for quadratic, cubic or higher-degree liquidus and solidus lines have also been obtained. The equation of normal freezing for ideal ternary liquid solutions solidified into ideal solid solutions of the pseudobinary type is given. Sample computations with the use of this new equation were made and are given for the Ga-Al-As system.

Li, C. H.↗
NASA NTRS1992

Computational simulation of coupled material degradation processes for probabilistic lifetime strength of aerospace materials

The research included ongoing development of methodology that provides probabilistic lifetime strength of aerospace materials via computational simulation. A probabilistic material strength degradation model, in the form of a randomized multifactor interaction equation, is postulated for strength degradation of structural components of aerospace propulsion systems subjected to a number of effects or primative variables. These primative variable may include high temperature, fatigue or creep. In most cases, strength is reduced as a result of the action of a variable. This multifactor interaction strength degradation equation has been randomized and is included in the computer program, PROMISS. Also included in the research is the development of methodology to calibrate the above described constitutive equation using actual experimental materials data together with linear regression of that data, thereby predicting values for the empirical material constraints for each effect or primative variable. This regression methodology is included in the computer program, PROMISC. Actual experimental materials data were obtained from the open literature for materials typically of interest to those studying aerospace propulsion system components. Material data for Inconel 718 was analyzed using the developed methodology.

Boyce, Lola↗
NASA NTRS1972

Analytic evaluation of display requirements for approach to landing

A pilot-vehicle-display model is used to study information and display requirements and the effects on system performance and reliability of pilot-induced randomness, wind gusts, configurational changes, etc. A brief description of a control theoretic systems model is given and its use and validity are demonstrated by applying it in a piloted approach to landing situation. The analysis procedure assumes that the vehicle dynamics are represented by linearized equations of motion.

Kleinman, D. L.↗
NASA NTRS1974

Limiting performance of ground transportation vehicles subject to transient loading

A computational approach to determining the limiting performance of vehicles subject to transient disturbances based on response variable criteria is set forth. For the purposes of a limiting performance study, the transportation system dynamics are described using second- or first-order equations in which the sought for quantity is a vector of time-varying functions called control or isolator forces that have replaced portions of the physical system. Computations are performed as a linear programming problem. As an example of the limiting performance of vehicles in protecting passengers or cargo under crash conditions, the problem of lading damage of a rail vehicle that is struck by another vehicle is considered.

Pilkey, W. D.↗
NASA NTRS1980

Linear stochastic optimal control and estimation problem

Problem involves design of controls for linear time-invariant system disturbed by white noise. Solution is Kalman filter coupled through set of optimal regulator gains to produce desired control signal. Key to solution is solving matrix Riccati differential equation. LSOCE effectively solves problem for wide range of practical applications. Program is written in FORTRAN IV for batch execution and has been implemented on IBM 360.

Geyser, L. C.↗