Source author record

Antoine Laurain

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

ResearcherUnclaimed source record

math.OC math.AP math.NA Numerical Analysis

Catalog footprint

What is connected

6works

4topics

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.

preprint2026arXiv

FormOpt: A FEniCSx toolbox for level set-based shape optimization supporting parallel computing

This article presents the toolbox FormOpt for two- and three-dimensional shape optimization with parallel computing capabilities, built on the FEniCSx software framework. We introduce fundamental concepts of shape sensitivity analysis and their numerical applications, mainly for educational purposes, while also emphasizing computational efficiency via parallelism for practitioners. We adopt an optimize-then-discretize strategy based on the distributed shape derivative and its tensor representation, following the approach of \cite{MR3843884} and extending it in several directions. The numerical shape modeling relies on a level set method, whose evolution is driven by a descent direction computed from the shape derivative. Geometric constraints are treated accurately through a Proximal-Perturbed Lagrangian approach. FormOpt leverages the powerful features of FEniCSx, particularly its support for weak formulations of partial differential equations, diverse finite element types, and scalable parallelism. The implementation supports three different parallel computing modes: data parallelism, task parallelism, and a mixed mode. Data parallelism exploits FEniCSx's mesh partitioning features, and we implement a task parallelism mode which is useful for problems governed by a set of partial differential equations with varying parameters. The mixed mode conveniently combines both strategies to achieve efficient utilization of computational resources.

preprint2022arXiv

A shape optimization approach for electrical impedance tomography with point measurements

Working within the class of piecewise constant conductivities, the inverse problem of electrical impedance tomography can be recast as a shape optimization problem where the discontinuity interface is the unknown. Using Gröger's $W^{1}_p$-estimates for mixed boundary value problems, the averaged adjoint method is extended to the case of Banach spaces, which allows to compute the derivative of shape functionals involving point evaluations. We compute the corresponding distributed expression of the shape derivative and show that it may contain Dirac measures in addition to the usual domain integrals. We use this distributed shape derivative to devise a numerical algorithm, show various numerical results supporting the method, and based on these results we discuss the influence of the point measurements patterns on the quality of the reconstructions.

preprint2020arXiv

Shape optimization for superconductors governed by H(curl)-elliptic variational inequalities

This paper is devoted to the theoretical and numerical study of an optimal design problem in high-temperature superconductivity (HTS). The shape optimization problem is to find an optimal superconductor shape which minimizes a certain cost functional under a given target on the electric field over a specific domain of interest. For the governing PDE-model, we consider an elliptic curl-curl variational inequality (VI) of the second kind with an L1-type nonlinearity. In particular, the non-smooth VI character and the involved H(curl)-structure make the corresponding shape sensitivity analysis challenging. To tackle the non-smoothness, a penalized dual VI formulation is proposed, leading to the G{â}teaux differentiability of the corresponding dual variable mapping. This property allows us to derive the distributed shape derivative of the cost functional through rigorous shape calculus on the basis of the averaged adjoint method. The developed shape derivative turns out to be uniformly stable with respect to the penalization parameter, and strong convergence of the penalized problem is guaranteed. Based on the achieved theoretical findings, we propose 3D numerical solutions, realized using a level set algorithm and a Newton method with the Nedelec edge element discretization. Numerical results indicate a favourable and efficient performance of the proposed approach for a specific HTS application in superconducting shielding.

preprint2016arXiv

Properties of optimizers of the principal eigenvalue with indefinite weight and Robin conditions

In this paper, we are interested in the analysis of a well-known free boundary/shape optimization problem motivated by some issues arising in population dynamics. The question is to determine optimal spatial arrangements of favorable and unfavorable regions for a species to survive. The mathematical formulation of the model leads to an indefinite weight linear eigenvalueproblem in a fixed box $Ω$ and we consider the general case of Robin boundary conditions on $\partialΩ$. It is well known that it suffices to consider {\it bang-bang} weights taking two values of different signs, that can be parametrized by the characteristic function of the subset $E$ of $Ω$ on which resources are located. Therefore, the optimal spatial arrangement is obtained by minimizing the positive principal eigenvalue with respect to $E$, under a volume constraint. By using symmetrization techniques, as well as necessary optimality conditions, we prove new qualitative results on the solutions. Namely, we completely solve the problem in dimension 1, we prove the counter-intuitive result that the ball is almost never a solution in dimension 2 or higher, despite what suggest the numerical simulations. We also introduce a new rearrangement in the ball allowing to get a better candidate than the ball for optimality when Neumann boundary conditions are imposed. We also provide numerical illustrations of our results and of the optimal configurations.

preprint2015arXiv

Distributed shape derivative via averaged adjoint method and applications

The structure theorem of Hadamard-Zolésio states that the derivative of a shape functional is a distribution on the boundary of the domain depending only on the normal perturbations of a smooth enough boundary. Actually the domain representation, also known as distributed shape derivative, is more general than the boundary expression as it is well-defined for shapes having a lower regularity. It is customary in the shape optimization literature to assume regularity of the domains and use the boundary expression of the shape derivative for numerical algorithms. In this paper we describe several advantages of the distributed shape derivative in terms of generality, easiness of computation and numerical implementation. We identify a tensor representation of the distributed shape derivative, study its properties and show how it allows to recover the boundary expression directly. We use a novel Lagrangian approach, which is applicable to a large class of shape optimization problems, to compute the distributed shape derivative. We also apply the technique to retrieve the distributed shape derivative for electrical impedance tomography. Finally we explain how to adapt the level set method to the distributed shape derivative framework and present numerical results.

preprint2010arXiv

On a Bernoulli problem with geometric constraints

A Bernoulli free boundary problem with geometrical constraints is studied. The domain $\Om$ is constrained to lie in the half space determined by $x_1\geq 0$ and its boundary to contain a segment of the hyperplane $\{x_1=0\}$ where non-homogeneous Dirichlet conditions are imposed. We are then looking for the solution of a partial differential equation satisfying a Dirichlet and a Neumann boundary condition simultaneously on the free boundary. The existence and uniqueness of a solution have already been addressed and this paper is devoted first to the study of geometric and asymptotic properties of the solution and then to the numerical treatment of the problem using a shape optimization formulation. The major difficulty and originality of this paper lies in the treatment of the geometric constraints.