Researcher profile

Sebastian Schöps

Sebastian Schöps contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
31works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

31 published item(s)

preprint2026arXiv

Energy-based modeling for field-circuit coupling

This paper presents a generalized energy-based modeling framework extending recent formulations tailored for differential-algebraic equations. The proposed structure, inspired by the port-Hamiltonian formalism, ensures passivity, preserves the power balance, and facilitates the consistent interconnection of subsystems. A particular focus is put on low-frequency power applications in electrical engineering. Stranded, solid, and foil conductor models are investigated in the context of the eddy current problem. Each conductor model is shown to fit into the generalized energy-based structure, which allows their structure-preserving coupling with electrical circuits described by modified nodal analysis. Theoretical developments are validated through a numerical simulation of an oscillator circuit, demonstrating energy conservation in lossless scenarios and controlled dissipation when eddy currents are present. The applicability of the methodology towards engineering applications is studied through a numerical simulation of a nonlinear three-phase transformer.

preprint2026arXiv

Learning electromagnetic fields based on finite element basis functions

Parametric surrogate models of electric machines are widely used for efficient design optimization and operational monitoring. Addressing geometry variations, spline-based computer-aided design representations play a pivotal role. In this study, we propose a novel approach that combines isogeometric analysis, proper orthogonal decomposition and deep learning to enable rapid and physically consistent predictions by directly learning spline basis coefficients. The effectiveness of this method is demonstrated using a parametric nonlinear magnetostatic model of a permanent magnet synchronous machine.

preprint2026arXiv

Multi-Level Monte Carlo sampling with Parallel-in-Time Integration for Uncertainty Quantification in Electric Machine Simulation

While generally considered computationally expensive, Uncertainty Quantification using Monte Carlo sampling remains beneficial for applications with uncertainties of high dimension. As an extension of the naive Monte Carlo method, the Multi-Level Monte Carlo method reduces the overall computational effort, but is unable to reduce the time to solution in a sufficiently parallel computing environment. In this work, we propose a Uncertainty Quantification method combining Multi-Level Monte Carlo sampling and Parallel-in-Time integration for select samples, exploiting remaining parallel computing capacity to accelerate the computation. While effective at reducing the time-to-solution, Parallel-in-Time integration methods greatly increase the total computational effort. We investigate the tradeoff between time-to-solution and total computational effort of the combined method, starting from theoretical considerations and comparing our findings to two numerical examples. There, a speedup of 12 - 45% compared to Multi-Level Monte Carlo sampling is observed, with an increase of 15 - 18% in computational effort.

preprint2025arXiv

Foil Conductor Model for Efficient Simulation of HTS Coils in Large Scale Applications

Homogenization techniques are an appealing approach to reduce computational complexity in systems containing coils with large numbers of high temperature superconductor (HTS) tapes. Resolving all the coated conductor layers and turns in coils is often computationally prohibitive. In this paper, we extend the foil conductor model, well-known in normal conducting applications, to applications with insulated HTS coils. To enhance the numerical performance of the model, the conventional formulation based on A-V is extended to J-A-V. The model is verified to be suitable for simulations of superconductors and to accelerate the calculations compared to resolving all the individual layers. The performance of both the A-V and J-A-V formulated models is examined, and the J-A-V variant is concluded to be advantageous.

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

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

Freeform shape optimization of a compact DC photo-electron gun using isogeometric analysis

Compact DC high-voltage photo-electron guns are able to meet the sophisticated demands of high-current applications such as energy recovery linacs. A main design parameter for such sources is the electric field strength, which depends on the electrode geometry and is limited by the field emission threshold of the electrode material. In order to minimize the maximum field strength for optimal gun operation, isogeometric analysis (IGA) can be used to exploit the axisymmetric geometry and describe its cross section by non-uniform rational B-splines, the control points of which are the parameters to be optimized. This computationally efficient method is capable of describing CAD-generated geometries using open source software (GeoPDEs, NLopt, Octave) and it can simplify the step from design to simulation. We will present the mathematical formulation, the software workflow, and the results of an IGA-based shape optimization for a planned high-voltage upgrade of the DC photogun teststand Photo-CATCH at TU Darmstadt. The software builds on a general framework for isogeometric analysis and allows for easy adaptations to other geometries or quantities of interest. Simulations assuming a bias voltage of -300 kV yielded maximum field gradients of 9.06 MV/m on the surface of an inverted insulator electrode and below 3 MV/m on the surface of the photocathode.

preprint2022arXiv

Hermite-type modifications of BOBYQA for optimization with some partial derivatives

In this work we propose two Hermite-type optimization methods, Hermite least squares and Hermite BOBYQA, specialized for the case that some partial derivatives of the objective function are available and others are not. The main objective is to reduce the number of objective function calls by maintaining the convergence properties. Both methods are modifications of Powell's derivative-free BOBYQA algorithm. But instead of (underdetermined) interpolation for building the quadratic subproblem in each iteration, the training data is enriched with first and -- if possible -- second order derivatives and then (weighted) least squares regression is used. Proofs for global convergence are discussed and numerical results are presented. Further, the applicability is verified for a realistic test case in the context of yield optimization. Numerical tests show that the Hermite least squares approach outperforms classic BOBYQA if half or more partial derivatives are available. In addition, the Hermite-type approaches achieve more robustness and thus better performance in case of noisy objective functions.

preprint2021arXiv

Deep Learning-based Prediction of Key Performance Indicators for Electrical Machine

The design of an electrical machine can be quantified and evaluated by Key Performance Indicators (KPIs) such as maximum torque, critical field strength, costs of active parts, sound power, etc. Generally, cross-domain tool-chains are used to optimize all the KPIs from different domains (multi-objective optimization) by varying the given input parameters in the largest possible design space. This optimization process involves magneto-static finite element simulation to obtain these decisive KPIs. It makes the whole process a vehemently time-consuming computational task that counts on the availability of resources with the involvement of high computational cost. In this paper, a data-aided, deep learning-based meta-model is employed to predict the KPIs of an electrical machine quickly and with high accuracy to accelerate the full optimization process and reduce its computational costs. The focus is on analyzing various forms of input data that serve as a geometry representation of the machine. Namely, these are the cross-section image of the electrical machine that allows a very general description of the geometry relating to different topologies and the classical way with scalar parametrization of geometry. The impact of the resolution of the image is studied in detail. The results show a high prediction accuracy and proof that the validity of a deep learning-based meta-model to minimize the optimization time. The results also indicate that the prediction quality of an image-based approach can be made comparable to the classical way based on scalar parameters.

preprint2021arXiv

Numerical Analysis of the Screening Current-Induced Magnetic Field in the HTS Insert Dipole Magnet Feather-M2.1-2

Screening currents are field-induced dynamic phenomena which occur in superconducting materials, leading to persistent magnetization. Such currents are of importance in ReBCO tapes, where the large size of the superconducting filaments gives rise to strong magnetization phenomena. In consequence, superconducting accelerator magnets based on ReBCO tapes might experience a relevant degradation of the magnetic field quality in the magnet aperture, eventually leading to particle beam instabilities. Thus, persistent magnetization phenomena need to be accurately evaluated. In this paper, the 2D finite element model of the Feather-M2.1-2 magnet is presented. The model is used to analyze the influence of the screening current-induced magnetic field on the field quality in the magnet aperture. The model relies on a coupled field formulation for eddy current problems in time-domain. The formulation is introduced and verified against theoretical references. Then, the numerical model of the Feather-M2.1-2 magnet is detailed, highlighting the key assumptions and simplifications. The numerical results are discussed and validated with available magnetic measurements. A satisfactory agreement is found, showing the capability of the numerical tool in providing accurate analysis of the dynamic behavior of the Feather-M2.1-2 magnet.

preprint2021arXiv

On torque computation in electric machine simulation by harmonic mortar methods

The use of trigonometric polynomials as Lagrange multipliers in the harmonic mortar method enables an efficient and elegant treatment of relative motion in the stator-rotor coupling of electric machine simulation. Explicit formulas for the torque computation are derived by energetic considerations, and their realization by harmonic mortar finite element and isogemetric analysis discretizations is discussed. Numerical tests are presented to illustrate the theoretical results and demonstrate the potential of harmonic mortar methods for the evaluation of torque ripples.

preprint2021arXiv

Parallel-in-Time Solution of Eddy Current Problems Using Implicit and Explicit Time-stepping Methods

The time domain analysis of eddy current problems often requires the simulation of long time intervals, e.g. until a steady state is reached. Fast-switching excitations e.g. in pulsedwidth modulated signals require in addition very small time step sizes that significantly increase computation time. To speed up the simulation, parallel-in-time methods can be used. In this paper, we investigate the combination of explicit and implicit time integration methods in the context of the parallel-in-time method Parareal and using a simplified model for the coarse propagator.

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

A Darwin Time Domain Scheme for the Simulation of Transient Quasistatic Electromagnetic Fields Including Resistive, Capacitive and Inductive Effects

The Darwin field model addresses an approximation to Maxwell's equations where radiation effects are neglected. It allows to describe general quasistatic electromagnetic field phenomena including inductive, resistive and capacitive effects. A Darwin formulation based on the Darwin-Ampère equation and the implicitly included Darwin-continuity equation yields a non-symmetric and ill-conditioned algebraic systems of equations received from applying a geometric spatial discretization scheme and the implicit backward differentiation time integration method. A two-step solution scheme is presented where the underlying block-Gauss-Seidel method is shown to change the initially chosen gauge condition and the resulting scheme only requires to solve a weakly coupled electro-quasistatic and a magneto-quasistatic discrete field formulation consecutively in each time step. Results of numerical test problems validate the chosen approach.

preprint2020arXiv

A multilevel Monte Carlo method for high-dimensional uncertainty quantification of low-frequency electromagnetic devices

This work addresses uncertainty quantification of electromagnetic devices determined by the eddy current problem. The multilevel Monte Carlo (MLMC) method is used for the treatment of uncertain parameters while the devices are discretized in space by the finite element method. Both methods yield numerical approximations such that the total errors is split into stochastic and spatial contributions. We propose a particular implementation where the spatial error is controlled based on a Richardson-extrapolation-based error indicator. The stochastic error in turn is efficiently reduced in the MLMC approach by distributing the samples on multiple grids. The method is applied to a toy problem with closed-form solution and a permanent magnet synchronous machine with uncertainties. The uncertainties under consideration are related to the material properties in the stator and the magnets in the rotor. The examples show that the error indicator works reliably, the meshes used for the different levels do not have to be nested and, most importantly, MLMC reduces the computational cost by at least one order of magnitude compared to standard Monte Carlo.

preprint2020arXiv

Efficient Parallel-in-Time Solution of Time-Periodic Problems Using a Multi-Harmonic Coarse Grid Correction

This paper presents a highly-parallelizable parallel-in-time algorithm for efficient solution of nonlinear time-periodic problems. It is based on the time-periodic extension of the Parareal method, known to accelerate sequential computations via parallelization on the fine grid. The proposed approach reduces the complexity of the periodic Parareal solution by introducing a simplified Newton algorithm, which allows an additional parallelization on the coarse grid. In particular, at each Newton iteration a multi-harmonic correction is performed, which converts the block-cyclic periodic system in the time domain into a block-diagonal system in the frequency domain, thereby solving for each frequency component in parallel. The convergence analysis of the method is discussed for a one-dimensional model problem. The introduced algorithm and several existing solution approaches are compared via their application to the eddy current problem for both linear and nonlinear models of a coaxial cable. Performance of the considered methods is also illustrated for a three-dimensional transformer model.

preprint2020arXiv

Efficient Simulation of Field/Circuit Coupled Systems with Parallelised Waveform Relaxation

This paper proposes an efficient parallelised computation of field/circuit coupled systems co-simulated with the Waveform Relaxation (WR) technique. The main idea of the introduced approach lies in application of the parallel-in-time method parareal to the WR framework. Acceleration obtained by the time-parallelisation is further increased in the context of micro-macro parareal. Here, the field system is replaced by a lumped model in the circuit environment for the sequential computations of parareal. The introduced algorithm is tested with a model of a single-phase isolation transformer coupled to a rectifier circuit.

preprint2020arXiv

Enhanced adaptive surrogate models with applications in uncertainty quantification for nanoplasmonics

We propose an efficient surrogate modeling technique for uncertainty quantification. The method is based on a well-known 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 non-trivial 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.

preprint2020arXiv

On the stability of harmonic mortar methods with application to electric machines

Harmonic stator-rotor coupling offers a promising approach for the interconnection of rotating subsystems in the simulation of electric machines. This paper studies the stability of discretization schemes based on harmonic coupling in the framework of mortar methods for Poisson-like problems. A general criterion is derived that allows to ensure the relevant inf-sup stability condition for a variety of specific discretization approaches, including finite-element methods and isogeometric analysis with harmonic mortar coupling. The validity and sharpness of the theoretical results is demonstrated by numerical tests.

preprint2020arXiv

Parallel-in-Time Simulation of an Electrical Machine using MGRIT

We apply the multigrid-reduction-in-time (MGRIT) algorithm to an eddy current simulation of a two-dimensional induction machine supplied by a pulse-width-modulation signal. To resolve the fast-switching excitations, small time steps are needed, such that parallelization in time becomes highly relevant for reducing the simulation time. The MGRIT algorithm is well suited for introducing time parallelism in the simulation of electrical machines using existing application codes, as MGRIT is a non-intrusive approach that essentially uses the same time integrator as a traditional time-stepping algorithm. We investigate effects of spatial coarsening on MGRIT convergence when applied to two numerical models of an induction machine, one with linear material laws and a full nonlinear model. Parallel results demonstrate significant speedup in the simulation time compared to sequential time stepping, even for moderate numbers of processors.

preprint2020arXiv

Parallel-in-Time Simulation of Power Converters Using Multirate PDEs

This paper presents an efficient numerical algorithm for the simulation of pulse-width modulated power converters via parallelization in time domain. The method applies the multirate partial differential equation approach on the coarse grid of the (two-grid) parallel-in-time algorithm Parareal. Performance of the proposed approach is illustrated via its application to a DC-DC converter.

preprint2020arXiv

Towards a Parallel-in-Time Calculation of Time-Periodic Solutions with Unknown Period

This paper presents a novel parallel-in-time algorithm able to compute time-periodic solutions of problems where the period is not given. Exploiting the idea of the multiple shooting method, the proposed approach calculates the initial values at each subinterval as well as the corresponding period iteratively. As in the Parareal method, parallelization in the time domain is performed using discretization on a two-level grid. A special linearization of the time-periodic system on the coarse grid is introduced to speed up the computations. The iterative algorithm is verified via its application to the Colpitt oscillator model.

preprint2020arXiv

Waveform relaxation for low frequency coupled field/circuit differential-algebraic models of index 2

Motivated by the task to design quench protection systems for superconducting magnets in particle accelerators we address a coupled field/circuit simulation based on a magneto-quasistatic field modeling. We investigate how a waveform relaxation of Gauß-Seidel type performs for a coupled simulation when circuit solving packages are used that describe the circuit by the modified nodal analysis. We present sufficient convergence criteria for the coupled simulation of FEM discretised field models and circuit models formed by a differential-algebraic equation (DAE) system of index 2. In particular, we demonstrate by a simple benchmark system the drastic influence of the circuit topology on the convergence behavior of the coupled simulation.

preprint2020arXiv

Yield Optimization using Hybrid Gaussian Process Regression and a Genetic Multi-Objective Approach

Quantification and minimization of uncertainty is an important task in the design of electromagnetic devices, which comes with high computational effort. We propose a hybrid approach combining the reliability and accuracy of a Monte Carlo analysis with the efficiency of a surrogate model based on Gaussian Process Regression. We present two optimization approaches. An adaptive Newton-MC to reduce the impact of uncertainty and a genetic multi-objective approach to optimize performance and robustness at the same time. For a dielectrical waveguide, used as a benchmark problem, the proposed methods outperform classic approaches.

preprint2019arXiv

3D Field Simulation Model for Bond Wire On-Chip Inductors Validated by Measurements

This paper proposes 3D field simulation models for different designs of integrated bond wire on-chip inductors. To validate the simulation models, prototypes for three designs with air and ferrite cores are manufactured and measured. For air core inductors, high agreement between simulation and measurement is obtained. For ferrite core inductors, accurate models require an exact characterization of the ferrite material. These models enable the prediction of magnetic field influences on underlying integrated circuits.

preprint2019arXiv

Accelerated Steady-State Torque Computation for Induction Machines using Parallel-In-Time Algorithms

This paper focuses on efficient steady-state computations of induction machines. In particular, the periodic Parareal algorithm with initial-value coarse problem (PP-IC) is considered for acceleration of classical time-stepping simulations via non-intrusive parallelization in time domain, i.e., existing implementations can be reused. Superiority of this parallel-in-time method is in its direct applicability to time-periodic problems, compared to, e.g, the standard Parareal method, which only solves an initial-value problem, starting from a prescribed initial value. PP-IC is exploited here to obtain the steady state of several operating points of an induction motor, developed by Robert Bosch GmbH. Numerical experiments show that acceleration up to several dozens of times can be obtained, depending on availability of parallel processing units. Comparison of PP-IC with existing time-periodic explicit error correction method highlights better robustness and efficiency of the considered time-parallel approach.

preprint2019arXiv

IsoGeometric Approximations for Electromagnetic Problems in Axisymmetric Domains

We propose a numerical method for the solution of electromagnetic problems on axisymmetric domains, based on a combination of a spectral Fourier approximation in the azimuthal direction with an IsoGeometric Analysis (IGA) approach in the radial and axial directions. This combination allows to blend the flexibility and accuracy of IGA approaches with the advantages of a Fourier representation on axisymmetric domains. It also allows to reduce significantly the computational cost by decoupling of the computations required for each Fourier mode. We prove that the discrete approximation spaces employed functional space constitute a closed and exact de Rham sequence. Numerical simulations of relevant benchmarks confirm the high order convergence and other computational advantages of the proposed method.

preprint2019arXiv

Uncertainty Modeling and Analysis of the European X-ray Free Electron Laser Cavities Manufacturing Process

This paper reports on comprehensive efforts on uncertainty quantification and global sensitivity analysis for accelerator cavity design. As a case study object the TESLA shaped superconducting cavities, as produced for the European X-ray Free Electron Laser (EXFEL), are selected. The choice for these cavities is explained by the available measurement data that can be leveraged to substantiate the simulation model. Each step of the manufacturing chain is documented together with the involved uncertainties. Several of these steps are mimicked on the simulation side, e.g. by introducing a random eigenvalue problem. The uncertainties are then quantified numerically and in particular the sensitivities give valuable insight into the systems behavior. We also compare these findings to purely statistical studies carried out for the manufactured cavities. More advanced, adaptive, surrogate modeling techniques are adopted, which are crucial to incorporate a large number of uncertain parameters. The main contribution is the detailed comparison and fusion of measurement results for the EXFEL cavities on the one hand and simulation based uncertainty studies on the other hand. After introducing the quantities of physical interest for accelerator cavities and the Maxwell eigenvalue problem, the details on the manufacturing of the EXFEL cavities and measurements are reported. This is followed by uncertainty modeling with quantification studies.

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.

preprint2017arXiv

A 2-D Finite-Element Model for Electro-Thermal Transients in Accelerator Magnets

Superconducting accelerator magnets require sophisticated monitoring and means of protection due to the large energy stored in the magnetic field. Numerical simulations play a crucial role in understanding transient phenomena occurring within the magnet, and can, therefore, help to prevent disruptive consequences. We present a 2-D FEM model for the simulation of electro-thermal transients occurring in superconducting accelerator magnets. The magnetoquasistatic problem is solved with a modified magnetic vector potential formulation, where the cable eddy currents are resolved in terms of their equivalent magnetization. The heat balance equation is then investigated, and the relevant heat sources are discussed. The model implements a two-port component interface and is resolved, as part of an electrical circuit, in a cooperative simulation scheme with a lumped-parameter network.