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.
Similar content being viewed by others
References
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.
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).
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).
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).
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.
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.
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).
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).
K.N. Lee: Current status of environmental barrier coatings for Si-based ceramics. Surf. Coat. Technol. 133–134, 1–7 (2000).
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).
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.
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).
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).
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).
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).
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).
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).
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).
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).
N. Ashgriz and J.Y. Poo: Coalescence and separation in binary collisions of liquid drops. J. Fluid Mech. 221, 183–204 (1990).
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).
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).
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).
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).
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).
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).
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..
G. Tryggvason, R. Scardovelli, and S. Zaleski: Direct Numerical Simulations of Gas–Liquid Multiphase Flows (Cambridge University Press, New York, NY, USA, 2011).
J.U. Brackbill, D.B. Kothe, and C. Zemach: A continuum method for modeling surface tension. J. Comput. Phys. 100, 335–354 (1992).
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.
M. Gorokhovski and M. Herrmann: Modeling primary atomization. Annu. Rev. Fluid Mech. 40, 343–366 (2008).
S. Popinet: An accurate adaptive solver for surface-tension-driven interfacial flows. J. Comput. Phys. 228, 5838–5866 (2009).
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).
J. Qian and C.K. Law: Regimes of coalescence and separation in droplet collision. J. Fluid Mech. 331, 59–80 (1997).
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
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
Corresponding author
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:
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
-
(i)
level 6 at liquid/gas interface, level 5 for the droplet interior, and level 3 for the rest of the domain—L6
-
(ii)
level 7 at liquid/gas interface, level 6 for the droplet interior, and level 3 for the rest of the domain—L7
-
(iii)
level 8 at liquid/gas interface, level 7 for the droplet interior, and level 4 for the rest of the domain—L8
-
(iv)
level 9 at liquid/gas interface, level 8 for the droplet interior, and level 5 for the rest of the domain—L9
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.
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
About this article
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
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1557/jmr.2020.153