103 Assyr Abdulle · Simone Deparis Daniel Kressner · Fabio Nobile Marco Picasso Editors Numerical Mathematics and Advanced Applications ENUMATH 2013 Editorial Board T. J.Barth M.Griebel D.E.Keyes R.M.Nieminen D.Roose T.Schlick Lecture Notes in Computational Science 103 and Engineering Editors: TimothyJ.Barth MichaelGriebel DavidE.Keyes RistoM.Nieminen DirkRoose TamarSchlick Moreinformationaboutthisseriesat http://www.springer.com/series/3527 Assyr Abdulle • Simone Deparis (cid:129) Daniel Kressner (cid:129) Fabio Nobile (cid:129) Marco Picasso Editors Numerical Mathematics and Advanced Applications ENUMATH 2013 Proceedings of ENUMATH 2013, the 10th European Conference on Numerical Mathematics and Advanced Applications, Lausanne, August 2013 123 Editors AssyrAbdulle SimoneDeparis DanielKressner FabioNobile MarcoPicasso EcolePolytechniqueFédéraledeLausanne Lausanne Switzerland ISSN1439-7358 ISSN2197-7100(electronic) ISBN978-3-319-10704-2 ISBN978-3-319-10705-9(eBook) DOI10.1007/978-3-319-10705-9 SpringerChamHeidelbergNewYorkDordrechtLondon LibraryofCongressControlNumber:2014956263 MathematicsSubjectClassification(2010):65XX,74XX,76XX,78XX ©SpringerInternationalPublishingSwitzerland2015 Thisworkissubjecttocopyright.AllrightsarereservedbythePublisher,whetherthewholeorpartof thematerialisconcerned,specificallytherightsoftranslation,reprinting,reuseofillustrations,recitation, broadcasting,reproductiononmicrofilmsorinanyotherphysicalway,andtransmissionorinformation storageandretrieval,electronicadaptation,computersoftware,orbysimilarordissimilarmethodology nowknownorhereafterdeveloped.Exemptedfromthislegalreservationarebriefexcerptsinconnection with reviews or scholarly analysis or material supplied specifically for the purpose of being entered and executed on a computer system, for exclusive use by the purchaser of the work. Duplication of this publication or parts thereof is permitted only under the provisions of the Copyright Law of the Publisher’slocation,initscurrentversion,andpermissionforusemustalwaysbeobtainedfromSpringer. PermissionsforusemaybeobtainedthroughRightsLinkattheCopyrightClearanceCenter.Violations areliabletoprosecutionundertherespectiveCopyrightLaw. Theuseofgeneraldescriptivenames,registerednames,trademarks,servicemarks,etc.inthispublication doesnotimply,evenintheabsenceofaspecificstatement,thatsuchnamesareexemptfromtherelevant protectivelawsandregulationsandthereforefreeforgeneraluse. While the advice and information in this book are believed to be true and accurate at the date of publication,neithertheauthorsnortheeditorsnorthepublishercanacceptanylegalresponsibilityfor anyerrorsoromissionsthatmaybemade.Thepublishermakesnowarranty,expressorimplied,with respecttothematerialcontainedherein. CoverPhoto:CampusofÉcolepolytechniquefédéraledeLausanne.BycourtesyofAlainHerzog/EPFL. Printedonacid-freepaper SpringerispartofSpringerScience+BusinessMedia(www.springer.com) Preface TheEuropeanConferenceonNumericalMathematicsandAdvancedApplications (ENUMATH) is a series of conferences held every 2 years to provide a forum for discussion on recent aspects of numerical mathematics and scientific and industrial applications. The previous ENUMATH meetings took place in Paris (1995),Heidelberg(1997),Jyvaskyla(1999),Ischia(2001),Prague(2003),Santiago deCompostela(2005),Graz(2007),Uppsala(2009),andLeicester(2011). This book containsa selection of invited and contributedlectures of the ENU- MATH2013conferenceorganisedbytheMathematicalInstituteofComputational ScienceandEngineering(MATHICSE),EPFL,andheldinLausanne,Switzerland, August 26–30, 2013. It gives an overview of recent developments in numerical analysis, computational mathematics, and applications by leading experts in the field. The conference attracted around 400 participants from around the world including11invitedtalksby: (cid:129) Ruth Elizabeth Baker (University of Oxford, UK), on “Developing multiscale modelsforexploringbiologicalphenomena” (cid:129) Eric Cancès (CERMICS, Ecole des Ponts ParisTech, France), on “Electronic structurecalculation” (cid:129) Omar Ghattas (ICES, University of Texas at Austin, USA), on “Stochastic NewtonMCMCmethodsforBayesianinverseproblemswithapplicationtoice sheetdynamics” (cid:129) Ernst Hairer (Université de Genève, Switzerland), on “Long-term analysis of numericalandanalyticaloscillations” (cid:129) Jan Hesthaven (Brown University, USA), on “High-order reduced basis multi- scalefiniteelementmethods” (cid:129) Petr Knobloch(UniverzitaKarlova,CzechRepublic),on“Finiteelementmeth- odsforconvectiondominatedproblems” (cid:129) DmitriKuzmin(FriedrichAlexanderUniversitätErlangen-Nürnberg,Germany), on“Vertex-basedlimitersforcontinuousanddiscontinuousGalerkinmethods” (cid:129) Ilaria Perugia (Università deglistudi di Pavia, Italy), on “Trefftz-discontinuous Galerkinmethodsfortime-harmonicwaveproblems” v vi Preface (cid:129) RolfStenberg(AaltoUniversity,Finland),on“Mixedfiniteelementmethodsfor elasticity” (cid:129) Martin Vetterli(EPFL, Switzerland),public lectureon “Inverseproblemsregu- larizedbysparsity” (cid:129) Barbara Wohlmuth (TU München, Germany), on “Interfaces, corner and point sources” There were 24 minisymposia and numerous contributed talks covering a broad spectrum of numerical mathematics. This ENUMATH 2013 proceeding will be useful for a wide range of readers giving them a state-of-the-art overview of advancedtechniques,algorithms,andresultsin numericalmathematicsandscien- tific computing. Advances on finite element methods, time integrators, multiscale methods,numericallinearalgebra,anddiscretisationtechniquesforfluidmechanics andopticsarepresented.Thisbookcontainsaselectionof79papersbytheinvited speakers and from the minisymposia as well as the contributed sessions. It is organisedin11partsasfollows: PartISpaceDiscretisationMethodsforPDEs PartIITimeIntegrationSchemes PartIIIAPosterioriErrorEstimationandAdaptiveMethods PartIVNumericalLinearAlgebra PartVMultiscaleModelingandSimulation PartVIReducedOrderModeling PartVIIIUncertainty,StochasticModeling,andApplications PartIXSolvers,HighPerformanceComputing,andSoftwareLibraries PartXComputationalFluidandStructuralMechanics PartXIComputationalElectromagnetics We would like to thank all the participants for their valuable contributions and scientific discussions during the conferenceand to the minisymposiumorganisers forhelpingtoshapethecorestructureofthemeeting.ThemembersoftheScientific Committee have helped us tremendously in reviewing the contributions to this proceedings. This conference would not have been possible without all the work and guidance provided by the Programme Committee: Franco Brezzi, Miloslav Feistauer, Roland Glowinski, Gunilla Kreiss, Yuri Kuznetsov, Jacques Periaux, Alfio Quarteroni, Rolf Rannacher, Endre Süli. We also thank our sponsors for theirgeneroussupport:theSchoolofBasicSciencesandtheCentreInterfacultaire Bernoulli from EPFL, the Center for Advanced Modeling Science (CADMOS), MathWorks and Springer. Last but not least we would like to acknowledge the tireless effort of Virginie Ledouble leading the administration tasks, Corinne Craman who coordinated the edition of this proceedings, all the secretaries of MATHICSEfortheirtremendoushelpinorganisingthisconference,andourPhDs andPost-Docsthathavehelpedusinmanyways. Preface vii We hope that this volume reflects the inspiring talks and lively scientific exchangesthattookplaceatEPFLduringtheENUMATH2013meeting. Lausanne,Switzerland AssyrAbdulle August2014 SimoneDeparis DanielKressner FabioNobile MarcoPicasso AlfioQuarteroni Contents PartI SpaceDiscretisationMethodsforPDEs WeaklySymmetricMixedFiniteElementsforLinearElasticity........... 3 RolfStenberg Energy-CorrectedFiniteElementMethods forScalarEllipticProblems..................................................... 19 ThomasHorger,MarkusHuber,UlrichRüde,ChristianWaluga, andBarbaraWohlmuth StabilizedGalerkinforLinearAdvectionofVectorFields.................. 37 HolgerHeumannandRalfHiptmair FiniteElement-BoundaryElement MethodsforDielectric RelaxationSpectroscopy ........................................................ 47 StephanC.KramerandGertLube OntheLocalMeshSizeofNitsche’sMethodforDiscontinuous MaterialParameters............................................................. 57 MikaJuntunen RobustLocalFluxReconstructionforVarious FiniteElementMethods......................................................... 65 RolandBecker,DanielaCapatina,andRobertLuce Onthe Use ofReconstructionOperatorsinDiscontinuous GalerkinSchemes................................................................ 75 VáclavKucˇera Adaptive DiscontinuousGalerkinMethodsfor Nonlinear Diffusion-Convection-ReactionEquations .................................... 85 BulentKarasözen,MuratUzunca,andMuratManguogˇlu ix