high performance computing on graphics processing units: hgpu.org

We present a generalization of the inertial coupling (IC) [Usabiaga et al. J. Comp. Phys. 2013] which permits the resolution of radiation forces on small particles with arbitrary acoustic contrast factor. The IC method is based on a Eulerian-Lagrangian approach: particles move in continuum space while the fluid equations are solved in a regular mesh (here we use the finite volume method). Thermal fluctuations in the fluid stress, important below the micron scale, are also taken into account following the Landau-Lifshitz fluid description. Each particle is described by a minimal cost resolution which consists on a single small kernel (bell-shaped function) concomitant to the particle. The main role of the particle kernel is to interpolate fluid properties and spread particle forces. Here, we extend the kernel functionality to allow for an arbitrary particle compressibility. The particle-fluid force is obtained from an imposed no-slip constraint which enforces similar particle and kernel fluid velocities. This coupling is instantaneous and permits to capture the fast, non-linear effects underlying the radiation forces on particles. Acoustic forces arise either because an excess in particle compressibility (monopolar term) or in mass (dipolar contribution) over the fluid values. Comparison with theoretical expressions show that the present generalization of the IC method correctly reproduces both contributions. Due to its low computational cost, the present method allows for simulations with many particles using a standard Graphical Processor Unit (GPU).