期刊
MULTISCALE MODELING & SIMULATION
卷 8, 期 3, 页码 871-911出版社
SIAM PUBLICATIONS
DOI: 10.1137/090774501
关键词
fluctuating hydrodynamics; hybrid methods; direct simulation Monte Carlo; adiabatic piston
A previously developed hybrid particle-continuum method [J. B. Bell, A. Garcia, and S. A. Williams, Multiscale Model. Simul., 6 (2008), pp. 1256-1280] is generalized to dense fluids and two-and three-dimensional flows. The scheme couples an explicit fluctuating compressible Navier Stokes solver with the isotropic direct simulation Monte Carlo (DSMC) particle method [A. Donev, A. L. Garcia, and B. J. Alder, J. Stat. Mech. Theory Exp., 2009 (2009), article P11008]. To achieve bidirectional dynamic coupling between the particle (microscale) and continuum (macroscale) regions, the continuum solver provides state-based boundary conditions to the particle subdomain, while the particle solver provides flux-based boundary conditions for the continuum subdomain. This type of coupling ensures both state and flux continuity across the particle-continuum interface analogous to coupling approaches for deterministic parabolic partial differential equations; here, when fluctuations are included, a small (< 1%) mismatch is expected and observed in the mean density and temperature across the interface. By calculating the dynamic structure factor for both a bulk (periodic) and a finite system, it is verified that the hybrid algorithm accurately captures the propagation of spontaneous thermal fluctuations across the particle-continuum interface. The equilibrium diffusive (Brownian) motion of a large spherical bead suspended in a particle fluid is examined, demonstrating that the hybrid method correctly reproduces the velocity autocorrelation function of the bead but only if thermal fluctuations are included in the continuum solver. Finally, the hybrid is applied to the well-known adiabatic piston problem, and it is found that the hybrid correctly reproduces the slow nonequilibrium relaxation of the piston toward thermodynamic equilibrium but, again, only if the continuum solver includes stochastic (white-noise) flux terms. These examples clearly demonstrate the need to include fluctuations in continuum solvers employed in hybrid multiscale methods.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据