SEARCH · Engineering Papers

Results for “Reduce order modeling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35

All sources NASA NTRS DOE OSTI

NASA NTRS2017

Transient Thermal Analyses of Passive Systems on SCEPTOR X-57

As efficiency, emissions, and noise become increasingly prominent considerations in aircraft design, turning to an electric propulsion system is a desirable solution. Achieving the intended benefits of distributed electric propulsion (DEP) requires thermally demanding high power systems, presenting a different set of challenges compared to traditional aircraft propulsion. The embedded nature of these heat sources often preclude the use of traditional thermal management systems in order to maximize performance, with less opportunity to exhaust waste heat to the surrounding environment. This paper summarizes the thermal analyses of X-57 vehicle subsystems that don't employ externally air-cooled heat sinks. The high-power battery, wires, high-lift motors, and aircraft outer surface are subjected to heat loads with stringent thermal constraints. The temperature of these components are tracked transiently, since they never reach a steady-state equilibrium. Through analysis and testing, this report demonstrates that properly characterizing the material properties is key to accurately modeling peak temperature of these systems, with less concern for spatial thermal gradients. Experimentally validated results show the thermal profile of these systems can be sufficiently estimated using reduced order approximations.

Thermal Analyses↗

NASA NTRS1976

Low-perigee aerodynamic heating during orbital flight of an atmosphere Explorer

An extensive, low-perigee orbital aerodynamic heating study was undertaken in support of the Atmosphere Explorer-C Temperature Alarm. State of the art of low-density, high-speed flows, some models of the earth's atmosphere, external flow-field definition, thermodynamic and transport properties of atmospheric gases, the accommodation coefficient orbital thermal environment, and correlation of theory and measurements are discussed. Aerodynamic heating rates are determined for eight selected orbits by means of a reduced, analytical model verified by both ground test and flight data. These heating rates are compared with classical free-molecule and first-order collision regime values.

Caruso, P. S., Jr.↗

NASA NTRS1976

Theoretical and empirical low perigee aerodynamic heating during orbital flight of an atmosphere explorer

This document presents the results of an extensive, low perigee, orbital aerodynamic heating study undertaken in support of the Atmosphere Explorer-C Temperature Alarm. Based upon in-flight orbital temperature data from the Temperature Alarm tungsten resistance wire thermometer, aerodynamic heating rates have been determined for eight selected orbits by means of a reduced thermal analytical model verified by both ground test and flight data. These heating rates are compared with the classical free molecular and first order collision regime values. It has been concluded that, for engineering purposes, the aerodynamic heating rate of atmospheric gases at perigee altitudes between 170 and 135 km on pure tungsten wire is 30 to 60% of the value set by the classical free molecular limit. Relative to the more usual orbital thermal input attributable to direct solar radiation, the aerodynamic heating rate at the lowest altitude attempted with the spacecraft despun (135 km) is the equivalent of about 1.2 solar constants incident on a tungsten wire with a solar absorptivity of 0.85.

Caruso, P. S., Jr.↗

NASA NTRS1978

Aircraft model prototypes which have specified handling-quality time histories

Several techniques for obtaining linear constant-coefficient airplane models from specified handling-quality time histories are discussed. The pseudodata method solves the basic problem, yields specified eigenvalues, and accommodates state-variable transfer-function zero suppression. The algebraic equations to be solved are bilinear, at worst. The disadvantages are reduced generality and no assurance that the resulting model will be airplane like in detail. The method is fully illustrated for a fourth-order stability-axis small motion model with three lateral handling quality time histories specified. The FORTRAN program which obtains and verifies the model is included and fully documented.

Johnson, S. H.↗

NASA NTRS1982

Evaluation and reduction of errors induced by the Guyan transformation

The Guyan reduction refers to a method used to reduce the number of degrees of freedom in a structural model for dynamic analysis. Experience shows that, if the method is properly employed, then this reduction method does in fact provide a "reasonably" accurate approximation of the dynamic characteristics of the unreduced model. A method that provides an accurate estimate of the error in each mode of the reduced eigenvalue problem and, in addition, gives first order correction terms that greatly improve the accuracy of the eigenvectors and eigenvalues is presented. The method is demonstrated by standard COSMIC NASTRAN DAMP alters to rigid format 3, real normal mode analysis. Comparative computer time of modal extraction versus error analysis is discussed for the VAX-11/780 version of COSMIC NASTRAN.

Fox, G. L.↗

NASA NTRS1992

Adding dynamic rules to self-organizing fuzzy systems

This paper develops a Dynamic Self-Organizing Fuzzy System (DSOFS) capable of adding, removing, and/or adapting the fuzzy rules and the fuzzy reference sets. The DSOFS background consists of a self-organizing neural structure with neuron relocation features which will develop a map of the input-output behavior. The relocation algorithm extends the topological ordering concept. Fuzzy rules (neurons) are dynamically added or released while the neural structure learns the pattern. The DSOFS advantages are the automatic synthesis and the possibility of parallel implementation. A high adaptation speed and a reduced number of neurons is needed in order to keep errors under some limits. The computer simulation results are presented in a nonlinear systems modelling application.

Buhusi, Catalin V.↗

NASA NTRS1994

Linear quadratic Gaussian and feedforward controllers for the DSS-13 antenna

The controller development and the tracking performance evaluation for the DSS-13 antenna are presented. A trajectory preprocessor, linear quadratic Gaussian (LQG) controller, feedforward controller, and their combination were designed, built, analyzed, and tested. The antenna exhibits nonlinear behavior when the input to the antenna and/or the derivative of this input exceeds the imposed limits; for slewing and acquisition commands, these limits are typically violated. A trajectory preprocessor was designed to ensure that the antenna behaves linearly, just to prevent nonlinear limit cycling. The estimator model for the LQG controller was identified from the data obtained from the field test. Based on an LQG balanced representation, a reduced-order LQG controller was obtained. The feedforward controller and the combination of the LQG and feedforward controller were also investigated. The performance of the controllers was evaluated with the tracking errors (due to following a trajectory) and the disturbance errors (due to the disturbances acting on the antenna). The LQG controller has good disturbance rejection properties and satisfactory tracking errors. The feedforward controller has small tracking errors but poor disturbance rejection properties. The combined LQG and feedforward controller exhibits small tracking errors as well as good disturbance rejection properties. However, the cost for this performance is the complexity of the controller.

Gawronski, W. K.↗

NASA NTRS1995

Structural Analysis of Composite Test panels

An analytical study was conducted for a 1/4 scale model skin-stringer compression panel using STructural Analysis of General Shells (STAGS) code for finite element analysis. This was a study to investigate scaling methods for composite structures using reduced ply thicknesses for 1/2 and 1/4 scale models. A computer model for a curved panel in a pressure-box was also constructed in order to simulate test conditions in the pressure box test machine. The pressure box is a test facility at NASA Langley Research Center used to simulate flight conditions in fuselage panels. This model will be used to investigate modifications to the pressure-box for curved panels with different radii.

Berhan, Sean↗

NASA NTRS2014

Analytical-HZETRN Model for Rapid Assessment of Active Magnetic Radiation Shielding

The use of active radiation shielding designs has the potential to reduce the radiation exposure received by astronauts on deep-space missions at a significantly lower mass penalty than designs utilizing only passive shielding. Unfortunately, the determination of the radiation exposure inside these shielded environments often involves lengthy and computationally intensive Monte Carlo analysis. In order to evaluate the large trade space of design parameters associated with a magnetic radiation shield design, an analytical model was developed for the determination of flux inside a solenoid magnetic field due to the Galactic Cosmic Radiation (GCR) radiation environment. This analytical model was then coupled with NASA's radiation transport code, HZETRN, to account for the effects of passive/structural shielding mass. The resulting model can rapidly obtain results for a given configuration and can therefore be used to analyze an entire trade space of potential variables in less time than is required for even a single Monte Carlo run. Analyzing this trade space for a solenoid magnetic shield design indicates that active shield bending powers greater than 15 Tm and passive/structural shielding thicknesses greater than 40 g/cm2 have a limited impact on reducing dose equivalent values. Also, it is shown that higher magnetic field strengths are more effective than thicker magnetic fields at reducing dose equivalent.

Washburn, S. A.↗

NASA NTRS2003

Task II: Three-dimensional Rotating Stall Inception and Effects of Rotating Tip Clearance Asymmetry in Axial Compressors

The effects of two types of flow non-uniformity on stall inception behavior were assessed with linearized stability analyses of two compressor flow models. Response to rotating tip clearance asymmetries induced by a whirling rotor shaft or rotor height variations were investigated with a two-dimensional flow model. A 3-D compressor model was also developed to study the stability of both full-span and part-span rotating stall modes in annular geometries with radial flow variations. The studies focussed on (1) understanding what compressor designs were sensitive to these types of circumferential and spanwise flow non-uniformities, and (2) situations where 2-D stability theories were inadequate because of 3-D flow effects. Rotating tip clearance non-uniformity caused the greatest performance loss for shafts whirling at the rotating stall frequency. A whirling shaft displacement of 1 percent chord caused the stalling mass flow to rise by as much as 10 percent and the peak pressure rise to decrease by 6 percent. These changes were an order of magnitude larger than for equivalent-sized stationary or rotor-locked clearance asymmetries. Spanwise flow non-uniformity always destabilized the compressor, so that 2-D models over-predicted that stall margin compared to 3-D theory. The difference increased for compressors with larger spanwise variations of characteristic slope and reduced characteristic curvature near the peak. Differences between 2-D and 3-D stall point predictions were generally unacceptable (2 - 4 percent of flow coefficient) for single-stage configurations, but were less than 1 percent for multistage compressors. 2-D analyses predicted the wrong stall mode for specific cases of radial inlet flow distortion, mismatching and annulus area contraction, where higher-order radial modes led to stall. The stability behavior of flows with circumferential or radial non-uniformity was unified through a single stability criterion. The stall point for both cases was set by the integral around the annulus of the pressure rise characteristic slope, weighted by the amplitude of the mode shape. For the case of steady circumferential variations, this criterion reduced to the integrated mean slope (IMS) condition associated with steady inlet distortions. The rotating tip clearance asymmetry model was also used to demonstrate the feasibility of actively controlling the shaft position to suppress rotating stall. In axisymmetric mean flow, this method only stabilized the first harmonic mode, increasing the operating range until surge or higher harmonic modes became unstable.

Suder, Kenneth↗

NASA NTRS2014

Comparison Between THOR Anthropomorphic Test Device and THOR Finite Element Model

Extended time spent in reduced gravity can cause physiologic deconditioning of astronauts, reducing their ability to sustain excessive forces during dynamic phases of spaceflight such as landing. To make certain that the crew is safe during these phases, NASA must take caution when determining what types of landings are acceptable based on the accelerations applied to the astronaut. In order to test acceptable landings, various trials have been run accelerating humans, cadavers, and Anthropomorphic Test Devices (ATDs), or crash test dummies, at different acceleration and velocity rates on a sled testing platform. Using these tests, risks of injury will be created and metrics will be developed for the likelihood of injuries due to the acceleration. A finite element model (FEM) of the Test Device for Human Occupant Restraint (THOR) ATD has been developed that can simulate these test trials and others (Putnam, 2014), reducing the need for human and ATD testing. Additionally, this will give researchers a more effective way to test the accelerations and orientations encountered during spaceflight landings during design of new space vehicles for crewed missions. However, the FEM has not been proven and must be validated by comparing the forces, accelerations, and other measurements of all parts of the body between the physical tests already completed and computer simulated trials. The purpose of my research was to validate the FEM for the ATD using previously run trials with the physical THOR ATD.

Moore, Erik↗

NASA NTRS1974

Decomposition-aggregation stability analysis of the spinning Skylab

Stability of an 11-th order linear model of the spinning Skylab is determined by the decomposition-aggregation method based upon the comparison principle and vector Liapunov functions. To reduce the inherent conservativeness of the method an optimization problem is formulated and resolved producing the optimum comparison system. The system provides the best estimate of the stability region of the important structural parameter - asymmetry in the boom settings.

Cuk, S. M.↗

NASA NTRS1977

The effect of blade sweep on propeller performance

The advanced propeller developed for high Mach number cruise incorporates swept blades to reduce compressibility losses. In order to evaluate the induced flow-field vortex lattice methods are applied to a swept propeller blade. The blade is modeled by a radial distribution of helical horseshoe vortices with a single swept bound vortex at the quarter chord and the control point at the three-quarter chord of each radial section. The results of numerical calculations show that the power coefficient decreases as the blade is swept and the power loading distribution shifts inboard.

Sullivan, J. P.↗

NASA NTRS1993

Controllability and observability measures for Craig-Bampton substructure representations

Measures of controllability and observability are derived for a Craig-Bampton substructure representation for which fixed interface modes and constraint modes are computed relative to control actuator locations. The measures facilitate modal ordering and act as an index of dynamic completeness for modal truncation. Example problems are given to demonstrate that the controllable and observable reduced substructure models have a greater accuracy in predicting closed loop pole locations, transfer functions and actuation commands than conventional normal mode models.

Triller, Michael J.↗

NASA NTRS1996

Low-Dimensional Dynamical Models of Thermal Convection

A low-dimensional dynamic model for transitional buoyancy-driven flow in a differentially heated tall enclosure is presented. The full governing partial differential equations with the associated boundary conditions are solved by a spectral element method for a cavity of aspect ratio A=20. Proper orthogonal decomposition is applied to the oscillatory solution at Prandtl number Pr=P tau (omega) = 0.71 and Grashof number G tau (omega) = 3.2 x 10 (exp 4) to construct empirical eigenfunctions. Using the four most energetic empirical eigenfunctions for the velocity and temperature as basis functions and applying Galerkin's method, a reduced model consisting of eight nonlinear ordinary differential equations is obtained. Close to the 'design' conditions (P tau(omega) G tau(omega)), the low-order model (LOM) predictions are in excellent agreement with the predictions of the full model. In particular, the critical Grashof number at the onset of the first temporal flow instability (Hopf bifurcation) was well as the frequency and amplitude of oscillations at supercritical conditions are in excellent agreement with the predictions of the full model. Far from the 'design' conditions, the LOM predicts the existence of multiple stable steady solutions at large values of G tau, and a unique stable steady solution at small values of G tau, and exhibits hysteretic behavior that is qualitatively similar to that observed in direct numerical simulations based on the full model.

Liakopoulos, Anthony↗

NASA NTRS1997

Parameterizing Grid-Averaged Longwave Fluxes for Inhomogeneous Marine Boundary Layer Clouds

This paper examines the relative impacts on grid-averaged longwave flux transmittance (emittance) for Marine Boundary Layer (MBL) cloud fields arising from horizontal variability of optical depth tau and cloud sides, First, using fields of Landsat-inferred tau and a Monte Carlo photon transport algorithm, it is demonstrated that mean all-sky transmittances for 3D variable MBL clouds can be computed accurately by the conventional method of linearly weighting clear and cloudy transmittances by their respective sky fractions. Then, the approximations of decoupling cloud and radiative properties and assuming independent columns are shown to be adequate for computation of mean flux transmittance. Since real clouds have nonzero geometric thicknesses, cloud fractions A'(sub c) presented to isotropic beams usually exceed the more familiar vertically projected cloud fractions A(sub c). It is shown, however, that when A(sub c)less than or equal to 0.9, biases for all-sky transmittance stemming from use of A(sub c) as opposed to A'(sub c) are roughly 2-5 times smaller than, and opposite in sign to, biases due to neglect of horizontal variability of tau. By neglecting variable tau, all-sky transmittances are underestimated often by more than 0.1 for A(sub c) near 0.75 and this translates into relative errors that can exceed 40% (corresponding errors for all-sky emittance are about 20% for most values of A(sub c). Thus, priority should be given to development of General Circulation Model (GCM) parameterizations that account for the effects of horizontal variations in unresolved tau, effects of cloud sides are of secondary importance. On this note, an efficient stochastic model for computing grid-averaged cloudy-sky flux transmittances is furnished that assumes that distributions of tau, for regions comparable in size to GCM grid cells, can be described adequately by gamma distribution functions. While the plane-parallel, homogeneous model underestimates cloud transmittance by about an order of magnitude when 3D variable cloud transmittances are less than or equal to 0.2 and by approx. 20% to 100% otherwise, the stochastic model reduces these biases often by more than 80%.

Barker, Howard W.↗

NASA NTRS2001

Dynamic regulation of heart rate during acute hypotension: new insight into baroreflex function

To examine the dynamic properties of baroreflex function, we measured beat-to-beat changes in arterial blood pressure (ABP) and heart rate (HR) during acute hypotension induced by thigh cuff deflation in 10 healthy subjects under supine resting conditions and during progressive lower body negative pressure (LBNP). The quantitative, temporal relationship between ABP and HR was fitted by a second-order autoregressive (AR) model. The frequency response was evaluated by transfer function analysis. Results: HR changes during acute hypotension appear to be controlled by an ABP error signal between baseline and induced hypotension. The quantitative relationship between changes in ABP and HR is characterized by a second-order AR model with a pure time delay of 0.75 s containing low-pass filter properties. During LBNP, the change in HR/change in ABP during induced hypotension significantly decreased, as did the numerator coefficients of the AR model and transfer function gain. Conclusions: 1) Beat-to-beat HR responses to dynamic changes in ABP may be controlled by an error signal rather than directional changes in pressure, suggesting a "set point" mechanism in short-term ABP control. 2) The quantitative relationship between dynamic changes in ABP and HR can be described by a second-order AR model with a pure time delay. 3) The ability of the baroreflex to evoke a HR response to transient changes in pressure was reduced during LBNP, which was due primarily to a reduction of the static gain of the baroreflex.

Clinical Trial↗

NASA NTRS2013

GRAIL Refinements to Lunar Seismic Structure

Joint interpretation of disparate geophysical datasets helps reduce drawbacks that can result from analyzing them individually. The Apollo seismic network was situated on the lunar nearside surface in a roughly equilateral triangle having sides approximately 1000 km long, with stations 12/14 nearly co-located at one corner. Due to this limited geographical extent, near-surface ray coverage from moonquakes is low, but increase with depth. In comparison, gravity surveys and their resulting gravity anomaly maps have traditionally offered optimal resolution at crustal depths. Gravimetric maps and seismic data sets are therefor well suited to joint inversion, since the complementary information reduces inherent model ambiguity. We will perform a joint inversion of Apollo seismic delay times and gravity data collected by GRAIL lunar gravity mission, in order to recover seismic velocity and density as a function of latitude, longitude and depth within the Moon. We will relate density (rho) to seismic velocity (v) using a linear relationship that is allowed to be depth-dependent. The corresponding coefficient (B) can reflect a variety of material properties that vary with depth, including temperature and composition. The inversion seeks to recover the set of rho, v, and B perturbations that minimize (in a least-squares sense) the difference between the observed and calculated data.

Weber, Renee C.↗