A first Integral of Navier­Stokes equations for two-dimensional flows

As wellknown, Bernoulli's equation is obtained as the first integral of Euler's equations in the absence of vorticity. Even in case of non-vanishing vorticity, a first integral from Euler's equations is obtained by using the so called Clebsch transformation [1] for inviscid flows. In contrast to this, a generalisation of this procedure towards viscous flows has not been established so far. In the present paper a first integral of Navier-Stokes equations is constructed in the case of two-dimensional flow by making use of an alternative representation of the fields in terms of complex coordinates and introducing a potential representation for the pressure. The associated boundary conditions are also considered. The first integral is a suitable tool for the development of new analytical methods and numerical codes in fluid dynamics. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The Bernoulli Integral for a Certain Class of Non‐Stationary Viscous Vortical Flows of Incompressible Fluid

It has been shown in our paper [1] that there is a wide class of 3D motions of incompressible viscous fluid which can be described by one scalar function dabbed the quasi‐potential. This class of fluid flows is characterized by three‐component velocity field having two‐component vorticity field; both these fields can depend of all three spatial variables and time, in general. Governing equations for the quasi‐potential have been derived and simple illustrative example of 3D flow has been presented. Here, we derive the Bernoulli integral for that class of flows and compare it against the known Bernoulli integrals for the potential flows or 2D stationary vortical flows of inviscid fluid. We show that the Bernoulli integral for this class of fluid motion possesses unusual features: it is valid for the vortical nonstationary motions of a viscous incompressible fluid. We present a new very nontrivial analytical example of 3D flow with two‐component vorticity which hardly can be obtained by any of known methods. In the last section, we suggest a generalization of the developed concept which allows one to describe a certain class of 3D flows with the 3D vorticity.     

A High-Order Kernel-Free Boundary Integral Method for Incompressible Flow Equations in Two Space Dimensions

This paper presents a fourth-order kernel-free boundary integral method for the time-dependent, incompressible Stokes and Navier-Stokes equations defined on irregular bounded domains. By the stream function-vorticity formulation, the incompressible flow equations are interpreted as vorticity evolution equations. Time discretization methods for the evolution equations lead to a modified Helmholtz equation for the vorticity, or alternatively, a modified biharmonic equation for the stream function with two clamped boundary conditions. The resulting fourth-order elliptic boundary value problem is solved by a fourth-order kernel-free boundary integral method, with which integrals in the reformulated boundary integral equation are evaluated by solving corresponding equivalent interface problems, regardless of the exact expression of the involved Green's function. To solve the unsteady Stokes equations, a four-stage composite backward differential formula of the same order accuracy is employed for time integration. For the Navier-Stokes equations, a three-stage third-order semi-implicit Runge-Kutta method is utilized to guarantee the global numerical solution has at least third-order convergence rate. Numerical results for the unsteady Stokes equations and the Navier-Stokes equations are presented to validate efficiency and accuracy of the proposed method.     

Effective velocity in compressible Navier-Stokes equations with third-order derivatives

A formulation of certain barotropic compressible Navier-Stokes equations with third-order derivatives as a viscous Euler system is proposed by using an effective velocity variable. The equations model, for instance, viscous Korteweg or quantum Navier-Stokes flows. The formulation in the new variable allows for the derivation of an entropy identity, which is known as the BD (Bresch-Desjardins) entropy equation. As a consequence of this estimate, a new global-in-time existence result for the one-dimensional quantum Navier-Stokes equations with strictly positive particle densities is proved.     

The Solution of Two-Dimensional Oseen Flow Problems Using Integral Conditions

A general method of solving Oseen's linearized equations fortwo-dimensional steady flow of a viscous incompressible fluidpast a cylinder in an unbounded field is developed. The analysisis developed in terms of the scalar vorticity and stream functionand it is shown that the vorticity for Oseen flow problems canbe obtained separately from the stream function. The determinationof the vorticity can be effected using conditions of an integralcharacter deduced from the no-slip condition at the cylindersurface together with the conditions at large distances. Theindependent determination of the vorticity seems to be a newstep in Oseen theory. The method enables one to obtain manyproperties of the flow in terms ofthe Reynolds number by usingonly the vorticity without the necessity of finding the streamfunction. The use of integral conditions makes the detailedcalculations straightforward, systematic, and elementary. Themethod is tested by applying it to the case of uniform flowpast an elliptic cylinder at an arbitrary angle of incidenceand also to cases of symmetrical and asymmetrical flows pastcircular cylinders. The leading approximation for small Reynoldsnumber is obtained where possible. In the case of flow pasta rotating cylinder, the only possible solution is the Oseensolution for the nonrotating case with the addition of a potentialvortex.     

On a potential-velocity formulation of incompressible Navier-Stokes equations

Computational fluid dynamics has emerged as an essential investigative tool in nearly every field of technology. Despite a well-developed mathematical theory and the existence of readily available commercial software codes, computing solutions to the governing equations of fluid motion remains challenging, especially due to the non-linearity involved. Additionally, in the case of free surface film flows the dynamic boundary condition at the free surface complicates the mathematical treatment notably. Recently, by introduction of an auxiliary potential field, a first integral of the two-dimensional incompressible Navier-Stokes equations has been constructed leading to a set of equations, the differential order of which is lower than that of the original equations [1]. A useful application to free surface simulation was found in [2]. Moreover the new formulation is naturally extendible to three dimensions via tensor calculus, involving a non-unique symmetric tensor potential. The corresponding degrees of freedom can be used in order to achieve a numerically convenient representation. Finally an efficient staggered-grid finite difference scheme is applied to a Stokes flow problem in a 3D lid-driven cavity to demonstrate the capabilities of the new approach. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the convergence of multi-grid methods with transforming smoothers

Summary In the present paper we give a convergence theory for multi-grid methods with transforming smoothers as introduced in [31] applied to a general system of partial differential equations. The theory follows Hackbusch's approach for scalar pde and allows a convergence proof for some well-known multi-grid methods for Stokes- and Navier-Stokes equations as DGS by Brandt-Dinar, [5], TILU from [31] and the SIMPLE-methods by Patankar-Spalding, [23].This work was supported in part by Deutsche Forschungsgemeinschaft     

High-order ADI method for stream-function vorticity equations

Many recent works have demonstrated the efficiency of high-order compact (HOC) difference schemes on the stream-function and vorticity formulation of 2-D incompressible Navier-Stokes equations. HOC discretizations induce cross spatial derivatives which are treated explicitly in most ADI schemes. Recently, Karaa and Zhang proposed a fourth-order ADImethod for solving convection-diffusion problems efficiently. In this work, we extend this method to the solution of incompressible Navier-Stokes equations. The driven flow in a square cavity is used as a model problem and numerical results are compared with other results. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Meshless solution of two-dimensional incompressible flow problems using the radial basis integral equation method

The two-dimensional incompressible fluid flow problems governed by the velocity–vorticity formulation of the Navier–Stokes equations were solved using the radial basis integral (RBIE) equation method. The RBIE is a meshless method based on the multi-domain boundary element method with overlapping subdomains. It solves at each node for the potential and its spatial derivatives. This feature of the RBIE is advantageous in solving the velocity–vorticity formulation of the Navier–Stokes equations since the calculated velocity gradients can be used to compute the vorticity that is prescribed as a boundary condition to the vorticity transport equation. The accuracy of the numerical solution was examined by solving the test problem with known analytical solution. Two benchmark problems, i.e. the lid driven cavity flow and the thermally driven cavity flow were also solved. The numerical results obtained using the RBIE showed very good agreement with the benchmark solutions.     

Adapted kussmaul formulation of acoustic scattering problems: an overview

Problems of exterior acoustic scattering may be conveniently formulated by means of boundary integral equations. The problem seeks to find a wave function which gives velocity potential profile, pressure density profile, etc. of the acoustic wave at points in space. At the background of the formulations are two theories viz. (Helmholtz) Potential theory and the Green's representation formula. Potential theory gives rise to the so-called indirect formulation and the Green's representation formula to the direct formulations. Classical boundary integral formulations fail at the eigenfrequencies of the interior domain. That is, if a solution is sought of the exterior problem by first solving a homogeneous boundary integral equation, one is inevitably led to the conclusion that these homogeneous boundary equations have nontrivial solutions at certain wave-numbers which are the eigenvalues of the corresponding interior problem. At lower wave-numbers, these eigenfrequencies are thinly distributed but the higher the wave-number, the denser it becomes. This is a well-known drawback for both time-harmonic acoustics and elastodynamics. This is not a physical difficulty but arises entirely as a result of a deficiency in the integral equation is representation. Why then use It? The use has many advantages notably in that the meshing region is reduced from the infinite domain exterior to the body to its finite surface. This created the need for some robust formulations. A proof of the Kussmaul [1] formulation is presented. The formulation has a hypersingular kernel in the integral operator, which creates a havoc in computation (e.g., ill conditioning). The hyper-singularity can be avoided [2], as a result a new formulation is proposed. This paper presents a broad overview of the Adapted Kussmaul Formulation (AKF).     

Mean-Square Approximation of Navier-Stokes Equations with Additive Noise in Vorticity-Velocity Formulation

We consider a time discretization of incompressible Navier-Stokes equations with spatial periodic boundary conditions and additive noise in the vorticity-velocity formulation. The approximation is based on freezing the velocity on time subintervals resulting in a linear stochastic parabolic equation for vorticity. At each time step, the velocity is expressed via vorticity using a formula corresponding to the Biot-Savart-type law. We prove the first mean-square convergence order of the vorticity approximation.     

The method of lines for the computation of incompressible viscous flows

A common approach to finding numerical solutions of the time-dependent incompressible Navier-Stokes equations is considered within the method of lines framework [1]. For the determination of viscous incompressible flows the stream-function formulation for the fourth-order equation [2, 3], an artificial compressibility method [4], and a modified velocity correction method [5] are designed. Some improved and extended results of numerical simulations obtained by the author in the previous works are presented. Test calculations have been done for various flows inside square, triangular and semicircular cavities with one moving wall, the backward-facing step, double bent channels and for the flow around an aerofoil at large angle of attack. An alternative and practical methodology for resolving the Navier-Stokes equations in arbitrarily complex geometries using Cartesian meshes is proposed. Some of complex geometrical configurations can be decomposed into a set of subdomains. The simplest approach for specifying boundary conditions near curved or irregular boundaries is to transfer all the variables from the boundaries to the nearest grid knots. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On Some Model Equations of Euler and Navier-Stokes Equations

The author proposes a two-dimensional generalization of Constantin-Lax-Majda model. Some results about singular solutions are given. This model might be the first step toward the singular solutions of the Euler equations. Along the same line (vorticity formulation), the author presents some further model equations. He possibly models various aspects of difficulties related with the singular solutions of the Euler and Navier-Stokes equations. Some discussions on the possible connection between turbulence and the singular solutions of the Navier-Stokes equations are made.     

The Fourier Pseudospectral Method for Two-Dimensional Vorticity Equations

In this paper we develop a Fourier pseudospectral method forsolving two-dimensional vorticity equations. We prove the generalizedstability of the schemes and give convergence estimations dependingon the smoothness of the solution of the vorticity equations. Spectral methods have been applied widely to the partial differentialequations of fluid dynamics [4–11]. Guo Ben-yu proposeda technique to estimate strictly the error of the spectral schemesfor the K.D.V.-Burgers equation, the two-dimensional vorticityequations, and the Navier-Stokes equations [5,6,8]. On the otherhand, the authors [7,10] developed a pseudospectral method byusing Riesz spherical means to get better results. In this paper,we generalize this method to two-dimensional vorticity equations.The generalized stability and the convergence are proved. Thenumerical results show the advantage of such a method.     

Local mass conservative Hermite interpolation for the solution of flow problems by a multi-domain boundary element approach

This paper presents a local Hermite radial basis function interpolation scheme for the velocity and pressure fields. The interpolation for velocity satisfies the continuity equation (mass conservative interpolation) while the pressure interpolation obeys the pressure equation. Additionally, the Dual Reciprocity Boundary Element method (DRBEM) is applied to obtain an integral representation of the Navier-Stokes equations. Then, the proposed local interpolation is used to obtain the values of the field variables and their partial derivatives at the boundary of the sub-domains. This interpolation allows one to obtain the boundary values needed for the integral formulas for velocity and pressure at some nodes within the sub-domains. In the proposed approach the boundary elements are merely used to parameterize the geometry, but not for the evaluation of the integrals as it is usually done. The presented multi-domain approach is different from the traditional ones in boundary elements because the resulting integral equations are non singular and the boundary data needed for the boundary integrals are approximated using a local interpolation. Some accurate results for simple Stokes problems and for the Navier-Stokes equations at low Reynolds numbers up to Re = 400 were obtained.     

Analysis and finite element approximation of a Ladyzhenskaya model for viscous flow in streamfunction form

In this paper we consider a model for the motion of incompressible viscous flows proposed by Ladyzhenskaya. The Ladyzhenskaya model is written in terms of the velocity and pressure while the studied model is written in terms of the streamfunction only. We derived the streamfunction equation of the Ladyzhenskaya model and present a weak formulation and show that this formulation is equivalent to the velocity–pressure formulation. We also present some existence and uniqueness results for the model. Finite element approximation procedures are presented. The discrete problem is proposed to be well posed and stable. Some error estimates are derived. We consider the 2D driven cavity flow problem and provide graphs which illustrate differences between the approximation procedure presented here and the approximation for the streamfunction form of the Navier–Stokes equations. Streamfunction contours are also displayed showing the main features of the flow.     

The Clebsch transformation and its capabilities towards fluid and solid mechanics

In fluid dynamics, Clebsch made use of the representation for the velocity field in terms of three potentials Φ, α, β in order to construct a first integral of the equations of motion in case of an inviscid flow with vortices. Apart from this, he received a self-adjoint form of the equations allowing for deriving them from a variational formulation. In latter times the Clebsch transformation has been applied to different physical problems, for instance to baroclinic flow, Maxwell equations in classical electrodynamics [1], in Magnetohydrodynamics and even quantum theory within the context of a quantization of vortex tubes. Viscous flow, however, has not yet been formulated in terms of Clebsch variables to our best knowledge. It is the aim of this paper to demonstrate how Clebsch variables can be applied to viscous flow on the one hand, leading to a first integral of Navier-Stokes equations as a first example. As a second example, solid mechanics is considered: by making use of an analogy between vortices in fluid flow on the one hand and dislocations in crystals on the other hand, a dynamic theory of dislocations can be established by using a certain modification of the Clebsch transformation. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Theoretical study of a Bénard-Marangoni problem

In this paper we prove the existence of strong solutions for the stationary Bénard-Marangoni problem in a finite domain flat on the top, bifurcating from the basic heat conductive state. The Bénard-Marangoni problem is a physical phenomenon of thermal convection in which the effects of buoyancy and surface tension are taken into account. This problem is modelled with a system of partial differential equations of the type Navier-Stokes and heat equation. The boundary conditions include crossed boundary conditions involving tangential derivatives of the temperature and normal derivatives of the velocity field. To define tangential derivatives at the boundary, intended in the trace sense, it is necessary order two derivatives in the interior of the domain and thus the boundary term contains as high derivatives as the interior term. We overcome this difficulty by considering the weak formulation, and transforming the boundary integral into an equivalent integral defined in the whole domain. This allows us to reformulate the weak problem with a temperature having only order one weak derivatives. Concerning regularity results, we obtain strong solutions for the stationary Bénard-Marangoni problem.     

Vorticity-velocity formulation for the Navier-Stokes equations: The three-dimensional case

In this article, we propose a mixed method for the vorticity-velocity formulation of the stationary Stokes and Navier-Stokes equations in space dimension three, the unknowns being the vorticity and the velocity of the fluid. We give a similar variational formulation for the nonstationary Stokes equations in the vorticity-velocity variables. © 1997 John Wiley & Sons, Inc.     

Mixed formulation of the two‐layer quasi‐geostrophic equations of the ocean

In this article, we propose a mixed variational formulation for the streamfunction vorticity potential form for the two‐layer quasi‐geostrophic model of the ocean. We prove the existence and uniqueness of solutions of the mixed variational problem. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 489–502, 1999