A coupled mathematical system of four quasi-linear partial differential equations and the initial-boundary value conditions is presented to interpret transient behavior of three dimensional semiconductor device with heat conduction. The electric potential is defined by an elliptic equation, the electron and hole concentrations are determined by convection-dominated diffusion equations and the temperature is interpreted by a heat conduction equation. A mixed finite element approximation is used to get the electric field potential and one order of computational accuracy is improved. Two concentration equations and the heat conduction equation are solved by a fractional step scheme modified by a second-order upwind difference method, which can overcome numerical oscillation, dispersion and computational complexity. This changes the computation of a three dimensional problem into three successive computations of one-dimensional problem where the method of speedup is used and the computational work is greatly shortened. An optimal second-order error estimate in L2 norm is derived by prior estimate theory and other special techniques of partial differential equations. This type of parallel method is important in numerical analysis and is most valuable in numerical application of semiconductor device and it can successfully solve this international famous problem.

The miscible displacement of one incompressible fluid by another in a porous medium is governed by a system of two equations. One is elliptic form equation for the pressure and the other is parabolic form equation for the concentration of one of the fluids. Since only the velocity and not the pressure appears explicitly in the concentration equation, we use a mixed finite element method for the approximation of the pressure equation and mixed finite element method with characteristics for the concentration equation. To linearize the mixed-method equations, we use a two-grid algorithm based on the Newton iteration method for this full discrete scheme problems. First, we solve the original nonlinear equations on the coarse grid, then, we solve the linearized problem on the fine grid used Newton iteration once. It is shown that the coarse grid can be much coarser than the fine grid and achieve asymptotically optimal approximation as long as the mesh sizes satisfy h = H2 in this paper. Finally, numerical experiment indicates that two-grid algorithm is very effective.

In this paper, we present and study a mixed variational method in order to approximate,with the finite element method, a Stokes problem with Tresca friction boundary conditions.These non-linear boundary conditions arise in the modeling of mold filling process bypolymer melt, which can slip on a solid wall. The mixed formulation is based on adualization of the non-differentiable term which define the slip conditions. Existence anduniqueness of both continuous and discrete solutions of these problems is guaranteed bymeans of continuous and discrete inf-sup conditions that are proved. Velocity and pressureare approximated by P1 bubble-P1 finite element and piecewise linearelements are used to discretize the Lagrange multiplier associated to the shear stress onthe friction boundary. Optimal a priori error estimates are derived usingclassical tools of finite element analysis and two uncoupled discrete inf-sup conditionsfor the pressure and the Lagrange multiplier associated to the fluid shear stress.

A kind of compressible miscible displacement problems which include molecular diffusion and dispersion in porous media are investigated. The mixed finite element method is applied to the flow equation, and the transport one is solved by the symmetric interior penalty discontinuous Galerkin method. Based on a duality argument, employing projection estimates and approximation properties, a posteriori residual-type hp error estimates for the coupled system are presented, which is often used for guiding adaptivity. Comparing with the error analysis carried out by Yang (Int. J. Numer. Meth. Fluids, 65(7) (2011), pp. 781–797), the current work is more complicated and challenging.

We consider a non-standard mixed method for the Stokes problem in ℝn, n Є {2,3}, with Dirichlet boundary conditions, in which, after using the incompressibility condition to eliminate the pressure, the pseudostress tensor σ and the velocity vector u become the only unknowns. Then, we apply the Babuška-Brezzi theory to prove the well-posedness of the corresponding continuous and discrete formulations. In particular, we show that Raviart-Thomas elements of order k≥0 for σ and piecewise polynomials of degree k for u ensure unique solvability and stability of the associated Galerkin scheme. In addition, we introduce and analyze an augmented approach for our pseudostress-velocity formulation. The methodology employed is based on the introduction of the Galerkin least-squares type terms arising from the constitutive and equilibrium equations, and the Dirichlet boundary condition for the velocity, all of them multiplied by suitable stabilization parameters. We show that these parameters can be chosen so that the resulting augmented variational formulation is defined by a strongly coercive bilinear form, whence the associated Galerkin scheme becomes well posed for any choice of finite element subspaces. For instance, Raviart-Thomas elements of order k≥0 for σ and continuous piecewise polynomials of degree k+1 for u become a feasible choice in this case. Finally, extensive numerical experiments illustrating the good performance of the methods and comparing them with other procedures available in the literature, are provided.

In this paper, we develop a multiscale mortar multipoint flux mixed finite element method for second order elliptic problems. The equations in the coarse elements (or subdomains) are discretized on a fine grid scale by a multipoint flux mixed finite element method that reduces to cell-centered finite differences on irregular grids. The subdomain grids do not have to match across the interfaces. Continuity of flux between coarse elements is imposed via a mortar finite element space on a coarse grid scale. With an appropriate choice of polynomial degree of the mortar space, we derive optimal order convergence on the fine scale for both the multiscale pressure and velocity, as well as the coarse scale mortar pressure. Some superconvergence results are also derived. The algebraic system is reduced via a non-overlapping domain decomposition to a coarse scale mortar interface problem that is solved using a multiscale flux basis. Numerical experiments are presented to confirm the theory and illustrate the efficiency and flexibility of the method.

We consider an incompressible flow problem in a N-dimensional fracturedporous domain (Darcy’s problem). The fracture is represented by a(N − 1)-dimensional interface, exchanging fluid with the surroundingmedia. In this paper we consider the lowest-order(ℝ T0, ℙ0) Raviart-Thomas mixed finite elementmethod for the approximation of the coupled Darcy’s flows in the porous media and withinthe fracture, with independent meshes for the respective domains. This is achieved thanksto an enrichment with discontinuous basis functions on triangles crossed by the fractureand a weak imposition of interface conditions. First, we study the stability andconvergence properties of the resulting numerical scheme in the uncoupled case, when theknown solution of the fracture problem provides an immersed boundary condition. We detailthe implementation issues and discuss the algebraic properties of the associated linearsystem. Next, we focus on the coupled problem and propose an iterative porousdomain/fracture domain iterative method to solve for fluid flow in both the porous mediaand the fracture and compare the results with those of a traditional monolithic approach.Numerical results are provided confirming convergence rates and algebraic propertiespredicted by the theory. In particular, we discuss preconditioning and equilibrationtechniques to make the condition number of the discrete problem independent of theposition of the immersed interface. Finally, two and three dimensional simulations ofDarcy’s flow in different configurations (highly and poorly permeable fracture) areanalyzed and discussed.

In this paper, a general method to derive asymptotic error expansion formulas for the mixed finite element approximations of the Maxwell eigenvalue problem is established. Abstract lemmas for the error of the eigenvalue approximations are obtained. Based on the asymptotic error expansion formulas, the Richardson extrapolation method is employed to improve the accuracy of the approximations for the eigenvalues of the Maxwell system from to when applying the lowest order Nédélec mixed finite element and a nonconforming mixed finite element. To our best knowledge, this is the first superconvergence result of the Maxwell eigenvalue problem by the extrapolation of the mixed finite element approximation. Numerical experiments are provided to demonstrate the theoretical results.

In this paper, we investigate the superconvergence results for optimal control problems governed by parabolic equations with semidiscrete mixed finite element approximation. We use the lowest order mixed finite element spaces to discrete the state and costate variables while use piecewise constant function to discrete the control variable. Superconvergence estimates for both the state variable and its gradient variable are obtained.

We state a superconvergence result for the lowest order Raviart-Thomas approximation of eigenvalue problems. It is known that a similar superconvergence result holds for the mixed approximation of Laplace problem; here we introduce a new proof, since the one given for the source problem cannot be generalized in a straightforward way to the eigenvalue problem. Numerical experiments confirm thesuperconvergence property and suggest that it also holds for the lowest orderBrezzi-Douglas-Marini approximation.

In this paper we develop a residual based a posteriori error analysis for an augmentedmixed finite element method applied to the problem of linear elasticity in the plane.More precisely, we derive a reliable and efficient a posteriori error estimator for thecase of pure Dirichlet boundary conditions. In addition, several numericalexperiments confirming the theoretical properties of the estimator, andillustrating the capability of the corresponding adaptive algorithm to localize the singularities and the large stress regions of the solution, are also reported.

Our aim here is to study the thermal diffusion phenomenon in a forced convective flow. A system of nonlinear parabolic equations governs the evolution of the mass fractions in multicomponent mixtures. Some existence and uniqueness results are given under suitable conditions onstate functions. Then, we present a numerical scheme based on a "mixed finite element"method adapted to a finite volume scheme, of which we give numerical analysis. In a last part, we apply an homogenization technique to the studied equations in order to obtain an efficient modelling of Soret effect and adsorption in a porous medium at a macroscopic scale.

In this paper, we construct a combined upwinding and mixed finite element method for the numerical solution of a two-dimensional mean field model of superconducting vortices. An advantage of our method is that it works for any unstructured regular triangulation. A simple convergence analysis is given without resorting to the discrete maximum principle. Numerical examples are also presented.

Error estimates in L∞(0,T;L 2(Ω)), L∞(0,T;L 2(Ω)2),L∞(0,T;L ∞(Ω)), L∞(0,T;L∞ (Ω)2), Ω in ${\mathbb R}^2$ , are derived for a mixed finite element method for the initial-boundary value problem for integro-differential equation $$u_t={\rm div}\{a\bigtriangledown u+\int^t_0b_1\bigtriangledown u{\rm d}\tau+\int^t_0{\bf c}u{\rm d}\tau\}+f$$ based on the Raviart-Thomas space V h x Wh ⊂ H(div;Ω) x L2(Ω). Optimal order estimates are obtained for theapproximation of u,ut in L∞(0,T;L 2(Ω)) and the associated velocity p in L∞(0,T;L 2(Ω)2), divp in L∞(0,T;L 2(Ω)). Quasi-optimal order estimates are obtained for the approximation of u in L∞ (0,T;L ∞(Ω)) and p in L∞(0,T;L ∞(Ω)2.

In this paper, a multi-parameter error resolutiontechnique is applied into a mixed finite element method for theStokes problem. By using this technique and establishing a multi-parameter asymptotic error expansion for the mixed finite element method, an approximation of higheraccuracy is obtained by multi-processor computers in parallel.