We consider a mathematical model which describes the contact between a linearly elastic body and an obstacle, the so-called foundation. The process is static and the contact is bilateral, i.e., there is no loss of contact. The friction is modeled with a nonmotonone law. The purpose of this work is to provide an error estimate for the Galerkin method as well as to present and compare two numerical methods for solving the resulting nonsmooth and nonconvex frictional contact problem. The first approach...

Isogeometric analysis has been developed recently to use basis functions resulting from the CAO description of the computational domain for the finite element spaces. The goal of this study is to develop an axisymmetric Finite Element PIC code in which specific spline Finite Elements are used to solve the Maxwell equations and the same spline functions serve as shape function for the particles. The computational domain itself is defined using splines...

The aim of this paper is to analyze a formulation of the eddy current problem in terms of a time-primitive of the electric field in a bounded domain with input current intensities or voltage drops as source data. To this end, we introduce a Lagrange multiplier to impose the divergence-free condition in the dielectric domain. Thus, we obtain a time-dependent weak mixed formulation leading to a degenerate parabolic problem which we prove is well-posed. We propose a finite element method for space...

This paper extends previous results on nonlinear Schwarz preconditioning (Cai and Keyes 2002) to unstructured finite element elliptic problems exploiting now nonlocal (but small) subspaces. The nonlocal finite element subspaces are associated with subdomains obtained from a non-overlapping element partitioning of the original set of elements and are coarse outside the prescribed element subdomain. The coarsening is based on a modification of the agglomeration based AMGe method proposed in Jones...

An energy analysis is carried out for the usual semidiscrete Galerkin method for the semilinear equation in the region$\Omega$ (E) ${\left(t{u}_t\right)}_t={\sum }_{i,j=1}{\left(a_{ij}\left(x\right)u_{x_i}\right)}_{x_j}-a_0\left(x\right)u+f\left(u\right)$, subject to the initial and boundary conditions, $u=0$ on $\partial \Omega$ and $u(x,0)=u_0$. (E) is degenerate at $t=0$ and thus, even in the case $f\equiv 0$, time derivatives of $u$ will blow up as $t\to 0$. Also, in the case where $f$ is locally Lipschitz, solutions of (E) can blow up for $t>0$ in finite time. Stability and convergence of the scheme in $W^{2,1}$ is shown in the linear case without assuming $u_{tt}$ (which can blow up as $t\to 0$ is...

Singularly perturbed reaction-diffusion problems exhibit in general solutions with anisotropic features, e.g. strong boundary and/or interior layers. This anisotropy is reflected in a discretization by using meshes with anisotropic elements. The quality of the numerical solution rests on the robustness of the a posteriori error estimator with respect to both, the perturbation parameters of the problem and the anisotropy of the mesh. The equilibrated residual method has been shown to provide one...

The tetrahedral stress element is introduced and two different types of a finite piecewise linear approximation of the dual elasticity problem are investigated on a polyhedral domain. Fot both types a priori error estimates $O\left(h^2\right)$ in $L_2$-norm and $O\left(h^{1/2}\right)$ in $L_{\infty }$-norm are established, provided the solution is smooth enough. These estimates are based on the fact that for any polyhedron there exists a strongly regular family of decomprositions into tetrahedra, which is proved in the paper, too.

The existence of a continuous right inverse of the divergence operator in $W{₀}^{1,p}\left(\Omega \right)ⁿ$, 1 < p < ∞, is a well known result which is basic in the analysis of the Stokes equations. The object of this paper is to show that the continuity also holds for some weighted norms. Our results are valid for Ω ⊂ ℝⁿ a bounded domain which is star-shaped with respect to a ball B ⊂ Ω. The continuity results are obtained by using an explicit solution of the divergence equation and the classical theory of singular integrals...

In this paper, we discuss an hp-discontinuous Galerkin finite element method (hp-DGFEM) for the laser surface hardening of steel, which is a constrained optimal control problem governed by a system of differential equations, consisting of an ordinary differential equation for austenite formation and a semi-linear parabolic differential equation for temperature evolution. The space discretization of the state variable is done using an hp-DGFEM, time and control discretizations are based on a discontinuous Galerkin...

In this paper, we discuss an hp-discontinuous Galerkin finite element method (hp-DGFEM) for the laser surface hardening of steel, which is a constrained optimal control problem governed by a system of differential equations, consisting of an ordinary differential equation for austenite formation and a semi-linear parabolic differential equation for temperature evolution. The space discretization of the state variable is done using an hp-DGFEM, time and control discretizations are based on a discontinuous Galerkin...

In this paper we study the finite element approximations to the parabolic and hyperbolic integrodifferential equations and present an immediate analysis for global superconvergence for these problems, without using the Ritz projection or its modified forms.

In this article, we consider the initial value problem which is obtained after a space discretization (with space step $h$) of the equations governing the solidification process of a multicomponent alloy. We propose a numerical scheme to solve numerically this initial value problem. We prove an error estimate which is not affected by the step size $h$ chosen in the space discretization. Consequently, our scheme provides global convergence without any stability condition between $h$ and the time step size...