胶体有效粘滞系数的理论研究

胶体是自然界和工业应用中的常见对象。胶体领域的一个中心理论问题是如何根据胶体系统的微结构来确定胶体的流变性质。由于处理多粒子系统边值问题的困难,现有的胶体理论都局限于低颗粒浓度。该文中发展了变换场方法,用该方法可以计算含胶体颗粒的不可压缩粘滞流体的有效粘滞系数,颗粒可以是固体也可以是流体。在低颗粒浓度,该理论预测与爱因斯坦关于悬浮体的公式以及Taylor关于乳浊液的公式完全吻合。在高颗粒浓度,该文的结果与Nunan和Keller的结果相一致,其方法可以用于预测非球形颗粒悬浮体的有效粘滞系数。     

On an effective model for ferronematic liquid crystals

The tools of homogenization theory are used to develop an effective model for ferronematics–colloidal suspensions of small ferromagnetic particles in a nematic liquid crystalline medium. The results are presented for dilute suspensions of identical prolate spheroidal particles. For needle-like particles of large eccentricity, the model reduces to a known expression [5]. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A continuum theory of dense suspensions

A continuum theory is introduced for viscous fluids carrying dense suspensions (such as blood) or emulsions of arbitrary shape and inertia. Suspended particles possess microinertia that make the mixture an anisotropic fluid whose viscosity changes with motion and orientation of suspensions. The microinertia balance law coupled with the equations of motion of an anisotropic fluid govern the ultimate outcome. By means of the second law of thermodynamics, constitutive equations are obtained in terms of the frame-independent tensors. In a special case, a theory of bar-like suspensions is obtained. The field equations, boundary and initial conditions are given for both the arbitrarily-shaped suspensions and the bar-like suspensions. The theory is demonstrated with the solution of the channel flow problem. The mean viscosity of the fluid with suspensions is determined. The motions of suspensions down flow are demonstrated.     

A continuum theory of dense suspensions

A continuum theory is introduced for viscous fluids carrying dense suspensions (such as blood) or emulsions of arbitrary shape and inertia. Suspended particles possess microinertia that make the mixture an anisotropic fluid whose viscosity changes with motion and orientation of suspensions. The microinertia balance law coupled with the equations of motion of an anisotropic fluid govern the ultimate outcome. By means of the second law of thermodynamics, constitutive equations are obtained in terms of the frame-independent tensors. In a special case, a theory of bar-like suspensions is obtained. The field equations, boundary and initial conditions are given for both the arbitrarily-shaped suspensions and the bar-like suspensions. The theory is demonstrated with the solution of the channel flow problem. The mean viscosity of the fluid with suspensions is determined. The motions of suspensions down flow are demonstrated.     

A numerical model for magnetic chromatography

In this paper, we present details of a mathematical model for magnetic chromatography (MC) systems where strong distorted magnetic fields are used to separate particles from a colloidal mixture. The model simulates the effect of magnetic field gradients on particle motion, and includes calculation of the fluid flow, magnetic field, and particle concentration field. It is based on the finite-volume method (FVM) and uses an expanding-grid technique to handle domains with large aspect ratios. The model has been validated against the results from an analytical model. The numerical model has been used to simulate the performance of a real MC system under various operating conditions.     

基于内变量理论的电流变液本构关系

研究了电流变液的微结构本构关系.其理论框架是基于内变量理论和机理的分析.电流变液是由高介电常数的颗粒悬浮在某种液体中组成的.在电场作用下,极化的颗粒将沿着电场方向聚集在一起形成链状结构.颗粒聚集体的大小和方向将随外加电场和应变率的变化进行调整,因而可以通过建立起能量守恒方程和力平衡方程来确定颗粒聚集体的大小和方向的变化.那么,一个三维的、清晰的本构关系可以由相互作用能和系统的耗散能导出.具体考虑和讨论了在简单剪切载荷作用下的系统响应,发现电流变液的切变剪薄粘滞系数同系统Mason数之间近似于幂指数∝(Mn)-082的关系.     

A general theoretical model for the magnetohydrodynamic flow of micropolar magnetic fluids. Application to Stokes flow

Many practical applications, which have an inherent interest of physical and mathematical nature, involve the hydrodynamic flow in the presence of a magnetic field. Magnetic fluids comprise a novel class of engineering materials, where the coexistence of liquid and magnetic properties provides us with the opportunity to solve problems with high mathematical and technical complexity. Here, our purpose is to examine the micropolar magnetohydrodynamic flow of magnetic fluids by considering a colloidal suspension of ferromagnetic material (usually non‐conductive) in a carrier magnetic liquid, which is in general electrically conductive. In this case, the ferromagnetic particles behave as rigid magnetic dipoles. Thus, the application of an external magnetic field, apart from the creation of an induced magnetic field of minor significance, will prevent the rotation of each particle, increasing the effective viscosity of the fluid and will cause the appearance of an additional magnetic pressure. Despite the fact that the general consideration consists of rigid particles of arbitrary shape, the assumption of spherical geometry is a very good approximation as a consequence of their small size. Our goal is to develop a general three‐dimensional theoretical model that conforms to physical reality and at the same time permits the analytical investigation of the partial differential equations, which govern the micropolar hydrodynamic flow in such magnetic liquids. Furthermore, in the aim of establishing the consistency of our proposed model with the principles of both ferrohydrodynamics and magnetohydrodynamics, we take into account both magnetization and electrical conductivity of the fluid, respectively. Under this consideration, we perform an analytical treatment of these equations in order to obtain the three‐dimensional effective viscosity and total pressure in terms of the velocity field, the total (applied and induced) magnetic field and the hydrodynamic and magnetic properties of the fluid, independently of the geometry of the flow. Moreover, we demonstrate the usefulness of our analytical approach by assuming a degenerate case of the aforementioned method, which is based on the reduction of the partial differential equations to a simpler shape that is similar to Stokes flow for the creeping motion of magnetic fluids. In view of this aim, we use the potential representation theory to construct a new complete and unique differential representation of magnetic Stokes flow, valid for non‐axisymmetric geometries, which provides the velocity and total pressure fields in terms of easy‐to‐find potentials, via an analytical fashion. Copyright © 2009 John Wiley & Sons, Ltd.     

Modeling,Mathematical and Numerical Analysis of Electrorheological Fluids

Many electrorheological fluids are suspensions consisting of solid particles and a carrier oil. If such a suspension is exposed to a strong electric field the effective viscosity increases dramatically. In this paper we first derive a model which captures this behaviour. For the resulting system of equations we then prove local in time existence of strong solutions for large data. For these solutions we finally derive error estimates for a fully implicit time-discretization.     

Mathematical modeling for colloidal dispersion undergoing Brownian motion

An averaged motion approach for modeling Brownian dynamics for suspension systems of electrically charged particles in liquid is developed. The continuum model for the motion of particles consists of a system of integral equations coupled with a degenerate parabolic equation. Existence and uniqueness of global solution for the coupled system are established, and numerical results for the non-Newtonian viscosity of the mixture in terms of shear rate or Pechlet number are obtained. The model reveals some non-Newtonian properties such as the well-known shear thinning phenomenon for the viscosity of colloidal dispersions.     

Application of the asymptotic solution to EM field scattering problem for creation of media with prescribed permeability

Scattering of electromagnetic (EM) waves by many small impedance particles (bodies), embedded in a homogeneous medium, is studied. Physical properties of the particles are described by their boundary impedances. The limiting equation is obtained for the effective EM field in the limiting medium, in the limit a→0, where a is the characteristic size of a particle and the number M(a) of the particles tends to infinity at a suitable rate. The proposed theory allows one to create a medium with a desirable spatially inhomogeneous permeability. The main new physical result is the explicit analytical formula for the permeability μ(x) of the limiting medium. The computational results confirm a possibility to create the media with various distributions of μ(x).     

在两相粘性跨音速喷管流动中簿层方程的一种隐式求解方法*

本文略去沿流动方向的粘性,将任意曲线坐标系中无量纲化的N-S方程简化为薄层方程.采用隐式近似因子分解法求解气相控制方程,采用特征线法跟踪颗粒,然后获得两相跨音速湍流充分耦合的数值方法.其中,颗粒尺寸是分级的,用参考平面中的拟特征线法处理喷管的粘性亚音速进口边界条件,湍流采用代数模型.该计算方法应用于火箭喷管两相粘流计算,并预估了固体火箭发动机的推力和比冲,计算与试验结果吻合很好.文中还讨论了不同颗粒尺寸、不同颗粒质量百分数和颗粒尺寸分级等对流场的影响,分析了颗粒、二维径向分速和粘性对发动机比冲的影响.本文的方法具有节省机时的优点,尤其是对颗粒尺寸分级的计算,效果更为显着.     

两球形颗粒间横向毛细力的格子Boltzmann研究

采用以Shan-Chen多组分模型为基础的格子Boltzmann-伪固体模型对两颗粒间的浮体和浸润横向毛细力展开数值模拟研究,其中流体-固体间的相互作用及颗粒润湿性质在介观层次上采用简单形式得以充分考虑．三维测试表明,与已有理论解相比,成功再现了横向毛细力与颗粒间距的“1/L”关系,并确认了浸润横向毛细力与表面张力间的线性关系．这表明可进一步应用该模型研究横向毛细力作用下的颗粒自聚集等现象.     

NSM solution for unsteady MHD Couette flow of a dusty conducting fluid with variable viscosity and electric conductivity

The effects of dependence on temperature of the viscosity and electric conductivity, Reynolds number and particle concentration on the unsteady MHD flow and heat transfer of a dusty, electrically conducting fluid between parallel plates in the presence of an external uniform magnetic field have been investigated using the network simulation method (NSM) and the electric circuit simulation program Pspice. The fluid is acted upon by a constant pressure gradient and an external uniform magnetic field perpendicular is applied to the plates. We solved the steady-state and transient problems of flow and heat transfer for both the fluid and dust particles. With this method, only discretization of the spatial co-ordinates is necessary, while time remains as a real continuous variable. Velocity and temperature are studied for different values of the viscosity and magnetic field parameters and for different particle concentration and upper wall velocity.     

Separation of magnetic particles in channel flows by BEM

In-line separation of suspensions can become difficult in case of particles with comparable values of densities. For flows in micro devices in such cases gravitational settling is inefficient, and other separation techniques must be applied. In case of magneto active particles, the action of Kelvin magnetic force in a non-uniform magnetic field could be used in order to achieve a higher degree of particles separation. The contribution therefore deals with Euler-Lagrangian formulation of dilute two-phase flows. The Boundary element based computational algorithm solves the incompressible Navier-Stokes equations written in velocity-vorticity formulation. The non-uniform magnetic field is defined analytically for the case of a set of long thin wires. The particle trajectories are computed by applying the 4th order Runge-Kutta method. The computed test case consists of a narrow channel with laminar flow of suspension under Re = 1 − 10. Particle trajectories under the influence of a non-uniform magnetic field are computed for the case of magnetite and aluminium particles suspended in water. The efficiency of separation on basis of particle trajectories for different values of Re number and magnetic field strength is performed, clearly indicating superior separation of magneto active particles. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Comparison of non-cohesive resolved and coarse grain DEM models for gas flow through particle beds

The Discrete Element Method (DEM) is a widely used approach for modelling granular systems. Currently, the number of particles which can be tractably modelled using DEM is several orders of magnitude lower than the number of particles present in common large-scale industrial systems. Practical approaches to modelling such industrial system therefore usually involve modelling over a limited domain, or with larger particle diameters and a corresponding assumption of scale invariance. These assumption are, however, problematic in systems where granular material interacts with gas flow, as the dynamics of the system depends heavily on the number of particles. This has led to a number of suggested modifications for coupled gas–grain DEM to effectively increase the number of particles being simulated. One such approach is for each simulated particle to represent a cluster of smaller particles and to re-formulate DEM based on these clusters. This, known as a representative or ‘coarse grain’ method, potentially allows the number of virtual DEM particles to be approximately the same as the real number of particles at relatively low computational cost. We summarise the current approaches to coarse grain models in the literature, with emphasis on discussion of limitations and assumptions inherent in such approaches. The effectiveness of the method is investigated for gas flow through particle beds using resolved and coarse grain models with the same effective particle numbers. The pressure drop, as well as the pre and post fluidisation characteristics in the beds are measured and compared, and the relative saving in computational cost is weighed against the effectiveness of the coarse grain approach. In general, the method is found perform reasonably well, with a considerable saving of computational time, but to deviate from empirical predictions at large coarse grain ratios.     

悬浮流中柱状粒子取向分布函数的数值模拟

推导并数值模拟了悬浮流中柱状粒子取向的分布函数,并将其用于楔形流场计算,得到了取向分布函数与粒子方向角的关系。说明与短时间内形成的最可几角分布相比,稳定状态的分布变化幅度不大,主要变化在于右上区域的逆时针转动,速度梯度值大的地方,转动角度也越大;粒子最可几方向接近于流线方向:在流线方向上,随着极径减小,最可几角增大,同时它与流线夹角减小。     

Mathematical modeling of the movement of suspended particles subjected to acoustic and flow fields

When particles are subjected to an acoustic field particle trajectories depend on the particle and fluid compressibility and density values. Hence a combination of acoustic and flow fields on particles can be used to deflect and trap, or to segregate and/or fractionate fine particles in fluid suspensions. Using particle physics in an acoustic field, a mathematical model was developed to calculate trajectories of deflected particles due to the application of acoustic standing waves. The resulting second order ordinary differential equation was quite stiff and hence difficult to solve numerically and did not have a closed form solution. The analysis of the above equation showed that the basic problems with numerical solutions could not be ameliorated through the use of standard rescaling techniques. A combination of phase space and asymptotic analysis turns out to be far more useful in obtaining approximate solutions. An approximate solution was derived which enabled the calculation of the particle trajectories and concentration at collection planes in the acoustic field. Analysis of the solution showed that all the particles move toward the pressure node to which the particles are supposed to move. Particles with 2 μm diameter took approximately 20 s to reach that node. Then at the bench scale, the above technology was implemented by building a flow chamber with two transducers at opposite ends to generate an acoustic standing wave. SiC particle trajectories were tracked using captured digital images from a high-resolution microscope. The displacements of SiC particles due to an acoustic force were compared with the mathematical model predictions. For input power levels between 3.0 and 5.0 W, the experimental data were comparable to mathematical model predictions. Hence it was concluded that the proposed approximate solution was both quantitatively and qualitatively closer to experimental results than the simplified form ignoring the second order term reported in the literature.     

On some simple examples of mechanical systems with hyperbolic chaos

Examples of mechanical systems with hyperbolic chaos are discussed, including the Thurston–Weeks–Hunt–MacKay hinge mechanism, in which conservative Anosov dynamics is realized, and dissipative systems with Smale–Williams type attractors (a particle on a plane under periodic kicks, interacting particles sliding on two alternately rotating disks, and a string with parametric excitation by modulated pump). The examples considered in the paper are interesting from the viewpoint of filling hyperbolic theory, as a well-developed field of the mathematical theory of dynamical systems, with physical content. The results of computer tests for hyperbolicity of the systems are presented that are based on the analysis of the statistics of intersection angles of stable and unstable manifolds.     

What is a particle from the standpoint of systems theory?

Our theory of collective motions of open system [2,3] is generalized to the case of systems distributed inm-dimensional space. We propose a rigorous definition of a particle in terms of systems theory. A concrete image of a particle, conceived as a localized collective entity in Space, is obtained. It turns out that wave-corpuscle duality is nothing else but a simple consequence of general equations of collective motions: particles appear at the output as collective manifestations of inner states, and waves appear at the input as guiding waves of particles.The David and Helena Zlotowski chair in Operator Theory and Systems Theory.     

The Green formula in the theory of potential scattering and corrections to the eikonal scattering amplitude

A method based on the Green formula is developed for calculating the scattering amplitude of fast charged particles in an external field. The scattering amplitude is representable as an integral over an arbitrary closed surface enveloping the domain of influence of the external field on the particle. Corrections to the eikonal scattering amplitude are simply derived without using the specific form of the potential. The resulting formulas can be used to investigate the interaction between particles and fields of complex configuration. Translated from Teoreticheskaya i Matematicheskaya Fizika. Vol. No. 2, pp. 280–288, May, 1998.