Investigation of Boundary-Value Problem for Stationary System of Equations of Viscous Non-Isothermal Fluid

In the Stokes approximation at small Reynolds and Peclet numbers, we obtain a solution to the boundary-value problem of flow around of particles of spherical shape for stationary system of equations of a viscous non-isothermal fluid comprising a linearized by speed Navier–Stokes equation system and the equation of heat transfer given an exponential-power law of dependence of viscosity of fluid on temperature.     

The quasistationary approximation in the problem of the motion of an isolated volume of a viscous incompressible capillary liquid

The equations of the quasistationary approximation in the problem of the motion of an isolated volume of a viscous incompressible capillary liquid are derived from the exact equations using an expansion in a small quasistationary parameter, which is equal to the ratio of the Stokes time to the capillary time. The problem contains yet another dimensionless parameter, which is proportional to the modulus of the conserved angular momentum of the liquid volume, which is also assumed to be small. Depending on the relation between these parameters, three versions of the limiting problem are obtained: the traditional version and two new versions. Asymptotic solutions of the problems which arise when the quasistationary parameter tends to zero are constructed.     

解任意形状扁轴对称体Stokes流动的连续奇点分布法

本文以Sampson球形无穷级数作为基本奇点,采用分段等强度和分段二次抛物分布两种体内连续分布法解任意形状扁轴对称体的Stokes流动.通过扁球的无界绕流问题,对这两种方法的收敛性,精度和适用范围做了检验和比较.结果表明,在一定的范围内,无论是阻力系数或压力分布,它们的计算结果都和精确解符合得很好,而且,随着分布函数逼近程度的提高,其收敛性得到改善,适用范围也随之扩大.作为一般算例,分别用这两种方法解决了卡西尼扁卵形体的绕流问题,得到了一致的结果.最后,用分段二次连续分布法计算了具有一定生理意义的红细胞体的Stokes流动,首次得到了它的阻力系数和表面压力分布.     

Longitudinal waves in partially saturated porous media: the effect of gas bubbles

The problem of the propagation of longitudinal waves in a liquid-saturated porous medium when there are gas bubbles present is considered. The decay factor and the phase velocity of Frenkel–Biot waves of the first and second kind are found as a function of the frequency in the linear approximation. It is shown that, in the neighbourhood of the resonance frequency of the bubbles, longitudinal Frenkel–Biot waves change their form. A wave of the first kind is transformed from a fast wave at low frequencies into a slow wave at high frequencies. The dispersion curve of a wave of the second kind consists of two branches – a “low-frequency” branch, the oscillations of which possess the classical properties, and a “high-frequency” branch, which is a weakly decaying high-velocity mode. The frequency dependences of the ratio of the mass velocities of a gas-liquid mixture and of a porous matrix, and also of the perturbations of the stress in the matrix and the pressure in the mixture, are constructed. It is shown that the “high-frequency” branch of a wave of the second kind is characterized by the in phase motion of the gas-liquid mixture and of the porous matrix, while their mass velocities are close, which explains the weak decay of this mode of oscillations. An analytical expression is obtained for the “boundary frequency”, which determines the offset of the “high-frequency” branch of the dispersion curve of the wave of the second kind.     

A new approach to the classical Stokes flow problem: Part I Methodology and first-order analytical results

The problem of determining the axisymmetric Stokes flow past an arbitrary body, the boundary shape of which can be represented by an analytic function, is examined by developing an exact method. An appropriate nonorthogonal coordinate system is introduced, and it is shown that the Hilbert space to which the stream function belongs is spanned by the set of Gegenbauer polynomials based on the physical argument that the drag on a body should be finite. The partial differential equation of the original problem is then reduced to two simultaneous vector differential equations. By the truncation of this infinite-dimensional system to the one-dimensional subspace, an explicit analytic solution to the Stokes equation valid for all bodies in question is obtained as a first approximation.     

A defect correction method for the time‐dependent Navier‐Stokes equations

A method for solving the time dependent Navier‐Stokes equations, aiming at higher Reynolds' number, is presented. The direct numerical simulation of flows with high Reynolds' number is computationally expensive. The method presented is unconditionally stable, computationally cheap, and gives an accurate approximation to the quantities sought. In the defect step, the artificial viscosity parameter is added to the inverse Reynolds number as a stability factor, and the system is antidiffused in the correction step. Stability of the method is proven, and the error estimations for velocity and pressure are derived for the one‐ and two‐step defect‐correction methods. The spacial error is O(h) for the one‐step defect‐correction method, and O(h²) for the two‐step method, where h is the diameter of the mesh. The method is compared to an alternative approach, and both methods are applied to a singularly perturbed convection–diffusion problem. The numerical results are given, which demonstrate the advantage (stability, no oscillations) of the method presented. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Surface waves due to blasts on and above inviscid liquids of finite depth

For the problem of waves due to an explosion above the surface of a homogeneous ocean of finite depth, asymptotic expressions of the velocity potential and the surface displacement are determined for large times and distances from the pressure area produced by the incident shock. It is shown that the first item in Sakurai's approximation scheme for the pressure field inside the, blast wave as well as the results of Taylor's point blast theory can be used to yield realistic expressions of surface displacement. Some interesting features of the wave motion in general are described. Finally some numerical calculations for the surface elevation were performed and included as a particular case.     

Fluid flow in an inhomogeneous granular medium

The seepage of a compressible fluid in an inhomogeneous undeformable granular medium is investigated. It is assumed that the fluid flow in a porous space is described by the Navier–Stokes equations. It is shown that, in the case of an inhomogeneous velocity field, a tensor of additional effective stresses occurs in connection with the transfer of fluid particles in a transverse direction when flow occurs around the granules of the medium in a longitudinal direction. Using the fundamental propositions of Reynolds’ averaging theory and Prandtl's mixing path, the structure of the effective viscosity coefficient is determined and hypotheses are formulated which enable it to be assumed to be independent of the flow velocity. It is established by comparison with experimental data that the effective viscosity coefficient can exceed the viscosity coefficient of the flowing fluid by an order of magnitude. The equations of average motion are obtained, which in the case of an incompressible fluid have the form of the Navier–Stokes equations with body forces proportional to the velocity. It is established that, in addition to the well-known dimensionless flow numbers, there is a new number which characterizes the ratio of the Darcy porous drag forces to the effective viscosity forces. The proposed equations are extended to the case of the flow of an aerated fluid. The components of the angular momentum vector are used as the required functions instead of the components of the velocity vector. This enables a solving system of equations to be obtained, which, apart from the notation, is identical with the similar equations for the case of an incompressible fluid. The solution of a new problem of the fluid flow in a plane channel with permeable walls is presented using three models: Darcy's law for an incompressible and aerated fluid, and also of an aerated fluid taking the effective viscosity into account. It is established that, for the same pressure drop, the maximum flow rate corresponds to Darcy's law. Compressibility leads to its reduction, but by simultaneously taking into account the compressibility and the effective viscosity one obtains minimum values of the flow rate. The effective viscosity and aeration of the fluid has a considerable effect on the flow parameters.     

A multilevel correction method for Stokes eigenvalue problems and its applications

In this paper, a new multilevel correction scheme is proposed to solve Stokes eigenvalue problems by the finite element method. This new scheme contains a series of correction steps, and the accuracy of eigenpair approximation can be improved after each step. In each correction step, we only need to solve a Stokes problem on the corresponding fine finite element space and a Stokes eigenvalue problem on the coarsest finite element space. This correction scheme can improve the efficiency of solving Stokes eigenvalue problems by the finite element method. As applications of this multilevel correction method, a multigrid method and an adaptive finite element technique are introduced for Stokes eigenvalue problems. Some numerical results are given to validate our schemes. Copyright © 2013 John Wiley & Sons, Ltd.     

The self-propulsion of a body with moving internal masses in a viscous fluid

An investigation of the characteristics of motion of a rigid body with variable internal mass distribution in a viscous fluid is carried out on the basis of a joint numerical solution of the Navier — Stokes equations and equations of motion for a rigid body. A nonstationary three-dimensional solution to the problem is found. The motion of a sphere and a drop-shaped body in a viscous fluid in a gravitational field, which is caused by the motion of internal material points, is explored. The possibility of self-propulsion of a body in an arbitrary given direction is shown.     

Exact and approximate expressions for resistance coefficients of a colloidal sphere approaching a solid surface at intermediate Reynolds numbers

A mathematical model has been developed to describe the force of liquid flow acting on a colloidal spherical particle as it approaches a solid surface at intermediate-Reynolds-number-flow regime. The model has incorporated bispherical coordinates to determine a stream function for the flow disturbed by the sphere. The stream function was then used to derive the flow force on the particle as a function of the inter-surface separation distance. The force equation was related to the modified Stokes equation to obtain an exact analytical expression for the correction factor to the Stokes law. Finally, a rational approximation is presented, which is in good agreement with the exact numerical result, and can be readily applied to more general particle–surface interactions involving short-range hydrodynamics associated with colloidal particles in the near vicinity of a large solid collector surface at intermediate Reynolds number of the supporting flow.     

The generation of beams of three-dimensional periodic internal waves in an exponentially stratified fluid

The spectral method is used to construct an exact solution of the linearized problem of the generation of disturbances by localized sources that execute arbitrary periodic motions in a viscous exponentially stratified fluid. The expressions obtained do not contain any adjusting parameters and describe conical beams of three-dimensional periodic internal waves and two types of boundary layers, the spatial scale of which is given by the kinematic viscosity and the buoyancy frequency of the medium. The thickness of one of them, which is analogous to Stokes periodic flow in a homogeneous viscous fluid, is specified by the kinematic viscosity and the wave frequency, that is, it additionally depends on a ratio of the wave and buoyancy frequencies. The thickness of the specific internal boundary layer also depends on the geometry of the problem. In the approximation of weak stratification and low viscosity, asymptotic estimates of the expressions obtained are presented for two types of generators, namely, in the form of a plane inclined rectangle that vibrates along its surface (a frictional source) and along the normal to it (a piston source) in the non-degenerate case when the wave cone does not touch the radiating plane. In limiting cases the analytical expressions obtained agree with known exact solutions of the problem of generating axially symmetric and two-dimensional periodic internal waves.     

A self-consistent field method applied to the dynamics of viscous suspensions

A method for the approximate solution of the problem of many bodies of spherical form in a viscous fluid is developed in the Stokes approximation. Using a purely hydrodynamic approach, based on the use of the concept of a self-consistent field, the classical boundary value problem is reduced to a formal procedure for solving a linear system of algebraic equations in the tensor coefficients, which occur in the solution obtained for the velocity field and pressure of the liquid. A procedure for the approximate solution of this system of equations is constructed for the case of dilute suspensions, when the ratio of the size of the dispersed particles to the characteristic distance between them is a small parameter. Finally, the initial boundary value problem is reduced to solving a recurrent system of equations, in which each subsequent approximation for all the required quantities depends solely on the previous approximations. The system of recurrent equations obtained can be solved analytically in any specified approximation with respect to a small parameter. It is shown that this system of equations contains in itself all possible physical formulations of the problems, and, within the frameworks of the mathematical procedure constructed, they are distinguished solely by a set of specified and required functions. The practical possibilities of the method are in no way limited by the number of dispersed particles in the fluid.     

Exact controllability of a fluid–rigid body system

This paper is devoted to the controllability of a 2D fluid–structure system. The fluid is viscous and incompressible and its motion is modelled by the Navier–Stokes equations whereas the structure is a rigid ball which satisfies Newton's laws. We prove the local null controllability for the velocities of the fluid and of the rigid body and the exact controllability for the position of the rigid body. An important part of the proof relies on a new Carleman inequality for an auxiliary linear system coupling the Stokes equations with some ordinary differential equations.     

Numerical analysis of a Picard multilevel stabilization of mixed finite volume method for the 2D/3D incompressible flow with large data

In this article, we develop a branch of nonsingular solutions of a Picard multilevel stabilization of mixed finite volume method for the 2D/3D stationary Navier‐Stokes equations without relying on the unique solution condition. The method presented consists of capturing almost all information of initial problem (the nonlinear problems) on the coarsest mesh and then performs one Picard defect correction (the linear problems) on each subsequent mesh based on previous information thus only solving one large linear systems. What is more, the method presented can results in a better coefficient matrix in the model presented with small viscosity. Theoretical results show that the method presented is derived with the convergence rate of the same order as the corresponding finite volume method/finite element method solving the stationary Navier‐Stokes equations on a fine mesh. Therefore, the method presented is definitely more efficient than the standard finite volume method/finite element method. Finally, numerical experiments clearly show the efficiency of the method presented for solving the stationary Navier‐Stokes equations.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 34: 30–50, 2018     

The motion of a cluster of spherical particles in an ideal fluid

The motion of a cluster of an arbitrary number of spherical particles, attached to one another, in an ideal incompressible fluid is considered. Using the previously developed self-consistent field method, an expression is obtained for the virtual mass of the cluster in the form of an explicit function of the coordinates of all the particles. It is shown that, for special cases, the solution obtained is identical with the corresponding results known in the literature. For a statistically uniform particle distribution, using an averaging procedure over all the different possible configurations in the fluid inside a spherical volume, a simple analytic relation is obtained for the average value of the virtual mass of the spherical cluster as a function of its radius in the first approximation with respect to the volume fraction of the particles.     

A note on the suspension viscosity of porous spherical particles

The problem of the motion of a porous sphere undergoing a shearing motion is studied. An expression for the effective viscosity of the suspension is obtained and it is found that the viscosity, as compared to solid spheres, is smaller in the case of porous spheres.     

An iterative solution method for the Stokes problem with variable viscosity

An iterative method for efficient solution of the Stokes problem with a variable viscosity is considered. A preconditioner for the Shur complement is constructed taking into account the variable viscosity. The efficiency analysis is given. An application of the preconditioner for solving one problem of the mantle convection modeling is considered.     

The virtual mass of a sphere in a suspension of spherical particles

The problem of the virtual mass of a sphere, moving in an ideal incompressible fluid when there are other identical spherical particles of arbitrary mass present is considered. A solution is constructed for the velocity potential of the fluid in the form of the superposition of perturbation fields, introduced into the flow by each of the particles. The perturbation fields are obtained in the form of functional series, the coefficients of which are mutually consistent by a defined system of equations. An explicit expression is obtained for the hydrodynamic force acting on the sphere in the form of a function of the coordinates of all the particles. A simple analytical dependence of the mean value of the force and the virtual mass of the sphere on the particle-to-fluid density ratio in a first approximation of the volume fraction of the dispersed phase is obtained for a statistically uniform distribution of the dispersed particles in the suspension, using the procedure of averaging over their different possible configurations in space.     

Effect of homogeneous and heterogeneous chemical reactions on peristaltic transport of an MHD micropolar fluid with wall effects

In this paper, the influence of magnetic field on the dispersion of a solute in peristaltic flow of an incompressible micropolar fluid is studied as a model of fluid transport in the human intestinal system with wall properties. Long wavelength approximation, Taylor's limiting condition, and dynamic boundary conditions at the flexible walls are used to obtain the average effective dispersion coefficient in the presence of combined homogeneous and heterogeneous chemical reactions. The effects of various pertinent parameters on the effective dispersion coefficient are discussed. Average effective dispersion coefficient increases with amplitude ratio, which implies that dispersion is more in the presence of peristalsis. It also increases with the cross‐viscosity coefficient, heterogeneous chemical reaction rate, and wall parameters. Further, dispersion decreases with micropolar parameter, magnetic parameter, and homogeneous chemical reaction rates. Copyright © 2015 John Wiley & Sons, Ltd.