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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

A Three-Dimensional Thermodynamically Based Function for the Progressive Failure of Unidirectional Composites

A thermodynamically-based work potential theory for modeling progressive damage for laminated, unidirectional composites assuming plane stress (2D Schaperystheory) is extended to three dimensioanl (3D). An internal state variable, S, is defined to account for the dissipated energy due to damage evolution in the form of microstructure changes in the matrix. With the stationary of the total work potential with respect to the internal state variable, a thermodynamically-consistent set of evolution equationsis derived. The internal state variable is related to the transverse and shear modulithrough microdamage functions. In the first part of this work, coupon specimensare prepared to conduct experiments to characterize the relations between the internalstate variable and the transverse modulus as well as shear modulus. The informationis subsequent used for the prediction of three point bending test. In the second partof this work, objectivity is studied. Three separate methods utilizing different definitions of a reduced internal state variable or of the order of the polynomials are used to represent the matrix micro-damage functions are employed. The three methods are implemented in a user defined subroutine within a commercial finite element method software package. Results from numerical simulations of a center-notched composites panel are compared. The agreement in the maximum stress predictions among the three methods indicates that objectivity, with respect to the functional form of themicrodamage functions, is satisfied.

carbon fiber reinforced polymer composites↗
NASA NTRS2019

Cloud Giovanni: Reining in Costs and Improving Performance with Analytical Data Stores Using Scalable Serverless Architecture

Giovanni is the Geospatial Interactive Online Visualization ANd aNalysis Infrastructure developed at NASA GES DISC which provides a simple and intuitive way to visualize, analyze, and access vast amounts of Earth science data. It receives large number of user requests each day for a variety of analysis and visualization services, which leads to the big data challenge of serving gradually increasing large data volumes with diverse statistical algorithms. We hereby propose a multi-dimensional accumulation method which provides fast and cost-efficient cloud analysis for diverse services including both area averaging and time averaging. This method involves the weighted volume integration over multiple variable dimensions (time and space), and is implemented in AWS using Athena providing serverless and highly scalable data analysis. Compared to the standard method, this approach dramatically reduced the computational time by order of magnitude with a minimal AWS cost incurred. For example, for a benchmark of 10-year area averaging over the 1x1 degree daily variable, the computational time was reduced from minutes to seconds, and the Athena cost is only $5 for 100,000 requests.

Zhang, Hailiang↗
NASA NTRS2021

Design of Solar Sailing Trajectories Resilient to Safe Mode Events

Solar sails are an enabling technology for stand-alone small-satellite deep-space exploration. However, their always-on nature, combined with the time of flight required for deep-space missions, makes them particularly susceptible to safe mode events. Unlike missions using electric propulsion, exclusively solar sail- ing missions cannot carry extra propellant to make up for a safe mode event, and more powerful methods like expected availability cannot be directly applied. This work extends the expected availability and duty cycle approaches to solar sailing. Through an application to NASA’s Solar Cruiser mission, a 46 % increase in tra- jectory resilience is obtained at negligible change in the mission’s time of flight. Additionally, it is shown that when a safe mode event stops the spacecraft from reaching its target orbit, the developed methods reduce the expected final error by approximately an order of magnitude.

Johnson, Les↗
NASA NTRSDate not supplied

Flutter Prediction Report in Support of the High Angle Working Group at the Third Aeroelastic Prediction Workshop

This paper presents a summary of the computational flutter results associated with the AIAA third Aeroelastic Prediction Workshop (AePW-3) High Angle Working Group (HAWG). The computational results are compared against the experimental data collected during the Pitch and Plunge Apparatus (PAPA) Benchmark Supercritical Wing (BSCW) test campaign conducted in the Transonic Dynamics Tunnel (TDT) at NASA Langley Research Center thirty years ago. During that test several flutter points were identified at transonic conditions. One of these points, specifically near Mach 0.8 and 5° angle of attack, became a focal point of the computational challenge with in HAWG. Various fidelity time-domain, Reduced Order Model (ROM) and Linearized Frequency Domain (LFD) methods were used by seven participating teams. While there are encouraging trends in the computational results, the range of the predicted flutter dynamic pressure is still quite large due to the stall flutter mechanism. A description of each participating team’s software and methods is included.

Pawel Chwalowski↗
NASA NTRS1993

The theoretical accuracy of Runge-Kutta time discretizations for the initial boundary value problem: A careful study of the boundary error

The conventional method of imposing time dependent boundary conditions for Runge-Kutta (RK) time advancement reduces the formal accuracy of the space-time method to first order locally, and second order globally, independently of the spatial operator. This counter intuitive result is analyzed in this paper. Two methods of eliminating this problem are proposed for the linear constant coefficient case: (1) impose the exact boundary condition only at the end of the complete RK cycle, (2) impose consistent intermediate boundary conditions derived from the physical boundary condition and its derivatives. The first method, while retaining the RK accuracy in all cases, results in a scheme with much reduced CFL condition, rendering the RK scheme less attractive. The second method retains the same allowable time step as the periodic problem. However it is a general remedy only for the linear case. For non-linear hyperbolic equations the second method is effective only for for RK schemes of third order accuracy or less. Numerical studies are presented to verify the efficacy of each approach.

Carpenter, Mark H.↗
NASA NTRS1990

An efficient design sensitivity analysis of eigenvectors

Subspace iteration has been a major advance in solving large eigen problems when only a subset of eigen-pairs is required. The essence of this method is a transformation from displacement coordinates of an n-th order eigensystem to generalized coordinates of a smaller q-th order. The eigenvalue problem is then solved in the reduced space. The method was first developed by Clint and Jennings for real symmetric systems and was then called 'simultaneous iteration'. The success of the method prompted further research along this line and there have been many improved algorithms developed. This approach has been widely used by structural engineers for extracting the most useful natural frequencies and mode shapes of large-scale dynamic systems. This paper exploits into a new direction which is in the form of iterative process for simultaneously calculating eigenvector derivatives of many eigenvectors with respect to multi-variables. The method fully uses all the available information from preceding eigenvalue solution and, thus, effectively economizes computational efforts. It iterates through two equations derived from the first variation of the two fundamental equations used in subspace iteration method. There is no expensive large matrix decomposition required and the process converges to acceptable solution in a finite number of iterations. Therefore, the procedure increases its efficiency superiority over the others as the system size or the number of interested eigenvectors become larger and larger.

Ting, T.↗
NASA NTRS1990

Projection and assembly method for multibody component model reduction

The problem addressed is that of obtaining reduced-order component models for use in simulating the dynamics of a multibody system. In certain cases, nonlinear system models may be constructed using linear dynamic models for each component, but allowing large angle motion between components. Without some form of model reduction, system models constructed in this manner may be too large for use in control system design and simulation trades. This paper analyzes one method of component model reduction that allows systems level requirements (e.g., capturing the effect of body 1 reaction wheel noise on body 2 camera pointing) to aid in the selection of the reduced-order component models. Briefly stated, important modes are selected at the system level and projected onto the components, and reduced-order components are then assembled into a reduced-order system model that retains the projected modes.

Bernard, Douglas E.↗
NASA NTRS1986

On model reduction

Three model reduction methods are described. These are the discrete balanced realizations of Mullis and Roberts (1976) where a characterization of the reduction error is given and a previously unknown L(infinity) norm bound on the reduction error, is obtained. Another method is a new model reduction technique for discrete time systems which has the advantage that the reduced order model is balanced and has an L(infinity) norm bound on the reduction error. The last method derived is a frequency weighting technique for continuous and discrete systems where it is possible to specify the approximation accuracy with frequency and also, for this method, an L(infinity) norm on the weighted reduction error is obtained.

Al-Saggaf, Ubaid M.↗
NASA NTRS1989

The application of the reduced order model Kalman filter to motion estimation of degraded image sequences

Motion estimation is a field of great interest because of its many applications in areas such as robotics and image coding. The optic flow method is one such scheme which, although fairly accurate, is prone to error in the presence of noise. This thesis describes the use of the reduced order model Kalman filter (ROMKF) in reducing errors in displacement estimation due to degradation of the sequence. The implementation of filtering and motion estimation algorithms on the SUN workstation is also discussed. Results from preliminary testing were used to determine the degrees of freedom available for the ROMKF in the SUN software. The tests indicated that increasing the state to the left leads to slight improvement over the minimum state case. Therefore, the software uses the minimum model, with the option of adding states to the left only. The ROMKF was then used in conjunction with a hierarchical pel recursive motion estimation algorithm. Applying the ROMKF to the degraded displacements themselves generally yielded slight improvements in cases with noise degradation and noise plus blur. Filtering the images of the degraded sequence prior to motion estimation was less effective in these cases. Both methods performed badly in the case of blur alone, resulting in increased displacement errors. This is thought to be due in part to filter artifacts. Some improvements were obtained by varying the filter parameters when filtering the displacements directly. This result suggests that further study in varying filter parameters may lead to better results. The results of this thesis indicate that the ROMKF can play a part in reducing motion estimation errors from degraded sequences. However, more work needs to be done before the use of the ROMKF can be a practical solution.

Simpson, Elizabeth C.↗
NASA NTRS2001

Reduced-Order Modeling: New Approaches for Computational Physics

In this paper, we review the development of new reduced-order modeling techniques and discuss their applicability to various problems in computational physics. Emphasis is given to methods ba'sed on Volterra series representations and the proper orthogonal decomposition. Results are reported for different nonlinear systems to provide clear examples of the construction and use of reduced-order models, particularly in the multi-disciplinary field of computational aeroelasticity. Unsteady aerodynamic and aeroelastic behaviors of two- dimensional and three-dimensional geometries are described. Large increases in computational efficiency are obtained through the use of reduced-order models, thereby justifying the initial computational expense of constructing these models and inotivatim,- their use for multi-disciplinary design analysis.

Beran, Philip S.↗
NASA NTRS1989

Component model reduction via the projection and assembly method

The problem of acquiring a simple but sufficiently accurate model of a dynamic system is made more difficult when the dynamic system of interest is a multibody system comprised of several components. A low order system model may be created by reducing the order of the component models and making use of various available multibody dynamics programs to assemble them into a system model. The difficulty is in choosing the reduced order component models to meet system level requirements. The projection and assembly method, proposed originally by Eke, solves this difficulty by forming the full order system model, performing model reduction at the the system level using system level requirements, and then projecting the desired modes onto the components for component level model reduction. The projection and assembly method is analyzed to show the conditions under which the desired modes are captured exactly; to the numerical precision of the algorithm.

Bernard, Douglas E.↗
NASA NTRS1973

Singular perturbation equations for flexible satellites

The dynamic model of a flexible satellite with n generalized structural position coordinates requires the solution of a set of N coupled nonlinear ordinary or partial differential equations. For single composite body satellites, N is equal to (n + 3). For dual-spin systems, N is equal to (n + 9). These equations usually involve time derivatives up to the second order. For large values of n, linearization of the system has so far been the only practicable way of solution. The present study shows a method of avoiding this linearization by reducing the N equations to three or nine nonlinear, coupled, first-order ordinary differential equations involving only the angular velocities of the composite bodies. The solutions for these angular velocities lead to linear equations in the n generalized structural position coordinates, which can then be solved by known methods.

Huang, T. C.↗
NASA NTRS2024

Iterated Sigma-Point Kalman Filtering for Trajectory Reconstruction

Two second-order Kalman filtering methods are added to the New Statistical Trajec-tory Estimation Program (NewSTEP). The two filters are of the class of Sigma-Point Kalman Filters, namely the Unscented Kalman Filter and the Divided Difference Filter. Both filters are implemented with an iterative measurement update step. Results from simulation test case and a flight data test case are provided to show the filter performance. The results indicate that the new second-order filters pro-duce trajectory estimates with reduced uncertainties when compared to the existing first-order filtering method.

Sophia Vedvik↗
NASA NTRS1989

Baseline estimation with semidynamic and geometric satellite methods

The potential of a semidynamic method for an estimation of satellite position was investigated in a real data environment, using observations of Lageos collected during the MERIT campaign. In this method, quasi-simultaneous laser range observations of Lageos by pairs of stations were transformed to simultaneous range differences (SRDs), in order to reduce the effects of orbital and observational residual biases. Baselines obtained via the SRD method were compared to those obtained by the geometric method. It was found that, for baselines of regional extent, the SRD method is very efficient and at least as accurate as the more complex dynamic methods.

Dedes, George C.↗
NASA NTRSDate not supplied

Explore Information Content Efficiently from Current and Future Hyperspectral Satellite Missions using a Spectral Fingerprinting Method

Hyperspectral remote sensors from current and future missions provide measurements of the Top of Atmosphere (TOA) radiance or reflectance spectra with high information content. For example, the Atmospheric Infrared Sounder (AIRS), together with the Cross-track Infrared Sounder (CrIS), and the Infrared Atmospheric Sounding Interferometer (IASI) have provided more than 20 years radiance measurements with thousands of spectral channels. These measurements will be continued for the next two decades with the same or more advanced hyperspectral sensors. The upcoming missions such as NASA’s CLARREO Pathfinder (CPF) and ESA’s TRUTHS will provide unprecedented accurate TOA hyperspectral radiance measurements in solar spectral region. Traditional ways to derive Climate Data Records (CDRs) from these measurements are performing spatial and temporal averages of the retrieved Level-2 products. However, it is a time-consuming process to generate decades of Level-2 data from Level-1 data. Furthermore, the differences in Level-2 algorithms used for different satellite sensors will introduce errors in derived CDRs. In this presentation, we will describe a spectral fingerprinting method to generate high-quality CDRs directly from spatiotemporally averaged Level-1 data. By using consistent radiative kernels which contain the spectral information of various atmospheric and surface CDRs, we can reduce the errors due to algorithm inconsistency. Additionally, the spectral fingerprinting method reduces the time needed to generate CDRs by more than three orders of magnitude. This makes it easy to reprocess CDRs once the Level-1 data from different satellites have been improved via either re-calibrations or inter-satellite calibrations. We will present results of applying spectral fingerprinting method to 20-years of AIRS and CrIS data. The resulting CDRs include: 1) vertical profiles of atmospheric temperature and water vapor, 2) cloud properties such as optical depth, effective size, and height, 3) vertical profiles or column amounts for atmospheric trace gases such as O3 and CO, and 4) surface emissivity spectra and skin temperatures. These CDRs will be publicly available at NASA GES DISC in late 2024. The same fingerprinting method is planned to be applied to future CPF and TRUTHS data for solar spectral region.

hyperspectral remote sensing↗
NASA NTRS1976

An efficient numerical method for solving the time-dependent compressible Navier-Stokes equations at high Reynolds number

A fine-mesh method incorporating two new operators, which drastically reduces the computation time, has been developed for solving the time-dependent Navier-Stokes equations at flight Reynolds numbers. The approach time-splits the equations into a hyperbolic part and a parabolic part, solves the hyperbolic part by a new explicit numerical method based on characteristics theory, and solves the parabolic part by a new efficient implicit parabolic method. The method has reduced the computation time by one and two orders of magnitude from that required previously to solve for the interaction of a shock wave with a boundary layer on a flat plate.

Maccormack, R. W.↗
NASA NTRS1984

The efficient numerical calculation of condensational cloud drop growth

A modified Runge-Kutta integration technique is applied to condensational cloud droplet growth calculations, using the device of limiting dr, the growth of a droplet of radius r over the time step dt to some function of the absolute value of r(eq)-r, where r(eq) denotes the equilibrium value. The use of variable time steps to obtain an efficient choice of dt is also discussed. Results using these techniques are compared with those obtained with the standard fourth-order Runge-Kutta method, showing that the modified technique reduces the requisite computer time by a factor of about one hundred.

Robinson, N. F.↗
NASA NTRS1991

A hierarchy for modeling high speed propulsion systems

General research efforts on reduced order propulsion models for control systems design are overviewed. Methods for modeling high speed propulsion systems are discussed including internal flow propulsion systems that do not contain rotating machinery, such as inlets, ramjets, and scramjets. The discussion is separated into four areas: (1) computational fluid dynamics models for the entire nonlinear system or high order nonlinear models; (2) high order linearized models derived from fundamental physics; (3) low order linear models obtained from the other high order models; and (4) low order nonlinear models (order here refers to the number of dynamic states). Included in the discussion are any special considerations based on the relevant control system designs. The methods discussed are for the quasi-one-dimensional Euler equations of gasdynamic flow. The essential nonlinear features represented are large amplitude nonlinear waves, including moving normal shocks, hammershocks, simple subsonic combustion via heat addition, temperature dependent gases, detonations, and thermal choking. The report also contains a comprehensive list of papers and theses generated by this grant.

Hartley, Tom T.↗