95 research outputs found
The ROMES method for statistical modeling of reduced-order-model error
This work presents a technique for statistically modeling errors introduced
by reduced-order models. The method employs Gaussian-process regression to
construct a mapping from a small number of computationally inexpensive `error
indicators' to a distribution over the true error. The variance of this
distribution can be interpreted as the (epistemic) uncertainty introduced by
the reduced-order model. To model normed errors, the method employs existing
rigorous error bounds and residual norms as indicators; numerical experiments
show that the method leads to a near-optimal expected effectivity in contrast
to typical error bounds. To model errors in general outputs, the method uses
dual-weighted residuals---which are amenable to uncertainty control---as
indicators. Experiments illustrate that correcting the reduced-order-model
output with this surrogate can improve prediction accuracy by an order of
magnitude; this contrasts with existing `multifidelity correction' approaches,
which often fail for reduced-order models and suffer from the curse of
dimensionality. The proposed error surrogates also lead to a notion of
`probabilistic rigor', i.e., the surrogate bounds the error with specified
probability
Model Reduction for Multiscale Lithium-Ion Battery Simulation
In this contribution we are concerned with efficient model reduction for
multiscale problems arising in lithium-ion battery modeling with spatially
resolved porous electrodes. We present new results on the application of the
reduced basis method to the resulting instationary 3D battery model that
involves strong non-linearities due to Buttler-Volmer kinetics. Empirical
operator interpolation is used to efficiently deal with this issue.
Furthermore, we present the localized reduced basis multiscale method for
parabolic problems applied to a thermal model of batteries with resolved porous
electrodes. Numerical experiments are given that demonstrate the reduction
capabilities of the presented approaches for these real world applications
Model order reduction approaches for infinite horizon optimal control problems via the HJB equation
We investigate feedback control for infinite horizon optimal control problems
for partial differential equations. The method is based on the coupling between
Hamilton-Jacobi-Bellman (HJB) equations and model reduction techniques. It is
well-known that HJB equations suffer the so called curse of dimensionality and,
therefore, a reduction of the dimension of the system is mandatory. In this
report we focus on the infinite horizon optimal control problem with quadratic
cost functionals. We compare several model reduction methods such as Proper
Orthogonal Decomposition, Balanced Truncation and a new algebraic Riccati
equation based approach. Finally, we present numerical examples and discuss
several features of the different methods analyzing advantages and
disadvantages of the reduction methods
The GNAT method for nonlinear model reduction: effective implementation and application to computational fluid dynamics and turbulent flows
The Gauss--Newton with approximated tensors (GNAT) method is a nonlinear
model reduction method that operates on fully discretized computational models.
It achieves dimension reduction by a Petrov--Galerkin projection associated
with residual minimization; it delivers computational efficency by a
hyper-reduction procedure based on the `gappy POD' technique. Originally
presented in Ref. [1], where it was applied to implicit nonlinear
structural-dynamics models, this method is further developed here and applied
to the solution of a benchmark turbulent viscous flow problem. To begin, this
paper develops global state-space error bounds that justify the method's design
and highlight its advantages in terms of minimizing components of these error
bounds. Next, the paper introduces a `sample mesh' concept that enables a
distributed, computationally efficient implementation of the GNAT method in
finite-volume-based computational-fluid-dynamics (CFD) codes. The suitability
of GNAT for parameterized problems is highlighted with the solution of an
academic problem featuring moving discontinuities. Finally, the capability of
this method to reduce by orders of magnitude the core-hours required for
large-scale CFD computations, while preserving accuracy, is demonstrated with
the simulation of turbulent flow over the Ahmed body. For an instance of this
benchmark problem with over 17 million degrees of freedom, GNAT outperforms
several other nonlinear model-reduction methods, reduces the required
computational resources by more than two orders of magnitude, and delivers a
solution that differs by less than 1% from its high-dimensional counterpart
POD\u2013Galerkin reduced order methods for combined Navier\u2013Stokes transport equations based on a hybrid FV-FE solver
The purpose of this work is to introduce a novel POD\u2013Galerkin strategy for the semi-implicit hybrid high order finite volume/finite element solver introduced in Berm\ufadez et al. (2014) and Busto et al. (2018). The interest is into the incompressible Navier\u2013Stokes equations coupled with an additional transport equation. The full order model employed in this article makes use of staggered meshes. This feature will be conveyed to the reduced order model leading to the definition of reduced basis spaces in both meshes. The reduced order model presented herein accounts for velocity, pressure, and a transport-related variable. The pressure term at both the full order and the reduced order level is reconstructed making use of a projection method. More precisely, a Poisson equation for pressure is considered within the reduced order model. Results are verified against three-dimensional manufactured test cases. Moreover a modified version of the classical cavity test benchmark including the transport of a species is analysed
Non-intrusive polynomial chaos method applied to full-order and reduced problems in computational fluid dynamics: A comparison and perspectives
In this work, Uncertainty Quantification (UQ) based on non-intrusive Polynomial Chaos Expansion (PCE) is applied to the CFD problem of the flow past an airfoil with parameterized angle of attack and inflow velocity. To limit the computational cost associated with each of the simulations required by the non-intrusive UQ algorithm used, we resort to a Reduced Order Model (ROM) based on Proper Orthogonal Decomposition (POD)-Galerkin approach. A first set of results is presented to characterize the accuracy of the POD-Galerkin ROM developed approach with respect to the Full Order Model (FOM) solver (OpenFOAM). A further analysis is then presented to assess how the UQ results are affected by substituting the FOM predictions with the surrogate ROM ones
Finite volume POD-Galerkin stabilised reduced order methods for the parametrised incompressible Navier\u2013Stokes equations
In this work a stabilised and reduced Galerkin projection of the incompressible unsteady Navier\u2013Stokes equations for moderate Reynolds number is presented. The full-order model, on which the Galerkin projection is applied, is based on a finite volumes approximation. The reduced basis spaces are constructed with a POD approach. Two different pressure stabilisation strategies are proposed and compared: the former one is based on the supremizer enrichment of the velocity space, and the latter one is based on a pressure Poisson equation approach
Model Reduction Opportunities in Detailed Simulations of Combustion Dynamics
Rocket and gas turbine combustion dynamics involves a confluence of diverse physics and interaction across a number of system components. Any comprehensive, self-consistent numerical model is burdened by a very large computational mesh, stiff unsteady processes which limit the permissible time step, and the need to perform tedious, repeated calculations for a broad parametric range. Predictive CFD models rely on very large scale simulations and advanced hardware. Reduced Basis Methods (RBM) have grown in usage during the past decade, as promising new techniques in making large simulations more accessible. These methods create models with far fewer unknown quantities than the original system, by generating “proper” fundamental solutions and their Galerkin projections, while guaranteeing accuracy and computational efficiency. RBMs seek to reproduce full CFD solutions, rather than solutions to a simplified or linearized set of equations. We present here some recent work in this area, focusing on approaches to model large scale combustor systems. The maturation of methods leading to LES-based turbulent combustion modeling is discussed, and model reduction goals and strategies are explored from the perspective of applicability in real life problems in both gas turbine, as well as rocket engines
- …