Hostname: page-component-cd9895bd7-lnqnp Total loading time: 0 Render date: 2024-12-27T11:20:50.418Z Has data issue: false hasContentIssue false

Simulating Three-Dimensional Free Surface Viscoelastic Flows using Moving Finite Difference Schemes

Published online by Cambridge University Press:  28 May 2015

Yubo Zhang
Affiliation:
Department of Computer Science, 2063 Kemper Hall, University of California at Davis, One Shields Avenue, Davis, CA 95616-8562, USA
Tao Tang
Affiliation:
Department of Mathematics, Hong Kong Baptist University, Kowloon Tong, Kowloon, Hong Kong
Get access

Abstract

An efficient finite difference framework based on moving meshes methods is developed for the three-dimensional free surface viscoelastic flows. The basic model equations are based on the incompressible Navier-Stokes equations and the Oldroyd-B constitutive model for viscoelastic flows is adopted. A logical domain semi-Lagrangian scheme is designed for moving-mesh solution interpolation and convection. Numerical results show that harmonic map based moving mesh methods can achieve better accuracy for viscoelastic flows with free boundaries while using much less memory and computational time compared to the uniform mesh simulations.

Type
Research Article
Copyright
Copyright © Global Science Press Limited 2011

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1] anderson, D.M., McFadden, G.B., and Wheeler, A.A.. Diffuse-interface methods in fluid mechanics. Annual Reviews in Fluid Mechanics, 30(1):139165, 1998.CrossRefGoogle Scholar
[2] Brackbill, J.U., Kothe, D.B., and Zemach, C.. A continuum method for modeling surface tension. Journal of Computational Physics, 100(2) :335354, 1992.Google Scholar
[3] Ceniceros, H.D. and Hou, T.Y.. An efficient dynamically adaptive mesh for potentially singular solutions. Journal of Computational Physics, 172(2):609639, 2001.Google Scholar
[4] Chorin, A.J.. Numerical solution of the Navier-Stokes equations. Mathematics of Computation, 22(104) :745762, 1968.Google Scholar
[5] Courant, R., Isaacson, E., and Rees, M.. On the solution of nonlinear hyperbolic differential equations by finite differences. Commun. Pure Appl. Math., 5:243255, 1952.CrossRefGoogle Scholar
[6] Di, Y., Li, R., and Tang, T.. A general moving mesh framework in 3D and its application for simulating the mixture of multi-phase flows. Commun. Comput. Phys, 3:582602, 2008.Google Scholar
[7] Di, Y., Li, R., Tang, T., and Zhang, P.. Moving mesh finite element methods for the incompressible Navier-Stokes equations. SIAM J. Sci. Comput., 26(3):10361056, 2005.CrossRefGoogle Scholar
[8] Y Di, Li, R., Tang, T., and Zhang, P.. Level set calculations for incompressible two-phase flows on a dynamically adaptive grid. J. Sci. Comput, 31(1-2):7598, 2007.Google Scholar
[9] Dupont, T.F. and Liu., Y Back and forth error compensation and correction methods for removing errors induced by uneven gradients of the level set function. Journal of Computational Physics, 190(1):311324, 2003.CrossRefGoogle Scholar
[10] Enright, D., Fedkiw, R., Ferziger, J., and Mitchell, I.. A hybrid particle level set method for improved interface capturing. Journal of Computational Physics, 183(1):83116, 2002.CrossRefGoogle Scholar
[11] Feng, W. M., Yu, P., Hu, S. Y., Liu, Z. K., Du, Q., and Chen, L. Q.. A fourier spectral moving mesh method for the cahn-hilliard equation with elasticity. Commun. Comput. Phys, 5:582599, 2009.Google Scholar
[12] Gibou, F., Fedkiw, R.P., Cheng, L.T, and Kang, M.. A second-order-accurate symmetric discretization of the Poisson equation on irregular domains. Journal of Computational Physics, 176(1) :205227, 2002.Google Scholar
[13] Hirt, C.W. and Nichols, B.D.. Volume of fluid/VOF/ method for the dynamics of free boundaries. Journal of Computational Physics, 39:201225, 1981.CrossRefGoogle Scholar
[14] Hu, X.-L., Li, R., and Tang, T.. A multi-mesh adaptive finite element approximation to phase field models. Commun. Comput. Phys, 5:10121029, 2009.Google Scholar
[15] Jacqmin, D.. Calculation of two-phase Navier-Stokes flows using phase-field modeling. Journal of Computational Physics, 155(1):96127, 1999.CrossRefGoogle Scholar
[16] Li, R., Tang, T., and Zhang, P.. A moving mesh finite element algorithm for singular problems in two and three space dimensions. Journal of Computational Physics, 177(2) :365393, 2002.Google Scholar
[17] Li, R., Zhang, P., and Tang, T.. Moving mesh methods in multiple dimensions based on harmonic maps. Journal of Computational Physics, 170(2) :562588, 2001.Google Scholar
[18] Oldroyd, J.G.. Non-Newtonian flow of liquids and solids. Rheology: Theory and Applications, 1:653682, 1956.Google Scholar
[19] Oran, E.S. and Boris, J.P.. Numerical Simulation of Reactive Flow. Cambridge University Press, 2001.Google Scholar
[20] Osher, S. and Fedkiw, R.. Level Set Methods and Dynamic Implicit Surfaces. Applied Mathematical Sciences, 2003.Google Scholar
[21] Osher, S. and Sethian, J.. Fronts propagating with curvature-dependent speed- Algorithms based on Hamilton-Jacobi formulations. Journal of Computational Physics, 79:1249, 1988.CrossRefGoogle Scholar
[22] Qiao, Z.-H.. Numerical investigations of the dynamical behaviors and instabilities for the gierer-meinhardt system. Communications in Computational Physics, 3:406426, 2008.Google Scholar
[23] Scardovelli, R. and Zaleski, S.. Direct numerical simulation of free-surface and interfacial flow. Annual Review of Fluid Mechanics, 31(1):567603, 1999.Google Scholar
[24] Selle, A., Fedkiw, R., Kim, B.M., Liu, Y., and Rossignac, J.. An unconditionally stable MacCormack method. Journal of Scientific Computing, 35(2):350371, 2008.Google Scholar
[25] Sethian, J.A.. Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science. Cambridge University Press, 1999.Google Scholar
[26] Sussman, M. and Puckett, E.G.. A coupled level set and volume-of-fluid method for computing 3D and axisymmetric incompressible two-phase flows. Journal of Computational Physics, 162(2):301337, 2000.Google Scholar
[27] Sussman, M., Smereka, P., and Osher, S.. A level set approach for computing solutions to incompressible two-phase flow. Journal of Computational Physics, 114(1):146159, 1994.Google Scholar
[28] Tang, H.-Z. and Tang, T.. Adaptive mesh methods for one-and two-dimensional hyperbolic conservation laws. SIAM Journal on Numerical Analysis, 41(2):487515, 2004.Google Scholar
[29] Trebotich, D., Colella, P., and Miller, G.H.. A stable and convergent scheme for viscoelastic flow in contraction channels. Journal of Computational Physics, 205(1):315–342, 2005.CrossRefGoogle Scholar
[30] Dam, A. van and Zegeling, P.A.. Balanced monitoring of flow phenomena in moving mesh methods. Commun. Comput. Phys, 7:138170, 2010.Google Scholar
[31] Wang, H.-Y. and Li, R.. Mesh sensitivity for numerical solutions of phase-field equations using r-adaptive finite element methods. Communications in Computational Physics, 3:357375, 2008.Google Scholar
[32] Wang, H.-Y., Li, R., and Tang, T.. Efficient computation of dendritic growth with r-adaptive finite element methods. Journal of Computational Physics, 227(12):59846000, 2008.Google Scholar
[33] Yokoi, K.. A numerical method for free-surface flows and its application to droplet impact on a thin liquid layer. Journal of Scientific Computing, 35(2):372396, 2008.Google Scholar
[34] Yue, P., Zhou, C., Feng, J.J., Ollivier-Gooch, C.F., and Hu, H.H.. Phase-field simulations of interfacial dynamics in viscoelastic fluids using finite elements with adaptive meshing. Journal of Computational Physics, 219(1):4767, 2006.Google Scholar