Variational solution to the lattice Boltzmann method for Couette flow.

Literature studies of the lattice Boltzmann method (LBM) demonstrate hydrodynamics beyond the continuum limit. This includes exact analytical solutions to the LBM, for the bulk velocity and shear stress of Couette flow under diffuse reflection at the walls through the solution of equivalent moment equations. We prove that the bulk velocity and shear stress of Couette flow with Maxwell-type boundary conditions at the walls, as specified by two-dimensional isothermal lattice Boltzmann models, are inherently linear in Mach number. Our finding enables a systematic variational approach to be formulated that exhibits superior computational efficiency than the previously reported moment method. Specifically, the number of partial differential equations (PDEs) in the variational method grows linearly with quadrature order while the number of moment method PDEs grows quadratically. The variational method directly yields a system of linear PDEs that provide exact analytical solutions to the LBM bulk velocity field and shear stress for Couette flow with Maxwell-type boundary conditions. It is anticipated that this variational approach will find utility in calculating analytical solutions for novel lattice Boltzmann quadrature schemes and other flows.

Highly Spherical Nanoparticles Probe Gigahertz Viscoelastic Flows of Simple Liquids Without the No-Slip Condition.

Simple liquids are conventionally described by Newtonian fluid mechanics, based on the assumption that relaxation processes in the flow occur much faster than the rate at which the fluid is driven. Nanoscale solids, however, have characteristic mechanical response times on the picosecond scale, which are comparable to mechanical relaxation times in simple liquids; as a result, viscoelastic effects in the liquid must be considered. These effects have been observed using time-resolved optical measurements of vibrating nanoparticles, but interpretation has often been complicated by finite velocity slip at the liquid-solid interface. Here, we use highly spherical gold nanoparticles to drive flows that are theoretically modeled without the use of the no-slip boundary condition at the particle surface. We obtain excellent agreement with this analytical theory that considers both the compression and shear relaxation properties of the liquid.

Influence of polydispersity on micromechanics of granular materials.

We study the effect of polydispersity on the macroscopic physical properties of granular packings in two and three dimensions. A mean-field approach is developed to approximate the macroscale quantities as functions of the microscopic ones. We show that the trace of the fabric and stress tensors are proportional to the mean packing properties (e.g., packing fraction, average coordination number, and average normal force) and dimensionless correction factors, which depend only on the moments of the particle-size distribution. Similar results are obtained for the elements of the stiffness tensor of isotropic packings in the linear affine response regime. Our theoretical predictions are in good agreement with the simulation results.

Morphology and linear-elastic moduli of random network solids.

The effective linear-elastic moduli of disordered network solids are analyzed by voxel-based finite element calculations. We analyze network solids given by Poisson-Voronoi processes and by the structure of collagen fiber networks imaged by confocal microscopy. The solid volume fraction ϕ is varied by adjusting the fiber radius, while keeping the structural mesh or pore size of the underlying network fixed. For intermediate ϕ, the bulk and shear modulus are approximated by empirical power-laws K(phi)proptophin and G(phi)proptophim with n≈1.4 and m≈1.7. The exponents for the collagen and the Poisson-Voronoi network solids are similar, and are close to the values n=1.22 and m=2.11 found in a previous voxel-based finite element study of Poisson-Voronoi systems with different boundary conditions. However, the exponents of these empirical power-laws are at odds with the analytic values of n=1 and m=2, valid for low-density cellular structures in the limit of thin beams. We propose a functional form for K(ϕ) that models the cross-over from a power-law at low densities to a porous solid at high densities; a fit of the data to this functional form yields the asymptotic exponent n≈1.00, as expected. Further, both the intensity of the Poisson-Voronoi process and the collagen concentration in the samples, both of which alter the typical pore or mesh size, affect the effective moduli only by the resulting change of the solid volume fraction. These findings suggest that a network solid with the structure of the collagen networks can be modeled in quantitative agreement by a Poisson-Voronoi process.

Lattice Boltzmann simulation of fluid flow in fracture networks with rough, self-affine surfaces.

Using the lattice Boltzmann method, we study fluid flow in a two-dimensional (2D) model of fracture network of rock. Each fracture in a square network is represented by a 2D channel with rough, self-affine internal surfaces. Various parameters of the model, such as the connectivity and the apertures of the fractures, the roughness profile of their surface, as well as the Reynolds number for flow of the fluid, are systematically varied in order to assess their effect on the effective permeability of the fracture network. The distribution of the fractures' apertures is approximated well by a log-normal distribution, which is consistent with experimental data. Due to the roughness of the fractures' surfaces, and the finite size of the networks that can be used in the simulations, the fracture network is anisotropic. The anisotropy increases as the connectivity of the network decreases and approaches the percolation threshold. The effective permeability K of the network follows the power law K approximately (beta), where is the average aperture of the fractures in the network and the exponent beta may depend on the roughness exponent. A crossover from linear to nonlinear flow regime is obtained at a Reynolds number Re approximately O(1), but the precise numerical value of the crossover Re depends on the roughness of the fractures' surfaces.