TY - JOUR
T1 - Inertial coupling method for particles in an incompressible fluctuating fluid
AU - Balboa Usabiaga, Florencio
AU - Delgado-Buscalioni, Rafael
AU - Griffith, Boyce E.
AU - Donev, Aleksandar
N1 - Funding Information:
We thank Alejandro Garcia, Charles Peskin, Paul Atzberger, Eric Vanden-Eijnden, Martin Maxey, Tony Ladd and Burkhard Dünweg for informative and inspiring discussions. A. Donev was supported in part by the Air Force Office of Scientific Research under Grant No. FA9550-12-1-0356. B. Griffith acknowledges research support from the National Science Foundation under awards OCI 1047734 and DMS 1016554. R. Delgado-Buscalioni and F. Balboa acknowledge funding from the Spanish government FIS2010-22047-C05 and from the Comunidad de Madrid MODELICO-CM (S2009/ESP-1691). Collaboration between A. Donev and R. Delgado-Buscalioni was fostered at the Kavli Institute for Theoretical Physics in Santa Barbara, California, and supported in part by the National Science Foundation under Grant No. NSF PHY05-51164.
PY - 2014/2/1
Y1 - 2014/2/1
N2 - We develop an inertial coupling method for modeling the dynamics of point-like "blob" particles immersed in an incompressible fluid, generalizing previous work for compressible fluids (Balboa Usabiaga et al., 2013 [42]). The coupling consistently includes excess (positive or negative) inertia of the particles relative to the displaced fluid, and accounts for thermal fluctuations in the fluid momentum equation. The coupling between the fluid and the blob is based on a no-slip constraint equating the particle velocity with the local average of the fluid velocity, and conserves momentum and energy. We demonstrate that the formulation obeys a fluctuation-dissipation balance, owing to the non-dissipative nature of the no-slip coupling. We develop a spatio-temporal discretization that preserves, as best as possible, these properties of the continuum formulation. In the spatial discretization, the local averaging and spreading operations are accomplished using compact kernels commonly used in immersed boundary methods. We find that the special properties of these kernels allow the blob to provide an effective model of a particle; specifically, the volume, mass, and hydrodynamic properties of the blob are remarkably grid-independent. We develop a second-order semi-implicit temporal integrator that maintains discrete fluctuation-dissipation balance, and is not limited in stability by viscosity. Furthermore, the temporal scheme requires only constant-coefficient Poisson and Helmholtz linear solvers, enabling a very efficient and simple FFT-based implementation on GPUs. We numerically investigate the performance of the method on several standard test problems. In the deterministic setting, we find the blob to be a remarkably robust approximation to a rigid sphere, at both low and high Reynolds numbers. In the stochastic setting, we study in detail the short and long-time behavior of the velocity autocorrelation function and observe agreement with all of the known behavior for rigid sphere immersed in a fluctuating fluid. The proposed inertial coupling method provides a low-cost coarse-grained (minimal resolution) model of particulate flows over a wide range of time-scales ranging from Brownian to convection-driven motion.
AB - We develop an inertial coupling method for modeling the dynamics of point-like "blob" particles immersed in an incompressible fluid, generalizing previous work for compressible fluids (Balboa Usabiaga et al., 2013 [42]). The coupling consistently includes excess (positive or negative) inertia of the particles relative to the displaced fluid, and accounts for thermal fluctuations in the fluid momentum equation. The coupling between the fluid and the blob is based on a no-slip constraint equating the particle velocity with the local average of the fluid velocity, and conserves momentum and energy. We demonstrate that the formulation obeys a fluctuation-dissipation balance, owing to the non-dissipative nature of the no-slip coupling. We develop a spatio-temporal discretization that preserves, as best as possible, these properties of the continuum formulation. In the spatial discretization, the local averaging and spreading operations are accomplished using compact kernels commonly used in immersed boundary methods. We find that the special properties of these kernels allow the blob to provide an effective model of a particle; specifically, the volume, mass, and hydrodynamic properties of the blob are remarkably grid-independent. We develop a second-order semi-implicit temporal integrator that maintains discrete fluctuation-dissipation balance, and is not limited in stability by viscosity. Furthermore, the temporal scheme requires only constant-coefficient Poisson and Helmholtz linear solvers, enabling a very efficient and simple FFT-based implementation on GPUs. We numerically investigate the performance of the method on several standard test problems. In the deterministic setting, we find the blob to be a remarkably robust approximation to a rigid sphere, at both low and high Reynolds numbers. In the stochastic setting, we study in detail the short and long-time behavior of the velocity autocorrelation function and observe agreement with all of the known behavior for rigid sphere immersed in a fluctuating fluid. The proposed inertial coupling method provides a low-cost coarse-grained (minimal resolution) model of particulate flows over a wide range of time-scales ranging from Brownian to convection-driven motion.
KW - Brownian motion
KW - Fluctuating hydrodynamics
KW - Immersed-boundary method
KW - Inertial coupling
KW - Minimally-resolved particulate flows
UR - http://www.scopus.com/inward/record.url?scp=84888260181&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84888260181&partnerID=8YFLogxK
U2 - 10.1016/j.cma.2013.10.029
DO - 10.1016/j.cma.2013.10.029
M3 - Article
AN - SCOPUS:84888260181
SN - 0045-7825
VL - 269
SP - 139
EP - 172
JO - Computer Methods in Applied Mechanics and Engineering
JF - Computer Methods in Applied Mechanics and Engineering
ER -