The Numerical Scheme for the Basset Type Integro-Differential Equation in Hydrodynamics

A high-order precision numerical scheme was proposed for the Basset kernel integrals and the Basset type integro-differential equations that often arise in hydrodynamics. The scheme was tested on a model integral and on a real problem of particle focusing in a standing acoustic wave in liquid. The scheme showed around 2 orders higher speed of integral estimation compared with the existing analogs. The obtained results show that the integral can be approximated with high order of precision and the real problem is simulated well by the proposed scheme. A variable step technique was used to increase the precision of integro-differential equation simulation even more by eliminating the discrepancy at the very start of the simulation. The proposed numerical scheme may found its applications in many biological and medical problems related to acoustophoresis or in particle sedimentation simulations.