Source author record

Peter Arbenz

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

ResearcherUnclaimed source record

physics.acc-ph math-ph math.MP math.NA math.OC math.PR physics.comp-ph

Catalog footprint

What is connected

5works

7topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2019arXiv

Constrained multi-objective shape optimization of superconducting RF cavities considering robustness against geometric perturbations

High current storage rings, such as the Z-pole operating mode of the FCC-ee, require accelerating cavities that are optimized with respect to both the fundamental mode and the higher order modes. Furthermore, the cavity shape needs to be robust against geometric perturbations which could, for example, arise from manufacturing inaccuracies or harsh operating conditions at cryogenic temperatures. This leads to a constrained multi-objective shape optimization problem which is computationally expensive even for axisymmetric cavity shapes. In order to decrease the computation cost, a global sensitivity analysis is performed and its results are used to reduce the search space and redefine the objective functions. A massively parallel implementation of an evolutionary algorithm, combined with a fast axisymmetric Maxwell eigensolver and a frequency-tuning method is used to find an approximation of the Pareto front. The computed Pareto front approximation and a cavity shape with desired properties are shown. Further, the approach is generalized and applied to another type of cavity.

preprint2015arXiv

A Homotopy Method for Large-Scale Multi-Objective Optimization

A homotopy method for multi-objective optimization that produces uniformly sampled Pareto fronts by construction is presented. While the algorithm is general, of particular interest is application to simulation-based engineering optimization problems where economy of function evaluations, smoothness of result, and time-to-solution are critical. The presented algorithm achieves an order of magnitude improvement over other geometrically motivated methods, like Normal Boundary Intersection and Normal Constraint, with respect to solution evenness for similar computational expense. Furthermore, the resulting uniformity of solutions extends even to more difficult problems, such as those appearing in common Evolutionary Algorithm test cases.

preprint2014arXiv

Multilevel Monte Carlo for the Feynman-Kac Formula for the Laplace Equation

Since its formulation in the late 1940s, the Feynman-Kac formula has proven to be an effective tool for both theoretical reformulations and practical simulations of differential equations. The link it establishes between such equations and stochastic processes can be exploited to develop Monte Carlo sampling methods that are effective, especially in high dimensions. There exist many techniques of improving standard Monte Carlo sampling methods, a relatively new development being the so-called Multilevel Monte Carlo method. This paper investigates the applicability of multilevel ideas to the stochastic representation of partial differential equations by the Feynman- Kac formula, using the Walk on Sphere algorithm to generate the required random paths. We focus on the Laplace equation, the simplest elliptic PDE, while mentioning some extension possibilities.

preprint2012arXiv

A novel adaptive time stepping variant of the Boris-Buneman integrator for the simulation of particle accelerators with space charge

We show that adaptive time stepping in particle accelerator simulation is an enhancement for certain problems. The new algorithm has been implemented in the OPAL (Object Oriented Parallel Accelerator Library) framework, and is compared to the existing code. The idea is to adjust the frequency of costly self field calculations, which are needed to model Coulomb interaction (space charge) effects. In analogy to a Kepler orbit simulation that requires a higher time step resolution at the close encounter, we propose to choose the time step based on the magnitude of the space charge forces. Inspired by geometric integration techniques, our algorithm chooses the time step proportional to a function of the current phase space state instead of calculating a local error estimate like a conventional adaptive procedure. In this paper we build up on first observations made in recent work. A more profound argument is given on how exactly the time step should be chosen. An intermediate algorithm, initially built to allow a clearer analysis by introducing separate time steps for external field and self field integration, turned out to be useful in itself already for a large class of problems.

preprint2011arXiv

Perfectly Matched Layers in a Divergence Preserving ADI Scheme for Electromagnetics

For numerical simulations of highly relativistic and transversely accelerated charged particles including radiation fast algorithms are needed. While the radiation in particle accelerators has wavelengths in the order of 100 um the computational domain has dimensions roughly 5 orders of magnitude larger resulting in very large mesh sizes. The particles are confined to a small area of this domain only. To resolve the smallest scales close to the particles subgrids are envisioned. For reasons of stability the alternating direction implicit (ADI) scheme by D. N. Smithe et al. (J. Comput. Phys. 228 (2009) pp.7289-7299) for Maxwell equations has been adopted. At the boundary of the domain absorbing boundary conditions have to be employed to prevent reflection of the radiation. In this paper we show how the divergence preserving ADI scheme has to be formulated in perfectly matched layers (PML) and compare the performance in several scenarios.