Magnetic forces produced by rectangular permanent magnets in static microsystems

Finite element numerical simulations were carried out in 2D geometries to map the magnetic field and force distribution produced by rectangular permanent magnets as a function of their size and position with respect to a microchannel. A single magnet, two magnets placed in attraction and in repulsion have been considered. The goal of this work is to show where magnetic beads are preferentially captured in a microchannel. These simulations were qualitatively corroborated, in one geometrical case, by microscopic visualizations of magnetic bead plug formation in a capillary. The results show that the number of plugs is configuration dependent with: in attraction, one plug in the middle of the magnets; in repulsion, two plugs near the edges of the magnets; and with a single magnet, a plug close to the center of the magnet. The geometry of the magnets (h and l are the height and length of the magnets respectively) and their relative spacing s has a significant impact on the magnetic flux density. Its value inside a magnet increases with the h/l ratio. Consequently, bar magnets produce larger and more uniform values than flat magnets. The l/s ratio also influences the magnetic force value in the microchannel, both increasing concomitantly for all the configurations. In addition, a zero force zone in the middle appears in the attraction configuration as the l/s ratio increases, while with a single magnet, the number of maxima and minima goes from one to two, producing two focusing zones instead of only one.

Modélisation, analyse mathématique et simulation numérique de la dynamique des glaciers

Guillaume Jouvet

We address the free boundary problem that consists in finding the shape of a three dimensional glacier over a given period and under given climatic conditions. Glacier surface moves by sliding, internal shear and external exchange of mass. Ice is modelled as a non Newtonian fluid. Given the shape of the glacier, the velocity of ice is obtained by solving a stationary non-linear Stokes problem with a sliding law along the bedrock-ice interface. The shape of the glacier is updated by computing a Volume Of Fluid (VOF) function, which satisfies a transport equation. Climatic effects (accumulation and ablation of ice) are taken into account in the source term of this equation. A decoupling algorithm with a two-grid method allows the velocity of ice and the VOF to be computed using different numerical techniques, such that a Finite Element Method (FEM) and a characteristics method. On a theoretical level, we prove the well-posedness of the non-linear Stokes problem. A priori estimates for the convergence of the FEM are established by using a quasi-norm technique. Eventually, convergence of the linearisation schemes, such that a fixed point method and a Newton method, is proved. Several applications demonstrate the potential of the numerical method to simulate the motion of a glacier during a long period. The first one consists in the simulation of Rhone et Aletsch glacier from 1880 to 2100 by using climatic data provided by glaciologists. The glacier reconstructions over the last 120 years are validated against measurements. Afterwards, several different climatic scenarios are investigated in order to predict the shape the glaciers until 2100. A dramatic retreat during the 21st century is anticipated for both glaciers. The second application is an inverse problem. It aims to find a climate parametrization allowing a glacier to fit some of its moraines. Two other aspects of glaciology are also addressed in this thesis. The first one consists in modeling and in simulating ice collapse during the calving process. The previous ice flow model is supplemented by a Damage variable which describes the presence of micro crack in ice. An additional numerical scheme allows the Damage field to be solved and a two dimensional simulation of calving to be performed. The second problem aims to prove the existence of stationary ice sheet when considering shallow ice model and a simplified geometry. Numerical investigation confirms the theoretical result and shows physical properties of the solution.

EPFL2010

Adaptive analysis-aware defeaturing

Ondine Gabrielle Chanon

Removing geometrical details from a complex domain is a classical operation in computer aided design for simulation and manufacturing. This procedure simplifies the meshing process, and it enables faster simulations with less memory requirements. However, depending on the partial differential equation that one wants to solve in the geometrical model of interest, removing some important geometrical features may greatly impact the solution accuracy. For instance, in solid mechanics simulations, such features can be holes or fillets near stress concentration regions. Unfortunately, the effect of geometrical simplification on the accuracy of the problem solution is often neglected, because its analysis is a time-consuming task that is often performed manually, based on the expertise of engineers. It is therefore important to have a better understanding of the effect of geometrical model simplification, also called defeaturing, to improve our control on the simulation accuracy along the design and analysis phase.In this thesis, we formalize the process of defeaturing, and we analyze its impact on the accuracy of solutions of some partial differential problems. To achieve this goal, we first precisely define the error between the problem solution defined in the exact geometry, and the one defined in the simplified geometry. Then, we introduce an a posteriori estimator of the energy norm of this error. This allows us to reliably and efficiently control the error coming from the addition or the removal of geometrical features. We subsequently consider a finite element approximation of the defeatured problem, and the induced numerical error is integrated to the proposed defeaturing error estimator. In particular, we address the special case of isogeometric analysis based on (truncated) hierarchical B-splines, in possibly trimmed and multipatch geometries. In this framework, we derive a reliable a posteriori estimator of the overall error, i.e., of the error between the exact solution defined in the exact geometry, and the numerical solution defined in the defeatured geometry.We then propose a two-fold adaptive strategy for analysis-aware defeaturing, which starts by considering a coarse mesh on a fully-defeatured computational domain. On the one hand, the algorithm performs classical finite element mesh refinements in a (partially) defeatured geometry. On the other hand, the strategy also allows for geometrical refinement. That is, at each iteration, the algorithm is able to choose which missing geometrical features should be added to the simplified geometrical model, in order to obtain a more accurate solution.Throughout the thesis, we validate the presented theory, the properties of the aforementioned estimators and the proposed adaptive strategies, thanks to an extensive set of numerical experiments.

EPFL2022

Adaptive analysis-aware defeaturing

Ondine Gabrielle Chanon

EPFL2022

Modélisation, analyse mathématique et simulation numérique de la dynamique des glaciers

Guillaume Jouvet

EPFL2010