High order and energy preserving discontinuous Galerkin methods for the Vlasov-Poisson system

We present a computational study for a family of discontinuous Galerkin methods for the one dimensional Vlasov-Poisson system that has been recently introduced. We introduce a slight modification of the methods to allow for feasible computations while preserving the properties of the original methods. We study numerically the verification of the theoretical and convergence analysis, discussing also the conservation properties of the schemes. The methods are validated through their application to some of the benchmarks in the simulation of plasma physics.Comment: 44 pages, 28 figure

A BEM based on the Bézier/Bernstein polynomial for acoustic waveguide modelization

42nd International Conference on Boundary Elements and other Mesh Reduction Methods, BEM/MRM 2019; ITeCons-University of CoimbraCoimbra; Portugal; 2 July 2019 through 4 July 2019; Code 155806. Publicado en WIT Transactions on Engineering Sciences, Vol 126This paper proposes a novel boundary element approach formulated on the Bézier–Bernstein basis to yield a geometry-independent field approximation. The proposed method is geometrically based on both computer aided design (CAD) and isogeometric analysis (IGA), but field variables are independently approximated from the geometry. This approach allows the appropriate approximation functions for the geometry and variable field to be chosen. We use the Bézier–Bernstein form of a polynomial as an approximation basis to represent both geometry and field variables. The solution of the element interpolation problem in the Bézier–Bernstein space defines generalised Lagrange interpolation functions that are used as element shape functions. The resulting Bernstein–Vandermonde matrix related to the Bézier–Bernstein interpolation problem is inverted using the Newton–Bernstein algorithm. The applicability of the proposed method is demonstrated by solving the Helmholtz equation over an unbounded region in a two-and-a-half dimensional (2.5D) domain.Ministerio de Economía y Competitividad BIA2016-75042-C2-1-

Conformally Mapped Polynomial Chaos Expansions for Maxwell's Source Problem with Random Input Data

Generalized Polynomial Chaos (gPC) expansions are well established for forward uncertainty propagation in many application areas. Although the associated computational effort may be reduced in comparison to Monte Carlo techniques, for instance, further convergence acceleration may be important to tackle problems with high parametric sensitivities. In this work, we propose the use of conformal maps to construct a transformed gPC basis, in order to enhance the convergence order. The proposed basis still features orthogonality properties and hence, facilitates the computation of many statistical properties such as sensitivities and moments. The corresponding surrogate models are computed by pseudo-spectral projection using mapped quadrature rules, which leads to an improved cost accuracy ratio. We apply the methodology to Maxwell's source problem with random input data. In particular, numerical results for a parametric finite element model of an optical grating coupler are given