We consider the problem of quantifying uncertainty regarding the output of an electromagnetic field problem in the presence of a large number of uncertain input parameters. In order to reduce the growth in complexity with the number of dimensions, we employ a dimensionadaptive stochastic collocation method based on nested univariate nodes. We examine the accuracy and performance of collocation schemes based on Clenshaw-Curtis and Leja rules, for the cases of uniform and bounded, non-uniform random inputs, respectively. Based on numerical experiments with an academic electromagnetic field model, we compare the two rules in both the univariate and multivariate case and for both quadrature and interpolation purposes. Results for a real-world electromagnetic field application featuring high-dimensional input uncertainty are also presented.keywordsdimension adaptivity, Clenshaw-Curtis, computational electromagnetics, Leja, sparse grids, stochastic collocation, uncertainty quantification.

This work presents a data-driven magnetostatic finite-element solver that is specifically well suited to cope with strongly nonlinear material responses. The data-driven computing framework is essentially a multiobjective optimization procedure matching the material operation points as closely as possible to given material data while obeying Maxwell's equations. Here, the framework is extended with heterogeneous (local) weighting factors-one per finite element-equilibrating the goal function locally according to the material behavior. This modification allows the data-driven solver to cope with unbalanced measurement data sets, that is, data sets suffering from unbalanced space filling. This occurs particularly in the case of strongly nonlinear materials, which constitute problematic cases that hinder the efficiency and accuracy of standard data-driven solvers with a homogeneous (global) weighting factor. The local weighting factors are embedded in the distance-minimizing data-driven algorithm used for noiseless data, likewise for the maximum entropy data-driven algorithm used for noisy data. Numerical experiments based on a quadrupole magnet model with a soft magnetic material show that the proposed modification results in major improvements in terms of solution accuracy and solver efficiency. For the case of noiseless data, local weighting factors improve the convergence of the data-driven solver by orders of magnitude. When noisy data are considered, the convergence rate of the data-driven solver is doubled.

This paper develops a data-driven magnetostatic finite-element (FE) solver which directly exploits measured material data instead of a material curve constructed from it. The distances between the field solution and the measurement points are minimized while enforcing Maxwell's equations. The minimization problem is solved by employing the Lagrange multiplier approach. The procedure wraps the FE method within an outer data-driven iteration. The method is capable of considering anisotropic materials and is adapted to deal with models featuring a combination of exact material knowledge and measured material data. Thereto, three approaches with an increasing level of intrusivity according to the FE formulation are proposed. The numerical results for a quadrupole-magnet model show that data-driven field simulation is feasible and affordable and overcomes the need of modeling the material law.

We propose an efficient surrogate modeling technique for uncertainty quantification. The method is based on a wellknown dimension-adaptive collocation scheme. We improve the scheme by enhancing sparse polynomial surrogates with conformal maps and adjoint error correction. The methodology is applied to Maxwell's source problem with random input data. This setting comprises many applications of current interest from computational nanoplasmonics, such as grating couplers or optical waveguides. Using a nontrivial benchmark model, we show the benefits and drawbacks of using enhanced surrogate models through various numerical studies. The proposed strategy allows us to conduct a thorough uncertainty analysis, taking into account a moderately large number of random parameters.

The design of intrinsically flat two-dimensional optical components, i.e., metasurfaces, generally requires an extensive parameter search to target the appropriate scattering properties of their constituting building blocks. Such design methodologies neglect important near-field interaction effects, playing an essential role in limiting the device performance. Optimization of transmission, phase-addressing and broadband performances of metasurfaces require new numerical tools. Additionally, uncertainties and systematic fabrication errors should be analysed. These estimations, of critical importance in the case of large production of metaoptics components, are useful to further project their deployment in industrial applications. Here, we report on a computational methodology to optimize metasurface designs. We complement this computational methodology by quantifying the impact of fabrication uncertainties on the experimentally characterized components. This analysis provides general perspectives on the overall metaoptics performances, giving an idea of the expected average behavior of a large number of devices.