Flow simulation and hemolysis modeling for a blood centrifuge device

A common blood separation technique is to accelerate the sedimentation of blood constituents by centrifugation and then collect these constituents separately. When subjected to conditions of excessive stress, hemolysis can however occur, whereby the membrane of red blood cells ruptures and free hemoglobin is released into the plasma. In this paper, results of numerical simulations of the complex turbulent flow in the collection region of a blood centrifuge device are presented. These simulation results provide a detailed qualitative evaluation of the complex physical phenomena causing hemolysis. For a quantitative analysis, a 5-parameter model is proposed for the production of hemolysis. This model has been calibrated for hemolysis occurring during the collection process using experimentally measured values for different collector geometry, hematocrit level and rotational speed of the blood centrifuge device. (C) 2013 Elsevier Ltd. All rights reserved.

Approximation numérique des écoulements turbulents dans des cuves d'électrolyse de l'aluminium

Jonathan Rochat

Aluminium is a metal extracted from bauxite ore using electrolysis process in cells of big size. A huge electric current goes in the cell through an electrolytic bath and aluminium liquid. These currents generate strong magnetic forces that allow the bath and the aluminium to move. A good knowledge of these turbulent flows is very important to optimize the process. The purpose of this thesis is to study and simulate turbulent flows in the aluminium smelting process. These flows are solved numerically with a finite element method. In particular, the Navier-Stokes equations for bifluid flows with free moving interface are solved numerically. In the first part of this work, we develop some mixing-length models that take into account the effects of the wall. A theoretical mathematical study shows the validity of these models and we give some recommendation on the choice of the parameters of the computation. In the second part we study the resolution of the Navier-Stokes equations. The study focuses on algorithms that decouple the computation of the speed and pressure, commonly called projection method or Chorin-Temam algorithm. The final section provides answers on the relevance of wall modelling and projection methods in numerical simulation of turbulent flows in the aluminium smelting process. In particular, we obtain a numerical model that produces a realistic flow with a reasonable CPU time and we discuss the choice of certain parameters involved in the different models.

EPFL2016

Numerical modeling of cohesive sediment transport and bed morphology in estuaries

Célestin Leupi

Two major lines of investigation have been pursued in this thesis: (1) More efficient, robust and realistic numerical techniques are designed for the simulation of complex turbulent fluid flows; (2) A new algorithm and its analysis is performed in the context of multiphasic fluid flow, for a cohesive fine-grained sediment (fluid mud) transport in estuaries. Estuaries exist between marine and freshwater system where waters of different physical, chemical and biological composition meet, combine and disperse primarily due to tidal influences. In the present thesis, the behavior of cohesive sediment in estuaries is reviewed based on the existing literature. Basic theories and recent developments are introduced to describe the formation of fluid mud from a very dilute suspension and its motion down a natural river bed with complex bathymetry. The present work contributes to the numerical simulation of complex turbulent multiphasic fluid flows encountered in estuarine channels, with the aim of the better understanding of the underlying physical processes as well as predicting realistically the cohesive sediment transport and bed morphology in such a zone. The model is based on the mass preserving method by using the so-called Raviart-Thomas finite element on the unstructured mesh in the horizontal plane. In the vertical, the computational domain is divided into number of layers at predefined heights and the method uses a conventional conforming P1 finite element scheme, with the advantage that the lowermost and uppermost layers variable height allow a faithful representation of the time-varying bed and free surface, respectively. Concerning the modeling of turbulence, the research effort focuses on the turbulence two-equation k - ε closure for the vertical parameterization of eddy viscosity. More precisely, a robust up-to-date algorithm is used for this issue. The new methodology is developed with the aim to account for more general relevant effects in the closure. The proposed model offers the capability to cope with the stiffness problem introduced by the large difference between the solid phase flow time scale and the hydrodynamic one, by using a sub-cycling strategy, whereas the splitting scheme is adopted with the aim of stability and the positivity of the relevant turbulent variables. The flexibility of the model and its performance are evaluated on several free-surface flow configurations with increasing complexity : homogeneous unsteady non-uniform flows in plane open channel flows, U-shaped (193°) curved open channel flow. Concerning the cohesive sediment transport, most of the existing models in the literature assume the analogous transport characteristics with that of the coarse sediment and adopt the relevant developed sediment transport for the latter to treat the former. Moreover, these existing models do not account for the consolidation of the mud-bed. The present research effort focused on a novel methodology based on the realistic empirical relationships, which accounts for the mutually exclusive processes for re-suspension and/or erosion and deposition of fine sediment, whereas only a limited range of bed shear stresses is allowed for simultaneous erosion and deposition to occur. Hence, on this basis, the new model investigated the bed morphology evolution by taking into account of the fluidization and/or consolidation of the fluid mud, which was handled by modeling the bed in three layers: (i) the mud-bed layer, (ii) the partially consolidated bed and (iii) the fully consolidated bed. The prediction of deposition/re-suspension using these two different methods lead to a non negligible difference in the results. Therefore, investigation of the true mechanism of erosion/deposition processes for cohesive sediments and their implementation in the numerical model is very important. This suggests that a realistic prediction must account for the fresh mud-bed re-suspension once deposited, as well as the consolidation and/or fluidization of the mud-bed deposits. Finally, the capability and improvements of the model are demonstrated, and its predicting performance is successfully evaluated by applying it to the simulation of the Po River Estuary (PRE) in Italy, which is the main source of river water discharge into the Northern Adriatic Sea. The analysis showed that the consolidation/fluidization process at the bed may have important influence on the prediction of bed morphology evolution. The three-layer approach used in this thesis is a first attempt to model these processes in detail within a numerical model.

EPFL2005

Simulation of shear-driven flows

Roland Bouffanais

The research work reported in the present dissertation is aimed at the analysis of complex physical phenomena involving instabilities and nonlinearities occurring in fluids through state-of-the-art numerical modeling. Solutions of intricate fluid physics problems are devised in two particularly arduous situations: fluid domains with moving boundaries and the high-Reynolds-number regime dominated by nonlinear convective effects. Shear-driven flows of incompressible Newtonian fluids enclosed in cavities of varying geometries are thoroughly investigated in the two following frameworks: transition with a free surface and confined turbulence. The physical system we consider is made of an incompressible Newtonian fluid filling a bounded, or partially bounded cavity. A series of shear-driven flows are easily generated by setting in motion some part of the container boundary. These driven-cavity flows are not only technologically important, they are of great scientific interest because they display almost all physical fluid phenomena that can possibly occur in incompressible flows, and this in the simplest geometrical settings. Thus corner eddies, secondary flows, longitudinal vortices, complex three-dimensional patterns, chaotic particle motions, nonuniqueness, transition, and turbulence all occur naturally and can be studied in the same geometry. This facilitates the comparison of results from experiments, analysis, and computation over the whole range of Reynolds numbers. The flows under consideration are part of a larger class of confined flows driven by linear or angular momentum gradients. This dissertation reports a detailed study of a novel numerical method developed for the simulation of an unsteady free-surface flow in three-space-dimensions. This method relies on a moving-grid technique to solve the Navier-Stokes equations expressed in the arbitrary Lagrangian-Eulerian (ALE) kinematics and discretized by the spectral element method. A comprehensive analysis of the continuous and discretized formulations of the general problem in the ALE frame, with nonlinear, non-homogeneous and unsteady boundary conditions is presented. In this dissertation, we first consider in the internal turbulent flow of a fluid enclosed in a bounded cubical cavity driven by the constant translation of its lid. The solution of this flow relied on large-eddy simulations, which served to improve our physical understanding of this complex flow dynamics. Subsequently, a novel subgrid model based on approximate deconvolution methods coupled with a dynamic mixed scale model was devised. The large-eddy simulation of the lid-driven cubical cavity flow based on this novel subgrid model has shown improvements over traditional subgrid-viscosity type of models. Finally a new interpretation of approximate deconvolution models when used with implicit filtering as a way to approximate the projective grid filter was given. This led to the introduction of the grid filter models. Through the use of a newly-developed method of numerical simulation, in this dissertation we solve unsteady flows with a flat and moving free-surface in the transitional regime. These flows are the incompressible flow of a viscous fluid enclosed in a cylindrical container with an open top surface and driven by the steady rotation of the bottom wall. New flow states are investigated based on the fully three-dimensional solution of the Navier-Stokes equations for these free-surface cylindrical swirling flows, without resorting to any symmetry properties unlike all other results available in the literature. To our knowledge, this study delivers the most general available results for this free-surface problem due to its original mathematical treatment. This second part of the dissertation is a basic research task directed at increasing our understanding of the influence of the presence of a free surface on the intricate transitional flow dynamics of shear-driven flows.