4.7 Article

Velocity-vorticity-helicity formulation and a solver for the Navier-Stokes equations

Journal

JOURNAL OF COMPUTATIONAL PHYSICS
Volume 229, Issue 11, Pages 4291-4303

Publisher

ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/j.jcp.2010.02.012

Keywords

Incompressible Navier-Stokes equations; Finite element method; Vorticity; Helicity

Funding

  1. RAS [01.2.00104588]
  2. RFBR [08-01-00415]
  3. National Science Foundation [DMS 0914478]

Ask authors/readers for more resources

For the three-dimensional incompressible Navier-Stokes equations, we present a formulation featuring velocity, vorticity and helical density as independent variables. We find the helical density can be observed as a Lagrange multiplier corresponding to the divergence-free constraint on the vorticity variable, similar to the pressure in the case of the incompressibility condition for velocity. As one possible practical application of this new formulation, we consider a time-splitting numerical scheme based on an alternating procedure between vorticity-helical density and velocity-Bernoulli pressure systems of equations. Results of numerical experiments include a comparison with some well-known schemes based on pressure-velocity formulation and illustrate the competitiveness on the new scheme as well as the soundness of the new formulation. (C) 2010 Elsevier Inc. All rights reserved.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.7
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

Article Mathematics, Applied

INF-SUP STABILITY OF THE TRACE P2-P1 TAYLOR-HOOD ELEMENTS FOR SURFACE PDES

Maxim A. Olshanskii, Arnold Reusken, Alexander Zhiliakov

Summary: The paper introduces a trace finite element method for solving the Stokes system on a closed smooth surface. It achieves inf-sup stability through volume normal derivative stabilization and demonstrates optimal convergence and interpolation properties. Numerical examples illustrate the method's effectiveness.

MATHEMATICS OF COMPUTATION (2021)

Article Mathematics, Applied

Error analysis of higher order trace finite element methods for the surface Stokes equation

Thomas Jankuhn, Maxim A. Olshanskii, Arnold Reusken, Alexander Zhiliakov

Summary: The paper introduces a higher order unfitted finite element method for the Stokes system on a surface in R-3, utilizing parametric P-k-Pk-1 finite element pairs on a tetrahedral bulk mesh. It proves stability and optimal order convergence results, including a quantification of geometric errors from approximate parametric representation of the surface. Numerical experiments demonstrate the method's effectiveness.

JOURNAL OF NUMERICAL MATHEMATICS (2021)

Article Biochemistry & Molecular Biology

Experimental validation of a phase-field model to predict coarsening dynamics of lipid domains in multicomponent membranes

A. Zhiliakov, Y. Wang, A. Quaini, M. Olshanskii, S. Majd

Summary: Membrane phase-separation is a mechanism used by biological membranes to concentrate specific lipid species for organizing membrane processes. The computational approach based on the surface Cahn-Hilliard phase-field model complements experimental investigations in designing patchy liposomes by providing both qualitative and accurate quantitative information about membrane organization dynamics. The computational model informed by experiments has the potential to assist in designing liposomes with spatially organized surfaces, reducing cost and time in the design process.

BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES (2021)

Article Mathematics, Applied

A finite element method for Allen-Cahn equation on deforming surface

Maxim Olshanskii, Xianmin Xu, Vladimir Yushutin

Summary: The paper studies an Allen-Cahn-type equation defined on a time-dependent surface as a model of phase separation with order-disorder transition in a thin material layer. It shows that the limiting behavior of the solution is a geodesic mean curvature type flow in reference coordinates through formal inner-outer expansion. A geometrically unfitted finite element method, known as trace FEM, is considered for numerical solution with full stability and convergence analysis accounting for interpolation errors and approximate geometry recovery.

COMPUTERS & MATHEMATICS WITH APPLICATIONS (2021)

Article Engineering, Multidisciplinary

A decoupled, stable, and linear FEM for a phase-field model of variable density two-phase incompressible surface flow

Yerbol Palzhanov, Alexander Zhiliakov, Annalisa Quaini, Maxim Olshanskii

Summary: This paper presents a thermodynamically consistent phase-field model of two-phase flow of incompressible viscous fluids which allows for a non-linear dependence of the fluid density on the phase-field order parameter. An unfitted finite element method is applied to discretize the system, and a fully discrete time-stepping scheme is introduced to ensure the stability of the numerical solution. Numerical examples demonstrate the stability, accuracy, and overall efficiency of the approach, revealing interesting dependencies of flow statistics on geometry in two-phase surface flows.

COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING (2021)

Article Mathematics, Applied

An Unfitted Finite Element Method for Two-Phase Stokes Problems with Slip Between Phases

Maxim Olshanskii, Annalisa Quaini, Qi Sun

Summary: We propose an isoparametric unfitted finite element approach for simulating two-phase Stokes problems with slip, demonstrating stability and optimal error estimates independent of various factors. Numerical results in two and three dimensions support the theoretical findings and showcase the robustness of the approach.

JOURNAL OF SCIENTIFIC COMPUTING (2021)

Article Mathematics, Applied

Recycling augmented Lagrangian preconditioner in an incompressible fluid solver

Maxim A. Olshanskii, Alexander Zhiliakov

Summary: The paper discusses the reuse of matrix factorization as a building block in augmented Lagrangian and modified AL preconditioners for nonsymmetric saddle point linear algebraic systems. The strategy is applied to efficiently solve two-dimensional incompressible fluid problems independent of Reynolds number, and is tested on simulating surface fluid motion motivated by lateral fluidity of inextensible viscous membranes. New eigenvalue estimates for the AL preconditioner are derived from numerical examples including Kelvin-Helmholtz instability problems on the sphere and torus.

NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS (2022)

Article Mathematics

A Comparison of Cahn-Hilliard and Navier-Stokes-Cahn-Hilliard Models on Manifolds

Maxim Olshanskii, Yerbol Palzhanov, Annalisa Quaini

Summary: This paper investigates phase-field models for the numerical simulation of lateral phase separation and coarsening in lipid membranes, with and without lateral flow. An unfitted finite element method is used for the numerical solution of these models, allowing for complex and evolving shapes without explicit surface parametrization. The effect of lateral flow on phase evolution is examined through several numerical tests, focusing on the impact of variable line tension, viscosity, membrane composition, and surface shape on pattern formation.

VIETNAM JOURNAL OF MATHEMATICS (2022)

Article Biochemistry & Molecular Biology

Lipid domain coarsening and fluidity in multicomponent lipid vesicles: A continuum based model and its experimental validation

Y. Wang, Y. Palzhanov, A. Quaini, M. Olshanskii, S. Majd

Summary: Researchers have proposed a computational platform based on the principles of continuum mechanics and thermodynamics to model the membrane coarsening dynamics of liposomes. The platform has been quantitatively validated and shown to be a valuable tool in experimental practice.

BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES (2022)

Article Engineering, Multidisciplinary

Interpolatory tensorial reduced order models for parametric dynamical systems

Alexander V. Mamonov, Maxim A. Olshanskii

Summary: This paper presents a reduced order model (ROM) for numerical integration of a dynamical system with multiple parameters. The ROM utilizes compressed tensor formats to find a low rank representation for high-fidelity snapshots of the system state. The computational cost of the online phase depends only on tensor compression ranks, making it efficient and accurate.

COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING (2022)

Article Mathematics, Applied

A finite element scheme for the numerical solution of the Navier-Stokes/Biot coupled problem

Alexander Lozovskiy, Maxim A. Olshanskii, Yuri V. Vassilevski

Summary: This paper investigates a finite element method for a fluid-porous structure interaction problem, with a corrected balance of stresses on the fluid-structure interface. The deformations of the elastic medium are not necessarily small and are modeled using the Saint Venant-Kirchhoff (SVK) constitutive relation. The stability of the method is proved through an energy bound for the finite element solution.

RUSSIAN JOURNAL OF NUMERICAL ANALYSIS AND MATHEMATICAL MODELLING (2022)

Article Mathematics, Applied

Tangential Navier-Stokes equations on evolving surfaces: Analysis and simulations

Maxim A. Olshanskii, Arnold Reusken, Alexander Zhiliakov

Summary: This paper studies the lateral flow of a Boussinesq-Scriven fluid on a passively evolving surface embedded in Double-struck capital R-3, and introduces a well-posed weak formulation and a numerical solution method.

MATHEMATICAL MODELS & METHODS IN APPLIED SCIENCES (2022)

Article Mechanics

On equilibrium states of fluid membranes

Maxim A. A. Olshanskii

Summary: The paper explores equilibrium configurations of inextensible elastic membranes that exhibit lateral fluidity. The mechanical equilibrium is described by differential equations derived from a continuum description of membrane motions using the surface Navier-Stokes equations with bending forces. The equilibrium conditions are independent of lateral viscosity and relate tension, pressure, and tangential velocity of the fluid. Only surfaces with Killing vector fields, such as axisymmetric shapes, can support non-zero stationary mass flow.

PHYSICS OF FLUIDS (2023)

Article Mathematics, Applied

A cutFEM divergence-free discretization for the stokes problem

Haoran Liu, Michael Neilan, Maxim Olshanskii

Summary: We propose a CutFEM discretization for the Stokes problem based on the Scott-Vogelius pair, where discrete piecewise polynomial spaces are defined on non-fitted macro-element triangulations. Boundary conditions are imposed through Nitsche-type discretization, and stability is ensured by adding local ghost penalty stabilization terms. The scheme exhibits stability and a divergence-free property of the discrete velocity outside an O(h) neighborhood of the boundary. Additionally, a local grad-div stabilization is introduced to mitigate the error caused by violation of the divergence-free condition. Error analysis shows optimal order error estimates with a grad-div parameter scaling like O(h(-1)).

ESAIM-MATHEMATICAL MODELLING AND NUMERICAL ANALYSIS (2023)

Article Mechanics

Speed-direction description of turbulent flows

Maxim A. Olshanskii

PHYSICS OF FLUIDS (2020)

Article Computer Science, Interdisciplinary Applications

A new type of non-polynomial based TENO scheme for hyperbolic conservation laws

Tian Liang, Lin Fu

Summary: In this work, a new shock-capturing framework is proposed based on a new candidate stencil arrangement and the combination of infinitely differentiable non-polynomial RBF-based reconstruction in smooth regions with jump-like non-polynomial interpolation for genuine discontinuities. The resulting scheme achieves high order accuracy and resolves genuine discontinuities with sub-cell resolution.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

A high-order residual-based viscosity finite element method for incompressible variable density flow

Lukas Lundgren, Murtazo Nazarov

Summary: In this paper, a high-order accurate finite element method for incompressible variable density flow is introduced. The method addresses the issues of saddle point system and stability problem through Schur complement preconditioning and artificial compressibility approaches, and it is validated to have high-order accuracy for smooth problems and accurately resolve discontinuities.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Convergence analysis and optimization of a Robin Schwarz waveform relaxation method for time-periodic parabolic optimal control problems

Gabriele Ciaramella, Laurence Halpern, Luca Mechelli

Summary: This paper presents a novel convergence analysis of the optimized Schwarz waveform relaxation method for solving optimal control problems governed by periodic parabolic PDEs. The analysis is based on a Fourier-type technique applied to a semidiscrete-in-time form of the optimality condition, which enables a precise characterization of the convergence factor at the semidiscrete level. The behavior of the optimal transmission condition parameter is also analyzed in detail as the time discretization approaches zero.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Data-driven Whitney forms for structure-preserving control volume analysis

Jonas A. Actor, Xiaozhe Hu, Andy Huang, Scott A. Roberts, Nathaniel Trask

Summary: This article introduces a scientific machine learning framework that uses a partition of unity architecture to model physics through control volume analysis. The framework can extract reduced models from full field data while preserving the physics. It is applicable to manifolds in arbitrary dimension and has been demonstrated effective in specific problems.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Higher-continuity s-version of finite element method with B-spline functions

Nozomi Magome, Naoki Morita, Shigeki Kaneko, Naoto Mitsume

Summary: This paper proposes a novel strategy called B-spline based SFEM to fundamentally solve the problems of the conventional SFEM. It uses different basis functions and cubic B-spline basis functions with C-2-continuity to improve the accuracy of numerical integration and avoid matrix singularity. Numerical results show that the proposed method is superior to conventional methods in terms of accuracy and convergence.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

A cell-centred Eulerian volume-of-fluid method for compressible multi-material flows

Timothy R. Law, Philip T. Barton

Summary: This paper presents a practical cell-centred volume-of-fluid method for simulating compressible solid-fluid problems within a pure Eulerian setting. The method incorporates a mixed-cell update to maintain sharp interfaces, and can be easily extended to include other coupled physics. Various challenging test problems are used to validate the method, and its robustness and application in a multi-physics context are demonstrated.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Two-step multi-resolution reconstruction-based compact gas-kinetic scheme on tetrahedral mesh

Xing Ji, Fengxiang Zhao, Wei Shyy, Kun Xu

Summary: This paper presents the development of a third-order compact gas-kinetic scheme for compressible Euler and Navier-Stokes solutions, constructed particularly for an unstructured tetrahedral mesh. The scheme demonstrates robustness in high-speed flow computation and exhibits excellent adaptability to meshes with complex geometrical configurations.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Multiscale sampling for the inverse modeling of partial differential equations

Alsadig Ali, Abdullah Al-Mamun, Felipe Pereira, Arunasalam Rahunanthan

Summary: This paper presents a novel Bayesian statistical framework for the characterization of natural subsurface formations, and introduces the concept of multiscale sampling to localize the search in the stochastic space. The results show that the proposed framework performs well in solving inverse problems related to porous media flows.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Constrained optimized dynamic mode decomposition with control for physically stable systems with exogeneous inputs

Jacob Rains, Yi Wang, Alec House, Andrew L. Kaminsky, Nathan A. Tison, Vamshi M. Korivi

Summary: This paper presents a novel method called constrained optimized DMD with Control (cOptDMDc), which extends the optimized DMD method to systems with exogenous inputs and can enforce the stability of the resulting reduced order model (ROM). The proposed method optimally places eigenvalues within the stable region, thus mitigating spurious eigenvalue issues. Comparative studies show that cOptDMDc achieves high accuracy and robustness.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

A hybridizable discontinuous Galerkin formulation for the Euler-Maxwell plasma model

Andrea La Spina, Jacob Fish

Summary: This work introduces a hybridizable discontinuous Galerkin formulation for simulating ideal plasmas. The proposed method couples the fluid and electromagnetic subproblems monolithically based on source and employs a fully implicit time integration scheme. The approach also utilizes a projection-based divergence correction method to enforce the Gauss laws in challenging scenarios. Numerical examples demonstrate the high-order accuracy, efficiency, and robustness of the proposed formulation.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Numerical solution of the cavity scattering problem for flexural waves on thin plates: Linear finite element methods

Junhong Yue, Peijun Li

Summary: This paper proposes two numerical methods (IP-FEM and BP-FEM) to study the flexural wave scattering problem of an arbitrary-shaped cavity on an infinite thin plate. These methods successfully decompose the fourth-order plate wave equation into the Helmholtz and modified Helmholtz equations with coupled conditions on the cavity boundary, providing an effective solution to this challenging problem.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Fast and scalable computation of shape-morphing nonlinear solutions with application to evolutional neural networks

William Anderson, Mohammad Farazmand

Summary: We develop fast and scalable methods, called RONS, for computing reduced-order nonlinear solutions. These methods have been proven to be highly effective in tackling challenging problems, but become computationally prohibitive as the number of parameters grows. To address this issue, three separate methods are proposed and their efficacy is demonstrated through examples. The application of RONS to neural networks is also discussed.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

A second order directional split exponential integrator for systems of advection-diffusion-reaction equations

Marco Caliari, Fabio Cassini

Summary: In this paper, a second order exponential scheme for stiff evolutionary advection-diffusion-reaction equations is proposed. The scheme is based on a directional splitting approach and uses computation of small sized exponential-like functions and tensor-matrix products for efficient implementation. Numerical examples demonstrate the advantage of the proposed approach over state-of-the-art techniques.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

A conservative semi-Lagrangian method for inhomogeneous Boltzmann equation

Sebastiano Boscarino, Seung Yeon Cho, Giovanni Russo

Summary: This work proposes a high order conservative semi-Lagrangian method for the inhomogeneous Boltzmann equation of rarefied gas dynamics. The method combines a semi-Lagrangian scheme for the convection term, a fast spectral method for computation of the collision operator, and a high order conservative reconstruction and a weighted optimization technique to preserve conservative quantities. Numerical tests demonstrate the accuracy and efficiency of the proposed method.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Uniqueness and numerical scheme for spherical shell-structured sources from the far field patterns with at most two frequencies

Jialei Li, Xiaodong Liu, Qingxiang Shi

Summary: This study shows that the number, centers, scattering strengths, inner and outer diameters of spherical shell-structured sources can be uniquely determined from the far field patterns. A numerical scheme is proposed for reconstructing the spherical shell-structured sources, which includes a migration series method for locating the centers and an iterative method for computing the inner and outer diameters without computing derivatives.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)