skip to main content

Explore Research Products in the NSF-PAR

Title: An evolve‐then‐filter regularized reduced order model for convection‐dominated flows
Summary

In this paper, we propose a new evolve‐then‐filter reduced order model (EF‐ROM). This is a regularized ROM (Reg‐ROM), which aims to add numerical stabilization to proper orthogonal decomposition (POD) ROMs for convection‐dominated flows. We also consider the Leray ROM (L‐ROM). These two Reg‐ROMs use explicit ROM spatial filtering to smooth (regularize) various terms in the ROMs. Two spatial filters are used: a POD projection onto a POD subspace (Proj) and a POD differential filter (DF). The four Reg‐ROM/filter combinations are tested in the numerical simulation of the three‐dimensional flow past a circular cylinder at a Reynolds numberRe=1000. Overall, the most accurate Reg‐ROM/filter combination is EF‐ROM‐DF. Furthermore, the spatial filter has a higher impact on the Reg‐ROM than the regularization used. Indeed, the DF generally yields better results than Proj for both the EF‐ROM and L‐ROM. Finally, the CPU times of the four Reg‐ROM/filter combinations are orders of magnitude lower than the CPU time of the DNS. Copyright © 2017 John Wiley & Sons, Ltd.

  more » « less
NSF-PAR ID:
10237982
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
International Journal for Numerical Methods in Fluids
Volume:
84
Issue:
10
ISSN:
0271-2091
Format(s):
Medium: X Size: p. 598-615
Size(s):
["p. 598-615"]
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ( , Fluids)
    null (Ed.)
    Developing accurate, efficient, and robust closure models is essential in the construction of reduced order models (ROMs) for realistic nonlinear systems, which generally require drastic ROM mode truncations. We propose a deep residual neural network (ResNet) closure learning framework for ROMs of nonlinear systems. The novel ResNet-ROM framework consists of two steps: (i) In the first step, we use ROM projection to filter the given nonlinear system and construct a spatially filtered ROM. This filtered ROM is low-dimensional, but is not closed. (ii) In the second step, we use ResNet to close the filtered ROM, i.e., to model the interaction between the resolved and unresolved ROM modes. We emphasize that in the new ResNet-ROM framework, data is used only to complement classical physical modeling (i.e., only in the closure modeling component), not to completely replace it. We also note that the new ResNet-ROM is built on general ideas of spatial filtering and deep learning and is independent of (restrictive) phenomenological arguments, e.g., of eddy viscosity type. The numerical experiments for the 1D Burgers equation show that the ResNet-ROM is significantly more accurate than the standard projection ROM. The new ResNet-ROM is also more accurate and significantly more efficient than other modern ROM closure models. 
    more » « less
  2. ; ( , AIAA Journal)

    Low-dimensional and computationally less-expensive reduced-order models (ROMs) have been widely used to capture the dominant behaviors of high-4dimensional systems. An ROM can be obtained, using the well-known proper orthogonal decomposition (POD), by projecting the full-order model to a subspace spanned by modal basis modes that are learned from experimental, simulated, or observational data, i.e., training data. However, the optimal basis can change with the parameter settings. When an ROM, constructed using the POD basis obtained from training data, is applied to new parameter settings, the model often lacks robustness against the change of parameters in design, control, and other real-time operation problems. This paper proposes to use regression trees on Grassmann manifold to learn the mapping between parameters and POD bases that span the low-dimensional subspaces onto which full-order models are projected. Motivated by the observation that a subspace spanned by a POD basis can be viewed as a point in the Grassmann manifold, we propose to grow a tree by repeatedly splitting the tree node to maximize the Riemannian distance between the two subspaces spanned by the predicted POD bases on the left and right daughter nodes. Five numerical examples are presented to comprehensively demonstrate the performance of the proposed method, and compare the proposed tree-based method to the existing interpolation method for POD basis and the use of global POD basis. The results show that the proposed tree-based method is capable of establishing the mapping between parameters and POD bases, and thus adapt ROMs for new parameters.

     
    more » « less
  3. ; ( , International Journal of Robust and Nonlinear Control)
    Summary

    A new approach to model order reduction of nonlinear control systems is aimed at developing persistentreduced order models(ROMs) that are robust to the changes in system's energy level. A multivariate analysis method calledsmooth orthogonal decomposition(SOD) is used to identify the dynamically relevant modal structures of the control system. The identified SOD subspaces are used to develop persistent ROMs. Performance of the resultant SOD‐based ROM is compared withproper orthogonal decomposition(POD)–based ROM by evaluating their robustness to the changes in system's energy level. Results show that SOD‐based ROMs are valid for a relatively wider range of the nonlinear control system's energy when compared with POD‐based models. In addition, the SOD‐based ROMs show considerably faster computations compared to the POD‐based ROMs of same order. For the considered dynamic system, SOD provides more effective reduction in dimension and complexity compared to POD.

     
    more » « less
  4. ; ; ( , Fluids)
    We present a stabilized POD–Galerkin reduced order method (ROM) for a Leray model. For the implementation of the model, we combine a two-step algorithm called Evolve-Filter (EF) with a computationally efficient finite volume method. In both steps of the EF algorithm, velocity and pressure fields are approximated using different POD basis and coefficients. To achieve pressure stabilization, we consider and compare two strategies: the pressure Poisson equation and the supremizer enrichment of the velocity space. We show that the evolve and filtered velocity spaces have to be enriched with the supremizer solutions related to both evolve and filter pressure fields in order to obtain stable and accurate solutions with the supremizer enrichment method. We test our ROM approach on a 2D unsteady flow past a cylinder at Reynolds number 0≤Re≤100. We find that both stabilization strategies produce comparable errors in the reconstruction of the lift and drag coefficients, with the pressure Poisson equation method being more computationally efficient. 
    more » « less
  5. ; ; ( , Fluids)
    In this paper, we introduce the evolve-then-filter (EF) regularization method for reduced order modeling of convection-dominated stochastic systems. The standard Galerkin projection reduced order model (G-ROM) yield numerical oscillations in a convection-dominated regime. The evolve-then-filter reduced order model (EF-ROM) aims at the numerical stabilization of the standard G-ROM, which uses explicit ROM spatial filter to regularize various terms in the reduced order model (ROM). Our numerical results are based on a stochastic Burgers equation with linear multiplicative noise. The numerical result shows that the EF-ROM is significantly better than G-ROM. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email