Mixed FEM and BEM coupling for the three‐dimensional magnetostatic problem

We present two new mixed finite element methods coupled with a boundary method for the three dimensional magnetostatic problem. Such formulations are obtained by coupling a finite element method inside a bounded domain with a boundary integral method involving either the Calderon equations or the inverse of Dirichlet Neumann operator to treat the exterior domain. First, we present the formulations and then prove that our mixed formulations are well posed and that they lead to a convergent Galerkin method. Finally, we give numerical results for a sphere immersed in a homogeneous (source) field in the two formulations. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 443–462, 2003     

Asymptotic expansions for the voltage potentials with two‐dimensional and three‐dimensional thin interfaces

We derive asymptotic formulae for two‐dimensional and three‐dimensional steady state voltage potentials associated with thin conductivity imperfections having no uniform thickness. These formulae recover highly conducting inclusions and those with interfacial resistance. Our calculations are rigorous and based on layer potential techniques. Copyright © 2011 John Wiley & Sons, Ltd.     

Quasi‐local multitrace boundary integral formulations

In the context of time harmonic wave scattering by piecewise homogeneous penetrable objects, we present a new variant of the multitrace boundary integral formulation, that differs from the local multitrace approach of (Jerez‐Hanckes and Hiptmair, 2012) [13] by the presence of regularization terms involving boundary integral operators and localized at junctions i.e. points where at least three subdomains abut. We prove well‐posedness and quasioptimal convergence of conforming Galerkin discretizations for this new formulation, and present numerical results. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 2043–2062, 2015     

On the coupling of LDG-FEM and BEM methods for the three dimensional magnetostatic problem

We present two new coupling models for the three dimensional magnetostatic problem. In the first model, we propose a new coupled formulation, prove that it is well posed and solves Maxwell’s equations in the whole space. In the second, we propose a new coupled formulation for the Local Discontinuous Galerkin method, the finite element method and the boundary element method. This formulation is obtained by coupling the LDG method inside a bounded domain Ω₁ with the FEM method inside a layer where Ω is a bounded domain which is made up of material of permeability μ and such that , and with a boundary element method involving Calderon’s equations. We prove that our formulation is consistent and well posed and we present some a priori error estimates for the method.     

Vorticity‐pressure formulations for the Brinkman‐Darcy coupled problem

We introduce a new variational formulation for the Brinkman‐Darcy equations formulated in terms of the scaled Brinkman vorticity and the global pressure. The velocities in each subdomain are fully decoupled through the momentum equations, and can be later recovered from the principal unknowns. A new finite element method is also proposed, consisting in equal‐order Nédélec and piecewise continuous elements, for vorticity and pressure, respectively. The error analysis for the scheme is carried out in the natural norms, with bounds independent of the fluid viscosity. An adequate modification of the formulation and analysis permits us to specify the presentation to the case of axisymmetric configurations. We provide a set of numerical examples illustrating the robustness, accuracy, and efficiency of the proposed discretization.     

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     

Variational formulations for scattering in a three‐dimensional acoustic waveguide

Variational formulations for direct time‐harmonic scattering problems in a three‐dimensional waveguide are formulated and analyzed. We prove that the operators defined by the corresponding forms satisfy a Gårding inequality in adequately chosen spaces of test and trial functions and depend analytically on the wavenumber except at the modal numbers of the waveguide. It is also shown that these operators are strictly coercive if the wavenumber is small enough. It follows that these scattering problems are uniquely solvable except possibly for an infinite series of exceptional values of the wavenumber with no finite accumulation point. Furthermore, two geometric conditions for an obstacle are given, under which uniqueness of solution always holds in the case of a Dirichlet problem. Copyright © 2007 John Wiley & Sons, Ltd.     

Mixed boundary element solution for three-dimensional convection-diffusion problem with a velocity profile

A mixed boundary element formulation is presented for convection-diffusion problems with a velocity profile. In this formulation the convection-diffusion equation is considered as a nonlinear diffusion equation with inhomogeneous terms in which the convective term is involved additionally, because the spatial distribution of the drift velocity cannot be straightforwardly expressed in boundary integral form. Accordingly, a corresponding boundary integral equation may be described usually in the form of a so-called hybrid-type boundary integral equation.

In the present paper, mixed boundary elements are employed in a discrete model of the original convection-diffusion system. In the mixed element, potentials are approximated linearly, and their normal derivatives to boundaries are assumed constant. A simple iterative scheme is adopted in order to solve hybrid-type mixed boundary element equations. Simple three-dimensional models are dealt with in numerical experiments. The proposed approach gives more accurate and stable solutions compared with constant boundary elements which have been reported.     

A space‐time mixed‐hybrid finite element method for the damped wave equation

A space‐time finite element method is introduced to solve the linear damped wave equation. The scheme is constructed in the framework of the mixed‐hybrid finite element methods, and where an original conforming approximation of H(div;Ω) is used, the latter permits us to obtain an upwind scheme in time. We establish the link between the nonstandard finite difference scheme recently introduced by Mickens and Jordan and the scheme proposed. In this regard, two approaches are considered and in particular we employ a formulation allowing the solution to be marched in time, i.e., one only needs to consider one time increment at a time. Numerical results are presented and compared with the analytical solution illustrating good performance of the present method. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

Non - Classical Mixed Interface Problems for Anisotropic Bodies

The paper deals with the three-dimensional mathematical problems of the elasticity theory of anisotropic piece-wise homogeneous bodies. Non - classical mixed type boundary value problems are studied when on one part (S₁) of the interface the rigid contact conditions (jumps of displacement and stress vectors) are given, while conditions of the non - classical interface Problem H or Problem G are imposed on the remaining part (S₂) of the interface, i. e., in both cases jumps of the normal components of displacement and stress vectors are known on S₂ and, in addition, in the first one (Problem H) the tangent components of the displacement vector and in the second one (Problem G) the tangent components of the stress vector are given from the both sides on S₂. The investigation is carried out by the potential method and the theory of pseudodifferential equations on manifolds with boundary.     

Exponential meshes and three‐dimensional computation of a magnetic field

We describe the simulation of an exterior problem using a magnetic field deriving from magnetostatics, with a numerical method mixing the approaches of C. I. Goldstein and L.‐A. Ying. This method is based on the use of a graded mesh obtained by gluing homothetic layers in the exterior domain. On this mesh, we use an edge elements discretization and a recently proposed mixed formulation. In this paper, we provide both a theoretical and numerical study of the method. We establish an error estimate, describe the implementation, propose some preconditioning techniques and show the numerical results. We also compare these results with those obtained from an equivalent boundary elements approach. In this way, we retain that our method leads to a practical numerical implementation, a saving of storage, and turns out to be an alternative to the classical boundary elements method. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 595–637, 2003     

The method of mixed volume element‐characteristic mixed volume element and its numerical analysis for three‐dimensional slightly compressible two‐phase displacement

Numerical simulation of oil‐water two‐phase displacement is a fundamental problem in energy mathematics. The mathematical model for the compressible case is defined by a nonlinear system of two partial differential equations: (1) a parabolic equation for pressure and (2) a convection‐diffusion equation for saturation. The pressure appears within the saturation equation, and the Darcy velocity controls the saturation. The flow equation is solved by the conservative mixed volume element method. The order of the accuracy is improved by the Darcy velocity. The conservative mixed volume element with characteristics is applied to compute the saturation, that is, the diffusion is discretized by the mixed volume element and convection is computed by the method of characteristics. The method of characteristics has strong computational stability at sharp fronts and avoids numerical dispersion and nonphysical oscillation. Small time truncation error and accuracy are obtained through this method. The mixed volume element simulates diffusion, saturation, and the adjoint vector function simultaneously. By using the theory and technique of a priori estimates of differential equations, convergence of the optimal second order in norm is obtained. Numerical examples are provided to show the effectiveness and viability of this method. This method provides a powerful tool for solving challenging benchmark problems.     

A combined hybrid mixed element method for incompressible miscible displacement problem with local discontinuous Galerkin procedure

In this article, we propose a combined hybrid discontinuous mixed finite element method for miscible displacement problem with local discontinuous Galerkin method. Here, to obtain more accurate approximation and deal with the discontinuous case, we use the hybrid mixed element method to approximate the pressure and velocity, and use the local discontinuous Galerkin finite element method for the concentration. Compared with other combined methods, this method can improve the efficiency of computation, deal with the discontinuous problem well and keep local mass balance. We study the convergence of this method and give the corresponding optimal error estimates in L^∞(L²) for velocity and concentration and the super convergence in L^∞(H¹) for pressure. Finally, we also present some numerical examples to confirm our theoretical analysis.     

Combined boundary integral equations for acoustic scattering‐resonance problems

In this paper, boundary integral formulations for a time‐harmonic acoustic scattering‐resonance problem are analyzed. The eigenvalues of eigenvalue problems resulting from boundary integral formulations for scattering‐resonance problems split in general into two parts. One part consists of scattering‐resonances, and the other one corresponds to eigenvalues of some Laplacian eigenvalue problem for the interior of the scatterer. The proposed combined boundary integral formulations enable a better separation of the unwanted spectrum from the scattering‐resonances, which allows in practical computations a reliable and simple identification of the scattering‐resonances in particular for non‐convex domains. The convergence of conforming Galerkin boundary element approximations for the combined boundary integral formulations of the resonance problem is shown in canonical trace spaces. Numerical experiments confirm the theoretical results. Copyright © 2016 John Wiley & Sons, Ltd.     

Initial boundary value problem for the evolution system of MHD type describing geophysical flow in three‐dimensional domains

The initial boundary value problem for the evolution system describing geophysical flow in three‐dimensional domains was considered. The existence and uniqueness of global strong solution to the evolution system were proved under assumption on smallness of data. Moreover, solvable compatibility conditions of initial data and boundary values which guarantee the existence and uniqueness of global strong solution were discussed. Copyright © 2003 John Wiley & Sons, Ltd.     

One approach to a magnetostatic problem for bodies with inclusions in an inhomogeneous external field

A direct magnetostatic problem for magnets with a finite-size inclusion is considered in an integrodifferential form. An approach is used that, under certain conditions, reduces the problem to a single integral equation on a two-dimensional manifold-the inclusion surface. As an important illustrative example, finite formulas are derived to compute the resulting field of a magnetic half-space with a spherical cavity in an arbitrary external field.     

The boundary integral equation approach for numerical solution of the one‐dimensional Sine‐Gordon equation

This article describes a numerical method based on the boundary integral equation and dual reciprocity method for solving the one‐dimensional Sine‐Gordon (SG) equation. The time derivative is approximated by the time‐stepping method and a predictor–corrector scheme is employed to deal with the nonlinearity which appears in the problem. Numerical results are presented for some problems to demonstrate the usefulness and accuracy of this approach. In addition, the conservation of energy in SG equation is investigated. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

Mixed formulation of the two‐layer quasi‐geostrophic equations of the ocean

In this article, we propose a mixed variational formulation for the streamfunction vorticity potential form for the two‐layer quasi‐geostrophic model of the ocean. We prove the existence and uniqueness of solutions of the mixed variational problem. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 489–502, 1999     

Finite element solutions for three‐dimensional elliptic boundary value problems on unbounded domains

The finite element (FE) solutions of a general elliptic equation ?div([a_ij] ??u) + u = f in an exterior domain Ω, which is the complement of a bounded subset of R ³, is considered. The most common approach to deal with exterior domain problems is truncating an unbounded subdomain Ω_∞, so that the remaining part Ω_B = Ω\Ω _∞ is bounded, and imposing an artificial boundary condition on the resulted artificial boundary Γ_a = Ω _∞ ∩ Ω _B. In this article, instead of discarding an unbounded subdomain Ω_∞ and introducing an artificial boundary condition, the unbounded domain is mapped to a unit ball by an auxiliary mapping. Then, a similar technique to the method of auxiliary mapping, introduced by Babu?ka and Oh for handling the domain singularities, is applied to obtain an accurate FE solution of this problem at low cost. This method thus does have neither artificial boundary nor any restrictions on f. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006