A high-order non-conforming discontinuous Galerkin method for time-domain electromagnetics
|
|
- Lawrence Elliott
- 6 years ago
- Views:
Transcription
1 A high-order non-conforming discontinuous Galerkin method for time-domain electromagnetics Hassan Fahs, Stephane Lanteri To cite this version: Hassan Fahs, Stephane Lanteri. A high-order non-conforming discontinuous Galerkin method for time-domain electromagnetics. Journal of Computational and Applied Mathematics, Elsevier, 010, 34 (4), pp < /j.cam >. <hal > HAL Id: hal Submitted on 10 Jul 013 HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers. L archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.
2 A high-order non-conforming discontinuous Galerkin method for time-domain electromagnetics Hassan Fahs a,, Stéphane Lanteri a a INRIA, 004 Route des Lucioles, BP 93, F-0690 Sophia Antipolis Cedex, France Abstract In this paper, we discuss the formulation, stability and validation of a high-order non-dissipative discontinuous Galerkin (DG) method for solving Maxwell s equations on non-conforming simplex meshes. The proposed method combines a centered approximation for the numerical fluxes at inter element boundaries, with either a second-order or a fourth-order leap-frog time integration scheme. Moreover, the interpolation degree is defined at the element level and the mesh is refined locally in a non-conforming way resulting in arbitrary-level hanging nodes. The method is proved to be stable and conserves a discrete counterpart of the electromagnetic energy for metallic cavities. Numerical experiments with high-order elements show the potential of the method. Key words: computational electromagnetism, time-domain Maxwell s equations, discontinuous Galerkin method, explicit time integration, non-conforming meshes 1 Introduction Time domain solutions of Maxwell s equations find applications in the applied sciences and engineering problems such as the design and optimization of antennas and radars, the design of emerging technologies (high speed electronics, integrated optics, etc.), the study of human exposure to electromagnetic waves [8], to name a few. These problems require high fidelity approximate solutions with a rigorous control of the numerical errors. Even for linear problems such conditions force one to look beyond standard computational techniques and Corresponding author. address: Hassan.Fahs@gmail.com (Hassan Fahs). Article published in J. Comput. Appl. Math. 34 (010)
3 seek new numerical frameworks enabling the accurate, efficient, and robust modeling of wave phenomena over long simulation times in settings of realistic geometrical complexity. The finite difference time-domain (FDTD) method, first introduced by Yee in 1966 [13] and later developed by Taflove and others [10], has been used for a broad range of applications in computational electromagnetics. In spite of its flexibility and second-order accuracy in a homogeneous medium, the Yee scheme suffers from serious accuracy degradation when used to model complex geometries. In recent years, a number of efforts aimed at addressing the shortcomings of the classical FDTD scheme, e.g. embedding schemes to overcome staircasing [1], high-order finite difference schemes [10]-[14], nonconforming orthogonal FDTD methods []. Most of these methods, however, have not really penetrated into main stream user community, partly due to their complicated nature and partly because these methods themselves often introduce other complications. The discontinuous Galerkin methods enjoy an impressive favor nowadays and are now used in various applications. Being higher order versions of traditional finite volume methods[6], discontinuous Galerkin time-domain(dgtd) methods based on discontinuous finite element spaces, easily handle elements of various types and shapes, irregular non-conforming meshes [4], and even locally varying polynomial degree. They hence offer great flexibility in the mesh design, but also lead to (block-) diagonal mass matrices and therefore yield fully explicit, inherently parallel methods when coupled with explicit time stepping [1]. Moreover, continuity is weakly enforced across mesh interfaces by adding suitable bilinear forms (so-called numerical fluxes) to the standard variational formulations. Whereas high-order discontinuous Galerkin time-domain methods have been developed on conforming hexahedral [3] and tetrahedral [5] meshes, the design of non-conforming discontinuous Galerkin time-domain methods is still in its infancy. In practice, the non-conformity can result from a local refinement of the mesh (i.e. h-refinement), of the interpolation degree (i.e. p-enrichment) or of both of them (i.e. hp-refinement). In this paper, we present a high-order DGTD method on non-conforming simplicial meshes. It is an extension of the DG formulation recently studied in [4]. One of the most important properties which should be aimed at is the conservation of a discrete counterpart of the electromagnetic energy on a general non-conforming simplex mesh with arbitrary level hanging nodes, including hp-type refinement. This cannot be obtained with DG methods based on upwind fluxes [7]. The rest of the paper is organized as follows. In section, we introduce the high-order non-conforming DGTD method for solving the firstorder Maxwell equations, based on totally centered fluxes and a high-order leap-frog time integration scheme. We prove the stability of the resulting fully discretized scheme and its energy conservation properties in section 3. The stability result is more general than the ones obtained in [4]-[5]. Numerical
4 results are presented in section 4. Finally, section 5 concludes this paper and states future research directions. Discontinuous Galerkin Time-Domain method We consider the Maxwell equations in three space dimensions for heterogeneous anisotropic linear media with no source. The electric permittivity tensor ǫ(x) and the magnetic permeability tensor µ(x) are varying in space, time-invariant and both symmetric positive definite. The electric field E = t (E x,e y,e z ) and the magnetic field H = t (H x,h y,h z ) verify: ǫ t E = curl H, µ t H = curl E, (1) where the symbol t denotes a time derivative. These equations are set and solved on a bounded polyhedral domain Ω of R 3. For the sake of simplicity, a metallic boundary condition is set everywhere on the domain boundary Ω, i.e. n E = 0 (where n denotes the unitary outwards normal). We consider a partition Ω h of Ω into a set of tetrahedra τ i of size h i = diam(τ i ) with boundaries τ i such that h = max τi Ω h h i. To each τ i Ω h we assign an integer p i 0 (the local interpolation order) and we collect the p i in the vector p = {p i : τ i Ω h }. Of course, if p i is uniform in all element τ i of the mesh, we have p = p i. Within this construction we admit meshes with possibly hanging nodes i.e. by allowing non-conforming (or irregular) meshes where element vertices can lie in the interior of faces of other elements. Each tetrahedron τ i is assumed to be the image, under a smooth bijective(diffeomorphic) mapping, of a fixed reference tetrahedron ˆτ = {ˆx,ŷ,ẑ ˆx,ŷ,ẑ 0;ˆx+ŷ + ẑ 1}. For each τ i, V i denotes its volume, and ǫ i and µ i are respectively the local electric permittivity and magnetic permeability tensors of the medium, which could be varying inside the element τ i. For two distinct tetrahedra τ i and τ k in Ω h, theintersection τ i τ k is atrianglea ik which we will call interface, withunitary normal vector n ik, oriented from τ i towards τ k. For the boundary interfaces, the index k corresponds to a fictitious element outside the domain. Finally, we denote by V i the set of indices of the elements which are neighbors of τ i (having an interface in common). We also define the perimeter P i of τ i by P i = k V i s ik. We have the following geometrical property for all elements: k V i s ik n ik = 0. In the following, for a given partition Ω h and vector p, we seek approximate solutions to (1) in the finite dimensional subspace V p (Ω h ) := { v L (Ω) 3 : v τi P pi (τ i ), τ i Ω h }, where P pi (τ i ) denotes thespace of nodal polynomials 3
5 of degree at most p i inside the element τ i. Note that the polynomial degree, p i, may vary from element to element in the mesh. By non-conforming interface we mean an interface a ik for which at least one of its vertices is a hanging node or/and such that p i aik p k aik. According to the discontinuous Galerkin approach, the electric and magnetic fields inside each finite element are linear combinations ( E i, H i ) of linearly independent basis vector fields ϕ ij, 1 j d i, where d i denotes the local number of degrees of freedom (DOF) inside τ i. We denote by P i = Span( ϕ ij, 1 j d i ). The approximate fields ( E h, H h ), defined by ( i, E h τi = E i, H h τi = H i ) are allowed to be completely discontinuous across element boundaries. For such a discontinuous field U h, we define its average { U h } ik through any internal interface a ik, as { U h } ik = ( U i aik + U k aik )/. Note that for any internal interface a ik, { U h } ki = { U h } ik. Because of this discontinuity, a global variational formulation cannot be obtained. However, dot-multiplying (1) by any given vector function ϕ P i, integrating over each single element τ i and integrating by parts, yields: ϕ ǫ i te = curl ϕ τi H ϕ ( H n), τ i τ i ϕ µ i th = curl ϕ τi E+ () ϕ ( E n). τ i τ i In equations (), we now replace the exact fields E and H by the approximate fields E h and H h in order to evaluate volume integrals. For integrals over τ i, a specific treatment must be introduced since the approximate fields are discontinuous through element faces. We choose to use completely centered fluxes,i.e. i, k V i, E aik { E h } ik, H aik { H h } ik.themetallicboundary condition on a boundary interface a ik (k in the element index of the fictitious neighboring element) is dealt with weakly, in the sense that traces of fictitious fields E k and H k are used for the computation of numerical fluxes for the boundary element τ i. In the present case, where all boundaries are metallic, wesimplytake E k aik = E i aik andh k aik = H i aik.replacingsurfaceintegrals using centered fluxes in () and re-integrating by parts yields: ϕ ǫ i tei = 1 (curl ϕ τ i H i +curlh i ϕ) τ i 1 ϕ ( H k n ik ), k V a ik i ϕ µ i thi = 1 (3) (curl ϕ τ i E i +curle i ϕ) τ i + 1 ϕ ( E k n ik ). a ik k V i 4
6 We can rewrite this formulation in terms of scalar unknowns. Inside each element, the fields are recomposed according to E i = 1 j d i E ij ϕ ij, H i = 1 j d i H ij ϕ ij. Let us denote by E i and H i respectively the column vectors (E il ) 1 l di and (H il ) 1 l di. Equations (3) can be rewritten as: M i ǫ te i = K i H i S ik H k, k V i M µ i t H i = K i E i + S ik E k, k V i (4) where the symmetric positive definite mass matrices Mi σ (σ stands for ǫ or µ), and the symmetric stiffness matrix K i (all of size d i ) are given by : (Mi σ ) jl = t ϕ ij σ i ϕ il and (K i ) jl = 1 t ϕ ij curl ϕ il + t ϕ il curl ϕ ij. For any interface τ i τ i a ik, the d i d k rectangular matrix S ik is given by: (S ik ) jl = 1 a ik t ϕ ij ( ϕ kl n ik ), 1 j d i, 1 l d k. (5) Concerning the time discretization, we propose to use a leap-frog (LF N, N =,4) scheme. This kind of time scheme has both advantages to be explicit and to be non-dissipative. In the sequel, superscripts refer to time stations and t is the fixed time-step. The unknowns related to the electric field are approximated at integer time-stations t n = n t and are denoted by E n i. The unknowns related to the magnetic field are approximated at half-integer timestations t n+1/ = (n+1/) t and are denoted by H n+1/ i. The LF N (N =,4) integrator is constructed as follows [15]-[9]: T 1 = t(mi) ǫ 1 curlh n+1 i, T 1 = t(m µ i ) 1 curle n+1 i, T = t(m µ i ) 1 curlt 1, T = t(mǫ i ) 1 curlt 1, T 3 = t(mi ǫ) 1 curlt, T 3 = t(mµ i ) 1 curlt. E n+1 i = E n i LF : +T 1, H n+3 i = H n+1 i +T 1. E n+1 i = E n i LF 4 : +T 1 +T 3 /4, H n+3 i = H n+1 i +T 1 +T 3/4. (6) For the treatment of the boundary condition on an interface a ik, we use: E n k a ik = E n i a ik and H n+1 k a ik = H n+1 i a ik. (7) 5
7 3 Stability of the discontinuous Galerkin method We aim at giving and proving a sufficient condition for the L -stability of the proposed discontinuous Galerkin method with only metallic boundary conditions. We use the same kind of energy approach as in [5], where a quadratic form plays the role of a Lyapunov function of the whole set of numerical unknowns. To this end, we suppose that all electric (resp. magnetic) unknowns are gathered in a column vector E (resp. H) of size d = id i, then the space discretized system (4) can be rewritten as: M ǫ t E = KH AH BH, M µ t H = KE+AE BE, (8) where we have the following definitions and properties: M ǫ,m µ and K are d d block diagonal matrices with diagonal blocks equal tom ǫ i,mµ i andk i respectively. ThereforeM ǫ andm µ aresymmetricpositive definite matrices, and K is a symmetric matrix. A is also a d d block sparse matrix, whose non-zero blocks are equal to S ik when a ik is an internal interface of the mesh. Since n ki = n ik, it can be checked from (5) that (S ik ) jl = (S ki ) lj and then S ki = t S ik ; thus A is a symmetric matrix. B is a d d block diagonal matrix, whose non-zero blocks are equal to S ik when a ik is a metallic boundary interface of the mesh. In that case, (S ik ) jl = (S ik ) lj, and S ik = t S ik ; thus B is a skew-symmetric matrix. ThediscontinuousGalerkinDGTD-P pi methodusingcenteredfluxescombined with Nth order leap-frog (LF N ) time scheme and arbitrary local accuracy and basis functions can be written, in function of the matrix S = K A B, in the general form: M ǫen+1 E n t M µhn+3 H n+1 t = S N H n+1, = t S N E n+1, (9) where the matrix S N (N being the order of the leap-frog scheme) verifies: S if N =, S N = S(I t 4 M µ t SM ǫ S) if N = 4. (10) 6
8 We now define the following discrete version of the electromagnetic energy. Definition 1 We consider the following electromagnetic energies inside each tetrahedron τ i and in the whole domain Ω: the local energy : i, E n i = 1 ( t E n i Mǫ i En i + t H n 1 i M µ i H n+1 i ), (11) the global energy : E n = 1 ( t E n M ǫ E n + t H n 1 M µ H n+1 ). (1) In the following, we shall prove that the global energy(1) is conserved through a time step and that it is a positive definite quadratic form of all unknowns under a CFL-like condition on the time-step t. Lemma 1 Using the DGTD-P pi method (9)-(10) for solving (1) with metallic boundaries only, the global discrete energy (1) is exactly conserved, i.e. E n+1 E n = 0, n. Proof. We denote by E n+1 = En+1 +E n. We have : E n+1 E n = t E n+1 M ǫ (E n+1 E n )+ 1 t H n+1 M µ (H n+3 H n 1 ) = t t 1 E n+1 SN H n+1 t t H n+ 1 ( t S N E n+1 + t S N E n ) = t t H n+1 ( t S N t S N )E n+1 = 0. This concludes the proof. Lemma Using the DGTD-P pi method (9)-(10), the global discrete electromagnetic energy E n (1) is a positive definite quadratic form of all unknowns if: t, with d N = M µ t S N M ǫ, (13) d N where. denotes a matrix norm, and the matrix M σ is the inverse square root of M σ. Also, for a given mesh, the stability limit of the LF 4 scheme is roughly.85 times larger than that of the LF scheme. Proof. The mass matrices M ǫ and M µ are symmetric positive definite and we can construct in a simple way their square root (also symmetric positive definite) denoted by M ǫ and M µ respectively. Using the scheme (9) to develop H n+1 in function of E n and H n 1, yields: 7
9 E n = 1 t E n M ǫ E n + 1 t H n 1 M µ H n+1 = 1 t E n M ǫ E n + 1 t H n 1 M µ t H n 1 1 Mǫ E n + 1 Mµ H n 1 t t H n 1 t S N E n t H n 1 t S N M ǫ 1 Mǫ E n + 1 Mµ H n 1 d N t M µ H n 1 ǫ M E n. µ µ M M M ǫ E n Atthispoint,wechoosetouseanupperboundfortheterm M µ H n 1 M ǫ E n which might lead to sub-optimal lower bounds for the energy (and then to a slightly too severe stability limit for the scheme). Anyway, this stability limit is only sufficient, and not really close to necessary. We use the inequality: M µ H n 1 M ǫ E n 1 ( Mµ H n 1 + M ǫ E n ). We then sum up the lower bounds for the E n to obtain: E n 1 (1 d N t ) M ǫ E n + 1 (1 d N t ) M µ H n 1. Then, under the condition proposed in Lemma, the electromagnetic energy E n is a positive definite quadratic form of all unknowns. Moreover, for a given mesh, using the definition (10) of S N, the LF 4 scheme is stable if: t M µ t S 4 M ǫ, t M µ t (S t 4 S M µ t S M ǫ S )M ǫ, td t3 4 d3. This inequality is verified if and only if d t ( ) (.847). This concludes the proof. Now, our objective is to give an explicit CFL condition on t under which the local energy (11) is a positive definite quadratic form of the numerical unknowns E n i and H n 1 i. We first need some classical definitions. Definition We assume that the tensors ǫ i and µ i are piecewise constant, i.e. ǫ i = ǫ i and µ i = µ i. We denote by c i = 1/ ǫ i µ i the propagation speed in 8
10 the finite element τ i. We also assume that there exist dimensionless constants α i and β ik (k V i ) such that: X P i, curl X τi α ip i V i X τi, X a ik β iks ik V i X τ i, (14) where X τi and X aik denote the L -norm of the vector field X over τ i and the interface a ik respectively. Lemma 3 Using the LF scheme (4)-(6)-(7), under assumptions of Definition, the local discrete energy E n i (11) is a positive definite quadratic form of all unknowns (E n i,hn 1 i ) and the scheme is stable if the time step t is such that: i, k V i, c i t[α i +β ik ] < 4V i, (15) P i (with the convention that, in the above formula, k should be replaced by i for a metallic boundary interface a ik ). in the defi- Proof. Using the scheme (3) to replace the occurrences of H n+1 i nition of E i, and using the boundary fluxes given in (7), we get: Ei n = ǫ i En i τ i + µ 1 i Hn i τ i t 4 Xn i, with ( X n i = curlh n 1 i E n i +curl E n i H n 1 i ) k Vi τ i ( H n i a ik 1 E n k ) n ik. In the remainder of this proof, we omit the superscripts n and n-1/ respectively in the electric and magnetic variables. We have the following identities: X n i curl H i τi E i τi + curl E i τi H i τi + 1 k V i ( ) µi H ǫ i ǫi a ik + E i µ k a ik i α ip i V i H i τi E i τi + 1 k V i ( µi ǫ i β ik s ik V i H i τ i + ) ǫi β ki s ik E µ i V k τ k. k Noticing that H i τi E i τi c i (µ i H i τ i + ǫ i E i τ i ), gathering all lower bounds for terms in the expression of E n i and using P i = k V i s ik leads to: 9
11 Ei n ( 1 s ik α ic i t k V i P i 4V i t 8 k V i s ik ) (ǫ i E i τ i +µ i H i τ i ) ( µi β ik H ǫ i V i ǫi β ki τ i + E i µ i V k τ k k ). Then, summing upthese inequalities inorder toobtainalower boundfor ie i leads to an expression that we reorganize as sum over interfaces. We find that ie i a ik s ik W ik with: ( W ik = ǫ i E 1 i τ i α ic i t β ) ikc i t + P i 4V i 8V i ( µ i H 1 i τ i α ic i t β ) ikc i t + P i 4V i 8V i ( ǫ k E 1 k τ k α kc k t β ) kic k t + P k 4V k 8V k ( µ k H 1 k τ k α kc k t β ) kic k t. P k 4V k 8V k UndertheconditionsproposedinLemma3,W ik isapositivedefinitequadratic form of all unknowns and so is the local energy. This concludes the proof. Note that, the existence of the constants α i and β ik (k V i ) is always ensured. The values of α i only depend on the local polynomial order p i while the values of β ik depend on p i and on the number of hanging nodes on the interface a ik. Forinstance, fororthogonalpolynomials onad-simplex β ik = (p i +1)(p i +d)/d (see [11]), and for arbitrary basis functions these values are given by: ( α i P i ; β iks ik ) = ( M 1/ S Vi 1 M 1/ ; M 1/ S M 1/ ), V i where M is the mass matrix without material parameter, S = S ik, and S 1 = curl ϕ ij curl ϕ il, 1 j,l d i. τ i 4 Numerical experiments We consider here the Maxwell equations in two space dimensions and in the TM-polarization; i.e. we solve for (H x,h y,e z ). We validate the theory by considering the propagation of an eigenmode which is a standing wave of frequency f = 1 MHz and wavelength λ = 1.4 m in a unitary metallic cavity 10
12 with ǫ = µ = 1 in normalized units. Owing to the existence of an exact analytical solution, this problem allows us to appreciate the numerical results at any point and time in the cavity. Numerical simulations make use of triangular meshes of the square [0,1] [0,1] and a non-conforming mesh is obtained by locally refining (two refinement levels) the square zone [0.5, 0.75] [0.5, 0.75] of a coarse conforming mesh as shown on Fig. 1. The resulting non-conforming mesh consists of 78 triangles and 44 nodes (36 of them are hanging nodes). For this non-conforming mesh, we assign to coarse (i.e. non refined) elements a high polynomial degree p 1 and to the refined region a low polynomial degree p. The resulting scheme is referred to as DGTD-P (p1,p ). If p 1 = p = p, the scheme is simply called DGTD-P p. Note that, for a conforming interface a ik, Fig. 1. Non-conforming locally refined triangular mesh. the matrix S ik defined in (5) can be evaluated in a direct way once and for all. However, for a non-conforming interface, we cannot calculate this matrix with an exact formula because it depends on the number of hanging nodes on the interface a ik. For that, and only for non-conforming interfaces, we calculate the matrix S ik by using a Gaussian quadrature formula. All simulations are carried out for time t = 150 which corresponds to 106 periods. In Tab. 1, we summarize the CFL values of the LF DGTD-P p method. If p 1 p, the DGTD-P (p1,p ) method has the same stability limit as the DGTD-P min(p1,p ) method, as long as the mesh is actually refined. We plot on Fig. the time evolution of the overall L error of the DGTD-P p and DGTD-P (p1,p ) methods using the LF and LF 4 schemes. Tab. gives the L error, the number of degrees of freedom and the CPU time to reach time t = 150. It can be observed from Fig. that the gain in the L error is noticeable when the accuracy in space and time is increased. Moreover, it is clear from (6) and Lemma that, for the same non-conforming mesh, each time step of LF 4 requires times more memory than the LF time step, but its stability limit is almost.85 times less restrictive. Then, LF 4 requires almost 1.5 times less CPU time and is roughly 15 times more efficient than LF. Furthermore, for a given accuracy, the LF 4 DGTD-P (p1,p ) method requires less CPU time than the LF 4 DGTD-P p method. Fig. 3 illustrates the numerical convergence of the DGTD- 11
13 P p anddgtd-p (p1,p ) methods. Corresponding asymptotic convergence orders are summarized in Tab. 3. As it could be expected from the use of the Nth accurate time integration scheme, the asymptotic convergence order is bounded by N independently of the interpolation degree. Table 1 The CFL values of the LF DGTD methods. DGTD-P p method, p = CFL(LF ) DGTD-P (p1,p ) method, (p 1,p ) = (3,) (4,) (4,3) (5,3) (5,4) CFL(LF ) e-04 Error 1e-05 1e-06 DGTD-P, LF DGTD-P, LF4 DGTD-P3, LF DGTD-P3, LF4 DGTD-P4, LF DGTD-P4, LF4 1e Time e-04 Error 1e-05 1e-06 DGTD-P(3,), LF DGTD-P(3,), LF4 DGTD-P(4,3), LF DGTD-P(4,3), LF4 DGTD-P(5,4), LF DGTD-P(5,4), LF4 1e Fig.. Time evolution of the L error. DGTD-P p (top) and DGTD-P (p1,p ) (bottom) methods. Time 1
14 Table #DOF, L errorsandcputimeinminutesusingthelf andlf 4 DGTDmethods. DGTD-P p method LF LF 4 p # DOF Error CPU (min) Error CPU (min) E E E E E E E E DGTD-P (p1,p ) method LF LF 4 (p 1,p ) # DOF Error CPU (min) Error CPU (min) (3,) E E-05 1 (4,) E E (4,3) E E (5,4) E E DGTD P p method with LF time scheme DGTD P, LF DGTD P3, LF DGTD P4, LF DGTD P p method with LF4 time scheme DGTD P, LF4 DGTD P3, LF4 DGTD P4, LF L error 10 4 L error (DOF) 1/ (DOF) 1/ 10 DGTD P(p1,p) method with LF time scheme DGTD P(3,), LF DGTD P(4,), LF DGTD P(4,3), LF DGTD P(5,3), LF DGTD P(5,4), LF DGTD P(p1,p) method with LF4 time scheme DGTD P(3,), LF4 DGTD P(4,), LF4 DGTD P(4,3), LF4 DGTD P(5,3), LF4 DGTD P(5,4), LF L error 10 4 L error (DOF) 1/ (DOF) 1/ Fig. 3. Numerical convergence of the DGTD-P p and DGTD-P (p1,p ) methods. 13
15 Table 3 Asymptotic convergence orders of the LF and LF 4 DGTD methods. DGTD-P p method, p = 3 4 LF scheme LF 4 scheme DGTD-P (p1,p ) method, (p 1,p ) = (3,) (4,) (4,3) (5,3) (5,4) LF scheme LF 4 scheme Concluding remarks In this paper, we have studied a high-order discontinuous Galerkin method for the discretization of the time-domain Maxwell equations on non-conforming simplicial meshes. We proved that the method conserves a discrete equivalent of the electromagnetic energy and it is stable under some CFL-type stability condition. Numerical simulations were performed by considering an eigenmode problem in two space dimensions. We have shown that, for a given non-conforming mesh, the DGTD methods coupled to the LF 4 scheme is at least 15 times more accurate and requires roughly 1.5 times less CPU time than the LF DGTD methods. Concerning future works, our objective is to design a truly hp-adaptive method through the construction of an appropriate error estimator. References [1] M. Bernacki, L. Fezoui, S. Lanteri, and S. Piperno. Parallel unstructured mesh solvers for heterogeneous wave propagation problems. Appl. Math. Model., 30: , 006. [] F. Collino, T. Fouquet, and P. Joly. Conservative space-time mesh refinement methods for the FDTD solution of Maxwell s equations. J. Comput. Phys., 11:9 35, 006. [3] G. Cohen, X. Ferrieres, and S. Pernet. A spatial high-order hexahedral discontinuous Galerkin method to solve Maxwell s equations in time domain. J. Comput. Phys., 17: , 006. [4] H. Fahs. Development of a hp-like discontinuous Galerkin time-domain method on non-conforming simplicial meshes for electromagnetic wave propagation. Int. J. Numer. Anal. Model., 6:193 16, 009. [5] L. Fezoui, S. Lanteri, S. Lohrengel, and S. Piperno. Convergence and stability of a discontinuous Galerkin time-domain method for the heterogeneous Maxwell 14
16 equations on unstructured meshes. ESAIM: Math. Model. and Numer. Anal., 39: , 005. [6] J.S. Hesthaven and T. Warburton. Nodal high-order methods on unstructured grids. I. Time-domain solution of Maxwell s equations. J. Comput. Phys., 181:186 1, 00. [7] D. Kopriva, S.L. Woodruff, and M.Y. Hussaini. Discontinuous spectral element approximation of Maxwell s equations. In B. Cockburn, G.E. Karniadakis, and C.W. Shu, editors, Discontinuous Galerkin Methods: Theory, Computation and Applications, vol. 11, Springer-Verlag, 000, pp [8] G. Scarella, O. Clatz, S. Lanteri, G. Beaume, S. Oudot, J.-P. Pons, S. Piperno, P. Joly, and J. Wiart. Realistic numerical modelling of human head tissue exposure to electromagnetic waves from cellular phones. Comptes Rendus Physique, 7: , 006. [9] H. Spachmann, R. Schuhmann, and T. Weiland. High order explicit time integration schemes for Maxwell s equations. Int. J. Numer. Model., 15: , 00. [10] A. Taflove. Advances in computational electrodynamics, the finite-difference time-domain method. Artech House, Boston, London, [11] T. Warburton and J.S. Hesthaven. On the constants in hp-finite element trace inverse inequalities. Comput. Methods Appl. Mech. Engrg., 19: , 003. [1] T. Xiao and Q.H. Liu. A staggered upwind embedded boundary (SUEB) method to eliminate the FDTD staircasing error. IEEE Trans. Antennas and Propagat., 5: , 004. [13] K.S. Yee. Numerical solution of initial boundary value problems involving Maxwell s equations in isotropic media. IEEE Trans. Antennas and Propagat., 14:30 307, [14] A. Yefet and P.G. Petropoulos. A staggered fourth-order accurate explicit finite difference scheme for the time-domain Maxwell s equations. J. Comput. Phys., 168:86 315, 001. [15] J.L. Young. High-order, leapfrog methodology for the temporally dependent Maxwell s equations. Radio Science, 36:9 17,
Approximation SEM-DG pour les problèmes d ondes elasto-acoustiques
Approximation SEM-DG pour les problèmes d ondes elasto-acoustiques Helene Barucq, Henri Calandra, Aurélien Citrain, Julien Diaz, Christian Gout To cite this version: Helene Barucq, Henri Calandra, Aurélien
More informationNodal and divergence-conforming boundary-element methods applied to electromagnetic scattering problems
Nodal and divergence-conforming boundary-element methods applied to electromagnetic scattering problems M. Afonso, Joao Vasconcelos, Renato Mesquita, Christian Vollaire, Laurent Nicolas To cite this version:
More informationFinite volume method for nonlinear transmission problems
Finite volume method for nonlinear transmission problems Franck Boyer, Florence Hubert To cite this version: Franck Boyer, Florence Hubert. Finite volume method for nonlinear transmission problems. Proceedings
More informationReduced Vlasov-Maxwell modeling
Reduced Vlasov-Maxwell modeling Philippe Helluy, Michel Massaro, Laurent Navoret, Nhung Pham, Thomas Strub To cite this version: Philippe Helluy, Michel Massaro, Laurent Navoret, Nhung Pham, Thomas Strub.
More informationTemporal convergence analysis of a locally implicit discontinuous galerkin time domain method for electromagnetic wave propagation in dispersive media
Temporal convergence analysis of a locally implicit discontinuous galerkin time domain method for electromagnetic wave propagation in dispersive media Stéphane Descombes, Stéphane Lanteri, Ludovic Moya
More informationOptimized Schwarz Methods for Maxwell Equations with Discontinuous Coefficients
Optimized Schwarz Methods for Maxwell Equations with Discontinuous Coefficients Victorita Dolean, Martin Gander, Erwin Veneros To cite this version: Victorita Dolean, Martin Gander, Erwin Veneros. Optimized
More informationPalindromic Discontinuous Galerkin Method
Palindromic Discontinuous Galerkin Method David Coulette, Emmanuel Franck, Philippe Helluy, Michel Mehrenberger, Laurent Navoret To cite this version: David Coulette, Emmanuel Franck, Philippe Helluy,
More informationA FV Scheme for Maxwell s equations
A FV Scheme for Maxwell s equations Convergence Analysis on unstructured meshes Stephanie Lohrengel * Malika Remaki ** *Laboratoire J.A. Dieudonné (UMR CNRS 6621), Université de Nice Sophia Antipolis,
More informationA non-commutative algorithm for multiplying (7 7) matrices using 250 multiplications
A non-commutative algorithm for multiplying (7 7) matrices using 250 multiplications Alexandre Sedoglavic To cite this version: Alexandre Sedoglavic. A non-commutative algorithm for multiplying (7 7) matrices
More informationAntipodal radiation pattern of a patch antenna combined with superstrate using transformation electromagnetics
Antipodal radiation pattern of a patch antenna combined with superstrate using transformation electromagnetics Mark Clemente Arenas, Anne-Claire Lepage, Xavier Begaud To cite this version: Mark Clemente
More informationThe magnetic field diffusion equation including dynamic, hysteresis: A linear formulation of the problem
The magnetic field diffusion equation including dynamic, hysteresis: A linear formulation of the problem Marie-Ange Raulet, Benjamin Ducharne, Jean-Pierre Masson, G. Bayada To cite this version: Marie-Ange
More informationCase report on the article Water nanoelectrolysis: A simple model, Journal of Applied Physics (2017) 122,
Case report on the article Water nanoelectrolysis: A simple model, Journal of Applied Physics (2017) 122, 244902 Juan Olives, Zoubida Hammadi, Roger Morin, Laurent Lapena To cite this version: Juan Olives,
More informationFull-order observers for linear systems with unknown inputs
Full-order observers for linear systems with unknown inputs Mohamed Darouach, Michel Zasadzinski, Shi Jie Xu To cite this version: Mohamed Darouach, Michel Zasadzinski, Shi Jie Xu. Full-order observers
More informationOn Symmetric Norm Inequalities And Hermitian Block-Matrices
On Symmetric Norm Inequalities And Hermitian lock-matrices Antoine Mhanna To cite this version: Antoine Mhanna On Symmetric Norm Inequalities And Hermitian lock-matrices 016 HAL Id: hal-0131860
More informationInteractions of an eddy current sensor and a multilayered structure
Interactions of an eddy current sensor and a multilayered structure Thanh Long Cung, Pierre-Yves Joubert, Eric Vourc H, Pascal Larzabal To cite this version: Thanh Long Cung, Pierre-Yves Joubert, Eric
More informationOn Symmetric Norm Inequalities And Hermitian Block-Matrices
On Symmetric Norm Inequalities And Hermitian lock-matrices Antoine Mhanna To cite this version: Antoine Mhanna On Symmetric Norm Inequalities And Hermitian lock-matrices 015 HAL Id: hal-0131860
More informationThe sound power output of a monopole source in a cylindrical pipe containing area discontinuities
The sound power output of a monopole source in a cylindrical pipe containing area discontinuities Wenbo Duan, Ray Kirby To cite this version: Wenbo Duan, Ray Kirby. The sound power output of a monopole
More informationA new simple recursive algorithm for finding prime numbers using Rosser s theorem
A new simple recursive algorithm for finding prime numbers using Rosser s theorem Rédoane Daoudi To cite this version: Rédoane Daoudi. A new simple recursive algorithm for finding prime numbers using Rosser
More informationInfluence of a Rough Thin Layer on the Potential
Influence of a Rough Thin Layer on the Potential Ionel Ciuperca, Ronan Perrussel, Clair Poignard To cite this version: Ionel Ciuperca, Ronan Perrussel, Clair Poignard. Influence of a Rough Thin Layer on
More informationOn sl3 KZ equations and W3 null-vector equations
On sl3 KZ equations and W3 null-vector equations Sylvain Ribault To cite this version: Sylvain Ribault. On sl3 KZ equations and W3 null-vector equations. Conformal Field Theory, Integrable Models, and
More informationNonlocal computational methods applied to composites structures
Nonlocal computational methods applied to composites structures Norbert Germain, Frédéric Feyel, Jacques Besson To cite this version: Norbert Germain, Frédéric Feyel, Jacques Besson. Nonlocal computational
More informationOn the longest path in a recursively partitionable graph
On the longest path in a recursively partitionable graph Julien Bensmail To cite this version: Julien Bensmail. On the longest path in a recursively partitionable graph. 2012. HAL Id:
More informationA non-commutative algorithm for multiplying (7 7) matrices using 250 multiplications
A non-commutative algorithm for multiplying (7 7) matrices using 250 multiplications Alexandre Sedoglavic To cite this version: Alexandre Sedoglavic. A non-commutative algorithm for multiplying (7 7) matrices
More informationON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS
ON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS Abdelhafid Younsi To cite this version: Abdelhafid Younsi. ON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS. 4 pages. 212. HAL Id:
More informationA proximal approach to the inversion of ill-conditioned matrices
A proximal approach to the inversion of ill-conditioned matrices Pierre Maréchal, Aude Rondepierre To cite this version: Pierre Maréchal, Aude Rondepierre. A proximal approach to the inversion of ill-conditioned
More informationOn Newton-Raphson iteration for multiplicative inverses modulo prime powers
On Newton-Raphson iteration for multiplicative inverses modulo prime powers Jean-Guillaume Dumas To cite this version: Jean-Guillaume Dumas. On Newton-Raphson iteration for multiplicative inverses modulo
More informationA Simple Proof of P versus NP
A Simple Proof of P versus NP Frank Vega To cite this version: Frank Vega. A Simple Proof of P versus NP. 2016. HAL Id: hal-01281254 https://hal.archives-ouvertes.fr/hal-01281254 Submitted
More informationSolution to Sylvester equation associated to linear descriptor systems
Solution to Sylvester equation associated to linear descriptor systems Mohamed Darouach To cite this version: Mohamed Darouach. Solution to Sylvester equation associated to linear descriptor systems. Systems
More informationAll Associated Stirling Numbers are Arithmetical Triangles
All Associated Stirling Numbers are Arithmetical Triangles Khaled Ben Letaïef To cite this version: Khaled Ben Letaïef. All Associated Stirling Numbers are Arithmetical Triangles. 2017.
More informationOn The Exact Solution of Newell-Whitehead-Segel Equation Using the Homotopy Perturbation Method
On The Exact Solution of Newell-Whitehead-Segel Equation Using the Homotopy Perturbation Method S. Salman Nourazar, Mohsen Soori, Akbar Nazari-Golshan To cite this version: S. Salman Nourazar, Mohsen Soori,
More informationNumerical modification of atmospheric models to include the feedback of oceanic currents on air-sea fluxes in ocean-atmosphere coupled models
Numerical modification of atmospheric models to include the feedback of oceanic currents on air-sea fluxes in ocean-atmosphere coupled models Florian Lemarié To cite this version: Florian Lemarié. Numerical
More informationRHEOLOGICAL INTERPRETATION OF RAYLEIGH DAMPING
RHEOLOGICAL INTERPRETATION OF RAYLEIGH DAMPING Jean-François Semblat To cite this version: Jean-François Semblat. RHEOLOGICAL INTERPRETATION OF RAYLEIGH DAMPING. Journal of Sound and Vibration, Elsevier,
More informationEddy-Current Effects in Circuit Breakers During Arc Displacement Phase
Eddy-Current Effects in Circuit Breakers During Arc Displacement Phase Olivier Chadebec, Gerard Meunier, V. Mazauric, Yann Le Floch, Patrice Labie To cite this version: Olivier Chadebec, Gerard Meunier,
More informationOptical component modelling and circuit simulation using SERENADE suite
Optical component modelling and circuit simulation using SERENADE suite Laurent Guilloton, Smail Tedjini, Tan-Phu Vuong To cite this version: Laurent Guilloton, Smail Tedjini, Tan-Phu Vuong. Optical component
More informationFinite Difference Solution of Maxwell s Equations
Chapter 1 Finite Difference Solution of Maxwell s Equations 1.1 Maxwell s Equations The principles of electromagnetism have been deduced from experimental observations. These principles are Faraday s law,
More informationPseudospectral and High-Order Time-Domain Forward Solvers
Pseudospectral and High-Order Time-Domain Forward Solvers Qing H. Liu G. Zhao, T. Xiao, Y. Zeng Department of Electrical and Computer Engineering Duke University DARPA/ARO MURI Review, August 15, 2003
More informationSmart Bolometer: Toward Monolithic Bolometer with Smart Functions
Smart Bolometer: Toward Monolithic Bolometer with Smart Functions Matthieu Denoual, Gilles Allègre, Patrick Attia, Olivier De Sagazan To cite this version: Matthieu Denoual, Gilles Allègre, Patrick Attia,
More informationMethylation-associated PHOX2B gene silencing is a rare event in human neuroblastoma.
Methylation-associated PHOX2B gene silencing is a rare event in human neuroblastoma. Loïc De Pontual, Delphine Trochet, Franck Bourdeaut, Sophie Thomas, Heather Etchevers, Agnes Chompret, Véronique Minard,
More informationOn size, radius and minimum degree
On size, radius and minimum degree Simon Mukwembi To cite this version: Simon Mukwembi. On size, radius and minimum degree. Discrete Mathematics and Theoretical Computer Science, DMTCS, 2014, Vol. 16 no.
More informationJoseph Charles. INRIA Sophia Antipolis-Méditerranée Nachos project-team
Amélioration des performances de méthodes Galerkin discontinues d ordre élevé pour la résolution numérique des équations de Maxwell instationnaires sur des maillages simplexes Joseph Charles INRIA Sophia
More informationA New Integral Formulation for Eddy Current Computation in Thin Conductive Shells
A New Integral Formulation for Eddy Current Computation in Thin Conductive Shells Tung Le Duc, Gérard Meunier, Olivier Chadebec, Jean-Michel Guichon To cite this version: Tung Le Duc, Gérard Meunier, Olivier
More informationClassification of high dimensional data: High Dimensional Discriminant Analysis
Classification of high dimensional data: High Dimensional Discriminant Analysis Charles Bouveyron, Stephane Girard, Cordelia Schmid To cite this version: Charles Bouveyron, Stephane Girard, Cordelia Schmid.
More informationSome approaches to modeling of the effective properties for thermoelastic composites
Some approaches to modeling of the ective properties for thermoelastic composites Anna Nasedkina Andrey Nasedkin Vladimir Remizov To cite this version: Anna Nasedkina Andrey Nasedkin Vladimir Remizov.
More informationFORMAL TREATMENT OF RADIATION FIELD FLUCTUATIONS IN VACUUM
FORMAL TREATMENT OF RADIATION FIELD FLUCTUATIONS IN VACUUM Frederic Schuller, Renaud Savalle, Michael Neumann-Spallart To cite this version: Frederic Schuller, Renaud Savalle, Michael Neumann-Spallart.
More informationOverview of FFT-based homogenization techniques from the Galerkin point of view (slides)
Overview of FFT-based homogenization techniques from the Galerkin point of view (slides) Sébastien Brisard To cite this version: Sébastien Brisard. Overview of FFT-based homogenization techniques from
More informationA Context free language associated with interval maps
A Context free language associated with interval maps M Archana, V Kannan To cite this version: M Archana, V Kannan. A Context free language associated with interval maps. Discrete Mathematics and Theoretical
More informationInductive thermography nondestructive testing applied to carbon composite materials: multiphysics and multiscale modeling
Inductive thermography nondestructive testing applied to carbon composite materials: multiphysics and multiscale modeling H.K. Bui, G Berthiau, D Trichet, G Wasselynck To cite this version: H.K. Bui, G
More informationQuasi-periodic solutions of the 2D Euler equation
Quasi-periodic solutions of the 2D Euler equation Nicolas Crouseilles, Erwan Faou To cite this version: Nicolas Crouseilles, Erwan Faou. Quasi-periodic solutions of the 2D Euler equation. Asymptotic Analysis,
More informationQuantum efficiency and metastable lifetime measurements in ruby ( Cr 3+ : Al2O3) via lock-in rate-window photothermal radiometry
Quantum efficiency and metastable lifetime measurements in ruby ( Cr 3+ : Al2O3) via lock-in rate-window photothermal radiometry A. Mandelis, Z. Chen, R. Bleiss To cite this version: A. Mandelis, Z. Chen,
More informationThermodynamic form of the equation of motion for perfect fluids of grade n
Thermodynamic form of the equation of motion for perfect fluids of grade n Henri Gouin To cite this version: Henri Gouin. Thermodynamic form of the equation of motion for perfect fluids of grade n. Comptes
More informationAbout the implementation of the finite element method for computer aided education in electrical engineering
About the implementation of the finite element method for computer aided education in electrical engineering François Buret, Olivier Fabrègue, Daniel Muller, Alain Nicolas, Laurent Nicolas, Frédéric Thollon
More informationNumerical modeling of diffusion within composite media
Numerical modeling of diffusion within composite media Miljan Milosevic To cite this version: Miljan Milosevic. Numerical modeling of diffusion within composite media. 2nd ECCOMAS Young Investigators Conference
More informationCutwidth and degeneracy of graphs
Cutwidth and degeneracy of graphs Benoit Kloeckner To cite this version: Benoit Kloeckner. Cutwidth and degeneracy of graphs. IF_PREPUB. 2009. HAL Id: hal-00408210 https://hal.archives-ouvertes.fr/hal-00408210v1
More informationComputation and Experimental Measurements of the Magnetic Fields between Filamentary Circular Coils
Computation and Experimental Measurements of the Magnetic Fields between Filamentary Circular Coils Ao Anele, Y Hamam, L Chassagne, J Linares, Y Alayli, Karim Djouani To cite this version: Ao Anele, Y
More informationAnalyzing large-scale spike trains data with spatio-temporal constraints
Analyzing large-scale spike trains data with spatio-temporal constraints Hassan Nasser, Olivier Marre, Selim Kraria, Thierry Viéville, Bruno Cessac To cite this version: Hassan Nasser, Olivier Marre, Selim
More informationb-chromatic number of cacti
b-chromatic number of cacti Victor Campos, Claudia Linhares Sales, Frédéric Maffray, Ana Silva To cite this version: Victor Campos, Claudia Linhares Sales, Frédéric Maffray, Ana Silva. b-chromatic number
More informationNew estimates for the div-curl-grad operators and elliptic problems with L1-data in the half-space
New estimates for the div-curl-grad operators and elliptic problems with L1-data in the half-space Chérif Amrouche, Huy Hoang Nguyen To cite this version: Chérif Amrouche, Huy Hoang Nguyen. New estimates
More informationOn the Griesmer bound for nonlinear codes
On the Griesmer bound for nonlinear codes Emanuele Bellini, Alessio Meneghetti To cite this version: Emanuele Bellini, Alessio Meneghetti. On the Griesmer bound for nonlinear codes. Pascale Charpin, Nicolas
More informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
More informationDiurnal variation of tropospheric temperature at a tropical station
Diurnal variation of tropospheric temperature at a tropical station K. Revathy, S. R. Prabhakaran Nayar, B. V. Krishna Murthy To cite this version: K. Revathy, S. R. Prabhakaran Nayar, B. V. Krishna Murthy.
More informationMODal ENergy Analysis
MODal ENergy Analysis Nicolas Totaro, Jean-Louis Guyader To cite this version: Nicolas Totaro, Jean-Louis Guyader. MODal ENergy Analysis. RASD, Jul 2013, Pise, Italy. 2013. HAL Id: hal-00841467
More informationEaster bracelets for years
Easter bracelets for 5700000 years Denis Roegel To cite this version: Denis Roegel. Easter bracelets for 5700000 years. [Research Report] 2014. HAL Id: hal-01009457 https://hal.inria.fr/hal-01009457
More informationDG discretization of optimized Schwarz methods for Maxwell s equations
DG discretization of optimized Schwarz methods for Maxwell s equations Mohamed El Bouajaji, Victorita Dolean,3, Martin J. Gander 3, Stéphane Lanteri, and Ronan Perrussel 4 Introduction In the last decades,
More informationA Simple Model for Cavitation with Non-condensable Gases
A Simple Model for Cavitation with Non-condensable Gases Mathieu Bachmann, Siegfried Müller, Philippe Helluy, Hélène Mathis To cite this version: Mathieu Bachmann, Siegfried Müller, Philippe Helluy, Hélène
More informationCan we reduce health inequalities? An analysis of the English strategy ( )
Can we reduce health inequalities? An analysis of the English strategy (1997-2010) Johan P Mackenbach To cite this version: Johan P Mackenbach. Can we reduce health inequalities? An analysis of the English
More informationPositive mass theorem for the Paneitz-Branson operator
Positive mass theorem for the Paneitz-Branson operator Emmanuel Humbert, Simon Raulot To cite this version: Emmanuel Humbert, Simon Raulot. Positive mass theorem for the Paneitz-Branson operator. Calculus
More informationOn Solving Aircraft Conflict Avoidance Using Deterministic Global Optimization (sbb) Codes
On Solving Aircraft Conflict Avoidance Using Deterministic Global Optimization (sbb) Codes Sonia Cafieri, Frédéric Messine, Ahmed Touhami To cite this version: Sonia Cafieri, Frédéric Messine, Ahmed Touhami.
More informationAnalysis of Boyer and Moore s MJRTY algorithm
Analysis of Boyer and Moore s MJRTY algorithm Laurent Alonso, Edward M. Reingold To cite this version: Laurent Alonso, Edward M. Reingold. Analysis of Boyer and Moore s MJRTY algorithm. Information Processing
More informationFinite Element Method (FEM)
Finite Element Method (FEM) The finite element method (FEM) is the oldest numerical technique applied to engineering problems. FEM itself is not rigorous, but when combined with integral equation techniques
More informationHook lengths and shifted parts of partitions
Hook lengths and shifted parts of partitions Guo-Niu Han To cite this version: Guo-Niu Han Hook lengths and shifted parts of partitions The Ramanujan Journal, 009, 9 p HAL Id: hal-00395690
More informationOn the Earth s magnetic field and the Hall effect
On the Earth s magnetic field and the Hall effect J. E. Allen To cite this version: J. E. Allen. On the Earth s magnetic field and the Hall effect. Nonlinear rocesses in Geophysics, European Geosciences
More informationStickelberger s congruences for absolute norms of relative discriminants
Stickelberger s congruences for absolute norms of relative discriminants Georges Gras To cite this version: Georges Gras. Stickelberger s congruences for absolute norms of relative discriminants. Journal
More informationMultiple sensor fault detection in heat exchanger system
Multiple sensor fault detection in heat exchanger system Abdel Aïtouche, Didier Maquin, Frédéric Busson To cite this version: Abdel Aïtouche, Didier Maquin, Frédéric Busson. Multiple sensor fault detection
More informationImpedance Transmission Conditions for the Electric Potential across a Highly Conductive Casing
Impedance Transmission Conditions for the Electric Potential across a Highly Conductive Casing Hélène Barucq, Aralar Erdozain, David Pardo, Victor Péron To cite this version: Hélène Barucq, Aralar Erdozain,
More informationTropical Graph Signal Processing
Tropical Graph Signal Processing Vincent Gripon To cite this version: Vincent Gripon. Tropical Graph Signal Processing. 2017. HAL Id: hal-01527695 https://hal.archives-ouvertes.fr/hal-01527695v2
More informationVibro-acoustic simulation of a car window
Vibro-acoustic simulation of a car window Christophe Barras To cite this version: Christophe Barras. Vibro-acoustic simulation of a car window. Société Française d Acoustique. Acoustics 12, Apr 12, Nantes,
More informationOn path partitions of the divisor graph
On path partitions of the divisor graph Paul Melotti, Eric Saias To cite this version: Paul Melotti, Eric Saias On path partitions of the divisor graph 018 HAL Id: hal-0184801 https://halarchives-ouvertesfr/hal-0184801
More informationOn infinite permutations
On infinite permutations Dmitri G. Fon-Der-Flaass, Anna E. Frid To cite this version: Dmitri G. Fon-Der-Flaass, Anna E. Frid. On infinite permutations. Stefan Felsner. 2005 European Conference on Combinatorics,
More informationSOLAR RADIATION ESTIMATION AND PREDICTION USING MEASURED AND PREDICTED AEROSOL OPTICAL DEPTH
SOLAR RADIATION ESTIMATION AND PREDICTION USING MEASURED AND PREDICTED AEROSOL OPTICAL DEPTH Carlos M. Fernández-Peruchena, Martín Gastón, Maria V Guisado, Ana Bernardos, Íñigo Pagola, Lourdes Ramírez
More informationNorm Inequalities of Positive Semi-Definite Matrices
Norm Inequalities of Positive Semi-Definite Matrices Antoine Mhanna To cite this version: Antoine Mhanna Norm Inequalities of Positive Semi-Definite Matrices 15 HAL Id: hal-11844 https://halinriafr/hal-11844v1
More informationDGTD methods using modal basis functions and symplectic local time-stepping Application to wave propagation problems
DGTD methods using modal basis functions and symplectic local time-stepping Application to wave propagation problems Serge Piperno CERMICS ENPC 6-8 av. Blaise Pascal - Cité Descartes - Champs-sur-Marne
More informationTowards an active anechoic room
Towards an active anechoic room Dominique Habault, Philippe Herzog, Emmanuel Friot, Cédric Pinhède To cite this version: Dominique Habault, Philippe Herzog, Emmanuel Friot, Cédric Pinhède. Towards an active
More informationRENORMALISATION ON THE PENROSE LATTICE
RENORMALISATION ON THE PENROSE LATTICE C. Godreche, Henri Orland To cite this version: C. Godreche, Henri Orland. RENORMALISATION ON THE PENROSE LATTICE. Journal de Physique Colloques, 1986, 47 (C3), pp.c3-197-c3-203.
More informationA new approach of the concept of prime number
A new approach of the concept of prime number Jamel Ghannouchi To cite this version: Jamel Ghannouchi. A new approach of the concept of prime number. 4 pages. 24. HAL Id: hal-3943 https://hal.archives-ouvertes.fr/hal-3943
More informationA note on the acyclic 3-choosability of some planar graphs
A note on the acyclic 3-choosability of some planar graphs Hervé Hocquard, Mickael Montassier, André Raspaud To cite this version: Hervé Hocquard, Mickael Montassier, André Raspaud. A note on the acyclic
More informationGENERALIZED OPTICAL BISTABILITY AND CHAOS IN A LASER WITH A SATURABLE ABSORBER
GENERALIZED OPTICAL BISTABILITY AND CHAOS IN A LASER WITH A SATURABLE ABSORBER E. Arimondo, F. De Tomasi, B. Zambon, F. Papoff, D. Hennequin To cite this version: E. Arimondo, F. De Tomasi, B. Zambon,
More informationOn Poincare-Wirtinger inequalities in spaces of functions of bounded variation
On Poincare-Wirtinger inequalities in spaces of functions of bounded variation Maïtine Bergounioux To cite this version: Maïtine Bergounioux. On Poincare-Wirtinger inequalities in spaces of functions of
More informationDerivation and Analysis of Piecewise Constant Conservative Approximation for Anisotropic Diffusion Problems
Derivation Analysis of Piecewise Constant Conservative Approximation for Anisotropic Diffusion Problems A. Agouzal, Naïma Debit o cite this version: A. Agouzal, Naïma Debit. Derivation Analysis of Piecewise
More informationNumerical Modeling of Eddy Current Nondestructive Evaluation of Ferromagnetic Tubes via an Integral. Equation Approach
Numerical Modeling of Eddy Current Nondestructive Evaluation of Ferromagnetic Tubes via an Integral Equation Approach Anastassios Skarlatos, Grégoire Pichenot, Dominique Lesselier, Marc Lambert, Bernard
More informationA remark on a theorem of A. E. Ingham.
A remark on a theorem of A. E. Ingham. K G Bhat, K Ramachandra To cite this version: K G Bhat, K Ramachandra. A remark on a theorem of A. E. Ingham.. Hardy-Ramanujan Journal, Hardy-Ramanujan Society, 2006,
More informationThomas Lugand. To cite this version: HAL Id: tel
Contribution à la Modélisation et à l Optimisation de la Machine Asynchrone Double Alimentation pour des Applications Hydrauliques de Pompage Turbinage Thomas Lugand To cite this version: Thomas Lugand.
More informationSub-domain finite element method for efficiently considering strong skin and proximity effects
Sub-domain finite element method for efficiently considering strong skin and proximity effects Patrick Dular, Ruth Sabariego, Johan Gyselinck, Laurent Krähenbühl To cite this version: Patrick Dular, Ruth
More informationEntropies and fractal dimensions
Entropies and fractal dimensions Amelia Carolina Sparavigna To cite this version: Amelia Carolina Sparavigna. Entropies and fractal dimensions. Philica, Philica, 2016. HAL Id: hal-01377975
More informationParticle-in-cell simulations of high energy electron production by intense laser pulses in underdense plasmas
Particle-in-cell simulations of high energy electron production by intense laser pulses in underdense plasmas Susumu Kato, Eisuke Miura, Mitsumori Tanimoto, Masahiro Adachi, Kazuyoshi Koyama To cite this
More informationLinear Quadratic Zero-Sum Two-Person Differential Games
Linear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard To cite this version: Pierre Bernhard. Linear Quadratic Zero-Sum Two-Person Differential Games. Encyclopaedia of Systems and Control,
More informationThe Accelerated Euclidean Algorithm
The Accelerated Euclidean Algorithm Sidi Mohamed Sedjelmaci To cite this version: Sidi Mohamed Sedjelmaci The Accelerated Euclidean Algorithm Laureano Gonzales-Vega and Thomas Recio Eds 2004, University
More informationAnalysis of Defect in Extreme UV Lithography Mask Using a Modal Method Based on Nodal B-Spline Expansion
Analysis of Defect in Extreme UV Lithography Mas Using a Modal Method Based on odal B-Spline Expansion Kofi Edee, Patric Schiavone, Gérard Granet To cite this version: Kofi Edee, Patric Schiavone, Gérard
More informationNumerical Exploration of the Compacted Associated Stirling Numbers
Numerical Exploration of the Compacted Associated Stirling Numbers Khaled Ben Letaïef To cite this version: Khaled Ben Letaïef. Numerical Exploration of the Compacted Associated Stirling Numbers. 2017.
More informationVoltage Stability of Multiple Distributed Generators in Distribution Networks
oltage Stability of Multiple Distributed Generators in Distribution Networks Andi Wang, Chongxin Liu, Hervé Guéguen, Zhenquan Sun To cite this version: Andi Wang, Chongxin Liu, Hervé Guéguen, Zhenquan
More informationSome diophantine problems concerning equal sums of integers and their cubes
Some diophantine problems concerning equal sums of integers and their cubes Ajai Choudhry To cite this version: Ajai Choudhry. Some diophantine problems concerning equal sums of integers and their cubes.
More information