Next Article in Journal
Distributed Adaptive Path-Following Control for Distance-Based Formation of Fixed-Wing UAVs under Input Saturation
Next Article in Special Issue
Using the Impulse Method to Determine High-Pressure Dynamic Burning Rate of Solid Propellants
Previous Article in Journal
Numerical Investigation of Hypersonic Flat-Plate Boundary Layer Transition Subjected to Bi-Frequency Synthetic Jet
Previous Article in Special Issue
Numerical Investigation on the Effect of Ammonium Perchlorate Content and Position on the Combustion Characteristics of an Ammonium Perchlorate/Hydroxyl-Terminated Polybutadiene Propellant
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Phenomenological Model for the Unsteady Combustion of Solid Propellants from a Zel’dovich-Novzhilov Approach

Science and Technology on Combustion, Internal Flow and Thermo-Structure Laboratory, Northwestern Polytechnical University, Xi’an 710072, China
*
Author to whom correspondence should be addressed.
Aerospace 2023, 10(9), 767; https://doi.org/10.3390/aerospace10090767
Submission received: 20 July 2023 / Revised: 25 August 2023 / Accepted: 28 August 2023 / Published: 29 August 2023
(This article belongs to the Special Issue Combustion Evaluation and Control of Solid Rocket Motors)

Abstract

:
Solid rocket motors are prone to combustion instabilities, which may lead to various problems for the rockets, from unexpected oscillations, precision decreasing, to explosion. The unsteady combustion dynamics of the propellants play a crucial role in most solid rocket motors experiencing combustion instabilities. A modeling framework for the unsteady combustion of the solid propellant is constructed via the Zel’dovich-Novozhilov (ZN) phenomenological perspective. The overall unsteady combustion features of a quasi-steady homogeneous one-dimensional (QSHOD) model are investigated. The phenomenological ZN parameters are then calculated. Compared with the traditional ZN-QSHOD linear equivalence relation, the new calculated system yields better results for the pressure coupling response, especially in the non-linear regime. The proposed phenomenological modeling provides a new methodology for the model reduction of the complex flame models.

1. Introduction

With the development of modern propulsion technology, both the energy efficiency and the energy density reach a considerably high level. Combustion instability has been an unsolved problem since the early application of both liquid rocket engines and solid rocket motors (SRMs) [1,2].
The sources of the combustion instabilities in SRMs are believed to be the interaction of the flow field, the acoustic fields and the unsteady combustion inside the chamber [1]. The coupling between the unsteady combustion and the acoustic fields, or the thermoacoustic mechanism, plays a central role in various SRMs, especially in small and medium SRMs [3]. While for the large scale SRMs with high length-to-diameter ratio, the hydrodynamical instability governs.
The thermoacoustic instability in SRMs has been attracting most of the researchers in the community since the 1960s. According to Culick [4], a reliable modeling of the thermoacoustic instability should properly include the non-linear effects in the acoustic system of the chamber and the combustion system of the propellants. Actually, the three combinations of linear acoustics with linear combustion response, the linear acoustics with non-linear combustion response and the non-linear acoustics with linear combustion response, are not able to model the important phenomenon of bi-stability and triggering. Following Culick’s approach [5,6], the non-linear acoustics in the chamber can be well-modeled. Nevertheless, on the other side, the modeling of the unsteady combustion remains difficult until now. Traditional approaches are based on different overall response functions [4]. The major limit of this approach is the historical effects of temperature distribution inside the solid propellants. On account of this, Mariappan and Sujith [7] improved the model by replacing the overall response function by a spatial position specified unsteady combustion model, yielding more detailed non-normal and non-linear features. In their model, the spatial temperature distribution in the solid propellant is taken into consideration, which leads to the historical effects that are included. Indeed, under the same pressure disturbance, the solid phase will respond differently with different temperature distribution, which is determined by the combustion history of the propellant. The unsteady combustion model used by Mariappan and Sujith [7] is the quasi-steady homogeneous one-dimensional (QSHOD) model equipped with a flame model, named the Krier-T’ien-Sirignano-Summerfield (KTSS) model [8], where the flame is in the gas phase, which is assumed to be a uniform one in a region above the burning surface. The author of the present paper finds that the Zel’dovich-Novozhilov(ZN) model serves better, and investigates the triggering instability and bi-stable region thoroughly in a previous work [9,10].
Unlike the steady combustion models, the unsteady combustion model is designed to be used in the combustion instability study of SRMs. The detailed modeling of the various physical and chemical processes in the steady model is no longer appropriate due to its mathematical complexity. Thus, the unsteady model should not be too simplified so as to include the major dynamics related to the combustion instability, and also not too complicated to be mathematically intractable [1]. The unsteady combustion models of solid propellants were developed separately by the researchers in the USA and the former USSR/Russia. Researchers in the USA adopt the QSHOD approach, where the processes in the condensed phase, the gas phase and the interface between them are detailed, modeled and then combined together [1]. Meanwhile, the Russian researchers follow the phenomenological approach, originally proposed by Zel’dovich, where the burning rates and surface temperature are related to the ambient pressure and temperature directly [11]. The comparison of the two classes of models have been discussed thoroughly since the 1990s, and this is summarized in [12], where the linear and non-linear responses of the two models under pressure oscillations with relative high amplitude are compared. It is well known that the QSHOD model and ZN model are equivalent in the linear region mathematically. In [12], the non-linear responses are compared based on this linear equivalence, which leads to the inevitable drawback that the equivalence is derived under linear assumption while the non-linear responses are investigated.
On account of this, the present paper returns to the original phenomenological construction of the models. The steady combustion features of solid propellants based on the QSHOD model are thoroughly calculated in a wide range of ambient pressure and temperature. The results are then used to determine the phenomenological parameters of the ZN model, yielding a set of ZN parameters for the wide range of ambient conditions. The responses of the ZN model with the parameters calculated from traditionally linear equivalence and the new phenomenological approach are compared, from the linear to non-linear regimes.
After the present introductory part, the mathematical formulations of the QSHOD model and the ZN model, as well as the relations between them, are discussed in Section 2. Based on the mathematical characteristics of the two models, Section 3 discusses the numerical algorithm used to investigate the equations. Using this algorithm, a phenomenological modeling of a QSHOD model is constructed, and the responses of both models are illustrated and discussed in Section 4. Finally, the concluding remarks are given in Section 5.

2. Mathematical Formulations

The mathematical formulations of the QSHOD model and ZN model share most parts in common, which has been discussed thoroughly in the literature [12]. In the present section, the governing equations and the boundary conditions are discussed, as well as their physical significance. However, not all of the derivations and assumptions are included here.

2.1. QSHOD Model with KTSS Flame

To start with, the non-dimensionalized governing equation of the model is [8]
θ τ R θ X 2 θ X 2 = 0 , X [ 0 , ) , τ [ 0 , ) R = θ s m
where τ , X are the non-dimensional temporal and spatial variables, θ , P are the non-dimensional temperature and pressure, and R is the non-dimensional burning rate. It should be noted that the exponential dependence of the burning rate on the surface temperature is an approximation of the exponential Arrhenius’ Law. Theses variables are defined as follows.
τ = t / α r 0 2 ,   X = x / α r 0 , R = r r 0 , P = p p 0 ,   θ = T T a T s , 0 T a ,   H = Q s c ( T s , 0 T a ) ,
where H represents a non-dimensional energy comparing the total heat release of the combustion Q s and the heat capacity of the propellant. The condition under investigation is determined by the ambient pressure p 0 and temperature T a . Under this condition, the steady burning rate and the burning surface temperature are r 0 , T s , 0 , respectively. The spatial and temporal variables are non-dimensionalized by r 0 and the heat diffusivity of the propellant α .
The QSHOD model is in a class of similar models with the same governing equations [12]. The QSHOD-KTSS model is a typical one of the class, where the above governing equation for the condensed phase is equipped with a flame in the gas phase, assumed to be uniform in a region. Under this assumption, the corresponding boundary condition can be derived as [8]
θ s X = P 2 n ( P n / m H ) R H R ,
where the dependence on pressure is introduced. The boundary condition is derived from the mass and energy balance at the phase interface [8]. It should be noticed that the governing Equation (1) is a partial differential equation (PDE) with strong coupling with its boundary condition (2). The equation is derived from a heat conduction problem with the regressing boundary, leading to a term with a coefficient R determined by the boundary value of temperature. Thus, the equation governs the temperature field on a moving domain, and the velocity of the moving boundary, which is also the burning rate R, which is determined by the temperature field. This complicated interaction situation precludes the possibility of an analytical solution. The previous study of the model focus on the analytical solution of its linearization system or on the numerical solution of the full system.
From the equations above, the three major parameters of the QSHOD-KTSS model are n, m, and H. The former two determine the dependence of the burning rate on the ambient pressure and burning surface temperature. H describes the energy released by the combustion of the propellant to the heat capacity.

2.2. ZN Model

The ZN model is a phenomenological approach without detailed modeling of the flame following the notion of Zel’dovich [11,13,14]. The governing equation for the solid phase is the same as the QSHOD model,
θ τ R θ X 2 θ X 2 = 0 , X [ 0 , ) , τ [ 0 , ) ,
where the non-dimensional variables are defined as the same as the QSHOD model.
According to Zel’dovich’s notion, there is no need to model the processes in the gas phase in detail. Alternatively, a phenomenological approach in the macroscopic perspective is adopted. To start with, there are two major features of steady combustion, the steady burning rate r 0 and the steady burning surface temperature T s , 0 , and two major ambient conditions, ambient pressure p 0 and temperature T a . In steady combustion conditions, there exists dependence relations between them.
r 0 = F r ( T a , p 0 ) , T s , 0 = F s ( T a , p 0 )
The specific forms of the functions F r , F s are not known, and they are not necessary to be known. Usually, a derived form of the above relations is used, based on the steady distribution of temperature fields in the solid propellant [14]. It should be noticed that the major feature of the above relation can be described by the following partial derivatives.
k = ( T s , 0 T a ) ln r T a p 0 ,   r Z N = T s T a p 0 , n = ln r ln p 0 T a ,   μ = 1 T s , 0 T a T s ln p 0 T a .
The central assumption of the ZN model is Zel’dovich’s notion that the chemical and gas dynamical processes in the gas phase are much faster than the heat conduction processes inside the condensed phase. Under this condition, the relaxation times of the gas phase are much smaller than those in the condensed phase, and the gas phase reaches equilibrium immediately, as the temperature in the condensed phase varies. Mathematically, the Equation (4) can be extended to the unsteady case [11]. For example, the burning rate r under a pressure disturbance p can be expressed as r = F r ( T a , p 0 + p ) . In this manner, the burning rate and the burning surface temperature of the unsteady combustion can be expressed by the parameters k , r Z N , n , μ , defined in Equation (5). The burning rate and boundary condition of the ZN model can be derived as
R = P n exp k θ s φ R , R = exp k ( θ s 1 ) r Z N θ ( τ , 0 ) X = R k ( n ln P + k θ s ln R ) = φ ,
where the heat flux at the burning surface is φ . The readers can refer to the literature [11,14] for detailed derivation processes.
Thus, the ZN model, consisting of Equations (3) and (6) and the QSHOD model, consisting of Equations (1) and (2) have the same governing equation as the two models share the same solid phase modeling. The burning rate expressions and boundary conditions in the two models are completely different, as the modeling methods for the gas phase and phase interface are different. Despite the obvious difference, it is well known that the linearization of the two models yields the equivalent results [1,14].
In the previous work [12], the two models are also compared based on a linear equivalence. The model variables of the QSHOD model are n , m , H , while the model variables of the ZN model are n , k , r Z N . Except for the pressure exponent n, the equivalent relation of the other two variables is
H = 2 k 1 + r Z N 2 k , m = k r Z N ,
which can be derived in the sense that the derivative of the burning rate on the surface temperature is the same [12]. Obviously, this linear equivalence holds near the reference state, where the derivatives are taken, which may not valid for a wide parameter range. Thus, a new approach is proposed in the present paper.

2.3. The Phenomenological Reconstruction of the QSHOD Model

From the discussion in the previous subsection, the linear equivalence holds accurately only on one state point. In the present paper, the calculations are taken in a wide state range, and the ZN parameters n , k , r Z N , μ are evaluated based on the information throughout the range comprehensively.
The method adopted in the present paper is based on the phenomenological notion of the ZN model. The steady state combustion states are evaluated by the QSHOD model with varying ambient pressures and temperatures, and the ZN parameters are then calculated according to the definition (5). In the QSHOD-KTSS model, the two major coefficients determining the main feature of the model are H and m. Thus, with the fixed H , m , the steady combustion can be solved by evolving the following equation to the steady state,
T τ R T X 2 T X 2 = 0 , x [ 0 , ) , t [ 0 , ) R = T s T a T s , 0 T a m T s X = ( T s , 0 T a ) P 2 n ( P n / m H ) R + H R ,
where only the temperature is reverted to the dimensional form to take the ambient temperature variation into consideration.
For each ambient state, ( P i , T a , j ) , one can solve a steady state from the above equations, characterized by the burning rate and surface temperature R ( i , j ) and T s , ( i , j ) . Based on the two matrices of { R } and { T s } , the parameters can be evaluated by the rows and columns. For example, k can be evaluated by each row of the matrix { R } , where the information in row { R ( i , · ) } will yield a k at fixed ambient pressure P i with varying ambient temperature T a , j , , T a , N p by a regression algorithm. In the present paper, an N p × N p mesh of ambient states are adopted, which can be adjusted for specific ranges as needed. Solving the above system for N p different P and N p different T a yields an N p × N p matrix of { R } N p × N p and an N p × N p matrix of { T s } N p × N p , the four coefficients in Equation (5) can then be calculated by regression in each row and column (the total number for each parameter is N p ). In this manner, a ZN model based on the results of the QSHOD-KTSS model is derived. The detailed modeling of the flame model is replaced by the phenomenological relation Equation (4) or Equation (5).

3. Numerical Methods

As mentioned earlier in Section 2, the governing PDE of the two models is difficult to solve due to the existence of the coupling term. The numerical algorithm to deal with the PDE is the spectral method proposed by Novozhilov in [13,14], which will be introduced briefly in this section.

3.1. The Expansion by the Lagurre Polynomials

The temperature in Equation (1) or (3) is expanded by the Laguerre polynomials as
θ ( τ , X ) e X i = 0 N u i ( τ ) L i ( X ) , L i ( X ) = e X i ! d i d X i e X X i ,
where L i ( X ) is the i-th order Laguerre polynomial, and the total number of Laguerre polynomials is set to be N. The derivatives of the polynomial has the relation
d L i ( X ) d X = s = 0 N 1 L s ( X ) , d 2 L i ( X ) d X 2 = s = 0 N 2 ( N s 1 ) L s ( X ) .
By substituting the expansion Equation (9) into the governing Equation (1) or (3), and using the above derivatives, the PDE can be reduced to a system of ordinary differential equations
u ˙ 0 = φ R θ s , u ˙ 1 = φ θ s + R ( u 0 θ s ) , u ˙ i = φ i θ s + j = 0 i 2 ( i j 1 ) u j + R j = 0 i 1 u j θ s , i = 2 , 3 , , N ,
which is an ( N + 1 ) -th order system of ODEs of u i . Additionally, the surface temperature θ s and the surface heat flux φ have a simple expression following the definition of the Laguerre polynomial,
θ s = i = 0 N u i , φ = i = 0 N ( i + 1 ) u i .
With the above equations, the system (1) or (3) is spatially discretized, leading to a explicit form, which can be easily integrated by a standard time-forwarding algorithm. The validation of the algorithm is included in [12], where an appropriate order of the polynomials used in the algorithm is shown to be N = 15 .

3.2. The Phenomenological Parameters

As discussed in Section 2, the phenomenological reconstruction of the QSHOD-KTSS model is based on the steady combustion results under different ambient pressures and temperatures. This can be complete by a temporal forwarding of the unsteady Equation (8) by a similar expansion,
T ( τ , X ) e X i = 0 N u i ( τ ) L i ( X ) .
A convergence criterion is used to check if the steady combustion is reached.
By conducting the steady combustion calculation in a range of ambient pressure P 1 , P 2 , , P N p and temperature T a , 1 , T a , 2 , , T a , N p , the steady combustion in N p × N p conditions are evaluated.
To calculate the phenomenological parameters in Equation (5), a straightforward method is to approximate the partial derivatives by finite differences. However, to collect the information in each state points together, the linear regression is used. Take the pressure exponent as example. The pressure exponent is the derivative of ln R to the ln P at fixed ambient temperature. Thus, for a selected temperature T a , j , there is a vector of pressure P 1 , P 2 , , P N p and a row of corresponding burning rates R 1 , j , R 2 , j , , R N p , j . After conducting the linear regression for the logarithm of both series, the slope is n j , which is the pressure exponent in this ambient temperature T a , j . In this manner, the N p pressure exponent will be calculated. Repeating the processes for the other parameters, N p values for each parameter can be calculated.

4. Results and Discussion

In this section, the unsteady combustion under varying pressure oscillations of the models will be discussed. The models in the comparison are the QSHOD model and the two equivalent models, in the sense of linear equivalence and the phenomenological approach, respectively.

4.1. Model Parameters

The basic parameters for the solid propellant under investigation are collected in Table 1. The calculations are conducted in the non-dimensional form, according to (8), and the results are reverted back and summarized in the dimensional form.

4.2. The Equivalent ZN Models

Traditionally, the equivalent counterpart of a QSHOD model in the ZN framework should be based on the linear equivalence [12,15]. In the present model, the corresponding ZN parameters can be calculated by Equation (7)
k L = 1.8293 ,   r Z N L = 0.3049 ,
where the superscript denotes the linear equivalence.
In the present paper, a new phenomenological approach is shown in Figure 1, where the logarithmic coordinates are used when the logarithm of the quantity is used in Equation (5). In each subfigure, five typical conditions are selected to draw the variation. It can be seen from the figure that the lines are of good linearity, and the slope of each line is almost the same. This indicates that, in the selected parameter range, the ZN parameters are close to constant. However, this does not imply that the above linear equivalence method is valid.
In the present paper, the ZN parameters are calculated according to the algorithm in Section 3.2. Under the N p × N p mesh, there are N p results for each parameter, and their mean value is used for the final results,
k = 2.034 , r Z N = 0.6873 , n = 0.3775 , μ = 0.0577 .
It is noticeable that the parameters generated from the new approach are quite different from the traditional linear equivalence results k L , r Z N L , and the pressure exponent is a bit larger.
Comparing the two methods, the linear equivalence is derived from the derivatives calculated at the reference state, while in the phenomenological method, the parameters are calculated from the comprehensive effects throughout the parameter range. The difference originates from two aspects. The first one is the effects of different parameters, which are coupled together, and not separated. The pressure exponent n is an example. Although the pressure exponent is introduced into the QSHOD-KTSS model explicitly in the boundary condition (8), its influence is not straightforward, and the pressure exponent of the full system differs from the original set value. Another important aspect is the wide state range. The four parameters are not constant throughout the state range and the dependences of r , T s on T a , p are not strictly linear. The results illustrate the discrepancies between the linear equivalence at a single state point and the overall variations in the state range.

4.3. The Unsteady Combustion of the QSHOD-KTSS Model

The comparison of the unsteady combustion of the QSHOD model and the ZN model has been discussed in [12] under various conditions. In the present paper, the harmonic excitation with increasing amplitude is discussed.
In the harmonic excitation problem, the unsteady combustion is driven by a periodically varying ambient pressure, which is the case in the chamber of an SRM under combustion instability. The driving signal is
P = 1 + Δ P sin ( 2 π τ ) ,
where the non-dimensional form is used. According to Equation (11), the pressure is oscillating around the reference pressure at an amplitude Δ P .
It has been shown in the literature that the non-linear dynamics of the unsteady combustion under the increasing driving amplitude is a series of period-doubling bifurcations to chaos [11,14]. In the present paper, the comparison is also made in this problem. In each calculation, time is evolved to τ = 500 and the range τ [ 400 , 500 ] is plotted to cut off the transient processes. The reason for cutting off the transient processes is the non-linear dynamics of the system, which varies mainly on the final states, which is the attractor of the system. Indeed, it will be shown that the combustion system will be attracted to a series of geometrically different attractors with increasing amplitude.
To start with, at Δ P = 0 , which corresponds to a steady constant pressure condition, the burning rate response is also steady at unity. With an increasing Δ P , there will be a burning rate oscillation, correspondingly. When the amplitude is small, the linear simplification applies, and a harmonic pressure perturbation will lead to a harmonic burning rate response with the same frequency. The burning rate responds to the pressure oscillation in a “linear” manner when Δ P < 0.28 , as shown in Figure 2a. In the three subfigures, the top one is the time series of the burning rate, the middle one is the phase portrait by phase reconstruction and the bottom one is the FFT of the time series. Actually, the response in this region has begun to deviate from the linear response, but the period of the burning rate is the same as the driving pressure. Mathematically, the attractor of the system is qualitatively the same for 0 < Δ P < 0.28 .
At the point Δ P = 0.28 , a period-doubling bifurcation occurs, where the limit cycle becomes self-intersected and the period of the cycle doubles, as shown in Figure 2b. As the attractor becomes geometrically different, the final state of the system is qualitatively different from the Δ P < 0.28 case, and a bifurcation occurs. With the further increasing of the Δ P , there will be a series of period-doubling to period-4 ( Δ P = 0.72 ), period-8, period- 2 n and finally to chaos ( Δ P = 0.86 ), as shown in Figure 2c,d.
It is noticeable that the burning rate response becomes extremely high as Δ P > 0.2 , which is unphysical. This is due to the over-simplification of the flame structure, which may lead to wrong results when applied in the combustion instability study of the SRMs.

4.4. The Unsteady Combustion of the Two ZN Models

Based on the linear equivalency, the ZN parameter should be k L = 1.8293 , r Z N L = 0.3049 , n = 0.3 . The corresponding unsteady combustion responses are collected in Figure 3. The overall bifurcation scenario is the same as the QSHOD model, qualitatively. Nevertheless, from the quantitative perspective, there are considerable differences. The bifurcation points are smaller, or each period-doubling becomes earlier. Additionally, the burning rate responses are confined in a realistic range. Thus, the ZN model should be more appropriate for the combustion instability study than the QSHOD model, as discussed in [12]. However, although the ZN model by linear equivalence is equivalent to the QSHOD model according to non-linear dynamics, in the sense of the bifurcation scenario, the quantitative features are quite different; especially the exact position of the bifurcation points.
In the new phenomenological approach proposed in the present paper, the ZN parameter should be k = 2.034 ,   r Z N = 0.6873 ,   n = 0.3775 . The corresponding unsteady combustion responses are shown in Figure 4. Again, the overall bifurcation scenario is the same as the previous two cases. But, the bifurcation points of the new model become more close to the original QSHOD model than the ZN model based on the linear equivalency.
Indeed, the typical bifurcation points of the three systems are collected in Table 2, and the ZN model from linear equivalence has much smaller bifurcation points. It should be noted that the unphysical large values of both bifurcation points and response amplitudes in the QSHOD-KTSS model are due to the simplifications in the modeling processes. For example, the burning rate R = θ s m is a simplification of the exponential function in Arrhenius’ Law, which is valid near the reference state or R = 1 . Admittedly, the accuracy of the QSHOD model is not the focus of the present paper. However, the results show that the unphysical late bifurcation effects are better captured by the phenomenological approach than the linear equivalence approach.
By the comparison of the two ZN models, there also exists the non-linear equivalency in the sense of the same bifurcation scenario. The ZN model from the phenomenological approach reconstructs the system better, as each bifurcation point is closer to the original QSHOD model than the linear equivalency. Indeed, the phenomenological approach makes full use of the information in a wide range of parameters, while the linear equivalency only use the information near the reference state.

5. Conclusions

A new unsteady combustion modeling framework is investigated based on the Zel’dovich phenomenological approach. The original unsteady combustion model based on flame modeling is reconstructed in the phenomenological form, with both good qualitative and quantitative agreements.
The spectral method is modified to investigate the combustion features of the QSHOD-KTSS model in a wide state range, yielding combustion characteristics under different ambient conditions. The results are then analyzed by linear regressions to calculate the four ZN parameters, in a phenomenological manner.
The unsteady burning rate responses of the the new constructed model are then compared with the traditional ZN model based on linear equivalence. The results show that the new phenomenological approach is able to reconstruct the original model better than the traditional linear equivalence model, in the sense that the bifurcation scenario is closer. This indicated that the phenomenological approach, in which information in a widely varying state range is taken into consideration, is able to capture dynamical properties of the original model better than the linear equivalence approach, in which information in only a single state is considered.
The proposed methodology can be used for more complicated combustion modeling for solid propellants, or even liquid propellants. Indeed, the modeling of unsteady combustion responses of different solid or liquid propellants focuses on the burning rate responses of various perturbations of ambient conditions. The phenomenological approach can be applied in these modeling problems.

Author Contributions

Conceptualization, Z.W.; methodology, Z.W.; software, Z.W.; validation, W.Z. and Y.L.; formal analysis, Z.W.; investigation, Z.W.; resources, Z.W.; writing—original draft preparation, Z.W.; writing—review and editing, Z.W., W.Z.; visualization, Z.W.; project administration, Z.W.; funding acquisition, Z.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Fundamental Research Funds for Central Universities (Grant No. G2021KY05108).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
SRMSolid rocket motor
QSHODquasi-steady homogeneous one-dimensional
KTSSKrier-T’ien-Sirignano-Summerfield
ZNZel’dovich-Novozhilov
PDEPartial differential equation
ODEOrdinary differential equation

References

  1. Culick, F. Combustion Instabilities in Solid Propellant Rocket Motors; Technical Report; California Institute of Technology Pasadena: Pasadena, CA, USA, 2004. [Google Scholar]
  2. Juniper, M.P.; Sujith, R. Sensitivity and nonlinearity of thermoacoustic oscillations. Annu. Rev. Fluid Mech. 2018, 50, 661–689. [Google Scholar] [CrossRef]
  3. Gallier, S.; Briquet, B.; Yiao, M. Aluminum Combustion Can Drive Instabilities in Solid Rocket Motors: T-Burner Study. J. Propuls. Power 2019, 35, 159–172. [Google Scholar] [CrossRef]
  4. Culick, F.E.C.; Burnley, V.; Swenson, G. Pulsed instabilities in solid-propellant rockets. J. Propuls. Power 1995, 11, 657–665. [Google Scholar] [CrossRef]
  5. Culick, F.E.C. Nonlinear behavior of acoustic waves in combustion chambers—I. Acta Astronaut. 1976, 3, 715–734. [Google Scholar] [CrossRef]
  6. Culick, F. Nonlinear behavior of acoustic waves in combustion chambers—II. Acta Astronaut. 1976, 3, 735–757. [Google Scholar] [CrossRef]
  7. Mariappan, S.; Sujith, R. Thermoacoustic instability in a solid rocket motor: Non-normality and nonlinear instabilities. J. Fluid Mech. 2010, 653, 1–33. [Google Scholar] [CrossRef]
  8. Krier, H.; Sirignano, W.; Summerfield, M.; T’ien, J. Nonsteady burning phenomena of solid propellants-Theory and experiments. AIAA J. 1968, 6, 278–285. [Google Scholar] [CrossRef]
  9. Wang, Z.; Liu, P.; Ao, W. A reduced-order model of thermoacoustic instability in solid rocket motors. Aerosp. Sci. Technol. 2020, 97, 105615. [Google Scholar] [CrossRef]
  10. Wang, Z.; Liu, P.; Jin, B.; Ao, W. Nonlinear characteristics of the triggering combustion instabilities in solid rocket motors. Acta Astronaut. 2020, 176, 371–382. [Google Scholar] [CrossRef]
  11. Novozhilov, V.B.; Novozhilov, B.V. Theory of Solid-Propellant Nonsteady Combustion; John Wiley & Sons: Hoboken, NJ, USA, 2020. [Google Scholar]
  12. Wang, Z.; Zhang, X.; Bai, Y.; Gan, X.; Wang, J.; Jin, B.; Liu, P. A Comparison of the Nonlinear Unsteady Responses of Two Classical Models for Solid Propellant. In Proceedings of the 2022 13th International Conference on Mechanical and Aerospace Engineering (ICMAE), Bratislava, Slovakia, 20–22 July 2022; pp. 258–263. [Google Scholar]
  13. Novozhilov, B.V. Randomization of the unsteady burning rate of gunpowder. Russ. J. Phys. Chem. 2004, 23, 68–74. [Google Scholar]
  14. Novozhilov, B. Combustion of energetic materials in an acoustic field. Combust. Explos. Shock Waves 2005, 41, 709–726. [Google Scholar] [CrossRef]
  15. Zarko, V.; Gusachenko, L. Critical review of phenomenological models for studying transient combustion of solid propellants. Int. J. Spray Combust. Dyn. 2010, 2, 151–167. [Google Scholar] [CrossRef]
Figure 1. The phenomenological calculation of ZN parameters: (a) k, (b) n, (c) r Z N , (d) μ .
Figure 1. The phenomenological calculation of ZN parameters: (a) k, (b) n, (c) r Z N , (d) μ .
Aerospace 10 00767 g001
Figure 2. The burning rate responses of the QSHOD model under harmonic excitation pressure with different amplitudes: (a) Δ P = 0.1 , (b) Δ P = 0.281 , (c) Δ P = 0.72 , (d) Δ P = 0.86 .
Figure 2. The burning rate responses of the QSHOD model under harmonic excitation pressure with different amplitudes: (a) Δ P = 0.1 , (b) Δ P = 0.281 , (c) Δ P = 0.72 , (d) Δ P = 0.86 .
Aerospace 10 00767 g002
Figure 3. The burning rate responses of the ZN model (linear equivalence) under harmonic excitation pressure with different amplitudes: (a) Δ P = 0.1 , (b) Δ P = 0.28 , (c) Δ P = 0.41 , (d) Δ P = 0.6 .
Figure 3. The burning rate responses of the ZN model (linear equivalence) under harmonic excitation pressure with different amplitudes: (a) Δ P = 0.1 , (b) Δ P = 0.28 , (c) Δ P = 0.41 , (d) Δ P = 0.6 .
Aerospace 10 00767 g003
Figure 4. The burning rate responses of the ZN model (phenomenological approach) under harmonic excitation pressure with different amplitudes: (a) Δ P = 0.1 , (b) Δ P = 0.28 , (c) Δ P = 0.62 , (d) Δ P = 0.86 .
Figure 4. The burning rate responses of the ZN model (phenomenological approach) under harmonic excitation pressure with different amplitudes: (a) Δ P = 0.1 , (b) Δ P = 0.28 , (c) Δ P = 0.62 , (d) Δ P = 0.86 .
Aerospace 10 00767 g004
Table 1. The physical and chemical parameters of the solid propellant.
Table 1. The physical and chemical parameters of the solid propellant.
ParameterValueUnit
Reference pressure, p 0 6MPa
Reference ambient temperature, T a , 0 300K
Reference burning surface temperature, T s , 0  11100K
Reference burning rate, r 0 6.871mm/s
Pressure range[3, 9]MPa
Temperature range[250, 350]K
Pressure exponent, n0.3
Energy constant, H0.81
Temperature exponent, m6
1 The burning surface temperature of the propellant at the reference state.
Table 2. The bifurcation points of different models.
Table 2. The bifurcation points of different models.
ModelPeriod-2 DoublingPeriod-4 DoublingChaos
QSHOD-KTSS0.2810.720.86
ZN (linear equivalence)0.280.410.6
ZN (phenomenological approach)0.280.620.86
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, Z.; Zhang, W.; Liu, Y. A Phenomenological Model for the Unsteady Combustion of Solid Propellants from a Zel’dovich-Novzhilov Approach. Aerospace 2023, 10, 767. https://doi.org/10.3390/aerospace10090767

AMA Style

Wang Z, Zhang W, Liu Y. A Phenomenological Model for the Unsteady Combustion of Solid Propellants from a Zel’dovich-Novzhilov Approach. Aerospace. 2023; 10(9):767. https://doi.org/10.3390/aerospace10090767

Chicago/Turabian Style

Wang, Zhuopu, Wenchao Zhang, and Yuanzhe Liu. 2023. "A Phenomenological Model for the Unsteady Combustion of Solid Propellants from a Zel’dovich-Novzhilov Approach" Aerospace 10, no. 9: 767. https://doi.org/10.3390/aerospace10090767

APA Style

Wang, Z., Zhang, W., & Liu, Y. (2023). A Phenomenological Model for the Unsteady Combustion of Solid Propellants from a Zel’dovich-Novzhilov Approach. Aerospace, 10(9), 767. https://doi.org/10.3390/aerospace10090767

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