TY - JOUR

T1 - A finite volume formulation for simulating drop motion in a viscoelastic two-phase system

AU - You, R.

AU - Borhan, A.

AU - Haj-Hariri, H.

N1 - Funding Information:
We wish to acknowledge the support of NASA Glenn Research Center. This research was funded by NASA Grant NAG3-2760 with Dr. R. Balasubramaniam as the technical monitor.

PY - 2008/8

Y1 - 2008/8

N2 - We present a numerical method for simulating the motion of a drop in a viscoelastic two-phase system. The axisymmetric incompressible flow equations combined with viscoelastic constitutive equations are solved on an orthogonal curvilinear coordinate system using the finite volume method. The artificial compressibility factor is introduced in the continuity equation to effectively increase the convergence rate. A flux-difference splitting scheme is utilized in both the momentum and constitutive equations. This scheme has been used in high-speed aerodynamics to automatically handle the change of type of equations between elliptic and hyperbolic. The combination of artificial compressibility and the flux-difference splitting scheme is capable of overcoming the numerical instability caused by high elasticity of the viscoelastic phase. The differential constitutive equations of the FENE-CR model are employed to model the polymer stress terms in the momentum equations if either phase of the two-fluid system is viscoelastic. The free-surface boundary condition consists of the dynamic and kinematic conditions. These are used to update the drop shape during its time evolution. Conforming to the deforming drop shape, a boundary-fitted grid is generated at each time step. The grid is generated by solving elliptic partial differential equations.

AB - We present a numerical method for simulating the motion of a drop in a viscoelastic two-phase system. The axisymmetric incompressible flow equations combined with viscoelastic constitutive equations are solved on an orthogonal curvilinear coordinate system using the finite volume method. The artificial compressibility factor is introduced in the continuity equation to effectively increase the convergence rate. A flux-difference splitting scheme is utilized in both the momentum and constitutive equations. This scheme has been used in high-speed aerodynamics to automatically handle the change of type of equations between elliptic and hyperbolic. The combination of artificial compressibility and the flux-difference splitting scheme is capable of overcoming the numerical instability caused by high elasticity of the viscoelastic phase. The differential constitutive equations of the FENE-CR model are employed to model the polymer stress terms in the momentum equations if either phase of the two-fluid system is viscoelastic. The free-surface boundary condition consists of the dynamic and kinematic conditions. These are used to update the drop shape during its time evolution. Conforming to the deforming drop shape, a boundary-fitted grid is generated at each time step. The grid is generated by solving elliptic partial differential equations.

UR - http://www.scopus.com/inward/record.url?scp=47049084041&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=47049084041&partnerID=8YFLogxK

U2 - 10.1016/j.jnnfm.2007.12.002

DO - 10.1016/j.jnnfm.2007.12.002

M3 - Article

AN - SCOPUS:47049084041

VL - 153

SP - 109

EP - 129

JO - Journal of Non-Newtonian Fluid Mechanics

JF - Journal of Non-Newtonian Fluid Mechanics

SN - 0377-0257

IS - 2-3

ER -