A hybrid particle-continuum method for hydrodynamics of complex fluids

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 (\u3c 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

Multiscale modeling of rapid granular flow with a hybrid discrete-continuum method

Both discrete and continuum models have been widely used to study rapid granular flow, discrete model is accurate but computationally expensive, whereas continuum model is computationally efficient but its accuracy is doubtful in many situations. Here we propose a hybrid discrete-continuum method to profit from the merits but discard the drawbacks of both discrete and continuum models. Continuum model is used in the regions where it is valid and discrete model is used in the regions where continuum description fails, they are coupled via dynamical exchange of parameters in the overlap regions. Simulation of granular channel flow demonstrates that the proposed hybrid discrete-continuum method is nearly as accurate as discrete model, with much less computational cost

Multiscale lattice Boltzmann approach to modeling gas flows

For multiscale gas flows, kinetic-continuum hybrid method is usually used to balance the computational accuracy and efficiency. However, the kinetic-continuum coupling is not straightforward since the coupled methods are based on different theoretical frameworks. In particular, it is not easy to recover the non-equilibrium information required by the kinetic method which is lost by the continuum model at the coupling interface. Therefore, we present a multiscale lattice Boltzmann (LB) method which deploys high-order LB models in highly rarefied flow regions and low-order ones in less rarefied regions. Since this multiscale approach is based on the same theoretical framework, the coupling precess becomes simple. The non-equilibrium information will not be lost at the interface as low-order LB models can also retain this information. The simulation results confirm that the present method can achieve model accuracy with reduced computational cost

Multiscale simulation strategies and mesoscale modelling of gas and liquid flows

This paper was presented at the 2nd Micro and Nano Flows Conference (MNF2009), which was held at Brunel University, West London, UK. The conference was organised by Brunel University and supported by the Institution of Mechanical Engineers, IPEM, the Italian Union of Thermofluid dynamics, the Process Intensification Network, HEXAG - the Heat Exchange Action Group and the Institute of Mathematics and its Applications.This paper presents a review of multiscale simulation strategies for the modelling of micro- and nanoscale flows. These have been developed in the last two decades in an attempt to bridge the application gap between molecular and continuum simulation methods preventing the simulation of many micro- and nanofluidic devices. The paper is focused on hybrid molecular-continuum methods and reviews different coupling strategies, including geometrical decomposition in conjunction with state- and flux coupling, pointwise coupling, the heterogeneous multiscale method and the equation free approach. The different applications of these methods are briefly discussed

Study on SPH Viscosity Term Formulations

For viscosity-dominated flows, the viscous effect plays a much more important role. Since the viscosity term in SPH-governing (Smoothed Particle Hydrodynamics) equations involves the discretization of a second-order derivative, its treatment could be much more challenging than that of a first-order derivative, such as the pressure gradient. The present paper summarizes a series of improved methods for modeling the second-order viscosity force term. By using a benchmark patch test, the numerical accuracy and efficiency of different approaches are evaluated under both uniform and non-uniform particle configurations. Then these viscosity force models are used to compute a documented lid-driven cavity flow and its interaction with a cylinder, from which the most recommended viscosity term formulation has been identified