Toplogical derivative for nonlinear magnetostatic problem

The topological derivative represents the sensitivity of a domain-dependent functional with respect to a local perturbation of the domain and is a valuable tool in topology optimization. Motivated by an application from electrical engineering, we derive the topological derivative for an optimization problem which is constrained by the quasilinear equation of two-dimensional magnetostatics. Here, the main ingredient is to establish a sufficiently fast decay of the variation of the direct state at scale 1 as $|x|\rightarrow \infty$. In order to apply the method in a bi-directional topology optimization algorithm, we derive both the sensitivity for introducing air inside ferromagnetic material and the sensitivity for introducing material inside an air region. We explicitly compute the arising polarization matrices and introduce a way to efficiently evaluate the obtained formulas. Finally, we employ the derived formulas in a level-set based topology optimization algorithm and apply it to the design optimization of an electric motor.Comment: 54 pages, 9 figure

Isogeometric Simulation and Shape Optimization with Applications to Electrical Machines

Future e-mobility calls for efficient electrical machines. For different areas of operation, these machines have to satisfy certain desired properties that often depend on their design. Here we investigate the use of multipatch Isogeometric Analysis (IgA) for the simulation and shape optimization of the electrical machines. In order to get fast simulation and optimization results, we use non-overlapping domain decomposition (DD) methods to solve the large systems of algebraic equations arising from the IgA discretization of underlying partial differential equations. The DD is naturally related to the multipatch representation of the computational domain, and provides the framework for the parallelization of the DD solvers

Asymptotic analysis and topological derivative for 3D quasi-linear magnetostatics

In this paper we study the asymptotic behaviour of the quasilinear $curl$-$curl$ equation of 3D magnetostatics with respect to a singular perturbation of the differential operator and prove the existence of the topological derivative using a Lagrangian approach. We follow the strategy proposed in our recent previous work (arXiv:1907.13420) where a systematic and concise way for the derivation of topological derivatives for quasi-linear elliptic problems in $H^1$ is introduced. In order to prove the asymptotics for the state equation we make use of an appropriate Helmholtz decomposition. The evaluation of the topological derivative at any spatial point requires the solution of a nonlinear transmission problem. We discuss an efficient way for the numerical evaluation of the topological derivative in the whole design domain using precomputation in an offline stage. This allows us to use the topological derivative for the design optimization of an electrical machine.Comment: 26 pages, 5 figure

A unified approach to shape and topological sensitivity analysis of discretized optimal design problems

We introduce a unified sensitivity concept for shape and topological perturbations and perform the sensitivity analysis for a discretized PDE-constrained design optimization problem in two space dimensions. We assume that the design is represented by a piecewise linear and globally continuous level set function on a fixed finite element mesh and relate perturbations of the level set function to perturbations of the shape or topology of the corresponding design. We illustrate the sensitivity analysis for a problem that is constrained by a reaction-diffusion equation and draw connections between our discrete sensitivities and the well-established continuous concepts of shape and topological derivatives. Finally, we verify our sensitivities and illustrate their application in a level-set-based design optimization algorithm where no distinction between shape and topological updates has to be made

Regularization and finite element error estimates for elliptic distributed optimal control problems with energy regularization and state or control constraints

In this paper we discuss the numerical solution of elliptic distributed optimal control problems with state or control constraints when the control is considered in the energy norm. As in the unconstrained case we can relate the regularization parameter and the finite element mesh size in order to ensure an optimal order of convergence which only depends on the regularity of the given target, also including discontinuous target functions. While in most cases, state or control constraints are discussed for the more common $L^2$ regularization, much less is known in the case of energy regularizations. But in this case, and for both control and state constraints, we can formulate first kind variational inequalities to determine the unknown state, from wich we can compute the control in a post processing step. Related variational inequalities also appear in obstacle problems, and are well established both from a mathematical and a numerical analysis point of view. Numerical results confirm the applicability and accuracy of the proposed approach

Shape Optimization of Rotating Electric Machines using Isogeometric Analysis and Harmonic Stator-Rotor Coupling

This work deals with shape optimization of electric machines using isogeometric analysis. Isogeometric analysis is particularly well suited for shape optimization as it allows to easily modify the geometry without remeshing the domain. A 6-pole permanent magnet synchronous machine is modeled using a multipatch isogeometric approach and rotation of the machine is realized by modeling the stator and rotor domain separately and coupling them at the interface using harmonic basis functions. Shape optimization is applied to the model minimizing the total harmonic distortion of the electromotive force as a goal functional

Semiclassical theory of cavity-assisted atom cooling

We present a systematic semiclassical model for the simulation of the dynamics of a single two-level atom strongly coupled to a driven high-finesse optical cavity. From the Fokker-Planck equation of the combined atom-field Wigner function we derive stochastic differential equations for the atomic motion and the cavity field. The corresponding noise sources exhibit strong correlations between the atomic momentum fluctuations and the noise in the phase quadrature of the cavity field. The model provides an effective tool to investigate localisation effects as well as cooling and trapping times. In addition, we can continuously study the transition from a few photon quantum field to the classical limit of a large coherent field amplitude.Comment: 10 pages, 8 figure

Large Observatory for X-ray Timing (LOFT-P): A Probe-Class Mission Concept Study

LOFT-P is a mission concept for a NASA Astrophysics Probe-Class (less than $1B) X-ray timing mission, based on the LOFT M-class concept originally proposed to ESA's M3 and M4 calls. LOFT-P requires very large collecting area, high time resolution, good spectral resolution, broadband spectral coverage (2-30 keV), highly flexible scheduling, and an ability to detect and respond promptly to time-critical targets of opportunity. Many of LOFTP's targets are bright, rapidly varying sources, so these measurements are synergistic to imaging and high-resolution spectroscopy instruments, addressing much smaller distance scales than are possible without very long baseline X-ray interferometry, and using complementary techniques to address the geometry and dynamics of emission regions. LOFT-P was presented as an example mission to the head of NASA's Astrophysics Division, to demonstrate the strong community support for creation of a probe-class, for missions costing between$500M and $1B. We submitted a white paper4 in response to NASA PhysPAG's call for white papers: Probe-class Mission Concepts, describing LOFT-P science and a simple extrapolation from the ESA study costs. The next step for probe-class missions will be input into the NASA Astrophysics Decadal Survey to encourage the creation of a probe-class opportunity. We report on a 2016 study by MSFC's Advanced Concepts Office of LOFT-P, a US-led probe-class LOFT concept