Source author record

Herbert De Gersem

Herbert De Gersem appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2023arXiv

A Stabilized Circuit-Consistent Foil Conductor Model

The magnetoquasistatic simulation of large power converters, in particular transformers, requires efficient models for their foils windings by means of homogenization techniques. In this article, the classical foil conductor model is derived and an inconsistency in terms of circuit theory is observed, which may lead to time-stepping instability. This can be related to the differential-algebraic nature of the resulting system of equations. It is shown how the foil conductor model can be adapted to mitigate this problem by a modified definition of the turn-by-turn conductance matrix. Numerical results are presented to demonstrate the instability and to verify the effectiveness of the new adapted foil conductor model.

preprint2022arXiv

A Neural Solver for Variational Problems on CAD Geometries with Application to Electric Machine Simulation

This work presents a deep learning-based framework for the solution of partial differential equations on complex computational domains described with computer-aided design tools. To account for the underlying distribution of the training data caused by spline-based projections from the reference to the physical domain, a variational neural solver equipped with an importance sampling scheme is developed, such that the loss function based on the discretized energy functional obtained after the weak formulation is modified according to the sample distribution. To tackle multi-patch domains possibly leading to solution discontinuities, the variational neural solver is additionally combined with a domain decomposition approach based on the Discontinuous Galerkin formulation. The proposed neural solver is verified on a toy problem and then applied to a real-world engineering test case, namely that of electric machine simulation. The numerical results show clearly that the neural solver produces physics-conforming solutions of significantly improved accuracy.

preprint2022arXiv

Assessing the Performance of Leja and Clenshaw-Curtis Collocation for Computational Electromagnetics with Random Input Data

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 dimension-adaptive 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.

preprint2022arXiv

Electromagnetic Quasistatic Field Formulations of Darwin Type

Electromagnetic quasistatic (EMQS) fields, where radiation effects are neglected, while Ohmic losses and electric and magnetic field energies are considered, can be modeled using Darwin-type field models as an approximation to the full Maxwell equations. Commonly formulated in terms of magnetic vector and electric scalar potentials, these EMQS formulations are not gauge invariant. Several EMQS formulations resulting from different gauge equations are considered and analyzed in terms of their structural properties and their modeling capabilities and limitations. Associated discrete field formulations in the context of the Maxwell-grid equations of the Finite Integration Technique are considered in frequency and time domain and are studied with respect to their algebraic properties. A comparison of numerical simulation results w.r.t. reference solutions obtained with established formulations for the full Maxwell equations are presented.

preprint2022arXiv

Power Module Heat Sink Design Optimization with Ensembles of Data-Driven Polynomial Chaos Surrogate Models

We consider the problem of optimizing the design of a heat sink used for cooling an insulated gate bipolar transistor (IGBT) power module. The thermal behavior of the heat sink is originally estimated using a high-fidelity computational fluid dynamics (CFD) simulation, which renders numerical optimization too computationally demanding. To enable optimization studies, we substitute the CFD simulation model with an inexpensive polynomial surrogate model that approximates the relation between the device's design features and a relevant thermal quantity of interest. The surrogate model of choice is a data-driven polynomial chaos expansion (DD-PCE), which learns the aforementioned relation by means of polynomial regression. Advantages of the DD-PCE include its applicability in small-data regimes and its easily adaptable model structure. To address the issue of model-form uncertainty and model robustness in view of limited training and test data, ensembles of DD-PCEs are generated based on data re-shuffling. Then, using the full ensemble of surrogate models, the surrogate-based predictions are accompanied by uncertainty metrics such as mean value and variance. Once trained and tested in terms of accuracy and robustness, the ensemble of DD-PCE surrogates replaces the high-fidelity simulation model in optimization algorithms aiming to identify heat sink designs that optimize the thermal behavior of the IGBT under geometrical and operational constraints. Optimized heat sink designs are obtained for a computational cost much smaller than utilizing the original model in the optimization procedure. Due to ensemble modeling, the optimization results can also be assessed in terms of uncertainty and robustness. Comparisons against alternative surrogate modeling techniques illustrate why the DD-PCE should be preferred in the considered setting.

preprint2020arXiv

A Coupled A-H Formulation for Magneto-Thermal Transients in High-Temperature Superconducting Magnets

The application of high-temperature superconductors to accelerator magnets for future particle colliders is under study. Numerical methods are crucial for an accurate evaluation of the complex dynamical behavior of the magnets, especially concerning the magnetic field quality and thermal behavior. We present a coupled $\textbf{A}$-$\textbf{H}$ field formulation for the analysis of magneto-thermal transients in accelerator magnets. The magnetic field strength $\textbf{H}$ accounts for the eddy current problem in the source regions containing the superconducting domains, while the magnetic vector potential $\textbf{A}$ represents the magnetoquasistatic problem in the normal and non-conducting domains. Furthermore, we include a slab approximation for the source regions, making the formulation suitable for large scale models composed of thousands of tapes. In this work, the relevant equations are derived and discussed, with emphasis on the coupling conditions. The weak formulation is derived, and numerical results are provided in order to both, verify the formulation and scale it to the size of an accelerator magnet.

preprint2020arXiv

Approximation and Uncertainty Quantification of Systems with Arbitrary Parameter Distributions using Weighted Leja Interpolation

Approximation and uncertainty quantification methods based on Lagrange interpolation are typically abandoned in cases where the probability distributions of one or more {system} parameters are not normal, uniform, or closely related {distributions}, due to the computational issues that arise when one wishes to define interpolation nodes for general distributions. This paper examines the use of the recently introduced weighted Leja nodes for that purpose. Weighted Leja interpolation rules are presented, along with a dimension-adaptive sparse interpolation algorithm, to be employed in the case of high-dimensional input uncertainty. The performance and reliability of the suggested approach is verified by four numerical experiments, where the respective models feature extreme value and truncated normal parameter distributions. Furthermore, the suggested approach is compared with a well-established polynomial chaos method and found to be either comparable or superior in terms of approximation and statistics estimation accuracy.

preprint2020arXiv

Magnetic Field Simulation with Data-Driven Material Modeling

This paper developes 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.

preprint2019arXiv

Robust Adaptive Least Squares Polynomial Chaos Expansions in High-Frequency Applications

We present an algorithm for computing sparse, least squares-based polynomial chaos expansions, incorporating both adaptive polynomial bases and sequential experimental designs. The algorithm is employed to approximate stochastic high-frequency electromagnetic models in a black-box way, in particular, given only a dataset of random parameter realizations and the corresponding observations regarding a quantity of interest, typically a scattering parameter. The construction of the polynomial basis is based on a greedy, adaptive, sensitivity-related method. The sequential expansion of the experimental design employs different optimality criteria, with respect to the algebraic form of the least squares problem. We investigate how different conditions affect the robustness of the derived surrogate models, that is, how much the approximation accuracy varies given different experimental designs. It is found that relatively optimistic criteria perform on average better than stricter ones, yielding superior approximation accuracies for equal dataset sizes. However, the results of strict criteria are significantly more robust, as reduced variations regarding the approximation accuracy are obtained, over a range of experimental designs. Two criteria are proposed for a good accuracy-robustness trade-off.

preprint2018arXiv

Coupled Simulation of Transient Heat Flow and Electric Currents in Thin Wires: Application to Bond Wires in Microelectronic Chip Packaging

This work addresses the simulation of heat flow and electric currents in thin wires. An important application is the use of bond wires in microelectronic chip packaging. The heat distribution is modeled by an electrothermal coupled problem, which poses numerical challenges due to the presence of different geometric scales. The necessity of very fine grids is relaxed by solving and embedding a 1D sub-problem along the wire into the surrounding 3D geometry. The arising singularities are described using de Rham currents. It is shown that the problem is related to fluid flow in porous 3D media with 1D fractures [C. D'Angelo, SIAM Journal on Numerical Analysis 50.1, pp. 194-215, 2012]. A careful formulation of the 1D-3D coupling condition is essential to obtain a stable scheme that yields a physical solution. Elliptic model problems are used to investigate the numerical errors and the corresponding convergence rates. Additionally, the transient electrothermal simulation of a simplified microelectronic chip package as used in industrial applications is presented.

preprint2016arXiv

Automatic calibration of damping layers in finite element time domain simulations

Matched layers are commonly used in numerical simulations of wave propagation to model (semi-)infinite domains. Attenuation functions describe the damping in layers, and provide a matching of the wave impedance at the interface between the domain of interest and the absorbing region. Selecting parameters in the attenuation functions is non-trivial. In this work, an optimisation procedure for automatically calibrating matched layers is presented. The procedure is based on solving optimisation problems constrained by partial differential equations with polynomial and piecewise-constant attenuation functions. We show experimentally that, for finite element time domain simulations, piecewise-constant attenuation function are at least as efficient as quadratic attenuation functions. This observation leads us to introduce consecutive matched layers as an alternative to perfectly matched layers, which can easily be employed for problems with arbitrary geometries. Moreover, the use of consecutive matched layers leads to a reduction in computational cost compared to perfectly matched layers. Examples are presented for acoustic, elastodynamic and electromagnetic problems. Numerical simulations are performed with the libraries FEniCS/DOLFIN and dolfin-adjoint, and the computer code to reproduce all numerical examples is made freely available.

preprint2016arXiv

Automatic Generation of Equivalent Electrothermal SPICE Netlists from 3D Electrothermal Field Models

Starting from a 3D electrothermal field problem discretized by the Finite Integration Technique, the equivalence to a circuit description is shown by exploiting the analogy to the Modified Nodal Analysis approach. Using this analogy, an algorithm for the automatic generation of a monolithic SPICE netlist is presented. Joule losses from the electrical circuit are included as heat sources in the thermal circuit. The thermal simulation yields nodal temperatures that influence the electrical conductivity. Apart from the used field discretization, this approach applies no further simplifications. An example 3D chip package is used to validate the algorithm.

preprint2016arXiv

Electrothermal Simulation of Bonding Wire Degradation under Uncertain Geometries

In this paper, electrothermal field phenomena in electronic components are considered. This coupling is tackled by multiphysical field simulations using the Finite Integration Technique (FIT). In particular, the design of bonding wires with respect to thermal degradation is investigated. Instead of resolving the wires by the computational grid, lumped element representations are introduced as point-to-point connections in the spatially distributed model. Fabrication tolerances lead to uncertainties of the wires' parameters and influence the operation and reliability of the final product. Based on geometric measurements, the resulting variability of the wire temperatures is determined using the stochastic electrothermal field-circuit model.

preprint2016arXiv

Isogeometric Simulation of Lorentz Detuning in Superconducting Accelerator Cavities

Cavities in linear accelerators suffer from eigenfrequency shifts due to mechanical deformation caused by the electromagnetic radiation pressure, a phenomenon known as Lorentz detuning. Estimating the frequency shift up to the needed accuracy by means of standard Finite Element Methods, is a complex task due to the non exact representation of the geometry and due to the necessity for mesh refinement when using low order basis functions. In this paper, we use Isogeometric Analysis for discretising both mechanical deformations and electromagnetic fields in a coupled multiphysics simulation approach. The combined high-order approximation of both leads to high accuracies at a substantially lower computational cost.