The location of the inception point of the air entrainment directly affects the energy dissipation ratio, the cavitation damage control, and the training wall height designs for a stepped spillway and a stilling basin...The location of the inception point of the air entrainment directly affects the energy dissipation ratio, the cavitation damage control, and the training wall height designs for a stepped spillway and a stilling basin. In this paper, the boundary layer theory of plates is used to predict the location of the inception point of the air entrainment over the stepped spillways by assuming the steps on the spillways as a kind of roughness. An empirical formula is presented based on the physical model experiments, with the maximum error less than 1% except at one point where the error is 1.6%, as compared to the experimental data. Meanwhile, it is shown that the location of the inception point of the air entrainment for the stepped spillway is much nearer to the top of the spillway than that for a smooth spillways, which explains why the high ratio of the energy dissipation is provided for the stepped spillway.展开更多
This paper is concerned with three-dimensional numerical simulation of a plunging liquid jet. The transient processes of forming an air cavity around the jet, capturing an initially large air bubble, and the break-up ...This paper is concerned with three-dimensional numerical simulation of a plunging liquid jet. The transient processes of forming an air cavity around the jet, capturing an initially large air bubble, and the break-up of this large toroidal-shaped bubble into smaller bubbles were analyzed. A stabilized finite element method (FEM) was employed under parallel numerical simulations based on adaptive, unstructured grid and coupled with a level-set method to track the interface between air and liquid. These simulations show that the inertia of the liquid jet initially depresses the pool's surface, forming an annular air cavity which surrounds the liquid jet. A toroidal liquid eddy which is subse- quently formed in the liquid pool results in air cavity collapse, and in turn entrains air into the liquid pool from the unstable annular air gap region around the liquid jet.展开更多
Trails induced by air entrainment of bow wave breaking can be observed clearly around an advancing ship,making it prong to be detected and causing reduction in hydrodynamic performance.In breaking bow wave region,diff...Trails induced by air entrainment of bow wave breaking can be observed clearly around an advancing ship,making it prong to be detected and causing reduction in hydrodynamic performance.In breaking bow wave region,different scale of coherent vortex structures related closely to the air entrainment are generated.The formation and evolution of bubble clouds can be accounted partly by the swirling vortex flow of the jet plunging.In this work,high-fidelity simulation of the bow wave breaking around a rectangular thin plate is performed with the in-house code BAMR-SJTU.A Liutex force model is constructed to investigate the interaction between the free surface and vortices,and to explore the possibility of alleviating or controlling the air entrainment.This study is the first attempt to manipulate vortices corresponding to the air entrainment of the bow wave breaking.Some typical vortex structures based on the Liutex vortex identification method and the distribution of the vortex force are presented.Comparison of the vortex structures and the entrapped bubbles is conducted to illustrate the efficiency of Liutex force model in air entrainment control.展开更多
In order to avoid the cavitation damage and the decrease of the energy dissipation of the stepped spillways with a large unit discharge, the air entrainments of the hydraulic jump aeration basin (HJAB) are theoretic...In order to avoid the cavitation damage and the decrease of the energy dissipation of the stepped spillways with a large unit discharge, the air entrainments of the hydraulic jump aeration basin (HJAB) are theoretically and experimentally investigated for the hydraulic-jump-stepped spillway developed by the authors. It is shown that the submerged degree of the hydraulic jump and the air concentration in the measuring section are all functions of the dimensionless discharge, the length and the end sill height of the H JAB. The submerged degree odecreases with the increase of the dimensionless discharge or the dimensionless length of the H JAB, but increases with the increase of the dimensionless height of the end sill of the HJAB. The flow regimes near the critical hydraulic jump, namely, at or= 1.0, have the best effect of the air entrainment for the flow at the measuring section and then that of the stepped spillway.展开更多
Based on the momentum integral equation of boundary layer on the steep slope and the equation of water depth deduced form both the equation of continuity including the effect of the boundary layer and the energy equat...Based on the momentum integral equation of boundary layer on the steep slope and the equation of water depth deduced form both the equation of continuity including the effect of the boundary layer and the energy equation outside the boundary layer on the steep slope, and with the aid of the function of velocity profile in the boundary layer, the water depth and boundary layer were numerically computed by Runge-Kutta method. The method determine parameters in the preceding equation were also discussed in detail for the laminar and turbulent flows. The water depth boundary layer and point of inception of air entrainment calculated by above method agree well with the data obtained by prototype.展开更多
Early in 1953 the experiments by Peterka proved that air entrainment has effects on decreasing cavitation damage. This technology has been widely used in the release works of high dams since the inception of air entra...Early in 1953 the experiments by Peterka proved that air entrainment has effects on decreasing cavitation damage. This technology has been widely used in the release works of high dams since the inception of air entrainment in the Grand Goulee Dam in 1960. Behavior, mechanism and application of air entrainment for cavitation damage control have been investigated for over half century. However, severe cavitation damage happened due to complex mechanism of air entrainment. The effects of air entrainment are related to many factors, including geometric parameters, hydraulic parameters and entrained air manners. In the present work an experimental set-up for air entrainment was specially designed, the behavior of reducing cavitation damage was experimentally investigated in the three aspects of entrained air pressure, air tube aera and air tube number. The results show that magnitude of reduction of cavitation damage is closely related to the entrained air tube number as well as entrained air pressure, air tube aera, and that the effect through three air tubes is larger than that through single air tube although the entrained air tubes have the same sum of tube aera, that is, 1 + 1 + 1 〉 3. Therefore, it is important to design an effective manner of air entrainment.展开更多
Physical modeling represents probably the oldest design tool in hydraulic engineering together with analytical approaches. In free surface flows, the similitude based upon a Froude similarity allows for a correct repr...Physical modeling represents probably the oldest design tool in hydraulic engineering together with analytical approaches. In free surface flows, the similitude based upon a Froude similarity allows for a correct representation of the dominant forces, namely gravity and inertia. As a result fluid flow properties such as the capillary forces and the viscous forces might be incorrectly reproduced, affecting the air entrainment and transport capacity of a high-speed model flow. Small physical models operating under a Froude similitude systematically underestimate the air entrainment rate and air-water interfacial properties. To limit scale effects, minimal values of Reynolds or Weber number have to be respected. The present article summarizes the physical background of such limitations and their combination in terms of the Morton number. Based upon a literature review, the existing limits are presented and discussed, resulting in a series of more conservative recommendations in terms of air concentration scaling. For other air-water flow parameters, the selection of the criteria to assess scale effects is critical because some parameters (e.g., bubble sizes, turbulent scales) can be affected by scale effects, even in relatively large laboratory models.展开更多
The cavitation erosion is one of the common damage modes in water engineering.The study of the interaction among cavitation bubbles,air bubbles and the wall is of great significance for understanding the mechanism of ...The cavitation erosion is one of the common damage modes in water engineering.The study of the interaction among cavitation bubbles,air bubbles and the wall is of great significance for understanding the mechanism of the air entrainment to alleviate the cavitation and to enhance the cavitation erosion mitigation effect of aeration.By using the high-speed camera,the regular patterns of the collapse of cavitation bubbles in the vicinity of the wall and the air bubbles are studied in this paper.It is shown that in the vicinity of air bubbles,the cavitation bubbles may only collapse towards or from air bubbles,while under the dual impacts of air bubbles and the wall,the direction of collapse of the cavitation bubbles depends on the combined vector of the impact forces of the air bubbles and the wall.When the air bubbles are very close to the cavitation bubbles,the air bubbles will be penetrated and stretched by the cavitation bubbles,when the distance between them is short enough,the cavitation bubbles and the air bubbles will connect to form gas-type cavitation bubbles,and the collapse strength will be significantly decreased,when there are two air bubbles near the cavitation bubbles,they may even be penetrated successively by the cavitation bubbles.Thus it can be concluded that during the air entrainment to alleviate the cavitation,the number of air bubbles is more important than the total air concentration.展开更多
The emergency transformation of various aspects of life and business these days requires prompt evaluation of autonomous vehicles.One of the primary reassessments deals with the applicability of the vehicle passive sa...The emergency transformation of various aspects of life and business these days requires prompt evaluation of autonomous vehicles.One of the primary reassessments deals with the applicability of the vehicle passive safety system to the protec-tion of arbitrarily positioned passengers.To mitigate possible risks caused by the simultaneous deployment of several big airbags,a new principle of their operation is required.Herein,the aspirated inflator for a driver airbag is developed that can provide 50L-airbag inflation within 30-40 ms.As a result,about 3/4 of the air is to be entrained into an airbag from the vehicle compartment.The process is initiated by a supersonic pulse jet(1/3 air volume)generated pyrotechnically.Then the Prandtl-Meyer problem formulation enables guiding linear and angular dimensions of the essential parts of the device.Accordingly,a family of experimental models of varied geometry is fabricated and tested to determine their operational effectiveness in a range of motive pressure within~3-7 MPa.Experiments are performed on a specially designed facility equipped with compressed-air tanks and a high-speed valve to mimic the inflator operation with the pyrotechnic gas generator.The aspirated inflator operability is characterized using multivariate measurements of pressure fields,high-speed video-recording of the airbag inflation process,and evaluation of aspiration(entrainment)ratio.The average volume aspiration ratio measured at 300 K is found to reach 2.8 and it’s expected to almost double at 1200 K.展开更多
Experimental study on aeration characteristics of various aeration devices was conducted in the spillway tunnel of the Pubugou hydropower project, Sichuan Province, China. It is shown by comparison that the new type o...Experimental study on aeration characteristics of various aeration devices was conducted in the spillway tunnel of the Pubugou hydropower project, Sichuan Province, China. It is shown by comparison that the new type of aeration device, namely, the aerator with a trapezoidal-shaped slot and a steep-slope section(ATSS), can avoid water accumulation in the cavity of the aeration device in the project, thus can effectively solve the backwater problems arising from this project and be used for a wide range of different water levels, without any drain facilities. Above the water level of 840 m, the water contained in the cavity can be eliminated completely, which means that the recommended new type of aerator can meet the aeration demands in the spillway of the project with low Froude number and may be of practical significance and of interest to other projects with similar types of aeration devices.展开更多
Air entrainment is an effective approach to protect release works from cavitation damage. The traditional method of aerator device designs is that, for given flow conditions, the geometries of the aerator device are d...Air entrainment is an effective approach to protect release works from cavitation damage. The traditional method of aerator device designs is that, for given flow conditions, the geometries of the aerator device are designed and then the effec(s are experimentally tested for cavitation damage control. The present paper proposes an inverse problem method of determining the bottom slopes in front of and behind an aerator if the requirements of air entrainment, flow conditions and some of aerator geometric parameters are given. An RBF neural network model is developed and the relevant bottom slopes are calculated in different conditions of flow and geometry on the basis of the data of 19 aerator devices from different discharge tunnels with safe operation. The case study shows that the methodology provides an effective way to design aerator devices under given target conditions.展开更多
This paper presents entrainment mechanism, and transported and diffusioncharacteristics at the point of entry of submerged jet. The profiles of both velocity andconcentration within the air-water mixing layer were the...This paper presents entrainment mechanism, and transported and diffusioncharacteristics at the point of entry of submerged jet. The profiles of both velocity andconcentration within the air-water mixing layer were theoretically deduced. And the comparisonsbetween theoretical values and measured data were made. Results show that the velocity profilewithin the air-water mixing layer exhibits a form of error function. The concentrations of airentrainment in the internal and external regions of air-water mixing layer correspond to Gaussiandistribution.展开更多
The interaction of multiple fires may lead to a higher flame height and more intense radiation flux than a single fire,which increases the possibility of flame spread and risks to the surroundings.Experiments were con...The interaction of multiple fires may lead to a higher flame height and more intense radiation flux than a single fire,which increases the possibility of flame spread and risks to the surroundings.Experiments were conducted using three burners with identical heat release rates(HRRs)and propane as the fuel at various spacings.The results show that flames change from non-merging to merging as the spacing decreases,which result in a complex evolution of flame height and merging point height.To facilitate the analysis,a novel merging criterion based on the dimensionless spacing S/z_(c) was proposed.For non-merging flames(S/z_(c)>0.368),the flame height is almost identical to a single fire;for merging flames(S/z_(c)≤0.368),based on the relationship between thermal buoyancy B and thrust P(the pressure difference between the inside and outside of the flame),a quantitative analysis of the flame height,merging point height,and air entrainment was formed,and the calculated merging flame heights show a good agreement with the measured experimental values.Moreover,the multi-point source model was further improved,and radiation fraction of propane was calculated.The data obtained in this study would play an important role in calculating the external radiation of propane fire.展开更多
The present paper attempts to solve equations in the initial stage and the two-phase flow regime of fuel spray penetration using the HPM-Padétechnique,which is a combination of the homotopy perturbation method(HP...The present paper attempts to solve equations in the initial stage and the two-phase flow regime of fuel spray penetration using the HPM-Padétechnique,which is a combination of the homotopy perturbation method(HPM)and Padéapproximation.At the initial stage,the effects of the droplet drag and the air entrainment were explained while in the two-phase flow stage,the spray droplets had the same velocities as the entrained air.The results for various injection pressures and ambient densities are presented graphically and then discussed upon.The obtained results for these two stages show a good agreement with previously obtained expressions via successive approximations in the available literature.The numerical result indicates that the proposed method is straight forward to implement,efficient and accurate for solving nonlinear equations of fuel spray.展开更多
基金Project supported by the National Natural Science Foundation of China(Grant No.51179114)
文摘The location of the inception point of the air entrainment directly affects the energy dissipation ratio, the cavitation damage control, and the training wall height designs for a stepped spillway and a stilling basin. In this paper, the boundary layer theory of plates is used to predict the location of the inception point of the air entrainment over the stepped spillways by assuming the steps on the spillways as a kind of roughness. An empirical formula is presented based on the physical model experiments, with the maximum error less than 1% except at one point where the error is 1.6%, as compared to the experimental data. Meanwhile, it is shown that the location of the inception point of the air entrainment for the stepped spillway is much nearer to the top of the spillway than that for a smooth spillways, which explains why the high ratio of the energy dissipation is provided for the stepped spillway.
基金supported by the Office of Naval Research(Grant ONRDC14292111)
文摘This paper is concerned with three-dimensional numerical simulation of a plunging liquid jet. The transient processes of forming an air cavity around the jet, capturing an initially large air bubble, and the break-up of this large toroidal-shaped bubble into smaller bubbles were analyzed. A stabilized finite element method (FEM) was employed under parallel numerical simulations based on adaptive, unstructured grid and coupled with a level-set method to track the interface between air and liquid. These simulations show that the inertia of the liquid jet initially depresses the pool's surface, forming an annular air cavity which surrounds the liquid jet. A toroidal liquid eddy which is subse- quently formed in the liquid pool results in air cavity collapse, and in turn entrains air into the liquid pool from the unstable annular air gap region around the liquid jet.
基金Project supported by the National Key Research and Development Program of China(Grant No.2019YFB1704200)the National Natural Science Foundation of China(Grant Nos.51879159,52131102).
文摘Trails induced by air entrainment of bow wave breaking can be observed clearly around an advancing ship,making it prong to be detected and causing reduction in hydrodynamic performance.In breaking bow wave region,different scale of coherent vortex structures related closely to the air entrainment are generated.The formation and evolution of bubble clouds can be accounted partly by the swirling vortex flow of the jet plunging.In this work,high-fidelity simulation of the bow wave breaking around a rectangular thin plate is performed with the in-house code BAMR-SJTU.A Liutex force model is constructed to investigate the interaction between the free surface and vortices,and to explore the possibility of alleviating or controlling the air entrainment.This study is the first attempt to manipulate vortices corresponding to the air entrainment of the bow wave breaking.Some typical vortex structures based on the Liutex vortex identification method and the distribution of the vortex force are presented.Comparison of the vortex structures and the entrapped bubbles is conducted to illustrate the efficiency of Liutex force model in air entrainment control.
基金Project supported by the National Natural Science Foundation of China(Grant No.51479057)the Research Innovation Program for College Graduates of Jiangsu Province(Grant No.2015B41714)
文摘In order to avoid the cavitation damage and the decrease of the energy dissipation of the stepped spillways with a large unit discharge, the air entrainments of the hydraulic jump aeration basin (HJAB) are theoretically and experimentally investigated for the hydraulic-jump-stepped spillway developed by the authors. It is shown that the submerged degree of the hydraulic jump and the air concentration in the measuring section are all functions of the dimensionless discharge, the length and the end sill height of the H JAB. The submerged degree odecreases with the increase of the dimensionless discharge or the dimensionless length of the H JAB, but increases with the increase of the dimensionless height of the end sill of the HJAB. The flow regimes near the critical hydraulic jump, namely, at or= 1.0, have the best effect of the air entrainment for the flow at the measuring section and then that of the stepped spillway.
文摘Based on the momentum integral equation of boundary layer on the steep slope and the equation of water depth deduced form both the equation of continuity including the effect of the boundary layer and the energy equation outside the boundary layer on the steep slope, and with the aid of the function of velocity profile in the boundary layer, the water depth and boundary layer were numerically computed by Runge-Kutta method. The method determine parameters in the preceding equation were also discussed in detail for the laminar and turbulent flows. The water depth boundary layer and point of inception of air entrainment calculated by above method agree well with the data obtained by prototype.
基金supported by the National Natural Science Foundation of China (Grant No.50879021)the National Science Fund for Distinguished Young Scholars (Grant No.50925932)the Ministry of Science and Technology of China (Grant No.2008BAB19B04)
文摘Early in 1953 the experiments by Peterka proved that air entrainment has effects on decreasing cavitation damage. This technology has been widely used in the release works of high dams since the inception of air entrainment in the Grand Goulee Dam in 1960. Behavior, mechanism and application of air entrainment for cavitation damage control have been investigated for over half century. However, severe cavitation damage happened due to complex mechanism of air entrainment. The effects of air entrainment are related to many factors, including geometric parameters, hydraulic parameters and entrained air manners. In the present work an experimental set-up for air entrainment was specially designed, the behavior of reducing cavitation damage was experimentally investigated in the three aspects of entrained air pressure, air tube aera and air tube number. The results show that magnitude of reduction of cavitation damage is closely related to the entrained air tube number as well as entrained air pressure, air tube aera, and that the effect through three air tubes is larger than that through single air tube although the entrained air tubes have the same sum of tube aera, that is, 1 + 1 + 1 〉 3. Therefore, it is important to design an effective manner of air entrainment.
基金supported by the Australian Re-search Council(Grant No.DP120100481)
文摘Physical modeling represents probably the oldest design tool in hydraulic engineering together with analytical approaches. In free surface flows, the similitude based upon a Froude similarity allows for a correct representation of the dominant forces, namely gravity and inertia. As a result fluid flow properties such as the capillary forces and the viscous forces might be incorrectly reproduced, affecting the air entrainment and transport capacity of a high-speed model flow. Small physical models operating under a Froude similitude systematically underestimate the air entrainment rate and air-water interfacial properties. To limit scale effects, minimal values of Reynolds or Weber number have to be respected. The present article summarizes the physical background of such limitations and their combination in terms of the Morton number. Based upon a literature review, the existing limits are presented and discussed, resulting in a series of more conservative recommendations in terms of air concentration scaling. For other air-water flow parameters, the selection of the criteria to assess scale effects is critical because some parameters (e.g., bubble sizes, turbulent scales) can be affected by scale effects, even in relatively large laboratory models.
基金supported by the National Key Research and Development Program of China(Grant No.2016YFC0401901)the National Natural Science Foundation of China(Grant No.51409180).
文摘The cavitation erosion is one of the common damage modes in water engineering.The study of the interaction among cavitation bubbles,air bubbles and the wall is of great significance for understanding the mechanism of the air entrainment to alleviate the cavitation and to enhance the cavitation erosion mitigation effect of aeration.By using the high-speed camera,the regular patterns of the collapse of cavitation bubbles in the vicinity of the wall and the air bubbles are studied in this paper.It is shown that in the vicinity of air bubbles,the cavitation bubbles may only collapse towards or from air bubbles,while under the dual impacts of air bubbles and the wall,the direction of collapse of the cavitation bubbles depends on the combined vector of the impact forces of the air bubbles and the wall.When the air bubbles are very close to the cavitation bubbles,the air bubbles will be penetrated and stretched by the cavitation bubbles,when the distance between them is short enough,the cavitation bubbles and the air bubbles will connect to form gas-type cavitation bubbles,and the collapse strength will be significantly decreased,when there are two air bubbles near the cavitation bubbles,they may even be penetrated successively by the cavitation bubbles.Thus it can be concluded that during the air entrainment to alleviate the cavitation,the number of air bubbles is more important than the total air concentration.
文摘The emergency transformation of various aspects of life and business these days requires prompt evaluation of autonomous vehicles.One of the primary reassessments deals with the applicability of the vehicle passive safety system to the protec-tion of arbitrarily positioned passengers.To mitigate possible risks caused by the simultaneous deployment of several big airbags,a new principle of their operation is required.Herein,the aspirated inflator for a driver airbag is developed that can provide 50L-airbag inflation within 30-40 ms.As a result,about 3/4 of the air is to be entrained into an airbag from the vehicle compartment.The process is initiated by a supersonic pulse jet(1/3 air volume)generated pyrotechnically.Then the Prandtl-Meyer problem formulation enables guiding linear and angular dimensions of the essential parts of the device.Accordingly,a family of experimental models of varied geometry is fabricated and tested to determine their operational effectiveness in a range of motive pressure within~3-7 MPa.Experiments are performed on a specially designed facility equipped with compressed-air tanks and a high-speed valve to mimic the inflator operation with the pyrotechnic gas generator.The aspirated inflator operability is characterized using multivariate measurements of pressure fields,high-speed video-recording of the airbag inflation process,and evaluation of aspiration(entrainment)ratio.The average volume aspiration ratio measured at 300 K is found to reach 2.8 and it’s expected to almost double at 1200 K.
文摘Experimental study on aeration characteristics of various aeration devices was conducted in the spillway tunnel of the Pubugou hydropower project, Sichuan Province, China. It is shown by comparison that the new type of aeration device, namely, the aerator with a trapezoidal-shaped slot and a steep-slope section(ATSS), can avoid water accumulation in the cavity of the aeration device in the project, thus can effectively solve the backwater problems arising from this project and be used for a wide range of different water levels, without any drain facilities. Above the water level of 840 m, the water contained in the cavity can be eliminated completely, which means that the recommended new type of aerator can meet the aeration demands in the spillway of the project with low Froude number and may be of practical significance and of interest to other projects with similar types of aeration devices.
基金Project supported by the National Natural Science Function of China(Grant No.51179114)
文摘Air entrainment is an effective approach to protect release works from cavitation damage. The traditional method of aerator device designs is that, for given flow conditions, the geometries of the aerator device are designed and then the effec(s are experimentally tested for cavitation damage control. The present paper proposes an inverse problem method of determining the bottom slopes in front of and behind an aerator if the requirements of air entrainment, flow conditions and some of aerator geometric parameters are given. An RBF neural network model is developed and the relevant bottom slopes are calculated in different conditions of flow and geometry on the basis of the data of 19 aerator devices from different discharge tunnels with safe operation. The case study shows that the methodology provides an effective way to design aerator devices under given target conditions.
文摘This paper presents entrainment mechanism, and transported and diffusioncharacteristics at the point of entry of submerged jet. The profiles of both velocity andconcentration within the air-water mixing layer were theoretically deduced. And the comparisonsbetween theoretical values and measured data were made. Results show that the velocity profilewithin the air-water mixing layer exhibits a form of error function. The concentrations of airentrainment in the internal and external regions of air-water mixing layer correspond to Gaussiandistribution.
基金supported by the National Natural Science Foundation of China(Grant Nos.52036009 and 51722605).
文摘The interaction of multiple fires may lead to a higher flame height and more intense radiation flux than a single fire,which increases the possibility of flame spread and risks to the surroundings.Experiments were conducted using three burners with identical heat release rates(HRRs)and propane as the fuel at various spacings.The results show that flames change from non-merging to merging as the spacing decreases,which result in a complex evolution of flame height and merging point height.To facilitate the analysis,a novel merging criterion based on the dimensionless spacing S/z_(c) was proposed.For non-merging flames(S/z_(c)>0.368),the flame height is almost identical to a single fire;for merging flames(S/z_(c)≤0.368),based on the relationship between thermal buoyancy B and thrust P(the pressure difference between the inside and outside of the flame),a quantitative analysis of the flame height,merging point height,and air entrainment was formed,and the calculated merging flame heights show a good agreement with the measured experimental values.Moreover,the multi-point source model was further improved,and radiation fraction of propane was calculated.The data obtained in this study would play an important role in calculating the external radiation of propane fire.
文摘The present paper attempts to solve equations in the initial stage and the two-phase flow regime of fuel spray penetration using the HPM-Padétechnique,which is a combination of the homotopy perturbation method(HPM)and Padéapproximation.At the initial stage,the effects of the droplet drag and the air entrainment were explained while in the two-phase flow stage,the spray droplets had the same velocities as the entrained air.The results for various injection pressures and ambient densities are presented graphically and then discussed upon.The obtained results for these two stages show a good agreement with previously obtained expressions via successive approximations in the available literature.The numerical result indicates that the proposed method is straight forward to implement,efficient and accurate for solving nonlinear equations of fuel spray.