Fourier spectral approximation to long-time behavior of three dimensional Ginzburg-Landau type equation*

In this paper, we consider a three dimensional Ginzburg–Landau type equation with a periodic initial value condition. A fully discrete Galerkin–Fourier spectral approximation scheme is constructed, and then the dynamical properties of the discrete system are analyzed. First, the existence and convergence of global attractors of the discrete system are proved by a priori estimates and error estimates of the discrete solution, and the numerical stability and convergence of the discrete scheme are proved. Furthermore, the long-time convergence and stability of the discrete scheme are proved. *This work was supported by the National Natural Science Foundation of China (No.: 10432010 and 10571010)     

A preconditioner for a kind of coupled FEM-BEM variational inequality

In this paper we are concerned with a kind of nonlinear transmission problem with Signorini contact conditions. This problem can be described by a coupled FEM-BEM variational inequality. We first develop a preconditioning gradient projection method for solving the variational inequality. Then we construct an effective domain decomposition preconditioner for the discrete system. The preconditioner makes the coupled inequality problem be decomposed into an equation problem and a “small” inequality problem, which can be solved in parallel. We give a complete analysis to the convergence speed of this iterative method.     

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.     

A dynamic viscoelastic contact problem with normal compliance,finite penetration and nonmonotone slip rate dependent friction

We consider a mathematical model which describes the dynamic evolution of a viscoelastic body in frictional contact with an obstacle. The contact is modelled with normal compliance and unilateral constraint, associated to a rate slip-dependent version of Coulomb’s law of dry friction. In order to approximate the contact conditions, we consider a regularized problem wherein the contact is modelled by a standard normal compliance condition without finite penetrations. For each problem, we derive a variational formulation and an existence result of the weak solution of the regularized problem is obtained. Next, we prove the convergence of the weak solution of the regularized problem to the weak solution of the initial nonregularized problem. Then, we introduce a fully discrete approximation of the variational problem based on a finite element method and on a second order time integration scheme. The solution of the resulting nonsmooth and nonconvex frictional contact problems is presented, based on approximation by a sequence of nonsmooth convex programming problems. Finally, some numerical simulations are provided in order to illustrate both the behaviour of the solution related to the frictional contact conditions and the convergence result.     

A Jacobi-collocation method for solving second kind Fredholm integral equations with weakly singular kernels

In this work, we propose a Jacobi-collocation method to solve the second kind linear Fredholm integral equations with weakly singular kernels. Particularly, we consider the case when the underlying solutions are sufficiently smooth. In this case, the proposed method leads to a fully discrete linear system. We show that the fully discrete integral operator is stable in both infinite and weighted square norms. Furthermore, we establish that the approximate solution arrives at an optimal convergence order under the two norms. Finally, we give some numerical examples, which confirm the theoretical prediction of the exponential rate of convergence.     

关于Stokes和线性Navier-Stokes方程的广义维数分裂迭代方法

本文研究了鞍点问题的迭代法.在Benzi等人提出的维数分裂(DS)迭代方法的基础上,提出了具有三个参数的广义维数分裂(GDS)迭代法,该方法包含了DS迭代法,理论分析表明该方法是无条件收敛的.通过对有限差分法和有限元法离散的Stokes问题及有限元法离散的Oseen问题的数值结果表明,本文所给方法是有效的.     

Convergence analysis of the FEM approximation of the first order projection method for incompressible flows with and without the inf-sup condition

In this paper we obtain convergence results for the fully discrete projection method for the numerical approximation of the incompressible Navier–Stokes equations using a finite element approximation for the space discretization. We consider two situations. In the first one, the analysis relies on the satisfaction of the inf-sup condition for the velocity-pressure finite element spaces. After that, we study a fully discrete fractional step method using a Poisson equation for the pressure. In this case the velocity-pressure interpolations do not need to accomplish the inf-sup condition and in fact we consider the case in which equal velocity-pressure interpolation is used. Optimal convergence results in time and space have been obtained in both cases.     

A TWO-LEVEL FINITE ELEMENT GALERKIN METHOD FOR THE NONSTATIONARY NAVIER-STOKES EQUATIONS II： TIME DISCRETIZATION

In this article we consider the fully discrete two-level finite element Galerkin method for the two-dimensional nonstationary incompressible Navier-Stokes equations. This method consists in dealing with the fully discrete nonlinear Navier-Stokes problem on a coarse mesh with width $H$ and the fully discrete linear generalized Stokes problem on a fine mesh with width $h << H$. Our results show that if we choose $H=O(h^{1/2}$) this method is as the same stability and convergence as the fully discrete standard finite element Galerkin method which needs dealing with the fully discrete nonlinear Navier-Stokes problem on a fine mesh with width $h$. However, our method is cheaper than the standard fully discrete finite element Galerkin method.     

Quadratic finite element approximation of the Signorini problem

Applying high order finite elements to unilateral contact variational inequalities may provide more accurate computed solutions, compared with linear finite elements. Up to now, there was no significant progress in the mathematical study of their performances. The main question is involved with the modeling of the nonpenetration Signorini condition on the discrete solution along the contact region. In this work we describe two nonconforming quadratic finite element approximations of the Poisson-Signorini problem, responding to the crucial practical concern of easy implementation, and we present the numerical analysis of their efficiency. By means of Falk's Lemma we prove optimal and quasi-optimal convergence rates according to the regularity of the exact solution.

     

Analysis of a coupled finite-infinite element method for exterior Helmholtz problems

Summary. This analysis of convergence of a coupled FEM-IEM is based on our previous work on the FEM and the IEM for exterior Helmholtz problems. The key idea is to represent both the exact and the numerical solution by the Dirichlet-to-Neumann operators that they induce on the coupling hypersurface in the exterior of an obstacle. The investigation of convergence can then be related to a spectral analysis of these DtN operators. We give a general outline of our method and then proceed to a detailed investigation of the case that the coupling surface is a sphere. Our main goal is to explore the convergence mechanism. In this context, we show well-posedness of both the continuous and the discrete models. We further show that the discrete inf-sup constants have a positive lower bound that does not depend on the number of DOF of the IEM. The proofs are based on lemmas on the spectra of the continuous and the discrete DtN operators, where the spectral characterization of the discrete DtN operator is given as a conjecture from numerical experiments. In our convergence analysis, we show algebraic (in terms of N) convergence of arbitrary order and generalize this result to exponential convergence. Received April 10, 1999 / Revised version received November 10, 1999 / Published online October 16, 2000     

Alternating direction method of multiplier for the unilateral contact problem with an automatic penalty parameter selection

We propose an alternating direction method of multiplier (ADMM) for the unilateral (frictionless) contact problem with an optimal parameter selection. We first introduce an auxiliary unknown to seprate the linear elasticity subproblem from the unilateral contact condition. Then an alternating direction is applied to the corresponding augmented Lagrangian. By eliminating the primal and auxiliary unknowns, at the discrete level, we derive a pure dual algorithm, starting point for the convergence analysis and the optimal parameter approximation. Numerical experiments are proposed to illustrate the efficiency of the proposed (optimal) penalty parameter selection method.     

A new reliable algorithm based on the sinc function for the time fractional diffusion equation

This paper deals with the numerical solution of time fractional diffusion equation. In this work, we consider the fractional derivative in the sense of Riemann-Liouville. At first, the time fractional derivative is discretized by integrating both sides of the equation with respect to the time variable and we arrive at a semi–discrete scheme. The stability and convergence of time discretized scheme are proven by using the energy method. Also we show that the convergence order of this scheme is O(τ2?α). Then we use the sinc collocation method to approximate the solution of semi–discrete scheme and show that the problem is reduced to a Sylvester matrix equation. Besides by performing some theorems, the exponential convergence rate of sinc method is illustrated. The numerical experiments are presented to show the excellent behavior and high accuracy of the proposed hybrid method in comparison with some other well known methods.     

The discrete Plateau Problem: Algorithm and numerics

We solve the problem of finding and justifying an optimal fully discrete finite element procedure for approximating minimal, including unstable, surfaces. In this paper we introduce the general framework and some preliminary estimates, develop the algorithm, and give the numerical results. In a subsequent paper we prove the convergence estimate. The algorithmic procedure is to find stationary points for the Dirichlet energy within the class of discrete harmonic maps from the discrete unit disc such that the boundary nodes are constrained to lie on a prescribed boundary curve. An integral normalisation condition is imposed, corresponding to the usual three point condition. Optimal convergence results are demonstrated numerically and theoretically for nondegenerate minimal surfaces, and the necessity for nondegeneracy is shown numerically.

     

Convergence of the modified discrete ordinates method for the anisotropic scattering transport equation

Criticality problem of nuclear tractors generally refers to an eigenvalue problem for the transport equations. In this paper, we deal with the eigenvalue of the anisotropic scattering transport equation in slab geometry. We propose a new discrete method which was called modified discrete ordinates method. It is constructed by redeveloping and improving discrete ordinates method in the space of L¹(X). Different from traditional methods, norm convergence of operator approximation is proved theoretically. Furthermore, convergence of eigenvalue approximation and the corresponding error estimation are obtained by analytical tools.     

An Algorithm for Continuous Type Optimal Location Problem

In this paper, we extend the ordinary discrete type facility location problems to continuous type ones. Unlike the discrete type facility location problem in which the objective function isn't everywhere differentiable, the objective function in the continuous type facility location problem is strictly convex and continuously differentiable. An algorithm without line search for solving the continuous type facility location problems is proposed and its global convergence, linear convergence rate is proved. Numerical experiments illustrate that the algorithm suggested in this paper have smaller amount of computation, quicker convergence rate than the gradient method and conjugate direction method in some sense.     

Analysis of two frictional viscoplastic contact problems with damage

In this work we study two quasistatic frictional contact problems arising in viscoplasticity including the mechanical damage of the material, caused by excessive stress or strain and modelled by an inclusion of parabolic type. The variational formulation is provided for both problems and the existence of a unique solution is proved for each of them. Then a fully discrete scheme is introduced using the finite element method to approximate the spatial domain and the Euler scheme to discretize the time derivatives. Error estimates are derived and, under suitable regularity assumptions, the linear convergence of the algorithm is deduced. Finally, some numerical examples are presented to show the performance of the method.     

NUMERICAL ANALYSIS OF A CONTACT PROBLEM IN RATE-TYPE VISCOPLASTICITY

In this paper, we consider numerical approximations of a contact problem in rate-type viscoplasticity. The contact conditions are described in term of a subdifferential and include as special cases some classical frictionless boundary conditions. The contact problem consists of an evolution equation coupled with a time-dependent variational inequality. Error estimates for both spatially semi-discrete and fully discrete solutions are derived and some convergence results are shown. Under appropriate regularity assumptions on the exact solution, error estimates are obtained.     

Galerkin finite element method for nonlinear fractional Schrödinger equations

In this paper, a class of nonlinear Riesz space-fractional Schrödinger equations are considered. Based on the standard Galerkin finite element method in space and Crank-Nicolson difference method in time, the semi-discrete and fully discrete systems are constructed. By Brouwer fixed point theorem and fractional Gagliardo-Nirenberg inequality, we prove the fully discrete system is uniquely solvable. Moreover, we focus on a rigorous analysis and consideration of the conservation and convergence properties for the semi-discrete and fully discrete systems. Finally, a linearized iterative finite element algorithm is introduced and some numerical examples are given to confirm the theoretical results.     

Numerical analysis of two frictionless elastic-piezoelectric contact problems

In this work, we consider two frictionless contact problems between an elastic-piezoelectric body and an obstacle. The linear elastic-piezoelectric constitutive law is employed to model the piezoelectric material and either the Signorini condition (if the obstacle is rigid) or the normal compliance condition (if the obstacle is deformable) are used to model the contact. The variational formulations are derived in a form of a coupled system for the displacement and electric potential fields. An existence and uniqueness result is recalled. Then, a discrete scheme is introduced based on the finite element method to approximate the spatial variable. Error estimates are derived on the approximate solutions and, as a consequence, the linear convergence of the algorithm is deduced under suitable regularity conditions. Finally, some two-dimensional examples are presented to demonstrate the performance of the algorithm.     

Convergence of a stabilized discontinuous Galerkin method for incompressible nonlinear elasticity

In this paper we present a discontinuous Galerkin method applied to incompressible nonlinear elastostatics in a total Lagrangian deformation-pressure formulation, for which a suitable interior penalty stabilization is applied. We prove that the proposed discrete formulation for the linearized problem is well-posed, asymptotically consistent and that it converges to the corresponding weak solution. The derived convergence rates are optimal and further confirmed by a set of numerical examples in two and three spatial dimensions.