Researcher profile

Ulrich Römer

Ulrich Römer contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
10topics
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

9 published item(s)

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.

preprint2021arXiv

Bayesian Dynamical System Identification With Unified Sparsity Priors And Model Uncertainty

This work is concerned with uncertainty quantification in reduced-order dynamical system identification. Reduced-order models for system dynamics are ubiquitous in design and control applications and recent efforts focus on their data-driven construction. Our starting point is the sparse-identification of nonlinear dynamics (SINDy) framework, which reformulates system identification as a regression problem, where unknown functions are approximated from a sparse subset of an underlying library. In this manuscript, we formulate this system identification method in a Bayesian framework to handle parameter and structural model uncertainties. We present a general approach to enforce sparsity, which builds on the recently introduced class of neuronized priors. We perform comparisons between different variants such as Lasso, horseshoe, and spike and slab priors, which are all obtained by modifying a single activation function. We also outline how state observation noise can be incorporated with a probabilistic state-space model. The resulting Bayesian regression framework is robust and simple to implement. We apply the method to two generic numerical applications, the pendulum and the Lorenz system, and one aerodynamic application using experimental measurements.

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

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

Goal-oriented a posteriori estimation of numerical errors in the solution of multiphysics systems

This paper develops a general methodology for a posteriori error estimation in time-dependent multiphysics numerical simulations. The methodology builds upon the generalized-structure additive Runge--Kutta (GARK) approach to time integration. GARK provides a unified formulation of multimethods that simulate complex systems by applying different discretization formulas and/or different time steps to individual components of the system. We derive discrete GARK adjoints and analyze their time accuracy. Based on the adjoint method, we establish computable a posteriori identities for the impacts of both temporal and spatial discretization errors on a given goal function. Numerical examples with reaction-diffusion systems illustrate the accuracy of the derived error measures. Local error decompositions are used to illustrate the power of this framework in adaptive refinements of both temporal and spatial meshes.

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

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

Solving Parameter Estimation Problems with Discrete Adjoint Exponential Integrators

The solution of inverse problems in a variational setting finds best estimates of the model parameters by minimizing a cost function that penalizes the mismatch between model outputs and observations. The gradients required by the numerical optimization process are computed using adjoint models. Exponential integrators are a promising family of time discretizations for evolutionary partial differential equations. In order to allow the use of these discretizations in the context of inverse problems adjoints of exponential integrators are required. This work derives the discrete adjoint formulae for a W-type exponential propagation iterative methods of Runge-Kutta type (EPIRK-W). These methods allow arbitrary approximations of the Jacobian while maintaining the overall accuracy of the forward integration. The use of Jacobian approximation matrices that do not depend on the model state avoids the complex calculation of Hessians in the discrete adjoint formulae, and allows efficient adjoint code generation via algorithmic differentiation. We use the discrete EPIRK-W adjoints to solve inverse problems with the Lorenz-96 model and a computational magnetics benchmark test. Numerical results validate our theoretical derivations.