A new methodology for anisotropic mesh refinement based upon error gradients

We introduce a new strategy for controlling the use of anisotropic mesh refinement based upon the gradients of an a posteriori approximation of the error in a computed finite element solution. The efficiency of this strategy is demonstrated using a simple anisotropic mesh adaption algorithm and the quality of a number of potential a posteriori error estimates is considered.     

Time‐harmonic and asymptotically linear Maxwell equations in anisotropic media

This paper is focused on following time‐harmonic Maxwell equation: where is a bounded Lipschitz domain, is the exterior normal, and ω is the frequency. The boundary condition holds when Ω is surrounded by a perfect conductor. Assuming that f is asymptotically linear as , we study the above equation by improving the generalized Nehari manifold method. For an anisotropic material with magnetic permeability tensor and permittivity tensor , ground state solutions are established in this paper. Applying the principle of symmetric criticality, we find 2 types of solutions with cylindrical symmetries in particular for the uniaxial material.     

A new formulation of regularized meshless method applied to interior and exterior anisotropic potential problems

This paper proposes a new formulation of regularized meshless method (RMM), which differs from the traditional RMM in that the traditional formulation generates the diagonal elements of influence matrix via null-field integral equations, while our new one directly employs the boundary integral equations at the domain point to evaluate the diagonal elements. We test the present RMM formulation to two-dimensional anisotropic potential problems in finite and infinite domains in comparison with the traditional RMM. Numerical results show that the present RMM sharply outperforms the traditional RMM in the solution of interior problems, while the latter is clearly superior for exterior problems. A rigorous theoretical analysis of circular domain case also corroborates such numerical experiment observations and is provided in the appendix of this paper.     

The No Response Test for the reconstruction of polyhedral objects in electromagnetics

We develop a No Response Test for the reconstruction of a polyhedral obstacle from two or few time-harmonic electromagnetic incident waves in electromagnetics. The basic idea of the test is to probe some region in space with waves which are small on some test domain and, thus, do not generate a response when the scatterer is inside of this test domain. We will prove that the No Response Test checks analytic continuability of a time-harmonic field from the far field pattern into the domain for a non-vibrating test domain B.We show that two incident waves, defined by one incident direction and two appropriately chosen directions of polarization, are enough to recover the convex hull of polyhedrals. Based on this uniqueness result, we build up the No Response Test and we prove convergence in the sense that it fully reconstructs a convex polyhedral scatterer D or the convex hull of an arbitrary polyhedral scatterer.Further, we will describe the algorithmic realization of the No Response Test and show the feasibility of the method by reconstruction of convex polyhedral objects in three dimensions. This is the first formulation of the No Response Test for electromagnetics.     

Finite element study of time‐dependent Maxwell's equations in dispersive media

In this article, we consider the time‐dependent Maxwell's equations in a bounded domain when dispersive media are involved. The Crank‐Nicolson scheme is developed to approximate the electric field equation by Nedelec edge elements and is proved to be optimal convergent in energy norm. The analysis is carried out for Debye medium, but the same results hold true for other dispersive media such as plasma and Lorentz medium. Furthermore, our analysis extends straightforward to cases when a dispersive medium and a simple medium (such as air) are coupled. Mathematics Subject Classification (2000): 65N30, 35L15, 78‐08. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

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     

On the difference approximations of overdetermined hyperbolic equations of classical mathematical physics

Issues concerning difference approximations of overdetermined systems of hyperbolic equations are examined. The formulations of extended overdetermined systems are given for hydrodynamics equations, magnetohydrodynamics equations, Maxwell equations, and elasticity equations. Some approaches to the construction of difference schemes are discussed for these systems.     

Hodge–Helmholtz decompositions of weighted Sobolev spaces in irregular exterior domains with inhomogeneous and anisotropic media

We study in detail Hodge–Helmholtz decompositions in nonsmooth exterior domains Ω??^N filled with inhomogeneous and anisotropic media. We show decompositions of alternating differential forms of rank q belonging to the weighted L²‐space L_s^{2, q}(Ω), s∈?, into irrotational and solenoidal q‐forms. These decompositions are essential tools, for example, in electro‐magnetic theory for exterior domains. To the best of our knowledge, these decompositions in exterior domains with nonsmooth boundaries and inhomogeneous and anisotropic media are fully new results. In the Appendix, we translate our results to the classical framework of vector analysis N=3 and q=1, 2. Copyright © 2008 John Wiley & Sons, Ltd.     

On electromagnetic waves in complex linear media in nonsmooth domains

Electromagnetic wave propagation in more complex linear materials such as bi‐anisotropic media have come to a considerable attention within the last 15–20 years. We shall propose a general framework to approach a class of highly complex materials. Such problems have been extensively studied mostly in the time‐harmonic case. In this paper, we focus on the time‐dependent case. A well‐posedness result for a large class of media is obtained. We also analyze Drude–Born–Fedorov type media in nonsmooth domains of arbitrary topological genus within this framework. Copyright © 2012 John Wiley & Sons, Ltd.     

Solving an initial value problem in inhomogeneous electrically and magnetically anisotropic uniaxial media

An initial value problem for a system of the second order partial differential equations, describing the electric wave propagation in vertically inhomogeneous electrically and magnetically anisotropic uniaxial media, is the main object of the study. The present paper suggests and justifies a new algorithm for solving this problem. This algorithm has several steps. On the first step the original initial value problem is written in terms of the Fourier images with respect to lateral variables. After that the obtained problem is transformed into an equivalent second kind vector integral equation of the Volterra type. A solution of this integral equation is constructed by successive approximations. At last, using the real Paley-Wiener theorem, a solution of the original initial value problem is found.     

A pseudo active kinematic constraint for a biological living soft tissue: An effect of the collagen network

Recent studies in mammalian hearts show that left ventricular wall thickening is an important mechanism for systolic ejection, and that during contraction the cardiac muscle develops significant stresses in the muscular cross-fiber direction. We suggested that the collagen network surrounding the muscular fibers could account for these mechanical behaviors. To test this hypothesis we develop a model for large deformation response of active, incompressible, nonlinear elastic and transversely isotropic living soft tissue (such as cardiac or arteries tissues) in which we include a coupling effect between the connective tissue and the muscular fibers. Then, a three-dimensional finite element formulation including this internal pseudo-active kinematic constraint is derived. Analytical and finite element solutions are in a very good agreement. The numerical results show this wall thickening effect with an order of magnitude compatible with the experimental observations.     

Reduced Finite Element Discretizations of the Stokes and Navier-Stokes Equations

ABSTRACT

If finite element spaces for the velocity and pressure do not satisfy the Babu?ka-Brezzi condition, a stable conforming discretization of the Stokes or Navier-Stokes equations can be obtained by enriching the velocity space by suitable functions. Writing any function from the enriched space as a sum of a function from the original space and a function from the supplementary space, the discretization will contain a number of additional terms compared with a conforming discretization for the original pair of spaces. We show that not all these terms are necessary for the solvability of the discrete problem and for optimal convergence properties of the discrete solutions, which is useful for saving computer memory and for establishing a connection to stabilized methods.     

On the convergence of finite element solutions to the interface problem for the Stokes system

The Stokes system with a discontinuous coefficient (Stokes interface problem) and its finite element approximations are considered. We firstly show a general error estimate. To derive explicit convergence rates, we introduce some appropriate assumptions on the regularity of exact solutions and on a geometric condition for the triangulation. We mainly deal with the MINI element approximation and then consider P1-iso-P2/P1 element approximation. Results are expected to give an instructive remark in numerical analysis for two-phase flow problems.     

Variational formulations for transmission problems of the electromagnetic field

The subject of this article is a review of all possible transmission problems for electromagnetic phenomena. In particular, we study the case of a perfect dielectric and a perfect conductor via a (formal) limit with conductivity approaching zero or infinity and discuss the expected regularity of the involved unknowns. Finally, we formulate equivalent variational formulations for each considered problem.     

Entire solutions of sublinear elliptic equations in anisotropic media

We study the nonlinear elliptic problem −Δu=ρ(x)f(u) in RN (N?3), lim|x|→∞u(x)=?, where ??0 is a real number, ρ(x) is a nonnegative potential belonging to a certain Kato class, and f(u) has a sublinear growth. We distinguish the cases ?>0 and ?=0 and prove existence and uniqueness results if the potential ρ(x) decays fast enough at infinity. Our arguments rely on comparison techniques and on a theorem of Brezis and Oswald for sublinear elliptic equations.     

FINITE ELEMENT SIMULATIONS FOR COMPRESSIBLE MISCIBLE DISPLACEMENT WITH MOLECULAR DISPERSION IN POROUS MEDIA

Abstract. We consider a nonlinear parabolic system describing compressible miscible displacement in a porous medium [5]. Continuous time and discrete time Galerkinmethods are introduced to approximate the solution snd optimal H1 error estimatesare obtained. One contribution of this paper is a demonstration of how moleculaxdispersion can be handled.     

Regularity of the Maxwell equations in heterogeneous media and Lipschitz domains

This note establishes regularity estimates for the solution of the Maxwell equations in Lipschitz domains with non-smooth coefficients and minimal regularity assumptions. The argumentation relies on elliptic regularity estimates for the Poisson problem with non-smooth coefficients.     

Global existence of weak solutions for the anisotropic compressible Stokes system

In this paper, we study the problem of global existence of weak solutions for the quasi-stationary compressible Stokes equations with an anisotropic viscous tensor. The key idea is a new identity that we obtain by comparing the limit of the equations of the energies associated to a sequence of weak-solutions with the energy equation associated to the system verified by the limit of the sequence of weak-solutions. In the context of stability of weak solutions, this allows us to construct a defect measure which is used to prove compactness for the density and therefore allowing us to identify the pressure in the limiting model. By doing so we avoid the use of the so-called effective flux. Using this new tool, we solve an open problem namely global existence of solutions à la Leray for such a system without assuming any restriction on the anisotropy amplitude. This provides a flexible and natural method to treat compressible quasilinear Stokes systems which are important for instance in biology, porous media, supra-conductivity or other applications in the low Reynolds number regime.     

Augmented high order finite volume element method for elliptic PDEs in non-smooth domains: Convergence study

The accuracy of a finite element numerical approximation of the solution of a partial differential equation can be spoiled significantly by singularities. This phenomenon is especially critical for high order methods. In this paper, we show that, if the PDE is linear and the singular basis functions are homogeneous solutions of the PDE, the augmentation of the trial function space for the Finite Volume Element Method (FVEM) can be done significantly simpler than for the Finite Element Method. When the trial function space is augmented for the FVEM, all the entries in the matrix originating from the singular basis functions in the discrete form of the PDE are zero, and the singular basis functions only appear in the boundary conditions. That is to say, there is no need to integrate the singular basis functions over the elements and the sparsity of the matrix is preserved without special care. FVEM numerical convergence studies on two-dimensional triangular grids are presented using basis functions of arbitrary high order, confirming the same order of convergence for singular solutions as for smooth solutions.