A velocity decomposition approach for moving interfaces in viscous fluids

JOURNAL OF COMPUTATIONAL PHYSICS

卷:228

Article

Beale, J. Thomas¹  Layton, Anita T.¹ 

[1] Duke Univ, Dept Math, Durham, NC 27708 USA

关键词: Navier-Stokes flow;    Stokes flow;    Boundary integral;    Stiff equations;    Fractional stepping;    Immersed interface;    Immersed boundary;   

DOI  :  10.1016/j.jcp.2009.01.023

来源: Elsevier

【 摘 要 】

We present a second-order accurate method for computing the coupled motion of a viscous fluid and an elastic material interface with zero thickness. The fluid flow is described by the Navier-Stokes equations, with a singular force due to the stretching of the moving interface. We decompose the velocity into a Stokes part and a regular part. The first part is determined by the Stokes equations and the singular interfacial force. The Stokes solution is obtained using the immersed interface method, which gives second-order accurate values by incorporating known jumps for the solution and its derivatives into a finite difference method. The regular part of the velocity is given by the Navier-Stokes equations with a body force resulting from the Stokes part. The regular velocity is obtained using a time-stepping method that combines the semi-Lagrangian method with the backward difference formula. Because the body force is continuous, jump conditions are not necessary. For problems with stiff boundary forces, the decomposition approach can be combined with fractional time-stepping, using a smaller time step to advance the interface quickly by Stokes flow, with the velocity computed using boundary integrals. The small time steps maintain numerical stability, while the overall solution is updated on a larger time step to reduce computational cost. (C) 2009 Elsevier Inc. All rights reserved.

【 授权许可】

Free