Related papers: Computing the shape gradient of stellarator coil c…
This thesis deals with shape optimization for contact mechanics. More specifically, the linear elasticity model is considered under the small deformations hypothesis, and the elastic body is assumed to be in contact (sliding or with Tresca…
The Physalis method is suitable for the simulation of flows with suspended spherical particles. It differs from standard immersed boundary methods due to the use of a local spectral representation of the solution in the neighborhood of each…
In this work, we utilize new coil objectives for stellarator optimization with autodifferentiation, including pointwise and net coil-coil forces and torques. We use these methods to perform the first large-scale optimization of planar…
One common approach to computing the magnetic field produced by a filamentary current-carrying coil is to approximate the coil as a series of straight segments. The Biot-Savart field from each straight segment is known analytically.…
The differentiable programming paradigm is a cornerstone of modern scientific computing. It refers to numerical methods for computing the gradient of a numerical model's output. Many scientific models are based on differential equations,…
A novel adaptive technique for electromagnetic Particle In Cell (PIC) plasma simulations is presented here. Two main issues are identified in designing adaptive techniques for PIC simulation: first, the choice of the size of the particle…
A finite difference scheme is used to develop a numerical method to solve the flow of an unbounded viscoelastic fluid with zero to moderate inertia around a prolate spheroidal particle. The equations are written in prolate spheroidal…
Shape optimization based on surface gradients and the Hadarmard-form is considered for a compressible viscous fluid. Special attention is given to the difference between the 'function composition' approach involving local shape derivatives…
In transient simulations of particulate Stokes flow, to accurately capture the interaction between the constituent particles and the confining wall, the discretization of the wall often needs to be locally refined in the region approached…
This paper presents a general framework for calculating the dimension of spline spaces over arbitrary rectilinear partitions using the smoothing cofactor method. The approach extends existing dimension theory for polynomial splines over…
This paper investigates, without any regularization or penalization procedure, a shape optimization problem involving a simplified friction phenomena modeled by a scalar Tresca friction law. Precisely, using tools from convex and…
We consider the inverse obstacle scattering problem of determining both the shape and the "equivalent impedance" from far field measurements at a fixed frequency. In this work, the surface impedance is represented by a second order surface…
Outside the core of the plasma, the plasma current and pressure rapidly transition to zero in a scrape-off or edge region or plasma-vacuum interface. However, existing tools for fixed-boundary magnetohydrodynamic equilibria in 2D and 3D…
This paper describes a method for solving smooth nonconvex minimization problems subject to bound constraints with good worst-case complexity guarantees and practical performance. The method contains elements of two existing methods: the…
The problem of finding of analytical gradients (derivatives over atoms coordinates) of solvation energies can be decomposed on two subtasks: at the first stage we search for parameters of the superficial devices (three coordinates, three…
We consider a problem relating to magnetic confinement devices known as stellarators. Plasma is confined by magnetic fields generated by current-carrying coils, and here we investigate how closely to the plasma they need to be positioned.…
Confining energetic ions such as alpha particles is a prime concern in the design of stellarators. However, directly measuring alpha confinement through numerical simulation of guiding-center trajectories has been considered to be too…
In this article, we consider the problem of optimal design of a compliant structure under a volume constraint, within the framework of linear elasticity. We introduce the pure displacement and the dual mixed formulations of the linear…
We develop a new optimisation technique that combines multiresolution subdivision surfaces for boundary description with immersed finite elements for the discretisation of the primal and adjoint problems of optimisation. Similar to wavelets…
Stellarators are magnetic confinement devices under active development to deliver steady-state carbon-free fusion energy. Their design involves a high-dimensional, constrained optimization problem that requires expensive physics simulations…