BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Traian Iliescu

Traian Iliescu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAphysics.flu-dynNumerical Analysisphysics.ao-phphysics.comp-phComputationMachine Learningmath.APmath.DSmath.OCnlin.PSphysics.geo-ph

Catalog footprint

What is connected

19works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2022arXiv

Physics Guided Machine Learning for Variational Multiscale Reduced Order Modeling

We propose a new physics guided machine learning (PGML) paradigm that leverages the variational multiscale (VMS) framework and available data to dramatically increase the accuracy of reduced order models (ROMs) at a modest computational cost. The hierarchical structure of the ROM basis and the VMS framework enable a natural separation of the resolved and unresolved ROM spatial scales. Modern PGML algorithms are used to construct novel models for the interaction among the resolved and unresolved ROM scales. Specifically, the new framework builds ROM operators that are closest to the true interaction terms in the VMS framework. Finally, machine learning is used to reduce the projection error and further increase the ROM accuracy. Our numerical experiments for a two-dimensional vorticity transport problem show that the novel PGML-VMS-ROM paradigm maintains the low computational cost of current ROMs, while significantly increasing the ROM accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Reduced Order Model Closures: A Brief Tutorial

In this paper, we present a brief tutorial on reduced order model (ROM) closures. First, we carefully motivate the need for ROM closure modeling in under-resolved simulations. Then, we construct step by step the ROM closure model by extending the classical Galerkin framework to the spaces of resolved and unresolved scales. Finally, we develop the data-driven variational multiscale ROM closure and then we test it in fluid flow simulations. Our tutorial on ROM closures is structured as a sequence of questions and answers, and is aimed at first year graduate students and advanced undergraduate students. Our goal is not to explain the "how," but the "why." That is, we carefully explain the principles used to develop ROM closures, without focusing on particular approaches. Furthermore, we try to keep the technical details to a minimum and describe the general ideas in broad terms while citing appropriate references for details.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stochastic Data-Driven Variational Multiscale Reduced Order Models

Trajectory-wise data-driven reduced order models (ROMs) tend to be sensitive to training data, and thus lack robustness. We propose to construct a robust stochastic ROM closure (S-ROM) from data consisting of multiple trajectories from random initial conditions. The S-ROM is a low-dimensional time series model for the coefficients of the dominating proper orthogonal decomposition (POD) modes inferred from data. Thus, it achieves reduction both space and time, leading to simulations orders of magnitude faster than the full order model. We show that both the estimated POD modes and parameters in the S-ROM converge when the number of trajectories increases. Thus, the S-ROM is robust when the training data size increases. We demonstrate the S-ROM on a 1D Burgers equation with a viscosity $ν= 0.002$ and with random initial conditions. The numerical results verify the convergence. Furthermore, the S-ROM makes accurate trajectory-wise predictions from new initial conditions and with a prediction time far beyond the training range, and it quantifies the spread of uncertainties due to the unresolved scales.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Multifidelity Ensemble Kalman Filter with Reduced Order Control Variates

This work develops a new multifidelity ensemble Kalman filter (MFEnKF) algorithm based on linear control variate framework. The approach allows for rigorous multifidelity extensions of the EnKF, where the uncertainty in coarser fidelities in the hierarchy of models represent control variates for the uncertainty in finer fidelities. Small ensembles of high fidelity model runs are complemented by larger ensembles of cheaper, lower fidelity runs, to obtain much improved analyses at only small additional computational costs. We investigate the use of reduced order models as coarse fidelity control variates in the MFEnKF, and provide analyses to quantify the improvements over the traditional ensemble Kalman filters. We apply these ideas to perform data assimilation with a quasi-geostrophic test problem, using direct numerical simulation and a corresponding POD-Galerkin reduced order model. Numerical results show that the two-fidelity MFEnKF provides better analyses than existing EnKF algorithms at comparable or reduced computational costs.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A long short-term memory embedding for hybrid uplifted reduced order models

In this paper, we introduce an uplifted reduced order modeling (UROM) approach through the integration of standard projection based methods with long short-term memory (LSTM) embedding. Our approach has three modeling layers or components. In the first layer, we utilize an intrusive projection approach to model dynamics represented by the largest modes. The second layer consists of an LSTM model to account for residuals beyond this truncation. This closure layer refers to the process of including the residual effect of the discarded modes into the dynamics of the largest scales. However, the feasibility of generating a low rank approximation tails off for higher Kolmogorov $n$-width systems due to the underlying nonlinear processes. The third uplifting layer, called super-resolution, addresses this limited representation issue by expanding the span into a larger number of modes utilizing the versatility of LSTM. Therefore, our model integrates a physics-based projection model with a memory embedded LSTM closure and an LSTM based super-resolution model. In several applications, we exploit the use of Grassmann manifold to construct UROM for unseen conditions. We performed numerical experiments by using the Burgers and Navier-Stokes equations with quadratic nonlinearity. Our results show robustness of the proposed approach in building reduced order models for parameterized systems and confirm the improved trade-off between accuracy and efficiency.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Data-Driven Correction Reduced Order Models for the Quasi-Geostrophic Equations: A Numerical Investigation

This paper investigates the recently introduced data-driven correction reduced order model (DDC-ROM) in the numerical simulation of the quasi-geostrophic equations. The DDC-ROM uses available data to model the correction term that is generally used to represent the missing information in low-dimensional ROMs. Physical constraints are added to the DDC-ROM to create the constrained data-driven correction reduced order model (CDDC-ROM) in order to further improve its accuracy and stability. Finally, the DDC-ROM is tested on time intervals that are longer than the time interval over which it was trained. The numerical investigation shows that, for low-dimensional ROMs, both the DDC-ROM and CDDC-ROM perform better than the standard Galerkin ROM (G-ROM) and the CDDC-ROM provides the best results.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An Evolve-Then-Filter Regularized Reduced Order Model

In this paper, we propose a new evolve-then-filter reduced order model (EF-ROM). This is a regularized ROM (Reg-ROM), which aims at the numerical stabilization of 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 new 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 number $Re=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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Approximate Deconvolution Reduced Order Modeling

This paper proposes a large eddy simulation reduced order model(LES-ROM) framework for the numerical simulation of realistic flows. In this LES-ROM framework, the proper orthogonal decomposition(POD) is used to define the ROM basis and a POD differential filter is used to define the large ROM structures. An approximate deconvolution(AD) approach is used to solve the ROM closure problem and develop a new AD-ROM. This AD-ROM is tested in the numerical simulation of the one-dimensional Burgers equation with a small diffusion coefficient(10^{-3})

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Disperse two-phase flows, with applications to geophysical problems

In this paper we study the motion of a fluid with several dispersed particles whose concentration is very small (smaller than $10^{-3}$), with possible applications to problems coming from geophysics, meteorology, and oceanography. We consider a very dilute suspension of heavy particles in a quasi-incompressible fluid (low Mach number). In our case the Stokes number is small and --as pointed out in the theory of multiphase turbulence-- we can use an Eulerian model instead of a Lagrangian one. The assumption of low concentration allows us to disregard particle--particle interactions, but we take into account the effect of particles on the fluid (two-way coupling). In this way we can study the physical effect of particle inertia (and not only passive tracers), with a model similar to the Boussinesq equations. The resulting model is used in both direct numerical simulations and large eddy simulations of a dam-break (lock-exchange) problem, which is a well-known academic test case. Keywords: Dilute suspensions, Eulerian models, direct and large eddy simulations, slightly compressible flows, dam-break (lock-exchange) problem.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Two-Level Finite Element Discretization of the Streamfunction Formulation of the Stationary Quasi-Geostrophic Equations of the Ocean

In this paper we proposed a two-level finite element discretization of the nonlinear stationary quasi-geostrophic equations, which model the wind driven large scale ocean circulation. Optimal error estimates for the two-level finite element discretization were derived. Numerical experiments for the two-level algorithm with the Argyris finite element were also carried out. The numerical results verified the theoretical error estimates and showed that, for the appropriate scaling between the coarse and fine mesh sizes, the two-level algorithm significantly decreases the computational time of the standard one-level algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Are the Snapshot Difference Quotients Needed in the Proper Orthogonal Decomposition?

This paper presents a theoretical and numerical investigation of the following practical question: Should the time difference quotients of the snapshots be used to generate the proper orthogonal decomposition basis functions? The answer to this question is important, since some published numerical studies use the time difference quotients, whereas other numerical studies do not. The criterion used in this paper to answer this question is the rate of convergence of the error of the reduced order model with respect to the number of proper orthogonal decomposition basis functions. Two cases are considered: the no_DQ case, in which the snapshot difference quotients are not used, and the DQ case, in which the snapshot difference quotients are used. The error estimates suggest that the convergence rates in the $C^0(L^2)$-norm and in the $C^0(H^1)$-norm are optimal for the DQ case, but suboptimal for the no_DQ case. The convergence rates in the $L^2(H^1)$-norm are optimal for both the DQ case and the no_DQ case. Numerical tests are conducted on the heat equation and on the Burgers equation. The numerical results support the conclusions drawn from the theoretical error estimates. Overall, the theoretical and numerical results strongly suggest that, in order to achieve optimal pointwise in time rates of convergence with respect to the number of proper orthogonal decomposition basis functions, one should use the snapshot difference quotients.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Variational Multiscale Proper Orthogonal Decomposition: Navier-Stokes Equations

We develop a variational multiscale proper orthogonal decomposition reduced-order model for turbulent incompressible Navier-Stokes equations. The error analysis of the full discretization of the model is presented. All error contributions are considered: the spatial discretization error (due to the finite element discretization), the temporal discretization error (due to the backward Euler method), and the proper orthogonal decomposition truncation error. Numerical tests for a three-dimensional turbulent flow past a cylinder at Reynolds number Re=1000 show the improved physical accuracy of the new model over the standard Galerkin and mixing-length proper orthogonal decomposition reduced-order models. The high computational efficiency of the new model is also showcased. Finally, the theoretical error estimates are confirmed by numerical simulations of a two-dimensional Navier-Stokes problem.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Finite Element Discretization of the Streamfunction Formulation of the Stationary Quasi-Geostrophic Equations of the Ocean

This paper presents a conforming finite element discretization of the streamfunction formulation of the one-layer stationary quasi-geostrophic equations, which are a commonly used model for the large scale wind- driven ocean circulation. Optimal error estimates for this finite element discretization with the Argyris element are derived. Numerical tests for the finite element discretization of the quasi-geostrophic equations and two of its standard simplifications (the linear Stommel model and the linear Stommel-Munk model) are carried out. By benchmarking the numerical results against those in the published literature, we conclude that our finite element discretization is accurate. Furthermore, the numerical results have the same convergence rates as those predicted by the theoretical error estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Approximate deconvolution large eddy simulation of a stratified two-layer quasigeostrophic ocean model

We present an approximate deconvolution (AD) large eddy simulation (LES) model for the two-layer quasigeostrophic equations. We applied the AD-LES model to mid-latitude two-layer square oceanic basins, which are standard prototypes of more realistic stratified ocean dynamics models. Two spatial filters were investigated in the AD-LES model: a tridiagonal filter and an elliptic differential filter. A sensitivity analysis of the AD-LES results with respect to changes in modeling parameters was performed. The results demonstrate that the AD-LES model used in conjunction with the tridiagonal or differential filters provides additional dissipation to the system, allowing the use of a smaller eddy viscosity coefficient. Changing the spatial filter makes a significant difference in characterizing the effective dissipation in the model. It was found that the tridiagonal filter introduces the least amount of numerical dissipation into the AD-LES model. The differential filter, however, added a significant amount of numerical dissipation to the AD-LES model for large values of the filter width. All AD-LES models reproduced the DNS results at a fraction of the cost within a reasonable level of accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Approximate deconvolution large eddy simulation of a barotropic ocean circulation model

This paper puts forth a new large eddy simulation closure modeling strategy for two-dimensional turbulent geophysical flows. This closure modeling approach utilizes approximate deconvolution, which is based solely on mathematical approximations and does not employ additional phenomenological arguments to the model. The new approximate deconvolution model is tested in the numerical simulation of the wind-driven circulation in a shallow ocean basin, a standard prototype of more realistic ocean dynamics. The model employs the barotropic vorticity equation driven by a symmetric double-gyre wind forcing, which yields a four-gyre circulation in the time mean. The approximate deconvolution model yields the correct four-gyre circulation structure predicted by a direct numerical simulation, on a coarser mesh but at a fraction of the computational cost. This first step in the numerical assessment of the new model shows that approximate deconvolution could represent a viable tool for under-resolved computations in the large eddy simulation of more realistic turbulent geophysical flows.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Proper Orthogonal Decomposition Closure Models For Turbulent Flows: A Numerical Comparison

This paper puts forth two new closure models for the proper orthogonal decomposition reduced-order modeling of structurally dominated turbulent flows: the dynamic subgrid-scale model and the variational multiscale model. These models, which are considered state-of-the-art in large eddy simulation, together with the mixing length and the Smagorinsky closure models, are tested in the numerical simulation of a 3D turbulent flow around a circular cylinder at Re = 1,000. Two criteria are used in judging the performance of the proper orthogonal decomposition reduced-order models: the kinetic energy spectrum and the time evolution of the POD coefficients. All the numerical results are benchmarked against a direct numerical simulation. Based on these numerical results, we conclude that the dynamic subgrid-scale and the variational multiscale models perform best.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Variational Multiscale Proper Orthogonal Decomposition: Convection-Dominated Convection-Diffusion Equations

We introduce a variational multiscale closure modeling strategy for the numerical stabilization of proper orthogonal decomposition reduced-order models of convection-dominated equations. As a first step, the new model is analyzed and tested for convection-dominated convection-diffusion equations. The numerical analysis of the finite element discretization of the model is presented. Numerical tests show the increased numerical accuracy over the standard reduced-order model and illustrate the theoretical convergence rates.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Impact of Boundary Conditions on Entrainment and Transport in Gravity Currents

Gravity currents have been studied numerically and experimentally both in the laboratory and in the ocean. The question of appropriate boundary conditions is still challenging for most complex flows. Gravity currents make no exception - appropriate, physically and mathematically sound boundary conditions are yet to be found. This task is further complicated by the technical limitations imposed by the current oceanographic techniques. In this paper, we make a first step toward a better understanding of the impact of boundary conditions on gravity currents. Specifically, we use direct numerical simulations to investigate the effect that the popular Neumann, and less popular Dirichlet boundary conditions on the bottom continental shelf have on the entrainment and transport of gravity currents. The finding is that gravity currents under these two different boundary conditions differ most in the way they transport heat from the top towards the bottom. This major difference occurs at medium temperature ranges. Entrainment and transport at high temperatures also show significant differences.

0 citations0 reviewsNo code linkedOpen access