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 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS2015

Coupled Inertial Navigation and Flush Air Data Sensing Algorithm for Atmosphere Estimation

This paper describes an algorithm for atmospheric state estimation that is based on a coupling between inertial navigation and flush air data sensing pressure measurements. In this approach, the full navigation state is used in the atmospheric estimation algorithm along with the pressure measurements and a model of the surface pressure distribution to directly estimate atmospheric winds and density using a nonlinear weighted least-squares algorithm. The approach uses a high fidelity model of atmosphere stored in table-look-up form, along with simplified models of that are propagated along the trajectory within the algorithm to provide prior estimates and covariances to aid the air data state solution. Thus, the method is essentially a reduced-order Kalman filter in which the inertial states are taken from the navigation solution and atmospheric states are estimated in the filter. The algorithm is applied to data from the Mars Science Laboratory entry, descent, and landing from August 2012. Reasonable estimates of the atmosphere and winds are produced by the algorithm. The observability of winds along the trajectory are examined using an index based on the discrete-time observability Gramian and the pressure measurement sensitivity matrix. The results indicate that bank reversals are responsible for adding information content to the system. The algorithm is then applied to the design of the pressure measurement system for the Mars 2020 mission. The pressure port layout is optimized to maximize the observability of atmospheric states along the trajectory. Linear covariance analysis is performed to assess estimator performance for a given pressure measurement uncertainty. The results indicate that the new tightly-coupled estimator can produce enhanced estimates of atmospheric states when compared with existing algorithms.

Karlgaard, Christopher D.↗

NASA NTRS1991

Structural control design based on reduced-order observer

An observer-based structural control design method is proposed in this paper. The method is a semi-inverse design procedure in that the control law is not designed before the observer system, but is a result that comes from the observer design. However, the observer design is not completely independent of the control design either, but seeks to yield a control law that is close to a prescribed control law. First, the observer design problem is considered as the reconstruction of a linear function of the state vector. The linear state function to be reconstructed is the given control law. Then, based on the derivation for linear state function observers, the observer design is formulated as a parameter optimization problem. The optimization objective is to generate a matrix that is close to the optimal feedback gain matrix. Based on that matrix, the form of the observer and a new control law can be determined. The semi-inverse design procedure can yield a reduced-order observer with dimension considerably smaller than that of the system. Two examples are used to demonstrate the proposed design procedure.

Su, Tzu-Jeng↗

NASA NTRS1979

Asymptotic continuation method for trajectory optimization

A continuation method is applied to a singular perturbation parameter to obtain a new numerical method for computing optimal trajectories. This method allows one to use simply calculated reduced order approximations as starting solutions and continue the perturbation parameter until the optimal full-order solution is obtained. The method does not require the calculation of higher order correction terms nor does it require the perturbation parameter to be small - thus, it has potentially superior convergence properties compared to conventional asymptotic expansions when the perturbation parameter is large. A simple trajectory optimization problem is considered to illustrate the method.

Washburn, R. B., Jr.↗

DOE OSTI2026

Lyapunov controlled counterdiabatic quantum optimization

We introduce a quantum algorithm that integrates counterdiabatic (CD) protocols with quantum Lyapunov control (QLC) to address combinatorial optimization problems. This approach offers versatility, allowing implementation as either a digital-analog or purely digital algorithm based on selected control strategies. By examining spin-glass Hamiltonians, we illustrate how the algorithm can explore alternative paths to enhance solution outcomes compared to conventional CD techniques. This method reduces dependence on extensive higher-order CD terms and on classical optimization techniques, making it more suitable for existing quantum computing platforms. The combination of digital compression via CD protocols and the adaptable nature of QLC methods positions this approach as a promising candidate for near-term quantum devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2004

The Application of Probabilistic Methods to the Mistuning Problem

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

Griffin, J. H.↗

DOE OSTI2024

Using intrusive approaches as a step towards accounting for stochasticity in wind turbine design

Current wind turbine design methods require tens of thousands of time-domain simulations and use different random seeds to account for the stochasticity of the environmental conditions. The account of stochasticity is nonintrusive because the sampling method calls a deterministic model multiple times without changing its underlying equations. In this work, we investigate and demonstrate using simple proof of concepts how intrusive approaches can be used to directly account for stochasticity in the equations representing a mechanical system. Our long term goal is to apply such methodology to the design of wind turbines without requiring an excessive number of simulations. Intrusive methods manipulate stochastic variables directly to provide the probability density functions (PDFs) of the states and outputs at any time as functions of the PDFs of the inputs. We illustrate how different methods can be used with a reduced-order model of a wind turbine with one degree of freedom and for linear and nonlinear models. We discuss how the methods can be extended and what it will take to apply them to a level of fidelity similar to current state-of-the-art wind turbine design tools.

17 WIND ENERGY↗

DOE OSTI2022

Proposal and application of ROM-Lasso method for sensitivity coefficient evaluation

We propose a novel method for evaluating sensitivity coefficients of neutronics parameters to cross sections, so-called the reduced-order modeling technique ROM-Lasso. In this method, cross sections of interest are randomly sampled, and corresponding perturbed core analyses are performed. Then, the sensitivity coefficient vector of the higher-level model is expanded via the active subspace bases obtained with the lower-level model whose dimensional complexity is smaller than that of the higher-level model, and the expansion coefficients are estimated by the Lasso regression. A unique feature of the ROM-Lasso method allows the use of different bases optimized for each neutronics parameter. We conducted a verification calculation for an accelerator-driven system and demonstrated that the ROM-Lasso method can reproduce the sensitivity coefficients with a much smaller number of forward calculations than the direct method. The proposed method can be used to practically evaluate sensitivity coefficients. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2024

Splitting tensile strength of shale cores: intact versus fractured and sealed with ureolysis-induced calcium carbonate precipitation (UICP)

Ureolysis-induced calcium carbonate precipitation (UICP) is a biomineral solution where the urease enzyme converts urea and calcium into calcium carbonate. The resulting biomineral can bridge gaps in fractured shale, reduce undesired fluid flow, limit fracture propagation, better store carbon dioxide, and potentially enhance well efficiency. The mechanical properties of shale cores were investigated using a modified Brazilian indirect tensile strength test. An investigation of intact shale using Eagle Ford and Wolfcamp cores was conducted at varying temperatures. Results show no significant difference between shale types (average tensile strength = 6.19 MPa). Eagle Ford displayed higher strength at elevated temperature, but temperature did not influence Wolfcamp. Comparatively, cores with a single, lengthwise heterogeneous fracture were sealed with UICP and further tested for tensile strength. UICP was delivered via a flow-through method which injected 20–30 sequential patterns of ureolytic microorganisms and UICP-promoting fluids into the fracture until permeability reduced by three orders of magnitude or with an immersion method which placed cores treated with guar gum and UICP-promoting fluids into a batch reactor, demonstrating that guar gum is a suitable inclusion and may reduce the number of flow-through injections required. Tensile results for both delivery methods were variable (0.15–8 MPa), and in some cores the biomineralized fracture split apart, possibly due to insufficient sealing and/or heterogeneity in the composite UICP-shale cores. Notably in other cores the biomineralized fracture remained intact, demonstrating more cohesion than the surrounding shale, indicating that UICP may produce a strong seal for subsurface application.

58 GEOSCIENCES↗

DOE OSTI2023

Tensor-based flow reconstruction from optimally located sensor measurements

Reconstructing high-resolution flow fields from sparse measurements is a major challenge in fluid dynamics. Existing methods often vectorize the flow by stacking different spatial directions on top of each other, hence confounding the information encoded in different dimensions. Here, we introduce a tensor-based sensor placement and flow reconstruction method which retains and exploits the inherent multidimensionality of the flow. We derive estimates for the flow reconstruction error, storage requirements and computational cost of our method. We show, with examples, that our tensor-based method is significantly more accurate than similar vectorized methods. Furthermore, the variance of the error is smaller when using our tensor-based method. While the computational cost of our method is comparable to similar vectorized methods, it reduces the storage cost by several orders of magnitude. The reduced storage cost becomes even more pronounced as the dimension of the flow increases. We demonstrate the efficacy of our method on three examples: a chaotic Kolmogorov flow, in situ and satellite measurements of the global sea surface temperature and three-dimensional unsteady simulated flow around a marine research vessel.

Mechanics↗

DOE OSTI2025

Reduced basis approximations of parameterized dynamical partial differential equations via neural networks

Projection-based reduced order models are effective at approximating parameter-dependent differential equations that are parametrically separable. When parametric separability is not satisfied, which occurs in both linear and nonlinear problems, projection-based methods fail to adequately reduce the computational complexity. Devising alternative reduced order models is crucial for obtaining efficient and accurate approximations to expensive high-fidelity models. In this work, we develop a timestepping procedure for dynamical parameter-dependent problems, in which a neural-network is trained to propagate the coefficients of a reduced basis expansion. This results in an online stage with a computational cost independent of the size of the underlying problem. Here, we demonstrate our method on several parabolic partial differential equations, including a problem that is not parametrically separable.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1991

Improved large perturbation propulsion models for control system design (1988-1989) and large perturbation models of high velocity propulsion systems (1989-1990) and reduced order propulsion models for control system design (1990-1991)

Methods for modeling high speed propulsion systems will be discussed. Included in this category are internal flow propulsion systems without rotating machinery, such as inlets, ramjets, and scramjets. Among the modeling topics discussed are modeling of linear isentropic flow, heat exchange, gasdynamics, lumped parameter systems, and infinite dimensional systems. Furthermore, a generalized overview of modeling high speed propulsion systems is presented in this collection of papers.

Hartley, Tom T.↗

NASA NTRS1986

Quasi-generalized variables

The numerical solution of a system of differential and algebraic equations is difficult, due to the appearance of numerical instabilities. A method is presented here which permits numerical solutions of such a system to be obtained which satisfy the algebraic constraint equations exactly without reducing the order of the differential equations. The method is demonstrated using examples from mechanics.

Baumgarten, J.↗

NASA NTRS1980

Singular perturbation equations for flexible satellites

Force equations of motion of the individual flexible elements of a satellite were obtained in a previous paper. Moment equations of motion of the composite bodies of a flexible satellite are to be developed using two sets of equations which form the basic system for any dynamic model of flexible satellites. This basic system consists of a set of N-coupled, nonlinear, ordinary, or partial differential equations, for a flexible satellite with n generalized, structural position coordinates. For single composite body satellites, N is equal to (n + 3); for dual-spin systems, N is equal to (n + 9). These equations involve time derivatives up to the second order. The study shows a method of avoiding this linearization by reducing the N equations to 3 or 9 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 be solved by known methods.

Huang, T. C.↗

NASA NTRS1990

Reduced order component models for flexible multibody dynamics simulations

Many flexible multibody dynamics simulation codes require some form of component description that properly characterizes the dynamic behavior of the system. A model reduction procedure for producing low order component models for flexible multibody simulation is described. Referred to as projection and assembly, the method is a Rayleigh-Ritz approach that uses partitions of the system modal matrix as component Ritz transformation matrices. It is shown that the projection and assembly method yields a reduced system model that preserves a specified set of the full order system modes. Unlike classical component mode synthesis methods, the exactness of the method described is obtained at the expense of having to compute the full order system modes. The paper provides a comprehensive description of the method, a proof of exactness, and numerical results demonstrating the method's effectiveness.

Tsuha, Walter S.↗

DOE OSTI2020

The Adjoint Petrov–Galerkin method for non-linear model reduction

Here, we formulate a new projection-based reduced-order modeling technique for non-linear dynamical systems. The proposed technique, which we refer to as the Adjoint Petrov–Galerkin (APG) method, is derived by decomposing the generalized coordinates of a dynamical system into a resolved coarse-scale set and an unresolved fine-scale set. A Markovian finite memory assumption within the Mori–Zwanzig formalism is then used to develop a reduced-order representation of the coarse scales. This procedure leads to a closed reduced-order model that displays commonalities with the adjoint stabilization method used in finite elements. The formulation is shown to be equivalent to a Petrov–Galerkin method with a non-linear, time-varying test basis, thus sharing some similarities with the Least-Squares Petrov–Galerkin method. Theoretical analysis examining a priori error bounds and computational cost is presented. Numerical experiments on the compressible Navier–Stokes equations demonstrate that the proposed method can lead to improvements in numerical accuracy, robustness, and computational efficiency over the Galerkin method on problems of practical interest. Improvements in numerical accuracy and computational efficiency over the Least-Squares Petrov–Galerkin method are observed in most cases.

42 ENGINEERING↗

NASA NTRS1992

An unsymmetric Lanczos algorithm for damped structural dynamics systems

A one-sided, unsymmetric block Lanczos algorithm is proposed for the model reduction of structural dynamics systems with unsymmetric damping and/or stiffness matrices. The algorithm is a three-term iteration scheme, which transforms the system matrix into an almost skew-symmetric, block-tridiagonal form. The Lanczos reduced-order model is guaranteed to be stable if the full-order system is stable. For unstable systems, a shifting method is available. Also, the algorithm offers flexibility in the choice of starting vectors and thus can yield more accurate reduced-order models. A linear system example and a plane truss structure example are used to show the efficacy of the proposed method.

Su, Tzu-Jeng↗

NASA NTRSDate not supplied

Stochastic Reconstruction of Thermal Protection Material Properties from Arc-Jet Experiments

Material response models are used to assess reliability using variances in the bond-line temperature predictions based on uncertainties in trajectory, aerothermal environment, and material properties. A key deficiency in the current approach is that input uncertainties are too often subjective, empirical, or ad-hoc, and are not rigorously linked to the arc-jet test data used to develop the TPS material model. While materials such as PICA are well understood, future missions may require more novel materials such as HEEET where unknown uncertainties have real consequences on the ability to assess reliability. A quantifiable estimate of reliability requires an iterative methodology where the parameters driving the variance in bond-line temperature (for example) are systematically identified. A test campaign to collect data or develop new models can then be identified to reduce those input uncertainties. A Bayesian inference loop defines these connections mathematically, i.e., prior knowledge about uncertainty is updated based on observation. While these concepts are well known (and often applied intuitively in a non-rigorous approach), only recent advances in reduced-order modelling have made them computationally viable methods for engineering. By replacing deterministic inverse methods with stochastic approaches, the hope is new materials proposed for future missions can more rapidly be developed with a greater understanding of the TPS material reliability. Two additional steps for the analysis of arc jet test data are discussed. The first is ability to construct a reduced-order model using material response simulations (Icarus/US3D) of the arc-jet test articles, and the second is the inclusion of this surrogate model in the Bayesian inversion process. Both capabilities will be demonstrated using prior PICA arc-jet test data. The quality of a surrogate model will be investigated and the variances on the calibrated material properties will be compared to our current understanding of the PICA material model.

Material response↗

DOE OSTI2022

Virtual Synchronous Generator Control of Multi-port Autonomous Reconfigurable Solar Plants (MARS)

Multi-port autonomous reconfigurable solar power plant (MARS) is an integrated concept for integration of solar photovoltaic (PV) and energy storage systems (ESS) to transmission alternating current (ac) grid and high-voltage direct current (HVdc) links. The integrated development incorporates advanced control methods to provide enhanced grid ancillary services such as fast frequency responses and dynamic voltage support. In this paper, the virtual synchronous generator (VSG) control algorithm of MARS is discussed. The proposed VSG-based control enables enhanced synthetic inertial response and flexible frequency response characteristics of hybrid PV and ESS power plants in accordance with grid ancillary service requirements. Multi-port power electronics interface allows controlled emergency power support from MARS to local transmission ac grid and remote grid through the HVdc link. The performance of VSG control is validated using a reduced-order model of MARS in Simulink. Methods for estimating grid ancillary service capabilities of MARS are also discussed.

Pan, Jiuping↗