Fluctuating hydrodynamics in periodic domains and heterogeneous adjacent multidomains: Thermal equilibrium

Xin Bian, Zhen Li, Mingge Deng, and George Em Karniadakis
Phys. Rev. E 92, 053302 – Published 10 November 2015

Abstract

We first study fluctuating hydrodynamics (FH) at equilibrium in periodic domains by use of the smoothed dissipative particle dynamics (SDPD) method. We examine the performance of SDPD by comparing it with the theory of FH. We find that the spatial correlation of particle velocity is always the Dirac δ function, irrespective of numerical resolution, in agreement with the theory. However, the spatial correlation of particle density has a finite range of rc, which is due to the kernel smoothing procedure for the density. Nevertheless, this finite range of correlation can be reduced to an arbitrarily small value by increasing the resolution, that is, reducing rc, similarly to how the smoothing kernel converges to the Dirac δ function. Moreover, we consider temporal correlation functions (CFs) of random field variables in Fourier space. For sufficient resolution, the CFs of SDPD simulations agree very well with analytical solutions of the linearized FH equations. This confirms that both the shear and sound modes are modeled accurately and that fluctuations are generated, transported, and dissipated in both thermodynamically and hydrodynamically consistent ways in SDPD. We also show that the CFs of the classical dissipative particle dynamics (DPD) method with proper parameters can recover very well the linearized solutions. As a reverse implication, the measurement of CFs provides an effective means of extracting viscosities and sound speed of a DPD system with a new set of input parameters. Subsequently, we study the FH in truncated domains in the context of multiscale coupling via the domain decomposition method, where a SDPD simulation in one subdomain is coupled with a Navier-Stokes (NS) solver in an adjacent subdomain with an overlapping region. At equilibrium, the mean values of the NS solution are known a priori and do not need to be extracted from actual simulations. To this end, we model a buffer region as an equilibrium boundary condition (EBC) at the truncated side of the SDPD simulation. In the EBC buffer, the velocity of particles is drawn from a known Gaussian distribution, that is, the Maxwell-Boltzmann distribution. Due to the finite range of spatial correlation, the density of particles in the EBC buffer must be drawn from a conditional Gaussian distribution, which takes into account the available density distribution of neighboring interior particles. We introduce a Kriging method to provide such a conditional distribution and hence preserve the spatial correlation of density. Spatial and temporal correlations of SDPD simulations in the truncated domain are compared to that in a single complete domain. We find that a gap region between the buffer and interior is important to reduce the extra dissipation generated by the artificial buffer at equilibrium, rendering more investigations necessary for thermal fluctuations in the multiscale coupling of nonequilibrium flows.

  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
21 More
  • Received 14 May 2015
  • Revised 8 September 2015

DOI:https://doi.org/10.1103/PhysRevE.92.053302

©2015 American Physical Society

Authors & Affiliations

Xin Bian*, Zhen Li, Mingge Deng, and George Em Karniadakis

  • Division of Applied Mathematics, Brown University, Providence, Rhode Island 02912, USA

  • *xin_bian@brown.edu
  • george_karniadakis@brown.edu

Article Text (Subscription Required)

Click to Expand

References (Subscription Required)

Click to Expand
Issue

Vol. 92, Iss. 5 — November 2015

Reuse & Permissions
Access Options
CHORUS

Article Available via CHORUS

Download Accepted Manuscript
Author publication services for translation and copyediting assistance advertisement

Authorization Required


×
×

Images

×

Sign up to receive regular email alerts from Physical Review E

Log In

Cancel
×

Search


Article Lookup

Paste a citation or DOI

Enter a citation
×