Love waves in a fluid-saturated porous layer under a rigid boundary and lying over an elastic half-space under gravity

In this paper, mathematical modeling of the propagation of Love waves in a fluid-saturated porous layer under a rigid boundary and lying over an elastic half-space under gravity has been considered. The equations of motion have been formulated separately for different media under suitable boundary conditions at the interface of porous layer, elastic half-space under gravity and rigid layer. Following Biot, the frequency equation has been derived which contain Whittaker’s function and its derivative that have been expanded asymptotically up to second term (for approximate result) for large argument due to small values of Biot’s gravity parameter (varying from 0 to 1). The effect of porosity and gravity of the layers in the propagation of Love waves has been studied. The effect of hydrostatic initial stress generated due to gravity in the half-space has also been shown in the phase velocity of Love waves. The phase velocity of Love waves for first two modes has been presented graphically. Frequency equations have also been derived for some particular cases, which are in perfect agreement with standard results. Subsequently the lower and upper bounds of Love wave speed have also been discussed.     

Error analysis of mixed finite element methods for wave propagation in double negative metamaterials

In this paper, we develop both semi-discrete and fully discrete mixed finite element methods for modeling wave propagation in three-dimensional double negative metamaterials. Optimal error estimates are proved for Nédélec spaces under the assumption of smooth solutions. To our best knowledge, this is the first error analysis obtained for Maxwell's equations when metamaterials are involved.     

Finite element approximation of coupled seismic and electromagnetic waves in fluid‐saturated poroviscoelastic media

This work presents a collection of global and iterative finite element procedures for the numerical approximation of coupled seismic and electromagnetic waves in 2D bounded fluid‐saturated porous media, with absorbing boundary conditions at the artificial boundaries. The equations being analyzed are the coupled Biot's equations of motion and Maxwell equations in the diffusive range. Both seismoelectric and electroseismic coupling are simultaneously included and analyzed in the model. The case of compressional and vertically polarized shear waves coupled with the transverse magnetic polarization (PSVTM‐mode) is analyzed in detail, including the derivation of a priori error estimates on the global finite element procedure and results on the convergence of a domain decomposition iterative algorithm. Later, the corresponding results for the case of horizontally polarized shear waves coupled with the transverse electric polarization (SHTE‐mode) are stated. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

A non-overlapping domain decomposition for low-frequency time-harmonic Maxwell’s equations in unbounded domains

In this paper, we are concerned with a non-overlapping domain decomposition method for solving the low-frequency time-harmonic Maxwell’s equations in unbounded domains. This method can be viewed as a coupling of finite elements and boundary elements in unbounded domains, which are decomposed into two subdomains with a spherical artificial boundary. We first introduce a discretization for the coupled variational problem by combining Nédélec edge elements of the lowest order and curvilinear elements. Then we design a D-N alternating method for solving the discrete problem. In the method, one needs only to solve the finite element problem (in a bounded domain) and calculate some boundary integrations, instead of solving a boundary integral equation. It will be shown that such iterative algorithm converges with a rate independent of the mesh size. The work of Qiya Hu was supported by Natural Science Foundation of China G10371129.     

Generation of higher-order polynomial basis of Nédélec H(curl) finite elements for Maxwell’s equations

The goal of this study is the automatic construction of a vectorial polynomial basis for Nédélec mixed finite elements, J.C. Nédélec (1980), in particular, the generation of finite elements without the expression of the polynomial basis functions, using symbolic calculus: the exhibition of basis functions has no practical interest.     

Boundary element method in Biot's linear consolidation

A boundary element method is formulated for the general theory of Biot's linear consolidation. Results for typical examples show good agreement with finite element solutions. As the cpu time in two dimensions is unsatisfactory an alternative method of calculating the stresses is suggested.     

Reconstruction of closely spaced small inhomogeneities via boundary measurements for the full time-dependent Maxwell’s equations

We consider for the full time-dependent Maxwell’s equations the inverse problem of identifying locations and certain properties of small electromagnetic inhomogeneities in a homogeneous background medium from dynamic boundary measurements on the boundary for a finite time interval.     

Low order Crouzeix-Raviart type nonconforming finite element methods for the 3D time-dependent Maxwell’s equations

Two Crouzeix-Raviart type nonconforming elements are used in a finite element scheme as well in a mixed finite element scheme for time-dependent Maxwell’s equations in three dimensions. The error estimates are obtained under anisotropic meshes, which are the same as those for conforming elements under regular meshes.     

Three-dimensional invisible cloaks with arbitrary shapes based on partial differential equation

The method of designing electromagnetic invisible cloak is usually based on the form-invariance of Maxwell’s equations in coordinate transformation. By solving the partial differential equations (PDEs) that describe how the coordinates transform, three-dimensional (3-D) electromagnetic and acoustic invisible cloaks with arbitrary shapes can be designed provided the boundary conditions of the cloaks can be determined by the corresponding transformation. Full wave simulations based on finite element method verify the designed cloaks. The proposed method can be easily used in designing other transformation media such as matter-wave cloaks.     

Dielectric waveguides of arbitrary cross sectional shape

Numerical guided mode solutions of arbitrary cross sectional shaped waveguides are obtained using a finite difference (FD) technique. The standard FD scheme is appropriately modified to capture all discontinuities, due to the change of the refractive index, across the waveguides’ interfaces taking into account the shape of each interface at the same time. The method is applied to the vector Helmholtz equation formulated to describe the electric or magnetic fields in the waveguide (one field is retrieved from the other through Maxwell’s equations). Computational cost is kept to a minimum by exploiting sparse matrix algebra. The waveguides under study have arbitrary cross sectional shape and arbitrary refractive index profile.     

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.     

3D dynamic Green’s functions in a multilayered poroelastic half-space

The complete 3D dynamic Green’s functions in the multilayered poroelastic media are presented in this study. A method of potentials in cylindrical coordinate system is applied first to decouple the Biot’s wave equations into four scalar Helmholtz equations, and then, general solutions to 3D wave propagation problems are obtained. After that, a three vector base and the propagator matrix method are introduced to treat 3D wave propagation problems in the stratified poroelastic half-space disturbed by buried sources. It is known that the original propagator algorithm has the loss-of-precision problem when the waves become evanescent. At present, an orthogonalization procedure is inserted into the matrix propagation loop to avoid the numerical difficulty of the original propagator algorithm. At last, the validity of the present approach for accurate and efficient calculating 3D dynamic Green’s functions of a multilayered poroelastic half-space is confirmed by comparing the numerical results with the known exact analytical solutions of a uniform poroelastic half-space.     

Numerical analysis of a PML model for time-dependent Maxwell’s equations

In this paper, we study a perfectly matched layer model for the three-dimensional time-dependent Maxwell’s equations. We develop both semi- and fully-discrete finite element methods for solving the truncated PML problem by Nedelec edge elements. Optimal convergence rates are proved for both semi- and fully-discrete schemes. To our knowledge, this is the first error analysis obtained for time domain finite element method for PML models.     

On Solution to an Optimal Shape Design Problem in 3-Dimensional Linear Magnetostatics

In this paper we present theoretical, computational, and practical aspects concerning 3-dimensional shape optimization governed by linear magnetostatics. The state solution is approximated by the finite element method using Nédélec elements on tetrahedra. Concerning optimization, the shape controls the interface between the air and the ferromagnetic parts while the whole domain is fixed. We prove the existence of an optimal shape. Then we state a finite element approximation to the optimization problem and prove the convergence of the approximated solutions. In the end, we solve the problem for the optimal shape of an electromagnet that arises in the research on magnetooptic effects and that was manufactured afterwards.     

Flux-splitting dispersion-relation-preserving dual-compact upwind scheme for solving the Maxwell’s equations on non-staggered grids

This study proposes a flux-splitting Maxwell’s equations solver for modeling electromagnetic waves in two-dimensional non-staggered grids. A fifth-order spatially accurate dual-compact upwind scheme was developed in a three-point grid stencil to approximate the first-order derivative term. The integrity of the proposed finite-difference time-domain method for solving TM-mode Maxwell’s equations verified using two-dimensional test problems. The benchmark Mie-scattering problem was also shown to be in good agreement with the semi-analytic result.     

Algebraic convergence for anisotropic edge elements in polyhedral domains

We study approximation errors for the h-version of Nédélec edge elements on anisotropically refined meshes in polyhedra. Both tetrahedral and hexahedral elements are considered, and the emphasis is on obtaining optimal convergence rates in the H(curl) norm for higher order elements. Two types of estimates are presented: First, interpolation error estimates for functions in anisotropic weighted Sobolev spaces. Here we consider not only the H(curl)-conforming Nédélec elements, but also the H(div)-conforming Raviart-Thomas elements which appear naturally in the discrete version of the de Rham complex. Our technique is to transport error estimates from the reference element to the physical element via highly anisotropic coordinate transformations. Second, Galerkin error estimates for the standard H(curl) approximation of time harmonic Maxwell equations. Here we use the anisotropic weighted Sobolev regularity of the solution on domains with three-dimensional edges and corners. We also prove the discrete compactness property needed for the convergence of the Maxwell eigenvalue problem. Our results generalize those of [40] to the case of polyhedral corners and higher order elements.     

Propagation of TM waves in a layer with arbitrary nonlinearity

A boundary value problem for Maxwell’s equations describing propagation of TM waves in a nonlinear dielectric layer with arbitrary nonlinearity is considered. The layer is located between two linear semi-infinite media. The problem is reduced to a nonlinear boundary eigenvalue problem for a system of second-order nonlinear ordinary differential equations. A dispersion equation for the eigenvalues of the problem (propagation constants) is derived. For a given nonlinearity function, the dispersion equation can be studied both analytically and numerically. A sufficient condition for the existence of at least one eigenvalue is formulated.     

A computer-assisted proof for photonic band gaps

We investigate photonic crystals, modeled by a spectral problem for Maxwell’s equations with periodic electric permittivity. Here, we specialize to a two-dimensional situation and to polarized waves. By Floquet–Bloch theory, the spectrum has band-gap structure, and the bands are characterized by families of eigenvalue problems on a periodicity cell, depending on a parameter k varying in the Brillouin zone K. We propose a computer-assisted method for proving the presence of band gaps: For k in a finite grid in K, we obtain eigenvalue enclosures by variational methods supported by finite element computations, and then capture all k ∈ K by a perturbation argument.     

Multilevel methods for mixed finite elements in three dimensions

In this paper we consider second order scalar elliptic boundary value problems posed over three–dimensional domains and their discretization by means of mixed Raviart–Thomas finite elements [18]. This leads to saddle point problems featuring a discrete flux vector field as additional unknown. Following Ewing and Wang [26], the proposed solution procedure is based on splitting the flux into divergence free components and a remainder. It leads to a variational problem involving solenoidal Raviart–Thomas vector fields. A fast iterative solution method for this problem is presented. It exploits the representation of divergence free vector fields as s of the –conforming finite element functions introduced by Nédélec [43]. We show that a nodal multilevel splitting of these finite element spaces gives rise to an optimal preconditioner for the solenoidal variational problem: Duality techniques in quotient spaces and modern algebraic multigrid theory [50, 10, 31] are the main tools for the proof. Received November 4, 1996 / Revised version received February 2, 1998     

A generalized thin-film equation in multidimensional space

This paper deals with the existence and nonnegativity of the weak periodic-domain solution for a degenerate fourth-order parabolic equation modeling the evolution of thin films. This study in the multidimensional domain follows the recent results related to the resolution in one-dimensional space [J.E. Rakotoson, J.M. Rakotoson, C. Verbeke, Generalized lubrification models blow-up and global existence, RACSAM 99 (2) (2005) 235–241; C. Verbeke, Quelques modèles d’équations d’évolution de surfaces: Explosion en temps fini et diverses propriétés qualitatives, Thèse de doctorat de l’Université de Poitiers (15 Décembre 2005); J.E. Rakotoson, J.M. Rakotoson, C. Verbeke, Higher-order equations related to thin film: Blow up and global existence, the influence of the initial data, (submitted for publication)].