TY - JOUR

T1 - Discrete exterior calculus discretization of incompressible Navier–Stokes equations over surface simplicial meshes

AU - Mohamed, Mamdouh S.

AU - Hirani, Anil N.

AU - Samtaney, Ravi

N1 - KAUST Repository Item: Exported on 2020-10-01

PY - 2016/2/11

Y1 - 2016/2/11

N2 - A conservative discretization of incompressible Navier–Stokes equations is developed based on discrete exterior calculus (DEC). A distinguishing feature of our method is the use of an algebraic discretization of the interior product operator and a combinatorial discretization of the wedge product. The governing equations are first rewritten using the exterior calculus notation, replacing vector calculus differential operators by the exterior derivative, Hodge star and wedge product operators. The discretization is then carried out by substituting with the corresponding discrete operators based on the DEC framework. Numerical experiments for flows over surfaces reveal a second order accuracy for the developed scheme when using structured-triangular meshes, and first order accuracy for otherwise unstructured meshes. By construction, the method is conservative in that both mass and vorticity are conserved up to machine precision. The relative error in kinetic energy for inviscid flow test cases converges in a second order fashion with both the mesh size and the time step.

AB - A conservative discretization of incompressible Navier–Stokes equations is developed based on discrete exterior calculus (DEC). A distinguishing feature of our method is the use of an algebraic discretization of the interior product operator and a combinatorial discretization of the wedge product. The governing equations are first rewritten using the exterior calculus notation, replacing vector calculus differential operators by the exterior derivative, Hodge star and wedge product operators. The discretization is then carried out by substituting with the corresponding discrete operators based on the DEC framework. Numerical experiments for flows over surfaces reveal a second order accuracy for the developed scheme when using structured-triangular meshes, and first order accuracy for otherwise unstructured meshes. By construction, the method is conservative in that both mass and vorticity are conserved up to machine precision. The relative error in kinetic energy for inviscid flow test cases converges in a second order fashion with both the mesh size and the time step.

UR - http://hdl.handle.net/10754/596175

UR - http://linkinghub.elsevier.com/retrieve/pii/S0021999116000929

UR - http://www.scopus.com/inward/record.url?scp=84959359127&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2016.02.028

DO - 10.1016/j.jcp.2016.02.028

M3 - Article

SN - 0021-9991

VL - 312

SP - 175

EP - 191

JO - Journal of Computational Physics

JF - Journal of Computational Physics

ER -