Viscosity of a Suspension with a Cubic Array of Spheres in a Shear Flow

Viscous fluid flow past an infinite periodic array of rigid spheres of the same radius is considered. A solution of the Stokes equations periodic in three variables is obtained for viscous incompressible flow with a linear velocity profile. The solution takes into account the hydrodynamic interaction of an infinite number of particles in the array. An expression for the effective viscosity of a suspension with a cubic array of particles is obtained.     

Symmetry breaking and preliminary results about a Hopf bifurcation for incompressible viscous flow in an expansion channel

This computational study shows, for the first time, a clear transition to two-dimensional Hopf bifurcation for laminar incompressible flows in symmetric plane expansion channels. Due to the well-known extreme sensitivity of this study on computational mesh, the critical Reynolds numbers for both the known symmetry-breaking (pitchfork) bifurcation and Hopf bifurcation were investigated for several layers of mesh refinement. It is found that under-refined meshes lead to an overestimation of the critical Reynolds number for the symmetry breaking and an underestimation of the critical Reynolds number for the Hopf bifurcation.     

Development of a doubly periodic system of cavities in viscous bodies

The plane unknown-boundary problem of the development of a doubly periodic system of cavities in viscous media in the presence of finite strains is considered. Under conditions of timedependent slow Newtonian viscous fluid flow the solution of the doubly periodic problem of the development of a system of identical cavities whose centers are arranged in square and triangular grids is obtained.     

Numerical simulation of plate autorotation in a viscous fluid flow

A general formulation of the plane coupled dynamical and aerodynamical problem of the motion of a rigid body with a rotational degree of freedom in a viscous incompressible fluid flow is given. A computation technique for solving the Navier-Stokes equations based on the meshless viscous vortex domain method is used. The autorotation of a single plate and a pair of plates is investigated. The effect of the reduced moment of inertia and the Reynolds number on the angular rotation velocity is determined. The time dependences of the hydrodynamic loads are compared with the corresponding instantaneous flow patterns. The increased the autorotation velocity of two plates in tandem is detected.     

黏性不可压缩流体流动前沿的数值模拟

提出了模拟注射成型中黏性、不可压缩流体流动前沿的新方法. 将Hele-Shaw流动应用于非 等温条件下的黏性、不可压缩流体，建立了流动分析模型，用充填因子的输运方程描述流动 前沿. 应用高阶Taylor展开式计算每一时间步长的充填因子，用Galerkin方法导出了计算 充填因子各阶导数的递推公式. 给出了时间增量的选取方法，证明了它的稳定性. 针对Han 设计的试验模具，用相同的材料及工艺条件模拟充填过程，比较了传统方法和该方法的模 拟结果与实验结果的差异. 算例分析表明，该方法可以有效地提高注射成型中流动前沿的 模拟精度和计算效率.     

The generation of three-dimensional internal waves and attendant boundary layers in a viscous continuously stratified fluid. Construction of an analytical solution

An analytical solution of a linearized problem of the emission of periodic internal waves by part of a plane which oscillates with a small amplitude in an arbitrary direction in a viscous exponentially stratified fluid is constructed. Solutions of the dispersion equation are given for all positions of the emitting surface (arbitrary, vertical, horizontal, and critical when one of the beam propagation directions is collinear with the emitting surface). The possibility of transition to the case of a uniform fluid, which is important for applications, is analyzed.     

DAMPING OF VERTICALLY EXCITED SURFACE WAVE IN WEAKLY VISCOUS FLUID

In a vertically oscillating circular cylindrical container, singular perturbation theory of two-time scale expansions is developed in weakly viscous fluids to investigate the motion of single free surface standing wave by linearizing the Navier-Stokes equation. The fluid field is divided into an outer potential flow region and an inner boundary layer region. The solutions of both two regions are obtained and a linear amplitude equation incorporating damping term and external excitation is derived. The condition to appear stable surface wave is obtained and the critical curve is determined. In addition, an analytical expression of damping coefficient is determined. Finally, the dispersion relation, which has been derived from the inviscid fluid approximation, is modified by adding linear damping. It is found that the modified results are reasonably closer to experimental results than former theory. Result shows that when forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when forcing frequency is high, the surface tension of the fluid is prominent.     

Onset of self-oscillations upon the loss of stability of spatially periodic two-dimensional viscous fluid flows relative to long-wave perturbations

The long-wave asymptotics of the secondary flow that arises after a steady, spatially periodic flow loses stability are studied when one of the periods tends to infinity and the rate of base flow along the longer period is equal to zero. It is shown that if certain non-degeneracy conditions are satisfied, then from the base solution a self-oscillatory regime branches off and both hard and soft stability loss is possible. For the leading terms of the asymptotics explicit formulas are obtained. Examples of self-oscillations calculated for specific flows are presented and the behavior of the fluid particle trajectories in the self-oscillatory regime branching off from the base flow is investigated.     

Bifurcation condition of crack pattern in the periodic rectangular array

Bifurcation condition of crack pattern in the periodic rectangular array plays an important role in determining the final failure pattern of rock material. An approximation for the critical crack size/spacing ratio is established for a uniformly growing periodic rectangular array yields to a non-uniform growing pattern of crack growth. Numerical results show that the critical crack size/spacing ratio λ_cr depends on the number of cracks, the crack spacing, the perpendicular distance between two adjacent rows, as well as the loading conditions. In general, λ_cr increases with the number of lines. It is observed that the critical crack size/spacing ratio λ_cr for the periodic rectangular array decreases with an increase in the perpendicular distance between two adjacent rows. It is clear that the critical crack size/spacing ratio λ_cr for the periodic rectangular array under shear stress increases with increasing the crack spacing.     

Simulation of interface and free surface flows in a viscous fluid using adapting quadtree grids

A new adaptive quadtree method for simulating laminar viscous fluid problems with free surfaces and interfaces is presented in this paper. The Navier–Stokes equations are solved with a SIMPLE‐type scheme coupled with the Compressive Interface Capturing Scheme for Arbitrary Meshes (CICSAM) (Numerical prediction of two fluid systems with sharp interfaces, Ph.D. Thesis, Imperial College of Science, Technology and Medicine, London, 1997) volume of fluid (VoF) method and PLIC reconstruction of the volume fraction field during refinement and derefinement processes. The method is demonstrated for interface advection cases in translating and shearing flow fields and found to provide high interface resolution at low computational cost. The new method is also applied to simulation of the collapse of a water column and the results are in excellent agreement with other published data. The quadtree grids adapt to follow the movement of the free surface, whilst maintaining a band of the smallest cells surrounding the surface. The calculation is made on uniform and adapting quadtree grids and the accuracy of the quadtree calculation is shown to be the same as that made on the equivalent uniform grid. Copyright © 2004 John Wiley & Sons, Ltd.     

MOTION AND DEFORMATION OF VISCOUS DROP IN STOKES FLOW NEAR RIGID WALL

A boundary integral method was developed for simulating the motion and deformation of a viscous drop in an axisymmetric ambient Stokes flow near a rigid wall and for direct calculating the stress on the wall. Numerical experiments by the method were performed for different initial stand-off distances of the drop to the wall, viscosity ratios, combined surface tension and buoyancy parameters and ambient flow parameters. Numerical results show that due to the action of ambient flow and buoyancy the drop is compressed and stretched respectively in axial and radial directions when time goes. When the ambient flow action is weaker than that of the buoyancy the drop raises and bends upward and the stress on the wall induced by drop motion decreases when time advances. When the ambient flow action is stronger than that of the buoyancy the drop descends and becomes flatter and flatter as time goes. In this case when the initial stand-off distance is large the stress on the wall increases as the drop evolutes but when the stand-off distance is small the stress on the wall decreases as a result of combined effects of ambient flow, buoyancy and the stronger wall action to the flow. The action of the stress on the wall induced by drop motion is restricted in an area near the symmetric axis, which increases when the initial stand-off distance increases. When the initial stand-off distance increases the stress induced by drop motion decreases substantially. The surface tension effects resist the deformation and smooth the profile of the drop surfaces. The drop viscosity will reduce the deformation and migration of the drop.     

Stability of a viscoelastic plate in fluid flow

The stability of an infinite viscoelastic plate on an elastic foundation in a viscous incompressible flow is studied. The Navier-Stokes system is linearized for an exponential velocity profile. The problem is reduced by a Fourier-Laplace transform to a system of ordinary differential equations, whose solution is found in the form of convergent series. The roots of the dispersion relation that characterize the stability of the system are found numerically. The effect of the viscosities of the fluid and the plate on the stability of the waves propagating upstream and downstream is studied. The results are compared with available data on the stability of a viscoelastic plate in an ideal fluid flow. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 4, pp. 66–74, July–August, 2006.     

Three-dimensional motion of a viscous incompressible fluid in a narrow tube

An approximate solution of the problem of unsteady motion of a viscous incompressible fluid in a long narrow deformable tube at low Reynolds numbers is obtained. Pressure oscillations and tube deformation are shown to be related by an integrodifferential equation. The solution obtained extends the Poiseuille solution in elliptic tubes to the case of comparatively arbitrary small deformations in terms of the tube length and angle. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 4, pp. 28–32, July–August, 2009.     

VELOCITY FIELD IN SHIP WAVES ON THE VISCOUS FLUID

ＩｎｔｒｏｄｕｃｔｉｏｎＷｈｅｎａｍｏｖｉｎｇｂｏｄｙ (ｆｏｒｅｘａｍｐｌｅ:ｓｈｉｐ)ｉｓｉｎｍｏｔｉｏｎｏｎｔｈｅｏｃｅａｎ ,ａＶ＿ｓｈａｐｅｄｓｈｉｐｗａｖｅａｐｐｅａｒｓｂｅｈｉｎｄｔｈｅｂｏｄｙ .Ｉｎｔｈｉｓａｒｔｉｃｌｅ,ｗｅｗｉｌｌｄｉｓｃｕｓｓｔｈｅｓｈａｐｅｏｆｔｈｅｖｅｌｏｃｉｔｙｆｉｅｌｄａｎｄｔｈｅＶ＿ｓｈａｐｅｄｓｈｉｐｗａｖｅ .Ａｔｆｉｒｓｔ,Ｒｅｆ.[1]ｓｔｕｄｉｅｄｔｈｅｑｕｅｓｔｉｏｎｓｏｆｓｈｉｐｗａｖｅｓｏｎａｎｉｎｃｏｍｐｒｅｓｓｉｂｌｅｉｎｖｉｓｃｉｄｆｌｕｉ…     

空腔内粘性流问题与哈密顿体系方法

以双板驱动空腔粘性流问题为研究对象,根据其特点建立了哈密顿体系下的对偶正则方程,将问题归结为辛体系下的本征值问题.利用辛本征解空间的完备性、正交性和展开理论,形成一套封闭的求解问题方法.算例的数值结果揭示了一些空腔流动的特点.同时这种方法也为研究其他问题提供了一条思路.     

Optimum shape of a cavity for the collection of a soil-saturating viscous fluid

The process of inflow of an oil-bearing fluid from an infinite region of porous soil into a cavity (reservoir) and its subsequent pumping into a well located at the center of the reservoir is investigated. The cavity is chosen from the class of flattened ellipsoids of revolution. Thus, a combined problem of the optimization of the shape of a cavity of given volume for which the maximum seepage flow is achieved and the approximate determination of the additional constraints on the cavity dimensions necessary to ensure the outflow of the oil-bearing fluid is considered.     

Fallopian tube assessment of the peristaltic-ciliary flow of a linearly viscous fluid in a finite narrow tube

The present theoretical assessment deals with the peristaltic-ciliary transport of a developing embryo within a fallopian tubal fluid in the human fallopian tube. A mathematical model of peristalsis-cilia induced flow of a linearly viscous fluid within a fallopian tubal fluid in a finite two-dimensional narrow tube is developed. The lubrication approximation theory is used to solve the resulting partial differential equation. The expressions for axial and radial velocities, pressure gradient, stream function, volume flow rate, and time mean volume flow rate are derived. Numerical integration is performed for the appropriate residue time over the wavelength and the pressure difference over the wavelength. Moreover, the plots of axial velocity, the appropriate residue time over wavelength, the vector, the pressure difference over wavelength, and the streamlines are displayed and discussed for emerging parameters and constants. Salient features of the pumping characteristics and the trapping phenomenon are discussed in detail. Furthermore, a comparison between the peristaltic flow and the peristaltic-ciliary flow is made as the special case. Relevance of the current results to the transport of a developing embryo within a fallopian tubal fluid from ampulla to the intramural in the fallopian tube is also explored. It reveals the fact that cilia along with peristalsis helps to complete the required mitotic divisions while transporting the developing embryo within a fallopian tubal fluid in the human fallopian tube.     

FLOW OF A TRAIN OF DEFORMABLE FLUID PARTICLES IN A TUBE

I.IntroductionWangandSkalakll]studiedforsolidspheresintubes.HymanandSkalak12--31studiedthemovementofacompactfluidspheres(non-deformable)intubes.Theyonlystudiedcompactdropswhoseradiusesarenotbiggerthan0.7.Pozrikidisl4]studiedthenlotiollofdeformabletluidparticlesintubes.Inthispaper,westudiedtwo-phasenowregime,consistingofaperiodictrainofequallyspacedliquidslugsmovingaxisymmetricallyinacontinuousliquidphasewhichisimmisciblewiththeslugphase(Fig.l(a)).Wereportnumericalsolutionsforthehydrodynamic…     

Asymptotic models of solidification in cooling thin-layer flows of a highly viscous fluid

Asymptotic models are constructed for the solidification process in a highly viscous film flow on the surface of a cone with a given mass supply at the cone apex. In the thin-layer approximation, the problem is reduced to two parabolic equations for the temperatures of the liquid and the solid coupled with an ordinary differential equation for the solidification front. For large Péclet numbers, an analytical steady-state solution for the solidification front is found. A nondimensional parameter which makes it possible to distinguish flows (i) without a solid crust, (ii) with a steady-state solid crust, and (iii) with complete solidification is determined. For finite Péclet numbers and large Stefan numbers, an analytical transient solution is found and the time of complete flow solidification is determined. In the general case, when all the governing parameters are of the order of unity, the original system of equations is studied numerically. The solutions obtained are qualitatively compared with the data of field observations for lava flows produced by extrusive volcanic eruptions.