Numerical Simulation of Orbitally Shaken Reactors

Mammalian cell cultures have become a major topic of research in the biopharmaceutical industry. This kind of cells requires specific conditions to grow. In this thesis, we study the hydrodynamics of orbitally shaken reactors (OSR), a recently introduced kind of bioreactors for mammalian cell cultures that represents a simple to operate and cheap alternative to commonly used reactors such as stirred tanks. OSRs can provide suitable conditions for small scale cell cultures, however a deeper understanding of the principles governing the OSRs is required to exploit their full potentiality and proceed with scaling up. This work aims at shedding light into the mechanisms of the OSRs through computational fluid dynamics. OSRs are only partially filled with liquid medium, the remaining space is occupied by air. When an OSR is agitated, the interface between the two phases moves and creates different shapes. This interface is at the heart of the simulation of OSRs: not only its location is part of the problem, but it can also carry singularities. In particular, the pressure has usually a low regularity in the vicinity of the interface and numerical methods might underperform if the singularity is not treated in an appropriate manner. This motivated the study of an elliptic problem in a medium with an internal interface carrying discontinuities. In this work, we devise a novel method called SESIC to solve this kind of problem. It uses the a priori knowledge to improve the numerical accuracy in the vicinity of the interface by removing the singularities. We prove that this method yields optimal orders of convergence in H1 and L2 norm. Numerical tests also show that optimal orders can be obtained in the L∞ norm in some cases. Regularized integration is also investigated with the perspective of further simplifying the scheme. It is found that, if the regularization bandwidth is suitably defined, good approximations can still be obtained, even if the convergence rate is decreased. We apply then the methodology of the SESIC method to the approximation of the two-phase Navier-Stokes equations, which amounts to correct the pressure. If adapted integration is used with it, the density and viscosity can be kept discontinuous across the interface without creating spurious velocity, as shown by numerical experiments. The sharp treatment of the discontinuities improves the accuracy of the simulations by retaining the physical meaning of the phenomena independently of the mesh size. We also pay attention to the boundary conditions used, which must be suitably chosen to allow the interface motion but still reproduce the wall friction. We show that imposing the zero normal component of the velocity yields the best results for the no-penetration condition and that it must be employed with a correction term to avoid spurious velocities. Robin-type conditions are used for the tangential components to recover the no-slip condition far from the contact line. Specific tests a