A smoothed particle hydrodynamics-based fluid model with a spatially dependent viscosity: application to flow of a suspension with a non-Newtonian fluid matrix

A smoothed particle hydrodynamics approach is utilized to model a non-Newtonian fluid with a spatially varying viscosity. In the limit of constant viscosity, this approach recovers an earlier model for Newtonian fluids of Español and Revenga (Phys Rev E 67:026705, 2003). Results are compared with numerical solutions of the general Navier–Strokes equation using the “regularized” Bingham model of Papanastasiou (J Rheol 31:385–404, 1987) that has a shear-rate-dependent viscosity. As an application of this model, the effect of having a non-Newtonian fluid matrix, with a shear-rate-dependent viscosity in a moderately dense suspension, is examined. Simulation results are then compared with experiments on mono-size silica spheres in a shear-thinning fluid and for sand in a calcium carbonate paste. Excellent agreement is found between simulation and experiment. These results indicate that measurements of the shear viscosity of simple shear-rate-dependent non-Newtonian fluids may be used in simulation to predict the viscosity of concentrated suspensions having the same matrix fluid.     

A numerical study of the flow of Bingham-like fluids in two-dimensional vane and cylinder rheometers using a smoothed particle hydrodynamics (SPH) based method

In this paper, a Lagrangian formulation of the Navier–Stokes equations, based on the smoothed particle hydrodynamics (SPH) approach, was applied to determine how well rheological parameters such as plastic viscosity can be determined from vane rheometer measurements. First, to validate this approach, a Bingham/Papanastasiou constitutive model was implemented into the SPH model and tests comparing simulation results to well established theoretical predictions were conducted. Numerical simulations for the flow of fluids in vane and coaxial cylinder rheometers were then performed. A comparison to experimental data was also made to verify the application of the SPH method in realistic flow geometries. Finally, results are presented from a parametric study of the flow of Bingham fluids with different yield stresses under various applied angular velocities of the outer cylindrical wall in the vane and coaxial cylinder rheometers. The stress, strain rate and velocity profiles, especially in the vicinity of the vane blades, were computed. By comparing the calculated stress and flow fields between the two rheometers, the validity of the assumption that the vane could be approximated as a cylinder for measuring the rheological properties of Bingham fluids at different shear rates was tested.     

A comparison of viscosity-concentration relationships for emulsions

Differential effective medium theory (D-EMT) has been used by a number of investigators to derive expressions for the shear viscosity of a colloidal suspension or an emulsion as a function of the volume fraction of the dispersed phase. Pal and Rhodes [R. Pal, E. Rhodes, J. Rheol. 33 (7) (1989) 1021-1045] used D-EMT to derive a viscosity-concentration expression for non-Newtonian emulsions, in which variations among different oil-water emulsions were accommodated by fitting the value of an empirical solvation factor by matching the volume fraction at which the ratio of each emulsion was experimentally observed to have a viscosity 100 times greater than that of the pure solvent. When the particles in suspension have occluded volume due to solvation or flocculation, we show that the application of D-EMT to the problem becomes more ambiguous than these investigators have indicated. In addition, the resulting equations either do not account for the limiting behavior near the critical concentration, that is, the concentration at which the viscosity diverges, or they incorporate this critical behavior in an ad hoc way. We suggest an alternative viscosity-concentration equation for emulsions, based on work by Bicerano and coworkers [J. Bicerano, J.F. Douglas, D.A. Brune, J. Macromol. Sci., Rev. Macromol. Chem. Phys. C 39 (4) (1999) 561-642]. This alternative equation has the advantages that (1) its parameters are more closely related to physical properties of the suspension and (2) it recovers the correct limiting behavior both in the dilute limit and near the critical concentration for rigid particles. In addition, the equation can account for the deformability of flexible particles in the semidilute regime. The proposed equation is compared to the equation proposed by Pal and Rhodes.     

Hydraulic radius and transport in reconstructed model three-dimensional porous media

Methods for reconstructing three-dimensional porous media from two-dimensional cross sections are evaluated in terms of the transport properties of the reconstructed systems. Two-dimensional slices are selected at random from model three-dimensional microstructures, based on penetrable spheres, and processed to create a reconstructed representation of the original system. Permeability, conductivity, and a critial pore diameter are computed for the original and reconstructed microstructures to assess the validity of the reconstruction technique. A surface curvature algorithm is utilized to further modify the reconstructed systems by matching the hydraulic radius of the reconstructed three-dimensional system to that of the two-dimensional slice. While having only minor effects on conductivity, this modification significantly improves the agreement between permeabilities and critical diameters of the original and reconstructed systems for porosities in the range of 25–40%. For lower porosities, critical pore diameter is unaffected by the curvature modification so that little improvement between original and reconstructed permeabilities is obtained by matching hydraulic radii.