arrow
Return

A second order virtual node algorithm for Navier-Stokes flow problems with interfacial forces and discontinuous material properties

delete2014-05-01
delete14
PRE
AI
C
Craig Schroeder *
A
Alexey Stomakhin
R
Russell Howes
J
Joseph Teran
DOI:10.1016/j.jcp.2014.01.051delete
deleteOriginal
deleteOriginal request for help
deleteShare
deleteSave
Abstract

Abstract

En 中文
We present a numerical method for the solution of the Navier-Stokes equations in three dimensions that handles interfacial discontinuities due to singular forces and discontinuous fluid properties such as viscosity and density. We show that this also allows for the enforcement of normal stress and velocity boundary conditions on irregular domains. The method improves on results in [1] (which solved the Stokes equations in two dimensions) by providing treatment of fluid inertia as well as a new discretization of jump and boundary conditions that accurately resolves null modes in both two and three dimensions. We discretize the equations using an embedded approach on a uniform MAC grid to yield discretely divergence-free velocities that are second order accurate. We maintain our interface using the level set method or, when more appropriate, the particle level set method. We show how to implement Dirichlet (known velocity), Neumann (known normal stress), and slip velocity boundary conditions as special cases of our interface representation. The method leads to a discrete, symmetric KKT system for velocities, pressures, and Lagrange multipliers. We also present a novel simplification to the standard combination of the second order semi-Lagrangian and BDF schemes for discretizing the inertial terms. Numerical results indicate second order spatial accuracy for the velocities (L-infinity and L-2) and first order for the pressure (in L-infinity, second order in L-2). Our temporal discretization is also second order accurate. (C) 2014 Elsevier Inc. All rights reserved.
Keywords:
Navier Stokes
Virtual node algorithms
Interface problems
Cartesian grids

Journal

Journal of Computational Physics cover
Journal of Computational Physics
IF:
3.8
Papers:
1.5W
Citations:
7.4W

Organization

University of California System cover
University of California System
Scholars:
37.5W
Papers: 33.7W
Citations: 6.6K