Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reduced order methods”

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 433 records · Page 24

Advanced Vibration Analysis Tool Developed for Robust Engine Rotor Designs

The primary objective of this research program is to develop vibration analysis tools, design tools, and design strategies to significantly improve the safety and robustness of turbine engine rotors. Bladed disks in turbine engines always feature small, random blade-to-blade differences, or mistuning. Mistuning can lead to a dramatic increase in blade forced-response amplitudes and stresses. Ultimately, this results in high-cycle fatigue, which is a major safety and cost concern. In this research program, the necessary steps will be taken to transform a state-of-the-art vibration analysis tool, the Turbo- Reduce forced-response prediction code, into an effective design tool by enhancing and extending the underlying modeling and analysis methods. Furthermore, novel techniques will be developed to assess the safety of a given design. In particular, a procedure will be established for using natural-frequency curve veerings to identify ranges of operating conditions (rotational speeds and engine orders) in which there is a great risk that the rotor blades will suffer high stresses. This work also will aid statistical studies of the forced response by reducing the necessary number of simulations. Finally, new strategies for improving the design of rotors will be pursued.

Min, James B.↗

The effect of a small initial curvature on the free vibration of clamped, rectangular plates

An analytical method of obtaining the natural frequencies and mode shapes of clamped, rectangular plates having a small initial curvature is presented. Specifically, the singular perturbation technique is used to reduce the fourth-order plate vibration problem to the simpler membrane problem with modified boundary conditions that account for the bending effects. The eigenfrequencies for plates with inverse aspect ratios varying between 0.1 and 1.0 and for the dimensionless normal prestress between 0.1 and 1.0 have been presented for values of epsilon, the normalized bending rigidity, ranging between 0.0010 and 0.2500. It is established that a small initial curvature has no effect on the frequency of vibration of the plate. However, its effect is manifested in the eigenmodes.

Adeniji-Fashola, A. A.↗

In Search of Grid Converged Solutions

Assessing solution error continues to be a formidable task when numerically solving practical flow problems. Currently, grid refinement is the primary method used for error assessment. The minimum grid spacing requirements to achieve design order accuracy for a structured-grid scheme are determined for several simple examples using truncation error evaluations on a sequence of meshes. For certain methods and classes of problems, obtaining design order may not be sufficient to guarantee low error. Furthermore, some schemes can require much finer meshes to obtain design order than would be needed to reduce the error to acceptable levels. Results are then presented from realistic problems that further demonstrate the challenges associated with using grid refinement studies to assess solution accuracy.

Lockard, David P.↗

Higher Order Modeling in Hybrid Approaches to the Computation of Electromagnetic Fields

Higher order geometry representations and interpolatory basis functions for computational electromagnetics are reviewed. Two types of vector-valued basis functions are described: curl-conforming bases, used primarily in finite element solutions, and divergence-conforming bases used primarily in integral equation formulations. Both sets satisfy Nedelec constraints, which optimally reduce the number of degrees of freedom required for a given order. Results are presented illustrating the improved accuracy and convergence properties of higher order representations for hybrid integral equation and finite element methods.

Wilton, Donald R.↗

Spin Testing for Durability Began on a Self-Tuning Impact Damper for Turbomachinery Blades

NASA and Pratt & Whitney will collaborate under a Space Act Agreement to perform spin testing of the impact damper to verify damping effectiveness and durability. Pratt & Whitney will provide the turbine blade and damper hardware for the tests. NASA will provide the facility and perform the tests. Effectiveness and durability will be investigated during and after sustained sweeps of rotor speed through resonance. Tests of a platform wedge damper are also planned to compare its effectiveness with that of the impact damper. Results from baseline tests without dampers will be used to measure damping effectiveness. The self-tuning impact damper combines two damping methods-the tuned mass damper and the impact damper. It consists of a ball located within a cavity in the blade. This ball rolls back and forth on a spherical trough under centrifugal load (tuned mass damper) and can strike the walls of the cavity (impact damper). The ball s rolling natural frequency is proportional to the rotor speed and can be designed to follow an engine-order line (integer multiple of rotor speed). Aerodynamic forcing frequencies typically follow these engineorder lines, and a damper tuned to the engine order will most effectively reduce blade vibrations when the resonant frequency equals the engine-order forcing frequency. This damper has been tested in flat plates and turbine blades in the Dynamic Spin Facility. During testing, a pair of plates or blades rotates in vacuum. Excitation is provided by one of three methods--eddy-current engine-order excitation (ECE), electromechanical shakers, and magnetic bearing excitation. The eddy-current system consists of magnets located circumferentially around the rotor. As a blade passes a magnet, a force is imparted on the blade. The number of magnets used can be varied to change the desired engine order of the excitation. The magnets are remotely raised or lowered to change the magnitude of the force on the blades. The other two methods apply force to the rotating shaft itself at frequencies independent of the rotor speed. During testing, blade vibration is monitored with strain gauges and laser displacement probes.

Duffy, Kirsten↗

A uniform geometrical optics and an extended uniform geometrical theory of diffraction for evaluating high frequency EM fields near smooth caustics and composite shadow boundaries

A uniform geometrical optics (UGO) and an extended uniform geometrical theory of diffraction (EUTD) are developed for evaluating high frequency electromagnetic (EM) fields within transition regions associated with a two and three dimensional smooth caustic of reflected rays and a composite shadow boundary formed by the caustic termination or the confluence of the caustic with the reflection shadow boundary (RSB). The UGO is a uniform version of the classic geometrical optics (GO). It retains the simple ray optical expressions of classic GO and employs a new set of uniform reflection coefficients. The UGO also includes a uniform version of the complex GO ray field that exists on the dark side of the smooth caustic. The EUTD is an extension of the classic uniform geometrical theory of diffraction (UTD) and accounts for the non-ray optical behavior of the UGO reflected field near caustics by using a two-variable transition function in the expressions for the edge diffraction coefficients. It also uniformly recovers the classic UTD behavior of the edge diffracted field outside the composite shadow boundary transition region. The approach employed for constructing the UGO/EUTD solution is based on a spatial domain physical optics (PO) radiation integral representation for the fields which is then reduced using uniform asymptotic procedures. The UGO/EUTD analysis is also employed to investigate the far-zone RCS problem of plane wave scattering from two and three dimensional polynomial defined surfaces, and uniform reflection, zero-curvature, and edge diffraction coefficients are derived. Numerical results for the scattering and diffraction from cubic and fourth order polynomial strips are also shown and the UGO/EUTD solution is validated by comparison to an independent moment method (MM) solution. The UGO/EUTD solution is also compared with the classic GO/UTD solution. The failure of the classic techniques near caustics and composite shadow boundaries is clearly demonstrated and it is shown that the UGO/EUTD results remain valid and uniformly reduce to the classic results away from the transition regions. Mathematical details on the asymptotic properties and efficient numerical evaluation of the canonical functions involved in the UGO/EUTD expressions are also provided.

Constantinides, E. D.↗

Modal cost analysis for linear matrix-second-order systems

Reduced models and reduced controllers for systems governed by matrix-second-order differential equations are obtained by retaining those modes which make the largest contributions to quadratic control objectives. Such contributions, expressed in terms of modal data, used as mode truncation criteria, allow the statement of the specific control objectives to influence the early model reduction from very high order models which are available, for example, from finite element methods. The relative importance of damping, frequency, and eigenvector in the mode truncation decisions are made explicit for each of these control objectives: attitude control, vibration suppression and figure control. The paper also shows that using modal cost analysis (MCA) on the closed loop modes of the optimally controlled system allows the construction of reduced control policies which feedback only those closed loop modal coordinates which are most critical to the quadratic control performance criterion. In this way, the modes which should be controlled (and hence the modes which must be observable by choice of measurements), are deduced from truncations of the optimal controller.

Skelton, R. E.↗

Application of global sensitivity equations in multidisciplinary aircraft synthesis

The present paper investigates the applicability of the Global Sensitivity Equation (GSE) method in the multidisciplinary synthesis of aeronautical vehicles. The GSE method provides an efficient approach for representing a large coupled system by smaller subsystems and accounts for the subsystem interactions by means of first-order behavior sensitivities. This approach was applied in an aircraft synthesis problem with performance constraints stemming from the disciplines of structures, aerodynamics, and flight mechanics. Approximation methods were considered in an attempt to reduce problem dimensionality and to improve the efficiency of the optimization process. The influence of efficient constraint representations, the choice of design variables, and design variable scaling on the conditioning of the system matrix was also investigated.

Hajela, P.↗

Linear and non-linear dynamic models of a geared rotor-bearing system

A three degree of freedom non-linear model of a geared rotor-bearing system with gear backlash and radial clearances in rolling element bearings is proposed here. This reduced order model can be used to describe the transverse-torsional motion of the system. It is justified by comparing the eigen solutions yielded by corresponding linear model with the finite element method results. Nature of nonlinearities in bearings is examined and two approximate nonlinear stiffness functions are proposed. These approximate bearing models are verified by comparing their frequency responses with the results given by the exact form of nonlinearity. The proposed nonlinear dynamic model of the geared rotor-bearing system can be used to investigate the dynamic behavior and chaos.

Kahraman, Ahmet↗

Computer Simulation Results and Analysis for a Root-Raised Cosine Filter Design using Canonical Signed Digits

Filters implemented using Canonical Signed-Digit (CSD) number representation for the coefficients have been studied for many years as an efficient way to increase the speed and reduce the hardware complexity. A simple and very effective method for designing a CSD filter with a specific set of filter parameters has been investigated and its simulation results are presented here. In order to optimize filter coefficients into the corresponding CSD numbers, the Minimum Mean Square Error (MMSE) criterion is used. Furthermore, an attempt is made to improve frequency response of the CSD filter by allocating an extra non-zero digit for normalized coefficients exceeding one-half. Due to limited filter aperture width in examples presented here, the frequency response of CSD filter as the number of non-zero digits in a CSD code increases is not affected much. A root-raised cosine filter model is employed as a base line for this design approach. Finally, the CSD filter simulation results and a hardware complexity comparison with a conventional filter are also shown along with the eye diagrams an Bit-Error-Rate (BER) performance curves.

Kim, Heechul↗

Modeling and Full-Scale Testing of an Aspirating Face Seal

A 36" diameter aspirating face seal for aircraft engine application has undergone extensive testing and analysis. Previous testing indicated that the seal tended to seek equilibrium at axial rotor clearances that were larger than expected. Parameter studies were conducted on several seal design parameters to evaluate effect on seal performance. Mixing of air flows from the air dam and air bearing regions of the seal was shown to have a significant impact on the seal's performance. Two methods of minimizing this flow interaction were studied both analytically and experimentally. The first method is to reduce the labyrinth tooth clearance, thereby limiting flow to the air dam itself. The second method involves utilizing a flow deflector between the air dam and air bearing regions of the seal in order to prevent radial flow from the air dam from disrupting the formation of a hydrostatic film at the air bearing. Both methods were shown to be effective design enhancements, allowing seal closure to be achieved. In both cases, the seal seeks an equilibrium position 0.0015" from the rotor surface, with corresponding low leakage rates.

Turnquist, Norman A.↗

Effects of a Rotating Aerodynamic Probe on the Flow Field of a Compressor Rotor

An investigation of distortions of the rotor exit flow field caused by an aerodynamic probe mounted in the rotor is described in this paper. A rotor total pressure Kiel probe, mounted on the rotor hub and extending up to the mid-span radius of a rotor blade channel, generates a wake that forms additional flow blockage. Three types of high-response aerodynamic probes were used to investigate the distorted flow field behind the rotor. These probes were: a split-fiber thermo-anemometric probe to measure velocity and flow direction, a total pressure probe, and a disk probe for in-flow static pressure measurement. The signals acquired from these high-response probes were reduced using an ensemble averaging method based on a once per rotor revolution signal. The rotor ensemble averages were combined to construct contour plots for each rotor channel of the rotor tested. In order to quantify the rotor probe effects, the contour plots for each individual rotor blade passage were averaged into a single value. The distribution of these average values along the rotor circumference is a measure of changes in the rotor exit flow field due to the presence of a probe in the rotor. These distributions were generated for axial flow velocity and for static pressure.

Lepicovsky, Jan↗

Subscale and Full-Scale Testing of Buckling-Critical Launch Vehicle Shell Structures

New analysis-based shell buckling design factors (aka knockdown factors), along with associated design and analysis technologies, are being developed by NASA for the design of launch vehicle structures. Preliminary design studies indicate that implementation of these new knockdown factors can enable significant reductions in mass and mass-growth in these vehicles and can help mitigate some of NASA s launch vehicle development and performance risks by reducing the reliance on testing, providing high-fidelity estimates of structural performance, reliability, robustness, and enable increased payload capability. However, in order to validate any new analysis-based design data or methods, a series of carefully designed and executed structural tests are required at both the subscale and full-scale level. This paper describes recent buckling test efforts at NASA on two different orthogrid-stiffened metallic cylindrical shell test articles. One of the test articles was an 8-ft-diameter orthogrid-stiffened cylinder and was subjected to an axial compression load. The second test article was a 27.5-ft-diameter Space Shuttle External Tank-derived cylinder and was subjected to combined internal pressure and axial compression.

Hilburger, Mark W.↗

Computing Sensitivities in Evolutionary Systems: A Real-time Reduced Order Modeling Strategy

We present a new methodology for computing sensitivities in evolutionary systems using a model-driven low-rank approximation. To this end, we formulate a variational principle that seeks to minimize the distance between the time derivative of the reduced approximation and sensitivity dynamics. The first order optimality condition of the variational principle leads to a system of closed form evolution equations for an orthonormal basis and corresponding sensitivity coefficients. This approach allows for the computation of sensitivities with respect to a large number of parameters in an accurate and tractable manner by extracting correlations between different sensitivities on the fly. The presented method requires solving forward evolution equations, sidestepping the restrictions imposed by the forward/backward workflow of adjoint sensitivities. For example, the presented method, unlike the adjoint equation, does not impose any input/output load and can be used in applications in which real-time sensitivities are of interest. We demonstrate the utility of the method for three test cases: (1) computing sensitivity with respect to model parameters in the Rössler system, (2) computing sensitivity with respect to an infinite-dimensional forcing parameter in the chaotic Kuramoto--Sivashinsky equation, and (3) computing sensitivity with respect to reaction parameters for species transport in a turbulent reacting flow.

Reduced order model↗

PSE-Based Aerodynamic Flow Transition Prediction Using Automated Unstructured CFD Integration

Accurate, robust, and efficient prediction of transition in viscous flows is a significant challenge in computational fluid dynamics. We present a coupled, high-fidelity, iterative approach that leverages the FUN3D flow solver and the LASTRAC stability code to predict transition in low-disturbance environments, initiated by the linear growth of boundary-layer instability modes. Our method integrates the ability of FUN3D to compute mixed laminar-transitional-turbulent mean flows via transition-sensitized Reynolds-Averaged Navier-Stokes equations with the ability of LASTRAC to perform linear stability analysis, all within an automated framework that requires no intermediate user involvement. Unlike conventional frameworks that rely on classical stability theory or reduced-order metamodels, our approach employs the parabolized stability equations to provide more accurate and reliable estimates of disturbance growth for multiple instability mechanisms, including Tollmien-Schlichting, Kelvin-Helmholtz, and crossflow modes. By accounting for the effects of mean-flow nonparallelism as well as the surface curvature, this approach lays the foundation for improved N-factor correlations for transition onset prediction in a broad class of flows. We apply this method to three distinct flow configurations: 1) flow over a zero-pressure-gradient flat plate, 2) the NLF-0416 airfoil with both natural and separation-induced transitions, and 3) a 6:1 prolate spheroid, where transition is primarily driven by crossflow instability. For the two-dimensional cases, a formulated intermittency distribution is used to model the transition zone in between the laminar and fully turbulent flows. The results include comparisons with experimental measurements, similar numerical approaches, and transport-equations-based models, demonstrating good agreement in surface pressure coefficients, transition onset locations, and skin-friction coefficients for all three configurations. Besides contributing a couple of new insights into boundary-layer transition in these canonical cases, this study provides a powerful tool for transition modeling in both research and design applications in aerodynamics.

Transition Prediction↗

Precision Asteroid Astrometry

Among the methods used to guide spacecraft to their destinations, Optical Navigation (OpNav) remains an effective option. OpNav makes use of star fields and small body ephemerides to precisely locate spacecraft. To facilitate accurate OpNav, the small body ephemerides must be constantly updated; asteroid orbits accumulate errors each year of a few milliarcseconds. Through extended exposures, taken with strategic offsets, a least-squares solution can be found that determines updated ephemeris data. This updated data can also be used by occultation astronomers to gain further information about the small bodies, including their size and shape. Using the 24-inch telescope at the Caltech Table Mountain Observatory (TMO), we capture two or more 180 second exposures of each target. These images, combined with a file for the predicted background star field and two reference files, are then processed through a series of scripts and programs. Starting with a prediction file and two to five exposures of the asteroid, the data is processed. This original data is about 32MB per observation. Once the data are reduced to only Right Ascension and Declination for each target, the data are ready for delivery. This consists of text only, and for each target takes about 80 bytes; this resulting data reduction is about five orders of magnitude. This method produces observed positions that are refined by about 12 milliarcseconds, a refinement that is accomplished almost nowhere else. The occultation observations that are facilitated by the ephemerides being refined also produce results that are not possible in any other way from ground-based observations.

Dial, Jason↗

Model reduction for discrete bilinear systems

A model reduction method for discrete bilinear systems is developed which matches q sets of Volterra and covariance parameters. These parameters are shown to represent both deterministic and stochastic attributes of the discrete bilinear system. A reduced order model which matches these q sets of parameters is defined to be a q-Volterra covariance equivalent realization (q-Volterra COVER). An algorithm is presented which constructs a class of q-Volterra COVERs parameterized by solutions to a Hermitian, quadratic, matrix equation. The algorithm is applied to a bilinear model of a robot manipulator.

King, A. M.↗

An Adaptive Semi-Implicit Scheme for Simulations of Unsteady Viscous Compressible Flows

A numerical scheme for simulation of unsteady, viscous, compressible flows is considered. The scheme employs an explicit discretization of the inviscid terms of the Navier-Stokes equations and an implicit discretization of the viscous terms. The discretization is second order accurate in both space and time. Under appropriate assumptions, the implicit system of equations can be decoupled into two linear systems of reduced rank. These are solved efficiently using a Gauss-Seidel method with multigrid convergence acceleration. When coupled with a solution-adaptive mesh refinement technique, the hybrid explicit-implicit scheme provides an effective methodology for accurate simulations of unsteady viscous flows. The methodology is demonstrated for both body-fitted structured grids and for rectangular (Cartesian) grids.

Steinthorsson, Erlendur↗