Méthodes numériques liées à la distribution d'alumine dans une cuve d'électrolyse d'aluminium

Metallic aluminium plays a key role in our modern economy. Primary metallic aluminium is produced by the transformation of aluminium oxide using the Hall-Héroult industrial process. This process, which requires enormous quantities of energy, consists in performing the electrolysis of an aluminium oxide solute in large pots and with hundreds of thousands of amperes of electrical current.

The topic of this thesis is the study of some selected aspects of the modelisation of the electrolysis process from the point of view of numerical simulation. This thesis is divided in two parts.

The first part is focused on the numerical modelisation of the alumina powder dissolution and transport in the electrolytic bath as a function of the bath temperature. We provide a mathematical model for the transport and dissolution of the alumina powder, followed by its time and space discretisation by means of a finite element method. Finally, we study the behavior of this numerical model in the case of an industrial electrolysis pot.

The second part is devoted to the development of a numerical scheme for the approximation of the fluid flows in an electrolysis pot. The scheme relies on a Fourier basis decomposition of the unknowns. The amplitude of each Fourier component satisfies a partial differential equation which is explicitely derived. The solution of this equation is approximated by means of a finite element method. Finally, the approximate fluid flow obtained with this new method is compared with the solution provided by the reference model in an industrial

The topic of this thesis is the study of some selected aspects of the modelisation of the electrolysis process from the point of view of numerical simulation. This thesis is divided in two parts.

Reduction techniques for PDEs built upon Reduced Basis and Domain Decomposition Methods with applications to hemodynamics

Luca Pegolotti

This thesis focuses on the development and validation of a reduced order technique for cardiovascular simulations. The method is based on the combined use of the Reduced Basis method and a Domain Decomposition approach and can be seen as a particular implementation of the Reduced Basis Element method. Our contributions include the application to the unsteady three-dimensional Navier--Stokes equations, the introduction of a reduced coupling between subdomains, and the reconstruction of arteries with deformed elementary building blocks. The technique is divided into two main stages: the offline and the online phases. In the offline phase, we define a library of reference building blocks (e.g., tubes and bifurcations) and associate with each of these a set of Reduced Basis functions for velocity and pressure. The set of Reduced Basis functions is obtained by Proper Orthogonal Decomposition of a large number of flow solutions called snapshots; this step is expensive in terms of computational time. In the online phase, the artery of interest is geometrically approximated as a composition of subdomains, which are obtained from the parametrized deformation of the aforementioned building blocks. The local solution in each subdomain is then found as a linear combination of the Reduced Basis functions defined in the corresponding building block. The strategy to couple the local solutions is of utmost importance. In this thesis, we devise a nonconforming method for the coupling of Partial Differential Equations that takes advantage of the definition of a small number of Lagrange multiplier basis functions on the interfaces. We show that this strategy allows us to preserve the h-convergence properties of the discretization method of choice for the primal variable even when a small number of Lagrange multiplier basis functions is employed. Moreover, we test the flexibility of the approach in scenarios in which different discretization algorithms are employed in the subdomains, and we also use it in a fluid-structure interaction benchmark. The introduction of the Lagrange multipliers, however, is associated with stability problems deriving from the saddle-point structure of the global system. In our Reduced Order Model, the stability is recovered by means of supremizers enrichment. In our numerical simulations, we specifically focus on the effects of the Reduced Basis and geometrical approximations on the quality of the results. We show that the Reduced Order Model performs similarly to the corresponding high-fidelity one in terms of accuracy. Compared to other popular models for cardiovascular simulations (namely 1D models), it also allows us to compute a local reconstruction of the Wall-Shear Stress on the vessel wall. The speedup with respect to the Finite Element method is substantial (at least one order of magnitude), although the current implementation presents bottlenecks that are addressed in depth throughout the thesis.

EPFL2021

A mixture model to take into account diluted gas in liquid flow: applications to aluminium electrolysis

Emile Tryphon Pierre Soutter

Aluminium is a metal sought in the industry because of its various physical properties. It is produced by an electrolysis reduction process in large cells. In these cells, a large electric current goes through the electrolytic bath and the liquid aluminium. This electric current generates electromagnetic forces that set the bath and the aluminium into motion. Moreover, large quantities of carbon dioxide gas are produced through chemical reactions in the electrolytic bath: the presence of these gases alleviates the density of the liquid bath and changes the dynamics of the flow. Accurate knowledge of this fluid flow is essential to improve the efficiency of the whole process.The purpose of this thesis is to study and approximate the interaction of carbon dioxide with the fluid flow in the aluminium electrolysis process.In the first chapter of this work, a mixture-averaged model is developed for mixtures of gas and liquid. The model is based on the conservation of mass and momentum equations of the two phases, liquid and gas. By combining these equations, a system is established that takes into account the velocity of the liquid-gas mixture, the pressure, the gas velocity and the local gas concentration as unknowns.In the second chapter, a simplified problem is studied theoretically. It is shown that under the assumption that the gas concentration is small, the problem is well-posed. Moreover, we prove a priori error estimates of a finite element approximation of this problem.In the third chapter, we compare this liquid-gas model with a water column reactor experiment. Finally, the last chapter shows that the fluid flow is changed in aluminium electrolysis cells when we take into account the density of the bath reduced by carbon dioxide. These changes are quantified as being of the order of 30% and explain partially the differences between previous models and observations from Rio Tinto Aluminium engineers.