Efficient Numerical Simulation of Biochemotaxis Phenomena in Fluid Environments.

Area of Science:

• Computational fluid dynamics
• Biochemical modeling
• Numerical analysis

Background:

• Biochemotaxis models describe cell movement in response to chemical gradients.
• Efficient numerical simulations are crucial for understanding complex biological processes.
• Coupled chemotaxis-fluid and Navier-Stokes equations present significant computational challenges.

Purpose of the Study:

• To propose a novel dimension splitting method for efficient numerical simulation of a biochemotaxis model.
• To reduce computational and storage costs associated with simulating coupled chemotaxis-fluid and Navier-Stokes equations.
• To analyze the evolution of cell plume shape and the impact of parameter changes on cell descent.

Main Methods:

• A second-order pressure correction method decouples velocity and pressure for Navier-Stokes equations.
• An alternating direction implicit scheme solves the velocity equation.
• Dimension splitting is applied to the pressure equation, replacing the traditional elliptic operator.
• Operator splitting and extrapolation techniques achieve second-order time accuracy for chemotactic equations.

Main Results:

• The proposed dimension splitting method effectively reduces computational and storage requirements.
• The method enables efficient simulation of the coupled biochemotaxis model.
• Experiments demonstrate the evolution of cell plume shape during descent.
• The study investigates the influence of parameter variations on velocity and plume dynamics.

Conclusions:

• The novel dimension splitting method offers an efficient approach for simulating biochemotaxis.
• This method significantly lowers computational and storage demands.
• The findings provide insights into cell plume dynamics and parameter effects during descent.