-
Approximating a branch of solutions to the Navier--Stokes equations by reduced-order modeling
Authors:
Maxim A. Olshanskii,
Leo G. Rebholz
Abstract:
This paper extends a low-rank tensor decomposition (LRTD) reduced order model (ROM) methodology to simulate viscous flows and in particular to predict a smooth branch of solutions for the incompressible Navier-Stokes equations. Additionally, it enhances the LRTD-ROM methodology by introducing a non-interpolatory variant, which demonstrates improved accuracy compared to the interpolatory method uti…
▽ More
This paper extends a low-rank tensor decomposition (LRTD) reduced order model (ROM) methodology to simulate viscous flows and in particular to predict a smooth branch of solutions for the incompressible Navier-Stokes equations. Additionally, it enhances the LRTD-ROM methodology by introducing a non-interpolatory variant, which demonstrates improved accuracy compared to the interpolatory method utilized in previous LRTD-ROM studies. After presenting both the interpolatory and non-interpolatory LRTD-ROM, we demonstrate that with snapshots from a few different viscosities, the proposed method is able to accurately predict flow statistics in the Reynolds number range $[25,400]$. This is a significantly wider and higher range than state of the art (and similar size) ROMs built for use on varying Reynolds number have been successful on. The paper also discusses how LRTD may offer new insights into the properties of parametric solutions.
△ Less
Submitted 6 May, 2024;
originally announced May 2024.
-
Analysis of a tensor POD-ROM for parameter dependent parabolic problems
Authors:
Alexander V. Mamonov,
Maxim A. Olshanskii
Abstract:
A space-time-parameters structure of the parametric parabolic PDEs motivates the application of tensor methods to define reduced order models (ROMs). Within a tensor-based ROM framework, the matrix SVD -- a traditional dimension reduction technique -- yields to a low-rank tensor decomposition (LRTD). Such tensor extension of the Galerkin proper orthogonal decomposition ROMs (POD-ROMs) benefits bot…
▽ More
A space-time-parameters structure of the parametric parabolic PDEs motivates the application of tensor methods to define reduced order models (ROMs). Within a tensor-based ROM framework, the matrix SVD -- a traditional dimension reduction technique -- yields to a low-rank tensor decomposition (LRTD). Such tensor extension of the Galerkin proper orthogonal decomposition ROMs (POD-ROMs) benefits both the practical efficiency of the ROM and its amenability for the rigorous error analysis when applied to parametric PDEs. The paper addresses the error analysis of the Galerkin LRTD-ROM for an abstract linear parabolic problem that depends on multiple physical parameters. An error estimate for the LRTD-ROM solution is proved, which is uniform with respect to problem parameters and extends to parameter values not in a sampling/training set. The estimate is given in terms of discretization and sampling mesh properties, and LRTD accuracy. The estimate depends on the smoothness rather than on the Kolmogorov n-widths of the parameterized manifold of solutions. Theoretical results are illustrated with several numerical experiments.
△ Less
Submitted 13 November, 2023;
originally announced November 2023.
-
On equilibrium states of fluid membranes
Authors:
Maxim A. Olshanskii
Abstract:
The paper studies the equilibrium configurations of inextensible elastic membranes exhibiting lateral fluidity. Using a continuum description of the membrane's motions based on the surface Navier--Stokes equations with bending forces, the paper derives differential equations governing the mechanical equilibrium. The equilibrium conditions are found to be independent of lateral viscosity and relate…
▽ More
The paper studies the equilibrium configurations of inextensible elastic membranes exhibiting lateral fluidity. Using a continuum description of the membrane's motions based on the surface Navier--Stokes equations with bending forces, the paper derives differential equations governing the mechanical equilibrium. The equilibrium conditions are found to be independent of lateral viscosity and relate tension, pressure, and tangential velocity of the fluid. These conditions suggest that either the lateral fluid motion ceases or non-decaying stationary flow of mass can only be supported by surfaces with Killing vector fields, such as axisymmetric shapes. A shape equation is derived that extends the classical Helfrich model with an area constraint to membranes of non-negligible mass. Furthermore, the paper suggests a simple numerical method to compute solutions of the shape equation. Numerical experiments conducted reveal a diverse family of equilibrium configurations. The stability of equilibrium states involving lateral flow of mass remains an unresolved question.
△ Less
Submitted 27 May, 2023; v1 submitted 29 March, 2023;
originally announced March 2023.
-
Interpolatory tensorial reduced order models for parametric dynamical systems
Authors:
Alexander V. Mamonov,
Maxim A. Olshanskii
Abstract:
The paper introduces a reduced order model (ROM) for numerical integration of a dynamical system which depends on multiple parameters. The ROM is a projection of the dynamical system on a low dimensional space that is both problem-dependent and parameter-specific. The ROM exploits compressed tensor formats to find a low rank representation for a sample of high-fidelity snapshots of the system stat…
▽ More
The paper introduces a reduced order model (ROM) for numerical integration of a dynamical system which depends on multiple parameters. The ROM is a projection of the dynamical system on a low dimensional space that is both problem-dependent and parameter-specific. The ROM exploits compressed tensor formats to find a low rank representation for a sample of high-fidelity snapshots of the system state. This tensorial representation provides ROM with an orthogonal basis in a universal space of all snapshots and encodes information about the state variation in parameter domain. During the online phase and for any incoming parameter, this information is used to find a reduced basis that spans a parameter-specific subspace in the universal space. The computational cost of the online phase then depends only on tensor compression ranks, but not on space or time resolution of high-fidelity computations. Moreover, certain compressed tensor formats enable to avoid the adverse effect of parameter space dimension on the online costs (known as the curse of dimension). The analysis of the approach includes an estimate for the representation power of the acquired ROM basis. We illustrate the performance and prediction properties of the ROM with several numerical experiments, where tensorial ROM's complexity and accuracy is compared to those of conventional POD-ROM.
△ Less
Submitted 12 May, 2022; v1 submitted 31 October, 2021;
originally announced November 2021.
-
A monolithic fluid-porous structure interaction finite element method
Authors:
Alexander Lozovskiy,
Maxim A. Olshanskii,
Yuri V. Vassilevski
Abstract:
The paper introduces a fully discrete quasi-Lagrangian finite element method for a monolithic formulation of a fluid-porous structure interaction problem. The method is second order in time and allows a standard $P_2-P_1$ (Taylor--Hood) finite element spaces for fluid problems in both fluid and porous domains. The performance of the method is illustrated on a series of numerical experiments.
The paper introduces a fully discrete quasi-Lagrangian finite element method for a monolithic formulation of a fluid-porous structure interaction problem. The method is second order in time and allows a standard $P_2-P_1$ (Taylor--Hood) finite element spaces for fluid problems in both fluid and porous domains. The performance of the method is illustrated on a series of numerical experiments.
△ Less
Submitted 12 May, 2021;
originally announced May 2021.
-
Speed-direction description of turbulent flows
Authors:
Maxim A. Olshanskii
Abstract:
In this note we introduce speed and direction variables to describe the motion of incompressible viscous flows. Fluid velocity ${\bf u}$ is decomposed into ${\bf u}=u{\bf r}$, with $u=|{\bf u}|$ and ${\bf r}={\bf u}/|{\bf u}|$. We consider a directional split of the Navier-Stokes equations into a coupled system of equations for $u$ and for ${\bf r}$. Equation for $u$ is particularly simple but sol…
▽ More
In this note we introduce speed and direction variables to describe the motion of incompressible viscous flows. Fluid velocity ${\bf u}$ is decomposed into ${\bf u}=u{\bf r}$, with $u=|{\bf u}|$ and ${\bf r}={\bf u}/|{\bf u}|$. We consider a directional split of the Navier-Stokes equations into a coupled system of equations for $u$ and for ${\bf r}$. Equation for $u$ is particularly simple but solely maintains the energy balance of the system. Under the assumption of a weak correlation between fluctuations in speed and direction in a developed turbulent flow, we further illustrate the application of $u$-${\bf r}$ variables to describe mean statistics of a shear turbulence. The standard (full) Reynolds stress tensor does not appear in a resulting equation for the mean flow profile.
△ Less
Submitted 18 November, 2020; v1 submitted 22 September, 2020;
originally announced September 2020.
-
A stable method for 4D CT-based CFD simulation in the right ventricle of a TGA patient
Authors:
Alexander Danilov,
Yushui Han,
Chun H. Lin,
Alexander Lozovskiy,
Maxim A. Olshanskii,
Victoria Yu. Salamatova,
Yuri V. Vassilevski
Abstract:
The paper discusses a stabilization of a finite element method for the equations of fluid motion in a time-dependent domain. After experimental convergence analysis, the method is applied to simulate a blood flow in the right ventricle of a post-surgery patient with the transposition of the great arteries disorder. The flow domain is reconstructed from a sequence of 4D CT images. The corresponding…
▽ More
The paper discusses a stabilization of a finite element method for the equations of fluid motion in a time-dependent domain. After experimental convergence analysis, the method is applied to simulate a blood flow in the right ventricle of a post-surgery patient with the transposition of the great arteries disorder. The flow domain is reconstructed from a sequence of 4D CT images. The corresponding segmentation and triangulation algorithms are also addressed in brief.
△ Less
Submitted 20 September, 2020;
originally announced September 2020.
-
Efficient discretizations for the EMAC formulation of the incompressible Navier-Stokes equations
Authors:
Sergey Charnyi,
Timo Heister,
Maxim A. Olshanskii,
Leo G. Rebholz
Abstract:
We study discretizations of the incompressible Navier-Stokes equations, written in the newly developed energy-momentum-angular momentum conserving (EMAC) formulation. We consider linearizations of the problem, which at each time step will reduce the computational cost, but can alter the conservation properties. We show that a skew-symmetrized linearization delivers the correct balance of (only) en…
▽ More
We study discretizations of the incompressible Navier-Stokes equations, written in the newly developed energy-momentum-angular momentum conserving (EMAC) formulation. We consider linearizations of the problem, which at each time step will reduce the computational cost, but can alter the conservation properties. We show that a skew-symmetrized linearization delivers the correct balance of (only) energy and that the Newton linearization conserves momentum and angular momentum, but conserves energy only up to the nonlinear residual. Numerical tests show that linearizing with 2 Newton steps at each time step is very effective at preserving all conservation laws at once, and giving accurate answers on long time intervals. The tests also show that the skew-symmetrized linearization is significantly less accurate. The tests also show that the Newton linearization of EMAC finite element formulation compares favorably to other traditionally used finite element formulation of the incompressible Navier-Stokes equations in primitive variables.
△ Less
Submitted 3 December, 2017;
originally announced December 2017.
-
Incompressible fluid problems on embedded surfaces: Modeling and variational formulations
Authors:
Thomas Jankuhn,
Maxim A. Olshanskii,
Arnold Reusken
Abstract:
Governing equations of motion for a viscous incompressible material surface are derived from the balance laws of continuum mechanics. The surface is treated as a time-dependent smooth orientable manifold of codimension one in an ambient Euclidian space. We use elementary tangential calculus to derive the governing equations in terms of exterior differential operators in Cartesian coordinates. The…
▽ More
Governing equations of motion for a viscous incompressible material surface are derived from the balance laws of continuum mechanics. The surface is treated as a time-dependent smooth orientable manifold of codimension one in an ambient Euclidian space. We use elementary tangential calculus to derive the governing equations in terms of exterior differential operators in Cartesian coordinates. The resulting equations can be seen as the Navier-Stokes equations posed on an evolving manifold. We consider a splitting of the surface Navier-Stokes system into coupled equations for the tangential and normal motions of the material surface. We then restrict ourselves to the case of a geometrically stationary manifold of codimension one embedded in $\Bbb{R}^n$. For this case, we present new well-posedness results for the simplified surface fluid model consisting of the surface Stokes equations. Finally, we propose and analyze several alternative variational formulations for this surface Stokes problem, including constrained and penalized formulations, which are convenient for Galerkin discretization methods.
△ Less
Submitted 9 October, 2018; v1 submitted 9 February, 2017;
originally announced February 2017.
-
A hybrid finite volume -- finite element method for bulk--surface coupled problems
Authors:
Alexey Y. Chernyshenko,
Maxim A. Olshanskii,
Yuri V. Vassilevski
Abstract:
The paper develops a hybrid method for solving a system of advection--diffusion equations in a bulk domain coupled to advection--diffusion equations on an embedded surface. A monotone nonlinear finite volume method for equations posed in the bulk is combined with a trace finite element method for equations posed on the surface. In our approach, the surface is not fitted by the mesh and is allowed…
▽ More
The paper develops a hybrid method for solving a system of advection--diffusion equations in a bulk domain coupled to advection--diffusion equations on an embedded surface. A monotone nonlinear finite volume method for equations posed in the bulk is combined with a trace finite element method for equations posed on the surface. In our approach, the surface is not fitted by the mesh and is allowed to cut through the background mesh in an arbitrary way. Moreover, a triangulation of the surface into regular shaped elements is not required. The background mesh is an octree grid with cubic cells. As an example of an application, we consider the modeling of contaminant transport in fractured porous media. One standard model leads to a coupled system of advection--diffusion equations in a bulk (matrix) and along a surface (fracture). A series of numerical experiments with both steady and unsteady problems and different embedded geometries illustrate the numerical properties of the hybrid approach. The method demonstrates great flexibility in handling curvilinear or branching lower dimensional embedded structures.
△ Less
Submitted 4 June, 2018; v1 submitted 31 December, 2016;
originally announced January 2017.
-
An adaptive numerical method for free surface flows passing rigidly mounted obstacles
Authors:
Kirill D. Nikitin,
Maxim A. Olshanskii,
Kirill M. Terekhov,
Yuri V. Vassilevski,
Ruslan Yanbarisov
Abstract:
The paper develops a method for the numerical simulation of a free-surface flow of incompressible viscous fluid around a streamlined body. The body is a rigid stationary construction partially submerged in the fluid. The application we are interested in the paper is a flow around a surface mounted offshore oil platform. The numerical method builds on a hybrid finite volume / finite difference disc…
▽ More
The paper develops a method for the numerical simulation of a free-surface flow of incompressible viscous fluid around a streamlined body. The body is a rigid stationary construction partially submerged in the fluid. The application we are interested in the paper is a flow around a surface mounted offshore oil platform. The numerical method builds on a hybrid finite volume / finite difference discretization using adaptive octree cubic meshes. The mesh is dynamically refined towards the free surface and the construction. Special care is taken to devise a discretization for the case of curvilinear boundaries and interfaces immersed in the octree Cartesian background computational mesh. To demonstrate the accuracy of the method, we show the results for two benchmark problems: the sloshing 3D container and the channel laminar flow passing the 3D cylinder of circular cross-section. Further, we simulate numerically a flow with surface waves around an offshore oil platform for the realistic set of geophysical data.
△ Less
Submitted 9 February, 2017; v1 submitted 18 September, 2016;
originally announced September 2016.
-
Finite stopping times for freely oscillating drop of a yield stress fluid
Authors:
Wanli Cheng,
Maxim A. Olshanskii
Abstract:
The paper addresses the question if there exists a finite stopping time for an unforced motion of a yield stress fluid with free surface. A variation inequality formulation is deduced for the problem of yield stress fluid dynamics with a free surface. Free surface is assumed to evolve with a normal velocity the flow. We also consider capillary forces acting along the free surface. Based on the var…
▽ More
The paper addresses the question if there exists a finite stopping time for an unforced motion of a yield stress fluid with free surface. A variation inequality formulation is deduced for the problem of yield stress fluid dynamics with a free surface. Free surface is assumed to evolve with a normal velocity the flow. We also consider capillary forces acting along the free surface. Based on the variational inequality formulation an energy equality is obtained, where kinetic and free energy rate of change is in a balance with the internal energy viscoplastic dissipation and the work of external forces. Further, the paper considers free small-amplitude oscillations of a droplet of Herschel-Bulkley fluid under the action of surface tension forces. Under certain assumptions it is shown that the finite stopping time $T_f$ of oscillations exists once the yield stress parameter is positive and the flow index $α$ satisfies ($α\ge1$). Results of several numerical experiments illustrate the analysis, reveal the dependence of $T_f$ on problem parameters and suggest an instantaneous transition of the whole drop from yielding state to the rigid one.
△ Less
Submitted 24 November, 2016; v1 submitted 31 May, 2016;
originally announced May 2016.
-
A connection between filter stabilization and eddy viscosity models
Authors:
Maxim A. Olshanskii,
Xin Xiong
Abstract:
Recently, a new approach for the stabilization of the incompressible Navier-Stokes equations for higher Reynolds numbers was introduced based on the nonlinear differential filtering of solutions on every time step of a discrete scheme. In this paper, the stabilization is shown to be equivalent to a certain eddy-viscosity model in LES. This allows a refined analysis and further understanding of des…
▽ More
Recently, a new approach for the stabilization of the incompressible Navier-Stokes equations for higher Reynolds numbers was introduced based on the nonlinear differential filtering of solutions on every time step of a discrete scheme. In this paper, the stabilization is shown to be equivalent to a certain eddy-viscosity model in LES. This allows a refined analysis and further understanding of desired filter properties. We also consider the application of the filtering in a projection (pressure correction) method, the standard splitting algorithm for time integration of the incompressible fluid equations. The paper proves an estimate on the convergence of the filtered numerical solution to the corresponding DNS solution.
△ Less
Submitted 13 April, 2013; v1 submitted 18 February, 2013;
originally announced February 2013.
-
Non-degenerate Eulerian finite element method for solving PDEs on surfaces
Authors:
Alexey Y. Chernyshenko,
Maxim A. Olshanskii
Abstract:
The paper studies a method for solving elliptic partial differential equations posed on hypersurfaces in $\mathbb{R}^N$, $N=2,3$. The method builds upon the formulation introduced in Bertalmio et al., J. Comput. Phys., 174 (2001), 759--780., where a surface equation is extended to a neighborhood of the surface. The resulting degenerate PDE is then solved in one dimension higher, but can be solved…
▽ More
The paper studies a method for solving elliptic partial differential equations posed on hypersurfaces in $\mathbb{R}^N$, $N=2,3$. The method builds upon the formulation introduced in Bertalmio et al., J. Comput. Phys., 174 (2001), 759--780., where a surface equation is extended to a neighborhood of the surface. The resulting degenerate PDE is then solved in one dimension higher, but can be solved on a mesh that is unaligned to the surface. We introduce another extended formulation, which leads to uniformly elliptic (non-degenerate) equations in a bulk domain containing the surface. We apply a finite element method to solve this extended PDE and prove the convergence of finite element solutions restricted to the surface to the solution of the original surface problem. Several numerical examples illustrate the properties of the method.
△ Less
Submitted 20 January, 2013;
originally announced January 2013.
-
A finite element solver and energy stable coupling for 3D and 1D fluid models
Authors:
Tatiana K. Dobroserdova,
Maxim A. Olshanskii
Abstract:
The paper develops a solver based on a conforming finite element method for a 3D--1D coupled incompressible flow problem. New coupling conditions are introduced to ensure a suitable bound for the cumulative energy of the model. We study the stability and accuracy of the discretization method, and the performance of some state-of-the-art linear algebraic solvers for such flow configurations. Motiva…
▽ More
The paper develops a solver based on a conforming finite element method for a 3D--1D coupled incompressible flow problem. New coupling conditions are introduced to ensure a suitable bound for the cumulative energy of the model. We study the stability and accuracy of the discretization method, and the performance of some state-of-the-art linear algebraic solvers for such flow configurations. Motivated by the simulation of the flow over inferior vena cava (IVC) filter, we consider the coupling of a 1D fluid model and a 3D fluid model posed in a domain with anisotropic inclusions. The relevance of our approach to realistic cardiovascular simulations is demonstrated by computing a blood flow over a model IVC filter.
△ Less
Submitted 20 March, 2013; v1 submitted 16 January, 2013;
originally announced January 2013.