Formulation of the theory of turbulence in an incompressible fluid

The theory of turbulence in an incompressible fluid is formulated using methods similar to those of quantum field theory. A systematic perturbation theory is set up, and the terms in the perturbation series are shown to be in one to one correspondence with certain diagrams analogous to Feynman diagrams. From a study of the diagrams it is shown that the perturbation series can be rearranged and partially summed in such a way as to reduce the problem to the solution of three simultaneous integral equations for three functions, one of which is the second order velocity correlation function. The equations have the form of infinite power series integral equations, and the first few terms in the power series are derived from an analysis of the diagrams to sixth order. Truncation of the integral equations at the lowest nontrivial order yields Chandrasekhar's equation, and truncation at a higher order yields the equations discussed by Kraichnan.     

Integro-differential equation method in the hydrodynamics of an incompressible viscous fluid

An integro-differential equation method is proposed to describe the motion of an incompressible viscous fluid. The method uses an analogy between the hydrodynamic equations for an incompressible viscous fluid and the magnetostatic equations. An analysis is made of the flow of an incompressible viscous fluid round an object as a specific application. The solution automatically satisfies the boundary conditions at the surface of the object and at infinity. Zh. éksp. Teor. Fiz. 112, 1332–1339 (October 1997)     

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.     

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.     

A matrix-free,implicit, incompressible fractional-step algorithm for fluid–structure interaction applications

In this paper we detail a fast, fully-coupled, partitioned fluid–structure interaction (FSI) scheme. For the incompressible fluid, new fractional-step algorithms are proposed which make possible the fully implicit, but matrix-free, parallel solution of the entire coupled fluid–solid system. These algorithms include artificial compressibility pressure-poisson solution in conjunction with upwind velocity stabilisation, as well as simplified pressure stabilisation for improved computational efficiency. A dual-timestepping approach is proposed where a Jacobi method is employed for the momentum equations while the pressures are concurrently solved via a matrix-free preconditioned GMRES methodology. This enables efficient sub-iteration level coupling between the fluid and solid domains. Parallelisation is effected for distributed-memory systems. The accuracy and efficiency of the developed technology is evaluated by application to benchmark problems from the literature. The new schemes are shown to be efficient and robust, with the developed preconditioned GMRES solver furnishing speed-ups ranging between 50 and 80.     

Wave motion of incompressible turbulent fluid

Evolution of small disturbances in a fully developed incompressible turbulent flow is considered on the base of the transport equation for the single-point probability density function (PDF) of velocity fluctuations. It is shown that at high frequencies this equation is similar to the Vlasov equation for charged plasma in a self-consistent electromagnetic field having longitudinal wave solutions for turbulent stresses similar to Langmuir waves. It is found that the longitudinal waves of turbulent stresses have a constant phase velocity and can be damped, neutral, or growing waves, depending on the type of undisturbed probability density function of velocity fluctuations. The obtained result differs from the previously published solutions to this problem using the statistical moments closures according to which the wave disturbances should be neutral or damped. The possibilities of experimental observation of longitudinal waves of turbulent stresses are analyzed.     

The Beltrami spectrum for incompressible fluid flows

Recently V. Yakhot, S. Orszag, and their co-workers have suggested that turbulent flows in various regions of space organize into a coherent hierarchy of weakly interacting superimposed approximate Beltrami flows. A mathematical framework is developed here to study organized Beltrami hierarchies in a systematic fashion. This framework is applied to several important classes of examples with universal Beltrami hierarchies. An analysis of the persistence of such Beltrami hierarchies is also presented for general solutions of the Navier-Stokes equations.Research partially supported by grant NSF DMS-860-2031. Sloan research fellowship gratefully acknowledgedThis research was partially supported by grants NSF DMS 86-11110 and DARPA — ONR N00014-86-K-0759     

Translation drag coefficient of a self-similar assembly of spheres immersed in an incompressible fluid.

On the basis of deterministic fractals and the Rotne-Prager hydrodynamic interaction tensor, we confirm the asymptotic as well as the finite size scaling of the friction coefficient lambda of a self-similar structure. The fractal assembly is made of N spheres with its dimension varying from D < 1 to D = 3. The number of spheres can be as high as N approximately O(10(4)). The asymptotic scaling behavior of the friction coefficient per sphere is lambda approximately N(1/D-1) for D > 1, lambda approximately (lnN)(-1) for D = 1, and lambda approximately N(0) for D < 1. The crossover behavior indicates that while in the regime of D > 1 the hydrodynamic screening effect grows with the size, for D<1 it is limited in a finite range, which decays with decreasing D.     

Interaction of Alfven waves in the linearized system of magnetohydrodynamics for an incompressible ideal fluid

Interaction of Alfven waves in the vicinity of multiple characteristic roots is considered. The asymptotic behavior for solutions of the Cauchy problem for the linearized equations of magnetohydrodynamics with highly oscillating initial data is obtained.     

The Hamiltonian description of incompressible fluid ellipsoids

We construct the noncanonical Poisson bracket associated with the phase space of first order moments of the velocity field and quadratic moments of the density of a fluid with a free-boundary, constrained by the condition of incompressibility. Two methods are used to obtain the bracket, both based on Dirac’s procedure for incorporating constraints. First, the Poisson bracket of moments of the unconstrained Euler equations is used to construct a Dirac bracket, with Casimir invariants corresponding to volume preservation and incompressibility. Second, the Dirac procedure is applied directly to the continuum, noncanonical Poisson bracket that describes the compressible Euler equations, and the moment reduction is applied to this bracket. When the Hamiltonian can be expressed exactly in terms of these moments, a closure is achieved and the resulting finite-dimensional Hamiltonian system provides exact solutions of Euler’s equations. This is shown to be the case for the classical, incompressible Riemann ellipsoids, which have velocities that vary linearly with position and have constant density within an ellipsoidal boundary. The incompressible, noncanonical Poisson bracket differs from its counterpart for the compressible case in that it is not of Lie-Poisson form.     

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.