Research Article  Open Access
Yanjie Ma, Futing Bao, Lin Sun, Yang Liu, Weihua Hui, "A New Erosive Burning Model of Solid Propellant Based on Heat Transfer Equilibrium at Propellant Surface", International Journal of Aerospace Engineering, vol. 2020, Article ID 8889333, 9 pages, 2020. https://doi.org/10.1155/2020/8889333
A New Erosive Burning Model of Solid Propellant Based on Heat Transfer Equilibrium at Propellant Surface
Abstract
Erosive burning refers to the augmentation of propellant burning rate appears when the velocity of combustion gas flowing parallel to the propellant surface is relatively high. Erosive burning can influence the total burning rate of propellant and performance of solid rocket motors dramatically. There have been many different models to evaluate erosive burning rate for now. Yet, due to the complication processes involving in propellant and solid rocket motor combustion, unknown constants often exist in these models. To use these models, trialanderror procedure must be implemented to determine the unknown constants firstly. This makes many models difficult to estimate erosive burning before plenty of experiments. In this paper, a new erosive burning rate model is proposed based on the assumption that the erosive burning rate is proportional to the heat flux at the propellant surface. With entrance effect, roughness, and transpiration considered, convective heat transfer coefficient correlation proposed in recent years is used to compute the heat flux. This allows the release of unknown constants, making the model universal and easy to implement. The computational data of the model are compared with different experimental and computational data from different models. Results show that good accuracy (10%) with experiments can be achieved by this model. It is concluded that the present model could be used universally for erosive burning rate evaluation of propellant and performance prediction of solid rocket motor as well.
1. Introduction
The performance of solid rocket motor (SRM) can be improved by raising the ratio of length to diameter. Yet, large length to diameter ratio would lead to higher port velocity, which in turn may cause substantial augment of the burning rate of the propellant. This augment of the burning rate is called erosive burning. When erosive burning occurs, the performance of SRMs becomes difficult to predict. Many mission failures are related to erosive burning. Thus, theory models are needed to explain the mechanism of erosive burning and predict the erosive burning rate of propellant. For decades, many researches on erosive burning are implemented. Some reviews are Ref. [1–5]. It is widely agreed that the enhancement of heat transfer from core gas in the port to the propellant surface is the reason why erosive burning occurs. Yet, because of the complexity of physical and chemical processes involved in the propellant combustion, the mechanism of erosive burning has not been revealed thoroughly. For now, empirical models are the most popular models in use.
The simplest series of models evaluating the erosive burning rate might be the linearcorrelation models. This type of models has the following formation [5]. where denotes the erosive burning ratio (the ratio between total burning rate to normal burning rate). The variable denotes different variables in different models, such as velocity of the core gas flow, Mach number, and mass flux of the port. The variable is the threshold value of , as the erosive burning occurs only when exceeds a threshold value. The constant is unknown and to be determine. Generally, and vary according to the propellant and the geometry of SRMs. They are often determined by trialanderror process, where computational data are compared to the experimental data repeatedly to obtain the best combination of and values.
LR model, a more precise and may be the most widely used model, is proposed by Lenoir and Robillard [6]. In LR model, the total burning rate is composed of two independent components. One is the normal burning rate depending on pressure according to Saint Robert’s law. The other component is the erosive burning rate depending on the gas flow condition. The authors postulated that the erosive burning rate is proportional to the heat transfer coefficient on the propellant surface and derived a correlation between erosive burning rate and gas parameters.
Different modifications of LR model were made by researchers. King [7] accounts the effect of heat feedback rate on the normal burning rate and proposed a model with comparable accuracy with LR correlation. Lawrence and coauthors [8] used port diameter as the characterized length in their model, rather than the length from the head in original LR model. Comparisons to the experimental data showed that Lawrence’s model can lead to more accurate results than original LR model does. This makes it popular in analyses of erosive burning, internal ballistic of SRMs, and ignition process [9–12].
Despite the enormous number of models trying to evaluate the erosive burning rate, there is an obstacle to overcome. In the models of linearcorrelation series and LR series (and many other models), there were often two or more unknown constants. These constants are dependent on propellant and geometry and sometimes other parameters. Trialanderror procedures are generally needed to obtain the values of these constants. This makes it difficult to predict erosive burning behavior of propellant before firing tests (that should have been one of objectives a model is established for). In 1997, Mukunda and Paul [13] developed a model for erosive burning. They modified their model [14] for threedimensional SRMs later, and the modified model was supported by Javed’s computation [15]. Although MP model was claimed to be universal with two constants determined, these constants are actually obtained by comparing with plenty of experimental data (another form of trialanderror procedure, indeed).
In present paper, the erosive burning rate is assumed to be proportional to the heat flux, rather than the heat transfer coefficient in LR model [6]. By computing the convection heat transfer coefficient using Gnielinski model [16] with the entrance effect, the surface roughness, and the transpiration from propellant surface considered, no unknown constants is introduced in this model. This makes the model easy to apply, universal, and having the potential to predict propellant erosive burning behavior and performance of SRMs before any firing tests. Comparison of computational results of this model show good agreement (10%) with experimental data and results from other models.
2. Model
As in Ref. [6], the total burning rate is postulated as the sum of normal burning rate and erosive burning rate in this paper.
Using Saint Robert’s law, the normal burning rate can be expressed as
The erosive burning rate is assumed to be proportional to the heat flux at the burning surface. where is proportionality coefficient, and is the heat flux. When the propellant is burning, the heat flux from gas to propellant surface results in the erosive burning rate. The energy balance at the propellant surface can be expressed as
In the righthand term of Equation (5), the heat release from the reaction on surface is ignored as in Ref. [17, 18]. The same equation was also used in LR model [6]. Combining Equation (4) and Equation (5), we get the expression of ,
Now, the total burning rate can be written as where the correlation is used.
2.1. Heat Transfer Coefficient with No Transpiration
As shown in Equation (7), the heat transfer coefficient at the propellant is needed for erosive burning rate computation. Erosive burning only occurs when the gas velocity is relatively high. In such situation, the convection takes the dominant role at the surface heat transfer. So, radiative and conductive heat transfer is neglected in this paper. Gnielinski [16] proposed an empirical correlation for convection. This correlation has good accuracy for both laminar flow and turbulent flow. During the burning process, the propellant surface temperature generally holds at a constant value. Thus, the boundary can be seen as a universal wall temperature condition (UWT). Under this condition, when the flow is laminar , Nusselt number can be computed by the following equation. where , , and .
When the flow is turbulent , Gnielinski [16] suggests
In Equation (9), entrance effect is considered. For gas, the property modification term can be written as , where for air. Darcy friction factor is computed with roughness considered (see Section 2.2).
When the flow is in transition state , linear interpolation of Equation (8) and Equation (9) should be used, as where the coefficient .
Now that is obtained, the heat transfer coefficient with no transpiration can be easily computed from
2.2. Darcy Friction Fraction with Effect of Roughness
Roughness of propellant surface has substantial effect on heat flux by influencing the value of Darcy friction fraction. For smooth tube, Prandtl correlation can be used to compute Darcy friction factor.
Colebrook [19] proposed an explicit form with good accuracy (less than 1%).
For rough tube, Von Karman correlation has good accuracy.
Haaland [20] combined Equation (13) and Equation (14) and proposed the following correlation for Darcy friction factor.
This correlation is validated for smooth tube, as well as rough tube.
2.3. Effect of Transpiration
When the propellant burns, the transpiration flow will decrease the heat transfer from core gas flow to the propellant surface. The relationship between the heat transfer coefficients with and without transpiration can be expressed as [18, 21]. with . The variable is heat transfer coefficient with no transpiration and computed using Equation (9), Equation (11), and Equation (15). The variable is the velocity of gas generated at propellant surface. Since the burning rate is neglectable comparing with , the equation and stands.
Now that heat transfer coefficient has been computed, Equation (7) can be resolved for total burning rate.
3. Validation of the Model
In this section, the total burning rate of propellant specimens at different erosive conditions and pressure history of SRMs with large lengthtodiameter ratio are computed using the erosive burning rate model described in Section 2. During the computation, some transport properties of the combustion gas such as viscosity and Prandtl number are needed for erosive burning computation. To obtain these parameters, software package Chemical Equilibrium with Applications (CEA) [22] is used. The results are compared with experimental data and numerical data from different researches.
3.1. Data of King
King [23–25] conducted a series of experiments on different composite propellants. In these experiments, effects of different factors, such as AP particle size, oxidizerfuel ratio, combustion gas temperature, and catalyst, were studied by experiments on different propellants. Properties of six different propellants of King’s paper are listed in Table 1. The total burning rate of these six different propellants of King’s paper is computed under different pressure and gas velocity conditions.
From the results shown in Figures 1–6, we can see that the computational data are in good agreement with King’s experimental data. As revealed by King’s papers (and many other researches, see Ref. [1–5]), propellants with comparable normal burning rate will have similar erosive burning characteristics. When a propellant has higher normal burning rate, no matter by decreasing AP particle diameter (comparison between propellant 4525 and 5051), using catalyst (propellant 4869), or changing the oxidizerfuel ratio (propellant 5542), the erosive burning becomes less substantial. Comparison of results from propellant 4525 and 5051 shows that as the AP particle diameter increases (from 20 μm of 4525 to 200 μm of 5051), the normal burning rate decreases, and the erosive burning ratio increases. At 6 MPa, the normal burning rate of 4525 is 0.83 cm/s, and that of 5051 is 0.51 cm/s. At the same pressure 6 MPa and gas velocity 183 m/s, the values of erosive burning ratio of 4525 and 5051 are 1.49 and 1.78, respectively. In propellant 4869, small amount of catalyst Fe_{2}O_{3} is used, leading to an increment of normal burning rate, while oxidizerfuel ratio and temperature remains almost the same as 4525. Results of Figures 1 and 3 show that at pressure 6 MPa and gas velocity 244 m/s, erosive burning ratios for 4525 and 4869 are 1.72 and 1.19, respectively. This indicates that with the normal burning rate increasing because of the catalyst, the erosive burning ratio decreases. Propellant 5542T has the same AP particle diameter with 4525 and larger oxidizerfuel ratio. This leads to higher flame temperature and higher normal burning rate (normal burning rate of 5542T is 1.23 cm/s at 6 MPa). The burning rate results of Figures 1 and 4 show that at 6 MPa and gas velocity 244 m/s, the erosive burning ratio of 5542T is 1.15, comparing with 1.72 of propellant 4525. The fact indicates that with oxidizerfuel ratio changes, normal burning rate and flame temperature would be affected, and again, the erosive burning ratio decreases with increasing normal burning rate. On the other hand, when normal burning rates of different propellants get close, the propellants would have comparable erosive burning behavior. Propellant 5565T differs from 4525 in oxidizerfuel ratio. Yet, by adjusting the particle size, 5565T has similar normal burning rate as 4525. As shown in Figure 5, at 6 MPa and gas velocity 244 m/s, the erosive burning ratio of 5565T is 1.51, not far from the value of 1.72 of 4525. The only difference between formulations 5555T and 5565T is the AP particle diameter distribution. With more small particles, 5555T has higher normal burning rate than 5565T. Again, results in Figures 5 and 6 show that the erosive burning ratio of 5555T is smaller than 5565T. In fact, erosive burning is negligibly small even when gas velocity achieves 305 m/s.
Another conclusion from King is that the temperature of core gas has very little effect on erosive burning. This was also mentioned by Marklund and Lake [26]. King, as well as Marklund and Lake, claimed that this phenomenon can prove that LR model’s flaw. The constant in LR model depends strongly on core gas temperature, and this would lead to obvious burning rate increment when core gas temperature increases, which is against to the experimental data. Yet, as Landsbaum [3] indicated, these researchers did not consider the reduction effect of temperature increment on mass flux. When core gas temperature increases, the mass flux decreases. According to LR model, lower mass flux would lead to lower erosive burning rate directly. As different factors would be affected by core gas temperature, it is a little difficult to determine the total effect of core gas temperature on burning rate precisely. Figure 7 depicts the computational data of erosive burning ratio of propellant 4525 at different core gas temperatures. When the velocity is not very high, the core gas temperature causes neglectable effect to erosive burning. At velocity 455 m/s and pressure 10 MPa, the erosive burning ratio at higher core gas temperature is 2.959, comparing with 2.832 at lower core gas temperature. This is consistent with experimental data from different researches [24, 26]. As the velocity gets higher, the effect of the core gas temperature gets more significant, following same trend as King’s model [24].
All the above results show that the model proposed in present paper can provide a good agreement with experiments. The regularity of erosive burning can be correctly revealed by the model.
3.2. Data of Marklund and Lake
Marklund and Lake [26] implemented experiments on different propellants, two of which, propellant A and propellant C, are computed in this paper. The parameters of these two propellants are listed in Table 1. Computational results of total burning rate at different pressure and flow velocity are shown in Figures 8 and 9, along with the experimental data from Ref [26]. It can be seen that the computational data are in good agreement with experimental data. Since they have different oxidizerfuel ratios, the flame temperature and normal burning rate are different for these two propellants. At 10 MPa, the normal burning rate of propellant A is 0.47 cm/s and that of propellant C is 0.85 cm/s. From the results of Figures 8 and 9, we can see that at 10 MPa and gas velocity 103 m/s, the total burning rate of propellant A is 0.95 cm/s, leading to an erosive burning ratio 2.02. At the same pressure 10 MPa and a little higher gas velocity 113 m/s, the total burning rate of propellant C is 1.09 cm/s, making the erosive burning ratio 1.28, much lower than that of propellant A. This again shows that the propellant with lower normal burning rate is likely to have larger erosive burning ratio.
3.3. Data of Hasegawa
Experiments on different SRMs with cylindrical grain were implemented by Hasegawa [27]. In these experiments, Xray was used to obtain the regression of propellant and burning rate at different positions. The headend pressure was recorded for each SRM. Parameters of the propellant (propellant 1 in Ref [27]) are shown in Table 1. The geometric parameters of SRMs are listed in Table 2.

In this paper, a onedimensional (1D) transient computational fluid dynamic (CFD) program is developed and used to compute the internal ballistic of the SRMs. SIMPLE scheme [28, 29] is used to solve the 1D NavierStokes equations. The total burning rate of propellant was computed using the model described in Section 2. Computational results of pressure history at headend of three SRMs are shown in Figures 10–12, along with the experimental and computational data from Hasegawa [27] and computational data from Javed [15]. It is worth noting that the computational results from Javed using MP model and from Hasegawa start with a relatively high pressure. This indicates that they ignored the ignition process of the SRM and might use quasitransient scheme, not transient scheme as this paper. Besides, for Hasegawa model, a complex trialanderror procedure to determine constants in the model has been implemented before the computation of SRMs.
Figure 10 depicts the pressure history curves of motor A of Ref [27]. Since the ignition process was ignored, the pressure peak at the beginning was missing in both Javed’s results and Hasegawa’s results. The computational results of this paper meet well with the experimental data in this stage. After ignition, a dramatic pressure drop occurs, where better results are obtained using MP model and Hasegawa model. As time goes, MP model seems has a trend to overpredict the total burning rate. Hasegawa model give the best results out of the three models. The precision of the model in this paper seems to be between the others until the moment just before tailing stage of the whole experiment. With ignition process excluded, the maximum difference between results of Javed and the experimental data is about 10%. The value from both of other two models is 3%.
Results of motor B of Ref [27] are shown in Figure 11. For this motor, all the three models have a deviation about 9% from experiment when the pressure drop occurs after the ignition. After that, computation of this paper gives the best results. The maximum deviation of Javed’s results is about 20%, occurring just before the tail begins, and that of Hasegawa’s computation is about 5% at the same time. For the results of this paper, the maximum deviation 3.6% after the pressure drop appears in the middle of the experiment.
All of the three models give good results for motor C of Ref [27]. The model of this paper shows the best results. MP model overpredicted the burning rate again, and this time, the overprediction appears at the beginning of the experiment, while results from the other models are close to experimental data and to each other. Reason for this might be the relatively low mass flux of this motor and low erosive burning accordingly, which makes the prediction of erosive burning a little easier. For this motor, the maximum deviations for MP model, Hasegawa model, and the model of this paper are 9%, 2.3%, and 1.4%, respectively.
From the comparisons with experimental data and results of different models, it can be seen that the model in this paper can provide good prediction for internal ballistics of SRMs. Considering the fact that the model is universal like MP model, not need to be adjusted by trialanderror procedure like Hasegawa model and many other models, it is a promising model, indeed.
4. Conclusion
A new erosive burning rate model is proposed based on the assumption that the erosive burning rate is proportional to the heat flux on the propellant surface. The heat flux is computed using convective heat transfer correlation developed in recent years. Unlike many other models (such as models with linear form and models of LR series), there are no unknown constants existing in the model, which makes the model universal and easy to implement. Using the model, it is possible to predict propellant burning rate and performance of SRMs before experiments. Comparison with experimental data and other models are implemented for propellant specimens and for real SRMs under different conditions. The computational results are in good agreement with experimental data. Results of propellant specimens also show that erosive burning regularities well accepted by researchers can be revealed correctly by this model. Results of motors indicate that the model can predict the pressure history with high accuracy, with deviation within 10%.
Nomenclature
:  Coefficient in Saint Robert’s law 
:  Specific heat of gas 
:  Specific heat of propellant 
:  Diameter of the port 
:  Darcy friction factor 
:  Proportionality coefficient in the model of this paper 
:  Proportionality coefficient in linear erosive burning model 
:  Heat transfer coefficient 
:  Heat transfer coefficient with no transpiration 
:  Proportionality coefficient 
:  Length of the port 
:  Exponent in Saint Robert’s law 
:  Nusselt number 
:  Pressure 
:  Prandtl number 
:  Heat flux at propellant surface 
:  Total burning rate 
:  Normal burning rate (burning rate with no erosive burning) 
:  Erosive burning rate 
:  Reynolds number 
:  Temperature 
:  Independent variable in linear erosive burning model 
:  Erosive burning ratio; roughness of the surface 
:  Thermal conductivity 
:  Velocity 
:  Density of gas phase 
:  Density of propellant. 
0:  Normal state 
:  Core gas flow 
:  Erosive burning 
:  Gas phase 
:  Initial state 
lam:  Laminar 
:  Interface between propellant and combustion gas 
turb:  Turbulent. 
Data Availability
All data, models, or code generated or used during the study are available from the corresponding author by request.
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
Acknowledgments
This work was financially supported by the Shanghai Academy of Spaceflight Technology Innovation Fund (Program no. SAST 2019107), the Fundamental Research Funds for the Central Universities, Northwestern Polytechnical University (Program no. 3102019HTXM006), and Natural Science Basic Research Program of Shaanxi (Program no. 2020JQ172).
References
 M. K. Razdan and K. K. Kuo, “Erosive burning of solid propellants,” Fundamentals of SolidPropellant Combustion, vol. 90, pp. 515–598, 1984. View at: Google Scholar
 M. K. King, “Erosive burning of solid propellants,” Journal of Propulsion and Power, vol. 9, no. 6, pp. 785–805, 1993. View at: Publisher Site  Google Scholar
 E. M. Landsbaum, “Erosive burning of solid rocket propellantsa revisit,” Journal of Propulsion and Power, vol. 21, no. 3, pp. 470–477, 2005. View at: Publisher Site  Google Scholar
 T. L. Jackson, “Modeling of heterogeneous propellant combustion: a survey,” AIAA Journal, vol. 50, no. 5, pp. 993–1006, 2012. View at: Publisher Site  Google Scholar
 V. A. Arkhipov, V. E. Zarko, I. K. Zharova et al., “Solid propellant combustion in a highvelocity crossflow of gases (review),” Combustion, Explosion, and Shock Waves, vol. 52, no. 5, pp. 497–513, 2016. View at: Publisher Site  Google Scholar
 J. M. Lenoir and G. Robillard, “A mathematical method to predict the effects of erosive burning in solidpropellant rockets,” Symposium (International) on Combustion, vol. 6, no. 1, pp. 663–667, 1957. View at: Publisher Site  Google Scholar
 M. K. King, “A modification of the composite propellant erosive burning model of Lenoir and Robillard,” Combustion and Flame, vol. 24, no. 24, pp. 365–368, 1975. View at: Publisher Site  Google Scholar
 L. Deverall, “The experimental and theoretical comparison of the erosive burning characteristics of composite propellants,” in Presented at the 3rd Solid Propulsion Conference, Atlantic City, NJ, U.S.A., 1968. View at: Google Scholar
 A. Peretz, K. K. Kuo, L. H. Caveny, and M. Summerfield, “Starting transient of solidpropellant rocket motors with high internal gas velocities,” AIAA Journal, vol. 11, no. 12, pp. 1719–1727, 1973. View at: Publisher Site  Google Scholar
 M. A. Willcox, M. Q. Brewster, K. C. Tang, D. S. Stewart, and I. Kuznetsov, “Solid rocket motor internal ballistics simulation using threedimensional grain burnback,” Journal of Propulsion and Power, vol. 23, no. 3, pp. 575–584, 2007. View at: Publisher Site  Google Scholar
 E. Cavallini, B. Favini, M. Di Giacinto, and F. Serraglia, “Internal ballistics simulation of a NAWC tactical SRM,” Journal of Applied Mechanics, vol. 78, no. 5, article 051018, 2011. View at: Publisher Site  Google Scholar
 A. K. Rettenmaier and S. D. Heister, “Experimental study of erosive and dynamic burning in polybutadienebased composite propellants,” Journal of Propulsion and Power, vol. 29, no. 1, pp. 87–94, 2013. View at: Publisher Site  Google Scholar
 H. S. Mukunda and P. J. Paul, “Universal behaviour in erosive burning of solid propellants,” Combustion and Flame, vol. 109, no. 1–2, pp. 224–236, 1997. View at: Publisher Site  Google Scholar
 H. S. Mukunda, P. J. Paul, A. Javed, and D. Chakraborty, “Extension of the universal erosive burning law to partly symmetric propellant grain geometries,” Acta Astronautica, vol. 93, pp. 176–181, 2014. View at: Publisher Site  Google Scholar
 A. Javed and D. Chakraborty, “Universal erosive burning model performance for solid rocket motor internal ballistics,” Aerospace Science and Technology, vol. 45, pp. 150–153, 2015. View at: Publisher Site  Google Scholar
 V. Gnielinski, “On heat transfer in tubes,” International Journal of Heat and Mass Transfer, vol. 63, pp. 134–140, 2013. View at: Publisher Site  Google Scholar
 D. R. Greatrix and J. J. Gottlieb, “Erosive burning model for compositepropellant rocket motors with large lengthtodiameter ratios,” Canadian Aeronautics and space journal, vol. 33, 1987. View at: Google Scholar
 D. R. Greatrix, “Regression rate estimation for standardflow hybrid rocket engines,” Aerospace Science and Technology, vol. 13, no. 7, pp. 358–363, 2009. View at: Publisher Site  Google Scholar
 C. F. Colebrook, T. Blench, H. Chatley, E. H. Essex, and J. R. Finniecome, “Turbulent flow in pipes, with particular reference to the transition region between the smooth and rough pipe LAWS,” Journal of the Institution of Civil Engineers, vol. 11, no. 4, pp. 133–156, 1939. View at: Publisher Site  Google Scholar
 S. E. Haaland, “Simple and explicit formulas for the friction factor in turbulent pipe flow,” Journal of Fluids Engineering, vol. 105, no. 1, pp. 8990, 1983. View at: Publisher Site  Google Scholar
 S. M. Ghiaasiaan, Convective Heat and Mass Transfer, Taylor & Francis, CRC Press, Boca Raton, 2018. View at: Publisher Site
 S. Gordon and B. J. McBride, Computer Program for Calculation of Complex Chemical Equilibrium Compositions and Applications. Part 1: Analysis, Publication NASARP1311, Cleveland, Ohio, 1994.
 M. K. King, Erosive Burning of Composite Propellants, Publication ADA033975. ATLANTIC RESEARCH CORP ALEXANDRIA, 1976.
 M. K. King, “Erosive burning of composite solid propellants: experimental and modeling studies,” Journal of Spacecraft and Rockets, vol. 16, no. 3, pp. 154–162, 1979. View at: Publisher Site  Google Scholar
 M. K. King, “Experimental and theoretical study of the effects of pressure and crossflow velocity on composite propellant burning rate,” Symposium (International) on Combustion, vol. 18, no. 1, pp. 207–216, 1981. View at: Publisher Site  Google Scholar
 T. Marklund and A. Lake, “Experimental investigation of propellant erosion,” ARS Journal, vol. 30, no. 2, pp. 173–178, 1960. View at: Publisher Site  Google Scholar
 H. Hasegawa, M. Hanzawa, S.I. Tokudome, and M. Kohno, “Erosive burning of aluminized composite propellants: Xray absorption measurement, correlation, and application,” Journal of Propulsion and Power, vol. 22, no. 5, pp. 975–983, 2006. View at: Publisher Site  Google Scholar
 S. V. Patankar and D. B. Spalding, “A calculation procedure for heat, mass and momentum transfer in threedimensional parabolic flows,” in Numerical Prediction of Flow, Heat Transfer, Turbulence and Combustion, pp. 54–73, Elsevier, 1983. View at: Publisher Site  Google Scholar
 F. Moukalled, L. Mangani, and M. Darwish, The Finite Volume Method in Computational Fluid Dynamics: An Advanced Introduction with Open FOAM® and Matlab, Springer International Publishing, Cham, 2016. View at: Publisher Site
Copyright
Copyright © 2020 Yanjie Ma et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.