بررسی عددی شدت فرسایش حباب زایی جریان اطراف هیدروفویل NACA0015 براساس دیدگاه اویلری-لاگرانژی

نوع مقاله : مقاله پژوهشی

نویسندگان

1 دانشکده ی مهندسی مکانیک، دانشگاه تربیت دبیر شهیدرجایی

2 دانشکده ی مهندسی مکانیک، دانشگاه تربیت دبیر شهید رجایی

چکیده

در این تحقیق برای بررسی شدت فرسایش ناشی از حباب‌زایی (کاویتاسیون)،رفتار حباب اطراف هیدروفویل دوبعدی ۰۰۱۵N‌A‌C‌A با استفاده ازدیدگاه
اویلری ـ لاگرانژی شبیه‌سازی شده است. بررسی درشت‌مقیاس (ماکروسکوپیک) جریان حباب‌زایی توسط مدل مخلوط همگن (روش اویلری) و مسیر حرکت حباب بر اساس نیروهای وارده با استفاده از قانون دوم نیوتن و توسعه کد عددی مشخص شد (روش لاگرانژی). برای ادغام نتایج به دست آمده و حل آنها از روش رانچ‌کوتای مرتبه ۴ با گام زمانی متغیر استفاده شده است، که سرعت حل داده‌ها را تا ۱۰ برابر افزایش داده است. نتایج به دست آمده نشان داد «شروع حباب‌زایی» برای جریان با هسته‌های بزرگ‌تر، بهتر قابل مشاهده است. با کاهش عدد حباب‌زایی نرخ رشد حباب بیشتر شده و با افزایش شعاع حباب، شدت فرسایش افزایش می‌یابد. مکان احتمالی فرسایش در محدوده‌ی انتهایی کویتی در سطح هیدروفویل قرار دارد. نتایج به دست آمده با دیگر کارهای منتشر شده مقایسه شد و از دقت قابل قبولی برخوردار بود.

کلیدواژه‌ها


عنوان مقاله [English]

Numerical Evaluation of Cavitation Erosion Intensity around NACA0015 Hydrofoil Based on Eulerian-Lagrangian Approach

نویسندگان [English]

  • A. katoozi 1
  • M. Mahdi 2
1 D‌e‌p‌t. o‌f M‌e‌c‌h‌a‌n‌i‌c‌a‌l E‌n‌g‌i‌n‌e‌e‌r‌i‌n‌g S‌h‌a‌h‌i‌d R‌a‌j‌a‌e‌e T‌e‌a‌c‌h‌e‌r T‌r‌a‌i‌n‌i‌n‌g U‌n‌i‌v‌e‌r‌s‌i‌t‌y
2 D‌e‌p‌t. o‌f M‌e‌c‌h‌a‌n‌i‌c‌a‌l E‌n‌g‌i‌n‌e‌e‌r‌i‌n‌g S‌h‌a‌h‌i‌d R‌a‌j‌a‌e‌e T‌e‌a‌c‌h‌e‌r T‌r‌a‌i‌n‌i‌n‌g U‌n‌i‌v‌e‌r‌s‌i‌t‌y
چکیده [English]

In this study, to investigate the intensity of cavitation-induced erosion, the bubble behavior around the NACA0015 2D hydrofoil was simulated from the Eulerian-Lagrangian perspective. Macroscopic examination of the cavitation flow was determined by a homogeneous mixture model (Eulerian method) and the trajectory of bubble motion based on the applied forces using Newton's second law and the development of numerical code (Lagrange method). One way to reduce the computational cost of the Lagrangian perspective is to use the Discrete Phase Model (DPM). In this method, the fluid is considered as a continuous environment, while the discrete phase is solved by tracing a large number of particles in the calculated flow field. The behavior of the bubble arises from the pressure gradient caused by the flow. Bubble oscillations were obtained from the modified Rayleigh-Plesset-Keller-Herring equation. This equation considers the compressible behavior of the bubble as the bubble collapse velocity approaches the speed of sound as well as the slip velocity between the bubble and the moving liquid. To pair the obtained results and solve them, the fourth-order Runge-Kutta method with variable time step was used, which increased the data solution speed up to 10 times. From the Keller& Kolodner relationship, a pressure wave emitted from the collapse of a spherical bubble, and the model of Soyama et al., the total energy of the cavitation-induced shocks, which is the result of the accumulation of all the shocks on each other, is obtained. The actual effects of flow for cavitation inception and erosion were investigated. Different cavitation numbers were used for cavitation inception with different radii. The results showed that the nucleation process occured in the cavitation inception numbers and the cavitation inception for flow with larger nuclei was visible better. As the cavitation number decreases, the bubble growth rate increases and as the bubble radius increases, the erosion intensity increases. At high cavitation numbers, the bubble oscillates around its initial radius; however, at the lowest cavitation number in this article, the number , we see an increase of nearly times the radius compared to the original radius. The erosion power of bubbles with an initial radius of is approximately times that of the erosion power of bubbles with an initial radius of and about times that of the initial bubbles of . The probable site of erosion is at the end of cavity at the hydrofoil level. As the bubbles increase in size, the number of collapses and their strength increase, and the dispersion of the distribution at the hydrofoil surface increases. The results were compared with other published works and had acceptable accuracy.

کلیدواژه‌ها [English]

  • Cavitation flow
  • Eulerian-Lagrangian method
  • Erosion Intensity
  • Bubble dynamics
  • Numerical prediction