A div-curl formulation of the stokes boundary value problem based on a harmonic representation formula

On a two-dimensional domain, we establish a div-curl formulation for the Stokes Dirichlet boundary value problem. The derivation of this formulation is based on a Harmonic representation formula given by Kratz. Existence and uniqueness of solutions for the div-curl formulation are proved.     

A Galerkin boundary element method based on interpolatory Hermite trigonometric wavelets

A Galerkin boundary element method based on interpolatory Hermite trigonometric wavelets is presented for solving 2-D potential problems defined inside or outside of a circular boundary in this paper. In this approach, an equivalent variational form of the corresponding boundary integral equation for the potential problem is used; the trigonometric wavelets are employed as trial and test functions of the variational formulation. The analytical formulae of the matrix entries indicate that most of the matrix entries are naturally zero without any truncation technique and the system matrix is a block diagonal matrix. Each block consists of four circular submatrices. Hence the memory spaces and computational complexity of the system matrix are linear scale. This approach could be easily coupled into domain decomposition method based on variational formulation. Finally, the error estimates of the approximation solutions are given and some test examples are presented.     

Geometrically nonlinear crystal plasticity implemented into a discontinuous Galerkin element formulation

Single crystal viscoplasticity, with a regularization technique for the power law, is presented and implemented into a discontinuous Galerkin (DG) framework. Although single crystal plasticity has been extensively studied, its examination with the regularization method in combination with a DG formulation leads to a numerically efficient and robust model. The performance of the DG framework in crystal viscoplasticity is shown by an example. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A discrete Galerkin method for a hyperslngular boundary integral equation

Received on 14 August 1995. Revised on 20 August 1996. Consider solving the interior Neumann problem with D a simply-connected planar region and S=D a smooth curve.A double-layer potential is used to represent the solution,and it leads to the problem of solving a hypersingular integralequation. This integral equation is reformulated as a Cauchysingular integral equation. A discrete Galerkin method withtrigonometric polynomials is then given for its solution. Anerror analysis is given, and numerical examples complete thepaper.     

Optimizations of a fast multipole symmetric Galerkin boundary element method code

Numerical Algorithms - This paper presents some optimizations of a fast multipole symmetric Galerkin boundary element method code. Except general optimizations, the code is specially sped up for...     

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.     

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.     

On a Kohn-Vogelius like formulation of free boundary problems

The present paper is concerned with the solution of a Bernoulli type free boundary problem by means of shape optimization. Two state functions are introduced, namely one which satisfies the mixed boundary value problem, whereas the second one satisfies the pure Dirichlet problem. The shape problem under consideration is the minimization of the L ²-distance of the gradients of the state functions. We compute the corresponding shape gradient and Hessian. By the investigation of sufficient second order conditions we prove algebraic ill-posedness of the present formulation. Our theoretical findings are supported by numerical experiments.     

Error estimates for a discretized Galerkin method for a boundary integral equation in two dimensions

We present a priori and a posteriori estimates for the error between the Galerkin and a discretized Galerkin method for the boundary integral equation for the single layer potential on the square plate. Using piecewise constant finite elements on a rectangular mesh we study the error coming from numerical integration. The crucial point of our analysis is the estimation of some error constants, and we demonstrate that this is necessary if our methods are to be used. After the determination of these constants we are in the position to prove invertibility and quasioptimal convergence results for our numerical scheme, if the chosen numerical integration formulas are sufficiently precise. © 1992 John Wiley & Sons, Inc.     

Hypersingular kernel integration in 3D Galerkin boundary element method

We consider Cauchy singular and Hypersingular boundary integral equations associated with 3D potential problems defined on polygonal domains, whose solutions are approximated with a Galerkin boundary element method, related to a given triangulation of the boundary. In particular, for constant and linear shape functions, the most frequently used basis functions, we give explicit results of the analytical inner integrations and suggest suitable quadrature schemes to evaluate the outer integrals required to form the Galerkin matrix elements. These numerical indications are given after an analysis of the singularities arising in the whole integration process, which is valid also for shape functions of higher degrees.     

An a posteriori error analysis of an augmented discontinuous Galerkin formulation for Darcy flow

In this paper we develop an a posteriori error analysis for an augmented discontinuous Garlerkin formulation applied to the Darcy flow. More precisely, we derive a reliable and efficient a posteriori error estimator, which consists of residual terms. Finally, we present several numerical experiments, showing the robustness of the method and the theoretical properties of the estimator, thus confirming the capability of the corresponding adaptive algorithms to localize the inner layers, the singularities and/or the large stress regions of the exact solution.     

A boundary element formulation in linear piezoelectric problems

A new boundary element formulation is presented for linear piezoelectric problems under the assumption of electrodynamic quasi-static approximation. The domain, its boundary and the time region are fully discretized into a finite number of boundary-volume-time elements, and then a set of linear equations are derived from the governing integral equations. The solving scheme for the discretized linear equations in which the displacement vector, the electric potential and the corresponding fluxes are included as unknowns, is discussed.

---

Zusammenfassung Eine neue Randelement-Formulierung wird für piezoelektrische Probleme vorgeschlagen. Dabei wird eine elektrodynamische quasistatische Näherung angenommen. Das Gebiet, der Rand und der Zeitbereich werden in eine endliche Zahl von Rand-Volumen-Zeit-Elementen diskretisiert, und dann wird aus den beherrschenden Integralgleichungen eine Reihe linearer Gleichungen hergeleitet. Das Lösungsverfahren für die diskretisierten linearen Gleichungen, in denen der Verschiebungsvektor, das elektrische Potential und die zugehörigen Flüsse als Unbekannte eingehen, wird diskutiert.

     

A boundary integral formulation for elastically deformable particles in a viscous fluid

This paper reports a formulation and implementation of a mixed (both direct and indirect) boundary element method using the double layer and its adjoint in a form suitable for solving Stokes flow problems involving elastically deformable particles. The formulation is essentially the Completed Double Layer Boundary Element Method for solving an exterior traction problem for the surrounding fluid or solid phase, followed by an interior displacement, and a mobility problem (if required) for the elastic particles. At the heart of the method is a deflation procedure that allows iterative solution strategies to be adopted, effectively opens the way for large-scale simulations of suspensions of deformable particles to be performed. Several problems are considered, to illustrate and benchmark the method. In particular, an analytical solution for an elastic sphere in an elongational flow is derived. The stresslet calculations for an elastic sphere in shear and elongational flows indicate that elasticity of the inclusions can potentially lead to positive second normal stress difference in shear flow, and an increase in the tensile resistance in elongational flow.This work is supported by a grant from the Australian Research Grant Council. X-J F wishes to acknowledge the support of the National Natural Science Foundation of China.     

Least‐squares mixed Galerkin formulation for variable‐coefficient fractional differential equations with D‐N boundary condition

We propose a least‐squares mixed variational formulation for variable‐coefficient fractional differential equations (FDEs) subject to general Dirichlet‐Neumann boundary condition by splitting the FDE as a system of variable‐coefficient integer‐order equation and constant‐coefficient FDE. The main contributions of this article are to establish a new regularity theory of the solution expressed in terms of the smoothness of the right‐hand side only and to develop a decoupled and optimally convergent finite element procedure for the unknown and intermediate variables. Numerical analysis and experiments are conducted to verify these findings.     

A boundary element Galerkin method for a hypersingular integral equation on open surfaces

A hypersingular boundary integral equation of the first kind on an open surface piece Γ is solved approximately using the Galerkin method. As boundary elements on rectangles we use continuous, piecewise bilinear functions which vanish on the boundary of Γ. We show how to compensate for the effect of the edge and corner singularities of the true solution of the integral equation by using an appropriately graded mesh and obtain the same convergence rate as for the case of a smooth solution. We also derive asymptotic error estimates in lower-order Sobolev norms via the Aubin–Nitsche trick. Numerical experiments for the Galerkin method with piecewise linear functions on triangles demonstrate the effect of graded meshes and show experimental rates of convergence which underline the theoretical results.