ICOSAHOM 2020
International Conference on Spectral and High Order Methods
12th  16th July 2021  Vienna, Austria
Conference Agenda
Overview and details of the sessions of this conference. Please select a date or location to show only sessions at that day or location. Please select a single session for detailed view (with abstracts and downloads if available).
Please note that all times are shown in the time zone of the conference. The current conference time is: 30th Nov 2022, 09:09:23pm CET

Session Overview 
Session  
Z9a: wave propagation problems
 
Presentations  
4:00pm  4:20pm
ComplexScaled Infinite Elements for Resonance Problems in Open Systems ^{1}ENSTA Paris; ^{2}TU Wien Complex scaling is a popular method to treat scattering and resonance problems in open domains. Thereby the unbounded domain is decomposed into a bounded interior and an unbounded exterior part. Subsequently the technique of complex scaling is applied to the exterior domain to obtain exponentially decreasing solutions. Finally, when using socalled perfectly matched layers (PMLs), the complexscaled exterior is usually truncated and discretized using finite elements or finite differences. Contrary to the use of PMLs we suggest the use of complexscaled infinite elements which are closely related to Hardy space infinite elements. For discretizing the exterior complexscaled problems we use a tensor product method based on the decomposition of the exterior domain into a surface and a generalized radial coordinate. The surface component is discretized using standard finite element basis functions while for the discretization of the radial component we use certain spectral ansatz functions. Due to this ansatz we evade introducing the error due to the truncation of the exterior domain and obtain superalgebraic approximation properties. Moreover, we are able to treat a class of inhomogeneous exterior domains which is not straightforward when using Hardy space infinite elements. 4:20pm  4:40pm
Asymptotic Green's function methods for wave propagations Iowa State University, United States of America We will present asymptotic methods for simulating high frequency wave propagation. The Huygens' principle or the Feymann's path integral is used as the time propagator for the wavefunction, where the Green's functions are approximated by asymptotic approximations. Upon obtaining analytic approximations for the phase and amplitude of the Green's functions, the resulting integral can be evaluated by fast Fourier transform after appropriate lowrank approximations. The perfectly matched layer method is incorporated to restrict the computation onto a bounded domain of interest. Numerical example will be presented to demonstrate the asymptotic Green's function methods. 4:40pm  5:00pm
A matrixfree highorder spectral implementation of the acoustic wave equation with a perfectly matched layer. ^{1}Universidade de São Paulo, Brazil; ^{2}Imperial College London This work studies the application of a matrixfree method, using sumfactorization, to a highorder spectral element implementation of the acoustic wave problem with a perfectly matched layer. Matrixfree methods have shown to be effective for increasing computational performance in highorder problems that require parallel scalability. Using spectral elements we discuss the computational cost and memory storage at different polynomial degrees and in different architectures. In this talk we will also compare the matrixfree highorder spectral element implementation in quadrilaterals with a matrixfree masslumped higherorder simplicial finite element wave propagation. The acoustic wave equation is solved using a bidimensional and tridimensional mesh and a perfectly matched layer is applied in order to reduce reflections at the boundary. The added PerfectlyMatching Layer increases the degrees of freedom of the problem by adding a vectorvalued auxiliary variable, in the twodimensional case, or a vectorvalued and scalar auxiliary variables, in the tridimensional case, therefore increasing the computational cost. The problem was solved using a open source coding package called spyro which utilizes the Firedrake framework and a discretization relevant to the forward part of Full Waveform Inversion in active source seismic imaging. This research was carried out in association with the ongoing R\&D project registered as ANP 207142, “Software technologies for modeling and inversion, with applications in seismic imaging” (University of São Paulo / Shell Brasil / ANP). 5:00pm  5:20pm
Divergence error based adaptivity in discontinuous Galerkin solution of timedomain Maxwell’s equations ^{1}Dept. of Aerospace Engineering, Indian Institute of Technology Bombay; ^{2}Dept. of Mechanical Engineering, Indian Institute of Technology Bombay A variety of physical phenomena are modelled as systems of partial differential equations that admit divergencefree solutions. In some of these, like the incompressible Euler and NavierStokes equations, this condition of the solution being divergencefree is enforced explicitly. In certain other systems like the timedomain Maxwell’s equations, the usual practice is to incorporate the solenoidal condition within the evolution equations, combined with the requirement that the initial conditions be solenoidal. It relies on the reasoning that if field variables are initially divergencefree, they remain so when evolved in time using the first order divcurl equations. Jiang et. al. [1] have challenged the arguments that this system is overdetermined and satisfying the solenoidal condition initially, ensures that it is satisfied at all times. Setting aside issues of formulation, not all numerical schemes satisfy the solenoidal condition. The finite difference time domain (FDTD) method proposed by Yee [2], satisfies the divergencefree condition by design. It uses a grid where different components of the field variables are computed at staggered spatial and temporal points. Advancements in higher order Godunov schemes for problems in computational electromagnetism (CEM) gave rise to the finite volume time domain (FVTD) [3] and discontinuous Galerkin time domain (DGTD) methods [4]. These schemes do not account for the divergence constraint. In literature, there are various approaches to meet the constraint imposed by Gauss’ law. A divergence cleaning step is often added that solves a Poisson equation for a correction potential. Assous [5] used a constrained variational formulation of Maxwell’s equations and applied a penalization technique. In [6], Munz et. al. reformulated the constrained Maxwell’s equations and introduced a coupling term into Gauss’ law, rendering a perfectly hyperbolic system of equations. This made for a natural extension of the explicit methods for Maxwell’s equations to a purely hyperbolic system. In DGTD, with standard piecewise polynomial spaces used and no dedicated measures for constraint preservation taken, it is observed that global divergence errors are kth order small when using polynomial bases of degree k to represent the solution [7]. Divergence errors accruing in conservative higherorder formulations do not significantly impact the overall accuracy of the solution [8] and are often disregarded in practice. In [8], Cioni et. al used a mixed finite volume/finite element method to show that divergence error, despite being linked to the accuracy of the solver and the underlying discretization, does not hamper the formal accuracy of the solution. In this paper, we propose another point of view, that of constructively using this error in divergence to improve spatial accuracy, rather than of either ignoring or eliminating the naturally occurring divergence error in conservative, non FDTD frameworks. We establish that relative divergence and relative truncation errors are related and propose that divergence error can be used as an effective truncation error indicator. Since solving the evolution equations does not decrease divergence error in computations, it constantly tracks truncation error. The proposed divergence based error indicator may be utilized to drive adaptive methods that assign spatial operators of varying accuracy in the computational domain, with the motive of achieving desired levels of accuracy using fewer degrees of freedom. We begin with the transverse magnetic (TM) mode of the timedomain Maxwell’s equations to formulate the formal relation between the relative truncation error in the residual appearing in the semidiscrete system and the propagation of the associated relative divergence errors. Here, relative refers to the difference between quantities computed using different discretizations [9]. We extend this definition to incorporate it in a padaptive DGTD framework. Different levels of discretizations are obtained by operators formed using polynomial bases of varying degrees. The formulated relation is applied on plane wave solutions and the resultant simplified expressions obtained, are verified by solving canonical problems using DGTD. The correlation between discretization error in solutions and their divergence is also shown with numerical examples. References [1] Bo Nan Jiang, Jie Wu, and L. A. Povinelli. “The origin of spurious solutions in computational electromagnetics”. In: Journal of Computational Physics 125.1 (1996), pp. 104–123. issn: 00219991. doi: 10.1006/jcph.1996.0082. [2] Kane Yee. “Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media”. In: IEEE Transactions on Antennas and Propagation 14.3 (1966), pp. 302–307. issn: 0018926X. doi: 10.1109/TAP.1966.1138693. url: http : / / ieeexplore.ieee.org/document/1138693/. [3] Vijaya Shankar, Alireza H Mohammadian, and William F Hall. “A Time Domain, FiniteVolume Treatment for the Maxwell Equations”. In: Electromagnetics 10.12 (1990), pp. 127–145. doi: 10.1080/02726349008908232. url: https://doi.org/10.1080/02726349008908232. [4] Jan S. Hesthaven and Tim Warburton. Nodal Discontinuous Galerkin Methods. Vol. 54. Texts in Applied Mathematics. New York, NY: Springer New York, 2008. isbn: 9780387720654. doi: 10.1007/9780387720678. url: http://link.springer.com/10.1007/9780387720678. [5] F. Assous et al. “On a FiniteElement Method for Solving the ThreeDimensional Maxwell Equations”. In: Journal of Computational Physics 109.2 (1993), pp. 222–237. issn: 00219991. doi: 10.1006/jcph.1993.1214. [6] C.D. Munz et al. “Divergence Correction Techniques for Maxwell Solvers Based on a Hyperbolic Model”. In: Journal of Computational Physics 161.2 (2000), pp. 484–511. issn: 00219991. doi: 10.1006/jcph.2000.6507. url: https://www.sciencedirect.com/science/article/pii/S0021999100965070. [7] Bernardo Cockburn, Fengyan Li, and Chi Wang Shu. “Locally divergencefree discontinuous Galerkin methods for the Maxwell equations”. In: Journal of Computational Physics 194.2 (2004), pp. 588–610. issn: 00219991. doi: 10.1016/j.jcp.2003.09.007. [8] J.P. Cioni, Loula Fezoui, and Herve Steve. “A Parallel TimeDomain Maxwell Solver Using Upwind Schemes and Triangular Meshes”. In: IMPACT of Computing in Science and Engineering 5.3 (1993), pp. 215–247. issn: 08998248. doi: 10.1006/icse.1993.1010. url: https://linkinghub.elsevier.com/retrieve/pii/S0899824883710104. [9] Avijit Chatterjee. “A Multilevel Numerical Approach with Application in TimeDomain Electromagnetics”. In: Communications in Computational Physics 17.3 (2015), pp. 703–720. issn: 18152406. doi: 10.4208/cicp.181113.271114a. url: https://www.cambridge.org/core/product/ identifier/S1815240615000158/type/journal{\_}article. 5:20pm  5:40pm
HighOrder Finite Element Simulations of Optical Fiber Amplifiers ^{1}The University of Texas at Austin, United States of America; ^{2}Air Force Research Laboratory, United States of America Fiber laser amplifiers are of interest in communication technology, medical applications, and military defense capabilities. Silica fiber amplifiers can achieve highpower operation with great efficiency. At high optical intensities, multimode amplifiers suffer from undesired thermal coupling effects which pose a major obstacle in powerscaling of such devices. In this talk, both modeling and computational advancements to a unique threedimensional finite element (FE) model for the simulation of laser amplification in a fiber amplifier are presented. This model is based on the timeharmonic Maxwell equations for two weakly coupled electromagnetic fields, and it incorporates thermal effects via coupling with the heat equation. The highfrequency nature of the wave propagation problem requires the use of highorder discretizations to effectively counter numerical pollution. The discontinuous PetrovGalerkin (DPG) FE method provides a stable discretization with a builtin error indicator. For simulating a significant fiber length, a scalable parallel implementation is critical. The model is parallelized with an MPI/OpenMP implementation, based on a parallel nested dissection solver suited for the DPG linear system. Weak scalability for the model is shown with up to 24,576 cores on manycore compute architectures, enabling the solution of the nonlinear 3D fiber model with more than 8,000 wavelengths. 5:40pm  6:00pm
An Energy Conservative High Order Method for Liouville’s Equation of Geometrical Optics ^{1}Eindhoven University of Technology, Netherlands, The; ^{2}Signify Traditionally, in geometrical optics the illuminance or luminous intensity at a target is computed using Monte Carlo ray tracing. This approach exhibits slow convergence, and in general does not obey energy conservation when using forward ray tracing. An alternative approach is based on Liouville’s equation which states the conservation of energy and describes the advection of basic luminance in phase space. The basic luminance is roughly the luminous flux per unit area, per unit solid angle, and phase space describes the set of positions and direction coordinates. From the basic luminance integrated quantities, such as the illuminance or luminous intensity, can be computed. Ray tracing can directly compute these integrated quantities, whereas for Liouville’s equation one has to compute the basic luminance on phase space first and subsequently integrate to obtain these quantities. Despite the increased dimensionality compared to ray tracing, Liouville’s equation can be solved using a high order numerical method and, therefore, potentially yield faster time to error. At an optical interface the laws of optics, such as specular reflection and Snell’s law, have to be applied. These laws describe a discontinuous change in the direction coordinate of a ray, i.e., a jump in phase space. Hence, the basic luminance is discontinuous across an optical interface. A suitable method to solve Liouville’s equation is the discontinuous Galerkin spectral element method (DGSEM). At an optical interface nonlocal boundary conditions arise as a consequence of the laws of optics. Moreover, these nonlocal boundary conditions must also satisfy energy conservation constraints. Due to these nonlocal boundary conditions the elements close to the optical interface are connected in a nontrivial manner. A method was developed to handle optical interfaces for twodimensional optics ensuring energy conservation up to machine precision. In an example numerical experiments show that the high convergence rate of the DGSEM is preserved. Moreover, the DGSEM is compared to ray tracing. The results show that the DGSEM significantly outperforms ray tracing for time to error, yielding a much faster convergence to accurate solutions. 
Contact and Legal Notice · Contact Address: Privacy Statement · Conference: ICOSAHOM2020 
Conference Software  ConfTool Pro 2.6.145+CC © 2001–2022 by Dr. H. Weinreich, Hamburg, Germany 