Error analysis of the L2 least‐squares finite element method for incompressible inviscid rotational flows

In this article we analyze the L² least‐squares finite element approximations to the incompressible inviscid rotational flow problem, which is recast into the velocity‐vorticity‐pressure formulation. The least‐squares functional is defined in terms of the sum of the squared L² norms of the residual equations over a suitable product function space. We first derive a coercivity type a priori estimate for the first‐order system problem that will play the crucial role in the error analysis. We then show that the method exhibits an optimal rate of convergence in the H¹ norm for velocity and pressure and a suboptimal rate of convergence in the L² norm for vorticity. A numerical example in two dimensions is presented, which confirms the theoretical error estimates. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004     

Discrete harmonics for stream function-vorticity Stokes problem

Summary. We consider the bidimensional Stokes problem for incompressible fluids in stream function-vorticity. For this problem, the classical finite element method of degree one converges only in for the norm of the vorticity. We propose to use harmonic functions to approach the vorticity along the boundary. Discrete harmonics are functions that are used in practice to derive a new numerical method. We prove that we obtain with this numerical scheme an error of order for the norm of the vorticity. Received January, 2000 / Revised version received May 15, 2001 / Published online December 18, 2001     

An implicit–explicit residual error estimator for the coupling of dual‐mixed finite elements and boundary elements in elastostatics

We consider the coupling of dual‐mixed finite elements and boundary elements to solve a mixed Dirichlet–Neumann problem of plane elasticity. We derive an a‐posteriori error estimate that is based on the solution of local Dirichlet problems and on a residual term defined on the coupling interface. The general error estimate does not make use of any special finite element or boundary element spaces. Here the residual term is given in a negative order Sobolev norm. In practical applications, where a certain boundary element subspace is used, this norm can be estimated by weighted local L²‐norms. Copyright © 2001 John Wiley & Sons, Ltd.     

Formulation of boundary conditions for vorticity in viscous incompressible flow problems

For the stream function-vorticity formulation of the Navier-Stokes equations, vorticity boundary conditions are required on the body surface and the far-field boundary. A two-parameter approximating formula is derived that relates the velocity and vorticity on the outer boundary of the computational domain. The formula is used to construct an algorithm for correcting the conventional far-field boundary conditions. Specifically, a soft boundary condition is set for the vorticity and a uniform flux is specified for the transversal velocity. A third-order accurate three-parameter formula for the vorticity on the wall is derived. The use of the formula does not degrade the convergence of the iterative process of finding the vorticity as compared with a previously derived and tested two-parameter formula.     

Chebyshev‐Legendre spectral method for solving the two‐dimensional vorticity equations with homogeneous Dirichlet conditions

The Chebyshev‐Legendre spectral method for the two‐dimensional vorticity equations is considered. The Legendre Galerkin Chebyshev collocation method is used with the Chebyshev‐Gauss collocation points. The numerical analysis results under the L²‐norm for the Chebyshev‐Legendre method of one‐dimensional case are generalized into that of the two‐dimensional case. The stability and optimal order convergence of the method are proved. Numerical results are given. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

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.     

An operator approach for an oblique derivative boundary‐transmission problem

We consider a boundary‐transmission problem for the Helmholtz equation, in a Bessel potential space setting, which arises within the context of wave diffraction theory. The boundary under consideration consists of a strip, and certain conditions are assumed on it in the form of oblique derivatives. Operator theoretical methods are used to deal with the problem and, as a consequence, several convolution type operators are constructed and associated to the problem. At the end, the well‐posedness of the problem is shown for a range of non‐critical regularity orders of the Bessel potential spaces, which include the finite energy norm space. In addition, an operator normalization method is applied to the critical orders case. Copyright © 2004 John Wiley & Sons, Ltd.     

Quasi‐Lipschitz condition in potential theory

The velocity of an incompressible flow in a bounded three‐dimensional domain is represented by its vorticity with the help of an apparently new representation formula. Using this formula we prove a quasi‐Lipschitz estimate for in dependence of the supremum norm of . Our quasi‐Lipschitz bound extends to the case where is represented by any continuous ≠ rot     

On the Inviscid Limit Problem of the Vorticity Equations for Viscous Incompressible Flows in the Half‐Plane

We consider the Navier‐Stokes equations for viscous incompressible flows in the half‐plane under the no‐slip boundary condition. By using the vorticity formulation we prove the local‐in‐time convergence of the Navier‐Stokes flows to the Euler flows outside a boundary layer and to the Prandtl flows in the boundary layer in the inviscid limit when the initial vorticity is located away from the boundary. © 2014 Wiley Periodicals, Inc.     

Enhanced algorithms for implementing the spectral discretization of the vorticity‐velocity‐pressure formulation of Stokes problem

In this paper, we deal with the numerical resolution of spectral discretization of the vorticity‐velocity‐pressure formulation of Stokes problem in a square or a cube provided with nonstandard boundary conditions, which involve the normal component of the velocity and the tangential components of the vorticity. Therefore, we propose two algorithms: the Uzawa algorithm and the global resolution. We implemented the two algorithms and compared their results. With global resolution, we obtained a very good accuracy with a small number of iteration.     

Two-parameter extremum problems of boundary control for stationary thermal convection equations

Two-parameter extremum problems of boundary control are formulated for the stationary thermal convection equations with Dirichlet boundary conditions for velocity and with mixed boundary conditions for temperature. The cost functional is defined as the root mean square integral deviation of the desired velocity (vorticity, or pressure) field from one given in some part of the flow region. Controls are the boundary functions involved in the Dirichlet condition for velocity on the boundary of the flow region and in the Neumann condition for temperature on part of the boundary. The uniqueness of the extremum problems is analyzed, and the stability of solutions with respect to certain perturbations in the cost functional and one of the functional parameters of the original model is estimated. Numerical results for a control problem associated with the minimization of the vorticity norm aimed at drag reduction are discussed.     

Spectral method for vorticity‐stream function form of Navier–Stokes equations with slip boundary conditions

In this paper, we propose a spectral method for the vorticity‐stream function form of the Navier–Stokes equations with slip boundary conditions. The numerical solutions fulfill the incompressibility and the physical boundary conditions automatically. The stability and convergence of the proposed methods are proven. Numeric results demonstrate the efficiency of suggested algorithm. Copyright © 2011 John Wiley & Sons, Ltd.     

Integral method for the Stokes problemMéthode intégrale pour le problème de Stokes

We consider the bidimensional Stokes problem for incompressible fluids in stream function-vorticity formulation. For this problem, the classical finite elements method of degree one converges only in $\text{O}\text{(}\text{h}\text{)}$ for the quadratic norm of the vorticity, if the domain is convex and the solution regular. We propose to use harmonic functions obtained by a simple layer potential to approach vorticity along the boundary. Numerical results are very satisfying and we prove that this new numerical scheme leads to an error of order $\text{O}\text{(h)}$ for the natural norm of the vorticity and under more regularity assumptions from $\text{O}\text{(h}{}^{\mathrm{3/2}}\text{)}$ to $\text{O}\text{(h}{}^2\text{)}$ for the quadratic norm of the vorticity. To cite this article: T. Abboud et al., C. R. Acad. Sci. Paris, Ser. I 334 (2002) 71–76     

A Crank‐Nicolson Petrov‐Galerkin method with quadrature for semi‐linear parabolic problems

We propose and analyze an application of a fully discrete C² spline quadrature Petrov‐Galerkin method for spatial discretization of semi‐linear parabolic initial‐boundary value problems on rectangular domains. We prove second order in time and optimal order H¹ norm convergence in space for the extrapolated Crank‐Nicolson quadrature Petrov‐Galerkin scheme. We demonstrate numerically both L² and H¹ norm optimal order convergence of the scheme even if the nonlinear source term is not smooth. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005.     

Least‐squares spectral method for velocity‐vorticity‐pressure form of the Stokes equations

The aim of this article is to present and analyze first‐order system least‐squares spectral method for the Stokes equations in two‐dimensional spaces. The Stokes equations are transformed into a first‐order system of equations by introducing vorticity as a new variable. The least‐squares functional is then defined by summing up the ‐ and ‐norms of the residual equations. The ‐norm in the least‐squares functional is replaced by suitable operator. Continuous and discrete homogeneous least‐squares functionals are shown to be equivalent to ‐norm of velocity and ‐norm of vorticity and pressure for spectral Galerkin and pseudospectral method. The spectral convergence of the proposed methods are given and the theory is validated by numerical experiment. Mass conservation is also briefly investigated. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 661–680, 2016     

The hp‐version of the BEM with quasi‐uniform meshes for a three‐dimensional crack problem: The case of a smooth crack having smooth boundary curve

The article considers a three‐dimensional crack problem in linear elasticity with Dirichlet boundary conditions. The crack in this model problem is assumed to be a smooth open surface with smooth boundary curve. The hp‐version of the boundary element method with weakly singular operator is applied to approximate the unknown jump of the traction which is not L₂‐regular due to strong edge singularities. Assuming quasi‐uniform meshes and uniform distributions of polynomial degrees, we prove an a priori error estimate in the energy norm. The estimate gives an upper bound for the error in terms of the mesh size h and the polynomial degree p. It is optimal in h for any given data and quasi‐optimal in p for sufficiently smooth data. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Stability and convergence of the mark and cell finite difference scheme for Darcy‐Stokes‐Brinkman equations on non‐uniform grids

In this paper, we consider the mark and cell (MAC) method for Darcy‐Stokes‐Brinkman equations and analyze the stability and convergence of the method on nonuniform grids. Firstly, to obtain the stability for both velocity and pressure, we establish the discrete inf‐sup condition. Then we introduce an auxiliary function depending on the velocity and discretizing parameters to analyze the super‐convergence. Finally, we obtain the second‐order convergence in L2 norm for both velocity and pressure for the MAC scheme, when the perturbation parameter ? is not approaching 0. We also obtain the second‐order convergence for some terms of ∥·∥_? norm of the velocity, and the other terms of ∥·∥_? norm are second‐order convergence on uniform grid. Numerical experiments are carried out to verify the theoretical results.     

More About Geršgorin‐type theorems

In the recent book of R.S. Varga, [3], one of two main recurring themes is that a nonsingular theorem for matrices gives rise to an equivalent eigenvalue inclusion set in the complex plane, and conversely. If such nonsingularity result can be extended via irreducibility, usually this can be used for obtaining more information about the boundary of the corresponding eigenvalue inclusion set. Here we will start with one of Geršgorin‐type theorem for eigenvalue inclusion, given in [1], (for which exists corresponding equivalent statement about nonsingularity of a particular class of matrices) and use it for proving necessary conditions for an eigenvalue to lie on the boundary of localization area. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Global strong solution to the three dimensional nonhomogeneous incompressible magnetohydrodynamic equations with density‐dependent viscosity and resistivity

In this paper, we consider an initial boundary value problem for the 3‐dimensional nonhomogeneous incompressible magnetohydrodynamic equations with density‐dependent viscosity and resistivity coefficients over a bounded smooth domain. Global in time unique strong solution is proved to exist when the L² norms of initial vorticity and current density are both suitably small with arbitrary large initial density, and the vacuum of initial density is also allowed. Finally, we revisit the Navier‐Stokes model without electromagnetic effect. We find that this initial boundary problem also admits a unique global strong solution under other conditions. In particular, we prove small kinetic‐energy strong solution exists globally in time, which extends the recent result of Huang and Wang.     

Isoperimetry in Two‐Dimensional Percolation

We study isoperimetric sets, i.e., sets with minimal boundary for a prescribed volume, on the unique infinite connected component of supercritical bond percolation on the square lattice. In the limit of the volume tending to infinity, properly scaled isoperimetric sets are shown to converge (in the Hausdorff metric) to the solution of an isoperimetric problem in ?² with respect to a particular norm. As part of the proof we also show that the anchored isoperimetric profile as well as the Cheeger constant of the giant component in finite boxes scale to deterministic quantities. This settles a conjecture of Itai Benjamini for the square lattice. © 2015 Wiley Periodicals, Inc.