Skip to main content
Log in

Binary collision of CMAS droplets—Part II: Unequal-sized droplets

  • Published:
Journal of Materials Research Aims and scope Submit manuscript

Abstract

The analysis presented in Part I of this study on the binary collision of equal molten calcium–magnesium–alumino–silicate (CMAS) droplets is extended to investigate the flow and interfacial dynamics of unequal CMAS droplet collision. Numerical investigations of head-on, off-center, and grazing collisions of two CMAS droplets of size 1 and 2 mm are conducted at pressure and temperature of 20 atm and 1548 K, respectively, that are representative of a gas-turbine combustor. At these conditions, the physical properties of CMAS are density, ρCMAS = 2690 kg/m3, surface tension between CMAS/air, σCMAS = 0.40 N/m, and viscosity, μCMAS = 11.0 N-s/m2. The primary difference between the CMAS and a fictitious fluid with viscosity 1/10 of CMAS was higher deformation for the lower viscosity case, leading to stretching and subsequent breakup of the liquid structure. These mechanisms are supported by the time evolution of surface, kinetic, and viscous dissipation energies.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Figure 1:
Figure 2:
Figure 3:
Figure 4:
Figure 5:
Figure 6:
Figure 7:
Figure 8:
Figure 9:
Figure 10:

Similar content being viewed by others

References

  1. H. Ganti, P. Khare, and L. Bravo: Binary collision of CMAS droplets—Part I: Equal sized droplets. J. Mater. Res. (Focus Issue - Sandphobic Thermal/Environmental Barrier Coatings for Gas Turbine Engines), 1–15 (2020). doi:10.1557/jmr.2020.138.

    Google Scholar 

  2. A. Ghoshal, M. Murugan, M.J. Walock, A. Nieto, B.D. Barnett, M.S. Pepi, J.J. Swab, D. Zhu, K.A. Kerner, C.R. Rowe, C.Y. Shiao, D.A. Hopkins, and G. A. Gazonas: Molten particulate impact on tailored thermal barrier coatings for gas turbine engine. J. Eng. Gas Turb. Power 140 (2018).

  3. M. Murugan, A. Ghoshal, A. Nieto, M. Walock, L. Bravo, N. Jain, M. Pepi, J. Swab, D. Zhu, R.T. Pegg, C. Rowe, A. Flatau, and K. Kerner: Prevention of molten sand attack on thermal barrier coatings for rotorcraft gas turbine blades—A round robin test evaluation. In Annual Forum Proceedings (AHS International, Fairfax, VA, USA, 2018).

  4. N. Libertowski, N. Plewacki, and J.P. Bons: The effect of temperature and melting relative to particle deposition in gas turbines. In AIAA Scitech 2019 Forum, Jan 7 to 11, San Diego, CA, USA (American Institute of Aeronautics and Astronautics, Reston, VA, USA, 2019).

  5. L.G. Bravo, Q. Xue, M. Murugan, A. Ghoshal, M. Walock, and A. Flatau: Particle transport analysis of sand ingestion in gas turbine jet engines. In 53rd AIAA/SAE/ASEE Joint Propulsion Conference, Jul 10 to Jul 12, Atlanta, GA, USA (AIAA, Reston, VA, USA, 2017); pp. 1–14.

  6. N. Jain, L. Bravo, S. Bose, D. Kim, M. Murugan, A. Ghoshal, and A. Flatau: Turbulent multiphase flow and particle deposition of sand ingestion for high-temperature turbine blades. In Proceedings of the Summer Program, Jun 24 to Jul 20, Stanford, CA, USA (Center for Turbulence Research, Stanford, CA, USA, 2018); pp. 35–44.

  7. I. Spitsberg, and J. Steibel: Thermal and environmental barrier coatings for SiC/SiC CMCs in aircraft engine applications*. Int. J. Appl. Ceram. Technol. 1(4), 291–301 (2004).

    Article  CAS  Google Scholar 

  8. K.N. Lee, D.S. Fox, J.I. Eldridge, D. Zhu, R.C. Robinson, N.P. Bansal, and R.A. Miller: Upper temperature limit of environmental barrier coatings based on mullite and BSAS. J. Am. Ceram. Soc. 86(8), 1299–1306 (2003).

    Article  CAS  Google Scholar 

  9. K.N. Lee: Current status of environmental barrier coatings for Si-based ceramics. Surf. Coat. Technol. 133–134, 1–7 (2000).

    Google Scholar 

  10. M. Murugan, A. Ghoshal, M. Walock, A. Nieto, L. Bravo, B. Barnett, M. Pepi, J. Swab, R.T. Pegg, C. Rowe, D. Zhu, and K. Kerner: Microstructure based material-sand particulate interactions and assessment of coatings for high temperature turbine blades. In Proceedings of ASME Turbo Expo (IGTI), Charlotte, NC USA (ASME, New York, NY, USA, 2017).

  11. D.K. Larry Fehrenbacher, J. Kutsch, I. Vesnovsky, E. Fehrenbacher, A. Ghoshal, M. Walock, M. Murugan, and A. Nieto: Advanced environmental barrier coatings for SiC CMCs. In Advances in Ceramics for Environmental, Functional, Structural, and Energy Applications II, Vol. 266, edited by M.M. Mahmoud, K. Sridharan, H. Colorado, A.S. Bhalla, J.P. Singh, S. Gupta, J. Langhorn, A. Jitianu, and N. Jose Manjooran (John Wiley & Sons, Hoboken, NJ, USA, 2019); pp 83–93.

    Google Scholar 

  12. W. Song, S. Yang, M. Fukumoto, Y. Lavallée, S. Lokachari, H. Guo, Y. You, and D.B. Dingwell: Impact interaction of in-flight high-energy molten volcanic ash droplets with jet engines. Acta Mater. 171, 119–131 (2019).

    Article  CAS  Google Scholar 

  13. S. Singh and D. Tafti: Particle deposition model for particulate flows at high temperatures in gas turbine components. Int. J. Heat Fluid Flow 52, 72–83 (2015).

    Article  Google Scholar 

  14. D. Pearson and R. Brooker: The accumulation of molten volcanic ash in jet engines, simulating the role of magma composition, ash particle size and thermal barrier coatings. J. Volcanol. Geotherm. Res. 389, 106707 (2020).

    Article  CAS  Google Scholar 

  15. L. Bravo, N. Jain, P. Khare, M. Murugan, A. Ghoshal, and A. Flatau: Physical aspects of particle dynamics and deposition in turboshaft engines. J. Mater. Res. (under review), (2020).

    Google Scholar 

  16. D. Liu, P. Zhang, C.K. Law, and Y. Guo: Collision dynamics and mixing of unequal-size droplets. Int. J. Heat Mass Transfer 57, 421–428 (2013).

    Article  Google Scholar 

  17. C. Tang, J. Zhao, P. Zhang, C.K. Law, and Z. Huang: Dynamics of internal jets in the merging of two droplets of unequal sizes. J. Fluid Mech. 795, 671–689 (2016).

    Article  CAS  Google Scholar 

  18. C. Tang, P. Zhang, and C.K. Law: Bouncing, coalescence, and separation in head-on collision of unequal-size droplets. Phys. Fluids 24, 22101–1–22101–15 (2012).

    Article  Google Scholar 

  19. K.-L. Pan, C.K. Law, and B. Zhou: Experimental and mechanistic description of merging and bouncing in head-on binary droplet collision. J. Appl. Phys. 103, 064901 (2008).

    Article  Google Scholar 

  20. N. Ashgriz and J.Y. Poo: Coalescence and separation in binary collisions of liquid drops. J. Fluid Mech. 221, 183–204 (1990).

    Article  CAS  Google Scholar 

  21. N. Nikolopoulos and G. Bergeles: The effect of gas and liquid properties and droplet size ratio on the central collision between two unequal-size droplets in the reflexive regime. Int. J. Heat Mass Transfer 54, 678–691 (2011).

    Article  CAS  Google Scholar 

  22. M. Yoshino, J. Sawada, and K. Suzuki: Numerical simulation of head-on collision dynamics of binary droplets with various diameter ratios by the two-phase lattice kinetic scheme. Comput. Fluids 168, 304–317 (2018).

    Article  Google Scholar 

  23. V. Notaro, P. Khare, and J.G. Lee: Mixing characteristics of non-Newtonian impinging jets at elevated pressures. Flow Turbul. Combust. 102, 355–372 (2019).

    Article  CAS  Google Scholar 

  24. P. Khare and V. Yang: Newtonian and non-newtonian liquid droplet breakup: Child droplet statistics. Proceedings of ICLASS 2015, 13th Triennial International Conference on Liquid Atomization and Spray Systems, Aug 23 to Aug 27, Tainan, Taiwan (2015).

  25. P. Khare and V. Yang: Breakup of non-newtonian liquid droplets. In 44th AIAA Fluid Dynamics Conference, Jan 16 to Jan 20, Atlanta, GA, USA (AIAA, Reston, VA, USA, 2014).

  26. J. Gamertsfelder, P. Khare, and L.G. Bravo: Investigation of atomization behavour of liquid monopropellants in pintle injectors. In ASME Turbo Expo, June 22–26, London, England (ASME, New York, NY, USA, 2020).

  27. X. Chen, D. Ma, P. Khare, and V. Yang: Energy and mass transfer during binary droplet collision. In 49th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, Jan 4 to Jan 7, Orlando, FL, USA (AIAA, Reston, VA, USA, 2011); pp. 1–14..

  28. G. Tryggvason, R. Scardovelli, and S. Zaleski: Direct Numerical Simulations of Gas–Liquid Multiphase Flows (Cambridge University Press, New York, NY, USA, 2011).

    Google Scholar 

  29. J.U. Brackbill, D.B. Kothe, and C. Zemach: A continuum method for modeling surface tension. J. Comput. Phys. 100, 335–354 (1992).

    Article  CAS  Google Scholar 

  30. O. Desjardins and V. Moureau: Methods for multiphase flows with high density ratio. In Proceedings of the Summer Program, Jun 27 to Jul 23, Stanford, CA, USA (Center for Turbulence Research, Stanford, CA, USA, 2010); pp. 313–322.

  31. M. Gorokhovski and M. Herrmann: Modeling primary atomization. Annu. Rev. Fluid Mech. 40, 343–366 (2008).

    Article  Google Scholar 

  32. S. Popinet: An accurate adaptive solver for surface-tension-driven interfacial flows. J. Comput. Phys. 228, 5838–5866 (2009).

    Article  CAS  Google Scholar 

  33. D. Gerlach, G. Tomar, G. Biswas, and F. Durst: Comparison of volume-of-fluid methods for surface tension-dominant two-phase flows. Int. J. Heat Mass Transfer 49, 740–754 (2006).

    Article  Google Scholar 

  34. J. Qian and C.K. Law: Regimes of coalescence and separation in droplet collision. J. Fluid Mech. 331, 59–80 (1997).

    Article  CAS  Google Scholar 

  35. H. Ganti, P. Khare, and L. Bravo: Binary collision of CMAS droplets-Part I: Equal sized droplets. J. Mater. Res. (Focus Issue–Sandphobic Thermal/Environmental Barrier Coatings for Gas Turbine Engines), 1–15 (2020). doi:10.1557/jmr.2020.138

    Google Scholar 

Download references

Acknowledgments

The authors gratefully acknowledge the Advanced Computing Center at the University of Cincinnati, and the DoD High Performance Computing Modernization Program (HPCMP) subproject number, ARLAP44862H70 for the computational resources. Luis Bravo was supported by the US Army Research Laboratory 6.1 basic research program in propulsion sciences.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Prashant Khare.

Appendix A: Computational Challenges, Model Validation, and Grid-Sensitivity Analysis

Appendix A: Computational Challenges, Model Validation, and Grid-Sensitivity Analysis

Irrespective of the numerical method, the challenges accompanying the numerical simulation of incompressible two-phase systems increase dramatically, as the density ratio increases [1, 2 ]. The time integration scheme used in the current approach involves a classical time-splitting projection method, which requires the solution of the Poisson equation to obtain the pressure field:

$$\nabla \cdot \left[{\displaystyle{{\Delta t} \over {\rho_{n + {1 \over 2}}}}\nabla p_{n + {1 \over 2}}} \right] = \nabla \cdot u_\ast.$$
((15))

Equation (A1) is solved using a standard multigrid V-cycle methodology, and for large density and viscosity ratios, its solution suffers from slow convergence rates. One of the ways to overcome this issue is by using high grid resolution to resolve the steep density and viscosity gradients at the interface to ensure consistency in the momentum equation. Another method of speeding up the convergence rate is to spatially filter the interface during reconstruction. Even though the current methodology performs very well for the current configuration of droplet interaction at high viscosity and density ratios, the convergence can seriously degrade, depending on the problem and interface topology [3 ], in comparison with other methods [4]. Therefore, for all the cases conducted as a part of this research effort, including the validation study described in the next section, we have used both the aforementioned strategies to ensure accuracy: high grid resolution and spatially filtering (at least once) to ensure numerical accuracy and adequate resolution of the gas–liquid interface.

As a first step, grid-sensitivity analysis which is conducted to ensure appropriate grid resolution is used to resolve the physics under consideration. The canonical configuration of equal CMAS droplets colliding head-on (B = 0.0) is selected for the grid-sensitivity study. Figure A1 shows a comparison of the liquid morphology for four different refinement levels described below

  1. (i)

    level 6 at liquid/gas interface, level 5 for the droplet interior, and level 3 for the rest of the domain—L6

  2. (ii)

    level 7 at liquid/gas interface, level 6 for the droplet interior, and level 3 for the rest of the domain—L7

  3. (iii)

    level 8 at liquid/gas interface, level 7 for the droplet interior, and level 4 for the rest of the domain—L8

  4. (iv)

    level 9 at liquid/gas interface, level 8 for the droplet interior, and level 5 for the rest of the domain—L9

Figure A1:
figure 11

Grid-sensitivity study based on levels 6, 7, 8, and 9. The figure on the right shows an overlay of the liquid interface for these levels.

As seen clearly, L6 is unable to refine the interface sufficiently. While L7 resolves the interface better, to ensure that the gas film when the droplets come closer to each other is resolved, L8 and L9 were investigated, which show almost identical results for interface deformation and evolution as well as the gas film. Therefore, for this study, L8 was selected as the grid resolution. While for high Weber number droplet collision, such as the current study, bouncing is not expected to be an outcome, in addition to gradient and value-based refinements, distance-based refinement is employed to ensure that the gas film is resolved accurately.

Next, model validation is conducted by simulating the experiments conducted by Qian and Law [5 ]. It should be noted that since no study in the past has investigated CMAS droplet collision, we have to resort to using data on tetradecane droplet collision for validation purposes. We chose a case that incorporates merging, retracting, and formation of satellite droplets to ensure that our framework can accurately model different aspects of the droplet collision phenomena. In this experiment, conducted at a pressure of 1 atm and 300 K, two droplets of diameter 336 μm with 2.48 m/s collide. The density and viscosity of air at these conditions are 1.18 kg/m3 and 1.79 × 10−4 N-s/m2, respectively. The density, viscosity, and surface tension of tetradecane are 785.88 kg/m3, 2.21× 10−3 N-s/m2, and 0.02656 N/m. Figure A2 shows the time evolution of events that take place, as these droplets collide. The left side shows the experimental measurements and the right side shows results from the current simulation, showing excellent comparison. All flow features, including droplet coalescence, ligament formation and elongation, subsequent separation by pinching, and satellite droplet formation, are accurately captured.

Figure A2:
figure 12

Time evolution of liquid interface when two tetradecane droplets collide at an impact factor, B = 0.06. (a) Experimental images of Qian and Law [5] shown on the left and (b) results from current simulations.

Details of grid-sensitivity and experimental validation can be found in the previous article for Binary Collision of CMAS Droplet–Part I: Equal Sized Droplets [6 ].

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ganti, H., Khare, P. & Bravo, L. Binary collision of CMAS droplets—Part II: Unequal-sized droplets. Journal of Materials Research 35, 2275–2287 (2020). https://doi.org/10.1557/jmr.2020.153

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1557/jmr.2020.153

Keywords

Navigation