Uniform convergence of multigrid V-cycle on adaptively refined finite element meshes for second order elliptic problems

In this paper we prove the uniform convergence of the standard multigrid V-cycle algorithm with the Gauss-Seidel relaxation performed only on the new nodes and their "immediate" neighbors for discrete elliptic problems on the adaptively refined finite element meshes using the newest vertex bisection algorithm. The proof depends on sharp estimates on the relationship of local mesh sizes and a new stability estimate for the space decomposition based on the Scott-Zhang interpolation operator. Extensive numerical results are reported, which confirm the theoretical analysis.     

Solution of double nonlinear problems in porous media by a combined finite volume–finite element algorithm

The combined finite volume–finite element scheme for a double nonlinear parabolic convection-dominated diffusion equation which models the variably saturated flow and contaminant transport problems in porous media is extended. Whereas the convection is approximated by a finite volume method (Multi-Point Flux Approximation), the diffusion is approximated by a finite element method. The scheme is fully implicit and involves a relaxation-regularized algorithm. Due to monotonicity and conservation properties of the approximated scheme and in view of the compactness theorem we show the convergence of the numerical scheme to the weak solution. Our scheme is applied for computing two dimensional examples with different degrees of complexity. The numerical results demonstrate that the proposed scheme gives good performance in convergence and accuracy.     

Finite volume element approximation and analysis for a kind of semiconductor device simulation

In this paper, we present a finite volume element scheme for a kind of two dimensional semiconductor device simulation. A general framework is developed for finite volume element approximation of the semiconductor problems. We construct a fully discrete finite volume element scheme based on triangulations with a piecewise linear finite element space and a general type of control volume. Optimal-order convergence in H ¹-norm is derived.     

The finite basis problem for infinite involution semigroups of triangular 2 $$\times $$ 2 matrices

Let \(T_n(\mathbb {F})\) and \(UT_n(\mathbb {F})\) be the semigroups of all upper triangular \(n\times n\) matrices and all upper triangular \(n\times n\) matrices with 0s and/or 1s on the main diagonal over a field \(\mathbb {F}\) with \(\mathsf {char}(\mathbb {F})=0\), respectively. In this paper, we address the finite basis problem for \(T_2(\mathbb {F})\) and \(UT_2(\mathbb {F})\) as involution semigroups under the skew transposition. By giving a sufficient condition under which an involution semigroup is nonfinitely based, we show that both \(T_2(\mathbb {F})\) and \(UT_2(\mathbb {F})\) are nonfinitely based, and that there is a continuum of nonfinitely based involution monoid varieties between the involution monoid variety \(\mathsf {var} UT_2(\mathbb {F})\) generated by \(UT_2(\mathbb {F})\) and the involution monoid variety \(\mathsf {var} T_2(\mathbb {F})\) generated by \(T_2(\mathbb {F})\). Moreover, \(\mathsf {var} UT_2(\mathbb {F})\) cannot be defined within \(\mathsf {var} T_2(\mathbb {F})\) by any finite set of identities.     

Fast calculation based on a spatial two-grid finite element algorithm for a nonlinear space–time fractional diffusion model

In this article, a spatial two-grid finite element (TGFE) algorithm is used to solve a two-dimensional nonlinear space–time fractional diffusion model and improve the computational efficiency. First, the second-order backward difference scheme is used to formulate the time approximation, where the time-fractional derivative is approximated by the weighted and shifted Grünwald difference operator. In order to reduce the computation time of the standard FE method, a TGFE algorithm is developed. The specific algorithm is to iteratively solve a nonlinear system on the coarse grid and then to solve a linear system on the fine grid. We prove the scheme stability of the TGFE algorithm and derive a priori error estimate with the convergence result O(Δt² + h^{r + 1 − η} + H^{2r + 2 − 2η}) . Finally, through a two-dimensional numerical calculation, we improve the computational efficiency and reduce the computation time by the TGFE algorithm.     

The critical groups of a family of graphs and elliptic curves over finite fields

Let q be a power of a prime, and E be an elliptic curve defined over  . Such curves have a classical group structure, and one can form an infinite tower of groups by considering E over field extensions for all k≥1. The critical group of a graph may be defined as the cokernel of L(G), the Laplacian matrix of G. In this paper, we compare elliptic curve groups with the critical groups of a certain family of graphs. This collection of critical groups also decomposes into towers of subgroups, and we highlight additional comparisons by using the Frobenius map of E over  . This work was partially supported by the NSF, grant DMS-0500557 during the author’s graduate school at the University of California, San Diego, and partially supported by an NSF Postdoctoral Fellowship.     

A nonconforming finite element method for a two-dimensional curl–curl and grad-div problem

A numerical method for a two-dimensional curl–curl and grad-div problem is studied in this paper. It is based on a discretization using weakly continuous P ₁ vector fields and includes two consistency terms involving the jumps of the vector fields across element boundaries. Optimal convergence rates (up to an arbitrary positive ) in both the energy norm and the L ₂ norm are established on graded meshes. The theoretical results are confirmed by numerical experiments. The work of the first author was supported in part by the National Science Foundation under Grant No. DMS-03-11790 and by the Humboldt Foundation through her Humboldt Research Award. The work of the third author was supported in part by the National Science Foundation under Grant No. DMS-06-52481.     

L ∞-error estimates of triangular mixed finite element methods for optimal control problems governed by semilinear elliptic equations

We look at L ^∞-error estimates for convex quadratic optimal control problems governed by nonlinear elliptic partial differential equations. In so doing, use is made of mixed finite element methods. The state and costate are approximated by the lowest order Raviart-Thomas mixed finite element spaces, and the control, by piecewise constant functions. L ^∞-error estimates of optimal order are derived for a mixed finite element approximation of a semilinear elliptic optimal control problem. Finally, numerical tests are presented which confirm our theoretical results.     

Stability Analysis of an Optimal Control Problem for a Hyperbolic Equation

In this paper, we derive some results concerning the continuous dependence on parameters of optimal solutions to an optimal control problem that involves a quasilinear hyperbolic differential equation with a boundary condition and a nonlinear integral functional of action; continuous dependence of such kind is sometimes referred to as stability or sensitivity. To present a sufficient condition for the continuous dependence, we use the Kuratowski–Painlevé upper limit of a sequence of sets. Also offered is a technical interpretation of the results.     

Existence and iterative algorithm of solutions for a system of generalized set-valued strongly nonlinear mixed variational-like inequalities problems in Banach spaces

In this paper, we introduce and study a new system of generalized set-valued strongly nonlinear mixed variational-like inequalities problems and its related auxiliary problems in reflexive Banach spaces. The auxiliary principle technique is applied to study the existence and iterative algorithm of solutions for the system of generalized set-valued strongly nonlinear mixed variational-like inequalities problems. Firstly, we prove the existence and uniqueness of solutions of the auxiliary problems for the system of generalized set-valued strongly nonlinear mixed variational-like inequalities problems. Secondly, an iterative algorithm for solving the system of generalized set-valued strongly nonlinear mixed variational-like inequalities problems is constructed by using this existence and uniqueness result. Finally, we show the existence of solutions of the system of generalized set-valued strongly nonlinear mixed variational-like inequalities problems and discuss the convergence analysis of this algorithm. These results improve, unify and generalize many corresponding known results given in literatures.     

The coupling of boundary integral and finite element methods for nonmonotone nonlinear problems ∗

In this paper, we apply the coupling of the boundary integral and finite element methods to study the weak solvability of certain nonmonotone nonlinear exterior boundary value problems. In order to convert the original exterior problem into an equivalent nonlocal boundary value problem on a finite region, we employ two different approaches based on the use of one and two integral equations on the coupling boundary. Existence of a solution for the associated weak formulation, and convergence properties of the corresponding Galerkin approximations are deduced from fundamental results in nonlinear functional analysis. Indeed, the main arguments of our proofs are based on a surjectivity theorem for mappings of type (S) and on the Fredholm alternative for nonlinear A-proper mappings.     

Mathematical Analysis and Pattern Formation for a Partial Immune System Modeling the Spread of an Epidemic Disease

Our motivation is a mathematical model describing the spatial propagation of an epidemic disease through a population. In this model, the pathogen diversity is structured into two clusters and then the population is divided into eight classes which permits to distinguish between the infected/uninfected population with respect to clusters. In this paper, we prove the weak and the global existence results of the solutions for the considered reaction-diffusion system with Neumann boundary. Next, mathematical Turing formulation and numerical simulations are introduced to show the pattern formation for such systems.     

Convergence analysis of finite element approximations of the Joule heating problem in three spatial dimensions

In this paper we present a finite element discretization of the Joule-heating problem. We prove existence of solution to the discrete formulation and strong convergence of the finite element solution to the weak solution, up to a sub-sequence. We also present numerical examples in three spatial dimensions. The first example demonstrates the convergence of the method in the second example we consider an engineering application.     

The analysis of a finite element method for the Navier–Stokes equations with compressibility

Summary. A finite element formulation is developed for the two dimensional nonlinear time dependent compressible Navier–Stokes equations on a bounded domain. The existence and uniqueness of the solution to the numerical formulation is proved. An error estimate for the numerical solution is obtained. Received September 9, 1997 / Revised version received August 12, 1999 / Published online July 12, 2000     

A priori error estimates for a mixed finite element discretization of the Richards’ equation

Summary. A mixed finite element discretization is applied to Richards equation, a nonlinear, possibly degenerate parabolic partial differential equation modeling water flow through porous medium. The equation is considered in its pressure formulation and includes both variably and fully saturated flow regime. Characteristic for such problems is the lack in regularity of the solution. To handle this we use a time-integrated scheme. We analyze the scheme and present error estimates showing its convergence.Mathematics Subject Classification (2000): 65M12, 65M60, 76S05, 35K65Acknowledgments. We would like to thank Markus Bause for very useful discussions and suggestions.     

Preconditioners for higher order edge finite element discretizations of Maxwell's equations

In this paper,we are concerned with the fast solvers for higher order edge finite element discretizations of Maxwell's equations.We present the preconditioners for the first family and second family of higher order N′ed′elec element equations,respectively.By combining the stable decompositions of two kinds of edge finite element spaces with the abstract theory of auxiliary space preconditioning,we prove that the corresponding condition numbers of our preconditioners are uniformly bounded on quasi-uniform grids.We also present some numerical experiments to demonstrate the theoretical results.