0
Research Papers: Gas Turbines: Structures and Dynamics

The Depth-Averaged Numerical Simulation of Laminar Thin-Film Flows With Capillary Waves

[+] Author and Article Information
Bruce Kakimpa

Gas Turbine and Transmissions Research
Centre (G2TRC),
University of Nottingham,
Nottingham, NG7 2RD, UK
e-mail: bruce.kakimpa@nottingham.ac.uk

Herve Morvan

Professor
Gas Turbine and Transmissions Research
Centre (G2TRC),
University of Nottingham,
Nottingham, NG7 2RD, UK
e-mail: herve.morvan@nottingham.ac.uk

Stephen Hibberd

Associate Professor
Gas Turbine and Transmissions Research
Centre (G2TRC),
University of Nottingham,
Nottingham, NG7 2RD, UK
e-mail: stephen.hibberd@nottingham.ac.uk

1Corresponding author.

Contributed by the Structures and Dynamics Committee of ASME for publication in the JOURNAL OF ENGINEERING FOR GAS TURBINES AND POWER. Manuscript received March 15, 2016; final manuscript received April 6, 2016; published online May 24, 2016. Editor: David Wisler.

J. Eng. Gas Turbines Power 138(11), 112501 (May 24, 2016) (10 pages) Paper No: GTP-16-1105; doi: 10.1115/1.4033471 History: Received March 15, 2016; Revised April 06, 2016

Thin-film flows encountered in engineering systems such as aero-engine bearing chambers often exhibit capillary waves and occur within a moderate to high Weber number range. Although the depth-averaged simulation of these thin-film flows is computationally efficient relative to traditional volume-of-fluid (VOF) methods, numerical challenges remain particularly for solutions involving capillary waves and in the higher Weber number, low surface tension range. A depth-averaged approximation of the Navier–Stokes equations has been used to explore the effect of surface tension, grid resolution, and inertia on thin-film rimming solution accuracy and numerical stability. In shock and pooling solutions where capillary ripples are present, solution stability, and accuracy are shown to be highly sensitive to surface tension. The common practice in analytical studies of enforcing unphysical low Weber number stability constraints is shown to stabilize the solution by artificially damping capillary oscillations. This approach, however, although providing stable solutions is shown to adversely affect solution accuracy. An alternative grid resolution-based stability criterion is demonstrated and used to obtain numerically stable shock and pooling solutions without recourse to unphysical surface tension values. This allows for the accurate simulation of thin-film flows with capillary waves within the constrained parameter space corresponding to physical material and flow properties. Results obtained using the proposed formulation and solution strategy show good agreement with available experimental data from literature for low Re coating flows and moderate to high Re falling wavy film flows.

FIGURES IN THIS ARTICLE
<>
Copyright © 2016 by ASME
Your Session has timed out. Please sign back in to continue.

References

O'Brien, S. B. G. , and Schwartz, L. W. , 2002, “ Theory and Modeling of Thin Film Flows,” Enc. Surf. Colloid Sci., 4(1), pp. 5283–5297.
Glahn, A. , and Wittig, S. , 1996, “ Two Phase Air/Oil Flow in Aero Engine Bearing Chambers: Characterization of Oil Film Flows,” ASME J. Eng. Gas Turbines Power, 118(3), pp. 578–583. [CrossRef]
Chandra, B. , Simmons, K. , Pickering, S. , Collicott, S. H. , and Wiedemann, N. , 2013, “ Study of Gas/Liquid Behavior Within an Aeroengine Bearing Chamber,” ASME J. Eng. Gas Turbines Power, 135(5), p. 051201. [CrossRef]
Hirt, C. , and Nichols, B. , 1981, “ Volume of Fluid (VOF) Method for the Dynamics of Free Boundaries,” J. Comput. Phys., 39(1), pp. 201–225. [CrossRef]
Gao, D. , Morley, N. , and Dhir, V. , 2003, “ Numerical Simulation of Wavy Falling Film Flow Using {VOF} Method,” J. Comput. Phys., 192(2), pp. 624–642. [CrossRef]
Miyara, A. , 2000, “ Numerical Simulation of Wavy Liquid Film Flowing Down on a Vertical Wall and an Inclined Wall,” Int. J. Therm. Sci., 39(9–11), pp. 1015–1027. [CrossRef]
Ashmore, J. , Hosoi, A. E. , and Stone, H. A. , 2003, “ The Effect of Surface Tension on Rimming Flows in a Partially Filled Rotating Cylinder,” J. Fluid Mech., 479, pp. 65–98. [CrossRef]
Benilov, E. S. , Lapin, V. N. , and O'Brien, S. B. G. , 2012, “ On Rimming Flows With Shocks,” J. Eng. Math., 75(1), pp. 49–62. [CrossRef]
Wang, C. , Morvan, H. P. , Hibberd, S. , Cliff, K. A. , Anderson, A. , and Jacobs, A. , 2012, “ Specifying and Benchmarking a Thin Film Model for Oil Systems Applications in ANSYS Fluent,” ASME Paper No. GT2012-68984.
Kay, E. D. , Hibberd, S. , and Power, H. , 2014, “ A Depth-Averaged Model for Non-Isothermal Thin-Film Rimming Flow,” Int. J. Heat Mass Transfer, 70, p. 10031015. [CrossRef]
Williams, J. , 2008, “ Thin Film Rimming Flows Subject to Droplet Impact at the Surface,” Ph.D. thesis, The University of Nottingham, Nottingham, UK.
Wang, C. , Morvan, H. P. , Hibberd, S. , and Cliffe, K. A. , 2011, “ Thin Film Modelling for Aero-Engine Bearing Chambers,” ASME Paper No. GT2011-46259.
Villegas-Diaz, M. , Power, H. , and Riley, D. S. , 2005, “ Analytical and Numerical Studies of the Stability of Thin-Film Rimming Flows Subject to Surface Shear,” J. Fluid Mech., 541, pp. 317–344. [CrossRef]
Taliby, S. , and Portalski, S. , 1961, “ The Optimum Concentration of Surface Active Agent for the Suppression of Ripples,” Trans. Inst. Chem. Eng., 39(a), pp. 328–336.
Tirumkudulu, M. , and Acrivos, A. , 2001, “ Coating Flows Within a Rotating Horizontal Cylinder: Lubrication Analysis, Numerical Computations, and Experimental Measurements,” Phys. Fluids, 13(1), pp. 14–19. [CrossRef]
Liu, J. , and Gollub, J. P. , 1994, “ Solitary Wave Dynamics of Film Flows,” Phys. Fluids, 6(5), pp. 1702–1712. [CrossRef]
MacCormack, R. , 2003, “ The Effect of Viscosity in Hypervelocity Impact Cratering,” J. Spacecr. Rockets, 40(5), pp. 757–763. [CrossRef]
Denner, F. , and van Wachem, B. G. , 2015, “ Numerical Time-Step Restrictions as a Result of Capillary Waves,” J. Comput. Phys., 285(1), pp. 24–40. [CrossRef]
Lamb, H. , 1994, Hydrodynamics, 6th ed., Cambridge University Press, Cambridge, UK. [PubMed] [PubMed]
Nosoko, T. , Yoshimura, P. , Nagata, T. , and Oyakawa, K. , 1996, “ Characteristics of Two-Dimensional Waves on a Falling Liquid Film,” Chem. Eng. Sci., 51(5), pp. 725–732. [CrossRef]

Figures

Grahic Jump Location
Fig. 1

Thin-film rimming flow geometry and coordinate system used

Grahic Jump Location
Fig. 2

Thin-film rimming flow solution classification into; (a) shear dominated smooth flow, (b) a transitional shock flow regime where shear and gravity are in balance, and (c) gravity dominated pool flow

Grahic Jump Location
Fig. 3

Thin-film flow dimensionless parameters

Grahic Jump Location
Fig. 4

(a) Film height profiles and (b) momentum source terms for smooth solutions in cases A1–A3

Grahic Jump Location
Fig. 5

(a) Film height profiles in cases B1–B3 and (b) momentum source terms for the stable shock solution in case B2

Grahic Jump Location
Fig. 6

(a) Film height profiles in cases C1 and C2 and (b) momentum source terms for the stable pool solutions in case C2

Grahic Jump Location
Fig. 7

Polar plot of film solution illustrating the artificial smoothing of pool solutions due to the use of very high surface tension coefficients

Grahic Jump Location
Fig. 8

Results of the grid sensitivity studies for (a) case B2 (Δx=1.15 mm), case B2f1 (Δx=0.58 mm), cases B2f2 (Δx=0.29 mm), case B2c (Δx=2.30 mm), and (b) case B1 (Δx=1.15 mm), case B1f3 (Δx=0.14 mm)

Grahic Jump Location
Fig. 9

Sensitivity of small-wavelength disturbances in stable cases B1 (σ=0.0245N/m) and B2 (σ=0.5N/m) to surface tension

Grahic Jump Location
Fig. 10

Effect of inertia representation on (a) smooth solution A1, (b) shock solution B2, and (c) pool solution C2

Grahic Jump Location
Fig. 11

Experimental measurements from Ref. [15] together with equivalent film thickness predictions for (a) case 1 and (b) case 2 as defined in Table 3

Grahic Jump Location
Fig. 12

ETFM predictions for solitary waves traveling down a 6.4 deg slope with Re = 29 and inlet forcing frequency of 1.5 Hz. ETFM model results for both simplified and full inertia representations are shown together with equivalent experimental measurements from Ref. [16].

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In