Equivalence of u–p and ζ–ψ formulations of the time-dependent Navier–Stokes equations

This paper deals with the non-stationary incompressible Navier--Stokes equations for two-dimensional flows expressed in terms of the velocity and pressure and of the vorticity and streamfunction. The equivalence of the two formulations is demonstrated, both formally and rigorously, by virtue of a condition of compatibility between the boundary and initial values of the normal component of velocity. This condition is shown to be the only compatibility condition necessary to allow for solutions of a minimal regularity, namely H¹ for the velocity, as in most current numerical schemes relying on spatial discretizations of local type.     

A systematic formula for the asymptotic expansion of singular integrals

Asymptotic theories like the lifting-line, the slender body or the slender ship lead to lineintegrals with singular kernels. Sometimes these integrals are improper, that is to say that they are defined only by their Finite Part. To find asymptotic expansions of these integrals, the Matched Asymptotic Expansion Method is widely used along with other more specific methods depending on the kernel type. The first method is laborious and not systematic, and the other methods are sometimes too much specific to treat general cases. Moreover, all of them are not well adapted to deal with Finite Part integrals.Here, a new method is proposed to avoid the previous difficulties. This method is systematic for homogeneous kernels and gives approximations up to any order, as long as the derivative of the weight function exists at this given order. Moreover the occurrence of logarithmic terms in the expansion is explained and easily predictable. An elliptic integral and the classical lifting-line theory are treated to illustrate the ease of this method.

---

Résumé Les théories asymptotiques telles que la ligne portante, le corps élancé ou le navire de grand allongement conduisent à des intégrales curvilignes à noyaux singuliers. Parfois, ces intégrales sont impropres c'est à dire qu'elles sont définies en Parties Finies. Différentes méthodes ont été mises au point pour trouver les développements asymptotiques de ces intégrales. Généralement elles dépendent fortement de la nature du noyau, et c'est finalement la méthode des développements raccordés qui est utilisées quand le noyau est trop compliqué. Cependant, cette méthode est laborieuse et comme les précèdentes non adaptée aux intégrales défines par leur Partie Finie.Une nouvelle méthode est proposée pour surmonter ces difficultés. Cette méthode est systématique pour les noyaux homogènes et donne les approximations à tout ordre pourvu que les dérivées de la fonction poids existent jusqu'à cet ordre. De plus la présence de termes logarithmiques dans le développement est expliquée et aisément prédictible.Une intégrale elliptique, ainsi que la fameuse théorie de la ligne portante sont traités pour illustrer les possibilités de la méthode.

---

Nomenclature D domain of integration - f(x) weight function - FP Finite Part - h(x) weight function - I ,I _o bounded intervals - j, J integers - K(x, ) singular kernel - L, L integers - M integer defining the approximation order - P _k (x) Legendre polynomial - R set of real numbers - R(ß) equals 1 if is an integer and 0 if not - R _{f, J} ,R _{K, L} remainders of Taylor developments - S () equals either 1 or the sign function:sgn() - t, u, v, x variable of integration - , real numbers - homogeneity order of the kernel - F () Euler's integral (gamma function) - small parameter - [.] integer part of     

On the error estimates for the rotational pressure-correction projection methods

In this paper we study the rotational form of the pressure-correction method that was proposed by Timmermans, Minev, and Van De Vosse. We show that the rotational form of the algorithm provides better accuracy in terms of the -norm of the velocity and of the -norm of the pressure than the standard form.

     

Remarques sur les m\'ethodes de projection pour l'approximation des \'equations de Navier--Stokes

Summary. L'objet de cet article est de montrer que les estimations de convergence sur la pression pour les m\'ethodes de projection d\'ecrites dans \cite{Shen1} et \cite{Shen2} ne sont pas obtenues correctement car elles sont toutes bas\'ees sur une in\'egalit\'e fausse. Il semble qu'on ait besoin d'une convergence en de la vitesse dans pour d\'emontrer les estimations de convergence sur la pression en . La question de savoir si la m\'ethode de projection a un taux de convergence pour la pression plus \'elev\'e que reste ouverte. Received June 1, 1993     

A robust SN-DG-approximation for radiation transport in optically thick and diffusive regimes

We introduce a new discontinuous Galerkin (DG) method with reduced upwind stabilization for the linear Boltzmann equation applied to particle transport. The asymptotic analysis demonstrates that the new formulation does not suffer from the limitations of standard upwind methods in the thick diffusive regime; in particular, the new method yields the correct diffusion limit for any approximation order, including piecewise constant discontinuous finite elements. Numerical tests on well-established benchmark problems demonstrate the superiority of the new method. The improvement is particularly significant when employing piecewise constant DG approximation for which standard upwinding is known to perform poorly in the thick diffusion limit.     

<Emphasis Type="Italic">L</Emphasis><Superscript>1</Superscript>-minimization methods for Hamilton–Jacobi equations: the one-dimensional case

A new approximation technique based on L ¹-minimization is introduced. It is proven that the approximate solution converges to the viscosity solution in the case of one-dimensional stationary Hamilton–Jacobi equation with convex Hamiltonian. This material is based upon work supported by the National Science Foundation grant DMS-0510650. J.-L. Guermond is on leave from LIMSI, UPRR 3251 CNRS, BP 133, 91403 Orsay Cedex, France.     

Start‐up flow in a three‐dimensional lid‐driven cavity by means of a massively parallel direction splitting algorithm

The purpose of this paper is to validate a new highly parallelizable direction splitting algorithm. The parallelization capabilities of this algorithm are illustrated by providing a highly accurate solution for the start‐up flow in a three‐dimensional impulsively started lid‐driven cavity of aspect ratio 1 × 1 × 2 at Reynolds numbers 1000 and 5000. The computations are done in parallel (up to 1024 processors) on adapted grids of up to 2 billion nodes in three space dimensions. Velocity profiles are given at dimensionless times t = 4, 8, and 12; at least four digits are expected to be correct at Re = 1000. Copyright © 2011 John Wiley & Sons, Ltd.     

A new class of fractional step techniques for the incompressible Navier–Stokes equations using direction splitting

A new direction-splitting-based fractional time stepping is introduced for solving the incompressible Navier–Stokes equations. The main originality of the method is that the pressure correction is computed by solving a sequence of one-dimensional elliptic problems in each spatial direction. The method is very simple to program in parallel, very fast, and has exactly the same stability and convergence properties as the Poisson-based pressure-correction technique, either in standard or rotational form.     

Vorticity‐velocity formulations of the Stokes problem in 3D

This work studies the three‐dimensional Stokes problem expressed in terms of vorticity and velocity variables. We make general assumptions on the regularity and the topological structure of the flow domain: the boundary is Lipschitz and possibly non‐connected and the flow domain may be multiply connected. Upon introducing a new variational space for the vorticity, five weak formulations of the Stokes problem are obtained. All the formulations are shown to lead to well‐posed problems and to be equivalent to the primitive variable formulation. The various formulations are discussed by interpreting the test functions for the vorticity (resp. velocity) equation as vector potentials for the velocity (resp. vorticity). Of the five sets of boundary conditions derived in the paper, three are already known, but only for domains with a trivial topological structure, while the remaining two are new. Copyright © 1999 John Wiley & Sons, Ltd.     

Subgrid stabilization of Galerkin approximations of linear contraction semi‐groups of class C0 in Hilbert spaces

This article presents a stabilized Galerkin technique for approximating linear contraction semi‐groups of class C⁰ in a Hilbert space. The main result of this article is that this technique yields an optimal approximation estimate in the graph norm. The key idea is two‐fold. First, it consists in introducing an approximation space that is broken up into resolved scales and subgrid scales, so that the bilinear form associated with the generator of the semi‐group satisfies a uniform inf‐sup condition with respect to this decomposition. Second, the Galerkin approximation is slightly modified by introducing an artificial diffusion on the subgrid scales. Numerical tests show that the method applies also to nonlinear semi‐groups. © 2001 John Wiley & Sons., Inc. Numer Methods Partial Differential Eq 17: 1–25, 2001