Modular hp-fem System HERMES and Its Application to the Maxwell s Equations 1
|
|
- Chrystal Simon
- 5 years ago
- Views:
Transcription
1 Modular hp-fem System HERMES and Its Application to the Maxwell s Equations 1 Tomáš Vejchodský a Pavel Šolín b Martin Zítka b a Mathematical Institute, Academy of Sciences, Žitná 25, Praha 1, Czech Republic b Department of Mathematical Sciences, University of Texas at El Paso, El Paso, Texas , USA Abstract In this paper we introduce the high-performance modular finite element system HERMES, a multi-physics hp-fem solver based on a novel approach where the finite element technology (mesh processing and adaptation, numerical quadrature, assembling and solution of the discrete problems, a-posteriori error estimation, etc.) is fully separated from the physics of the solved problems. The physics is represented via simple modules containing PDE-dependent parameters as well as hierarchic higher-order finite elements satisfying the conformity requirements imposed by the PDE. After describing briefly the modular structure of HERMES and some of its functionality, we focus on its application to the time-harmonic Maxwell s equations. We present numerical results which illustrate the capability of the hp-fem to reduce both the number of degrees of freedom and the CPU time dramatically compared to standard lowest-order FEM. Key words: hp-fem, time-harmonic Maxwell s equations, hierarchic higher-order edge elements 2000 MSC classification: 65N30 addresses: vejchod@math.cas.cz (Tomáš Vejchodský), solin@utep.edu (Pavel Šolín), zitka@math.utep.edu (Martin Zítka). 1 This work was supported by the Grant Agency of the Czech Republic, grants No. 201/04/P021 and No. 102/05/0629, and by the Academy of Sciences of the Czech Republic, Institutional Research Plan No. AV0Z Preprint submitted to Elsevier Science 5 September 2005
2 1 Introduction The hp-fem is a sophisticated version of the finite element method (FEM) which varies both the diameter and polynomial degree of elements in order to maximize the convergence rates. Its advantage over other numerical methods is an unconditional exponential convergence, even for problems with singular solutions. The method was first introduced in the 1980s [4,5,2,3], and its theoretical foundations are well established today. However, the understanding of practical aspects of the method (such as the design of optimal data structures and algorithms, error estimates, automatic hp-adaptivity, etc.), has not yet attained a sufficient level of maturity. This can be illustrated on the fact that the hp-fem has not yet become a standard in commercial engineering codes. Nowadays there are several groups both in the U.S. [7,8,11,13,14,17 19,24,23] and in Europe [1,10,15,20] dealing with the hp-fem. In all cases, the investigation of the practical aspects is a substantial part of their research effort. 2 The hp-fem system HERMES The multi-physics modular finite element system HERMES has been developed by the authors over the last several years. The system is based on a novel modular approach where the finite element technology is completely separated from the physics of the problems represented by concrete PDEs. The modular structure of the system is illustrated in Fig. 1. FEM/hp FEM Module Continous Problem Geometry PDE Bdy. Cond. Mesh Generator Weak Formulation Edge elements Continuous elements Taylor Hood elements FEM/hp FEM Kernel Other types of elements Output of Results Visualization, etc. Discrete Problem smatrix Interface Solution stiffness matrix load vector Trilinos PETSc UMFPACK Other Algebraic Module solution coefficient vector Fig. 1. The modular structure of HERMES. 2
3 The heart of the solver is the FEM/hp-FEM Module which contains all PDEindependent algorithms, such as processing and adaptation of finite element meshes, numerical quadrature, assembling procedures, solution of systems of linear and nonlinear algebraic equations, a-posteriori error estimation, etc. The physics of the problem is represented via smaller PDE modules which are independent of the FEM/hp-FEM Module. The PDE modules contain not only the definition of concrete PDEs along with related physical parameters, but also scalar or vector-valued finite elements which are appropriate for their discretization. The PDE modules are strictly separated from the FEM/hp- FEM Module at each level of the hierarchic data structures. Last we would like to mention the module smatrix which serves as a universal interface between the FEM/hp-FEM Module and the Algebraic Module. The latter contains solver packages for sparse systems of linear and nonlinear algebraic equations, such as Trilinos [12], PETSc [6], and UMFPACK [9]. When these modules are not installed on the system where HERMES is run, smatrix can use a set of standard ILU-preconditioned matrix solvers. It is not technically difficult to add new solvers to the Algebraic Module. 3 Application to time-harmonic Maxwell s equations In addition to a module for systems of nonlinear elliptic equations and a module for the incompressible Navier Stokes equations, HERMES contains a time-harmonic Maxwell s equations module. Let us describe this module in more detail. Consider a bounded polygonal domain Ω R 2 whose boundary Ω is split into two relatively open disjoint parts Γ P and Γ I representing standard perfect conducting and impedance boundary conditions, respectively. The classical (strong) formulation of the Maxwell s equations reads: curl ( µ 1 r curl E ) κ 2 ɛ r E = F in Ω, E τ = 0 on Γ P, (1) µ 1 r curl E iκλe τ = g τ on Γ I. Here curl c = ( c/ x 2, c/ x 1 ) and curl b = b 2 / x 1 b 1 / x 2 are the standard vector and scalar curl operators, τ = ( ν 2, ν 1 ) is the positivelyoriented unit tangent vector to Ω, ν = (ν 1, ν 2 ) is the unit outer normal vector to Ω, and the symbol i stands for the imaginary unit. Some related notations are presented in Table 1. 3
4 E = E(x) C 2 phasor of the electric field strength (unknown) µ r = µ r (x) R relative permeability ɛ r = ɛ r (x) C 2 2 relative permittivity κ = const. R wave number λ = λ(x) > 0 impedance F = F(x) C 2 right-hand side of the PDE g = g(x) C 2 right-hand side of the impedance boundary condition Table 1 Quantities used in the time-harmonic Maxwell s equations. Problem (1) is formulated in the weak sense as follows: Find E V = {E H(curl, Ω) : E τ = 0 on Γ P } such that a(e, Φ) = F(Φ) Φ V, where the sesquilinear form a and the antilinear functional F are given by a(e, Φ) = F(Φ) = Ω Ω µ 1 r curl E curl Φ dx κ 2 F Φ dx + (g τ)(φ τ) ds. Γ I Ω (ɛ r E) Φ dx iκ Γ I λ(e τ)(φ τ) ds, We assume a triangulation T h,p of the domain Ω and assign a polynomial degree p j 0 to every element K j T h,p. This hp-mesh defines the following piecewise polynomial subspace of V : V h,p = { E h,p V : E h,p Kj P p j (K j ) and the tangent componet E h,p τ k is continuous on each edge e k E h,p }, where E h,p is the set of all edges in T h,p, τ k is the tangent vector to the edge e k and P p j (K j ) denotes the space of vector valued polynomials defined on a triangle K j with both components being of degree p j. With the finite element space V h,p problem: Find E h,p V h,p such that in hand, we can formulate the discrete a(e h,p, Φ h,p ) = F(Φ h,p ) Φ h,p V h,p. (2) The solution of this system of linear algebraic equations yields the approximate solution E h,p. For more details on the hp-fem including edge elements we refer to [24,22]. 4
5 4 Numerical experiments Let us use HERMES to illustrate the superior performance of the hp-fem over the standard low-order FEM. Consider the time-harmonic Maxwell s equations in the L-shape domain Ω with the vertices B 1 = (0, 0), B 2 = (0, 1), B 3 = (1, 1), B 4 = (1, 1), B 5 = ( 1, 1), and B 6 = ( 1, 0). Perfect conducting boundary conditions are prescribed on the edges meeting at the origin (B 1 B 2 and B 6 B 1 ) and the rest of Ω is equipped with impedance boundary conditions. For simplicity, the equation coefficients are chosen constant: µ r = 1, ɛ r = I (2 2 identity matrix), κ = 1, λ = 1. The exact solution is given by E = u, where u = r 2 3 sin ((2θ + π)/3). Hence, E 1 = 2 3 r 1 3 cos ( π 6 + θ 3 ), E 2 = 2 3 r 1 3 sin ( π 6 + θ 3 The solution is shown in the left part of Fig. 2. Notice the singularity at the reentrant corner (truncated for visualization purposes). ). Fig. 2. Example 1: Left: Modulus of the electric field phasor E. Right: Geometry of the mesh. The right hand sides F and g τ are chosen to be compatible with the exact solution: F = E, and g τ = ie τ on Ω. The problem was solved by the lowest-order FEM (p = 0 in all elements), piecewise-linear FEM (p = 1 in all elements), and by the hp-fem. We used the mesh shown in Fig. 2. In the first two cases we had to refine this meshes uniformly in order to obtain an accuracy comparable to the hp-fem: Every element in the mesh was subdivided into 10,000 triangles in the case p = 0 and into 484 triangles in the case p = 1. The distribution of the polynomial degree in the hp-mesh is shown in Fig. 3. 5
6 Fig. 3. Example 1: Meshes used for the hp-fem with detailed views of the re-entrant corners. The polynomial degree varies from 2 (blue) to 4 (red) in the left part and from 2 (blue) to 7 (red) in the right part. Table 2 compares the lowest-order case (p = 0) with the hp-fem on the relative error level of approx %. One can see that the hp-fem performed approx times faster compared to the lowest-order method. The memory requirements of the hp-fem were about times less. DOFs CPU time rel. error p = min 26 s % hp s % Ratio Table 2 Example 1: performance of the lowest-order and hp-fem. The other experiment compares the hp-fem (see the right part of Fig. 3) to the first-order method (p = 1). The corresponding results are presented in Table 3. In this case, the hp-fem was about 100 times faster and used 50 times less computer memory. DOFs CPU time rel. error p = min 18 s % hp s % Ratio Table 3 Example 1: performance of the first-order and hp-fem. 5 Summary and outlook We presented a multi-physics modular hp-fem system HERMES based on a novel approach, where the finite element technology is fully separated from the physics of the solved problems. In this way, new techniques implemented in the FEM/hp-FEM kernel work automatically for various PDEs, such as systems of linear and nonlinear second-order elliptic equations, time-harmonic 6
7 Maxwell s equations, convection-diffusion equations, Stokes problem, incompressible Navier Stokes equations, etc. The modular structure cuts down the development cost dramatically, and it greatly facilitates team-work. Due to its modular structure, HERMES allows for performing various numerical experiments (such as testing various sets of shape functions, various a-posteriori error estimators, various hp-adaptive strategies) very easily. The development of HERMES is by no means finished. Currently the implementation of hierarchic higher-order Taylor Hood elements is in progress, as well as an implementation of a novel approach to constrained approximation with multiple constraint levels, automatic hp-adaptivity, parallelization of the code, and other tasks. References [1] M. Ainsworth, B. Senior: Aspects of an hp-adaptive finite element method: Adaptive strategy, conforming approximation and efficient solvers, Technical Report 1997/2, Department of Mathematics and Computer Science, University of Leicester, England, [2] I. Babuška, W. Gui: The h, p and hp-versions of the finite element method in 1 dimension - Part I. The error analysis of the p-version, Numer. Math. 49 (1986), [3] I. Babuška, W. Gui: The h, p and hp-versions of the finite element method in 1 dimension - Part II. The error analysis of the h and hp-versions, Numer. Math. 49 (1986), [4] I. Babuška, M. Suri: The hp-version of the finite element method with quasiuniform meshes, Model. Math. Anal. Numer. 21 (1987), [5] I. Babuška, B. Szabo, I.N. Katz: The p-version of the finite element method, SIAM J. Numer. Anal. 18 (1981), [6] S. Balay, K. Buschelman, V. Eijkhout, W.D. Gropp, D. Kaushik, M.G. Knepley, L.C. McInnes, B.F. Smith, H. Zhang: PETSc Users Manual, Tech. Report ANL- 95/11, Argonne National Laboratory, [7] A.C. Bauer, A.K. Patra: Performance of parallel preconditioners for adaptive hp-fem discretizations of incompressible flows, Commun. Numer. Meth. Engrg. 18 (2002), [8] A.C. Bauer, A.K. Patra: Robust and efficient domain decomposition preconditioners for adaptive hp finite element approximation for linear elasticity with and without discontinuous coefficients, Int. J. Numer. Meth. Engrg. 59 (2004),
8 [9] T.A. Davis: Algorithm 832: UMFPACK - an unsymmetric-pattern multifrontal method with a column pre-ordering strategy, ACM Trans. Math. Software 30 (2004), [10] I. Doležel, P. Šolín, M. Zítka: On the hp-fem for Singular Electrostatics Problems, In: Proceedings of ISEF, September 2005, Baiona, Spain, to appear. [11] S. Iqbal, G.F. Carey: Performance analysis of dynamic load balancing algorithms with variable number of processors, J. Parallel Distrib. Comput., to appear. [12] M.A. Heroux, R.A. Bartlett, V.E. Howle, R.J. Hoekstra, J.J. Hu, T.G. Kolda, R.B. Lehoucq, K.R. Long, R.P. Pawlowski, E.T. Phipps, A.G. Salinger, H.K. Thornquist, R.S. Tuminaro, J.M. Willenbring, A. Williams, K.S. Stanley, An Overview of the Trilinos Project, ACM Trans. Math. Software, accepted, December [13] G.E. Karniadakis, S.J. Sherwin: Spectral/hp Element Methods for CFD, Oxford University Press, Oxford, [14] A. Laszloffy, J. Long, A.K. Patra: Simple data management, scheduling and solution strategies for managing the irregularities in parallel adaptive hp finite element simulations, Parallel Computing 26 (2000), [15] J.M. Melenk: hp-finite Element Methods for Singular Perturbations, Lecture Notes in Math. 1796, Springer-Verlag, Berlin, [16] P. Monk: Finite element methods for Maxwell s equations. Oxford University Press, New York, [17] M. Paszynski, J. Kurtz, L. Demkowicz: Parallel, fully automatic hp-adaptive 2D finite element package, TICAM Report 04-07, The University of Texas at Austin, [18] A.K. Patra et al.: Parallel adaptive numerical simulation of dry avalanches over natural terrain, J. Volcan. Geotherm. Research 139 (2005), [19] W. Rachowicz, D. Pardo, L. Demkowicz: Fully automatic hp-adaptivity in three dimensions, ICES Report 04-22, The University of Texas at Austin, [20] Ch. Schwab: p- and hp-finite Element Methods, Clarendon Press, Oxford, [21] J.R. Shewchuk: Delaunay Refinement Algorithms for Triangular Mesh Generation. Comput. Geom. 22 (2002), [22] P. Šolín: Partial differential equations and the finite element method, John Wiley & Sons, [23] P. Šolín, L. Demkowicz: Goal-Oriented hp-adaptivity for Elliptic Problems, Comput. Methods Appl. Mech. Engrg. 193 (2004), [24] P. Šolín, K. Segeth, I. Doležel: Higher-order finite element methods. Chapman & Hall/CRC, Boca Raton, FL,
Discrete Maximum Principle for a 1D Problem with Piecewise-Constant Coefficients Solved by hp-fem
The University of Texas at El Paso Department of Mathematical Sciences Research Reports Series El Paso, Texas Research Report No. 2006-10 Discrete Maximum Principle for a 1D Problem with Piecewise-Constant
More informationInterval Finite Element Methods: New Directions
Interval Finite Element Methods: New Directions Rafi Muhanna 1, Vladik Kreinovich 2, Pavel Šolín2, Jack Chessa 2, Roberto Araiza 2, and Gang Xiang 2 1 Center for Reliable Engineering Computing (REC), Department
More informationPartial Differential Equations and the Finite Element Method
Partial Differential Equations and the Finite Element Method Pavel Solin The University of Texas at El Paso Academy of Sciences ofthe Czech Republic iwiley- INTERSCIENCE A JOHN WILEY & SONS, INC, PUBLICATION
More informationNUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL
Proceedings of ALGORITMY 212 pp. 29 218 NUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL MARTIN HADRAVA, MILOSLAV FEISTAUER, AND PETR SVÁČEK Abstract. The present paper
More informationA posteriori error estimates applied to flow in a channel with corners
Mathematics and Computers in Simulation 61 (2003) 375 383 A posteriori error estimates applied to flow in a channel with corners Pavel Burda a,, Jaroslav Novotný b, Bedřich Sousedík a a Department of Mathematics,
More informationGeneralized Finite Element Methods for Three Dimensional Structural Mechanics Problems. C. A. Duarte. I. Babuška and J. T. Oden
Generalized Finite Element Methods for Three Dimensional Structural Mechanics Problems C. A. Duarte COMCO, Inc., 7800 Shoal Creek Blvd. Suite 290E Austin, Texas, 78757, USA I. Babuška and J. T. Oden TICAM,
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationA Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationNumerical Solutions of Laplacian Problems over L-Shaped Domains and Calculations of the Generalized Stress Intensity Factors
WCCM V Fifth World Congress on Computational Mechanics July 7-2, 2002, Vienna, Austria Eds.: H.A. Mang, F.G. Rammerstorfer, J. Eberhardsteiner Numerical Solutions of Laplacian Problems over L-Shaped Domains
More informationPreconditioned Eigensolver LOBPCG in hypre and PETSc
Preconditioned Eigensolver LOBPCG in hypre and PETSc Ilya Lashuk, Merico Argentati, Evgueni Ovtchinnikov, and Andrew Knyazev Department of Mathematics, University of Colorado at Denver, P.O. Box 173364,
More informationElectrostatic microactuators for precise positioning and comparison of their parameters
Computer Applications in Electrical Engineering Electrostatic microactuators for precise positioning and comparison of their parameters Václav Štarman, Jan Kacerovský, Jindřich Jansa, Pavel Karban, Ivo
More informationPreprint Alexander Heinlein, Axel Klawonn, and Oliver Rheinbach Parallel Two-Level Overlapping Schwarz Methods in Fluid-Structure Interaction
Fakultät für Mathematik und Informatik Preprint 2015-15 Alexander Heinlein, Axel Klawonn, and Oliver Rheinbach Parallel Two-Level Overlapping Schwarz Methods in Fluid-Structure Interaction ISSN 1433-9307
More informationA Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation
A Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation Tao Zhao 1, Feng-Nan Hwang 2 and Xiao-Chuan Cai 3 Abstract In this paper, we develop an overlapping domain decomposition
More informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More informationEdwin van der Weide and Magnus Svärd. I. Background information for the SBP-SAT scheme
Edwin van der Weide and Magnus Svärd I. Background information for the SBP-SAT scheme As is well-known, stability of a numerical scheme is a key property for a robust and accurate numerical solution. Proving
More informationParallel scalability of a FETI DP mortar method for problems with discontinuous coefficients
Parallel scalability of a FETI DP mortar method for problems with discontinuous coefficients Nina Dokeva and Wlodek Proskurowski University of Southern California, Department of Mathematics Los Angeles,
More informationFrom the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes
www.oeaw.ac.at From the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes D. Copeland, U. Langer, D. Pusch RICAM-Report 2008-10 www.ricam.oeaw.ac.at From the Boundary Element
More informationAdaptive Time Space Discretization for Combustion Problems
Konrad-Zuse-Zentrum fu r Informationstechnik Berlin Takustraße 7 D-14195 Berlin-Dahlem Germany JENS LANG 1,BODO ERDMANN,RAINER ROITZSCH Adaptive Time Space Discretization for Combustion Problems 1 Talk
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES. A virtual overlapping Schwarz method for scalar elliptic problems in two dimensions
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES A virtual overlapping Schwarz method for scalar elliptic problems in two dimensions Juan Gabriel Calvo Preprint No. 25-2017 PRAHA 2017 A VIRTUAL
More informationKey words. Parallel iterative solvers, saddle-point linear systems, preconditioners, timeharmonic
PARALLEL NUMERICAL SOLUTION OF THE TIME-HARMONIC MAXWELL EQUATIONS IN MIXED FORM DAN LI, CHEN GREIF, AND DOMINIK SCHÖTZAU Numer. Linear Algebra Appl., Vol. 19, pp. 525 539, 2012 Abstract. We develop a
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationFrom the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes
From the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes Dylan Copeland 1, Ulrich Langer 2, and David Pusch 3 1 Institute of Computational Mathematics,
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationParallel Scalability of a FETI DP Mortar Method for Problems with Discontinuous Coefficients
Parallel Scalability of a FETI DP Mortar Method for Problems with Discontinuous Coefficients Nina Dokeva and Wlodek Proskurowski Department of Mathematics, University of Southern California, Los Angeles,
More informationOn angle conditions in the finite element method. Institute of Mathematics, Academy of Sciences Prague, Czech Republic
On angle conditions in the finite element method Michal Křížek Institute of Mathematics, Academy of Sciences Prague, Czech Republic Joint work with Jan Brandts (University of Amsterdam), Antti Hannukainen
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationConvergence Behavior of a Two-Level Optimized Schwarz Preconditioner
Convergence Behavior of a Two-Level Optimized Schwarz Preconditioner Olivier Dubois 1 and Martin J. Gander 2 1 IMA, University of Minnesota, 207 Church St. SE, Minneapolis, MN 55455 dubois@ima.umn.edu
More informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
More informationConvergence of an adaptive hp finite element strategy in higher space-dimensions
Convergence of an adaptive hp finite element strategy in higher space-dimensions M. Bürg W. Dörfler Preprint Nr. 1/6 INSTITUT FÜR WISSENSCHAFTLICHES RECHNEN UND MATHEMATISCHE MODELLBILDUNG Anschriften
More informationThe Plane Stress Problem
The Plane Stress Problem Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) February 2, 2010 Martin Kronbichler (TDB) The Plane Stress Problem February 2, 2010 1 / 24 Outline
More informationDiscretization of PDEs and Tools for the Parallel Solution of the Resulting Systems
Discretization of PDEs and Tools for the Parallel Solution of the Resulting Systems Stan Tomov Innovative Computing Laboratory Computer Science Department The University of Tennessee Wednesday April 4,
More information33 RASHO: A Restricted Additive Schwarz Preconditioner with Harmonic Overlap
Thirteenth International Conference on Domain Decomposition ethods Editors: N. Debit,.Garbey, R. Hoppe, J. Périaux, D. Keyes, Y. Kuznetsov c 001 DD.org 33 RASHO: A Restricted Additive Schwarz Preconditioner
More informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More informationSpace-time XFEM for two-phase mass transport
Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More informationInstitut de Recherche MAthématique de Rennes
LMS Durham Symposium: Computational methods for wave propagation in direct scattering. - July, Durham, UK The hp version of the Weighted Regularization Method for Maxwell Equations Martin COSTABEL & Monique
More informationComparative Analysis of Mesh Generators and MIC(0) Preconditioning of FEM Elasticity Systems
Comparative Analysis of Mesh Generators and MIC(0) Preconditioning of FEM Elasticity Systems Nikola Kosturski and Svetozar Margenov Institute for Parallel Processing, Bulgarian Academy of Sciences Abstract.
More informationNodal O(h 4 )-superconvergence of piecewise trilinear FE approximations
Preprint, Institute of Mathematics, AS CR, Prague. 2007-12-12 INSTITTE of MATHEMATICS Academy of Sciences Czech Republic Nodal O(h 4 )-superconvergence of piecewise trilinear FE approximations Antti Hannukainen
More informationFinal Ph.D. Progress Report. Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo
Final Ph.D. Progress Report Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo Dissertation Committee: I. Babuska, L. Demkowicz, C. Torres-Verdin, R. Van
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More informationSpace-time Finite Element Methods for Parabolic Evolution Problems
Space-time Finite Element Methods for Parabolic Evolution Problems with Variable Coefficients Ulrich Langer, Martin Neumüller, Andreas Schafelner Johannes Kepler University, Linz Doctoral Program Computational
More informationA High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations
A High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations Khosro Shahbazi 1, Paul F. Fischer 2 and C. Ross Ethier 1 1 University of Toronto and 2 Argonne National
More informationAdaptive methods for control problems with finite-dimensional control space
Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy
More informationNUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING TECHNIQUES SATISFYING THE DISCRETE MAXIMUM PRINCIPLE
Proceedings of the Czech Japanese Seminar in Applied Mathematics 2005 Kuju Training Center, Oita, Japan, September 15-18, 2005 pp. 69 76 NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationLecture 8: Boundary Integral Equations
CBMS Conference on Fast Direct Solvers Dartmouth College June 23 June 27, 2014 Lecture 8: Boundary Integral Equations Gunnar Martinsson The University of Colorado at Boulder Research support by: Consider
More informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
More informationOn discontinuity capturing methods for convection diffusion equations
On discontinuity capturing methods for convection diffusion equations Volker John 1 and Petr Knobloch 2 1 Universität des Saarlandes, Fachbereich 6.1 Mathematik, Postfach 15 11 50, 66041 Saarbrücken, Germany,
More informationSchwarz type solvers forhp-fem discretizations of mixed problems
Wegelerstraße 6 53115 Bonn Germany phone +49 228 73-3427 fax +49 228 73-7527 www.ins.uni-bonn.de S. Beuchler, M. Purrucker Schwarz type solvers forhp-fem discretizations of mixed problems INS Preprint
More informationPAijpam.eu NEW H 1 (Ω) CONFORMING FINITE ELEMENTS ON HEXAHEDRA
International Journal of Pure and Applied Mathematics Volume 109 No. 3 2016, 609-617 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu doi: 10.12732/ijpam.v109i3.10
More informationarxiv: v2 [math.na] 23 Apr 2016
Improved ZZ A Posteriori Error Estimators for Diffusion Problems: Conforming Linear Elements arxiv:508.009v2 [math.na] 23 Apr 206 Zhiqiang Cai Cuiyu He Shun Zhang May 2, 208 Abstract. In [8], we introduced
More informationEfficiency-based h- and hp-refinement strategies for finite element methods
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 7; : 5 [Version: /9/8 v.] Efficiency-based h- and hp-refinement strategies for finite element methods H. De Sterck, Thomas A. Manteuffel,
More informationA Posteriori Error Estimation Techniques for Finite Element Methods. Zhiqiang Cai Purdue University
A Posteriori Error Estimation Techniques for Finite Element Methods Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 16, 2017 Books Ainsworth & Oden, A posteriori
More informationDiscontinuous Galerkin Method for interface problem of coupling different order elliptic equations
Discontinuous Galerkin Method for interface problem of coupling different order elliptic equations Igor Mozolevski, Endre Süli Federal University of Santa Catarina, Brazil Oxford University Computing Laboratory,
More information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
More informationINTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR GRIDS. 1. Introduction
Trends in Mathematics Information Center for Mathematical Sciences Volume 9 Number 2 December 2006 Pages 0 INTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR
More informationComplementarity based a posteriori error estimates and their properties
Complementarity based a posteriori error estimates and their properties Tomáš Vejchodský November 3, 2009 Institute of Mathematics, Czech Academy of Sciences Žitná 25, CZ-115 67 Praha 1, Czech Republic
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY 11, USA Dan Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More informationFINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION
Proceedings of ALGORITMY pp. 9 3 FINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION JAN STEBEL Abstract. The paper deals with the numerical simulations of steady flows
More informationA Parallel Scalable PETSc-Based Jacobi-Davidson Polynomial Eigensolver with Application in Quantum Dot Simulation
A Parallel Scalable PETSc-Based Jacobi-Davidson Polynomial Eigensolver with Application in Quantum Dot Simulation Zih-Hao Wei 1, Feng-Nan Hwang 1, Tsung-Ming Huang 2, and Weichung Wang 3 1 Department of
More informationFinite difference computation of the permeability of textile reinforcements with a fast Stokes solver and new validation examples
Finite difference computation of the permeability of textile reinforcements with a fast Stokes solver and new validation examples B. Verleye, R. Croce, M. Griebel, M. Klitz, S.V. Lomov, I. Verpoest and
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationOn the equivalence of regularity criteria for triangular and tetrahedral finite element partitions
Computers and Mathematics with Applications 55 (2008) 2227 2233 www.elsevier.com/locate/camwa On the equivalence of regularity criteria for triangular and tetrahedral finite element partitions Jan Brandts
More informationConstruction of a New Domain Decomposition Method for the Stokes Equations
Construction of a New Domain Decomposition Method for the Stokes Equations Frédéric Nataf 1 and Gerd Rapin 2 1 CMAP, CNRS; UMR7641, Ecole Polytechnique, 91128 Palaiseau Cedex, France 2 Math. Dep., NAM,
More informationHIGHER-ORDER LINEARLY IMPLICIT ONE-STEP METHODS FOR THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
STUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume LIII, Number 1, March 2008 HIGHER-ORDER LINEARLY IMPLICIT ONE-STEP METHODS FOR THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IOAN TELEAGA AND JENS
More informationLinear Solvers. Andrew Hazel
Linear Solvers Andrew Hazel Introduction Thus far we have talked about the formulation and discretisation of physical problems...... and stopped when we got to a discrete linear system of equations. Introduction
More informationChapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction
Chapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction The a posteriori error estimation of finite element approximations of elliptic boundary value problems has reached
More information18. Balancing Neumann-Neumann for (In)Compressible Linear Elasticity and (Generalized) Stokes Parallel Implementation
Fourteenth nternational Conference on Domain Decomposition Methods Editors: smael Herrera, David E Keyes, Olof B Widlund, Robert Yates c 23 DDMorg 18 Balancing Neumann-Neumann for (n)compressible Linear
More informationTwo-scale Dirichlet-Neumann preconditioners for boundary refinements
Two-scale Dirichlet-Neumann preconditioners for boundary refinements Patrice Hauret 1 and Patrick Le Tallec 2 1 Graduate Aeronautical Laboratories, MS 25-45, California Institute of Technology Pasadena,
More informationTopology optimisation of passive coolers for light-emitting diode lamps
Downloaded from orbit.dtu.dk on: Aug 15, 2018 Topology optimisation of passive coolers for light-emitting diode lamps Alexandersen, Joe Published in: Proceedings of WCSMO-11 Publication date: 2015 Document
More informationInteraction of Incompressible Fluid and Moving Bodies
WDS'06 Proceedings of Contributed Papers, Part I, 53 58, 2006. ISBN 80-86732-84-3 MATFYZPRESS Interaction of Incompressible Fluid and Moving Bodies M. Růžička Charles University, Faculty of Mathematics
More informationExact a posteriori error analysis of the least squares nite element method 1
Applied Mathematics and Computation 116 (2000) 297±305 www.elsevier.com/locate/amc Exact a posteriori error analysis of the least squares nite element method 1 Jinn-Liang Liu * Department of Applied Mathematics,
More informationIn Proc. of the V European Conf. on Computational Fluid Dynamics (ECFD), Preprint
V European Conference on Computational Fluid Dynamics ECCOMAS CFD 2010 J. C. F. Pereira and A. Sequeira (Eds) Lisbon, Portugal, 14 17 June 2010 THE HIGH ORDER FINITE ELEMENT METHOD FOR STEADY CONVECTION-DIFFUSION-REACTION
More informationA NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION
A NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION JOHNNY GUZMÁN, ABNER J. SALGADO, AND FRANCISCO-JAVIER SAYAS Abstract. The analysis of finite-element-like Galerkin discretization techniques for the
More informationHETEROGENEOUS MULTISCALE METHOD IN EDDY CURRENTS MODELING
Proceedings of ALGORITMY 2009 pp. 219 225 HETEROGENEOUS MULTISCALE METHOD IN EDDY CURRENTS MODELING JÁN BUŠA, JR. AND VALDEMAR MELICHER Abstract. The induction of eddy currents in a conductive piece is
More informationParallel numerical solution of the time-harmonic Maxwell equations in mixed form
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. (2011) Published online in Wiley Online Library (wileyonlinelibrary.com)..782 Parallel numerical solution of the time-harmonic Maxwell
More informationDomain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions
Domain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions Bernhard Hientzsch Courant Institute of Mathematical Sciences, New York University, 51 Mercer Street, New
More informationFinite Element Superconvergence Approximation for One-Dimensional Singularly Perturbed Problems
Finite Element Superconvergence Approximation for One-Dimensional Singularly Perturbed Problems Zhimin Zhang Department of Mathematics Wayne State University Detroit, MI 48202 Received 19 July 2000; accepted
More informationOptimal multilevel preconditioning of strongly anisotropic problems.part II: non-conforming FEM. p. 1/36
Optimal multilevel preconditioning of strongly anisotropic problems. Part II: non-conforming FEM. Svetozar Margenov margenov@parallel.bas.bg Institute for Parallel Processing, Bulgarian Academy of Sciences,
More informationA Two-grid Method for Coupled Free Flow with Porous Media Flow
A Two-grid Method for Coupled Free Flow with Porous Media Flow Prince Chidyagwai a and Béatrice Rivière a, a Department of Computational and Applied Mathematics, Rice University, 600 Main Street, Houston,
More informationApplied Mathematics and Computation 169 (2005)
Applied Mathematics and Computation 169 (2005) 485 499 www.elsevier.com/locate/amc Solving Laplacian problems with boundary singularities: a comparison of a singular function boundary integral method with
More informationWeak Galerkin Finite Element Methods and Applications
Weak Galerkin Finite Element Methods and Applications Lin Mu mul1@ornl.gov Computational and Applied Mathematics Computationa Science and Mathematics Division Oak Ridge National Laboratory Georgia Institute
More informationNewton additive and multiplicative Schwarz iterative methods
IMA Journal of Numerical Analysis (2008) 28, 143 161 doi:10.1093/imanum/drm015 Advance Access publication on July 16, 2007 Newton additive and multiplicative Schwarz iterative methods JOSEP ARNAL, VIOLETA
More informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationANALYSIS OF A FINITE ELEMENT PML APPROXIMATION FOR THE THREE DIMENSIONAL TIME-HARMONIC MAXWELL PROBLEM
MATHEMATICS OF COMPUTATION Volume 77, Number 261, January 2008, Pages 1 10 S 0025-5718(07)02037-6 Article electronically published on September 18, 2007 ANALYSIS OF A FINITE ELEMENT PML APPROXIMATION FOR
More informationA Review of Preconditioning Techniques for Steady Incompressible Flow
Zeist 2009 p. 1/43 A Review of Preconditioning Techniques for Steady Incompressible Flow David Silvester School of Mathematics University of Manchester Zeist 2009 p. 2/43 PDEs Review : 1984 2005 Update
More informationA Posteriori Estimates for Cost Functionals of Optimal Control Problems
A Posteriori Estimates for Cost Functionals of Optimal Control Problems Alexandra Gaevskaya, Ronald H.W. Hoppe,2 and Sergey Repin 3 Institute of Mathematics, Universität Augsburg, D-8659 Augsburg, Germany
More informationEnergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations of the Navier-Stokes equations
INTRNATIONAL JOURNAL FOR NUMRICAL MTHODS IN FLUIDS Int. J. Numer. Meth. Fluids 19007; 1:1 [Version: 00/09/18 v1.01] nergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations
More informationNUMERICAL MODELING OF TRANSIENT ACOUSTIC FIELD USING FINITE ELEMENT METHOD
POZNAN UNIVE RSITY OF TE CHNOLOGY ACADE MIC JOURNALS No 73 Electrical Engineering 213 Lukáš KOUDELA* Jindřich JANSA* Pavel KARBAN* NUMERICAL MODELING OF TRANSIENT ACOUSTIC FIELD USING FINITE ELEMENT METHOD
More informationA Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers
Applied and Computational Mathematics 2017; 6(4): 202-207 http://www.sciencepublishinggroup.com/j/acm doi: 10.11648/j.acm.20170604.18 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Robust Preconditioned
More informationRESIDUAL BASED ERROR ESTIMATES FOR THE SPACE-TIME DISCONTINUOUS GALERKIN METHOD APPLIED TO NONLINEAR HYPERBOLIC EQUATIONS
Proceedings of ALGORITMY 2016 pp. 113 124 RESIDUAL BASED ERROR ESTIMATES FOR THE SPACE-TIME DISCONTINUOUS GALERKIN METHOD APPLIED TO NONLINEAR HYPERBOLIC EQUATIONS VÍT DOLEJŠÍ AND FILIP ROSKOVEC Abstract.
More informationAn a posteriori error estimator for the weak Galerkin least-squares finite-element method
An a posteriori error estimator for the weak Galerkin least-squares finite-element method James H. Adler a, Xiaozhe Hu a, Lin Mu b, Xiu Ye c a Department of Mathematics, Tufts University, Medford, MA 02155
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationLecture 8: Fast Linear Solvers (Part 7)
Lecture 8: Fast Linear Solvers (Part 7) 1 Modified Gram-Schmidt Process with Reorthogonalization Test Reorthogonalization If Av k 2 + δ v k+1 2 = Av k 2 to working precision. δ = 10 3 2 Householder Arnoldi
More informationDomain-decomposed Fully Coupled Implicit Methods for a Magnetohydrodynamics Problem
Domain-decomposed Fully Coupled Implicit Methods for a Magnetohydrodynamics Problem Serguei Ovtchinniov 1, Florin Dobrian 2, Xiao-Chuan Cai 3, David Keyes 4 1 University of Colorado at Boulder, Department
More information