Generalized H-η-accretive operators in Banach spaces with application to variational inclusions

In this paper, a new notion of a generalized H-η-accretive operator is introduced and studied, which provides a unifying framework for the generalized m-accretive operator and the H-η-monotone operator in Banach spaces. A resolvent operator associated with the generalized H-η-accretive operator is defined, and its Lipschitz continuity is shown. As an application, the solvability for a class of variational inclusions involving the generalized H-η-accretive operators in Banach spaces is considered. By using the technique of the resolvent mapping, an iterative algorithm for solving the variational inclusion in Banach spaces is constructed. Under some suitable conditions, it is proven that the solution for the variational inclusion and the convergence of the iterative sequence generated by the algorithm exist.     

Boundary-value problems for linear equations with a generalized invertible operator in a Banach space with basis

We consider linear boundary-value problems for operator equations with generalized invertible operators in Banach spaces that have bases. Using the technique of generalized inverse operators applied to generalized invertible operators in Banach spaces, we establish conditions for the solvability of linear boundary-value problems for these operator equations and obtain formulas for the representation of their solutions. We consider special cases of these boundary-value problems, namely, so-called n- and d-normally solvable boundary-value problems as well as normally solvable problems for Noetherian operator equations.     

Efficient algebraic multigrid solvers with elementary restriction and prolongation

An algebraic multigrid (AMG) scheme is presented for the efficient solution of large systems of coupled algebraic equations involving second-order discrete differentials. It is based on elementary (zero-order) intergrid transfer operators but exhibits convergence rates that are independent of the system bandwidth. Inconsistencies in the coarse-grid approximation are minimised using a global scaling approximation which requires no explicit geometrical information. Residual components of the error spectrum that remain poorly represented in the coarse-grid approximations are reduced by exploiting Krylof subspace methods. The scheme represents a robust, simple and cost-effective approach to the problem of slowly converging eigenmodes when low-order prolongation and restriction operators are used in multigrid algorithms. The algorithm investigated here uses a generalised conjugate residual (GCR) accelerator; it might also be described as an AMG preconditioned GCR method. It is applied to two test problems, one based on a solution of a discrete Poisson-type equation for nodal pressures in a pipe network, the other based on coupled solutions to the discrete Navier–Stokes equations for flows and pressures in a driven cavity. © 1998 John Wiley & Sons, Ltd.     

Stability of Discrete Stokes Operators in Fractional Sobolev Spaces

Using a general approximation setting having the generic properties of finite-elements, we prove uniform boundedness and stability estimates on the discrete Stokes operator in Sobolev spaces with fractional exponents. As an application, we construct approximations for the time-dependent Stokes equations with a source term in L ^p (0, T; L ^q (Ω)) and prove uniform estimates on the time derivative and discrete Laplacian of the discrete velocity that are similar to those in Sohr and von Wahl [20]. On long leave from LIMSI (CNRS-UPR 3251), BP 133, 91403, Orsay, France.     

Operator-splitting methods for the incompressible Navier-Stokes equations on non-staggered grids. Part 1: First-order schemes

New implicit finite difference schemes for solving the time-dependent incompressible Navier-Stokes equations using primitive variables and non-staggered grids are presented in this paper. A priori estimates for the discrete solution of the methods are obtained. Employing the operator approach, some requirements on the difference operators of the scheme are formulated in order to derive a scheme which is essentially consistent with the initial differential equations. The operators of the scheme inherit the fundamental properties of the corresponding differential operators and this allows a priori estimates for the discrete solution to be obtained. The estimate is similar to the corresponding one for the solution of the differential problem and guarantees boundedness of the solution. To derive the consistent scheme, special approximations for convective terms and div and grad operators are employed. Two variants of time discretization by the operator-splitting technique are considered and compared. It is shown that the derived scheme has a very weak restriction on the time step size. A lid-driven cavity flow has been predicted to examine the stability and accuracy of the schemes for Reynolds number up to 3200 on the sequence of grids with 21 × 21, 41 × 41, 81 × 81 and 161 × 161 grid points.     

A Nonlinear Galerkin/Petrov-Least Squares Mixed Element Method for the Stationary Navier-Stokes Equation

A nonlinear Galerkin/ Petrov- least squares mixed element (NGPLSME) method for the stationary Navier-Stokes equations is presented and analyzed. The scheme is that Petrov-least squares forms of residuals are added to the nonlinear Galerkin mixed element method so that it is stable for any combination of discrete velocity and pressure spaces without requiring the Babuska-Brezzi stability condition. The existence, uniqueness and convergence ( at optimal rate ) of the NGPLSME solution is proved in the case of sufficient viscosity ( or small data).     

A nonlinear Galerkin/Petrov-least squares mixed element method for the stationary Navier-Stokes equations

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｎｏｎｌｉｎｅａｒＧａｌｅｒｋｉｎｍｅｔｈｏｄｉｓａｍｕｌｔｉ＿ｌｅｖｅｌｓｃｈｅｍｅｔｏｆｉｎｄｔｈｅａｐｐｒｏｘｉｍａｔｅｓｏｌｕｔｉｏｎｆｏｒｔｈｅｄｉｓｓｉｐａｔｉｖｅＰＤＥ .ＴｈｉｓｍｅｔｈｏｄｈａｓｆｉｒｓｔｍａｉｎｌｙｂｅｅｎａｄｄｒｅｓｓｅｄｂｙＦｏｉａｓ＿Ｍａｎｌｅｙ＿Ｔｅｍａｍ[1],Ｍａｒｉｏｎ＿Ｔｅｍａｍ[2 ],Ｆｏｉａｓ＿Ｊｏｌｌｙ＿Ｋｅｖｒｅｋｉｄｉｓ＿Ｔｉｔｉ[3]ａｎｄＤｅｖｕｌｄｅｒ＿Ｍａｒｉｏｎ＿Ｔｉｔｉ[4 ]ｉｎｔｈｅｃａｓｅｏｆｓｐｅｃｔ…     

Element functions of discrete operator difference method

ＩｎｔｒｏｄｕｃｔｉｏｎＤｉｓｃｒｅｔｅｏｐｅｒａｔｏｒｗａｓｐｕｓｈｅｄｆｏｒｗａｒｄｉｎｐａｐｅｒｓ [1 ,2 ] ,ｗｈｉｃｈｔｒｉｅｄｔｏｕｎｉｆｙｆｉｎｉｔｅｅｌｅｍｅｎｔｍｅｔｈｏｄａｎｄｄｉｆｆｅｒｅｎｃｅｍｅｔｈｏｄｉｎｔｏｏｎｅｕｎｉｆｏｒｍｆｒａｍｅａｎｄｂｅｎｅｆｉｔｕｓｆｏｒｆｉｎｄｉｎｇｎｅｗｍｅｔｈｏｄｓ.ＰｒｏｆｅｓｓｏｒＬＩＲｏｎｇ＿ｈｕａｅｔａｌ.ｇａｖｅａｍｅｔｈｏｄ‘ｇｅｎｅｒａｔｅｄｉｆｆｅｒｅｎｃｅｍｅｔｈｏｄ’[3,4 ]ｉｓａｋｉｎｄｏｆｉｎｎｏｖａｔｉｏｎａｎｄｄ…     

On stabilized finite element methods for the Stokes problem in the small time step limit

Recent studies indicate that consistently stabilized methods for unsteady incompressible flows, obtained by a method of lines approach may experience difficulty when the time step is small relative to the spatial grid size. Using as a model problem the unsteady Stokes equations, we show that the semi‐discrete pressure operator associated with such methods is not uniformly coercive. We prove that for sufficiently large (relative to the square of the spatial grid size) time steps, implicit time discretizations contribute terms that stabilize this operator. However, we also prove that if the time step is sufficiently small, then the fully discrete problem necessarily leads to unstable pressure approximations. The semi‐discrete pressure operator studied in the paper also arises in pressure‐projection methods, thereby making our results potentially useful in other settings. Copyright © 2006 John Wiley & Sons, Ltd.     

A new full discrete stabilized viscosity method for transient Navier-Stokes equations

A new full discrete stabilized viscosity method for the transient Navier-Stokes equations with the high Reynolds number （small viscosity coefficient） is proposed based on the pressure projection and the extrapolated trapezoidal rule. The transient Navier-Stokes equations are fully discretized by the continuous equal-order finite elements in space and the reduced Crank-Nicolson scheme in time. The new stabilized method is stable and has many attractive properties. First, the system is stable for the equal-order combination of discrete continuous velocity and pressure spaces because of adding a pres- sure projection term. Second, the artifical viscosity parameter is added to the viscosity coefficient as a stability factor, so the system is antidiffusive. Finally, the method requires only the solution to a linear system at every time step. Stability and convergence of the method is proved. The error estimation results show that the method has a second-order accuracy, and the constant in the estimation is independent of the viscosity coefficient. The numerical results are given, which demonstrate the advantages of the method presented.     

Sensitivity analysis of generalized set-valued quasi-variational inclusion in Banach spaces

The sensitivity analysis for a class of generalized set-valued quasi-variational inclusion problems is investigated in the setting of Banach spaces. By using the resolvent operator technique, without assuming the differentiability and monotonicity of the given data, equivalence of these problems to the class of generalized resolvent equations is established.     

Interative approximation of fixed points for almost asymptotically nonexpansive type mappings in Banach spaces

IntroductionThroughoutthispaperweassumethatEisarealBanachspace ,E isthedualspaceofE ,DisanonemptysubsetofEandJ:E →2 E isthenormalizeddualitymappingdefinedbyJ(x) =f∈E :〈x ,f〉=‖x‖·‖f‖,‖f‖=‖x‖,　　 x∈E .　　Definition 1　LetT :D →Dbeamapping .1 )Tissaidtobeasymptoticallynonexpansive[1],ifthereexistsasequence kn [1 ,∞)withlimn→∞kn =1suchthat　　　　　‖Tnx-Tny‖≤kn‖x-y‖forall　　x ,y∈D ,n≥0 ;(1 )　　2 )Tissaidtobeofasymptoticallynonexpansivetype[2 ],if　　　　　lims…     

Mann and Ishikawa iterative approximation of solutions form-accretive operator equations

The purpose of this paper is to study the Mann and Ishikawa iterative approximation of solutions for m-accretive operator equations in Banach spaces. The results presented in this paper extend and improve some authors’ recent results.     

A Newton's method scheme for solving free-surface flow problems

A Newton's method scheme is described for solving the system of non-linear algebraic equations arising when finite difference approximations are applied to the Navier–Stokes equations and their associated boundary conditions. The problem studied here is the steady, buoyancy-driven motion of a deformable bubble, assumed to consist of an inviscid, incompressible gas. The linear Newton system is solved using both direct and iterative equation solvers. The numerical results are in excellent agreement with previous work, and the method achieves quadratic convergence.     

A class of fully third-order accurate projection methods for solving the incompressible Navier-Stokes equations

In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We then derive a sufficient condition for the continuous projection equations to be temporally third-order accurate approximations of the original Navier-Stokes equations by means of the local- truncation-error-analysis technique. The continuous projection equations are discretized temporally and spatially to third-order accuracy on the staggered grids, resulting in a fully third-order discrete projection scheme. The possibility to design higher-order projection methods is thus demonstrated in the present paper. A heuristic stability analysis is performed on this projection method showing the probability of its being stable. The stability of the present scheme is further verified through numerical tests. The third-order accuracy of the present projection method is validated by several numerical test cases.     

A class of fully third-order accurate projection methods for solving the incompressible Navier-Stokes equations

In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We then derive a sufficient condition for the continuous projection equations to be temporally third-order accurate approximations of the original Navier-Stokes equations by means of the local- truncation-error-analysis technique. The continuous projection equations are discretized temporally and spatially to third-order accuracy on the staggered grids, resulting in a fully third-order discrete projection scheme. The possibility to design higher-order projection methods is thus demonstrated in the present paper. A heuristic stability analysis is performed on this projection method showing the probability of its being stable. The stability of the present scheme is further verified through numerical tests. The third-order accuracy of the present projection method is validated by several numerical test cases. The project supported by the China NKBRSF (2001CB409604) The English text was polished by Yunming Chen     

Coagulation-Fragmentation Processes

We study the well-posedness of coagulation-fragmentation models with diffusion for large systems of particles. The continuous and the discrete case are considered simultaneously. In the discrete situation we are concerned with a countable system of coupled reaction-diffusion equations, whereas the continuous case amounts to an uncountable system of such equations. These problems can be handled by interpreting them as abstract vector-valued parabolic evolution equations, where the dependent variables take values in infinite-dimensional Banach spaces. Given suitable assumptions, we prove existence and uniqueness in the class of volume preserving solutions. We also derive sufficient conditions for global existence. Accepted: (August 18, 1999)     

Strong convergence theorems for nonexpansive semi-groups in Banach spaces

Some strong convergence theorems of explicit composite iteration scheme for nonexpansive semi-groups in the framework of Banach spaces are established.Results presented in the paper not only extend and improve the corresponding results of Shioji- Takahashi,Suzuki,Xu and Aleyner-Reich,but also give a partially affirmative answer to the open questions raised by Suzuki and Xu.     

A Class of Variational-Hemivariational Inequalities in Reflexive Banach Spaces

We study a new class of elliptic variational-hemivariational inequalities in reflexive Banach spaces. An inequality in the class is governed by a nonlinear operator, a convex set of constraints and two nondifferentiable functionals, among which at least one is convex. We deliver a result on existence and uniqueness of a solution to the inequality. Next, we show the continuous dependence of the solution on the data of the problem and we introduce a penalty method, for which we state and prove a convergence result. Finally, we consider a mathematical model which describes the equilibrium of an elastic body in unilateral contact with a foundation. The model leads to a variational-hemivariational inequality for the displacement field, that we analyse by using our abstract results.     

Analysis of a two-grid method for semiconductor device problem

The mathematical model of a semiconductor device is governed by a system of quasi-linear partial differential equations.The electric potential equation is approximated by a mixed finite element method,and the concentration equations are approximated by a standard Galerkin method.We estimate the error of the numerical solutions in the sense of the Lqnorm.To linearize the full discrete scheme of the problem,we present an efficient two-grid method based on the idea of Newton iteration.The main procedures are to solve the small scaled nonlinear equations on the coarse grid and then deal with the linear equations on the fine grid.Error estimation for the two-grid solutions is analyzed in detail.It is shown that this method still achieves asymptotically optimal approximations as long as a mesh size satisfies H=O(h^1/2).Numerical experiments are given to illustrate the efficiency of the two-grid method.