scholarly journals Comparison of multiphase models for computing shock-induced bubble collapse

2019 ◽  
Vol 30 (8) ◽  
pp. 3845-3877 ◽  
Author(s):  
Eric Goncalves Da Silva ◽  
Philippe Parnaudeau

Purpose The purpose of this paper is to quantify the relative importance of the multiphase model for the simulation of a gas bubble impacted by a normal shock wave in water. Both the free-field case and the collapse near a wall are investigated. Simulations are performed on both two- and three-dimensional configurations. The main phenomena involved in the bubble collapse are illustrated. A focus on the maximum pressure reached during the collapse is proposed. Design/methodology/approach Simulations are performed using an inviscid compressible homogeneous solver based on different systems of equations. It consists in solving different mixture or phasic conservation laws and a transport-equation for the gas volume fraction. Three-dimensional configurations are considered for which an efficient massively parallel strategy was developed. The code is based on a finite volume discretization for which numerical fluxes are computed with a Harten, Lax, Van Leer, Contact (HLLC) scheme. Findings The comparison of three multiphase models is proposed. It is shown that a simple four-equation model is well-suited to simulate such strong shock-bubble interaction. The three-dimensional collapse near a wall is investigated. It is shown that the intensity of pressure peaks on the wall is drastically increased (more than 200 per cent) in comparison with the cylindrical case. Research limitations/implications The study of bubble collapse is a key point to understand the physical mechanism involved in cavitation erosion. The bubble collapse close to the wall has been addressed as the fundamental mechanism producing damage. Its general behavior is characterized by the formation of a water jet that penetrates through the bubble and the generation of a blast wave during the induced collapse. Both the jet and the blast wave are possible damaging mechanisms. However, the high-speed dynamics, the small spatio-temporal scales and the complicated physics involved in these processes make any theoretical and experimental approach a challenge. Practical implications Cavitation erosion is a major problem for hydraulic and marine applications. It is a limiting point for the conception and design of such components. Originality/value Such a comparison of multiphase models in the case of a strong shock-induced bubble collapse is clearly original. Usually models are tested separately leading to a large dispersion of results. Moreover, simulations of a three-dimensional bubble collapse are scarce in the literature using such fine grids.

2021 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Xin Wang ◽  
Ting-Qiang Xie

Purpose Cavitation erosion has always been a common technical problem in a hydraulic discharging structure. This paper aims to investigate the cavitation erosion behavior of hydraulic concrete under high-speed flow. Design/methodology/approach A high-speed and high-pressure venturi cavitation erosion generator was used to simulate the strong cavitation. The characteristics of hydrodynamic loads of cavitation bubble collapse zone, the failure characteristics and the erosion development process of concrete were investigated. The main influencing factors of cavitation erosion were discussed. Findings The collapse of the cavitation bubble group produced a high frequency, continuous and unsteady pulse load on the wall of concrete, which was more likely to cause fatigue failure of concrete materials. The cavitation action position and the main frequency of impact load were greatly affected by the downstream pressure. A power exponential relationship between cavitation load, cavitation erosion and flow speed was observed. With the increase of concrete strength, the degree of damage of cavitation erosion was approximately linearly reduced. Originality/value After cavitation erosion, a skeleton structure was formed by the accumulation of granular particles, and the relatively independent bulk structure of the surface differed from the flake structure formed after abrasion.


Author(s):  
Alireza Rahimi ◽  
Aravindhan Surendar ◽  
Aygul Z. Ibatova ◽  
Abbas Kasaeipoor ◽  
Emad Hasani Malekshah

Purpose This paper aims to investigate the three-dimensional natural convection and entropy generation in the rectangular cuboid cavities included by chamfered triangular partition made by polypropylene. Design/methodology/approach The enclosure is filled by multi-walled carbon nanotubes (MWCNTs)-H2O nanofluid and air as two immiscible fluids. The finite volume approach is used for computation. The fluid flow and heat transfer are considered with combination of local entropy generation due to fluid friction and heat transfer. Moreover, a numerical method is developed based on three-dimensional solution of Navier–Stokes equations. Findings Effects of side ratio of triangular partitions (SR = 0.5, 1 and 2), Rayleigh number (103 < Ra < 105) and solid volume fraction (f = 0.002, 0.004 and 0.01 Vol.%) of nanofluid are investigated on both natural convection characteristic and volumetric entropy generation. The results show that the partitions can be a suitable method to control fluid flow and energy consumption, and three-dimensional solutions renders more accurate results. Originality/value The originality of this work is to study the three-dimensional natural convection and entropy generation of a stratified system.


Author(s):  
Alessandra Dal Cengio Leonardi ◽  
Nickolas Keane ◽  
Cynthia Bir ◽  
Pamela VandeVord

Studies on blast neurotrauma have focused on investigating the effects of exposure to free-field blast representing the simplest form of blast threat scenario without considering any reflecting surfaces. However, in reality personnel are often located within enclosures or nearby reflecting walls causing a complex blast environment, that is, involving shock reflections and/or compound waves from different directions. In fact, when a blast wave interacts with nearby structures, reflected shock waves are generated and complex three-dimensional shock waves are formed. Complex shock wave overpressure-time traces are significantly different from free-field profiles because reflections can cause super-positioning of shock waves resulting in increased pressure magnitudes and multiple pressure peaks. Very importantly, the shocks arrive from different directions which would invoke a different biomechanical response than a one-dimensional exposure. It has been reported that in complex wave environments, the extent of the injuries becomes a function of the location related to the surrounding structures rather than a function of the distance from the center of the explosion, as it is for free-field conditions (Yelverton et al. 1993; Mayorga 1997; Stuhmiller 1997). Furthermore, the resulting injuries when the individual is in confined spaces are noted to be more severe (Yelverton et al. 1993; Leibovici et al. 1996). The purpose of this study was to design a complex wave testing system and perform a preliminary investigation of the intracranial pressure (ICP) response of rats exposed to a complex blast wave environment. Furthermore, we explored the effects of head orientation in the same environment.


2020 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Soraya Torkaman ◽  
Ghasem Barid Loghmani ◽  
Mohammad Heydari ◽  
Abdul-Majid Wazwaz

Purpose The purpose of this paper is to investigate a three-dimensional boundary layer flow with considering heat and mass transfer on a nonlinearly stretching sheet by using a novel operational-matrix-based method. Design/methodology/approach The partial differential equations that governing the problem are converted into the system of nonlinear ordinary differential equations (ODEs) with considering suitable similarity transformations. A direct numerical method based on the operational matrices of integration and product for the linear barycentric rational basic functions is used to solve the nonlinear system of ODEs. Findings Graphical and tabular results are provided to illustrate the effect of various parameters involved in the problem on the velocity profiles, temperature distribution, nanoparticle volume fraction, Nusselt and Sherwood number and skin friction coefficient. Comparison between the obtained results, numerical results based on the Maple's dsolve (type = numeric) command and previous existing results affirms the efficiency and accuracy of the proposed method. Originality/value The motivation of the present study is to provide an effective computational method based on the operational matrices of the barycentric cardinal functions for solving the problem of three-dimensional nanofluid flow with heat and mass transfer. The convergence analysis of the presented scheme is discussed. The benefit of the proposed method (PM) is that, without using any collocation points, the governing equations are converted to the system of algebraic equations.


2013 ◽  
Vol 655-657 ◽  
pp. 449-456
Author(s):  
Hong Ming Zhang ◽  
Li Xiang Zhang

The paper presents numerical prediction of cavitation erosion on a Francis turbine runner using CFD code. The SST turbulence model is employed in the Reynolds averaged Navier–Stokes equations in this study. A mixture assumption and a finite rate mass transfer model were introduced. The computing domain is discretized with a full three-dimensional mesh system of unstructured tetrahedral shapes. The finite volume method is used to solve the governing equations of the mixture model and the pressure-velocity coupling is handled via a Pressure Implicit with Splitting of Operators(PISO) procedure. Comparison the numerical prediction results with a real runner with cavitation damage, the region of higher volume fraction by simulation with the region of runner cavitation damage is consistent.


Author(s):  
Junaid Ahmad Khan ◽  
M. Mustafa ◽  
T. Hayat ◽  
Mustafa Turkyilmazoglu ◽  
A. Alsaedi

Purpose The purpose of the present study is to explore a three-dimensional rotating flow of water-based nanofluids caused by an infinite rotating disk. Design/methodology/approach Mathematical formulation is performed using the well-known Buongiorno model which accounts for the combined influence of Brownian motion and thermophoresis. The recently suggested condition of passively controlled wall nanoparticle volume fraction has been adopted. Findings The results reveal that temperature decreases with an increase in thermophoresis parameter, whereas it is negligibly affected with a variation in the Brownian motion parameter. Axial velocity is negative because of the downward flow in the vertical direction. Originality/value Two- and three-dimensional streamlines are also sketched and discussed. The computations are found to be in very good agreement with the those of existing studies in the literature for pure fluid.


Author(s):  
Cesar Martin Venier ◽  
Andrés Reyes Urrutia ◽  
Juan Pablo Capossio ◽  
Jan Baeyens ◽  
Germán Mazza

Purpose The purpose of this study is to assess the performance of ANSYS Fluent® and OpenFOAM®, at their current state of development, to study the relevant bubbling fluidized bed (BFB) characteristics with Geldart A, B and D particles. Design/methodology/approach For typical Geldart B and D particles, both a three-dimensional cylindrical and a pseudo-two-dimensional arrangement were used to measure the bed pressure drop and solids volume fraction, the latter by digital image analysis techniques. For a typical Geldart A particle, specifically to examine bubbling and slugging phenomena, a 2 m high three-dimensional cylindrical arrangement of small internal diameter was used. The hydrodynamics of the experimentally investigated BFB cases were also simulated for identical geometries and operating conditions using OpenFOAM® v6.0 and ANSYS Fluent® v19.2 at identical mesh and numerical setups. Findings The comparison between experimental and simulated results showed that both ANSYS Fluent® and OpenFOAM® provide a fair qualitative prediction of the bubble sizes and solids fraction for freely-bubbling Geldart B and D particles. For Geldart A particles, operated in a slugging mode, the qualitative predictions are again quite fair, but numerical values of relevant slug characteristics (length, velocity and frequency) slightly favor the use of OpenFOAM®, despite some deviations of predicted slug velocities. Originality/value A useful comparison of computational fluid dynamics (CFD) software performance for different fluidized regimes is presented. The results are discussed and recommendations are formulated for the selection of the CFD software and models involved.


Author(s):  
Seyed Ali Atyabi ◽  
Ebrahim Afshari ◽  
Mohammad Yaghoub Abdollahzadeh Jamalabadi

Purpose In this paper, a single module of cross-flow membrane humidifier is evaluated as a three-dimensional multiphase model. The purpose of this paper is to analyze the effect of volume flow rate, dry temperature, dew point wet temperature and porosity of gas diffusion layer on the humidifier performance. Design/methodology/approach In this study, one set of coupled equations are continuity, momentum, species and energy conservation is considered. The numerical code is benchmarked by the comparison of numerical results with experimental data of Hwang et al. Findings The results reveal that the transfer rate of water vapor and dew point approach temperature (DPAT) increase by increasing the volume flow rate. Also, it is found that the water recovery ratio (WRR) and relative humidity (RH) decrease with increasing volume flow rate. In addition, all mixed results decrease with increasing dry side temperature especially at high volume flow rates and this trend in high volume flow rates is more sensible. Although the transfer rate of water vapor and DPAT increases with increasing the wet inlet temperature, WRR and RH reduce. Increasing dew point temperature effect is more sensible at the wet side is compared with the dry side. The humidification performance will be enhanced with increasing diffusion layer porosity by increasing the wet inlet dew point temperature, but has no meaningful effect on other operating parameters. The pressure drop along humidifier gas channels increases with rising flow rate, consequently, the required power of membrane humidifier will enhance. Originality/value According to previous studies, the three-dimensional numerical multiphase model of cross-flow membrane humidifier has not been developed.


2019 ◽  
Vol 29 (7) ◽  
pp. 2434-2449 ◽  
Author(s):  
Xiangda Cui ◽  
Ahmed Bakkar ◽  
Wagdi George Habashi

Purpose This paper aims to introduce a three-dimensional smoothed particle hydrodynamics (SPH) framework for simulating supercooled large droplets (SLD) dynamics at aeronautical speeds. Design/methodology/approach To include the effects of the surrounding air, a multiphase model capable of handling high density-ratio problems is adopted. A diffusive term is incorporated to smooth the density field and avoid numerical instabilities. Additionally, a particle shifting technique is used to eliminate anisotropic particle distributions. Findings The framework is validated against low-speed droplet impingement experimental results and then applied to the droplet impingement at high speeds typical of SLD conditions. Preliminary parametric studies are conducted to investigate the post-impact splashing. It is observed that a thicker water film can decrease the crown diameter and a smaller impact angle can suppress upward and forward splashing. Originality/value A three-dimensional multiphase SPH framework for SLD dynamics at a wide range of impact speed is developed and validated. The effects of particle resolution, water film thickness and impact angle on the post-impact crown evolution are investigated.


2019 ◽  
Vol 141 (6) ◽  
Author(s):  
Masoud Khojasteh-Manesh ◽  
Miralam Mahdi

In the present study, cavitation erosion is investigated by implementing an Eulerian–Lagrangian approach. Three-dimensional two-phase flow is simulated in a microscale nozzle using Reynolds-averaged Navier–Stokes (RANS) solver along with realizable k−ε turbulence model and Schnerr–Sauer cavitation model. The numerical results are in agreement with experimental observations. A modified form of Rayleigh–Plesset–Keller–Herring equation along with bubble motion equation is utilized to simulate bubble dynamics. Average values of mixture properties over bubble surface are used instead of bubble-center values in order to account for nonuniformities around the bubble. A one-way coupling method is used between Lagrangian analysis and RANS solution. The impact pressure resulted from bubble collapse is calculated for evaluation of erosion in diesel and soy methyl ester (SME) biodiesel in different situations. The results show that the initial size of the bubbles is an important factor for determining the intensity of erosion. So, the bubbles erosive power increases when their initial radius increases. It is also found that the intensity of erosion in diesel is much higher than that of biodiesel and this is because of the differences in fuels properties, especially in viscosity and vapor pressure. The effect of bubbles initial position on erosion intensity is also investigated in this study, and it is found that bubbles with the highest distance from sheet cavity termination have the highest contribution in erosion rate.


Sign in / Sign up

Export Citation Format

Share Document