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 163 records · Page 9

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

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 sections: (1) computational fluid dynamics model for the entire nonlinear system or high order nonlinear models; (2) high order linearized model derived from fundamental physics; (3) low order linear models obtained from other high order models; and (4) low order nonlinear models. Included are special considerations on any 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, moving normal shocks, hammershocks, subsonic combustion via heat addition, temperature dependent gases, detonation, and thermal choking.

Hartley, Tom T.↗

A decentralized linear quadratic control design method for flexible structures

A decentralized suboptimal linear quadratic control design procedure which combines substructural synthesis, model reduction, decentralized control design, subcontroller synthesis, and controller reduction is proposed for the design of reduced-order controllers for flexible structures. The procedure starts with a definition of the continuum structure to be controlled. An evaluation model of finite dimension is obtained by the finite element method. Then, the finite element model is decomposed into several substructures by using a natural decomposition called substructuring decomposition. Each substructure, at this point, still has too large a dimension and must be reduced to a size that is Riccati-solvable. Model reduction of each substructure can be performed by using any existing model reduction method, e.g., modal truncation, balanced reduction, Krylov model reduction, or mixed-mode method. Then, based on the reduced substructure model, a subcontroller is designed by an LQ optimal control method for each substructure independently. After all subcontrollers are designed, a controller synthesis method called substructural controller synthesis is employed to synthesize all subcontrollers into a global controller. The assembling scheme used is the same as that employed for the structure matrices. Finally, a controller reduction scheme, called the equivalent impulse response energy controller (EIREC) reduction algorithm, is used to reduce the global controller to a reasonable size for implementation. The EIREC reduced controller preserves the impulse response energy of the full-order controller and has the property of matching low-frequency moments and low-frequency power moments. An advantage of the substructural controller synthesis method is that it relieves the computational burden associated with dimensionality. Besides that, the SCS design scheme is also a highly adaptable controller synthesis method for structures with varying configuration, or varying mass and stiffness properties.

Su, Tzu-Jeng↗

Stability bounds for the control of large space structures

Balas (1977) has discussed the stability problem of reduced-order regulators and estimators in terms of control and observation 'spillover'. The term 'control spillover' was used to define that part of the feedback control which excites the uncontrolled (or residual) modes, and 'observation spillover' was used to define that part of the measurement which is contaminated by residual modes. In this paper, two sufficient conditions are derived via Lyapunov methods for asymptotic stability of large space structures using a class of reduced-order controllers. These conditions give allowable bounds on the spectral norms of control and observation spillover terms. The sufficient condition given by a specified inequality equation appears to be less conservative, and should be useful as a design tool for the control of large space structures.

Joshi, S. M.↗

Achieving High Efficiency in Reduced Order Modeling for Large Scale Polycrystal Plasticity Simulations

Reduced order models for the nonlinear response of heterogeneous microstructures typically require a construction (or training) stage to build the reduced order basis. In this manuscript, an efficient model construction strategy for the eigenstrain homogenization method (EHM) is presented. The proposed strategy relies on a parallel, element-by-element, conjugate gradient solver. Near linear scaling has been achieved with respect to the number of degrees of freedom used to resolve the microstructure. Linear scaling with respect to the number of pre-analyses required to construct the reduced order model (ROM) follows from the EHM formulation. Furthermore, a parallel implementation for fast evaluation of the constructed ROM has been developed using shared memory parallelization. It has been shown that for large microstructures with ≈ 10,000 grains, the total computational cost of evaluating the nonlinear response of a polycrystal could be reduced by approximately an order of magnitude using 32 cores with respect to serial ROM simulation. The present methodology has been verified using an additively manufactured polycrystalline microstructure of a nickel-based superalloy, Inconel 625. The capability of the developed framework to construct a ROM for such large microstructures, as well as the ability of the ROM to predict average and local quantities of interest has been demonstrated.

microscale↗

A method for high order linear system reduction and nonlinear system simplification

Least-squares-type algorithms for reducing the order of linear systems in the frequency domain and simplifying nonlinear systems in time domain are developed and demonstrated. The possible model structures are represented as nodes in a tree, and costs along the branches are assigned using the repeated-Gram-Schmidt orthogonalization procedure of Desrochers and Saridis (1980), permitting identification of the optimal n-term model by searching the tree to depth n, with no need for parameter identification. The efficiency and flexibility of the algorithms is shown in applications to the eighth-order linear system studied by Hsia (1972), a three-state eight-nonlinear-term aircraft-dynamics problem, and the related linear-controller problem (Garrard and Jordan, 1977).

Desrochers, A. A.↗

A Reduced-Order Model for Efficient Simulation of Synthetic Jet Actuators

A new reduced-order model of multidimensional synthetic jet actuators that combines the accuracy and conservation properties of full numerical simulation methods with the efficiency of simplified zero-order models is proposed. The multidimensional actuator is simulated by solving the time-dependent compressible quasi-1-D Euler equations, while the diaphragm is modeled as a moving boundary. The governing equations are approximated with a fourth-order finite difference scheme on a moving mesh such that one of the mesh boundaries coincides with the diaphragm. The reduced-order model of the actuator has several advantages. In contrast to the 3-D models, this approach provides conservation of mass, momentum, and energy. Furthermore, the new method is computationally much more efficient than the multidimensional Navier-Stokes simulation of the actuator cavity flow, while providing practically the same accuracy in the exterior flowfield. The most distinctive feature of the present model is its ability to predict the resonance characteristics of synthetic jet actuators; this is not practical when using the 3-D models because of the computational cost involved. Numerical results demonstrating the accuracy of the new reduced-order model and its limitations are presented.

Yamaleev, Nail K.↗

Application of Transfer Matrix Approach to Modeling and Decentralized Control of Lattice-Based Structures

This paper presents a modeling and control of aerostructure developed by lattice-based cellular materials/components. The proposed aerostructure concept leverages a building block strategy for lattice-based components which provide great adaptability to varying ight scenarios, the needs of which are essential for in- ight wing shaping control. A decentralized structural control design is proposed that utilizes discrete-time lumped mass transfer matrix method (DT-LM-TMM). The objective is to develop an e ective reduced order model through DT-LM-TMM that can be used to design a decentralized controller for the structural control of a wing. The proposed approach developed in this paper shows that, as far as the performance of overall structural system is concerned, the reduced order model can be as e ective as the full order model in designing an optimal stabilizing controller.

Transfer Matrix Method↗

Reduced-Order Kalman Filtering for Processing Relative Measurements

A study in Kalman-filter theory has led to a method of processing relative measurements to estimate the current state of a physical system, using less computation than has previously been thought necessary. As used here, relative measurements signifies measurements that yield information on the relationship between a later and an earlier state of the system. An important example of relative measurements arises in computer vision: Information on relative motion is extracted by comparing images taken at two different times. Relative measurements do not directly fit into standard Kalman filter theory, in which measurements are restricted to those indicative of only the current state of the system. One approach heretofore followed in utilizing relative measurements in Kalman filtering, denoted state augmentation, involves augmenting the state of the system at the earlier of two time instants and then propagating the state to the later time instant.While state augmentation is conceptually simple, it can also be computationally prohibitive because it doubles the number of states in the Kalman filter. When processing a relative measurement, if one were to follow the state-augmentation approach as practiced heretofore, one would find it necessary to propagate the full augmented state Kalman filter from the earlier time to the later time and then select out the reduced-order components. The main result of the study reported here is proof of a property called reduced-order equivalence (ROE). The main consequence of ROE is that it is not necessary to augment with the full state, but, rather, only the portion of the state that is explicitly used in the partial relative measurement. In other words, it suffices to select the reduced-order components first and then propagate the partial augmented state Kalman filter from the earlier time to the later time; the amount of computation needed to do this can be substantially less than that needed for propagating the full augmented Kalman state filter.

Bayard, David S.↗

Some Aspects of Essentially Nonoscillatory (ENO) Formulations for the Euler Equations, Part 3

An essentially nonoscillatory (ENO) formulation is described for hyperbolic systems of conservation laws. ENO approaches are based on smart interpolation to avoid spurious numerical oscillations. ENO schemes are a superset of Total Variation Diminishing (TVD) schemes. In the recent past, TVD formulations were used to construct shock capturing finite difference methods. At extremum points of the solution, TVD schemes automatically reduce to being first-order accurate discretizations locally, while away from extrema they can be constructed to be of higher order accuracy. The new framework helps construct essentially non-oscillatory finite difference methods without recourse to local reductions of accuracy to first order. Thus arbitrarily high orders of accuracy can be obtained. The basic general ideas of the new approach can be specialized in several ways and one specific implementation is described based on: (1) the integral form of the conservation laws; (2) reconstruction based on the primitive functions; (3) extension to multiple dimensions in a tensor product fashion; and (4) Runge-Kutta time integration. The resulting method is fourth-order accurate in time and space and is applicable to uniform Cartesian grids. The construction of such schemes for scalar equations and systems in one and two space dimensions is described along with several examples which illustrate interesting aspects of the new approach.

Chakravarthy, Sukumar R.↗

Active Control of High Frequency Combustion Instability in Aircraft Gas-Turbine Engines

Active control of high-frequency (greater than 500 Hz) combustion instability has been demonstrated in the NASA single-nozzle combustor rig at United Technologies Research Center. The combustor rig emulates an actual engine instability and has many of the complexities of a real engine combustor (i.e. actual fuel nozzle and swirler, dilution cooling, etc.) In order to demonstrate control, a high-frequency fuel valve capable of modulating the fuel flow at up to 1kHz was developed. Characterization of the fuel delivery system was accomplished in a custom dynamic flow rig developed for that purpose. Two instability control methods, one model-based and one based on adaptive phase-shifting, were developed and evaluated against reduced order models and a Sectored-1-dimensional model of the combustor rig. Open-loop fuel modulation testing in the rig demonstrated sufficient fuel modulation authority to proceed with closed-loop testing. During closed-loop testing, both control methods were able to identify the instability from the background noise and were shown to reduce the pressure oscillations at the instability frequency by 30%. This is the first known successful demonstration of high-frequency combustion instability suppression in a realistic aero-engine environment. Future plans are to carry these technologies forward to demonstration on an advanced low-emission combustor.

Corrigan, Bob↗

Alternative Modal Basis Selection Procedures for Nonlinear Random Response Simulation

Three procedures to guide selection of an efficient modal basis in a nonlinear random response analysis are examined. One method is based only on proper orthogonal decomposition, while the other two additionally involve smooth orthogonal decomposition. Acoustic random response problems are employed to assess the performance of the three modal basis selection approaches. A thermally post-buckled beam exhibiting snap-through behavior, a shallowly curved arch in the auto-parametric response regime and a plate structure are used as numerical test articles. The results of the three reduced-order analyses are compared with the results of the computationally taxing simulation in the physical degrees of freedom. For the cases considered, all three methods are shown to produce modal bases resulting in accurate and computationally efficient reduced-order nonlinear simulations.

Przekop, Adam↗

Calculations of Laminar Heat Transfer Around Cylinders of Arbitrary Cross Section and Transpiration-Cooled Walls with Application to Turbine Blade Cooling

An approximate method for the development of flow and thermal boundary layers in the laminar region on cylinders with arbitrary cross section and transpiration-cooled walls is obtained by the use of Karman's integrated momentum equation and an analogous heat-flow equation. Incompressible flow with constant property values throughout the boundary layer is assumed. The velocity and temperature profiles within the boundary layer are approximated by expressions composed of trigonometric functions. Shape parameters for these profiles and functions necessary for the solution of the boundary-layer equations are presented as graphs so that the calculation for any specific case is reduced to the solution of two first-order differential equations. This method is applied to determine local heat-transfer coefficients and surface temperatures in the laminar flow region of the transpiration-cooled turbine blades for a given coolant flow rate, or to calculate the coolant flow distribution which is necessary in order to keep the blade temperature uniform along the surface.

Eckert, E R G↗

Stabilization by modification of the Lagrangian

In order to reduce the error growth during a numerical integration, a method of stabilization of the differential equations of the Keplerian motion is offered. It is characterized by the use of the eccentric anomaly as an independent variable in such a way that the time transformation is given by a generalized Lagrange formalism. The control terms in the equations of motion obtained by this modified Lagrangian give immediately a completely Liapunov-stable set of differential equations. In contrast to other publications, here the equation of time integration is modified by a control term which leads to an integral which defined the time element for the perturbed Keplerian motion.

Baumgarte, J. W.↗

Stabilization by modification of the Lagrangian

In order to reduce the error growth during a numerical integration, a method for stabilizing the differential equations of Keplerian motion is offered. It is characterized by the use of the eccentric anomaly as independent variable in such a way that the time transformation is given by a generalized Lagrange formalism. The control terms in the equations of motion obtained by this modified Lagrangian give immediately a completely Lyapunov-stable set of differential equations. The equation of time integration is modified by a control term which leads to an integral which defined the time element for the perturbed Keplerian motion.

Baumgarte, J. W.↗

A case for nonlinear model simplification in the design of flight control systems

An algorithm is presented for obtaining a reduced order or simplified nonlinear model of a nonlinear plant. The method provides a systematic approach for quickly finding the optimal simplified model of the system. Controllers are then designed around this nonlinear model and are used to control the original plant. An example design, for the F-8 aircraft, illustrates that the simplified model does indeed enable the design of a linear controller whose performance is indistinguishable from higher order optimal controllers for the original nonlinear plant. The net result is a simplified design procedure for effectively controlling complex nonlinear systems.

Desrochers, A. A.↗

Some applicatiaons of Lanczos vectors in structural dynamics

This paper summarizes several applications of Lanczos vectors and Krylov vectors. Lanczos vectors or Krylov vectors for use in component synthesis, Lanczos vectors for systems with unsymmetric damping (unsymmetric Block-Lanczos method), and the use of Lanczos vectors to obtain reduced-order models for control of flexible structures. For each of the above applications the theoretical background is briefly summarized and a numerical example is given.

Craig, Roy R., Jr.↗

Backscattering from a randomly rough dielectric surface

A backscattering model for scattering from a randomly rough dielectric surface is developed based on an approximate solution of a pair of integral equations for the tangential surface fields. Both like and cross-polarized scattering coefficients are obtained. It is found that the like polarized scattering coefficients contain two types of terms: single scattering terms and multiple scattering terms. The single scattering terms in like polarized scattering are shown to reduce the first-order solutions derived from the small perturbation method when the roughness parameters satisfy the slightly rough conditions. When surface roughnesses are large but the surface slope is small, only a single scattering term corresponding to the standard Kirchhoff model is significant. If the surface slope is large, the multiple scattering term will also be significant. The cross-polarized backscattering coefficients satisfy reciprocity and contain only multiple scattering terms. The difference between vertical and horizontal scattering coefficients is found to increase with the dielectric constant and is generally smaller than that predicted by the first-order small perturbation model. Good agreements are obtained between this model and measurements from statistically known surfaces.

Fung, Adrian K.↗