The Fictitious Domain Method Based on Navier Slip Boundary Condition for Simulation of Flow-Particle Interaction

TL;DR

This paper introduces a least-squares/fictitious domain method for simulating fluid-particle interactions with Navier slip boundary conditions, combining explicit particle updates with implicit fluid and velocity computations.

Contribution

It develops a novel fictitious domain approach with a least-squares formulation for fluid-particle interaction under Navier slip conditions, including an explicit-implicit scheme for stability.

Findings

Method successfully simulates flow-particle interaction.

Numerical results verify accuracy and stability.

Approach handles Navier slip boundary conditions effectively.

Abstract

In this article, we develop a least--squares/fictitious domain method for direct simulation of fluid particle motion with Navier slip boundary condition at the fluid--particle interface. Let $\Omega$ and $B$ be two bounded domains of $\mathbb{R}^{d}$ such that $\overline{B} \subset \Omega$. The motion of solid particle $B$ is governed by Newton's equations. Our goal here is to develop a fictitious domain method where one solves a variant of the original problem on the full $\Omega$, followed by a well--chosen correction over $B$ and corrections related to translation velocity and angular velocity of the particle. This method is of the virtual control type and relies on a least--squares formulation making the problem solvable by a conjugate gradient algorithm operating in a well chosen control space. Since the fully explicit scheme to update the particle motion using Newton's equation is unstable, we propose and implement an explicit--implicit scheme in which, at each time step, the position of the particle is updated explicitly, and the solution of Navier-Stokes equations and particle velocities are solved by the the least--squares/fictitious domain method implicitly. Numerical results are given to verify our numerical method.