Orientation tensors in simple flows of dilute suspensions of non-Brownian rigid ellipsoids,comparison of analytical and approximate solutions

General analytical solutions are obtained for the planar orientation structure of rigid ellipsoid of revolutions subjected to an arbitrary homogeneous flow in a Newtonian fluid. Both finite and infinite aspect ratio particles are considered. The orientation structure is described in terms of two-dimensional, time-dependent tensors that are commonly employed in constitutive equations for anisotropic fluids such as fiber suspensions. The effect of particle aspect ratio on the evolution of orientation structure is studied in simple shear and planar elongational flows. With the availability of analytical solutions, accuracies of quadratic closure approximations used for nonhomogeneous flows are analyzed, avoiding numerical integration of orientation distribution function. In general, fourth-order orientation evolution equations with sixth-order quadratic closure approximations yield more accurate representations compared to the commonly used second-order evolution equations with fourth-order quadratic closure approximations. However, quadratic closure approximations of any order are found to give correct maximum orientation angle (i.e., preferred direction) results for all particle aspect ratios and flow cases.     

车辆主动悬架设计的约束H∞控制方法

针对线性离散时间系统,研究其时域硬约束下的H∞输出反馈控制问题。假定外界扰动能量有限,基于线性矩阵不等式处理方法,提出并证明了时域硬约束下H∞输出反馈控制器存在的充分条件,并应用于车辆主动悬架设计。四分之一车辆模型数值仿真结果表明,即使车辆模型参数存在不确定性,本文提出的控制器在提高车辆乘坐舒适性的同时,仍能很好地兼顾车辆的其他性能要求。     

Hierarchical modeling in multibody dynamics

Summary In this paper a hierarchical approach using several mechanical models with different complexities and modeling depths to describe a single engineering system is presented. The mechanical models are derived from (but not limited to) multibody dynamics. The computer power available and improvements in theoretical understanding allow today not only to perform analyses but also to attack the problem of multimodel synthesis. Therefore, hierarchical modeling is used as a basis to analyze simultaneously models with different complexities and different excitations, and to optimize the performance with the most appropriate model for an investigated mechanical effect. Since only one single engineering system is investigated, its different models must be coupled by shared parameters, and the different criteria have to be combined with multicriteria optimization algorithms in order to obtain a single feasible design. An example taken from vehicle dynamics demonstrates the application of the approach. Received 14 January 1997; accepted for publication 11 September 1997     

Hydromagnetic transverse instability of two highly viscous fluid-particle flows with finite ion Larmor radius corrections

A linear analysis of the combined effect of viscosity, finite ion Larmor radius and suspended particles on Kelvin-Helmholtz instability of two superposed incompressible fluids in the presence of a uniform magnetic field is carried out. The magnetic field is assumed to be transverse to the direction of streaming. A general dispersion relation for such a configuration has been obtained using appropriate boundary conditions. The stability analysis is discussed analytically, and the obtained results are numerically confirmed. Some special cases are recovered and corrected. The limiting cases of absence of suspended particles (or fluid velocities) and finite Larmor radius, absence of suspended particles are discussed in detail. In both cases, all other physical parameters are found to have stabilizing as well as destabilizing effects on the considered system. In the former case, the kinematic viscosity is found to has a stabilizing effect, while in the later case, the finite Larmor radius is found to has a stabilizing influence for a vortex sheet. It is shown also that both finite Larmor radius and kinematic viscosity stabilizations for interchange perturbations are similar to the stabilization effect due to a magnetic field for non-interchange perturbations. Received 13 January 2003 Published online 24 April 2003 RID="a" ID="a"Also at: Department of Mathematics, Faculty of Education, Ain Shams University, Roxy, Cairo, Egypt. e-mail: m.elsayed@uaeu.ac.ae     

Manipulation of particles using dielectrophoresis

A numerical scheme based on the distributed Lagrange multiplier method (DLM) is used to study the motion of particles of electrorheological suspensions subjected to non-uniform electric fields. At small Reynolds number, the time taken by the particles to collect at the minimums or maximums of the electric field is primarily determined by a parameter defined to be the ratio of the dielectrophoretic and viscous forces. Simulations show that in non-uniform electric fields the collection time is also influenced by a parameter defined by the ratio of the electrostatic particle–particle interaction and dielectrophoretic forces. The collection time decreases as this parameter decreases because when this parameter is less than one, particles move to the regions of high or low electric field regions individually. However, when this parameter is greater than one, particles regroup into chains which then move toward the electric field maximums or minimums without breaking. It is also shown that when the real part of the Clausius–Mosotti factor (β) is negative the positions of the local minimums of the electric field, and thus also the locations where particles collect, can be modified by changing the electric potential boundary conditions.     

Simulations of fibre orientation in dilute suspensions with front moving in the filling process of a rectangular channel using level-set method

The simulation of fibre orientation in dilute suspension with front moving is carried out using the projection and level-set methods. The motion of fibres is described using the Jeffery equation, and the contribution of fibres to the flow is accounted for by the configuration-field method. The dilute suspension of short fibres in Newtonian fluids is considered. The governing Navier–Stokes equation for the fluid flow is solved using the projection method with finite difference scheme, while the fibre-related equations are directly solved with the Runge–Kutta method. In the present study for fibres in dilute suspension flow for injection molding, the effects of various flow and material parameters on the fibre orientation, the velocity distributions and the shapes of the leading flow front are found and discussed. Our findings indicate that the presence of fibre motion has little influence on the front shape in the ranges of fibre parameters studied at the fixed Reynolds number. Influence of changing fibre parameters only causes variation of front shape in the region near the wall, and the front shape in the central core area does not vary much with the fibre parameters. On the other hand, the fibre motion has strong influence on the distributions of the streamwise and transverse velocities in the fountain flow. Fibre motion produces strong normal stress near the wall which leads to the reduction of transversal velocity as compared to the Newtonian flow without fibres, which in turn, leads to the increased streamwise velocity near the wall. Thus, the fibre addition to the flow weakens the strength of the fountain flow. The Reynolds number has also displayed significant influence on the distribution of the streamwise velocity behind the flow front for a given fibre concentration. It is also found that the fibre orientation is not always along the direction of the velocity vector in the process of mold filling. In the region of the fountain flow, the fibre near the centreline is more oriented across the streamwise direction compared to that in the region far behind the flow front. This leads to the fact that the fibre near the centreline in the region of fountain flow is more extended along the transverse direction. As the fibre orientation in the suspension flow and the shape of the flow front have great bearing on the quality of the product made from injection molding, this study has much implications for engineering applications. These results can also be useful in other fields dealing with fibre suspensions.     

Stability of some shear flows for concentrated suspensions

In this paper we investigate the stability of some viscometric flows for a concentrated suspension model which allows for the effects of shear-induced migration, including plane and circular Couette and Poiseulle flows, and unbounded and bounded torsional flows. In the bounded torsional flow, where its radial outer boundary is assumed frictionless, an exact closeform solution is given. With the exception of torsional flows, we find that a limit point for all the steady-state solutions can exist for certain range in the parameter values. In all cases, disturbances can persist for a long time, O (H ²/a ²), where H is a dimension of the flow field, and a is the particles' radius.     

悬浮颗粒运动的格子Boltzmann数值模拟

将固体颗粒的牛顿力学和格子Boltzmann方法相结合,研究不规则形状悬浮颗粒在流场中的运动。通过受力分析,精确求得其所受合力、合力矩、合力作用中心等。提出了跟随颗粒运动的动网格计算域技术和模拟悬浮颗粒转动运动的局部数组方法及Euler-Lagrange两套坐标技术。通过对椭圆颗粒运动的数值模拟和对照他人对矩形颗粒的研究,分析了其复杂运动规律,并提供了合理的物理解释。结果表明:运用格子Boltzmann方法和上述特殊技术可以得到与有限元方法相同的模拟精度,且具有计算速度快、对复杂形状边界处理方便灵活、程序简单及特别适合大规模并行计算等优点。     

On the behaviour of gelled fibre suspensions in steady shear

The shear rheological properties of suspensions of gelled agar fibres in a low viscosity Newtonian matrix fluid were investigated. Two classes of fibres, low aspect ratio fibres and high aspect ratio fibres with an aspect ratio of the order of 10 and 100 respectively were included in the investigations. For all fibre phase volumes investigated, from as low as 0.01 upwards, the flow curves are characterised by an apparent yield stress followed by shear-thinning which was independent of the fibre aspect ratio. Based on our analysis of the flow curves, we conclude that the high aspect ratio fibres behave like flexible threads in contrast to the low aspect ratio fibres whose high shear relative viscosity is successfully described by a relation for long rigid rods. These findings are supported by flow visualisation using an optical shearing stage coupled to a light microscope.