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.

SourceResearch connectionFrom yearThrough year

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

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

NASA NTRS2015

Projection-Based Reduced Order Modeling for Spacecraft Thermal Analysis

This paper presents a mathematically rigorous, subspace projection-based reduced order modeling (ROM) methodology and an integrated framework to automatically generate reduced order models for spacecraft thermal analysis. Two key steps in the reduced order modeling procedure are described: (1) the acquisition of a full-scale spacecraft model in the ordinary differential equation (ODE) and differential algebraic equation (DAE) form to resolve its dynamic thermal behavior; and (2) the ROM to markedly reduce the dimension of the full-scale model. Specifically, proper orthogonal decomposition (POD) in conjunction with discrete empirical interpolation method (DEIM) and trajectory piece-wise linear (TPWL) methods are developed to address the strong nonlinear thermal effects due to coupled conductive and radiative heat transfer in the spacecraft environment. Case studies using NASA-relevant satellite models are undertaken to verify the capability and to assess the computational performance of the ROM technique in terms of speed-up and error relative to the full-scale model. ROM exhibits excellent agreement in spatiotemporal thermal profiles (<0.5% relative error in pertinent time scales) along with salient computational acceleration (up to two orders of magnitude speed-up) over the full-scale analysis. These findings establish the feasibility of ROM to perform rational and computationally affordable thermal analysis, develop reliable thermal control strategies for spacecraft, and greatly reduce the development cycle times and costs.

0000↗

NASA NTRS1980

Optimal member damper controller design for large space structures

Consideration is given to the selection of velocity feedback gains for individual dampers for the members of a structurally controlled large flexible space structure. The problem is formulated as an optimal output feedback regulator problem, and necessary conditions are derived for minimizing a quadratic performance function. The diagonal nature of the gain matrix is taken into account, along with knowledge of noise covariances. It is pointed out that the method presented offers a systematic approach to the design of a class of controllers for enhancing structural damping, which have significant potential if used in conjunction with a reduced-order optimal controller for rigid-body modes and selected structural modes.

Joshi, S. M.↗

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 NTRS1997

Recovery of motion parameters from distortions in scanned images

Scanned images, such as those produced by the scanning-laser ophthalmoscope (SLO), show distortions when there is target motion. This is because pixels corresponding to different image regions are acquired sequentially, and so, in essence, are slices of different snapshots. While these distortions create problems for image registration algorithms, they are potentially useful for recovering target motion parameters at temporal frequencies above the frame rate. Stetter, Sendtner and Timberlake measured large distortions in SLO images to recover the time course of rapid horizontal saccadic eye movements. Here, this work is extended with the goal of automatically recovering small eye movements in two dimensions. Eye position during the frame interval is modeled using a low dimensional parametric description, which in turn is used to generate predicted distortions of a reference template. The input image is then registered to the distorted template using normalized cross correlation. The motion parameters are then varied, and the correlation recomputed, to find the motion which maximizes the peak value of the correlation. The location and value of the correlation maximum are determined with sub-pixel precision using biquadratic interpolation, yielding eye position resolution better than 1 arc minute. This method of motion parameter estimation is tested using actual SLO images as well as simulated images. Motion parameter estimation might also be applied to individual video lines in order to reduce pipeline delays for a near real-time system.

Mulligan, Jeffrey B.↗

NASA NTRS2005

Fundamental Mistuning Model for Probabilistic Analysis Studied Experimentally

The Fundamental Mistuning Model (FMM) is a reduced-order model for efficiently calculating the forced response of a mistuned bladed disk. FMM ID is a companion program that determines the mistuning in a particular rotor. Together, these methods provide a way to acquire mistuning data in a population of bladed disks and then simulate the forced response of the fleet. This process was tested experimentally at the NASA Glenn Research Center, and the simulated results were compared with laboratory measurements of a "fleet" of test rotors. The method was shown to work quite well. It was found that the accuracy of the results depends on two factors: (1) the quality of the statistical model used to characterize mistuning and (2) how sensitive the system is to errors in the statistical modeling.

Griffin, Jerry H.↗

NASA NTRS1975

The role of multiple scattering in one-dimensional radiative transfer

The usual methods of solving the radiative transfer equation yield answers which embrace all orders of scattering and thus shed little light on the underlying physical process. The present analysis examines the contributions of the various orders of scattering to the one dimensional transfer of radiation. In the one dimensional case an exact analytical solution exists and the problem reduces to that of expanding these exact solutions in powers of the albedo for single scattering. Formulas are given which permit the calculation of any order of scattering in an atmosphere of arbitrary optical thickness, particle albedo, and asymmetry parameter. The results should aid in identifying those physical situations where only the lowest orders of scattering play a significant role and where appropriate approximate methods might provide results of acceptable accuracy.

Adamson, D.↗

NASA NTRS2012

Space Launch System Mission Flexibility Assessment

The Space Launch System (SLS) is envisioned as a heavy lift vehicle that will provide the foundation for future beyond low Earth orbit (LEO) missions. While multiple assessments have been performed to determine the optimal configuration for the SLS, this effort was undertaken to evaluate the flexibility of various concepts for the range of missions that may be required of this system. These mission scenarios include single launch crew and/or cargo delivery to LEO, single launch cargo delivery missions to LEO in support of multi-launch mission campaigns, and single launch beyond LEO missions. Specifically, we assessed options for the single launch beyond LEO mission scenario using a variety of in-space stages and vehicle staging criteria. This was performed to determine the most flexible (and perhaps optimal) method of designing this particular type of mission. A specific mission opportunity to the Jovian system was further assessed to determine potential solutions that may meet currently envisioned mission objectives. This application sought to significantly reduce mission cost by allowing for a direct, faster transfer from Earth to Jupiter and to determine the order-of-magnitude mass margin that would be made available from utilization of the SLS. In general, smaller, existing stages provided comparable performance to larger, new stage developments when the mission scenario allowed for optimal LEO dropoff orbits (e.g. highly elliptical staging orbits). Initial results using this method with early SLS configurations and existing Upper Stages showed the potential of capturing Lunar flyby missions as well as providing significant mass delivery to a Jupiter transfer orbit.

Monk, Timothy↗

NASA NTRS1975

A steepest descents method for reentry optimization

A steepest descents optimization program is applied to the problem of a lifting vehicle entering the earth's atmosphere. The program employs penalty functions representing terminal conditions and inflight inequality constraints. During each iteration, it reduces a single performance measure which is the sum of the performance index and the penalty functions. Therefore, only one set of adjoint equations must be integrated per iteration. Values of weight factors, multiplying the penalty functions, are automatically adjusted before each iteration in order that the penalty functions will approach acceptable values. This method is shown to be a form of the classical Lagrange multiplier methods.

Andrus, J. F.↗

NASA NTRS1999

Methods of Viscosity Measurements in Sealed Ampoules

Viscosity of semiconductor and metallic melts is usually measured by oscillating cup method. This method utilizes the melts contained in vacuum sealed silica ampoules, thus the problems related to volatility, contamination, and high temperature and pressure can be alleviated. In a typical design, the time required for a single measurement is of the order of one hour. In order to reduce this time to a minute range, a high resolution (0.05 arc.sec) angular detection system is implemented in our design of the viscometer. Furthermore, an electromagnet generating a rotational magnetic field (RMF) is incorporated into the apparatus. This magnetic field can be used to remotely and non intrusively measure the electrical conductivity of the melt. It can also be used to induce a well controlled rotational flow in the system. The transient behavior of this flow can potentially yield the viscosity of the fluid. Based on RMF implementation, two novel viscometry methods are proposed in this work: a) the transient torque method, b) the resonance method. A unified theoretical approach to the three methods (oscillating cup, transient torque, and resonance) is presented along with the initial test results of the constructed apparatus. Advantages of each of the method are discussed.

Mazuruk, Konstantin↗

NASA NTRS1982

Pole placement and order reduction in two-time-scale control systems through Riccati iteration

A transformation of variables taken from singular perturbations may be applied to two-time-scale linear systems in state space form to reduce the system to block-diagonal form with slow and fast modes decoupled. The transformation is easily computed by applying the new Riccati iteration. The iteration yields a solution to the nonsymmetric algebraic Riccati equation obtained by partitioning the original system matrix A. The numerical procedure is initiated with the trivial iterate L(0) = 0, and is globally convergent to the desired unique time scale decoupling solution. After transformation, the decoupled system may be used in controller design to achieve exact closed loop pole placement in the slow subsystem without altering the poles of the fast subsystem. The decoupled form may also be used to reduce system order by wetting a small parameter to zero. Provided the fast subsystem is stable, the order reduction can be expected to yield a good approximation to the original system. These methods are demonstrated using the 16th order linear model of a turbofan engine.

Anderson, L. R.↗

NASA NTRS1995

Rainfall estimation from ESMIR-5 measurements and application to El Nino

An algorithm to estimate monthly 5 deg x 5 deg area-averaged rain rate over the oceans from January 1973 to December 1976 using single-channel microwave data from the Nimbus-5 satellite has been developed. This study extends the work of Shin et al. by including the full width of scan angles (from -50 deg to 50 deg) in order to reduce sampling error. The scan-angle dependence of the estimated rain rate due to variable antenna sidelobe effects, surface emissivity, and propagation pathlength is eliminated using a statistical method. A globally uniform beam-filling correction factor of 2.2 is applied in this study. Comparison with island station rainfall measurements over the Pacific shows a remarkably high correlation between two data in the equatorial dry zone and South Pacific convergence zone (SPCZ) but a low correlation in the extratropics and equatorial western Pacific. It is also proved that the retrieved rain rates are statistically significant. The rainfall deviations from non-El Nino years April 1973 to December 1975 reveal the temporal and spatial variations produced by the 1972-73 and 1976-77 El Nino episodes. We observe an increase of rainfall over the eastern and central equatorial Pacific Ocean and a decrease over the equatorial western Pacific Ocean and eastern Australia during these events. Consistent with previous work, the rainfall anomaly of the 1972-73 El Nino was much stronger than that of the 1976-77 El Nino.

Lim, Hyo-Suk↗

NASA NTRS1996

Near-Net Forging Technology Demonstration Program

Significant advantages in specific mechanical properties, when compared to conventional aluminum (Al) alloys, make aluminum-lithium (Al-Li) alloys attractive candidate materials for use in cryogenic propellant tanks and dry bay structures. However, the cost of Al-Li alloys is typically five times that of 2219 aluminum. If conventional fabrication processes are employed to fabricate launch vehicle structure, the material costs will restrict their utilization. In order to fully exploit the potential cost and performance benefits of Al-Li alloys, it is necessary that near-net manufacturing methods be developed to off-set or reduce raw material costs. Near-net forging is an advanced manufacturing method that uses elevated temperature metal movement (forging) to fabricate a single piece, near-net shape, structure. This process is termed 'near-net' because only a minimal amount of post-forge machining is required. The near-net forging process was developed to reduce the material scrap rate (buy-to-fly ratio) and fabrication costs associated with conventional manufacturing methods. The goal for the near-net forging process, when mature, is to achieve an overall cost reduction of approximately 50 percent compared with conventional manufacturing options for producing structures fabricated from Al-Li alloys. This NASA Marshall Space Flight Center (MSFC) sponsored program has been a part of a unique government / industry partnership, coordinated to develop and demonstrate near-net forging technology. The objective of this program was to demonstrate scale-up of the near-net forging process. This objective was successfully achieved by fabricating four integrally stiffened, 170- inch diameter by 20-inch tall, Al-Li alloy 2195, Y-ring adapters. Initially, two 2195 Al-Li ingots were converted and back extruded to produce four cylindrical blockers. Conventional ring rolling of the blockers was performed to produce ring preforms, which were then contour ring rolled to produce 'contour preforms'. All of the contour preforms on this first-of-a-kind effort were imperfect, and the ingot used to fabricate two of the preforms was of an earlier vintage. As lessons were learned throughout the program, the tooling and procedures evolved, and hence the preform quality. Two of the best contour preforms were near- net forged to produce a process pathfinder Y-ring adapter and a 'mechanical properties pathfinder' Y-ring adapter. At this point, Lockheed Martin Astronautics elected to procure additional 2195 aluminum-lithium ingot of the latest vintage, produce two additional preforms, and substitute them for older vintage material non-perfectly filled preforms already produced on this contract. The existing preforms could have been used to fulfill the requirements of the contract.

Hall, I. Keith↗

NASA NTRS1997

Simplifying silicon burning: Application of quasi-equilibrium to (alpha) network nucleosynthesis

While the need for accurate calculation of nucleosynthesis and the resulting rate of thermonuclear energy release within hydrodynamic models of stars and supernovae is clear, the computational expense of these nucleosynthesis calculations often force a compromise in accuracy to reduce the computational cost. To redress this trade-off of accuracy for speed, the authors present an improved nuclear network which takes advantage of quasi- equilibrium in order to reduce the number of independent nuclei, and hence the computational cost of nucleosynthesis, without significant reduction in accuracy. In this paper they will discuss the first application of this method, the further reduction in size of the minimal alpha network. The resultant QSE- reduced alpha network is twice as fast as the conventional alpha network it replaces and requires the tracking of half as many abundance variables, while accurately estimating the rate of energy generation. Such reduction in cost is particularly necessary for future generation of multi-dimensional models for supernovae.

Hix, W. R.↗

NASA NTRS2015

Methodology for Producing a Uniform Distribution of UO2 in a Tungsten Matrix

Current work at NASA's Marshall Space Flight Center (MSFC) is focused on the development CERMET fuel materials for Nuclear Thermal Propulsion (NTP). The CERMETs consist of uranium dioxide (UO2) fuel particles embedded in a tungsten (W) metal matrix. Initial testing of W-UO2 samples fabricated from fine angular powders performed reasonably well, but suffered from significant fuel loss during repeated thermal cycling due to agglomeration of the UO2 (1). The blended powder mixtures resulted in a non-uniform dispersion of the UO2 particles in the tungsten matrix, which allows rapid vaporization of the interconnected UO2 from the sample edges into the bulk material. Also, the angular powders create areas of stress concentrations due to thermal expansion mismatch, which eventually cracks the tungsten matrix. Evenly coating spherical UO2 particles with chemical vapor deposited (CVD) tungsten prior to consolidation was previously demonstrated to provide improved performance. However, the CVD processing technology is expensive and not currently available. In order to reduce cost and enhance performance, a powder coating process has been developed at MSFC to produce a uniform distribution of the spherical UO2 particles in a tungsten matrix. The method involves utilization of a polyethylene binder during mixing which leads to fine tungsten powders clinging to the larger UO2 spherical particles. This process was developed using HfO2 as a surrogate for UO2. Enough powder was mixed to make 8 discs (2cm diameter x 8mm thickness) using spark plasma sintering. A uniaxial pressure of 50 MPa was used at four different temperatures (2 samples at each temperature). The first two samples were heated to 1400C and 1500C respectively for 5 minutes. Densities for these samples were less than 85% of theoretical, so the time at temperature was increased to 20 minutes for the remaining samples. The highest densities were achieved for the two samples sintered at 1700C (approx. 92% of theoretical). Scanning electron microscopy (SEM) of the mixed powders and the sintered samples along with energy dispersive x-ray analysis was obtained. The SEM of the powders clearly show the fine W powder adhered to the larger HfO2 particles and a uniform distribution of HfO2 particles in a tungsten matrix upon densification. Vicker's Microhardness testing was also performed on all samples using 0.5, 1.0 and 2.0 kg loads. Five indents were made at each load level. All indents were placed in the tungsten matrix to assist as a proxy in measuring densification. The highest hardness value was obtained for the 1700C specimens. The hardness average for these samples was 312.14 MPa. This powder processing method has been applied to W/UO2 powders with the SEM of the powders appearing similar to the W/HfO2 powder images.

Tucker, Dennis S.↗

NASA NTRS2008

Accuracy of Gradient Reconstruction on Grids with High Aspect Ratio

Gradient approximation methods commonly used in unstructured-grid finite-volume schemes intended for solutions of high Reynolds number flow equations are studied comprehensively. The accuracy of gradients within cells and within faces is evaluated systematically for both node-centered and cell-centered formulations. Computational and analytical evaluations are made on a series of high-aspect-ratio grids with different primal elements, including quadrilateral, triangular, and mixed element grids, with and without random perturbations to the mesh. Both rectangular and cylindrical geometries are considered; the latter serves to study the effects of geometric curvature. The study shows that the accuracy of gradient reconstruction on high-aspect-ratio grids is determined by a combination of the grid and the solution. The contributors to the error are identified and approaches to reduce errors are given, including the addition of higher-order terms in the direction of larger mesh spacing. A parameter GAMMA characterizing accuracy on curved high-aspect-ratio grids is discussed and an approximate-mapped-least-square method using a commonly-available distance function is presented; the method provides accurate gradient reconstruction on general grids. The study is intended to be a reference guide accompanying the construction of accurate and efficient methods for high Reynolds number applications

Thomas, James↗

NASA NTRS2002

Three-Dimensional High-Order Spectral Finite Volume Method for Unstructured Grids

Many areas require a very high-order accurate numerical solution of conservation laws for complex shapes. This paper deals with the extension to three dimensions of the Spectral Finite Volume (SV) method for unstructured grids, which was developed to solve such problems. We first summarize the limitations of traditional methods such as finite-difference, and finite-volume for both structured and unstructured grids. We then describe the basic formulation of the spectral finite volume method. What distinguishes the SV method from conventional high-order finite-volume methods for unstructured triangular or tetrahedral grids is the data reconstruction. Instead of using a large stencil of neighboring cells to perform a high-order reconstruction, the stencil is constructed by partitioning each grid cell, called a spectral volume (SV), into 'structured' sub-cells, called control volumes (CVs). One can show that if all the SV cells are partitioned into polygonal or polyhedral CV sub-cells in a geometrically similar manner, the reconstructions for all the SVs become universal, irrespective of their shapes, sizes, orientations, or locations. It follows that the reconstruction is reduced to a weighted sum of unknowns involving just a few simple adds and multiplies, and those weights are universal and can be pre-determined once for all. The method is thus very efficient, accurate, and yet geometrically flexible. The most critical part of the SV method is the partitioning of the SV into CVs. In this paper we present the partitioning of a tetrahedral SV into polyhedral CVs with one free parameter for polynomial reconstructions up to degree of precision five. (Note that the order of accuracy of the method is one order higher than the reconstruction degree of precision.) The free parameter will be determined by minimizing the Lebesgue constant of the reconstruction matrix or similar criteria to obtain optimized partitions. The details of an efficient, parallelizable code to solve three-dimensional problems for any order of accuracy are then presented. Important aspects of the data structure are discussed. Comparisons with the Discontinuous Galerkin (DG) method are made. Numerical examples for wave propagation problems are presented.

Liu, Yen↗

NASA NTRS1987

Radial orbit error reduction and sea surface topography determination using satellite altimetry

A method is presented in satellite altimetry that attempts to simultaneously determine the geoid and sea surface topography with minimum wavelengths of about 500 km and to reduce the radial orbit error caused by geopotential errors. The modeling of the radial orbit error is made using the linearized Lagrangian perturbation theory. Secular and second order effects are also included. After a rather extensive validation of the linearized equations, alternative expressions of the radial orbit error are derived. Numerical estimates for the radial orbit error and geoid undulation error are computed using the differences of two geopotential models as potential coefficient errors, for a SEASAT orbit. To provide statistical estimates of the radial distances and the geoid, a covariance propagation is made based on the full geopotential covariance. Accuracy estimates for the SEASAT orbits are given which agree quite well with already published results. Observation equations are develped using sea surface heights and crossover discrepancies as observables. A minimum variance solution with prior information provides estimates of parameters representing the sea surface topography and corrections to the gravity field that is used for the orbit generation. The simulation results show that the method can be used to effectively reduce the radial orbit error and recover the sea surface topography.

Engelis, Theodossios↗

NASA NTRS2005

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