Dynamics of a rotating layer of an ideal electrically conducting incompressible fluid

A system of nonlinear partial differential equations is considered that models perturbations in a layer of an ideal electrically conducting rotating fluid bounded by spatially and temporally varying surfaces with allowance for inertial forces. The system is reduced to a scalar equation. The solvability of initial boundary value problems arising in the theory of waves in conducting rotating fluids can be established by analyzing this equation. Solutions to the scalar equation are constructed that describe small-amplitude wave propagation in an infinite horizontal layer and a long narrow channel.     

Finite element approximation of viscoelastic progressively incompressible flows

Summary To avoid any numerical locking in the finite element approximation of viscoelastic flow problems, we propose a three-field approximation of this problem. This approximation, which involves velocities, stresses, and pressures is proved to converge for all times. In the proof, we also obtain convergence results for the three-fields finite element approximation of incompressible elasticity problems.     

Error bounds for the finite element approximation of an incompressible material in an unbounded domain

Summary. In this paper we design high-order local artificial boundary conditions and present error bounds for the finite element approximation of an incompressible elastic material in an unbounded domain. The finite element approximation is formulated in a bounded computational domain using a nonlocal approximate artificial boundary condition or a local one. In fact there are a family of nonlocal approximate artificial boundary conditions with increasing accuracy (and computational cost) and a family of local ones for a given artificial boundary. Our error bounds indicate how the errors of the finite element approximations depend on the mesh size, the terms used in the approximate artificial boundary condition and the location of the artificial boundary. Numerical examples of an incompressible elastic material outside a circle in the plane is presented. Numerical results demonstrate the performance of our error bounds. Received August 31, 1998 / Revised version received November 6, 2001 / Published online March 8, 2002     

Dynamics of a rotating layer of an ideal electrically conducting incompressible inhomogeneous fluid in an equatorial region

The equations describing the three-dimensional equatorial dynamics of an ideal electrically conducting inhomogeneous rotating fluid are studied. The magnetic and velocity fields are represented as superpositions of unperturbed steady-state fields and those induced by wave motion. As a result, after introducing two auxiliary functions, the equations are reduced to a special scalar one. Based on the study of this equation, the solvability of initial-boundary value problems arising in the theory of waves propagating in a spherical layer of an electrically conducting density-inhomogeneous rotating fluid in an equatorial zone is analyzed. Particular solutions of the scalar equation are constructed that describe small-amplitude wave propagation.     

Stabilized finite element schemes with LBB-stable elements for incompressible flows

We study stabilized FE approximations of SUPG type to the incompressible Navier–Stokes problem. Revisiting the analysis for the linearized model, we show that for conforming LBB-stable elements the design of the stabilization parameters for many practical flows differs from that commonly suggested in literature and initially designed for the case of equal-order approximation. Then we analyze a reduced SUPG scheme often used in practice for LBB-stable elements. To provide the reduced scheme with appropriate stability estimates we introduce a modified LBB condition which is proved for a family of FE approximations. The analysis is given for the linearized equations. Numerical experiments for some linear and nonlinear benchmark problems support the theoretical results.     

Convergence analysis of an approximation to miscible fluid flows in porous media by combining mixed finite element and finite volume methods

This article deals with development and analysis of a numerical method for a coupled system describing miscible displacement of one incompressible fluid by another through heterogeneous porous media. A mixed finite element (MFE) method is employed to discretize the Darcy flow equation combined with a conservative finite volume (FV) method on unstructured grids for the concentration equation. It is shown that the FV scheme satisfies a discrete maximum principle. We derive L^∞ and BV estimates under an appropriate CFL condition. Then we prove convergence of the approximate solutions to a weak solution of the coupled system. Numerical results are presented to see the performance of the method in two space dimensions. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

A discrete solenoidal finite difference scheme for the numerical approximation of incompressible flows

Summary For a well known class of finite difference schemes for approximating incompressible flows it is shown that the condition of discrete incompressibility can be incorporated into the discrete space. This simplifies the structure of the linear or nonlinear discrete systems and reduces the number of unknowns.     

Stabilized finite element methods for the velocity-pressure-stress formulation of incompressible flows Nonlinear model

Formulated in terms of velocity, pressure and the extra stress tensor, the incompressible Navier-Stokes is discretized by stabilized finite element method. The stabilized method proposed is analyzed for the full nonlinear model, and is applicable to various combinations of interpolation functions, including the simplest equal-order linear and bilinear elements.     

On the flow of an electrically conducting,incompressible fluid near an accelerated plate in the presence of a parallel plate,under transverse magnetic field

This paper deals with hydromagnetic flow of an electrically conducting, incompressible viscous fluid near an accelerated flat, non-conducting plate, in the presence of another parallel plate, when there is a transversely applied magnetic field. Induced magnetic field is neglected in comparison with the applied magnetic field. Laplace transform techniques are used. The equations are integrated by applying residue principle, and expressions for velocity profiles and skin-friction at both plates are derived for different values of Hartmann number M. It is observed that, with the increase of the value of the Hartmann number M, the velocity profiles are flattened, the shear stress at the stationary plate decreases, as the value of the time T and Hartmann number M increases, but the shear stress at the accelerated plate increases directly in proportion with the increase in time and Hartmann number.     

Stabilized finite element method for the viscoelastic Oldroyd fluid flows

In this paper, a new stabilized finite element method based on two local Gauss integrations is considered for the two-dimensional viscoelastic fluid motion equations, arising from the Oldroyd model for the non-Newtonian fluid flows. This new stabilized method presents attractive features such as being parameter-free, or being defined for non-edge-based data structures. It confirms that the lowest equal-order P ₁???P ₁ triangle element and Q ₁???Q ₁ quadrilateral element are compatible. Moreover, the long time stabilities and error estimates for the velocity in H ¹-norm and for the pressure in L ²-norm are obtained. Finally, some numerical experiments are performed, which show that the new method is applied to this model successfully and can save lots of computational cost compared with the standard ones.     

An approximation of incompressible miscible displacement in porous media by mixed finite element method and characteristics-mixed finite element method

An approximation scheme is defined for incompressible miscible displacement in porous media. This scheme is constructed by using two methods. Standard mixed finite element is used for the Darcy velocity equation. A characteristics-mixed finite element method is presented for the concentration equation. Characteristic approximation is applied to handle the convection part of the concentration equation, and a lowest-order mixed finite element spatial approximation is adopted to deal with the diffusion part. Thus, the scalar unknown concentration and the diffusive flux can be approximated simultaneously. In order to derive the optimal L²$L^2$-norm error estimates, a post-processing step is included in the approximation to the scalar unknown concentration. This scheme conserves mass globally; in fact, on the discrete level, fluid is transported along the approximate characteristics. Numerical experiments are presented finally to validate the theoretical analysis.     

Superconvergence of a combined mixed finite element and discontinuous Galerkin approximation for an incompressible miscible displacement problem

A combined mixed finite element and discontinuous Galerkin approximation for an incompressible miscible displacement problem which includes molecular diffusion and dispersion in porous media is studied. That is to say, the mixed finite element method is applied to the flow equation, and the transport equation is solved by an interior penalty discontinuous Galerkin method. Convolution of the Darcy velocity approximation with the Bramble-Schatz kernel function and averaging are applied in the evaluation of the coefficients in the Galerkin procedure for the concentration. A superconvergence estimate is obtained. Numerical experimental results are presented to verify the theoretical analysis.     

Analysis and finite element approximation of an optimal control problem for the Oseen viscoelastic fluid flow

In this article we study a boundary control problem for an Oseen-type model of viscoelastic fluid flow. The existence of a unique optimal solution is proved and an optimality system is derived by the first-order necessary condition. We investigate finite element approximations to a solution of the optimality system, and a solution algorithm for the system based on the gradient method.     

On the quadratic finite element approximation to the obstacle problem

Summary. In this paper, we obtain the error bound for any , for the piecewise quadratic finite element approximation to the obstacle problem, without the hypothesis that the free boundary has finite length (see [3]). Received October 31, 2000 / Revised version received July 23, 2001 / Published online October 17, 2001 The project was supported by the National Natural Science Foundation of China     

On the finite element approximation of a cascade flow problem

Summary The finite element analysis of a cascade flow problem with a given velocity circulation round profiles is presented. The nonlinear problem for the stream function with nonstandard boundary conditions is discretized by conforming linear triangular elements. We deal with the properties of the discrete problem and study the convergence of the method both for polygonal and nonpolygonal domains, including the effect of numerical integration.     

On the interface boundary conditions between two interacting incompressible viscous fluid flows

We consider two incompressible viscous fluid flows interacting through thin non-Newtonian boundary layers of higher Reynolds? number. We study the asymptotic behaviour of the problem, with respect to the vanishing thickness of the layers, using Γ-convergence methods. We derive general interfacial boundary conditions between the two fluid flows. These boundary conditions are specified for some particular cases including periodic or fractal structures of layers.     

Review and complements on mixed-hybrid finite element methods for fluid flows

Mixed and hybrid finite element methods for the resolution of a wide range of linear and nonlinear boundary value problems (linear elasticity, Stokes problem, Navier–Stokes equations, Boussinesq equations, etc.) have known a great development in the last few years. These methods allow simultaneous computation of the original variable and its gradient, both of them being equally accurate. Moreover, they have local conservation properties (conservation of the mass and the momentum) as in the finite volume methods.The purpose of this paper is to give a review on some mixed finite elements developed recently for the resolution of Stokes and Navier–Stokes equations, and the linear elasticity problem. Further developments for a quasi-Newtonian flow obeying the power law are presented.     

On a decoupled algorithm for solving a finite element problem for the approximation of viscoelastic fluid flow

Summary. The aim of this work is to study a decoupled algorithm of a fixed point for solving a finite element (FE) problem for the approximation of viscoelastic fluid flow obeying an Oldroyd B differential model. The interest for this algorithm lies in its applications to numerical simulation and in the cost of computing. Furthermore it is easy to bring this algorithm into play. The unknowns are the viscoelastic part of the extra stress tensor, the velocity and the pressure. We suppose that the solution is sufficiently smooth and small. The approximation of stress, velocity and pressure are resp. discontinuous, continuous, continuous FE. Upwinding needed for convection of , is made by discontinuous FE. The method consists to solve alternatively a transport equation for the stress, and a Stokes like problem for velocity and pressure. Previously, results of existence of the solution for the approximate problem and error bounds have been obtained using fixed point techniques with coupled algorithm. In this paper we show that the mapping of the decoupled fixed point algorithm is locally (in a neighbourhood of ) contracting and we obtain existence, unicity (locally) of the solution of the approximate problem and error bounds. Received July 29, 1994 / Revised version received March 13, 1995     

On the optimal approximation rates for criss-cross finite element spaces

Let Δ denote the triangulation of the plane obtained by multi-integer translates of the four lines x=0, y=0, x=y and x=?y. By $\text{l}{}_{\mathrm{k,\; h}}{}^{\mathrm{\mu }}$ we mean the space of all piecewise polynomials of degree ?k with respect to the scaled triangulation hΔ having continuous partial derivatives of order $\text{?μ}\text{on}\text{R}{}^2$. We show that the approximation properties of $\text{l}{}_{\mathrm{k,\; h}}{}^{\mathrm{\mu }}$ are completely governed by those of the space spanned by the translates of all so called box splines contained in $\text{l}{}_{\mathrm{k,h}}{}^{\mathrm{\mu }}$. Combining this fact with Fourier analysis techniques allows us to determine the optimal controlled approximation rates for the above subspace of box splines where μ is the largest degree of smoothness for which these spaces are dense as h tends to zero. Furthermore, we study the question of local linear dependence of the translates of the box splines for the above criss-cross triangulations.