Next Article in Journal
An Intelligent Data-Driven Approach for Electrical Energy Load Management Using Machine Learning Algorithms
Previous Article in Journal
Analysis of Energy Loss on a Tunable Check Valve through the Numerical Simulation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Measurement of Gas Flow Rate at Gasification of Low-Melting Materials in a Flow-Through Gas Generator

by
Dmitry A. Vnuchkov
1,
Valery I. Zvegintsev
1,
Denis G. Nalivaichenko
1 and
Sergey M. Frolov
2,*
1
Khristianovich Institute of Theoretical and Applied Mechanics of the Siberian Branch of the Russian Academy of Sciences (ITAM SB RAS), 630090 Novosibirsk, Russia
2
Department of Combustion and Explosion, Semenov Federal Research Center for Chemical Physics of the Russian Academy of Sciences, 119991 Moscow, Russia
*
Author to whom correspondence should be addressed.
Energies 2022, 15(15), 5741; https://doi.org/10.3390/en15155741
Submission received: 6 July 2022 / Revised: 1 August 2022 / Accepted: 5 August 2022 / Published: 8 August 2022

Abstract

:
A semi-empirical method is proposed for determining the rate of gas production in a flow-through gas generator (GG) with the allocation of a part of the gas flow produced by gasification of a low-melting solid material (LSM) in the total gas flow rate through the GG. The method is verified by test fires with polypropylene sample gasification by hot air under conditions of incoming supersonic flow with Mach number 2.43, 2.94, and 3.81 and stagnation temperature 600–700 K. The mean flow rates of gasification products obtained in test fires were 0.08 kg/s at Mach 2.43, 0.10 kg/s at Mach 2.94, and 0.05–0.02 kg/s at Mach 3.81. For obtaining 1 kg of gasification products in the test fires there was a need of 1.61 to 2.86 kg of gasifying agent.

1. Introduction

Gas generators (GG) are designed for generating gases of desired composition, flow rate, temperature, and pressure through different physical and chemical processes such as liquid/solid charge vaporization, pyrolysis, gasification, combustion, etc. All GGs can be divided into two groups: closed-type GGs and flow-through GGs. In the former, a working medium is initially enclosed in a single housing, so such GGs are mostly single-action devices. In the latter, a working medium is continuously supplied to GG, which ensures long-term operation. The closed-type GGs include GGs deploying airbags in automobiles, injecting inhibitors in accidental fires, driving turbopumps in rocket engines, etc. [1,2,3,4,5,6,7,8].
The flow-through GGs are used in industry for production of various chemicals and energy, e.g., by converting crude oil, biomass, sewage sludge, municipal solid wastes, etc. to syngas, which is further used for production of hydrogen, synthetic natural gas, ammonia, methanol, Fisher–Tropsch diesel, and/or generating heat and electricity [9,10,11,12]. Given the widespread use of GGs, the methods for determining their characteristics at both design and operation stages is relevant.
For this purpose, methods of numerical simulation of physicochemical processes in GGs are now widely used [13,14,15,16,17,18]. However, at present the complexity of these processes does not allow predicting the main characteristic of GGs, i.e., the instantaneous rate of gas production, with sufficient accuracy. In view of this, experimental verification and refinement of numerical simulations remain important [19,20]. A good illustrative example is the GG of a hybrid rocket. Contrary to solid rockets, the direct calculation of the mass flow rate from chamber pressure and nozzle throat area is not possible for hybrid rockets because the characteristic exhaust velocity c* strongly depends on the oxidizer-to-fuel ratio. According to [21], the measurement of the instantaneous fuel mass flow rate for hybrid rockets is a great challenge, because it is a function of operation conditions, firing duration, port diameter, nozzle erosion conditions, etc. Therefore, the end-point averaging method based on the initial and final shapes of fuel charge is often used for estimating the average fuel mass flow rate. However, the fuel mass flow rate is usually not constant during motor firings. The approach based on short firing duration for the end-point averaging is also not reliable due to the uncertainties introduced by ignition and shutdown transients [22]. Several reconstruction techniques for estimating the instantaneous fuel mass flow rate based on the measured chamber pressure and oxidizer mass flow rate have been reported in the literature on hybrid rockets [23,24,25,26,27]. However, all these techniques imply a known value of the oxidizer mass flow rate, whereas in the flow-through GGs, the latter can be unknown, especially when the GG is placed in the free approaching air stream. In the latter case, the combustion/gasification process in the GG can affect the mass flow rate of air at the GG intake.
This manuscript presents a method for the semi-empirical determination of the instantaneous rate of gas production in a flow-through GG placed in a free air stream with the allocation of a part of the gas flow produced by gasification of a low-melting solid material (LSM) in the total gas flow rate at the GG exit. This method and its demonstration are the novel and distinctive features of the present work.

2. Experimental Setup

Figure 1 shows the general view of a free jet facility for testing a GG in supersonic flow of gasifying agent (air). The facility was based on the Model Aerodynamic Facility (MAF) at the ITAM SB RAS [28]. Compressed air from high-pressure cylinders passed through a prechamber and entered a fire heater operating on combustion of hydrogen to obtain a hot air. Mass flow rate of hydrogen to the fire heater varied from 2 to 4 g/s. Apart from air and hydrogen, oxygen was added to the air flow in an amount of 20 to 40 g/s for keeping the mass fraction of oxygen at 23%wt in the vitiated air outflow. Thus, the outflowing gasifying agent contained oxygen (23%), nitrogen (75%), and steam (2%). The change in thermophysical properties of air due to steam admixture was neglected.
Three profiled nozzles with an exit diameter of 30 mm designed for flow Mach numbers M0 = 2.5, 3.0 and 4.0 could be attached to the fire heater. Table 1 shows the calculated and measured parameters of the corresponding jet flow. Measurements of the real values of the flow Mach number M1 were performed using a Pitot pressure sensor (Pitot pressure P0′) installed at the nozzle exit. In each experiment, the total pressure P0(t) and stagnation temperature T0(t) were recorded for calculating the mass flow rate of gases, G0(t). A jet of hot gasifying agent from the facility nozzle partly entered the GG mounted downstream along the facility axis.

3. Gas Generator Design

Figure 2 shows the schematic of the GG. An axisymmetric cylindrical intake with a diameter of 15 mm and length of 46.6 mm was installed in the frontal part of GG. From the cylindrical intake, hot gasifying agent flowed through a conical diffuser to enter the LSM sample. In the diffuser, a total pressure sensor and thermocouple were installed, which measured the stagnation pressure P0,in and stagnation temperature T0,in in the gas flow at the inlet to the LSM sample. Around the diffuser, a 36-g pyrotechnic igniter with a calorific value of about 1.5 MJ/kg was placed. It was used to ignite combustion of LSM aimed at providing heat for endothermic physical and chemical processes accompanying sample gasification. The igniter was triggered at a preset time by a command of the synchronization system. Typical igniter burning time was 0.10–0.15 s. Preliminary estimates showed that the gas temperature at the inlet of the LSM sample could theoretically reach 1300–2600 K, however the actual measured gas temperature was on the level of 1300–1500 K, thus indicating large heat loss into the walls.
The LSM sample was assembled from 16 identical polypropylene (PP) blocks 80 mm in diameter and 23 mm long (Figure 3). Each block had 61 orifices with a diameter of 4 mm. When assembling the sample, the orifices of all blocks were aligned to form straight longitudinal channels. The total initial sample mass was 1.28 kg. The masses of each block and sample in assembly were measured by balances both before and after test fires with an error of 0.05 g.
A measuring nozzle with a throat of 13 or 20 mm in diameter was mounted at the GG outlet downstream of the LSM sample. In the nozzle, three thermocouples and total pressure sensor were installed to measure the stagnation pressure P0,out and stagnation temperature T0,out. To measure static pressure Pout and control the flow Mach number, a static pressure sensor was also mounted in the nozzle throat. After exiting the measuring nozzle, the exhaust jet was directed to the noise absorbing ventilation shaft at atmospheric pressure. For pressure measurements, pressure sensors RPD-I with a maximum pressure up to 1.0 and 2.5 MPa were used. The rated measurement error of the pressure sensors was 0.2%. The actual pressure measurement error was established by the results of numerous calibrations of the sensors in the expected pressure range. In this case, to control the set pressure, a PDE-020I reference pressure transducer was used with a basic relative error of ±0.02%. The resulting estimate of the actual pressure measurement error did not exceed 1%. For temperature measurements, tungsten-rhenium thermocouples were used. To convert the electrical signal of thermocouples into temperature readings, a standard calibration table was used. The temperature measurement error did not exceed 5%.
The data acquisition system was based on the National Instruments NI PCI-6255 board. The system had 80 differential measuring channels. The bit depth of the analog-to-digital converter (ADC) was 16 bits. The maximum sampling frequency was 1.25 MHz. In the experiments, a sampling frequency of 1000 samples for each channel was used.
The frequency characteristics of the measuring equipment including pressure and temperature sensors as well as the data acquisition system were measured by applying preset pulse signals. It was shown based on the results of such dynamic tests that the operation frequency of the measuring equipment used ranged from 0 to 100 Hz. Since the characteristic frequencies of the processes under investigation did not exceed 10 Hz, dynamic measurement errors were not considered and were not taken into account.

4. Test Results

All test fires were conducted with igniter triggering. However, combustion of LSM sample failed in some test fires. Figure 4, Figure 5 and Figure 6 show the typical pressure and temperature records in test fires without and with combustion of LSM samples. In test fires without combustion (tests 1, 4, and 7 in Figure 4, Figure 5 and Figure 6), after ignition triggering, all measured parameters returned to their original values. In this case, the measured sample mass before and after test fires showed virtually no LSM depletion. In test fires with combustion (tests 2, 3, 5, 6, 8, 9 in Figure 4, Figure 5 and Figure 6), the pressure at the LSM sample inlet and outlet after ignition triggering sharply increased by a factor of 1.5 to 2.0. The temperature at the sample outlet increased to 800–1000 K, although the temperature at the sample inlet corresponded to the temperature of the incoming gas flow and remained constant at 600–700 K. This indicated the existence of an exothermic oxidation process due to partial interaction of PP decomposition products with oxygen in the incoming flow of gasifying agent.

5. Processing of Test Results

5.1. Calculation of Gas Flow Rate at GG Intake

The flow rate of gasifying agent at the GG intake was calculated based on the stagnation parameters P0(t) and T0(t) and the approaching flow Mach number M1:
Gin(t) = φFinρ(t)V(t) = φFinM1P0(t)π(M1)[γRT0(t)τ(M1)]1/2[RT0(t)τ(M1)]−1 =
        φFinM1π(M1)τ(M1)−1/2(γ/R)P0(t)T0(t)−1/2
where t is time, ρ is gas density, V is gas velocity, Fin is the cross-sectional area of the GG intake, φ is the contraction ratio of GG intake, R = 287 kJ/kg/K is the gas constant for air, γ = 1.4 is the specific heat ratio for air, π(M1) and τ(M1) are the gas-dynamic functions.

5.2. Contraction Ratio of GG Intake

In test fires without combustion, a cylindrical GG intake cut out a stream of a cross-sectional area Fin with the Mach number M1 from the free jet emanating from the facility nozzle. In this case, the contraction ratio of the GG intake was φ = 1.0. Let us jointly consider Figure 2 and curves 2 (P0,in) in tests 1, 4, and 7 in Figure 4, Figure 5 and Figure 6. Downstream of the cylindrical section of the intake, the supersonic gas flow entered the expanding conical diffuser and accelerated to Mach number M2 (M2 > M1). Thereafter, due to geometric throttling in the LSM channels, the flow passed through a normal shock and became subsonic with a sharp decrease in the stagnation pressure to P0,in(M2).
A similar flow pattern was observed in test fires with combustion prior to ignition triggering. Figure 7 shows several examples of such test fires in terms of the measured time histories of P0,in/P0 ratio. Prior to ignition triggering, the P0,in/P0 ratio corresponded to the total pressure loss at elevated Mach number M2 > M1. This indicated that the normal shock was located in the conical diffuser upstream of the total pressure sensor. After ignition triggering and establishment of LSM sample combustion, thermal throttling of the flow caused the normal shock to move upstream and exit the GG intake. Starting from this time instant, the P0,in/P0 ratio corresponded to the total pressure loss in the normal shock at the freestream Mach number M1. A decrease in the Mach number ahead of the normal shock led, on the one hand, to a pressure rise in the flow entering the channels of LSM sample and thereby intensified the combustion and heat and mass transfer processes in sample channels. On the other hand, the head shock wave detached from the GG intake led to a decrease in the gas flow rate through the intake, i.e., to a decrease in the contraction ratio φ in Equation (1). This posed a problem of determining the realistic value of the contraction ratio φ < 1.0. The methodology for solving this problem by measuring the gas flow rate at the GG outlet is presented below. Once the contraction ratio is determined, Equation (1) can be used for calculating the mass flow rate of gasifying agent entering the GG intake.

5.3. Calculation of Gas Flow Rate at GG Outlet

When the hot gasifying agent flows through the channels in the LSM sample, the LSM is heated, melted, pyrolyzed, and gasified. The produced gasification products are mixed and chemically react with the gasifying agent, thus leading to sample combustion accompanied with heat release. The heat of combustion promotes further sample gasification, and the gas flow rate increases. The task is to determine the gas flow rate at GG outlet.
In the test fires, this gas flow rate was determined by mounting a sonic nozzle of cross-sectional area F* (asterisk * denotes sonic nozzle throat) at the GG outlet and measuring the values of stagnation pressure P0,out and stagnation temperature T0,out before the nozzle. The values of F*, P0,out, and T0,out were used to calculate Gout using the relation:
Gout(t) = mF*P0,out(t)T0,out(t)−1/2
where the dimensional coefficient m is a function of gas composition at the nozzle:
m = [(γ*/R*)(2/(γ* + 1))(γ* + 1)/(γ − 1)]1/2
The main difficulty in calculating Gout is then the determination of coefficient m, as the gas composition depends on the LSM combustion completeness, as well as local instantaneous pressure and temperature. Herein, the values of γ* and R* for the gas mixture were determined by iterations using the Astra 4 thermodynamic code [29].
As an example, Figure 8 shows the results of thermodynamic calculations for mixtures composed of PP and air with the PP-to-air mass ratio ranging from 0 to 50% at a constant pressure of 1 MPa. Note that the results of calculations in a pressure range between 0.1 and 1 MPa differ by less than 3% only.

5.4. Flow Rate of Polypropylene Gasification Products

Figure 9, Figure 10 and Figure 11 show the time histories of flow rates Gin and Gout calculated based on Equations (1) and (2), respectively, for different freestream Mach numbers M1. In addition to plots Gin(t) and Gout(t), the calculated difference between these flow rates, ΔG = Gout(t) − Gin(t) is also plotted in Figure 9 to Figure 11. Clearly, this difference corresponds to the total flow rate of product gases generated by LSM sample gasification.
As mentioned above, LSM combustion in tests 1, 4, and 7 was absent. Therefore, the flow rates Gin and Gout were virtually equal. In test fires with combustion, a noticeable LSM gasification accompanied with sample combustion was detected. In these latter test fires, the gas flow rates Gin and Gout differed markedly. According to Figure 9, Figure 10 and Figure 11 the flow rate of gasification products (ΔG = Gout(t) − Gin(t)) in test fires with combustion could be constant, like in Figure 9 and Figure 10, or variable in time, like in Figure 11. The mean values of the flow rates of gasification products in Figure 9, Figure 10 and Figure 11 took the values of 0.08 kg/s at M1 = 2.43, 0.10 kg/s at M1 = 2.94, and 0.05–0.02 kg/s at M1 = 3.81.
An important element of test data processing procedure was the determination of LSM sample mass before (W1) and after (W2) test fire. The difference between these masses (ΔW = W1W2) determines the amount of gasified sample material leaving the GG with exhaust gases. Based on the calculated flow rates Gin and Gout and measured sample masses W1 and W2, one can derive the following balance equation for determining the intake contraction ratio φ:
W1W2 = mF* t1ʃt2 P0,out(τ)T0,out(τ)−1/2dτ − φFinM1π(M1)τ(M1)−1/2(γ/R) t1ʃt2 P0(τ)T0(τ)−1/2
In Equation (4), time t1 corresponds to the start of gas temperature rise after igniter triggering and time t2 corresponds to the shutdown of air supply. An unknown value of the air intake contraction ratio φ can be now determined from Equation (4). Once φ is determined, the realistic value of gasifying agent flow rate at the inlet to LSM sample can be obtained. Furthermore, the value of coefficient m in Equation (3) can be refined.
Table 2 shows the test results thus obtained and refined. Now, the time integrals of functions ΔG(t) shown in Figure 9, Figure 10 and Figure 11 are equal to the corresponding values of ΔW in each test fire. As seen from Table 2, the ratio of total amounts of gasifying agent to gasification products, ʃGindt/ʃΔGdt, varied from 1.61 to 2.86 in the test fires. This means that for obtaining 1 kg of gasification products one consumes 1.61 to 2.86 kg of gasifying agent.

6. Amendment

The adopted technology of dividing the LSM sample into separate blocks and measuring their individual masses Wi before and after test fires allowed obtaining additional useful information on the zones of maximum LSM decomposition along the sample length. Figure 12 shows the distribution of masses Wi of individual blocks in a test sample, measured before and after test 3. The individual blocks are numbered from the inlet of the LSM sample. As seen, the zone of most intense LSM decomposition in this test was located between blocks 4 and 9. Figure 13 shows the photographs of all individual blocks of the LSM sample before and after test 3.

7. Conclusions

A novel semi-empirical method is proposed for evaluating the instantaneous mass flow rate of gasification products in the flow-through gas generators operating on air-assisted gasification of low-melting solid materials. The instantaneous mass flow rate of gasification products is derived based on determining the mass flow rate of air captured by the gas generator intake in a supersonic flow and the instantaneous mass flow rate of gases at the gas generator outlet, supplemented with measurements of sample mass before and after test firing.
To verify the methodology, a set of test fires with hot air as gasifying agent and polypropylene as a low-melting solid material were conducted under conditions of incident supersonic flow of gasifying agent with Mach numbers 2.43, 2.94, and 3.81 and stagnation temperatures 600–700 K. The mean flow rates of gasification products obtained experimentally were 0.08 kg/s at Mach 2.43, 0.10 kg/s at Mach 2.94, and 0.05–0.02 kg/s at Mach 3.81. In the test fires, the ratio of total amounts of air to gasification products varied from 1.61 to 2.86 (mass basis). This means that for obtaining 1 kg of gasification products one needs 1.61 to 2.86 kg of air.

Author Contributions

V.I.Z.: Conceptualization, methodology; D.A.V.: investigation; D.G.N.: investigation; S.M.F.: writing original draft and editing; project administration. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Ministry of Science and Higher Education of Russian Federation under state contract N13.1902.21.0014 (agreement N075-15-2020-806).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations and Nomenclature

ADCAnalog-to-digital converter
GGGas generator
LSMLow-melting solid material
MAFModel Aerodynamic Facility
ITAM SB RASKhristianovich Institute of Theoretical and Applied Mechanics of the Siberian Branch of the Russian Academy of Sciences
 
c*Characteristic exhaust velocitym/s
M0Design nozzle Mach number-
M1Air flow Mach number at the nozzle exit-
M2Mach number before the normal shock in diffuser-
P0Pitot pressure at the nozzle exitPa
P0Total pressure of the air flowPa
T0Total (stagnation) temperature of the flowK
G0Mass flow ratekg/s
P0,inTotal pressure at the entrance of LSM samplePa
T0,inTotal temperature at the entrance of LSM sampleK
P0,outTotal pressure at the exit of LSM samplePa
T0,outTotal temperature at the exit of LSM sampleK
PoutStatic pressure at the sound nozzlePa
FinArea of intake entrance cross-sectionm2
F*Area of sonic nozzle throat cross-sectionm2
GinAir flow rate through intakekg/s
GoutGas flow rate at the sonic nozzle throatkg/s
φContraction ratio of intake-
RGas constant for airJ/kg/K
R*Gas constant for gas mixtureJ/kg/K
γSpecific heat ratio for air-
γ*Specific heat ratio for gas mixture-
π(M1)Gas-dynamic function-
τ(M1)Gas-dynamic function-
mDimensional coefficient-
W1Sample mass before testkg
W2Sample mass after testkg
WiMasses of individual blocks in the LSM sampleg

References

  1. Schmitt, R.G.; Butler, P.B.; Freesmeier, J.J. Performance and CO Production of a Non-Azide Airbag Propellant in a Pre-Pressurized Gas Generator. Combust. Sci. Technol. 1997, 122, 305–330. [Google Scholar] [CrossRef]
  2. Freesmeier, J.J.; Butler, P.B. Analysis of a Hybrid Dual-Combustion-Chamber Solid-Propellant Gas Generator. J. Propuls. Power 1999, 15, 552–561. [Google Scholar] [CrossRef]
  3. Wu, W.T.; Hsieh, W.H.; Huang, C.H.; Wang, C.H. Theoretical Simulation of Combustion and Inflation Processes of Two-Stage Airbag Inflators. Combust. Sci. Technol. 2005, 177, 383–412. [Google Scholar] [CrossRef]
  4. Galfetti, L.; Merotto, L.; Boiocchi, M.; Maggi, F.; DeLuca, L.T. Experimental investigation of paraffin-based fuels for hybrid rocket propulsion. Prog. Propuls. Phys. 2013, 4, 59–74. [Google Scholar] [CrossRef] [Green Version]
  5. Wada, Y.; Kawabata, Y.; Kato, R.; Kato, N.; Hori, K. Observation of combustion behavior of low melting temperature fuel for a hybrid rocket using double slab motor. Int. J. Energetic Mater. Chem. Propuls. 2016, 15, 351–369. [Google Scholar] [CrossRef]
  6. Lee, D.; Lee, C. Hybrid Gas Generator for a Staged Hybrid Rocket Engine. J. Propuls. Power 2017, 33, 204–212. [Google Scholar] [CrossRef]
  7. Cheng, C.; Zhang, X.; Wang, C.; Wang, L. Numerical investigation on cooling performance of filter in a pyrotechnic gas generator. Def. Technol. 2020, 17, 343–351. [Google Scholar] [CrossRef]
  8. Frolov, S.; Shamshin, I.; Kazachenko, M.; Aksenov, V.; Bilera, I.; Ivanov, V.; Zvegintsev, V. Polyethylene Pyrolysis Products: Their Detonability in Air and Applicability to Solid-Fuel Detonation Ramjets. Energies 2021, 14, 820. [Google Scholar] [CrossRef]
  9. Mahinpey, N.; Gomez, A. Review of gasification fundamentals and new findings: Reactors, feedstock, and kinetic studies. Chem. Eng. Sci. 2016, 148, 14–31. [Google Scholar] [CrossRef]
  10. Zhang, Y.; Xu, P.; Liang, S.; Liu, B.; Shuai, Y.; Li, B. Exergy analysis of hydrogen production from steam gasification of biomass: A review. Int. J. Hydrogen Energy 2019, 44, 14290–14302. [Google Scholar] [CrossRef]
  11. Indrawan, N.; Kumar, A.; Moliere, M.; Sallam, K.A.; Huhnke, R.L. Distributed power generation via gasification of biomass and municipal solid waste: A review. J. Energy Inst. 2020, 93, 2293–2313. [Google Scholar] [CrossRef]
  12. Srivastava, D. Advances in biofuel technology: A review. IOSR J. Eng. 2021, 11, 32–41. [Google Scholar]
  13. Ku, X.; Li, T.; Løvås, T. CFD–DEM simulation of biomass gasification with steam in a fluidized bed reactor. Chem. Eng. Sci. 2015, 122, 270–283. [Google Scholar] [CrossRef] [Green Version]
  14. Liu, H.; Cattolica, R.J.; Seiser, R. CFD studies on biomass gasification in a pilot-scale dual fluidized-bed system. Int. J. Hydrogen Energy 2016, 41, 11974–11989. [Google Scholar] [CrossRef]
  15. Qi, T.; Lei, T.; Yan, B.; Chen, G.; Li, Z.; Fatehi, H.; Wang, Z.; Bai, X.-S. Biomass steam gasification in bubbling fluidized bed for higher-H2 syngas: CFD simulation with coarse grain model. Int. J. Hydrogen Energy 2019, 44, 6448–6460. [Google Scholar] [CrossRef]
  16. Yang, S.; Fan, F.; Wei, Y.; Hu, J.; Wang, H.; Wu, S. Three-dimensional MP-PIC simulation of the steam gasification of biomass in a spouted bed gasifier. Energy Convers. Manag. 2020, 210, 112689. [Google Scholar] [CrossRef]
  17. Rashkovskiy, S.; Yakush, S. Numerical simulation of low-melting temperature solid fuel regression in hybrid rocket engines. Acta Astronaut. 2020, 176, 710–716. [Google Scholar] [CrossRef]
  18. Salgansky, E.; Lutsenko, N. Effect of solid fuel characteristics on operating conditions of low-temperature gas generator for high-speed flying vehicle. Aerosp. Sci. Technol. 2020, 109, 106420. [Google Scholar] [CrossRef]
  19. Shiplyuk, A.; Zvegintsev, V.; Frolov, S.; Vnuchkov, D.; Kiseleva, T.; Kislovsky, V.; Lukashevich, S.; Melnikov, A.Y.; Nalivaychenko, D. Gasification of low-melting hydrocarbon material in the airflow heated by hydrogen combustion. Int. J. Hydrogen Energy 2020, 45, 9098–9112. [Google Scholar] [CrossRef]
  20. Shiplyuk, A.N.; Zvegintsev, V.I.; Frolov, S.M.; Vnuchkov, D.A.; Kislovsky, V.A.; Kiseleva, T.A.; Lukashevich, S.V.; Melnikov, A.Y.; Nalivaychenko, D.G. Gasification of Low-Melting Fuel in a High-Temperature Flow of Inert Gas. J. Propuls. Power 2021, 37, 20–28. [Google Scholar] [CrossRef]
  21. Kuo, K.K.; Chiaverini, M.J. Challenges of Hybrid Rocket Propulsion in the 21st Century. In Fundamentals of Hybrid Rocket Combustion and Propulsion; Chiaverini, M.J., Kuo, K.K., Eds.; Progress in astronautics and aeronautics; American Institute of Aeronautics and Astronautics: Reston, VA, USA, 2007; pp. 593–638. [Google Scholar]
  22. Frederick, R.A.; Greiner, B.E. Laboratory-scale hybrid rocket motor uncertainty analysis. J. Propuls. Power 1996, 12, 605–611. [Google Scholar] [CrossRef]
  23. Wernimont, E.J.; Heister, S.D. Reconstruction Technique for Reducing Hybrid-Rocket Combustion Test Data. J. Propuls. Power 1999, 15, 128–136. [Google Scholar] [CrossRef]
  24. George, P.; Krishnan, S.; Varkey, P.M.; Ravindran, M.; Ramachandran, L. Fuel regression rate in hydrox-yl-terminated-polybutadiene/gaseous-oxygen hybrid rocket motors. Propuls. Power 2001, 17, 35–42. [Google Scholar] [CrossRef]
  25. Nagata, H.; Nakayama, H.; Watanabe, M.; Wakita, M.; Totani, T. Accuracy and applicable range of a reconstruction technique for hybrid rockets. Adv. Aircr. Spacecr. Sci. 2014, 1, 273–289. [Google Scholar] [CrossRef] [Green Version]
  26. Carmicino, C.; Sorge, A.R. Role of injection in hybrid rockets regression rate behavior. Propuls. Power 2005, 21, 606–612. [Google Scholar] [CrossRef]
  27. Saito, Y.; Kamps, L.T.; Komizu, K.; Bianchi, D.; Nasuti, F.; Nagata, H. The accuracy of reconstruction techniques for determining hybrid rocket fuel regression rate. In Proceedings of the Joint Propulsion Conference AIAA Paper 2018–4923, AIAA Propulsion and Energy Forum, Cincinnati, OH, USA, 9–11 July 2018. [Google Scholar] [CrossRef]
  28. Zvegintsev, V.I. Short-Duration Gas-dynamic Facilities. In Part 1. Facilities for Scientific Research; Parallel Publishing: Novosibirsk, Russia, 2014. (In Russian) [Google Scholar]
  29. Trusov, B.G. Modeling Chemical and Phase Equilibria at High Temperatures “Astra 4.”; Bauman Moscow State Technical University: Moscow, Russia, 1991. (In Russian) [Google Scholar]
Figure 1. General view of the free jet test facility with GG.
Figure 1. General view of the free jet test facility with GG.
Energies 15 05741 g001
Figure 2. Schematic of flow-through GG.
Figure 2. Schematic of flow-through GG.
Energies 15 05741 g002
Figure 3. One of 16 identical polyethylene blocks for assembling a test sample; dimensions are given in millimeters.
Figure 3. One of 16 identical polyethylene blocks for assembling a test sample; dimensions are given in millimeters.
Energies 15 05741 g003
Figure 4. Test results at M1 = 2.43 (see Table 1): 1—P0; 2—P0,in; 3—P0,out; 4—Pout; 5—T0; 6—T0,in; 7—T0,out.
Figure 4. Test results at M1 = 2.43 (see Table 1): 1—P0; 2—P0,in; 3—P0,out; 4—Pout; 5—T0; 6—T0,in; 7—T0,out.
Energies 15 05741 g004
Figure 5. Test results at M1 = 2.94 (see Table 1): 1—P0; 2—P0,in; 3—P0,out; 4—Pout; 5—T0; 6—T0,in; 7—T0,out.
Figure 5. Test results at M1 = 2.94 (see Table 1): 1—P0; 2—P0,in; 3—P0,out; 4—Pout; 5—T0; 6—T0,in; 7—T0,out.
Energies 15 05741 g005
Figure 6. Test results at M1 = 3.81 (see Table 1): 1—P0; 2—P0,in; 3—P0,out; 4—Pout; 5—T0; 6—T0,in; 7—T0,out.
Figure 6. Test results at M1 = 3.81 (see Table 1): 1—P0; 2—P0,in; 3—P0,out; 4—Pout; 5—T0; 6—T0,in; 7—T0,out.
Energies 15 05741 g006
Figure 7. Formation of a diverted head shock at the GG intake in test fires with LSM sample combustion.
Figure 7. Formation of a diverted head shock at the GG intake in test fires with LSM sample combustion.
Energies 15 05741 g007
Figure 8. Calculated dependences of coefficient m on temperature for PP–air mixtures with different content of PP (mass basis) at pressure 1 MPa.
Figure 8. Calculated dependences of coefficient m on temperature for PP–air mixtures with different content of PP (mass basis) at pressure 1 MPa.
Energies 15 05741 g008
Figure 9. Calculated time histories of flow rates Gin (curve 8), Gout (9) and ΔG = Gout(t) − Gin(t) (10) at M1 = 2.43 in tests 1 to 3.
Figure 9. Calculated time histories of flow rates Gin (curve 8), Gout (9) and ΔG = Gout(t) − Gin(t) (10) at M1 = 2.43 in tests 1 to 3.
Energies 15 05741 g009
Figure 10. Calculated time histories of flow rates Gin (curve 8), Gout (9) and ΔG = Gout(t) − Gin(t) (10) at M1 = 2.94 in tests 4 to 6.
Figure 10. Calculated time histories of flow rates Gin (curve 8), Gout (9) and ΔG = Gout(t) − Gin(t) (10) at M1 = 2.94 in tests 4 to 6.
Energies 15 05741 g010
Figure 11. Calculated time histories of flow rates Gin (curve 8), Gout (9) and ΔG = Gout(t) − Gin(t) (10) at M1 = 3.81 in tests 7 to 9.
Figure 11. Calculated time histories of flow rates Gin (curve 8), Gout (9) and ΔG = Gout(t) − Gin(t) (10) at M1 = 3.81 in tests 7 to 9.
Energies 15 05741 g011
Figure 12. Masses of blocks composing LSM sample before (1) and after (2) test 3 with combustion. Burning time t2t1 = 3.46 – 1.15 = 2.31 s.
Figure 12. Masses of blocks composing LSM sample before (1) and after (2) test 3 with combustion. Burning time t2t1 = 3.46 – 1.15 = 2.31 s.
Energies 15 05741 g012
Figure 13. Photographs of (a) upstream and (b) downstream faces of blocks composing LSM sample after test 3 with combustion.
Figure 13. Photographs of (a) upstream and (b) downstream faces of blocks composing LSM sample after test 3 with combustion.
Energies 15 05741 g013
Table 1. Characteristics of the jet flow at the entrance of GG intake.
Table 1. Characteristics of the jet flow at the entrance of GG intake.
Designed
M0
Mean P0′/P0Measured
M1
P0, MPa T0, KG0, kg/s
2.50.5282.431.906650.85
3.00.3462.943.208300.79
4.00.1633.813.008500.32
Table 2. Results of tests with combustion.
Table 2. Results of tests with combustion.
TestφʃGindt,
kg
mΔW,
kg
ʃΔGdt,
kg
ʃGindt/ʃΔGdt
20.7510.3300.03060.1870.1861.77
30.7280.3260.03080.1760.1741.86
50.7740.4440.03110.2760.2731.61
60.8100.3720.03120.2190.2161.70
80.6620.1330.03440.0640.0642.10
90.7170.1440.03440.0500.0502.86
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vnuchkov, D.A.; Zvegintsev, V.I.; Nalivaichenko, D.G.; Frolov, S.M. Measurement of Gas Flow Rate at Gasification of Low-Melting Materials in a Flow-Through Gas Generator. Energies 2022, 15, 5741. https://doi.org/10.3390/en15155741

AMA Style

Vnuchkov DA, Zvegintsev VI, Nalivaichenko DG, Frolov SM. Measurement of Gas Flow Rate at Gasification of Low-Melting Materials in a Flow-Through Gas Generator. Energies. 2022; 15(15):5741. https://doi.org/10.3390/en15155741

Chicago/Turabian Style

Vnuchkov, Dmitry A., Valery I. Zvegintsev, Denis G. Nalivaichenko, and Sergey M. Frolov. 2022. "Measurement of Gas Flow Rate at Gasification of Low-Melting Materials in a Flow-Through Gas Generator" Energies 15, no. 15: 5741. https://doi.org/10.3390/en15155741

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop