A hybridized discontinuous Galerkin framework for high-order particle-mesh operator splitting of the incompressible Navier-Stokes equations

JOURNAL OF COMPUTATIONAL PHYSICS

卷:358

Article

Maljaars, Jakob M.¹  Labeur, Robert Jan¹  Moller, Matthias² 

[1] Delft Univ Technol, Fac Civil Engn & Geosci, Environm Fluid Mech, POB 5048,Stevinweg 1, NL-2600 GA Delft, Netherlands

[2] Delft Univ Technol, Dept Appl Math, Mekelweg 4, NL-2628 CD Delft, Netherlands

关键词: Incompressible Navier-Stokes equations;    Lagrangian-Eulerian;    Hybridized discontinuous Galerkin;    Finite elements;    Particle-in-cell;    Material point method;   

DOI  :  10.1016/j.jcp.2017.12.036

来源: Elsevier

【 摘 要 】

A generic particle-mesh method using a hybridized discontinuous Galerkin (HDG) frame-work is presented and validated for the solution of the incompressible Navier-Stokes equations. Building upon particle-in-cell concepts, the method is formulated in terms of an operator splitting technique in which Lagrangian particles are used to discretize an advection operator, and an Eulerian mesh-based HDG method is employed for the constitutive modeling to account for the inter-particle interactions. Key to the method is the variational framework provided by the HDG method. This allows to formulate the projections between the Lagrangian particle space and the Eulerian finite element space in terms of local (i.e. cellwise) l(2)-projections efficiently. Furthermore, exploiting the HDG framework for solving the constitutive equations results in velocity fields which excellently approach the incompressibility constraint in a local sense. By advecting the particles through these velocity fields, the particle distribution remains uniform over time, obviating the need for additional quality control. The presented methodology allows for a straightforward extension to arbitrary-order spatial accuracy on general meshes. A range of numerical examples shows that optimal convergence rates are obtained in space and, given the particular time stepping strategy, second-order accuracy is obtained in time. The model capabilities are further demonstrated by presenting results for the flow over a backward facing step and for the flow around a cylinder. (C) 2018 Elsevier Inc. All rights reserved.

【 授权许可】

Free