Skip to main content Accessibility help
×
Home

Application of the LS-STAG Immersed Boundary/Cut-Cell Method to Viscoelastic Flow Computations

  • Olivier Botella (a1), Yoann Cheny (a1), Farhad Nikfarjam (a1) and Marcela Stoica (a1)

Abstract

This paper presents the extension of a well-established Immersed Boundary (IB)/cut-cell method, the LS-STAG method (Y. Cheny & O. Botella, J. Comput. Phys. Vol. 229, 1043-1076, 2010), to viscoelastic flow computations in complex geometries. We recall that for Newtonian flows, the LS-STAG method is based on the finite-volume method on staggered grids, where the IB boundary is represented by its level-set function. The discretization in the cut-cells is achieved by requiring that global conservation properties equations be satisfied at the discrete level, resulting in a stable and accurate method and, thanks to the level-set representation of the IB boundary, at low computational costs.

In the present work, we consider a general viscoelastic tensorial equation whose particular cases recover well-known constitutive laws such as the Oldroyd-B, White-Metzner and Giesekus models. Based on the LS-STAG discretization of the Newtonian stresses in the cut-cells, we have achieved a compatible velocity-pressure-stress discretization that prevents spurious oscillations of the stress tensor. Applications to popular benchmarks for viscoelastic fluids are presented: the four-to-one abrupt planar contraction flows with sharp and rounded re-entrant corners, for which experimental and numerical results are available. The results show that the LS-STAG method demonstrates an accuracy and robustness comparable to body-fitted methods.

Copyright

Corresponding author

*Corresponding author. Email addresses: olivier.botella@univ-lorraine.fr (O. Botella), yoann.cheny@univ-lorraine.fr (Y. Cheny)

References

Hide All
[1] Mittal, R. and Iaccarino, G.. Immersed boundary methods. Annual Review of Fluid Mechanics, 37, 239261, 2005.
[2] Balaras, E.. Modeling complex boundaries using an external force field on fixed Cartesian grids in large-eddy simulations. Computers Fluids, 33, 375404, 2004.
[3] Parnaudeau, P., Carlier, J., Heitz, D., and Lamballais, E.. Experimental and numerical studies of the flow over a circular cylinder at reynolds number 3900. Physics of Fluids, 20 085101, 2008.
[4] Yang, J., Preidikman, S., and Balaras, E.. A strongly coupled, embedded-boundary method for fluid-structure interactions of elastically mounted rigid bodies. Journal of Fluids and Structures, 24(2), 167182, 2008.
[5] Monasse, L., Daru, V., Mariotti, C., and Piperno, S.. A conservative immersed boundary method for fluid-structure interaction in the compressible inviscid case. Computational Fluid Dynamics 2010, Kuzmin, A. (Ed.), Springer, 2010.
[6] De Tullio, M. D., Cristallo, A., Balaras, E., and Verzicco, R.. Direct numerical simulation of the pulsatile flow through an aortic bileaflet mechanical heart valve. Journal of Fluid Mechanics, 622, 259290, 2009.
[7] Bergmann, M. and Iollo, A.. Modeling and simulation of fish-like swimming. Journal of Computational Physics, 230, 329348, 2011.
[8] Teran, J., Fauci, L., and Shelley, M.. Peristaltic pumping and irreversibility of a Stokesian viscoelastic fluid. Physics of Fluids, 20, 073101, 2008.
[9] Nonaka, A., Trebotich, D., Miller, G. H., Graves, D. T., and Colella, P.. A higher-order upwind method for viscoelastic flow. Communications in Applied Mathematics and Computational Science, 4, 5783, 2009.
[10] Crochet, M. J., Davies, A. R., and Walters, K.. Numerical Simulation of Non-Newtonian Flow. Elsevier, Amsterdam, 1984.
[11] Walters, K. and Webster, M. F.. The distinctive CFD challenges of computational rheology. International Journal For Numerical Methods in Fluid, 43, 577596, 2003.
[12] Owens, R. G. and Phillips, T. N.. Computational Rheology. Imperial College Press, London, 2002.
[13] Kang, S., Iaccarino, G., and Moin, P.. Accurate and efficient immersed-boundary interpolations for viscous flows. In Center for Turbulence Research Briefs, NASA Ames/Stanford University, pp. 3143, 2004.
[14] Muldoon, F. and Acharya, S.. A divergence-free interpolation scheme for the immersed boundary method. International Journal for Numerical Methods in Fluids, 56, 18451884, 2008.
[15] Guy, R. D. and Hartenstine, D. A.. On the accuracy of direct forcing immersed boundary methods with projection methods. Journal of Computational Physics, 229(7), 24792496, 2010.
[16] Cheny, Y. and Botella, O.. An immersed boundary/level-set method for incompressible viscous flows in complex geometries with good conservation properties. European Journal of Mechanics, 18, 561587, 2009.
[17] Cheny, Y. and Botella, O.. The LS-STAG method: A new immersed boundary / level-set method for the computation of incompressible viscous flows in complex moving geometries with good conservation properties. Journal of Computational Physics, 229, 10431076, 2010.
[18] Botella, O., Ait-Messaoud, M., Pertat, A., Rigal, C., and Cheny, Y.. The LS-STAG immersed boundary method for non-Newtonian flows in irregular geometries: Flow of shear-thinning liquids between eccentric rotating cylinders. Theoretical and Computational Fluid Dynamics, 29, 93110, 2015.
[19] Bird, R. B., Armstrong, R. C., and Hassager, O.. Dynamics of Polymeric Liquids. Wiley-Interscience, New-York, 1987.
[20] Baaijens, F. P. T.. Mixed finite element methods for viscoelastic flow analysis: A review. Journal of Non-Newtonian Fluid Mechanics, 79, 361385, 1998.
[21] Darwish, M. S. and Whiteman, J. R.. Numerical modeling of viscoelastic liquids using a finite-volume method. Journal of Non-Newtonian Fluid Mechanics, 45, 311337, 1992.
[22] Gerritsma, M.. Time dependant numerical simulations of a viscoelastic fluid on a staggered grid. PhD thesis, University of Groningen, The Netherlands, 1996.
[23] Van Kemenade, V. and Deville, M. O.. Application of spectral elements to viscoelastic creeping flows. Journal of Non-Newtonian Fluid Mechanics, 51(3), 277308, 1994.
[24] Fattal, R. and Kupferman, R.. Time-dependent simulation of viscoelastic flows at high Weissenberg number using the log-conformation representation. Journal of Non-Newtonian Fluid Mechanics, 126, 2337, 2005.
[25] Hao, J., Pan, T.-W., Glowinski, R., and Joseph, D. D.. A fictitious domain/distributed Lagrange multiplier method for the particulate flow of Oldroyd-B fluids: A positive definiteness preserving approach. Journal of Non-Newtonian Fluid Mechanics, 156, 95111, 2009.
[26] Matallah, H., Townsend, P., and Webster, M. F.. Recovery and stress-splitting schemes for viscoelastic flows. Journal of Non-Newtonian Fluid Mechanics, 75, 139166, 1998.
[27] Alves, M. A., Pinho, F. T., and Oliveira, P. J.. Effect of a high-resolution differencing scheme on finite-volume predictions of viscoelastic flows. Journal of Non-Newtonian Fluid Mechanics, 93, 287314, 2000.
[28] Kim, J. M., Kim, C., Kim, J. H., Chung, C., Ahn, K. H., and Lee, S. J.. High-resolution finite element simulation of 4:1 planar contraction flow of viscoelastic fluid. Journal of Non-Newtonian Fluid Mechanics, 129, 2337, 2005.
[29] Dupret, F., Marchal, J. M., and Crochet, M. J.. On the consequence of discretization errors in the numerical calculation of viscoelastic flow. Journal of Non-Newtonian Fluid Mechanics, 18, 173186, 1985.
[30] Aboubacar, M., Matallah, H., Tamaddon-Jahromi, H. R., and Webster, M. F.. Numerical prediction of extensional flows in contraction geometries: Hybrid finite volume/element method. Journal of Non-Newtonian Fluid Mechanics, 104, 125164, 2002.
[31] Phillips, T. N. and Williams, A. J.. Viscoelastic flow through a planar contraction using a semi-Lagrangian finite volume method. Journal of Non-Newtonian Fluid Mechanics, 87, 215246, 1999.
[32] Yoo, J. Y. and Na, Y.. A numerical study of the planar contraction flow of a viscoelastic fluid using the SIMPLER algorithm. Journal of Non-Newtonian Fluid Mechanics, 39, 89106, 1991.
[33] Oliveira, P. J., Pinho, F. T., and Pinto, G. A.. Numerical simulation of non-linear elastic flows with a general collocated finite-volume method. Journal of Non-Newtonian Fluid Mechanics, 79, 143, 1998.
[34] Keiller, R. A.. Entry-flow calculations for the Oldroyd-B and FENE equations. Journal of Non-Newtonian Fluid Mechanics, 46, 143178, 1993.
[35] Quinzani, L. M., Amstrong, R. C., and Brown, R. A.. Birefringence and laser-Doppler velocimetry (LDV) studies of viscoelastic flow through planar contraction. Journal of Non-Newtonian Fluid Mechanics, 52, 136, 1994.
[36] Quinzani, L. M., Amstrong, R. C., and Brown, R. A.. Used of coupled birefringence and LDV studies of flow through planar contraction to test constitutive equations for concentrated polymer solutions. Journal of Rheology, 39, 12011228, 1995.
[37] Harlow, F. H. and Welch, J. E.. Numerical calculation of time-dependent viscous incompressible flow of fluid with free surfaces. Physics of Fluids, 8, 21812189, 1965.
[38] Phillips, T. N. and Williams, A. J.. Comparison of creeping and inertial flow of an Oldroyd-B fluid through planar and axisymmetric contractions. Journal of Non-Newtonian Fluid Mechanics, 108, 2547, 2002.
[39] Mompean, G. and Deville, M.. Unsteady finite volume simulation of Oldroyd-B fluid through a three-dimensional planar contraction. Journal of Non-Newtonian Fluid Mechanics, 872, 253279, 1997.
[40] Al Moatassime, H. and Jouron, C.. A multigrid method for solving steady viscoelastic fluid flow. Computer Methods in Applied Mechanics and Engineering, 190(31), 40614080, 2001.
[41] Sasmal, G. P.. A finite volume approach for calculation of viscoelastic flow through an abrupt axisymmetric contraction. Journal of Non-Newtonian Fluid Mechanics, 56(1), 1547, 1995.
[42] Xia, H., Tucker, P. G., and Dawes, W. N.. Level sets for CFD in aerospace engineering. Progress in Aerospace Sciences, 46, 274283, 2010.
[43] Verstappen, R. W. C. P. and Veldman, A. E. P.. Symmetry-preserving discretization of turbulent flow. Journal of Computational Physics, 187, 343368, 2003.
[44] van der Plas, P., van der Heiden, H. J. L., Veldman, A. E. P., Luppes, R., and Verstappen, R. W. C. P.. Efficiently simulating viscous flow effects by means of regularization turbulence modeling and local grid refinement. Seventh International Conference on Computational Fluid Dynamics (ICCFD7), Hawaii, paper ICCFD-2503, 2012.
[45] Puzikova, V. and Marchevsky, I.. Application of the LS-STAG immersed boundary method for numerical simulation in coupled aeroelastic problems. Proceedings of the 11th World Congress on ComputationalMechanics (WCCMXI), Oñate, E., Oliver, J. and Huerta, A. (Eds), 20-25 July 2014, Barcelona, Spain., 2014.
[46] Puzikova, V. and Marchevsky, I.. Extension of the LS-STAG cut-cell immersed boundary method for RANS-based turbulence models. Proceedings of the International Summer School-Conference “Advanced Problems in Mechanics”, June 30-July 5, St. Petersburg, Russia., 2014.
[47] Saramito, P.. Efficient simulation of nonlinear viscoelastic fluid flows. Journal of Non-Newtonian Fluid Mechanics, 60(2-3), 199223, 1995.
[48] Alves, M. A., Oliveira, P. J., and Pinho, F. T.. Benchmark solutions for the flow of Oldroyd-B and PTT fluids in planar contractions. Journal of Non-Newtonian Fluid Mechanics, 110, 4575, 2003.
[49] Edussuriya, S. S., Williams, A. J., and Bailey, C.. A cell-centred finite volume method for modelling viscoelastic flow. Journal of Non-Newtonian Fluid Mechanics, 117, 4761, 2004.
[50] Li, X., Han, X., and Wang, X.. Numerical modeling of viscoelastic flows using equal low-order finite elements. Computer Methods in Applied Mechanics and Engineering, 199(9-12), 570581, 2010.
[51] Azaiez, J., Guénette, R., and Aït-Kadi, A.. Entry flow calculations using multi-mode models. Journal of Non-Newtonian Fluid Mechanics, 66, 271281, 1996.
[52] Aboubacar, M. and Webster, M. F.. A cell-vertex finite volume/element method on triangles for abrupt contraction viscoelastic flows. Journal of Non-Newtonian Fluid Mechanics, 98, 83106, 2001.
[53] Orlanski, I.. A simple boundary condition for unbounded hyperbolic flows. Journal of Computational Physics, 21, 251269, 1976.
[54] Dieci, L. and Eirola, T.. Positive definiteness in the numerical solution of Riccati differential equations. Numerische Mathematik, 67, 303313, 1994.

Keywords

MSC classification

Application of the LS-STAG Immersed Boundary/Cut-Cell Method to Viscoelastic Flow Computations

  • Olivier Botella (a1), Yoann Cheny (a1), Farhad Nikfarjam (a1) and Marcela Stoica (a1)

Metrics

Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed