Focus on:
All days
Feb 4, 2016
Feb 5, 2016
All sessions
Social Dinner
Hide Contributions
Indico style
Indico style  inline minutes
Indico style  numbered
Indico style  numbered + minutes
Indico Weeks View
Back to Conference View
Choose Timezone
Use the event/category timezone
Specify a timezone
Africa/Abidjan
Africa/Accra
Africa/Addis_Ababa
Africa/Algiers
Africa/Asmara
Africa/Bamako
Africa/Bangui
Africa/Banjul
Africa/Bissau
Africa/Blantyre
Africa/Brazzaville
Africa/Bujumbura
Africa/Cairo
Africa/Casablanca
Africa/Ceuta
Africa/Conakry
Africa/Dakar
Africa/Dar_es_Salaam
Africa/Djibouti
Africa/Douala
Africa/El_Aaiun
Africa/Freetown
Africa/Gaborone
Africa/Harare
Africa/Johannesburg
Africa/Juba
Africa/Kampala
Africa/Khartoum
Africa/Kigali
Africa/Kinshasa
Africa/Lagos
Africa/Libreville
Africa/Lome
Africa/Luanda
Africa/Lubumbashi
Africa/Lusaka
Africa/Malabo
Africa/Maputo
Africa/Maseru
Africa/Mbabane
Africa/Mogadishu
Africa/Monrovia
Africa/Nairobi
Africa/Ndjamena
Africa/Niamey
Africa/Nouakchott
Africa/Ouagadougou
Africa/PortoNovo
Africa/Sao_Tome
Africa/Tripoli
Africa/Tunis
Africa/Windhoek
America/Adak
America/Anchorage
America/Anguilla
America/Antigua
America/Araguaina
America/Argentina/Buenos_Aires
America/Argentina/Catamarca
America/Argentina/Cordoba
America/Argentina/Jujuy
America/Argentina/La_Rioja
America/Argentina/Mendoza
America/Argentina/Rio_Gallegos
America/Argentina/Salta
America/Argentina/San_Juan
America/Argentina/San_Luis
America/Argentina/Tucuman
America/Argentina/Ushuaia
America/Aruba
America/Asuncion
America/Atikokan
America/Bahia
America/Bahia_Banderas
America/Barbados
America/Belem
America/Belize
America/BlancSablon
America/Boa_Vista
America/Bogota
America/Boise
America/Cambridge_Bay
America/Campo_Grande
America/Cancun
America/Caracas
America/Cayenne
America/Cayman
America/Chicago
America/Chihuahua
America/Costa_Rica
America/Creston
America/Cuiaba
America/Curacao
America/Danmarkshavn
America/Dawson
America/Dawson_Creek
America/Denver
America/Detroit
America/Dominica
America/Edmonton
America/Eirunepe
America/El_Salvador
America/Fort_Nelson
America/Fortaleza
America/Glace_Bay
America/Goose_Bay
America/Grand_Turk
America/Grenada
America/Guadeloupe
America/Guatemala
America/Guayaquil
America/Guyana
America/Halifax
America/Havana
America/Hermosillo
America/Indiana/Indianapolis
America/Indiana/Knox
America/Indiana/Marengo
America/Indiana/Petersburg
America/Indiana/Tell_City
America/Indiana/Vevay
America/Indiana/Vincennes
America/Indiana/Winamac
America/Inuvik
America/Iqaluit
America/Jamaica
America/Juneau
America/Kentucky/Louisville
America/Kentucky/Monticello
America/Kralendijk
America/La_Paz
America/Lima
America/Los_Angeles
America/Lower_Princes
America/Maceio
America/Managua
America/Manaus
America/Marigot
America/Martinique
America/Matamoros
America/Mazatlan
America/Menominee
America/Merida
America/Metlakatla
America/Mexico_City
America/Miquelon
America/Moncton
America/Monterrey
America/Montevideo
America/Montserrat
America/Nassau
America/New_York
America/Nome
America/Noronha
America/North_Dakota/Beulah
America/North_Dakota/Center
America/North_Dakota/New_Salem
America/Nuuk
America/Ojinaga
America/Panama
America/Pangnirtung
America/Paramaribo
America/Phoenix
America/PortauPrince
America/Port_of_Spain
America/Porto_Velho
America/Puerto_Rico
America/Punta_Arenas
America/Rankin_Inlet
America/Recife
America/Regina
America/Resolute
America/Rio_Branco
America/Santarem
America/Santiago
America/Santo_Domingo
America/Sao_Paulo
America/Scoresbysund
America/Sitka
America/St_Barthelemy
America/St_Johns
America/St_Kitts
America/St_Lucia
America/St_Thomas
America/St_Vincent
America/Swift_Current
America/Tegucigalpa
America/Thule
America/Tijuana
America/Toronto
America/Tortola
America/Vancouver
America/Whitehorse
America/Winnipeg
America/Yakutat
America/Yellowknife
Antarctica/Casey
Antarctica/Davis
Antarctica/DumontDUrville
Antarctica/Macquarie
Antarctica/Mawson
Antarctica/McMurdo
Antarctica/Palmer
Antarctica/Rothera
Antarctica/Syowa
Antarctica/Troll
Antarctica/Vostok
Arctic/Longyearbyen
Asia/Aden
Asia/Almaty
Asia/Amman
Asia/Anadyr
Asia/Aqtau
Asia/Aqtobe
Asia/Ashgabat
Asia/Atyrau
Asia/Baghdad
Asia/Bahrain
Asia/Baku
Asia/Bangkok
Asia/Barnaul
Asia/Beirut
Asia/Bishkek
Asia/Brunei
Asia/Chita
Asia/Choibalsan
Asia/Colombo
Asia/Damascus
Asia/Dhaka
Asia/Dili
Asia/Dubai
Asia/Dushanbe
Asia/Famagusta
Asia/Gaza
Asia/Hebron
Asia/Ho_Chi_Minh
Asia/Hong_Kong
Asia/Hovd
Asia/Irkutsk
Asia/Jakarta
Asia/Jayapura
Asia/Jerusalem
Asia/Kabul
Asia/Kamchatka
Asia/Karachi
Asia/Kathmandu
Asia/Khandyga
Asia/Kolkata
Asia/Krasnoyarsk
Asia/Kuala_Lumpur
Asia/Kuching
Asia/Kuwait
Asia/Macau
Asia/Magadan
Asia/Makassar
Asia/Manila
Asia/Muscat
Asia/Nicosia
Asia/Novokuznetsk
Asia/Novosibirsk
Asia/Omsk
Asia/Oral
Asia/Phnom_Penh
Asia/Pontianak
Asia/Pyongyang
Asia/Qatar
Asia/Qostanay
Asia/Qyzylorda
Asia/Riyadh
Asia/Sakhalin
Asia/Samarkand
Asia/Seoul
Asia/Shanghai
Asia/Singapore
Asia/Srednekolymsk
Asia/Taipei
Asia/Tashkent
Asia/Tbilisi
Asia/Tehran
Asia/Thimphu
Asia/Tokyo
Asia/Tomsk
Asia/Ulaanbaatar
Asia/Urumqi
Asia/UstNera
Asia/Vientiane
Asia/Vladivostok
Asia/Yakutsk
Asia/Yangon
Asia/Yekaterinburg
Asia/Yerevan
Atlantic/Azores
Atlantic/Bermuda
Atlantic/Canary
Atlantic/Cape_Verde
Atlantic/Faroe
Atlantic/Madeira
Atlantic/Reykjavik
Atlantic/South_Georgia
Atlantic/St_Helena
Atlantic/Stanley
Australia/Adelaide
Australia/Brisbane
Australia/Broken_Hill
Australia/Darwin
Australia/Eucla
Australia/Hobart
Australia/Lindeman
Australia/Lord_Howe
Australia/Melbourne
Australia/Perth
Australia/Sydney
Canada/Atlantic
Canada/Central
Canada/Eastern
Canada/Mountain
Canada/Newfoundland
Canada/Pacific
Europe/Amsterdam
Europe/Andorra
Europe/Astrakhan
Europe/Athens
Europe/Belgrade
Europe/Berlin
Europe/Bratislava
Europe/Brussels
Europe/Bucharest
Europe/Budapest
Europe/Busingen
Europe/Chisinau
Europe/Copenhagen
Europe/Dublin
Europe/Gibraltar
Europe/Guernsey
Europe/Helsinki
Europe/Isle_of_Man
Europe/Istanbul
Europe/Jersey
Europe/Kaliningrad
Europe/Kirov
Europe/Kyiv
Europe/Lisbon
Europe/Ljubljana
Europe/London
Europe/Luxembourg
Europe/Madrid
Europe/Malta
Europe/Mariehamn
Europe/Minsk
Europe/Monaco
Europe/Moscow
Europe/Oslo
Europe/Paris
Europe/Podgorica
Europe/Prague
Europe/Riga
Europe/Rome
Europe/Samara
Europe/San_Marino
Europe/Sarajevo
Europe/Saratov
Europe/Simferopol
Europe/Skopje
Europe/Sofia
Europe/Stockholm
Europe/Tallinn
Europe/Tirane
Europe/Ulyanovsk
Europe/Vaduz
Europe/Vatican
Europe/Vienna
Europe/Vilnius
Europe/Volgograd
Europe/Warsaw
Europe/Zagreb
Europe/Zurich
GMT
Indian/Antananarivo
Indian/Chagos
Indian/Christmas
Indian/Cocos
Indian/Comoro
Indian/Kerguelen
Indian/Mahe
Indian/Maldives
Indian/Mauritius
Indian/Mayotte
Indian/Reunion
Pacific/Apia
Pacific/Auckland
Pacific/Bougainville
Pacific/Chatham
Pacific/Chuuk
Pacific/Easter
Pacific/Efate
Pacific/Fakaofo
Pacific/Fiji
Pacific/Funafuti
Pacific/Galapagos
Pacific/Gambier
Pacific/Guadalcanal
Pacific/Guam
Pacific/Honolulu
Pacific/Kanton
Pacific/Kiritimati
Pacific/Kosrae
Pacific/Kwajalein
Pacific/Majuro
Pacific/Marquesas
Pacific/Midway
Pacific/Nauru
Pacific/Niue
Pacific/Norfolk
Pacific/Noumea
Pacific/Pago_Pago
Pacific/Palau
Pacific/Pitcairn
Pacific/Pohnpei
Pacific/Port_Moresby
Pacific/Rarotonga
Pacific/Saipan
Pacific/Tahiti
Pacific/Tarawa
Pacific/Tongatapu
Pacific/Wake
Pacific/Wallis
US/Alaska
US/Arizona
US/Central
US/Eastern
US/Hawaii
US/Mountain
US/Pacific
UTC
Save
Europe/Paris
English (United States)
Deutsch (Deutschland)
English (United Kingdom)
English (United States)
Español (España)
Français (France)
Polski (Polska)
Português (Brasil)
Türkçe (Türkiye)
Монгол (Монгол)
Українська (Україна)
中文 (中国)
Login
New Trends in Integral Equations
from
Thursday, February 4, 2016 (9:00 AM)
to
Friday, February 5, 2016 (6:00 PM)
Monday, February 1, 2016
Tuesday, February 2, 2016
Wednesday, February 3, 2016
Thursday, February 4, 2016
9:00 AM
Welcome meeting
Welcome meeting
9:00 AM  9:15 AM
Room: Amphi Becquerel
9:15 AM
Volume integral equation formulation for anisotropic elastodynamic scattering: solvability, application to smallinclusion asymptotics

Marc Bonnet
(
ENSTA
)
Volume integral equation formulation for anisotropic elastodynamic scattering: solvability, application to smallinclusion asymptotics
Marc Bonnet
(
ENSTA
)
9:15 AM  9:50 AM
Room: Amphi Becquerel
In contrast with the vast existing literature on the mathematical aspects of boundary integral equations and their application to scattering by impenetrable obstacles characterized by Dirichlet, Neumann or impedant boundary conditions, comparatively few studies are available regarding the mathematical properties of volume integral equation (VIE) formulations. This communication addresses the solvability of VIEs arising in elastodynamic scattering by penetrable obstacles. The elasticity tensor and mass density are allowed to be smoothly heterogeneous inside the obstacle and may be discontinuous across the backgroundobstacle interface, the background elastic material being homogeneous. Both materials may be anisotropic, within certain limitations for the background medium. The VIE associated with this problem is first derived, relying on known properties of the background fundamental tensor. To avoid difficulties associated with existing radiation conditions for anisotropic elastic media, we also propose a definition of the radiating character of transmission solutions. The unique solvability of the volume integral equation (and of the scattering problem) is established. For the important special case of isotropic background properties, our definition of a radiating solution is found to be equivalent to the SommerfeldKupradze radiation conditions. Moreover, usefulness of this result for the derivation and justification of smallinclusion asymptotic approximations is discussed.
9:50 AM
Hmatrix Solver for BEM: a taskbased approach

Guillaume Sylvand
(
Airbus Group
)
Hmatrix Solver for BEM: a taskbased approach
Guillaume Sylvand
(
Airbus Group
)
9:50 AM  10:25 AM
Room: Amphi Becquerel
For the numerical simulation of wave propagation in acoustics, Airbus Group Innovations relies on integral equations solved with the Boundary Elements Method. The advantages of this approach are well known: mainly accuracy, and simpler (surfacic) mesh. The main algorithm drawback is the need to cope with a dense matrix whose size can be quite large for wave propagation problems, where the mesh step is governed by the wavelength of the physical problem treated (in frequency domain). <br/> Since the late 90's, fast methods have been introduced to deal with these limitations. First, the Fast Multipole Method (FMM) allowed to compute fast matrixvector products (in O(n log2(n)) instead of O(n2) for the standard algorithm), and hence to design fast solvers using iterative methods. Lately, HMatrix methods have gained wide acceptance by introducing fast direct solvers, allowing to solve systems in O(n log2(n))  or less  without the hassle of using iterative solvers (unknown convergence rate and difficulty to find a goodpreconditionner). <br/> HMatrix [1] is a lossy, hierarchical storage scheme for matrices that, along with an associated arithmetic, provides a rich enough set of approximate operations to perform the matrix addition, multiplication, factorization (e.g. LU or LDLT ) and inversion. It relies on two core ideas : (a) nested clustering of the degrees of freedom, and of their products; and (b) adaptative compression of these clusters. Several choices exist in the literature for these two ingredients, the most common being Binary Space Partitioning for the clustering and Adaptative Cross Approximation for the compression. <br/> Together, they allow for the construction of a fast direct solver [2], which is especially important for BEM applications as it gracefully handles a large number of RightHand Sides (RHS). They also provide a kernelindependent fast solver, allowing one to use the method for dierent physics. <br/> Airbus Group Innovations has recently implemented the HMatrix arithmetic and successfully applied it to a wide range of industrial applications in electromagnetism and acoustics. Furthermore, these algorithms are hard to eciently parallelize, as the very scarce literature on the subject shows [3]. We developed a parallel solver that goes beyond the aforementioned reference, using innovative techniques on top of a stateoftheart runtime system StarPU [4][5]. This enables the solving of very large problems, with a very good eciency. In this presentation, we show some results on the accuracy of this method on several challenging applications, and its fast solving time and ecient use of resources. <br/> <br/> <br/> References <br/> <br/> [1] W. Hackbusch, A Sparse Matrix Arithmetic Based on HMatrices. Part I: Introduction to HMatrices, Computing, Volume 62, Issue 2 (1999).. <br/> <br/> [2] L. Grasedyck, W. Hackbusch, Construction and Arithmetics of HMatrices, Computing, Volume 70, Issue 4 (2003). <br/> <br/> [3] R. Kriemann, Parallel HMatrix Arithmetics on Shared Memory Systems, Computing, Volume 74, Issue 3 (2005). <br/> <br/> [4] B. Lize, Resolution directe rapide pour les elements nis de frontiere en electromagnetisme et acoustique : HMatrices. Parallelisme et applications industrielles, PhD thesis, Paris 13, 2014. <br/> <br/> [5] C. Augonnet, S. Thibault, R. Namyst, and P.A. Wacrenier, StarPU: A Unified Platform for Task Scheduling on Heterogeneous Multicore Architectures, Concurrency and Computation: Practice and Experience, Special Issue: EuroPar 2009, vol. 23, pp. 187198, Feb. 2011.
10:25 AM
Coffee break
Coffee break
10:25 AM  10:45 AM
Room: Amphi Becquerel
10:45 AM
Quasilocal multitrace formulations

Xavier Claeys
(
UPMC
)
Quasilocal multitrace formulations
Xavier Claeys
(
UPMC
)
10:45 AM  11:20 AM
Room: Amphi Becquerel
In the context of time harmonic wave scattering by piecewise homogeneous penetrable objects, we present a new variant of the multitrace boundary integral formulations (MTF). This new approach differs from the socalled "local" MTF by the presence of regularisation terms involving boundary integral operators and localised at junctions i.e. points where at least three subdomains abut. It lends itself to much more standard analysis: all operators are continuous in standard Sobolev trace spaces, a Garding inequality can be proved, which implies quasioptimal approximation property for conformal Galerkin discretisations. As regards numerical performances, this new formulation also appears slightly more accurate compared to preexisting MTF, while the speed of convergence of iterative solvers remains comparable. <br/> <br/> References: <br/> [1] X.Claeys, Quasilocal multitrace boundary integral formulations, Numer. Methods Partial Differential Equations, 31(6):2043–2062, 2015. <br/> [2] X.Claeys and R. Hiptmair and C. JerezHanckes, Multitrace boundary integral equations, chapter in Direct and Inverse Problems in Wave Propagation and Applications, 51–100, Radon Ser. Comput. Appl. Math., 14, De Gruyter, Berlin, 2013. <br/> <br/>
11:20 AM
Modelling tools for scattering and antennas problems

Jérôme Simon
(
Onera Palaiseau
)
Modelling tools for scattering and antennas problems
Jérôme Simon
(
Onera Palaiseau
)
11:20 AM  11:55 AM
Room: Amphi Becquerel
The Surface Integral Equation is one of the most used methods in the simulation of electromagnetic problems. Its implementation combined with a fast multipole algorithm (MLFMM) and iterative solvers leads efficient and accurate methods for the analysis of radar cross section (RCS) for a target. <br/> <br/> Different methods (MLFMM, ACA, Domain Decomposition Method …) and formulations EFIE (Electric Field Integral Equation), MFIE (Magnetic Field Integral Equation), CFIE (Combined Field Integral Equation) and GCSIE (integral formulations inherently wellconditioned using a regularizing operator ) were implemented at ONERA. <br/> <br/> In this talk, we will present various applications obtained with our tools like radar cross section (RCS), antenna radiation, ….
11:55 AM
Quasioptimal domain decomposition methods for timeharmonic acoustic and electromagnetic wave problems

Christophe Geuzaine
(
Université de Liège
)
Quasioptimal domain decomposition methods for timeharmonic acoustic and electromagnetic wave problems
Christophe Geuzaine
(
Université de Liège
)
11:55 AM  12:30 PM
Room: Amphi Becquerel
In terms of computational methods, solving threedimensional timeharmonic acoustic or electromagnetic wave problems is known to be challenging, especially in the high frequency regime and in the presence of inhomogenous media. The bruteforce application of the finite element method in this case leads to the solution of very large, complex and possibly indefinite linear systems. Direct sparse solvers do not scale well for such problems, and Krylov subspace iterative solvers can exhibit slow convergence, or even diverge. Domain decomposition methods provide an alternative, iterating between subproblems of smaller sizes, amenable to sparse direct solvers. In this talk I will present a class of nonoverlapping Schwarz domain decomposition methods that exhibit quasioptimal convergence properties, i.e., with a convergence that is optimal for the evanescent modes and significantly improved compared to competing approaches for the remaining modes [1, 2]. These improved properties result from a combination of an appropriate choice of transmission conditions and a suitable localization of the optimal, integral operators associated with the complementary of each subdomain [3]. The resulting algorithms are well suited for highperformance, large scale parallel computations in complex geometrical configurations when combined with appropriate preconditionners [4]. <br/><br/> [1] Y. Boubendir, X. Antoine and C. Geuzaine, A quasioptimal nonoverlapping domain decomposition algorithm for the Helmholtz equation. Journal of Computational Physics 231 (2), pp. 262280, 2012. <br/><br/> [2] M. El Bouajaji, B. Thierry, X. Antoine and C. Geuzaine. A quasioptimal domain decomposition algorithm for the timeharmonic Maxwell's equations. Journal of Computational Physics 294, pp. 3857, 2015. <br/><br/> [3] M. El Bouajaji, X. Antoine, C. Geuzaine. Approximate local magnetictoelectric surface operators for timeharmonic Maxwell's equations. Journal of Computational Physics 279 (15), 241260, 2014. <br/><br/> [4] A. Vion and C. Geuzaine. Double sweep preconditioner for optimized Schwarz methods applied to the Helmholtz problem. Journal of Computational Physics 266, 171190, 2014.
12:30 PM
Lunch
Lunch
12:30 PM  2:15 PM
Room: Club Magnan
2:15 PM
Preconditionning integral equations on the unit disc in 3d

JeanClaude Nédélec
(
Ecole Polytechnique
)
Preconditionning integral equations on the unit disc in 3d
JeanClaude Nédélec
(
Ecole Polytechnique
)
2:15 PM  2:50 PM
Room: Amphi Becquerel
2:50 PM
Analytical preconditioners for the solution of threedimensional surface scattering problems

Marion Darbas
(
Université de Picardie
)
Analytical preconditioners for the solution of threedimensional surface scattering problems
Marion Darbas
(
Université de Picardie
)
2:50 PM  3:25 PM
Room: Amphi Becquerel
The numerical solution of timeharmonic scattering problems remains challenging in the high frequency regime due to its specific computational bottlenecks. The techniques based on integral equations lead to the resolution of linear systems where the involved matrices are dense and usually badly conditioned. The improvement of these methods is a timely research area. One possibility to reduce the computational cost is to precondition iterative solvers (to speed up the convergence) and on the other hand to use fast methods to compute the matrixvector products needed at each iteration. <br/> <br/> We propose an analytical preconditioner taking inspiration of OnSurface Radiation Condition techniques. This preconditioner is an accurate approximation to the DirichlettoNeumann map. The associated integral equations are of the second kind. Moreover, the proposed preconditioner shows highly desirable advantages: sparse structure, ease of implementation and low additional computational cost. <br/> <br/> In this talk, we present first the principle of the method in the acoustic case. We show numerical simulations for various configurations. Next, we explain how to extend the approach to other types of waves, namely elastic waves (joint work with Stéphanie Chaillat et Frédérique Le Louër).
3:25 PM
Hybrid numericalasymptotic methods for wave scattering problems

Stephen Langdon
(
University of Reading
)
Hybrid numericalasymptotic methods for wave scattering problems
Stephen Langdon
(
University of Reading
)
3:25 PM  4:00 PM
Room: Amphi Becquerel
Linear wave scattering problems (e.g. for acoustic, electromagnetic and elastic waves) are ubiquitous in science and engineering applications. However, conventional numerical methods for such problems (e.g. FEM or BEM with piecewise polynomial basis functions) are prohibitively expensive when the wavelength of the scattered wave is small compared to typical lengthscales of the scatterer (the socalled "high frequency" regime). This is because the solution possesses rapid oscillations which are expensive to capture using conventional approximation spaces. In this talk we outline recent progress in the development of "hybrid numericalasymptotic" methods, which incur significantly reduced computational cost. These methods use approximation spaces containing oscillatory basis functions, carefully chosen to capture the high frequency asymptotic behaviour. In particular we discuss some of the interesting challenges arising from nonconvex, penetrable and threedimensional scatterers.
4:00 PM
Coffee break
Coffee break
4:00 PM  4:20 PM
Room: Amphi Becquerel
4:20 PM
Windowed Green function method for layeredmedia scattering

Catalin Turc
(
New Jersey Institute of Technology
)
Windowed Green function method for layeredmedia scattering
Catalin Turc
(
New Jersey Institute of Technology
)
4:20 PM  4:55 PM
Room: Amphi Becquerel
We present a new Windowed Green Function (WGF) method for the numerical integralequation solution of problems of electromagnetic scattering by obstacles in presence of dielectric or conducting halfplanes. The WGF method, which is based on use of integral kernels that can be expressed directly in terms of the freespace Green function, does not require evaluation of expensive Sommerfeld integrals. The proposed approach is fast, accurate, flexible and easy to implement. In particular straightforward modifications of existing solvers suffice to incorporate the WGF capability. The proposed method can be up to thousands of times faster, for a given accuracy, than a corresponding method based on the layerGreenfunction.
4:55 PM
Fast Boundary Elements Methods and applications

Matthieu Aussal
(
Ecole Polytechnique
)
Fast Boundary Elements Methods and applications
Matthieu Aussal
(
Ecole Polytechnique
)
4:55 PM  5:30 PM
Room: Amphi Becquerel
Fast convolution on unstructured grids have been developed for many applications (e.g. electrostatics, magnetostatics, acoustics, electro magnetics, etc.). The goal is to reduce the complexity of matrixvector products, from O(N^2) to O(N log N ). In this presentation, we describe a new efficient numerical method called Sparse Cardinal Sine Decomposition (SCSD), based on a suitable Fourier decomposition of the Green kernel, sparse quadrature formulae and TypeIII Non Uniform Fast Fourier Transforms (typeIII NUFFT). This talk summarizes this new way, provide comparisons between SCSD, FMM and HMatrix, and gives numerical results from our new Boundary Element solver, MyBEM.
7:00 PM
7:00 PM  9:00 PM
Room: Salon de marbre
Friday, February 5, 2016
9:00 AM
Connecting Integral Equation and Unified Transform Methods for Wave Scattering

Simon N. ChandlerWilde
(
University of Reading
)
Connecting Integral Equation and Unified Transform Methods for Wave Scattering
Simon N. ChandlerWilde
(
University of Reading
)
9:00 AM  9:35 AM
Room: Amphi Becquerel
In this talk we discuss the application of the unified transform method, due to A.S. Fokas and coauthors, to interior and exterior problems for time harmonic waves. Like integral equation methods, the method reduces to the solution of a problem on the boundary of the domain in question. We discuss numerical implementations, restricting the solution space to a finitedimensional subspace, and explain how some implementations can be interpreted as Galerkin methods, the convergence of which can be established by standard arguments. We focus particularly on implementations that are based on approximation by plane waves (and generalised, evanescent plane waves), and explain that these can be implemented in such a way that the numerical solution is precisely the best approximation from the plane wave subspace. In the case of scattering by diffraction gratings we note that this particular unified transform method is in fact precisely the SS* method proposed previously (though with a less precise analysis) in Arens, ChandlerWilde and De Santo (2006), where it is derived as a first kind integral equation formulation. Further details can be found in [1,2]. <br/> <br/> [1] Acoustic scattering: high frequency boundary element methods and unified transform methods. S N ChandlerWilde & S Langdon, in “Unified Transform for Boundary Value Problems: Applications and Advances”, A S Fokas & B Pelloni (editors), SIAM, 2015. <br/> <br/> [2] When all else fails, integrate by parts"  an overview of new and old variational formulations for linear elliptic PDEs. E A Spence in "Unified Transform Method for Boundary Value Problems: Applications and Advances", A S Fokas & B Pelloni (editors), SIAM, 2015.
9:35 AM
High Order Impedance Boundary Condition for 3D Integral Equations

Paul Soudais
(
Dassault Aviation
)
High Order Impedance Boundary Condition for 3D Integral Equations
Paul Soudais
(
Dassault Aviation
)
9:35 AM  10:10 AM
Room: Amphi Becquerel
Impedance Boundary Conditions (IBC) are widely used in computational electromagnetics to model thin coatings on perfectly conducting (PEC) objects. The IBC is used to reduce drastically the number of unknowns of the integral equations models and to obtain a better conditioned linear system that can be more efficiently iteratively solved. <br/> Integral equations models can be built that mix homogeneous regions, whose material characteristics are defined by their electrical permittivities and magnetic permeabilities, PEC regions and IBC. These models are very efficient for the modeling of complex structures. <br/> Usually, the impedance of a coating is assumed to be the same for all incidence angles and polarizations (standard IBC). This assumption is valid for coatings with high refractive index or significant losses. For coatings with smaller refractive index or lower losses, the accuracy of the IBC can be greatly improved by approximating it as a partial derivative equation ([1] D.J. Hoppe, Y. Rahmat Samii, Impedance Boundary Conditions in Electromagnetics, Taylor and Francis, 1995). The second order IBC can be written as: <br/> $E_t+b_1\nabla_\Gamma div E_tb_2 rot_\Gamma rot _\Gamma E_t = a_0 J + a_1 \nabla_\Gamma div J  a_2 rot_\Gamma rot_\Gamma J$ $E_t$ and $J$ are respectively the tangent electric field and the electric current. Integral equations formulations for structures with homogeneous regions, PEC regions and SIBC can be generalized to take into account second order IBC. In [1], the second order IBC is discretized with spline basis functions for 2D and body of revolution problems. Here, the second order IBC is applied to 3D problems. The currents $J$ and $M= Exn$ are discretized with the lower order HDiv functions that can be used to model objects both smooth or with sharp edges. <br/> A first set of validations will be presented that show the increased accuracy of high order IBC over standard IBC while the computational effort for solving the overall integral equations model remains similar.
10:10 AM
A new a posteriori error estimate for the BEM

Marc Bakry
Sébastien Pernet
(
Onera Toulouse
)
A new a posteriori error estimate for the BEM
Marc Bakry
Sébastien Pernet
(
Onera Toulouse
)
10:10 AM  10:45 AM
Room: Amphi Becquerel
A posteriori error estimates are tools which enable a measure of the numerical error. They are required to be equivalent to the norm of the error and to be locally computable. They are norms of values which can be computed from the numerical solution and the problem parameters. In the context of BEM, there is a loss of locality of the norms and therefore of the estimates. Standard localization techniques partially solve the issue but lead to a loss of control between the norm of the error and the estimate. <br/> We introduce a new localization technique based on the computation of the residual (the norm of the residual is an equivalent norm of the error). By applying a wellchosen isomorphism to the residual, we can "carry" it in some functional space where the norm is local (typically L2). We rst introduce the concept for the BEM in 3Dacoustics, then we give an introduction on what this estimation would look like for the EFIE.
10:45 AM
Coffee break
Coffee break
10:45 AM  11:05 AM
Room: Amphi Becquerel
11:05 AM
Solution schemes and approximation tools devoted to the simulation of electromagnetic testing by the boundary element method

Edouard Demaldent
(
CEA
)
Solution schemes and approximation tools devoted to the simulation of electromagnetic testing by the boundary element method
Edouard Demaldent
(
CEA
)
11:05 AM  11:40 AM
Room: Amphi Becquerel
Electromagnetic testing is widely used for the characterization of a medium as to the detection of defects. In particular, the eddy current nondestructive testing of tubes in steam generators is determinant to diagnose the integrity of heat exchangers in nuclear industry. A valuable support to the mastering of these processes is brought by modeling and the finite boundary element method (BEM) is an appropriate simulation tool to many inspection configurations. The department of imaging and simulation for the control at CEA LIST is developing a BEM code devoted to these applications, mostly for eddy current testing. Some of these tools are, or will be, integrated into the CIVA software platform, whose target users are experts in nondestructive testing (non numericians). <br/><br/> In this talk, we will present our research work and review the technical options related to this activity. We will start with an overview of a recent study carried out in collaboration with IRSN (the French public expert in nuclear and radiological risks) to illustrate our practical use of BEM. We will then introduce low frequency formulations that are studied in collaboration with the research group POEMS. We will in particular discuss a multistep algorithm for solving the transmission problem (known as PMCHWT), which is stable over a wide range of parameters that are relevant to electromagnetic testing, and in particular go beyond eddy currents. The third part of the talk will focus on the discretization tools developed for the BEM code at LIST. They are based on the use of basic interpolation techniques to simplify the construction of complex approximation spaces that meet our needs for light but accurate computations (such as the Helmholtz decomposition of highorder edge functions).
11:40 AM
A low dispersive Trefftz DG method based on a local BEM for the solution of the Helmholtz equation

Abderrahmane Bendali
(
INSA Toulouse
)
A low dispersive Trefftz DG method based on a local BEM for the solution of the Helmholtz equation
Abderrahmane Bendali
(
INSA Toulouse
)
11:40 AM  12:15 PM
Room: Amphi Becquerel
We first briefly present the dispersion phenomenon and show how it damages the numerical solution of wave problems over large distances of propagation. We then introduce a Trefftz Discontinuous Galerkin (TDG) symmetric formulation for the Helmholtz equation with piecewise constant coefficients. Recall that Trefftz methods are discretization processes based on the use of exact interior either local or global solutions. The Trefftz method considered in this work is based on using locally a Boundary Element Method (BEM), thus avoiding the restriction on the type of the local waves used in the usual TDG methods. We show then that accurate local approximations of the DirichlettoNeumann map have a direct impact on the reduction of the dispersion error. We then present some numerical tests bringing out that the obtained procedure can completely rub out the dispersion error contrary to the best Interior Penalty DG (IPDG) methods.
12:15 PM
Lunch
Lunch
12:15 PM  2:00 PM
Room: Club Magnan
2:00 PM
SecondKind Single Trace Boundary Integral Equations

Ralf Hiptmair
(
ETH Zürich
)
SecondKind Single Trace Boundary Integral Equations
Ralf Hiptmair
(
ETH Zürich
)
2:00 PM  2:35 PM
Room: Amphi Becquerel
For secondorder linear transmission problems involving a single closed interface separating two homogeneous materials, a wellposed secondkind boundary integral formulation has been known for a long time. It arises from a straightforward combination of interior and exterior Calderon identities. Apparently, this simple approach cannot be extended to "composite" settings involving more than two materials. <br/> <br/> The key observation is that the same secondkind boundary integral equations (BIE) can also be obtained through a multipotential representation formula. We can attach a potential to each boundary of a material subdomain, add them all up to a multi potential, and then we notice that, thanks to a nulleld property, the sum provides a representation of the field solution, when its traces a plugged into the potentials. Taking traces yields a BIE on the skeleton of the subdomain partition. The skeleton traces of the unknown field will solve it. <br/> <br/> Using the fact that multipotentials for a single homogeneous material must vanish, the BIE can be converted into secondorder form: for the scalar case (acoustics) its operator becomes a compact perturbation of the identity in $L^2$. Galerkin matrices arising from piecewise polynomial Galerkin boundary element (BEM) discretization will be intrinsically wellconditioned. <br/> <br/> The new secondkind boundary element method has been implemented both for acoustic and electromagnetic scattering at composite objects. Numerical tests confirm the excellent meshsize independent conditioning of the Galerkin BEM matrices and the resulting fast convergence of iterative solvers like GMRES. Furthermore, by simple postprocessing, we obtain discrete solutions of competitive accuracy compared to using BEM with the standard firstkind BIE. Wellposedness of the new secondkind formulations is an open problem, as is the compactness of the modulation of the identity in the case of Maxwell's equations. Reassuringly, computations have never hinted at a lack of stability. <br/> <br/> <br/> References <br/> <br/> [1] X. Claeys, R. Hiptmair, and E. Spindler, Secondkind boundary integral equa tions for scattering at composite partly impenetrable objects, Tech. Rep. 201519, Seminar for Applied Mathematics, ETH Zurich, Switzerland, 2015. Submitted to BIT. <br/> <br/> [2] X. Claeys, R. Hiptmair, and E. Spindler, A secondkind galerkin boundary element method for scattering at composite objects, BIT Numerical Mathematics, 55 (2015), pp. 3357. <br/> <br/>
2:35 PM
New results in BEM for Electromagnetic Compatibility

Toufic Abboud
(
IMACS
)
New results in BEM for Electromagnetic Compatibility
Toufic Abboud
(
IMACS
)
2:35 PM  3:10 PM
Room: Amphi Becquerel
Electromagnetic Compatibility (EMC) modelling is a major issue in aeronautic industry, with applications like Lightning Indirect Effects, HighIntensity Radiated Fields, antenna coupling, wifi… These simulations represent many challenges for the numerical solvers: very large frequency band (starting from DC), need for local models (wires, slots, equipment…), low levels of currents and fields inside cavities, presence of geometric and data singularities and the need of robust adaptive meshing procedures… and finally the model size requires efficient fast and parallel solvers. We present new results on these issues.
3:10 PM
Boundary Element Solutions of Electromagnetic Problems Based on Helmholtz Projectors

Francesco Andriulli
(
Télécom Bretagne
)
Boundary Element Solutions of Electromagnetic Problems Based on Helmholtz Projectors
Francesco Andriulli
(
Télécom Bretagne
)
3:10 PM  3:45 PM
Room: Amphi Becquerel
Integral equation solvers are widely used for simulating electromagnetic scattering and radiation from metallic and penetrable objects. Long popular in academic circles, these solvers have been in recent years incorporated into several commercial electromagnetic analysis and design tools, after the advent of fast multipole and related algorithms. The effectiveness of these solvers notwithstanding, boundary element methods are often plagued by several issues related to illconditioning and numerical instabilities both for low and high frequencies. <br/> This talk will focus on a new family of schemes which can effectively solve several of these issues. These solvers are based on implicit discrete Helmholtz decompositions obtained via suitably defined projectors. The use of these projectors, we will briefly delineate, allows to obtain electric/magnetic/combined Calderon equations with peculiarly favorable properties (both in frequency and time domain), efficient wavelet preconditioners, as well as Calderonlike approaches that do not require the use of dual boundary elements.
3:45 PM
Coffee break
Coffee break
3:45 PM  4:05 PM
Room: Amphi Becquerel
4:05 PM
Fast boundary element methods in industrial applications

Olaf Steinbach
Fast boundary element methods in industrial applications
Olaf Steinbach
4:05 PM  4:40 PM
Room: Amphi Becquerel
In this talk we consider direct and indirect boundary integral formulations for the solution of electro– and magnetostatic field computations in industrial applications. Particular interest is on the use of floating potentials in a multi–dielectric setting with high permittivity. <br/> This talk summarizes joint work with A. Blaszczyk and Z. Andjelic (ABB), and with D. Amann, G. Of, and P. Urthaler (TU Graz). <br/> <br/> <br/> References <br/> <br/> [1] Z. Andjelic, G. Of, O. Steinbach, P. Urthaler: Fast boundary element methods for industrial applications in magnetostatics. In: Fast Boundary Element Methods in Engineering and Industrial Applications (U. Langer, M. Schanz, O. Steinbach, W. L. Wendland eds.), Lecture Notes in Applied and Computational Mechanics, vol. 63, Springer, Heidelberg, pp. 111–143, 2012. <br/> <br/> [2] Z. Andjelic, G. Of, O. Steinbach, P. Urthaler: Boundary element methods for magnetostatic field problems: A critical view. Comput. Visual. Sci. 14 (2011) 117–130. <br/> <br/> [3] D. Amann, A. Blaszczyk, G. Of, O. Steinbach: Simulation of floating potentials in industrial applications by boundary element methods. J. Math. Ind. 4:13 (2014) 15p.
4:40 PM
An Hmatrix based direct solver for the Boundary Element Method in 3D elastodynamics

Stéphanie Chaillat
An Hmatrix based direct solver for the Boundary Element Method in 3D elastodynamics
Stéphanie Chaillat
4:40 PM  5:15 PM
Room: Amphi Becquerel
The main advantage of the Boundary Element Method (BEM) is that only the domain boundaries are discretized leading to a drastic reduction of the total number of degrees of freedom. <br/> In traditional BE implementation the dimensional advantage with respect to domain discretization methods is offset by the fullypopulated nature of the BEM coefficient matrix. In the present work, we propose a fast method to solve the BEM system in 3D frequencydomain elastodynamics. <br/> Using the Hmatrix arithmetic and lowrank approximations (performed with Adaptive Cross Approximation), we derive a fast direct solver. We assess the numerical efficiency and accuracy on the basis of numerical results obtained for problems having known solutions. In particular, we study the efficiency of lowrank approximations when the frequency is increased. The efficiency of the method is also illustrated to study seismic wave propagation in 3D domains. <br/><br/> This is a joint work with Patrick Ciarlet and Luca Desiderio