Discontinuous Galerkin methods for solving double obstacle problem

In this article the ideas in Wang et al. [SIAM J Numec Anal 48 (2010), 708–73] are extended to solve the double obstacle problem using discontinuous Galerkin methods. A priori error estimates are established for these methods, which reach optimal order for linear elements. We present a test example, and the numerical results on the convergence order match the theoretical prediction. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

Discontinuous Galerkin methods for solving a hyperbolic inequality

In this paper, we study spatially semi‐discrete and fully discrete schemes to numerically solve a hyperbolic variational inequality, with discontinuous Galerkin (DG) discretization in space and finite difference discretization in time. Under appropriate regularity assumptions on the solution, a unified error analysis is established for four DG schemes, which reaches the optimal convergence order for linear elements. A numerical example is presented, and the numerical results confirm the theoretical error estimates.     

Discontinuous Galerkin finite element methods for linear elasticity and quasistatic linear viscoelasticity

Summary. We consider a finite-element-in-space, and quadrature-in-time-discretization of a compressible linear quasistatic viscoelasticity problem. The spatial discretization uses a discontinous Galerkin finite element method based on polynomials of degree r—termed DG(r)—and the time discretization uses a trapezoidal-rectangle rule approximation to the Volterra (history) integral. Both semi- and fully-discrete a priori error estimates are derived without recourse to Gronwall's inequality, and therefore the error bounds do not show exponential growth in time. Moreover, the convergence rates are optimal in both h and r providing that the finite element space contains a globally continuous interpolant to the exact solution (e.g. when using the standard ^k polynomial basis on simplicies, or tensor product polynomials, ^k, on quadrilaterals). When this is not the case (e.g. using ^k on quadri-laterals) the convergence rate is suboptimal in r but remains optimal in h. We also consider a reduction of the problem to standard linear elasticity where similarly optimal a priori error estimates are derived for the DG(r) approximation. Mathematics Subject Classification (2000):65N36Shaw and Whiteman would like to acknowledge the support of the US Army Research Office, Grant #DAAD19-00-1-0421, and the UK EPSRC, Grant #GR/R10844/01. Whiteman would also like to acknowledge support from TICAM in the form of Visiting Research Fellowships.     

Discontinuous Galerkin method for solving the black‐oil problem in porous media

We introduce a new algorithm for solving the three‐component three‐phase flow problem in two‐dimensional and three‐dimensional heterogeneous media. The oil and gas components can be found in the liquid and vapor phases, whereas the aqueous phase is only composed of water component. The numerical scheme employs a sequential implicit formulation discretized with discontinuous finite elements. Capillarity and gravity effects are included. The method is shown to be accurate and robust for several test problems. It has been carefully designed so that calculation of appearance and disappearance of phases does not require additional steps.     

Analysis of a quasistatic contact problem for piezoelectric materials

We consider a mathematical model which describes the frictional contact between a piezoelectric body and an electrically conductive foundation. The process is quasistatic, the material behavior is modeled with an electro-viscoelastic constitutive law and the contact is described with subdifferential boundary conditions. We derive the variational formulation of the problem which is in the form of a system involving two history-dependent hemivariational inequalities in which the unknowns are the velocity and electric potential field. Then we prove the existence of a unique weak solution to the model. The proof is based on a recent result on history-dependent hemivariational inequalities obtained in Migórski et al. (submitted for publication) [16].     

The element‐free Galerkin method for a quasistatic contact problem with the Tresca friction in elastic materials

This paper is proposed for the error estimates of the element‐free Galerkin method for a quasistatic contact problem with the Tresca friction. The penalty method is used to impose the clamped boundary conditions. The duality algorithm is also given to deal with the non‐differentiable term in the quasistatic contact problem with the Tresca friction. The error estimates indicate that the convergence order is dependent on the nodal spacing, the time step, the largest degree of basis functions in the moving least‐squares approximation, and the penalty factor. Numerical examples demonstrate the effectiveness of the element‐free Galerkin method and verify the theoretical analysis.     

Discontinuous Galerkin finite element method for a nonlinear boundary value problem

In this paper, we investigate the a priori and a posteriori error estimates for the discontinuous Galerkin finite element approximation to a regularization version of the variational inequality of the second kind. We show the optimal error estimates in the DG-norm （stronger than the H1 norm） and the L2 norm, respectively. Furthermore, some residual-based a posteriori error estimators are established which provide global upper bounds and local lower bounds on the discretization error. These a posteriori analysis results can be applied to develop the adaptive DG methods.     

Discontinuous Galerkin methods for periodic boundary value problems

This article considers the extension of well‐known discontinuous Galerkin (DG) finite element formulations to elliptic problems with periodic boundary conditions. Such problems routinely appear in a number of applications, particularly in homogenization of composite materials. We propose an approach in which the periodicity constraint is incorporated weakly in the variational formulation of the problem. Both H¹ and L² error estimates are presented. A numerical example confirming theoretical estimates is shown. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2007     

A complementarity approach to a quasistatic multi-rigid-body contact problem

In this paper, we study the problem of predicting the quasistatic planar motion of a passive rigid body in frictional contact with a set of active rigid bodies. The active bodies can be thought of as the links of a mechanism or robot manipulator whose positions can be actively controlled by actuators. The passive body can be viewed as a grasped object, which moves only in response to contact forces and other external forces such as those due to gravity. We formulate this problem as a certain uncoupled complementarity problem, and show that it belongs to the class of NP-complete problems. Finally, numerical results of our proposed linear programming-based solution algorithm for this class of problems are presented and compared to the only other currently available solution algorithm.The research of this author was based on work supported by the National Science Foundation under grants DDM-9104078 and CCR-9213739.The research of this author was partially supported by the National Science Foundation under grant IRI-9304734, the Texas Advanced Technology Program under grant 999903-095, and the Texas Advanced Research Program under grant 999903-078.     

Discontinuous Galerkin methods for the chemotaxis and haptotaxis models

In this work, first we formulate and compare three different discontinuous Interior Penalty Galerkin methods for the 2D Keller–Segel chemotaxis model. Keller–Segel chemotaxis model is the important starting step in the modeling of the real biological system. We show in the numerical tests that two of the proposed methods fail to give accurate, oscillation-free solutions.     

Existence of a solution for a quasistatic problem of unilateral contact with local friction

The existence result in linear elasticity obtained for the quasistatic problem of unilateral contact with regularized Coulomb friction is extented to a local friction problem. After discretizing the implicit variational inequality with respect to time, we have to solve a sequence of variational inequalities similar to the one of the static problem. If the friction coefficient is small enough, we show the existence of the incremental solution. We construct a suitable sequence of functions converging towards a quasistatic solution of the problem.     

A quasistatic contact problem for viscoelastic materials with friction and damage

We consider a mathematical model which describes the frictional contact between a deformable body and a foundation. The process is quasistatic, the material is assumed to be viscoelastic with long memory and the frictional contact is modelled with subdifferential boundary conditions. The mechanical damage of the material is described by the damage function, which is modelled by a nonlinear partial differential equation. We derive the variational formulation of the problem, which is a coupled system of a hemivariational inequality and a parabolic equation. Then we prove the existence of a unique weak solution to the model. The proof is based on arguments of abstract stationary inclusion and a fixed point theorem.     

Discontinuous Galerkin finite element methods for dynamic linear solid viscoelasticity problems

We consider the usual linear elastodynamics equations augmented with evolution equations for viscoelastic internal stresses. A fully discrete approximation is defined, based on a spatially symmetric or non‐symmetric interior penalty discontinuous Galerkin finite element method, and a displacement‐velocity centred difference time discretisation. An a priori error estimate is given but only the main ideas in the proof of the error estimate are reported here due to the large number of (mostly technical) estimates that are required. The full details are referenced to a technical report. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

A quasistatic viscoplastic contact problem with normal compliance and friction

We consider a quasistatic contact problem between a viscoplasticbody and an obstacle, the so-called foundation. The contactis modelled with normal compliance and the associated versionof Coulomb's law of dry friction. We derive a variational formulationof the problem and, under a smallness assumption on the normalcompliance functions, we establish the existence of a weak solutionto the model. The proof is carried out in several steps. Itis based on a time-discretization method, arguments of monotonicityand compactness, Banach fixed point theorem and Schauder fixedpoint theorem.     

Galerkin methods for a semilinear parabolic problem with nonlocal boundary conditions

We formulate and analyze a Crank-Nicolson finite element Galerkin method and an algebraically-linear extrapolated Crank-Nicolson method for the numerical solution of a semilinear parabolic problem with nonlocal boundary conditions. For each method, optimal error estimates are derived in the maximum norm.Dedicated to Professor J. Crank on the occasion of his 80th birthdaySupported in part by the National Science Foundation grant CCR-9403461.Supported in part by project DGICYT PB95-0711.     

Discontinuous Galerkin approximations for elliptic problems

In this article, we analyze a discontinuous finite element method recently introduced by Bassi and Rebay for the approximation of elliptic problems. Stability and error estimates in various norms are proven. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 365–378, 2000     

A Galerkin method with modified piecewise polynomials for solving a second-order boundary value problem

Summary The classical Ritz-Galerkin method is applied to a linear, second-order, self-adjoint boundary value problem. The coefficient functions of the operator exhibit a piecewise smooth behaviour characteristic of some physical situations. A trial function is constructed using a modified quintic smooth Hermite space , in order to meet some desired regularity conditions for the approximate solution. A collocation technique is used to reduce the amount of computational work. Known convergence properties for the projection method are recalled which, in this particular case, are illustrated by a series of numerical experiments.     

Successive search methods for solving a canonical DC programming problem

In this article, we propose two successive search methods for solving a canonical DC programming problem constrained by the difference set between two compact convex sets in the case where the dimension number is greater than or equal to three. In order to find feasible solutions, the algorithms generate the directions based on a branch and bound procedure, successively. By exploring the provisional solutions throughout the intersection of the boundaries of two compact convex sets, both algorithms calculate an approximate solution.