The flow of fiber suspensions in complex geometries

A continuum theory for dilute suspensions of large-aspect-ratio particles is applied to the flow of fiber suspensions through contractions. The theory, which incorporates the statistical orientation distribution function into the stress equation, predicts that the flow of dilute suspensions will differ qualitatively from the flow of the suspending fluid. The theory is in excellent agreement with experiments on the flow of suspensions of chopped-glass fibers through axisymmetric contractions, where substantial enlargement of the recirculating corner vortex is observed at volume fractions of 0.1% and less.     

What do dry granular flows tell us about dense non-Brownian suspension rheology?

Recent years have seen significant progress in our understanding of the rheology of dry granular materials. A scale invariance of equations of motion in the rigid grain limit has helped identifying dimensionless quantities which govern flow. After a review of recent results on dry granular materials, we show how this same dimensional analysis carries over to the case of dense, non-Brownian suspensions. Our review is based on compiled data from various sets of numerical simulations, using both molecular and contact dynamics. It covers the breakdown of kinetic theory, which arises when contact forces dominate collisional forces in inertial flows, and the approach to dense flows up to the quasi-static limit. We show that simple invariance arguments permit to clarify the conditions of occurrence of viscous and inertial scaling in suspensions, and show in particular that both may occur in dense flows even in the presence of significant contact forces, up to the jamming limit. Some implications of the properties of steady uniform quasistatic granular flow under constant pressure for very dense suspensions near the maximum concentration are also discussed.     

Behavior of a net of fibers linked by viscous interactions: theory and mechanical properties

This paper presents an investigation of the macroscopic mechanical behavior of highly concentrated fiber suspensions for which the mechanical behavior is governed by local fiber-fiber interactions.The problem is approached by considering the case of a net of rigid fibers of uniform length, linked by viscous point interactions of power-law type. Those interactions may result in local forces and moments located at the contacting point between two fibers, and respectively power-law functions of the local linear and angular velocity at this point.Assuming the existence of an elementary representative volume which size is small compared to the size of the whole structure, the fiber net is regarded as a periodic assembly of identical cells. Macroscopic equilibrium and constitutive equations of the equivalent continuum are then obtained by the discrete and periodic media homogenization method, based on the use of asymptotic expansions.Depending on the order of magnitude of local translational viscosities and rotational viscosities, three types of the equivalent continua are proved to be possible. One of them leads to an effective Cosserat medium, the other ones being usual Cauchy media. Lastly, formulations that enable an effective computation of constitutive equations are detailed. They show that the equivalent continuum behaves like an anisotropic power-law fluid.     

Stable and contact-free time stepping for dense rigid particle suspensions

We consider suspensions of rigid bodies in a two-dimensional viscous fluid. Even with high-fidelity numerical methods, unphysical contact between particles occurs because of spatial and temporal discretization errors. We extend a time stepping method that avoids overlap by imposing a minimum separation distance between all pairs of bodies. In its original form, the method discretizes interactions between different particles explicitly. Therefore, to avoid stiffness, a large minimum separation distance is used. In this paper, we introduce a new implicit time stepping method that is able to simulate dense suspensions with large time step sizes and a small minimum separation distance. The method is tested on various unbounded and bounded flows, and rheological properties of the resulting suspensions are computed.     

Flow-driven orientation dynamics of semiflexible fiber systems

We study the flow-induced orientation dynamics of semiflexible fibers in dilute fiber suspensions. Starting from the equations of motion for a two-rod model of flexible fibers in Stokes flow, the Smoluchowski equation for a connected monomer orientation distribution function is derived. We then obtain a set of equations for the time dependence of the first and second moments of the orientation distribution function, thus extending the Folgar Tucker equations for short rigid fiber suspensions to flexible fiber suspensions. The resulting generalized equations for the orientation dynamics of a suspension of flexible fibers are solved for simple channel flow. It is shown that all qualitative effects of bending and straightening of fibers and their influence on the orientation of flexible fibers are captured within our model. A scalar measure for the distribution of bending in a flow is introduced, which allows to detect the degree of bending of fibers. Paper was presented at the 3rd Annual Rheology Conference, AERC 2006, April 27–29, 2006, Crete, Greece.     

理论力学研究性教学新探索: 刚体运动基点法公式与连续介质速度场分解

基础力学课程体系是现代工程教育的重要基础。理论力学作为力学课程体系的起点,提供了力学基本概念和原理。一方面,以离散体系(质点/刚体)为研究对象的理论力学与以连续体系(介质/变形微元体)为对象的连续介质力学,它们之间既存在密切联系又差异明显;另一方面,源于理论/经典力学的动力系统理论和非线性科学研究对应用力学学科产生了广泛影响。围绕这两个侧面,作者将在《理论力学研究性教学新探索》系列教研论文中阐释它们之间的联系。本篇探讨刚体运动基点法公式和连续介质速度场分解之间的关系,它们分别给出刚体模型和连续介质模型的速度分布规律,前者依赖刚体转动角速度矢量,而后者由速度梯度张量所刻画。本文将说明,两者存在对应理论关系,且刚体基点法公式是连续介质速度场分解的退化形式,即忽略变形效应。

     

高浓度固-液两相流紊流的动理学模型

采用分子动理学方法,基于固-液两相流液相分子或颗粒相颗粒的Boltzmann方程,对Boltzmann方程分别取零矩和一次矩,则得到高浓度固-液两相流紊流的连续方程和动量方程,再和较成熟的低浓度两相流连续方程和动量方程比较,取低浓度两相流控制方程中较成熟合理的有关项和高浓度时由动理学方法推导出的颗粒间碰撞项,则得到高浓度固-液两相流紊流的最终控制方程:连续方程和动量方程.     

Some unsteady parallel flows of particulate suspensions

Analytical and numerical solutions are obtained for five unsteady parallel flows of particulate suspensions which exhibit boundary-layer characteristics. The governing equations are based on a continuum theory of particulate suspension behavior. An attempt is made to infer some of the pertinent features of two-phase and multiphase boundary-layer flows from the solutions.     

A dynamic theory for magnetizable elastic solids with thermal and electrical conduction

A general dynamical theory of magnetizable, electrically and thermally conducting media is developed for soft ferromagnetic or paramagnetic materials in external electromagnetic fields. The general equations are linearized by assuming infinitesimal strains, linear constitutive equations and that all field variables may be divided into two parts: a "rigid body state" and a "perturbation state". The former is the same as the one in rigid body electrodynamics, and the latter which accounts for electromagnetic interaction with the deformable continuum is coupled with stress and strain through linearized field equations. The theory is developed for general anisotropy but specialized for materials with uniaxial, or higher, symmetry.     

The Fast Exact Closure for Jeffery’s equation with diffusion

Jeffery’s equation with diffusion is widely used to predict the motion of concentrated fiber suspensions in flows with low Reynold’s numbers. Unfortunately, the evaluation of the fiber orientation distribution can require excessive computation, which is often avoided by solving the related second order moment tensor equation. This approach requires a ‘closure’ that approximates the distribution function’s fourth order moment tensor from its second order moment tensor. This paper presents the Fast Exact Closure (FEC) which uses conversion tensors to obtain a pair of related ordinary differential equations; avoiding approximations of the higher order moment tensors altogether. The FEC is exact in that when there are no fiber interactions, it exactly solves Jeffery’s equation. Numerical examples for dense fiber suspensions are provided with both a Folgar–Tucker (1984) [3] diffusion term and the recent anisotropic rotary diffusion term proposed by Phelps and Tucker (2009) [9]. Computations demonstrate that the FEC exhibits improved accuracy with computational speeds equivalent to or better than existing closure approximations.     

格栅结构力学性能研究进展

格栅复合材料是一种新型轻质高强材料. 综述了格栅复合材料的周期构型特征和格栅结构的制备工艺. 归纳了二维周期格栅材料的等效刚度矩阵计算方法, 比较了不同构型格栅的基本力学性能, 介绍了胞元材料的微极弹性理论和格栅的强度与屈服面计算方法. 探讨了格栅的缺陷及其力学响应, 包括格栅的尺度效应、夹杂缺陷以及裂纹扩展特征, 介绍了波在格栅材料中传播机理的最新研究成果. 根据格栅材料在工程中的应用形式, 分类介绍了格栅板壳结构、格栅加筋板壳结构和格栅夹层结构的结构特点和破坏方式、设计优化准则和实验研究成果. 还归纳了作者所在研究小组近期在碳纤维格栅复合材料的制备、实验研究和理论分析等方面的最新工作进展.      

Contemporary perspectives in finite strain plasticity

From an internally consistent theory of finite strain plasticity in which are introduced a new stress tensor that includes rigid body rotation and a compatible finite strain tensor amenable to simple geometrical interpretation, a close correlation with experiment is obtained for stress paths of arbitrary composition and direction. Measured material constants are proportional to the elastic shear modulus, as indicated for dislocation theory. At large strain there exists a quantum mechanical structure for the continuum, in which are defined a series of relative reference configurations and second order transitions.     

A numerical study of shock wave reflections on low density foam

A continuum mixture theory is used to describe shock wave reflections on low density open-cell polyurethane foam. Numerical simulations are compared to the shock tube experiments of Skews (1991) and detailed wave fields are shown of a shock wave interacting with a layer of foam adjacent to a rigid wall boundary. These comparisons demonstrate that a continuum mixture theory describes well the shock interactions with low density foam.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

Comparative numerical study of two concentrated fiber suspension models

We consider two rheological models for concentrated fiber suspensions. In both models the equations for orientation and flow are fully coupled, i.e., the orientation influences the flow via a constitutive relation for the viscosity and the orientation of the fibers is determined by the flow field. The orientation state of the fibers is characterized by the Advani–Tucker orientation tensor. We are investigating suspensions of fibers in which the kinetic energies of the fibers are large compared to the thermal energies, i.e., the influence of Brownian motion may be neglected. The first model is the Folgar–Tucker model with backcoupling to the flow (FT model). The second model is an extension of Folgar–Tucker, which models phenomenologically the topological exclusion interaction in dense suspensions (FTMS model). As test cases for the simulation are considered channel flow, 8:1 contraction flow and flow around a cylinder.     

Kinetic theory and rheology of dilute polymer solutions

Summary A summary is given of some recent attempts to relate the results of the kinetic theory of rigid and flexible macromolecules to continuum mechanics results.With 1 table     

Direct simulation of flexible fibers

A numerical simulation of multiple flexible fibers in suspension in Newtonian simple shear flow is presented. The method used is similar to those of previous recent simulation works by Fan et al. [J. Non-Newtonian Fluid Mech. 74 (1998) 113] and Yamane et al. [J. Non-Newtonian Fluid Mech. 54 (1994) 405], however, the method has been modified to allow a small amount of bending and torsion in the fibers. A restoring moment acts to straighten the fibers as they interact in the flow.It is demonstrated that this simulation can be used to extract basic rheological information about the suspension including fiber orientations and suspension viscosity. The viscosity of semi-concentrated to concentrated flexible fiber suspensions are shown to increase by a magnitude of the order 7–10% greater than the equivalent rigid fiber suspension tested. This is in qualitative agreement with previous experimental work by Goto et al. [Rheologica Acta 25 (1986) 119] and Blakeney [J. Colloid Interface Sci. 22 (1966) 324]. The implication is that any constitutive relation involving particulate suspensions described by orientation vectors may quantitatively underestimate suspension viscosity, particularly for fibers of large aspect ratio, or low Young’s modulus, whereby the tendency to flex is greater [Rheologica Acta 25 (1986) 119]. If particulate deformation were accounted for (by whatever means) in the existing constitutive relationship, predictions of bulk suspension parameters such as viscosity should be noticeably improved. A method is developed to modify an existing rigid-fiber viscosity to an equivalent flexible fiber viscosity, hence improving viscosity prediction ability.     

Shear thickening in electrically-stabilized colloidal suspensions

A theory is presented for the onset of shear thickening in colloidal suspensions of particles, stabilized by an electrostatic repulsion. Based on an activation model, a critical shear stress can be derived for the onset of shear thickening in dense suspensions for a constant potential and a constant charge approach of the spheres. Unlike previous models, the total interaction potential is taken into account (sum of attraction and repulsion). The critical shear stress is related to the maximum of the total interaction potential scaled by the free volume per particle. A comparison with experimental investigations shows the applicability of the theory.     

Constitutive equations for fiber suspensions in viscoelastic media

The kinetic theory of elastic dumbells with a friction factor that depends on the fiber orientation is used to obtain constitutive equations for fiber suspensions in a polymer matrix. We followed the approach of Fan (X.J. Fan, in P. Moldenaers and R. Keunings (Eds.), Theoretical and Applied Rheology, Proceedings XIth International Congress on Rheology, Brussels, Elsevier, Amsterdam, 1992, pp. 850–852), and derived equations for polymer solutions based on the FENE-P, FENE-CR, and Giesekus models. Start-up and steady-state free shear flows are studied to explore the effects of the fiber-polymer coupling as well as the fiber volume fraction. Predictions based on different types of closure approximations for the fourth-order fiber orientation tensor are also discussed.