Projection conditions on the vorticity in viscous incompressible flows

The problem of establishing appropriate conditions for the vorticity transport equation is considered. It is shown that, in viscous incompressible flows, the boundary conditions on the velocity imply conditions of an integral type on the vorticity. These conditions determine a projection of the vorticity field on the linear manifold of the harmonic vector fields. Some computational consequences of the above result in two-dimensional calculations by means of the nonprimitive variables, stream function and vorticity, are examined. As an example of the application of the discrete analogue of the projection conditions, numerical solutions of the driven cavity problem are reported.     

A vorticity formulation for computational fluid dynamic and aeroelastic analyses of viscous flows

The paper deals with a novel computational formulation for the analysis of viscous flows past a solid body. The formulation is based upon a convenient decomposition of the flow field into potential and rotational velocity contributions, which has the distinguishing feature that the rotational velocity vanishes in much of, if not all, the region in which the vorticity is negligible. Contrary to related formulations implemented by the authors in the past, in the proposed approach, discontinuities of the potential and rotational velocity fields across a prescribed surface emanating from the trailing edge (such as the wake mid-surface) are eliminated, thereby facilitating numerical implementations. However, the main novelty is related to the application of the boundary condition: first, the expression for the velocity used for the condition on the body boundary is consistent with that for the velocity in the field; also—contrary to related formulations used by the authors in the past—in the proposed approach, the condition on the body boundary does not require the evaluation of the total vorticity (inside and outside the computational domain). The proposed algorithm, valid for three-dimensional compressible flows, is validated—as a first step—for the case of two-dimensional incompressible flows. Specifically, numerical results are presented for the aerodynamic analysis of two-dimensional incompressible viscous flows past a circular cylinder and past a Joukowski airfoil. In order to verify the desirable absence of artificial damping, we present also results pertaining to the flutter (i.e., dynamic aeroelastic) analysis of a spring-mounted circular cylinder in a viscous flow, free to move in a direction orthogonal to the unperturbed flow. In both cases (aerodynamics and aeroelasticity), the results are in good agreement with existing literature data.     

Multivariant finite elements for three-dimensional simulation of viscous incompressible flows

Within multivariant elements, which have restricted degrees of freedom at some nodes, different velocity components have different variations. Shape functions for the multivariant elements Q P_o and R P_o are developed. With such shape functions the value of a velocity component within a multivariant element is shown to depend upon all the independent components of velocity at the nodes of the element. The use of the Q₁ P₀ element to simulate flows with discontinuous boundary conditions generated disturbance throughout the flow domain, giving erroneous pressure and velocity distributions. The Q P_o element restricted the disturbance due to such discontinuities to a small region near the singular points, whereas the P P_o element completely eliminated the fluctuations. Flows with discontinuous boundary conditions were simulated with reasonable accuracy by partially relaxing the no-slip condition on the Q₁ P_o elements near the singular points.     

Parallel overlapping scheme for viscous incompressible flows

A 3D parallel overlapping scheme for viscous incompressible flow problems is presented that combines the finite element method, which is best suited for analysing flow in any arbitrarily shaped flow geometry, with the finite difference method, which is advantageous in terms of both computing time and computer storage. A modified ABMAC method is used as the solution algorithm, to which a sophisticated time integration scheme proposed by the present authors has been applied. Parallelization is based on the domain decomposition method. The RGB (recursive graph bisection) algorithm is used for the decomposition of the FEM mesh and simple slice decomposition is used for the FDM mesh. Some estimates of the parallel performance of FEM, FDM and overlapping computations are presented. © 1997 John Wiley & Sons, Ltd.     

A hybrid vertex-centered finite volume/element method for viscous incompressible flows on non-staggered unstructured meshes

This paper proposes a hybrid vertex-centered finite volume/finite element method for solution of the two dimensional (2D) incompressible Navier-Stokes equations on unstructured grids.An incremental pressure fractional step method is adopted to handle the velocity-pressure coupling.The velocity and the pressure are collocated at the node of the vertex-centered control volume which is formed by joining the centroid of cells sharing the common vertex.For the temporal integration of the momentum equations,an implicit second-order scheme is utilized to enhance the computational stability and eliminate the time step limit due to the diffusion term.The momentum equations are discretized by the vertex-centered finite volume method (FVM) and the pressure Poisson equation is solved by the Galerkin finite element method (FEM).The momentum interpolation is used to damp out the spurious pressure wiggles.The test case with analytical solutions demonstrates second-order accuracy of the current hybrid scheme in time and space for both velocity and pressure.The classic test cases,the lid-driven cavity flow,the skew cavity flow and the backward-facing step flow,show that numerical results are in good agreement with the published benchmark solutions.     

Unsteady flows of an inhomogeneous incompressible viscous fluid

This paper deals with a theoretical analysis of the transfer of reactive impurities by open and filtration flows of an incompressible viscous fluid. The first section of the paper studies the model of an inhomogeneous incompressible viscous fluid, which is widely used in meteorology and oceanology, with additional allowance for the drag of the magnetic field or porous medium. Another object of research in this paper is the model of filtration of an inhomogeneous incompressible fluid in porous media proposed by V. N. Monakhov (1977) (Section 2). In both models, hydrodynamic flows determine the motion of the mixture as a whole and the temperature and concentration distributions of the components of an inhomogeneous fluid are described by a common nonlinear system of equations of diffusive heat and mass transfer.Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 2, pp. 44–51, March–April, 2005.     

Segmented multigrid domain decomposition procedure for incompressible viscous flows

The application of grid stretching or grid adaptation is generally required in order to optimize the distribution of nodal points for fluid-dynamic simulation. This is necessitated by the presence of disjoint high gradient zones, that represent boundary or free shear layers, reversed flow or vortical flow regions, triple deck structures, etc. A domain decomposition method can be used in conjunction with an adaptive multigrid algorithm to provide an effective methodology for the development of optimal grids. In the present study, the Navier-Stokes (NS) equations are approximated with a reduced Navier-Stokes (RNS) system, that represents the lowest-order terms in an asymptotic Re expansion. This system allows for simplified boundary conditions, more generality in the location of the outflow boundary, and ensures mass conservation in all subdomain grid interfaces, as well as at the outflow boundary. The higher-order (NS) diffusion terms are included through a deferred corrector, in selected subdomains, when necessary. Adaptivity in the direction of refinement is achieved by grid splitting or domain decomposition in each level of the multigrid procedure. Normalized truncation error estimates of key derivatives are used to determine the boundaries of these subdomains. The refinement is optimized in two co-ordinate directions independently. Multidirectional adaptivity eliminates the need for grid stretching so that uniform grids are specified in each subdomain. The overall grid consists of multiple domains with different meshes and is, therefore, heavily graded. Results and computational efficiency are discussed for the laminar flow over a finite length plate and for the laminar internal flow in a backward-facing step channel.     

Numerical solution of a viscous incompressible flow problem through an orifice

A steady flow problem of a viscous, incompressible fluid through an orifice is widely applicable to many physical phenomena and has been studied previously by many researchers. A problem of such type has been solved by applying LAD method given by Roache [1]. The resulting system of linear equations is solved by Hockney's method [2].     

A boundary element method for steady viscous fluid flow using penalty function formulation

A new boundary element method is presented for steady incompressible flow at moderate and high Reynolds numbers. The whole domain is discretized into a number of eight-noded cells, for each of which the governing boundary integral equation is formulated exclusively in terms of velocities and tractions. The kernels used in this paper are the fundamental solutions of the linearized Navier–Stokes equations with artificial compressibility. Significant attention is given to the numerical evaluation of the integrals over quadratic boundary elements as well as over quadratic quadrilateral volume cells in order to ensure a high accuracy level at high Reynolds numbers. As an illustration, square driven cavity flows are considered for Reynolds numbers up to 1000. Numerical results demonstrate both the high convergence rate, even when using simple (direct) iterations, and the appropriate level of accuracy of the proposed method. Although the method yields a high level of accuracy in the primary vortex region, the secondary vortices are not properly resolved. © 1997 John Wiley & Sons, Ltd.     

Conjugation of the channel and filtration flows of a viscous incompressible fluid

Novosibirsk. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, No. 1, pp. 95–99, January–February, 1995.     

Non-singular direct formulation of boundary integral equations for potential flows

This paper presents a general direct integral formulation for potential flows. The singularities of Green's functions are desingularized theoretically, using a subtracting and adding back technique, so that Gaussian quadrature or any other numerical integration methods can be applied directly to evaluate all the integrals without any difficulty. When high-order quadrature formulas are applied globally, the number of unknowns can be reduced. Interpolation functions are not necessary for unknown variables in the present paper. Therefore, the present method is much simpler and more efficient than the conventional one. Several numerical examples are calculated and compared satisfactorily with analytical solutions or published results. © 1998 John Wiley & Sons, Ltd.     

A subdomain boundary element method for high‐Reynolds laminar flow using stream function‐vorticity formulation

The paper presents a new formulation of the integral boundary element method (BEM) using subdomain technique. A continuous approximation of the function and the function derivative in the direction normal to the boundary element (further ‘normal flux’) is introduced for solving the general form of a parabolic diffusion‐convective equation. Double nodes for normal flux approximation are used. The gradient continuity is required at the interior subdomain corners where compatibility and equilibrium interface conditions are prescribed. The obtained system matrix with more equations than unknowns is solved using the fast iterative linear least squares based solver. The robustness and stability of the developed formulation is shown on the cases of a backward‐facing step flow and a square‐driven cavity flow up to the Reynolds number value 50 000. Copyright © 2004 John Wiley & Sons, Ltd.     

Mixed boundary elements for laminar flows

This paper presents a mixed boundary element formulation of the boundary domain integral method (BDIM) for solving diffusion–convective transport problems. The basic idea of mixed elements is the use of a continuous interpolation polynomial for conservative field function approximation and a discontinuous interpolation polynomial for its normal derivative along the boundary element. In this way, the advantages of continuous field function approximation are retained and its conservation is preserved while the normal flux values are approximated by interpolation nodal points with a uniquely defined normal direction. Due to the use of mixed boundary elements, the final discretized matrix system is overdetermined and a special solver based on the least squares method is applied. Driven cavity, natural and forced convection in a closed cavity are studied. Driven cavity results at Re=100, 400 and 1000 agree better with the benchmark solution than Finite Element Method or Finite Volume Method results for the same grid density with 21×21 degrees of freedom. The average Nusselt number values for natural convection 10³≤Ra≤10⁶ agree better than 0.1% with benchmark solutions for maximal calculated grid densities 61×61 degrees of freedom. Copyright © 1999 John Wiley & Sons, Ltd.     

Fast boundary‐domain integral algorithm for the computation of incompressible fluid flow problems

The paper deals with the numerical solution of fluid dynamics using the boundary‐domain integral method (BDIM). A velocity–vorticity formulation of the Navier–Stokes equations is adopted, where the kinematic equation is written in its parabolic form. Computational aspects of the numerical simulation of two‐dimensional flows is described in detail. In order to lower the computational cost, the subdomain technique is applied. A preconditioned Krylov subspace method (PKSM) is used for the solution of systems of linear equations. Level‐based fill‐in incomplete lower upper decomposition (ILU) preconditioners are developed and their performance is examined. Scaling of stopping criteria is applied to minimize the number of iterations for the PKSM. The effectiveness of the proposed method is tested on several benchmark test problems. Copyright © 1999 John Wiley & Sons, Ltd.     

A new finite element formulation for both compressible and nearly incompressible fluid dynamics

A new finite element formulation designed for both compressible and nearly incompressible viscous flows is presented. The formulation combines conservative and non‐conservative dependent variables, namely, the mass–velocity (density * velocity), internal energy and pressure. The central feature of the method is the derivation of a discretized equation for pressure, where pressure contributions arising from the mass, momentum and energy balances are taken implicitly in the time discretization. The method is applied to the analysis of laminar flows governed by the Navier–Stokes equations in both compressible and nearly incompressible regimes. Numerical examples, covering a wide range of Mach number, demonstrate the robustness and versatility of the new method. Copyright © 2000 John Wiley & Sons, Ltd.     

A finite amplitude analysis of tunnel deformation by the finite element method with highly viscous fluid

A finite element method for highly viscous fluid is used to calculate the velocity and stress fields in the surrounding soft rock of a tunnel. In order to fit the calculated values with the measured displacement of tunnel wall, we inverted the boundary forces and the mechanical parameters of the surrounding rocks.     

Finite element analysis of the axially symmetric motion of an incompressible viscous fluid in a spherical annulus

This paper presents results obtained by employing a modified Galerkin finite element method to analyse the steady state flow of a fluid contained between two concentric, rotating spheres. The spheres are assumed to be rigid and the cavity region between the spheres is filled with an incompressible, viscous, Newtonian fluid. The inner sphere is constrained to rotate about a vertical axis with a prescribed angular velocity, while the outer sphere is fixed. Results for the circumferential function Ω, streamfunction ψ, vorticity function ζ and inner boundary torque T₁ are presented for Reynolds numbers Re ? 2000 and radius ratios 0.1 ? α ? 0.9. The method proved effective for obtaining results for a wide range of radius ratios (0.1 ? α ? 0.9) and Reynolds numbers (0 ? Re ? 2000). Previous investigators who employed the finite difference method experienced difficulties in obtaining results for cases with radius ratios α ? 0.2, except for small Reynolds numbers (Re ? 100). Results for Ω, Ψ, ζ and T₁ obtained in this study for radius ratios 0.8 ≤ α ≤ 0.9 verified the development of Taylor vortices reported by other investigators. The research indicates that the method may be useful for analysing other non-linear fluid flow problems.     

Gauge principle and variational formulation for flows of an ideal fluid

A gauge principle is applied to mass flows of an ideal compressible fluid subject to Galilei transformation. A free-field Lagrangian defined at the outset is invariant with respeet to global SO(3) gauge transformations as well as Galilei transformations. The action principle leads to the equation of potential flows under constraint of a continuity equation. However, the irrotational flow is not invariant with respect to local SO(3) gauge transformations. According to the gauge principle, a gauge-covariant derivative is defined by introducing a new gauge field. Galilei invariance of the derivative requires the gauge field to coincide with the vorticity, i.e. the curl of the velocity field. A full gauge-covariant variational formulation is proposed on the basis of the Hamilton‘‘s principle and an assoicated Lagrangian. By means of an isentropic material variation taking into account individual particle motion, the Euler‘‘s equation of motion is derived for isentropic flows by using the covariant derivative. Noether‘‘s law associated with global SO(3) gauge invariance leads to the conservation of total angular momentum. In addition, the Lagrangian has a local symmetry of particle permutation which results in local conservation law equivalent to the vorticity equation.     

A three-dimensional fictitious domain method for incompressible fluid flow problems

A new Galerkin finite element method for the solution of the Navier–Stokes equations in enclosures containing internal parts which may be moving is presented. Dubbed the virtual finite element method, it is based upon optimization techniques and belongs to the class of fictitious domain methods. Only one volumetric mesh representing the enclosure without its internal parts needs to be generated. These are rather discretized using control points on which kinematic constraints are enforced and introduced into the mathematical formulation by means of Lagrange multipliers. Consequently, the meshing of the computational domain is much easier than with classical finite element approaches. First, the methodology will be presented in detail. It will then be validated in the case of the two-dimensional Couette cylinder problem for which an analytical solution is available. Finally, the three-dimensional fluid flow inside a mechanically agitated vessel will be investigated. The accuracy of the numerical results will be assessed through a comparison with experimental data and results obtained with a standard finite element method. © 1997 John Wiley & Sons, Ltd.     

Variational formulation of three-dimensional viscous free-surface flows: Natural imposition of surface tension boundary conditions

We present a new surface-intrinsic linear form for the treatment of normal and tangential surface tension boundary conditions in C⁰-geometry variational discretizations of viscous incompressible free-surface flows in three space dimensions. The new approach is illustrated by a finite (spectral) element unsteady Navier-Stokes analysis of the stability of a falling liquid film.