Numerical simulation of orbitally shaken viscous fluids with free surface

被引:27
作者
Discacciati, Marco [3 ]
Hacker, David [2 ]
Quarteroni, Alfio [1 ,4 ]
Quinodoz, Samuel [1 ]
Tissot, Stephanie [2 ]
Wurm, Florian M. [2 ]
机构
[1] Ecole Polytech Fed Lausanne, Chair Modelling & Sci Comp, MATHICSE, CH-1015 Lausanne, Switzerland
[2] Ecole Polytech Fed Lausanne, Lab Cellular Biotechnol, CH-1015 Lausanne, Switzerland
[3] Univ Politecn Cataluna, UPC BarcelonaTech, ETSECCPB, Lab Calcul Numer LaCaN, E-08034 Barcelona, Spain
[4] Politecn Milan, MOX, I-20133 Milan, Italy
基金
瑞士国家科学基金会;
关键词
Navier-Stokes; multi-phase flows; free surface; finite elements; level set; FINITE-ELEMENT-METHOD; LEVEL SET; 2-PHASE FLOW; APPROXIMATION; CELLS;
D O I
10.1002/fld.3658
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
Orbitally shaken bioreactors are an emerging alternative to stirred-tank bioreactors for large-scale mammalian cell culture, but their fluid dynamics is still not well defined. Among the theoretical and practical issues that remain to be resolved, the characterization of the liquid free surface during orbital shaking remains a major challenge because it is an essential aspect of gas transfer and mixing in these reactors. To simulate the fluid behavior and the free surface shape, we developed a numerical method based on the finite element framework. We found that the large density ratio between the liquid and the gas phases induced unphysical results for the free surface shape. We therefore devised a new pressure correction scheme to deal with large density ratios. The simulations operated with this new scheme gave values of wave amplitude similar to the ones measured experimentally. These simulations were used to calculate the shear stress and to study the mixing principle in orbitally shaken bioreactors. Copyright (C) 2012 John Wiley & Sons, Ltd.
引用
收藏
页码:294 / 315
页数:22
相关论文
共 31 条
[1]  
[Anonymous], SAND0662 SAND NAT LA
[2]  
[Anonymous], 1999, Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science
[3]  
Belytschko T, 2001, INT J NUMER METH ENG, V50, P993, DOI 10.1002/1097-0207(20010210)50:4<993::AID-NME164>3.0.CO
[4]  
2-M
[5]   STREAMLINE UPWIND PETROV-GALERKIN FORMULATIONS FOR CONVECTION DOMINATED FLOWS WITH PARTICULAR EMPHASIS ON THE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS [J].
BROOKS, AN ;
HUGHES, TJR .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 1982, 32 (1-3) :199-259
[6]   Edge stabilization for Galerkin approximations of convection-diffusion-reaction problems [J].
Burman, E ;
Hansbo, P .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2004, 193 (15-16) :1437-1453
[7]   Continuous interior penalty finite element method for Oseen's equations [J].
Burman, Erik ;
Fernandez, Miguel A. ;
Hansbo, Peter .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2006, 44 (03) :1248-1274
[8]   An enriched finite element method and level sets for axisymmetric two-phase flow with surface tension [J].
Chessa, J ;
Belytschko, T .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 2003, 58 (13) :2041-2064
[9]   Improving Eulerian two-phase flow finite element approximation with discontinuous gradient pressure shape functions [J].
Coppola-Owen, AH ;
Codina, R .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, 2005, 49 (12) :1287-1304
[10]   Mass preserving finite element implementations of the level set method [J].
Di Pietro, Daniele A. ;
Lo Forte, Stefania ;
Parolini, Nicola .
APPLIED NUMERICAL MATHEMATICS, 2006, 56 (09) :1179-1195