Abstract:The Direction Splitting algorithm is a computational method that facilitates the solution of the Navier–Stokes equation by decomposing a multi-dimensional Poisson equation into a series of one-dimensional equations. This approach allows for efficient and scalable numerical computations. The incorporation of rigid bodies in the computational domain is made possible by modifying the diffusion and divergence stencils near the fluid-solid interface. In our recent study (Goyal and Wachs, 2023), we implemented stencil corrections for fixed non-spherical rigid bodies and demonstrated the scalability and accuracy of the DS algorithm for a wide range of test cases. This study presents an extension of the Direction Splitting solver to effectively compute the dynamics of moving non-spherical rigid bodies in the presence of fluid flow. The interaction between rigid bodies and fluid can be categorized as either one-way or two-way coupled. The stencils resolving the fluid-sloid interface are updated at every time-iteration to account for changes in particle position and orientation. We demonstrate that the adaptation of the diffusion and divergence stencil in Goyal and Wachs (2023); Keating and Minev (2013) is inadequate for accurately predicting the dynamics of moving rigid bodies in a two-way coupled simulation, especially at O(102) Reynolds number. Hence, we implement a conservative cut-cell method for the discretization of the divergence term, and we also take into account the proximity of the fluid-solid interface in the advection discretization. The results of multiple test cases demonstrate that the suggested stencil modifications effectively capture the dynamics of fluid flow in the presence of freely moving rigid bodies. However, the cut-cell discretization of the divergence term constitutes a significant computing overhead that impacts the fast computing property of the Direction Splitting algorithm.

Dimension Splitting Algorithms for Three Dimensional Steady/Unsteady Stokes Equations

A Dimension Split Method for the Incompressible Navier–Stokes Equations in Three Dimensions

Dimension splitting method for the three dimensional rotating Navier-Stokes equations

A Dimension Splitting Method for 3-D Incompressible Thermal Flow

A dimension splitting and characteristic projection method for three-dimensional incompressible flow

Parallel Finite Element Algorithms Based on Fully Overlapping Domain Decomposition for Time-dependent Navier-Stokes Equations

Block iterative methods for linear algebraic equation and domain decomposition method for incompressible viscous flow

Dimension Splitting Method for 3d Rotating Compressible Navier-Stokes Equations in the Turbomachinery

A Parallel Partition of Unity Method for the Unsteady Stokes Equations

Fast and stable dimension splitting simulations for the hydrodynamically coupled three-component conserved Allen-Cahn phase field model

A Parallel Fast Multipole Accelerated Integral Equation Scheme for 3d Stokes Equations

Local and Parallel Finite Element Algorithms Based on the Partition of Unity for the Stokes Problem

A PARALLEL SOLVING DESIGN FOR UNSTEADY NAVIER-STOKES EQUATION

A Dimension Split Method for the 3‐D Compressible Navier–Stokes Equations in Turbomachine

Parallel Finite Element Algorithm Based on Full Domain Partition for Stationary Stokes Equations

Efficient Numerical Simulation of Cahn-Hilliard Type Models by a Dimension Splitting Method

Parallel Domain Decomposition method for the fully-mixed Stokes-dual-permeability fluid flow model with Beavers-Joseph interface conditions

A Two-Level Consistent Splitting Scheme for the Navier-Stokes Equations

The Navier-Stokes Equations in Stream Layer or on Stream Surface and a Dimension Split Method

Decoupled Modified Characteristic Finite Element Method with Different Subdomain Time Steps for Nonstationary Dual-Porosity-Navier-Stokes Model

An accurate and scalable direction-splitting solver for flows laden with non-spherical rigid bodies—Part 2: Moving rigid bodies