On the coupling of regularization techniques and the boundary element method for a hemivariational inequality modelling a delamination problem

In this paper, we couple regularization techniques of nondifferentiable optimization with the h‐version of the boundary element method (h‐BEM) to solve nonsmooth variational problems arising in contact mechanics. As a model example, we consider the delamination problem. The variational formulation of this problem leads to a hemivariational inequality with a nonsmooth functional defined on the contact boundary. This problem is first regularized and then discretized by an h‐BEM. We prove convergence of the h‐BEM Galerkin solution of the regularized problem in the energy norm, provide an a priori error estimate and give a numerical examples. Copyright © 2016 John Wiley & Sons, Ltd.     

On the existence of incompressible current–vortex sheets: study of a linearized free boundary value problem

We study the initial boundary value problem resulting from the linearization of the equations of ideal incompressible magnetohydrodynamics and the jump conditions on the hypersurface of tangential discontinuity (current–vortex sheet) about an unsteady piecewise smooth solution. Under some assumptions on the unperturbed flow, we prove an energy a priori estimate for the linearized problem. Since the so‐called loss of derivatives in the normal direction to the boundary takes place even for the constant coefficients linearized problem, for the variable coefficients problem and non‐planar current–vortex sheets the natural functional setting is provided by the anisotropic weighted Sobolev space W₂^1,σ. The result of this paper is a necessary step to prove the local in time existence of solutions of the original non‐linear free boundary value problem. The uniqueness of the regular solution of this problem follows already from the a priori estimate we obtain for the linearized problem. Copyright © 2005 John Wiley & Sons, Ltd.     

The scattering of plane elastic waves by a one‐dimensional periodic surface

The two‐dimensional scattering problem for time‐harmonic plane waves in an isotropic elastic medium and an effectively infinite periodic surface is considered. A radiation condition for quasi‐periodic solutions similar to the condition utilized in the scattering of acoustic waves by one‐dimensional diffraction gratings is proposed. Under this condition, uniqueness of solution to the first and third boundary‐value problems is established. We then proceed by introducing a quasi‐periodic free field matrix of fundamental solutions for the Navier equation. The solution to the first boundary‐value problem is sought as a superposition of single‐ and double‐layer potentials defined utilizing this quasi‐periodic matrix. Existence of solution is established by showing the equivalence of the problem to a uniquely solvable second kind Fredholm integral equation. Copyright © 1999 John Wiley & Sons, Ltd.     

On the Multiplicity of Solutions Of a Differential Equation Arising in Chemical Reactor Theory

Consider the boundary value problem where β ? 0, τ ? 0. We are concerned with a mathematically rigorous numerical study of the number of solutions in any bounded portion of the positive quadrant (τ ? 0, β ? 0) of the τ, β plane. These correct computational results may then be matched with asymptotic (β→∞, τ ? 0) results developed earlier. These numerical results are based on the development of a posteriori error estimates for the numerical solution of an associated initial-value problem and a priori bounds on .     

On the solution of Maxwell's equations in polygonal domains

This paper is concerned with the structure of the singular and regular parts of the solution of time‐harmonic Maxwell's equations in polygonal plane domains and their effective numerical treatment. The asymptotic behaviour of the solution near corner points of the domain is studied by means of discrete Fourier transformation and it is proved that the solution of the boundary value problem does not belong locally to H² when the boundary of the domain has non‐acute angles. A splitting of the solution into a regular part belonging to the space H², and an explicitly described singular part is presented. For the numerical treatment of the boundary value problem, we propose a finite element discretization which combines local mesh grading and the singular field methods and derive a priori error estimates that show optimal convergence as known for the classical finite element method for problems with regular solutions. Copyright © 2006 John Wiley & Sons, Ltd.     

The hp‐version of the BEM with quasi‐uniform meshes for a three‐dimensional crack problem: The case of a smooth crack having smooth boundary curve

The article considers a three‐dimensional crack problem in linear elasticity with Dirichlet boundary conditions. The crack in this model problem is assumed to be a smooth open surface with smooth boundary curve. The hp‐version of the boundary element method with weakly singular operator is applied to approximate the unknown jump of the traction which is not L₂‐regular due to strong edge singularities. Assuming quasi‐uniform meshes and uniform distributions of polynomial degrees, we prove an a priori error estimate in the energy norm. The estimate gives an upper bound for the error in terms of the mesh size h and the polynomial degree p. It is optimal in h for any given data and quasi‐optimal in p for sufficiently smooth data. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

A singularly perturbed problem of photon transport in a time‐dependent region

The paper is devoted to the numerical study of a singularly perturbed transport linear integro‐differential equation, in a time‐dependent domain with slab geometry. After a brief summary on existence and uniqueness results for such a model, we test the error between the exact solution and its quasi‐static approximation, which satisfies a simpler equation. Copyright © 2002 John Wiley & Sons, Ltd.     

Web‐spline‐based finite element approximation of some quasi‐newtonian flows: Existence‐uniqueness and error bound

This article deals with the web‐spline‐based finite element approximation of quasi‐Newtonian flows. First, we consider the scalar elliptic p‐Laplace problem. Then, we consider quasi‐Newtonian flows where viscosity obeys power law or Carreau law. We prove well‐posedness at the continuous as well as the discrete level. We give some error bounds for the solution of quasi‐Newtonian flow problem based on the web‐spline method. Finally, we provide the numerical results for the p‐Laplace problem. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq31: 54–77, 2015     

The class of infinite dimensional neat reducts of quasi‐polyadic algebras is not axiomatizable

SC, CA, QA and QEA denote the classes of Pinter's substitution algebras, Tarski's cylindric algebras, Halmos' quasi‐polyadic algebras and quasi‐polyadic equality algebras, respectively. Let ω ≤ α < β and let K ∈ {SC,CA,QA,QEA}. We show that the class of α ‐dimensional neat reducts of algebras in K_β is not elementary. This solves a problem in [3]. Also our result generalizes results proved in [2] and [3]. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On a backward problem for fractional diffusion equation with Riemann-Liouville derivative

In the present paper, we study the initial inverse problem (backward problem) for a two-dimensional fractional differential equation with Riemann-Liouville derivative. Our model is considered in the random noise of the given data. We show that our problem is not well-posed in the sense of Hadamard. A truncated method is used to construct an approximate function for the solution (called the regularized solution). Furthermore, the error estimate of the regularized solution in L² and H^τ norms is considered and illustrated by numerical example.     

Natural hp‐BEM for the electric field integral equation with singular solutions

We apply the hp ‐version of the boundary element method (BEM) for the numerical solution of the electric field integral equation (EFIE) on a Lipschitz polyhedral surface Γ. The underlying meshes are supposed to be quasi‐uniform triangulations of Γ, and the approximations are based on either Raviart‐Thomas or Brezzi‐Douglas‐Marini families of surface elements. Nonsmoothness of Γ leads to singularities in the solution of the EFIE, severely affecting convergence rates of the BEM. However, the singular behavior of the solution can be explicitly specified using a finite set of functions (vertex‐, edge‐, and vertex‐edge singularities), which are the products of power functions and poly‐logarithmic terms. In this article, we use this fact to perform an a priori error analysis of the hp ‐BEM on quasi‐uniform meshes. We prove precise error estimates in terms of the polynomial degree p, the mesh size h, and the singularity exponents. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2012     

A uniform estimate for the ELLAM scheme for transport equations

We prove a priori error estimate in a weighted energy norm for the Eulerian‐Lagrangian localized adjoint method (ELLAM) for the transport equations, without any special refinement or numerical stabilization introduced. The estimate holds uniformly with respect to ?. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

A Priori Estimates for Solutions of Singular Perturbations with a Turning Point

A priori estimates are derived for the solutions of the boundary value problem εy″ + a(x)y′ + b(x)y = f(x), c ? x ? d, y(c) = α, y(d) = β. Here 0 < ε ? 1 is a small parameter and a(x) has a single simple zero in [c,d] (the turning point). It is shown that the solutions of this problem are uniformly bounded for ε→0 by the norms of f, α and β if and only if b(x)<0 at the turning point. However, in certain cases there are weak a priori estimates for the solutions even if this condition is not fulfilled.     

Error analysis of the L2 least‐squares finite element method for incompressible inviscid rotational flows

In this article we analyze the L² least‐squares finite element approximations to the incompressible inviscid rotational flow problem, which is recast into the velocity‐vorticity‐pressure formulation. The least‐squares functional is defined in terms of the sum of the squared L² norms of the residual equations over a suitable product function space. We first derive a coercivity type a priori estimate for the first‐order system problem that will play the crucial role in the error analysis. We then show that the method exhibits an optimal rate of convergence in the H¹ norm for velocity and pressure and a suboptimal rate of convergence in the L² norm for vorticity. A numerical example in two dimensions is presented, which confirms the theoretical error estimates. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004     

Exponential attractor for a planar shear‐thinning flow

We study the dynamics of an incompressible, homogeneous fluid of a power‐law type, with the stress tensor T = ν(1 + µ|Dv|)^p?2Dv, where Dv is a symmetric velocity gradient. We consider the two‐dimensional problem with periodic boundary conditions and p ∈ (1, 2). Under these assumptions, we estimate the fractal dimension of the exponential attractor, using the so‐called method of ??‐trajectories. Copyright © 2007 John Wiley & Sons, Ltd.     

Web‐compact spaces,Fréchet‐Urysohn groups and a Suslin closed graph theorem

A topological space X is strongly web‐compact if X admits a family {A_α: α ∈ ?^?} of relatively countably compact sets covering X and such that A_α ? A_β for α ≤ β. The main result of this paper states the following: Theorem A Let X and Y be topological groups and f a homomorphism between X and Y with closed graph. If X is Fréchet‐Urysohn and Baire and Y is strongly web‐compact, then f is continuous. This extends a result of Valdivia. We provide an example showing that the property of being strongly web‐compact is not productive. This applies to show that there are quasi‐Suslin spaces X whose product X × X is not quasi‐Suslin (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Finite volume method for the variational inequalities of first and second kinds

We propose and analyze the finite volume method for solving the variational inequalities of first and second kinds. The stability and convergence analysis are given for this method. For the elliptic obstacle problem, we derive the optimal error estimate in the H¹‐norm. For the simplified friction problem, we establish an abstract H¹‐error estimate, which implies the convergence if the exact solution u∈H¹(Ω) and the optimal error estimate if u∈H^{1 + α}(Ω),0 < α≤2. Copyright © 2015 John Wiley & Sons, Ltd.     

Exponential convergence and H‐c multiquadric collocation method for partial differential equations

The radial basis function (RBF) collocation method uses global shape functions to interpolate and collocate the approximate solution of PDEs. It is a truly meshless method as compared to some of the so‐called meshless or element‐free finite element methods. For the multiquadric and Gaussian RBFs, there are two ways to make the solution converge—either by refining the mesh size h, or by increasing the shape parameter c. While the h‐scheme requires the increase of computational cost, the c‐scheme is performed without extra effort. In this paper we establish by numerical experiment the exponential error estimate ? ～ O(λ^√c?h) where 0 < λ < 1. We also propose the use of residual error as an error indicator to optimize the selection of c. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 571–594, 2003     

BVPs on infinite intervals: A test problem,a nonstandard finite difference scheme and a posteriori error estimator

As far as the numerical solution of boundary value problems defined on an infinite interval is concerned, in this paper, we present a test problem for which the exact solution is known. Then we study an a posteriori estimator for the global error of a nonstandard finite difference scheme previously introduced by the authors. In particular, we show how Richardson extrapolation can be used to improve the numerical solution using the order of accuracy and numerical solutions from 2 nested quasi‐uniform grids. We observe that if the grids are sufficiently fine, the Richardson error estimate gives an upper bound of the global error.     

On a finite element method for three‐dimensional unsteady compressible viscous flows

In this article we analyze a finite element method for three‐dimensional unsteady compressible Navier‐Stokes equations. We prove the existence and uniqueness of the numerical solution, and obtain a priori error estimates uniform in time. Numerical computations are carried out to test the orders of accuracy in the error estimates. Blend function interpolations are applied in the calculation of numerical integrations. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 432–449, 2004.