A penalty method for grid matching in the mixed Herrmann-Miyoshi scheme

A penalty method for mixed finite element methods is formulated and studied. The Herrmann-Miyoshi scheme for the biharmonic equation is considered. The main idea is to build a perturbed problem with two parameters playing the role of penalties. The perturbed problem is constructed by replacing principal conditions in the mixed variational formulation at the interface by natural conditions containing parameters. Discretization of the perturbed problem is effected by a finite element method. Estimates for the norm of the difference between the solutions of a discrete perturbed problem and of an initial value problem are derived depending on the mesh size and penalties. Recommendations are given as to how to choose penalties so as to fit a mesh size.     

One-step smoothing Newton method for solving the mixed complementarity problem with a function

The mixed complementarity problem (denote by MCP(F)) can be reformulated as the solution of a smooth system of equations. In the paper, based on a perturbed mid function, we propose a new smoothing function, which has an important property, not satisfied by many other smoothing function. The existence and continuity of a smooth path for solving the mixed complementarity problem with a P₀ function are discussed. Then we presented a one-step smoothing Newton algorithm to solve the MCP with a P₀ function. The global convergence of the proposed algorithm is verified under mild conditions. And by using the smooth and semismooth technique, the rate of convergence of the method is proved under some suitable assumptions.     

分层网格上奇异摄动问题的一致NIPG分析

本文采用非对称内罚间断有限元方法(以下简称NIPG方法)求解一维对流扩散型奇异摄动问题.理论上证明了采用拉格朗日线性元的NIPG方法在分层网格上至多相差一个关于摄动参数对数因子的拟最优阶的一致收敛性,即在能量范数度量下其误差估计为O((log~2(1/e))/N),其中N为网格剖分中单元个数.数值算例验证了理论分析的正确性.     

MIXED FINITE ELEMENT METHODS FOR FRACTIONAL NAVIER-STOKES EQUATIONS

This paper gives the detailed numerical analysis of mixed finite element method for fractional Navier-Stokes equations.The proposed method is based on the mixed finite element method in space and a finite difference scheme in time.The stability analyses of semi-discretization scheme and fully discrete scheme are discussed in detail.Furthermore,We give the convergence analysis for both semidiscrete and flly discrete schemes and then prove that the numerical solution converges the exact one with order O(h2+k),where h and k:respectively denote the space step size and the time step size.Finally,numerical examples are presented to demonstrate the effectiveness of our numerical methods.     

A procedure for determining a spacewise dependent heat source and the initial temperature

The inverse problem of determining a spacewise dependent heat source, together with the initial temperature for the parabolic heat equation, using the usual conditions of the direct problem and information from two supplementary temperature measurements at different instants of time is studied. These spacewise dependent temperature measurements ensure that this inverse problem has a unique solution, despite the solution being unstable, hence the problem is ill-posed. We propose an iterative algorithm for the stable reconstruction of both the initial data and the source based on a sequence of well-posed direct problems for the parabolic heat equation, which are solved at each iteration step using the boundary element method. The instability is overcome by stopping the iterations at the first iteration for which the discrepancy principle is satisfied. Numerical results are presented for a typical benchmark test example, which has the input measured data perturbed by increasing amounts of random noise. The numerical results show that the proposed procedure gives accurate numerical approximations in relatively few iterations.     

Superconvergence analysis of a finite element method for a two-parameter singularly perturbed problem

We consider a singularly perturbed elliptic problem with two small independent parameters and its discretization by a finite element method using piecewise bilinear elements on a layer-adapted mesh. We analyze superconvergence property of the method as well as a postprocessing technique which yields more accurate discrete solution. Numerical tests confirm our theoretical results.     

Numerical solution of a mixed singularly perturbed parabolic-elliptic problem

A one-dimensional singularly perturbed problem of mixed type is considered. The domain under consideration is partitioned into two subdomains. In the first subdomain a parabolic reaction-diffusion problem is given and in the second one an elliptic convection-diffusion-reaction problem. The solution is decomposed into regular and singular components. The problem is discretized using an inverse-monotone finite volume method on condensed Shishkin meshes. We establish an almost second-order global pointwise convergence in the space variable, that is uniform with respect to the perturbation parameter.     

Determination of a spacewise dependent heat source

This paper investigates the inverse problem of determining a spacewise dependent heat source in the parabolic heat equation using the usual conditions of the direct problem and information from a supplementary temperature measurement at a given single instant of time. The spacewise dependent temperature measurement ensures that the inverse problem has a unique solution, but this solution is unstable, hence the problem is ill-posed. For this inverse problem, we propose an iterative algorithm based on a sequence of well-posed direct problems which are solved at each iteration step using the boundary element method (BEM). The instability is overcome by stopping the iterations at the first iteration for which the discrepancy principle is satisfied. Numerical results are presented for various typical benchmark test examples which have the input measured data perturbed by increasing amounts of random noise.     

Globally Convergent Broyden-Like Methods for Semismooth Equations and Applications to VIP, NCP and MCP

In this paper, we propose a general smoothing Broyden-like quasi-Newton method for solving a class of nonsmooth equations. Under appropriate conditions, the proposed method converges to a solution of the equation globally and superlinearly. In particular, the proposed method provides the possibility of developing a quasi-Newton method that enjoys superlinear convergence even if strict complementarity fails to hold. We pay particular attention to semismooth equations arising from nonlinear complementarity problems, mixed complementarity problems and variational inequality problems. We show that under certain conditions, the related methods based on the perturbed Fischer–Burmeister function, Chen–Harker–Kanzow–Smale smoothing function and the Gabriel–Moré class of smoothing functions converge globally and superlinearly.     

Finite difference/$H^1$-Galerkin MFE procedure for a fractional water wave model

In this article, an $H^1$-Galerkin mixed finite element (MFE) method for solving the time fractional water wave model is presented. First-order backward Euler difference method and $L1$ formula are applied to approximate integer derivative and Caputo fractional derivative with order $1/2$, respectively, and $H^1$-Galerkin mixed finite element method is used to approximate the spatial direction. The analysis of stability for fully discrete mixed finite element scheme is made and the optimal space-time orders of convergence for two unknown variables in both $H^1$-norm and $L^2$-norm are derived. Further, some computing results for a priori analysis and numerical figures based on four changed parameters in the studied problem are given to illustrate the effectiveness of the current method     

ON THE hp FINITE ELEMENT METHOD FOR THE ONE DIMENSIONAL SINGULARLY PERTURBED CONVECTION-DIFFUSION PROBLEMS

In this work, a singularly perturbed two-point boundary value problem of convection-diffusion type is considered. An hp version finite element method on a strongly graded piecewise uniform mesh of Shishkin type is used to solve the model problem. With the analytic assumption of the input data, it is shown that the method converges exponentially and the convergence is uniformly valid with respect to the singular perturbation parameter.     

MIXED FINITE ELEMENT METHOD FOR THE CONVECTION-DOMINATED DIFFUSION PROBLEMS WITH SMALL PARAMETER ε

In rids paper a mixed finite element method for the convection-dominated diffusion problems with small parameter ε is presented,the effect of the parameter ε on the approximation error is considered and a sufficient condition for optimal error estimates is derived. The paper also shows that under some conditions,the standard finite dement method only gives a hounded solution,however the mixed finite element method gives a convergent one.     

A finite element method for two–parameter singularly perturbed problems in 2D

We consider a singularly perturbed elliptic problem with two small parameters posed on the unit square. Based on a decomposition of the solution, we prove uniform convergence of a finite element method in an energy norm. The method uses piecewise bilinear functions on a layer-adapted Shishkin mesh. Numerical results confirm our theoretical analysis. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Convergence analysis of a hybrid Mann iterative scheme with perturbed mapping for variational inequalities and fixed point problems

The purpose of this article is to investigate the problem of finding a common element of the set of fixed points of a non-expansive mapping and the set of solutions of the variational inequality problem for a monotone, Lipschitz continuous mapping. We introduce a hybrid Mann iterative scheme with perturbed mapping which is based on the well-known Mann iteration method and hybrid (or outer approximation) method. We establish a strong convergence theorem for three sequences generated by this hybrid Mann iterative scheme with perturbed mapping. Utilizing this theorem, we also construct an iterative process for finding a common fixed point of two mappings, one of which is non-expansive and the other taken from the more general class of Lipschitz pseudocontractive mappings.     

On a problem with nonperiodic frequent alternation of boundary conditions imposed on fast oscillating sets

A singularly perturbed eigenvalue problem for the Laplacian in a cylinder is considered. The problem is characterized by frequent nonperiodic alternation of boundary conditions imposed on narrow strips lying on the cylinder’s lateral surface. The width of the strips is an arbitrary function of a small parameter and can oscillate rapidly, with the nature of the oscillations being arbitrary. Sharp estimates are derived for the convergence rate of the eigenvalues and eigenfunctions in the problem.     

On the boundary-value problem for a mixed-type equation: Asymptotic analysis and an effective numerical algorithm

A singularly perturbed boundary-value problem for an equation of mixed ellipticparabolic type is considered. The first part of this work is devoted to an asymptotic study of the solution to the problem. Modification of the boundary function method for mixed-type equations with small parameters at the highest derivatives is used. The second part is devoted to creating a numerical method that considers the structure of the solution for small parameter values. The idea of an approximate factorization of an elliptic operator into the product of two parabolic operators is employed. An effective numerical algorithm based on an asymptotic representation of the solution is developed.     

Error bounds for euler approximation of a state and control constrained optimal control problem 1

We examine convergence of the Euler approximation to a nonlinear optimal control problem subject to mixed state-control and pure state constraints. We prove that under smoothness, independence, controllability and coercivity conditions at a reference solution of the continuous problem, there exists a locally unique solution to the Euler approximation, for sufficiently fine discretization, which converges to the reference solution with rate proportional to the mesh size.     

A full-Newton step infeasible interior-point algorithm for monotone LCP based on a locally-kernel function

We propose a new full-Newton step infeasible interior-point algorithm for monotone linear complementarity problems based on a simple locally-kernel function. The algorithm uses the simple locally-kernel function to determine the search directions and define the neighborhood of central path. Two types of full-Newton steps are used, feasibility step and centering step. The algorithm starts from strictly feasible iterates of a perturbed problem, on its central path, and feasibility steps find strictly feasible iterates for the next perturbed problem. By using centering steps for the new perturbed problem, we obtain strictly feasible iterates close enough to the central path of the new perturbed problem. The procedure is repeated until an ?-approximate solution is found. We analyze the algorithm and obtain the complexity bound, which coincides with the best-known result for monotone linear complementarity problems.     

A Petrov-Galerkin mixed finite element method with exponential fitting

We discuss a Petrov-Galerkin mixed finite element formulation of the semiconductor continuity equations on a rectangular domain. We give error estimates for equations that are in principle degenerate in the singularly perturbed case. We give arguments that indicate that the method is also effective in the singularly perturbed case. We develop a discretization that gives a higher-order accurate solution for use in an a posteriori error estimator. © 1995 John Wiley & Sons, Inc.     

An elliptic singularly perturbed problem with two parameters II: Robust finite element solution

In this paper we consider a singularly perturbed elliptic model problem with two small parameters posed on the unit square. The problem is solved numerically by the finite element method using piecewise linear or bilinear elements on a layer-adapted Shishkin mesh. We prove that method with bilinear elements is uniformly convergent in an energy norm. Numerical results confirm our theoretical analysis.