Nested Cartesian grid method in incompressible viscous fluid flow

In this work, the local grid refinement procedure is focused by using a nested Cartesian grid formulation. The method is developed for simulating unsteady viscous incompressible flows with complex immersed boundaries. A finite-volume formulation based on globally second-order accurate central-difference schemes is adopted here in conjunction with a two-step fractional-step procedure. The key aspects that needed to be considered in developing such a nested grid solver are proper imposition of interface conditions on the nested-block boundaries, and accurate discretization of the governing equations in cells that are with block-interface as a control-surface. The interpolation procedure adopted in the study allows systematic development of a discretization scheme that preserves global second-order spatial accuracy of the underlying solver, and as a result high efficiency/accuracy nested grid discretization method is developed. Herein the proposed nested grid method has been widely tested through effective simulation of four different classes of unsteady incompressible viscous flows, thereby demonstrating its performance in the solution of various complex flow–structure interactions. The numerical examples include a lid-driven cavity flow and Pearson vortex problems, flow past a circular cylinder symmetrically installed in a channel, flow past an elliptic cylinder at an angle of attack, and flow past two tandem circular cylinders of unequal diameters. For the numerical simulations of flows past bluff bodies an immersed boundary (IB) method has been implemented in which the solid object is represented by a distributed body force in the Navier–Stokes equations. The main advantages of the implemented immersed boundary method are that the simulations could be performed on a regular Cartesian grid and applied to multiple nested-block (Cartesian) structured grids without any difficulty. Through the numerical experiments the strength of the solver in effectively/accurately simulating various complex flows past different forms of immersed boundaries is extensively demonstrated, in which the nested Cartesian grid method was suitably combined together with the fractional-step algorithm to speed up the solution procedure.     

A Lagrangian and canonical formalism in the hydrodynamics of an incompressible fluid

The geometrical formulation of the variational principle in the hydrodynamics of an incompressible fluid is given. The Lagrange-Euler equations in the language of differential forms are formulated. The canonical formalism in the hydrodynamics in the general framework of multi-phase spaces ([4]) is considered and the construction of the phase space for hydrodynamics is given.     

Heat transfer of an oscillating cylinder in a viscous incompressible fluid flow

Unsteady flow and heat transfer of the heated cylinder with transverse and longitudinal oscillations in unbounded flow of viscous incompressible heat conducting fluid are considered. Influence of amplitude, frequency and direction of harmonic oscillations of the cylinder on the near-wake structure and heat transfer are investigated.     

Stability of axisymmetric swirl flows of viscous incompressible fluid

A new method of solution to the problem of stability of the swirl flow of viscous incompressible fluid is developed. The method based on expansion of the required function into power series of radial coordinate allows an avoidance of difficulties related to numerical integration of the system of differential equations with a singular point. Stability of the Poiseuille flow in a rotating pipe is considered as an example.     

Singularities motion equations in 2-dimensional ideal hydrodynamics of incompressible fluid

In this Letter, we have obtained motion equations for a wide class of one-dimensional singularities in 2D ideal hydrodynamics. The simplest of them, are well known as point vortices. More complicated singularities correspond to vorticity point dipoles. It has been proved that point multipoles of a higher order (quadrupoles and more) are not the exact solutions of two-dimensional ideal hydrodynamics. The motion equations for a system of interacting point vortices and point dipoles have been obtained. It is shown that these equations are Hamiltonian ones and have three motion integrals in involution. It means the complete integrability of two-particle system, which has a point vortex and a point dipole.     

基于不可压缩光滑粒子动力学的黏性液滴变形过程仿真

应用基于投影算法的不可压缩光滑粒子动力学(incompressible smoothed particle hydrodynamics, ISPH)法对黏性液滴变形过程进行了数值仿真. 对于张力失稳导致的粒子非物理簇集问题, 采用粒子移位技术加以解决. 为了验证本文ISPH 算法的精度和稳定性, 分别模拟了圆形黏性液滴的拉伸变形过程以及方形液滴的旋转变形过程, 得到了不同时刻液滴内部的压力变化特征, 准确地捕捉了液滴自由面演化过程, 并将数值计算结果与文献中的解析解进行了比较.分析结果表明, 基于投影算法的不可压缩光滑粒子动力学方法结合粒子移位技术, 能够有效地模拟黏性液滴变形过程, 可以得到精确和稳定的结果. 关键词： 不可压缩光滑粒子动力学 黏性液滴 自由面流动 数值仿真     

A numerical method for incompressible viscous flow simulation

We describe a numerical scheme for computing time-dependent solutions of the incompressible Navier-Stokes equations in the primitive variable formulation. This scheme uses finite elements for the space discretization and operator splitting techniques for the time discretization. The resulting discrete equations are solved using specialized nonlinear optimization algorithms that are computationally efficient and have modest storage requirements. The basic numerical kernel is the preconditioned conjugate gradient method for symmetric, positive-definite, sparse matrix systems, which can be efficiently implemented on the architectures of vector and parallel processing supercomputers.     

On the motion of a sphere with arbitrary slip in a viscous incompressible fluid

An expression for the force on a sphere moving with a time-dependent velocity through an incompressible fluid in nonstationary, nonhomogeneous flow is obtained for the case of arbitrary slip on the surface of the sphere.     

Thermodynamic theory of incompressible hydrodynamics

The grand potential for open systems describes thermodynamics of fluid flows at low Mach numbers. A new system of reduced equations for the grand potential and the fluid momentum is derived from the compressible Navier-Stokes equations. The incompressible Navier-Stokes equations are the quasistationary solution to the new system. It is argued that the grand canonical ensemble is the unifying concept for the derivation of models and numerical methods for incompressible fluids, illustrated here with a simulation of a minimal Boltzmann model in a microflow setup.     

An analysis of the incompressible viscous flows problem by meshless method

Generally the incompressible viscous flow problem is described by the Navier--Stokes equation. Based on the weighted residual method the discrete formulation of element-free Galerkin is inferred in this paper. By the step-by-step computation in the field of time, and adopting the least-square estimation of the-same-order shift, this paper has calculated both velocity and pressure from the decoupling independent equations. Each time fraction Newton--Raphson iterative method is applied for the velocity and pressure. Finally, this paper puts the method into practice of the shear-drive cavity flow, verifying the validity, high accuracy and stability.     

New equation for the polarization fourier components of the small-scale velocity of an incompressible fluid in anisotropic turbulence

A new equation for the small-scale polarization Fourier components of the incompressible fluid velocity in the case of anisotropic turbulence is suggested. The principal invariant of the strain rate tensor for the large-scale velocity is found. This invariant is of most significance for the subgrid simulation of fully developed turbulence.     

Asymmetric self-similar flows of a viscous incompressible fluid along a right-angle corner

Symmetric and asymmetric self-similar flows of a viscous incompressible fluid along a semi-infinite right-angle dihedral corner with a preset streamwise pressure gradient have been considered. Equations describing such flows in the framework of boundary layer approximation have been derived. The asymptotic behavior of solutions of the derived equations far from the corner edge has been theoretically investigated. A new method of computation of these solutions has been developed. Solutions for two types of asymptotic behavior have been obtained.     

A solution-adaptive lattice Boltzmann method for two-dimensional incompressible viscous flows

A stencil adaptive lattice Boltzmann method (LBM) is developed in this paper. It incorporates the stencil adaptive algorithm developed by Ding and Shu [26] for the solution of Navier–Stokes (N–S) equations into the LBM calculation. Based on the uniform mesh, the stencil adaptive algorithm refines the mesh by two types of 5-points symmetric stencils, which are used in an alternating sequence for increased refinement levels. The two types of symmetric stencils can be easily combined to form a 9-points symmetric structure. Using the one-dimensional second-order interpolation recently developed by Wu and Shu [27] along the straight line and the D2Q9 model, the adaptive LBM calculation can be effectively carried out. Note that the interpolation coefficients are only related to the lattice velocity and stencil size. Hence, the simplicity of LBM is not broken down and the accuracy is maintained. Due to the use of adaptive technique, much less mesh points are required in the simulation as compared to the standard LBM. As a consequence, the computational efficiency is greatly enhanced. The numerical simulation of two dimensional lid-driven cavity flows is carried out. Accurate results and improved efficiency are reached. In addition, the steady and unsteady flows over a circular cylinder are simulated to demonstrate the capability of proposed method for handling problems with curved boundaries. The obtained results compare well with data in the literature.     

耦合不可压流场输运方程的格子Boltzmann方法研究

基于格子Boltzmann方法,提出了求解耦合不可压缩流场输运方程的一种改进数值方法. 该方法使用格子Boltzmann方法求解流场方程,并根据流场格子模型的密度分布函数构建了输运方程的二阶离散格式. 通过二维平板通道流场输运系统验证了该方法的有效性.数值结果表明,该方法可以有效地减少计算过程中出现的非物理耗散, 并克服了传统模型所需巨大存储量的缺点.