UNDERSTANDING of unsteady combustion of solid
|
|
- Wendy Cook
- 6 years ago
- Views:
Transcription
1 JOURNAL OF PROPULSION AND POWER Vol. 4, No. 5, September October 8 Analytical Solution or Pressure-Coupled Combustion Response Functions o Composite Solid Propellants Michael Shusser and Fred E. C. Culick Caliornia Institute o Technology, Pasadena, Caliornia 9115 and Norman S. Cohen Cohen Proessional Services, Los Angeles, Caliornia 935 DOI: 1.514/1.15 This paper etends the classical analytical solution or small perturbation analysis o the pressure-coupled response o a homogeneous propellant to any two-component composite propellant. The solution obtained is general and can be used with any particular model or propellant combustion. As an eample, the Cohen and Strand ammonium perchlorate propellant model or a single ammonium perchlorate particle size was used in this work. The results and their mechanistic signiicance are presented and discussed. It is shown that, or a two-component composite propellant, two orms o pressure eponents arise rom the analysis. The signiicance o the second eponent is that it enables the composite propellant to be viewed as a homogeneous propellant with a requencydependent eponent via the coupling coeicients. It is ound that the ammonium perchlorate is the main source o instability because o its condensed phase eothermicity and monopropellant lame kinetics. This will be a problem with energetic materials in general. The inert binder provides a stabilizing inluence because o its endothermicity and the diusion lame ormed with the ammonium perchlorate. Eects o ammonium perchlorate particle size and pressure stem rom the changing lame structure and its eect on burning rate. Nomenclature A = dimensionless parameter characteristic o surace decomposition B = dimensionless parameter characteristic o heat eedback law C = coupling parameter D = deined by Eqs. (3) and (4) F = classical response unction, Eq. (17) = requency o oscillations n = burn rate pressure eponent n = deined by Eq. (49) p = pressure R = response unction r = burning rate S = raction o propellant surace area occupied by a component T = temperature T = bulk temperature t = time = distance (negative into the solid) z = dimensionless distance = raction o propellant mass low rate rom ammonium perchlorate = = 1 = dimensionless temperature = thermal diusivity = ammonium perchlorate comple requency Presented as Paper 346 at the 37th AIAA/ASME/SAE/ASEE Joint Propulsion Conerence and Ehibit, Salt Lake City, UT, 8 11 July 1; received 1 December 5; revision received 18 May 8; accepted or publication 6 May 8. Copyright 8 by the American Institute o Aeronautics and Astronautics, Inc. All rights reserved. Copies o this paper may be made or personal or internal use, on condition that the copier pay the $1. per-copy ee to the Copyright Clearance Center, Inc., Rosewood Drive, Danvers, MA 193; include the code /8 $1. in correspondence with the CCC. Postdoctoral Scholar; currently Faculty o Mechanical Engineering, Technion, 3 Haia, Israel. Member AIAA. Proessor o Mechanical Engineering and Jet Propulsion. Fellow AIAA. Subcontractor. Associate Fellow AIAA. 158 = binder comple requency = density = dimensionless time = dimensionless requency Subscripts = binder Im = imaginary part prop = propellant Re = real part s = surace value = oidizer 1 = oidizer = binder Superscripts (Im) = imaginary part (Re) = real part = perturbations I. Introduction UNDERSTANDING o unsteady combustion o solid propellants is a basis or analyzing, predicting, and suppressing o combustion instabilities in solid rocket motors. The main cause o these instabilities is the response o the combustion process to pressure oscillations in the combustion chamber. Hence, the most important physical quantity characterizing unsteady combustion is the response o mass burning rate to pressure oscillations. The simplest models o the unsteady combustion o solid propellants [quasi-steady homogeneous one-dimensional (QSHOD) models [1]] assume that the propellant is homogeneous, the gas phase is quasi steady, and the solid phase is one-dimensional. This approimation results in the well-known classical analytical solution or small perturbation analysis o a homogeneous propellant []. I one deines the pressure-coupled response unction R p R p r r p (1) p
2 SHUSSER, CULICK, AND COHEN 159 where r, p are steady-state values (overbars denote steady-state values throughout the paper) o burning rate and pressure, and r, p are their perturbations, respectively, then nab R p A= 1 A1 B Here, is a nondimensional comple requency o oscillations, whereas each o the basic parameters n, A, and B characterize the dependence o the steady-state burning rate r on a certain physical quantity. Thus, the pressure eponent n characterizes the dependence o r on pressure, A (characteristic o the surace decomposition) is related to the dependence on the surace temperature, and B (characteristic o the heat eedback law) shows the dependence on the surace temperature gradient, that is, heat lu. The classical solution has been used etensively in combustion instability studies. It is, however, limited to homogeneous propellants, whereas most modern solid propellants are composite propellants. To model composite propellant combustion, one must account or the heterogeneity and two-dimensionality o the problem. Following the classical work o Nachbar [3], a widely used approach has been to describe these eects approimately by using a sandwich model. This model represents a composite structure as side-by-side slabs o uel and oidizer. Steady burning o such sandwiches has been the object o numerous studies [4]. Cohen and Strand perturbed their steady-state model [5] to describe the nonsteady combustion o a composite propellant [6,7]. They were able to calculate the linearized combustion response and to eplain eects o ammonium perchlorate (AP) particle size, pressure, and crosslow on combustion response properties. Shusser et al. [8] used a sandwich model to numerically compute nonlinear combustion response properties o a composite propellant by solving the heat conduction equation in each component, with the Cohen and Strand model [5] as a boundary condition. A similar approach was used by Rasmussen and Frederick [9] and Surzhikov et al. [1]. Recently, a combustion model based on calculating the burning o a three-dimensional array o randomly packed spheres was proposed [11,1]. This model is irst to enable detailed consideration o the microstructure o the propellant. However, as with simpler sandwich models, the predictive value o the model is limited by the approimate constitutive relations it uses. Moreover, the need or computationally ast models remains or uture coupling o the propellant combustion with the gas dynamics o the combustion chamber in comprehensive motor analyses. It is well recognized that closed-orm solutions are very useul, as they provide insights not otherwise available. A general analytic solution or the response unction o composite propellants has not yet been obtained. Previous attempts to use the classical solution () or AP composite propellants were unsuccessul. For eample, Brown and Muzzy [13], who tried to use solution () to match their eperiments, obtained physically unrealistic values o the parameters. The reason or this ailure is obvious. Although a total propellant burning rate eists, surace temperature and heat lu are dierent or the AP and the binder components o the propellant. There are no average or total surace temperatures or surace temperature gradients and, thereore, no A and B parameters or the propellant as a whole can be deined. In the present work, we obtain the irst closed-orm analytical solution or a heterogeneous unsteady model o composite propellants. Such solution can serve several purposes. First, it can give more insight into the physics o the process. It also provides a better way to see trends with changes in physical conditions, such as pressure or temperature. It can be very helpul or adjusting parameters to achieve better modeling. Finally, it serves as a benchmark or numerical computations. Though only pressure-coupled response is considered in this analytical paper, the model can be etended to include velocity coupling by including an appropriate dynamic erosive burning model, assuming that dynamic erosive burning is the essential mechanism. This was done in the numerical work [8]. () In the net section, we will demonstrate how to obtain such an analytical solution or a composite propellant or a single (monomodal) particle size. For such a coupled two-component system, each component has its own set o basic parameters analogous to A, B, n, whereas two additional coeicients are due to the coupling between the oidizer and the binder. This solution is general in the sense that it is not limited to any particular ingredient or model or propellant combustion. A choice o a speciic model inluences only the values o the parameters, as they are obtained by linearizing the model. In Sec. III, we calculate the coeicients or the Cohen and Strand AP propellant model [5 7], show how they change with pressure and particle size, and compare the analytical solution with numerical computations (the numerical scheme has been described elsewhere [8]). We then proceed with a parametric study o the general solution in Sec. IV. Intrinsic instability is considered in Sec. V. The pressure eponent o a composite propellant is studied in Sec. VI. II. Analytical Solution The physical model used or the analysis is the classical AP/binder sideways sandwich. The approach yields a system o two unsteady one-dimensional heat conduction equations coupled through their boundary conditions at the solid surace. Their linearized orm, which is appropriate or small perturbation analysis, is as @z r 1 r 1 e z 1 (3) r r e z (4) Here, the subscripts 1, relate to the AP and the binder, respectively, the overbar denotes steady-state values, and the primed quantities are their perturbations. The dimensionless temporal and spatial coordinates 1,, z 1, z, and temperature perturbations 1, are deined by z 1 r 1 1 ; 1 r 1 t 1 ; z r (5) r t (6) 1 T 1 T 1 T 1s T ; T T T s T (7) Here, is the distance rom the solid surace (negative into the solid), t is the time, 1 and are AP and binder thermal diusivities, T 1s and T s are AP and binder surace temperatures, and T is the temperature in the bulk o the solid. The classical theory [] uses the ollowing boundary condition at the solid surace o a homogeneous p z p 1 B 1 r (8) A r Following the classical ormulation [], we write the boundary conditions at the solid surace and the relationships between burning rate and temperature perturbations as z 1 n B p p 1 B 1 A r 1 r 1 C r r (9) We use AP to represent the oidizer, but the generality o the approach is applicable to any two-component system.
3 16 SHUSSER, CULICK, AND COHEN p B p C r1 1 B r 1 r (1) 1 A r r 1 r 1 A 1 j z 1 ; r r A j z (11) One sees that now both AP and binder have their own set o basic parameters n, A, B, describing the relations between corresponding surace temperatures, heat lues, and burning rates. The n and n arise rom the burning o the composite propellant: n is not the value or AP by itsel as a monopropellant, and or the inert binder by itsel, n would be meaningless. Their meaning will be discussed in Sec. VI. The coupling coeicients C, C characterize how each heat lu is inluenced by the other component s burning rate. Physically, they are caused mainly by the coupling through the multiple lame structure. It should be noted that, due to the coupling, each component (AP or binder) behaves dierently when in a propellant and when by itsel. Thus, the values o the AP parameters A, B will also be dierent rom those or pure AP [14], as we shall see in the net section. We emphasize that the values o the parameters cannot be obtained by considering each component separately. They are determined by the physics o unsteady combustion o a composite propellant and can be obtained rom its suitable model. To complete the statement o the problem, we add the boundary conditions in the bulk o the solid z 1! ; z!1 (1) We now impose on the system an oscillatory pressure perturbation o requency. Deining two dimensionless oscillatory requencies 1 1 r ; 1 r (13) and AP and binder response unctions R r 1 p r 1 p ; R r p r p (14) one obtains (the details can be ound elsewhere [15]) n B C F R (15) =A 1= 1B 1=A C C =n B F n B C F R (16) =A 1= 1B 1=A C C =n B F Here, F is shorthand or the classical response unction () n F A B A = 1 A 1 B (17) with a corresponding ormula or F. The comple requencies, are obtained rom i 1 ; i (18) One sees rom Eqs. (15) and (16) that, when one o the coupling coeicients vanishes, the corresponding response unction reduces to its homogeneous orm [Eq. ()]: R F or C ; R F or C (19) The real and imaginary parts o R and R are written out in the Appendi. III. Cohen and Strand Model or Constitutive Combustion Relations The solution obtained in the previous section depends on several parameters. The values o the parameters can be determined by a particular model or propellant combustion. We have used the Cohen and Strand model [5] or this purpose, as was veriied by incorporating it as the constituent combustion model in our response unction work and comparing numerical solutions o the present model with a large amily o response unction data in [8]. Thus, to estimate the physically reasonable values o the combustion response parameters, we are justiied in using the Cohen and Strand model [5] to irst calculate the coeicients. The way to derive the eight coeicients A, A, B, B, n, n, C, C is to linearize the nonlinear surace boundary condition as obtained rom the original model (see [5 8]) and to compare its linearized orm with Eqs. (9) and (1). The resulting detailed epressions are given elsewhere [15]. Figure 1 compares the analytical solution with numerical computations or the individual component AP and binder response unctions as they behave in the composite propellant or a particular case. Both real and imaginary parts are shown in the igure. One sees that the agreement is very good. We studied the inluence o absolute pressure and particle size on the coeicients. AP concentration was kept at 87% throughout the calculations. Propellant bulk temperature was 5 C. We show here only the main results. The details can be ound in [15]. The B parameter is shown in Fig.. We see immediately that B in a propellant is much higher than in the classical ormulation. Although the B parameter usually lies between.6 and 1, B can be almost 3. On the other hand, or high pressures and large particle sizes, B becomes less than one; the reason or this is that, under these conditions, C, as we shall see later, and hence the less stable AP is practically decoupled rom the binder. B decreases with increasing pressure mainly because o decreasing lame heights. The jump at about 8 atm or the m case is caused by a change in the competing lame structure, as relected in the F model parameter (see [5,16] or more details). The binder parameters A, B, n may not be compared with the classical homogeneous A, B, n because an inert binder does not burn by itsel. Here, we show only B in Fig. 3. One sees that B is even larger than B and always stays above.4. Thus, an inert binder is highly stabilizing to the AP. The calculations also demonstrate [15] that n decreases much aster than n with increasing pressure because the binder is heated only by the diusion lame. AP is also heated by its monopropellant lame, which is the destabilizing inluence, as can be inerred rom our results in [14]. Figures 4 and 5 are plots o the coupling coeicients C and C. C is always negative and increases with pressure and particle size. This relects the stabilizing inluence o the binder on the AP by the coupling. The binder decomposition is endothermic and urnishes a diusion lame. C vanishes or high pressures and large sizes, where the AP becomes decoupled rom the binder and behaves more like pure AP. The dependence o C on pressure and particle size is more complicated, as seen rom Fig. 5. It is positive ecept or m AP at low pressures. This means that the AP generally has a destabilizing inluence on the binder when the binder is in act burning. The AP monopropellant lame and eothermic condensed phase tend to be destabilizing. The analytical solution has served to clariy this or a composite propellant. IV. Parametric Study A parametric study o the analytical solution [Eqs. (15) and (16)] is not easy because o the large number o parameters. The solution depends on eight coeicients A, A, B, B, n, n, C, C, and also It is possible that B will decrease with decreasing pressure at very low pressures, atm, due to the importance o condensed phase eothermicity o the AP at low pressures [8,14]. However, the presence o the binder has a stabilizing inluence on the AP.
4 SHUSSER, CULICK, AND COHEN 161 Fig. B parameter or the Cohen and Strand model [5]. Fig. 3 B parameter or the Cohen and Strand model [5]. Fig. 1 Comparison o analytical and numerical solution or composite propellant response unction; 87% AP ( m)/hydroyl-terminated polybutadiene, 68 atm, 98 K: a) AP response unction, b) binder response unction. on the ratio o dimensionless requencies that we will denote as : = 1 () We thereore adopt the ollowing approach. We will choose reasonable basic values o each parameter rom Sec. III results and vary only one o them at a time. Only the real part o the AP response unction as it behaves in the propellant will be considered due to its being the most important driver or practical applications. Results o this eercise can provide guidance or more study in uture work. The basic parameters chosen are given as ollows: A 1, B 1:5, n 1, A 8, B 3, n :8, C :3, C :5, :. Fig. 4 Coupling coeicient C or the Cohen and Strand model [5].
5 16 SHUSSER, CULICK, AND COHEN Fig. 5 Coupling coeicient C or the Cohen and Strand model [5]. All o the parameters are naturally divided into our groups: 1) AP coeicients A, B, n ; ) binder coeicients A, B, n ; 3) coupling coeicients C, C ; and 4) requency ratio. The inluence o the AP coeicients turned out to be analogous to A, B, n in the classical solution. The increase in A increases the response and moves its maimum to higher requencies. Decreasing B makes the response unction curve more peaked, whereas larger n increases the value o the response but not its requency dependence. The inluence o the binder coeicients is shown in Figs We see that or A it is rather limited. One has to vary it considerably to achieve a small change. This is analogous to the small eect o binder decomposition kinetics on steady-state burning rate. Increasing B makes the binder heat eedback smaller and thereore the AP response larger. A certain shit o the peak requency to the right is also seen in Fig. 7. On the other hand, increasing n decreases the AP response as more pressure-sensitive heating goes to the binder. Figures 9 and 1 are plots o the response unction or several C and C. As epected, the AP response is larger or less negative C and smaller C (less coupling rom or less coupling to the binder). No change in the orms o the curves can be observed in the igures. Again, we see the stabilizing inluence o the binder. Finally, Fig. 11 shows the inluence o the requency ratio. What is interesting is that it is requency dependent. Although or low requencies ( 1 < 5), the response is larger or lower, or high requencies ( 1 > ), it is larger or higher. There is an intermediate region in between. To understand the reason or this, one needs to observe that increasing or given 1 makes the binder requency higher. I is small, it is shited toward the peak area, and so the binder response becomes larger and the AP response correspondingly smaller. On the other hand, i is higher than the peak requency, then the binder response decreases and the AP response increases. Frequency ratio relects dierences in mean regression rates and thermal diusivities between AP and binder. Again, we see the binder being more or less a stabilizing inluence on the system. Fig. 6 Real part o AP response unction or several A. Fig. 7 Real part o AP response unction or several B. V. Intrinsic Instability It is well known that QSHOD-type models o homogeneous propellants have intrinsic stability limits, that is, there are regions o its parameters in which no steady burning is possible because the combustion is inherently unstable []. The goal o this section is to consider the possibility o intrinsic instability or the present model. We shall derive suicient conditions to prevent intrinsic instability and use the Cohen and Strand model [5] to see i they are satisied in practice. Fig. 8 Real part o AP response unction or several n.
6 SHUSSER, CULICK, AND COHEN 163 The solutions (15) and (16) can be written as R n B D C n B D D C C (1) Here, R n B D C n B D D C C () D =A 1= 1 B 1=A (3) D =A 1= 1 B 1=A (4) For intrinsic instability, the denominator in Eqs. (1) and () must vanish. This is equivalent to Fig. 9 Real part o AP response unction or several C. D Re D Re D Im D Im C C (5) where D Re D, D., D Re D Re D Im D Im D Re (6), D Im, D Im are the real and imaginary parts o We consider the case where both D Re Then, to satisy Eq. (6), we would need, D Re are positive: D Re > ; D Re > (7) D Im D Im < (8) Comparing Eqs. (5) and (8) one sees that i Eq. (7) is satisied then no intrinsic stability is possible or C C <.IC C > then, to prevent intrinsic instability, one needs to impose the additional condition D Re D Re >C C (9) Fig. 1 Real part o AP response unction or several C. In Sec. III, we obtained that, or the Cohen and Strand model [5], C C is almost always negative (see Figs. 4 and 5). This is consistent with the physics o the problem, as one would epect a destabilizing inluence o the AP on the binder (C > ) and a stabilizing inluence o the inert binder on the AP (C < ). However, or ine AP and low pressure, C C > was obtained. Whether this is an artiact o the model or relects the real physical situation is at present unclear. Nevertheless, we will consider both C C < and C C > in our analysis. To obtain conditions or the parameters, we write the real parts eplicitly: r Re B (3) D Re A Fig. 11 Real part o AP response unction or several. r Re B (31) D Re A Here, Re, Re are the real parts o the comple requencies, (see the Appendi). Equation (18) was used in obtaining Eqs. (3) and (31). Calculating the minimum values o D Re 1 Re < 1, one sees that Eq. (7) is equivalent to A, D Re or 1 Re < 1, s B 1 > 1 A (3)
7 164 SHUSSER, CULICK, AND COHEN n ln r ln p ; ln ln p (36) However, the situation or composite propellants turns out to be more complicated. In this section, we will demonstrate the eistence o two types o pressure eponents or a composite propellant and discuss the physical reasons or this phenomenon. A. Two Types o Pressure Eponent The classical pressure eponent or a homogeneous propellant n can be deined in two ways. Besides the usual deinition [Eq. (35)], one can use the surace boundary condition (8) to deine it as (37) z Fig. 1 s B 1 > 1 A (33) A Equations (3) and (33) are the suicient conditions to prevent intrinsic instability or C C <. IC C >, another condition must be added: B s A Suicient condition to prevent intrinsic instability. 1 1 B A s A 1 1 A >C C (34) To see i Eqs. (3 34) are satisied in practice, we use the values o the coeicients or the Cohen and Strand model [5] calculated in Sec. IV (see Figs. 5). One sees rom these igures that Eq. (33) is always satisied and Eq. (34) is satisied or C C >. However, this is not the case or condition (3). The let-hand side o condition (3) is plotted vs pressure or our particle sizes in Fig. 1. The igure shows that condition (3) is violated or high pressure and coarse AP. For m AP, the limit value o pressure is slightly above 16 atm. In numerical computations o composite propellant combustion [8], intrinsic instability was ound at 17 atm. Hence, the criterion o condition (3) is not only a suicient condition but probably close to the necessary condition as well. It should be stressed that, or high pressure and coarse AP, C and AP is practically decoupled rom the binder and behaves like pure AP. Only under these conditions, the Cohen and Strand model [5] allows intrinsic instability. One can thereore reach a tentative conclusion that the stabilizing inluence o an inert binder prevents an intrinsic instability under most conditions o interest. We show again the importance o avoiding coarse AP at high pressures, this time eplained with mathematics. VI. Pressure Eponents in Composite Propellants Pressure eponents represent the pressure dependence o the steady-state burning rate. For a homogeneous propellant, it is deined ln r ln p (35) Similarly, or a composite propellant, one can deine AP and binder pressure eponents that represent the pressure dependence o the steady-state burning rates r 1 and r : The deinitions o Eqs. (35) and (37) describe two aspects o the same phenomena. The ormer shows the pressure dependence o the steady-state burning rate, whereas the latter relates heat eedback and pressure oscillations. Burning rate is determined by surace temperature, which in turn is obtained rom the heat lu boundary condition at the solid surace. For a composite propellant, we have introduced two parameters n and n in boundary conditions (9) and (1). They can be deined by etending Eq. (37) to composite propellants 1 (38) 1 z1 1 r (39) z 1 r However, there are other ways, besides Eqs. (38) and (39), to etend Eq. (37) to composite propellants. The deinitions o Eqs. (38) and (39) assume that one calculates the p derivative while keeping both r1 and r at zero. However, it is not necessary to keep both o them at zero. We can etend Eq. (37) to a composite propellant while keeping each component s corresponding steady-state burning rate, that n 1 lim (4) 1 z1 p n z r 1 r (41) The limits in Eqs. (4) and (41) are needed in accordance with the steady-state limit. To calculate n 1 and n, we use the surace boundary conditions (9) and (1) and the response unction deinition Eq. (14), to obtain rom Eqs. (4) and (41) n 1 n C B R j ; n n C B R j (4) As will be proven later [see Eq. (5)], it ollows rom Eq. (4) that n 1 n C B n ; n n C B n 1 (43) Substituting into Eq. (4) the zero-requency limits o both response unctions that can be retrieved rom Eqs. (15) and (16), one obtains n 1 B n B C n (44) B B C C
8 SHUSSER, CULICK, AND COHEN 165 n B n B C n B B C C (45) B. Method o Homogenization To understand the meaning o n and n, we will use what one might call a homogenization o the propellant. Consider the AP surace boundary condition p B 1 p 1 B 1 r 1 r C A r (46) 1 r By deinition o the binder response unction R, r r R p p (47) Hence, Eq. (46) can be written n B 1 p 1 B 1 r 1 (48) A r 1 Fig. 13 AP pressure eponent in a propellant. where n n 1 C R n B (49) Equation (48) looks eactly like the surace boundary condition or a homogeneous propellant, Eq. (8). O course, it is an unusual homogeneous propellant because its pressure eponent n depends on the oscillatory requency through R. Nevertheless, the classical relationships remain valid or this propellant. For eample, its response unction R is given by n R A B A = 1 A 1 B (5) Because o n being a unction o R, Eq. (5) yields a relation between AP and binder response unctions R and R. One can use Eqs. (15) and (16) to veriy Eq. (5). The same method can also be applied to the binder component, or which n can be similarly deined. Thus, we can make an important conclusion that each component in the composite propellant behaves like a homogeneous propellant with a requency-dependent pressure eponent. The n and n, thereore, are pressure eponents characteristic o the nonsteady coupling. For C, n n n 1. For C, n n n ; the binder can still burn because o the steady-state heat eedback but there is no nonsteady interaction. The steady-state pressure eponents are obtained by taking the zero-requency limit o n and n, which yields Eq. (4) or n 1 and or n. One sees rom Eq. (43) that, without coupling, n 1 n i C ; n n i C (51) This is consistent with the result that, when one o the unsteady coupling coeicients vanishes, the corresponding response unction reduces to its homogeneous orm. It can also be veriied that, as or homogeneous propellant, the response unction value or zero requency is equal to the corresponding pressure eponent: R n 1 or ; R n or (5) The method o homogenization can be used to calculate the response unction o a multicomponent propellant (see [15]). Fig. 14 Binder pressure eponent in a propellant. C. Practical Values To estimate the physically reasonable values o the pressure eponents, we calculated them or the Cohen and Strand model [5]. As in Sec. III, AP concentration was kept at 87% throughout the calculations and propellant bulk temperature was 5 C. Both pressure eponents n 1 and n are plotted in Figs. 13 and 14. One can see that the dependencies o n 1 and n are more or less like n and n. However, n 1 <n because C < and n >n because C > [see Figs. 4 and 5 and Eq. (43)]. One can also deine the pressure eponent or the whole propellant: n ln r ln p The total propellant burning rate is given by (53) r prop S r 1 S r prop (54)
9 166 SHUSSER, CULICK, AND COHEN Here,,, and prop are AP, binder, and propellant densities; S and S are ractions o propellant surace area occupied by AP and binder, respectively. This also yields the propellant response unction rom the component contributions. Dierentiating Eq. (54) with respect to ln p and using Eq. (36), we obtain n prop n n (55) where is the raction o propellant mass low rate rom AP. The calculation shows [15] that n prop is high or ine and coarse particles but can be considerably smaller or intermediate sizes where diusion control is important in the lame. VII. Conclusions We have obtained a general analytical solution or the pressurecoupled response o a monomodal composite propellant, applied to AP and inert binder. The solution depends on the AP coeicients A, B, and n, the binder coeicients A, B, and n, the coupling coeicients C and C, and the requency ratio. The inluence o these parameters on the real part o the AP response unction in the propellant has been investigated. The inluence o the AP coeicients is analogous to the A, B, and n parameters o the classical homogeneous propellant solution. The increase in A increases the response and moves its maimum to higher requencies. Decreasing B makes the response unction curve more peaked, whereas larger n increases the value o the response but not its requency dependence. The inluence o A was ound to be limited due to the small eect o binder decomposition kinetics on burning rate. Increasing B makes the binder heat eedback smaller and thereore the AP response larger, whereas increasing n decreases the AP response as more heat goes to the binder. Less coupling rom the binder (less negative C ) or less coupling to the binder (smaller C ) increases the response. The inluence o the requency ratio was ound to be requency dependent. The parameters A, A, B, B, n, n, C, and C were calculated or a particular case o the Cohen and Strand model [5]. The analytical solution was shown to be in ecellent agreement with numerical computation results, and has the advantage o revealing important mechanisms. Thus, we see the destabilizing inluence o the monopropellant AP and the stabilizing inluence o the inert binder in the system. The inluence o absolute pressure and AP particle size on the parameters has been studied. B in a propellant turned out to be much higher than in the classical ormulation due to the stabilizing inluence o the binder in the coupling. B decreases with increasing pressure mainly because o decreasing lame heights. The B parameter was ound to be even larger than B. This shows that an inert binder is highly stabilizing to the propellant system. Because o the binder being heated only by the diusion lame, n decreases much aster than n with increasing pressure. The AP coupling coeicient C is always negative. This relects the stabilizing inluence o the binder on the AP by the coupling. However, C becomes less negative with increasing pressure and particle size. C vanishes or high pressures and large sizes, where the AP becomes decoupled rom the binder and behaves more like pure AP. The binder coupling coeicient C is generally positive, showing that the AP has a destabilizing inluence on the binder. The AP monopropellant combustion is the destabilizing inluence on the system, whereas the diusion lame with binder is stabilizing. Suicient conditions preventing intrinsic instability in a composite propellant were obtained. Using the Cohen and Strand model [5], we demonstrated that these conditions are satisied in most practical situations. In agreement with previous numerical results, it was ound mathematically that intrinsic instability is possible only at high pressures and with coarse AP. At this condition, AP is practically decoupled rom the binder and behaves like pure AP. Thus, it appears that in most practical cases, the stabilizing inluence o an inert binder prevents intrinsic instability. To show this, it is necessary to treat the two-component system as it really is rather than as an averaged one-component system. We have shown that, or a two-component composite propellant, two orms o pressure eponents arise rom the analysis. Whereas the conventional n 1 and n show the pressure dependence o both AP and binder steady-state burning rates, n and n arise rom the oscillations o the unsteady heat eedback and unsteady coupling. The general epressions or n 1 and n in propellants have been obtained and their dependencies on pressure and AP particle size were studied or the Cohen and Strand model [5]. In addition, the pressure eponent o the aggregate composite propellant was deined. We also demonstrated that both AP and binder in the composite propellant system behave like a homogeneous propellant with a requency-dependent pressure eponent characteristic o the nonsteady coupling. The convenience o an analytical solution will acilitate urther studies o the response unction o composite propellants. Other advanced energetic ingredients may be more destabilizing than AP and can be evaluated or linearized analysis by studying the monopropellants themselves in the ramework o this theoretical approach, which began with the study o AP [14]. Appendi: Real and Imaginary Parts o Binder and Ammonium Perchlorate Response Functions The real and imaginary parts o the binder and AP response unctions R and R are as ollows: Here, R e 1e 5 e e 6 e 5 i e e 5 e 1 e 6 e 6 e 5 e 6 R e 3e 5 e 4 e 6 e 5 i e 4e 5 e 3 e 6 e 6 e 5 e 6 (A1) (A) e 1 n B Re A Re 1 B 1A n B C Im (A3) e n B Im A Re Im e 3 n B r A Re Im e 4 n B Im A Re Im (A4) 1 B 1A n B C e 5 Re A Re 1 B 1A Im Re A Re 1 B 1A Im Im Im A Re Im A Re C C Im (A5) (A6) (A7)
10 SHUSSER, CULICK, AND COHEN 167 e 6 Im A Re Im Re A Re 1 B 1A Im Im A Re Im Re 1 B A Re 1A (A8) Im The real and imaginary parts o the comple requencies and are given by Re 1 p = 1 (A9) Im 1 p Re 1 1 1= p 1 16 = Im 1 p = (A1) (A11) (A1) Acknowledgments This work was sponsored partly by the Caliornia Institute o Technology and partly by the U.S. Oice o Naval Research Multidisciplinary University Research Initiative under Grant No. N14-95-I Judah Goldwasser was the Navy Program Manager. The supplemental support under Multidisciplinary University Research Initiative Grant No. N14-95-I-1339, through the courtesy o Herman Krier, University o Illinois at Urbana-Champaign, is grateully acknowledged. Reerences [1] Price, E. W., and Flandro, G. A., Status and Prospects or Future Developments, Nonsteady Burning and Combustion Stability o Solid Propellants, edited by L. De Luca, E. W. Price, and M. Summerield, Vol. 143, Progress in Astronautics and Aeronautics, AIAA, Washington, D.C., 199, pp [] Culick, F. E. C., Review o Calculations or Unsteady Burning o a Solid Propellant, AIAA Journal, Vol. 6, No. 1, 1968, pp doi:1.514/3.498 [3] Nachbar, W., Theoretical Study o the Burning o a Solid Propellant Sandwich, Solid Propellant Rocket Research, edited by M. Summerield, Vol. 1, ARS Progress in Astronautics and Rocketry, Academic Press, New York, 196, pp [4] Price, E. W., Sambamurthi, J. K., Sigman, R. K., and Panyam, R. R., Combustion o Ammonium Perchlorate Polymer Sandwiches, Combustion and Flame, Vol. 63, No. 3, 1986, pp doi:1.116/1-18(86)97-6 [5] Cohen, N. S., and Strand, L. D., Improved Model or the Combustion o AP Composite Propellants, AIAA Journal, Vol., No. 1, 198, pp doi:1.514/3.813 [6] Cohen, N. S., and Strand, L. D., Combustion Response to Compositional Fluctuations, AIAA Journal, Vol. 3, No. 5, 1985, pp doi:1.514/ [7] Cohen, N. S., and Strand, L. D., Eect o AP Particle Size on Combustion Response to Crosslow, AIAA Journal, Vol. 3, No. 5, 1985, pp doi:1.514/ [8] Shusser, M., Culick, F. E. C., and Cohen, N. S., Combustion Response o Ammonium Perchlorate Composite Propellants, Journal o Propulsion and Power, Vol. 18, No. 5,, pp [9] Rasmussen, B., and Frederick, R. A., Jr., Nonlinear Hterogeneous Model o Composite Solid Propellant Combustion, Journal o Propulsion and Power, Vol. 18, No. 5,, pp [1] Surzhikov, S. T., Murphy, J. J., and Krier, H., Two-Dimensional Model or Unsteady Burning Heterogeneous AP/Binder Solid Propellants, AIAA Paper -3573, July. [11] Jackson, T. L., and Buckmaster, J., Heterogeneous Propellant Combustion, AIAA Journal, Vol. 4, No. 6,, pp [1] Massa, L., Jackson, T. L., Buckmaster, J., and Campbell, M., Three- Dimensional Heterogeneous Propellant Combustion, Proceedings o the Combustion Institute, Vol. 9, No.,, pp doi:1.116/s () [13] Brown, R. S., and Muzzy, R. J., Linear and Nonlinear Pressure Coupled Combustion Instability o Solid Propellants, AIAA Journal, Vol. 8, No. 8, 197, pp doi:1.514/3.599 [14] Shusser, M., Culick, F. E. C., and Cohen, N. S., Combustion Response o Ammonium Perchlorate, AIAA Journal, Vol. 4, No. 4,, pp [15] Shusser, M., and Culick, F. E. C., Analytical Solution or Composite Solid Propellant Response Function, Guggenheim Jet Propulsion Center, Caliornia Inst. o Technology, Rept. CI 1-, May 1. [16] Beckstead, M. W., Derr, R. L., and Price, C. F., Model o Composite Solid-Propellant Combustion Based on Multiple Flames, AIAA Journal, Vol. 8, No. 1, 197, pp. 7. doi:1.514/3.687 V. Yang Editor-in-Chie
RATIONAL FUNCTIONS. Finding Asymptotes..347 The Domain Finding Intercepts Graphing Rational Functions
RATIONAL FUNCTIONS Finding Asymptotes..347 The Domain....350 Finding Intercepts.....35 Graphing Rational Functions... 35 345 Objectives The ollowing is a list o objectives or this section o the workbook.
More informationNumerical Methods - Lecture 2. Numerical Methods. Lecture 2. Analysis of errors in numerical methods
Numerical Methods - Lecture 1 Numerical Methods Lecture. Analysis o errors in numerical methods Numerical Methods - Lecture Why represent numbers in loating point ormat? Eample 1. How a number 56.78 can
More information( 1) ( 2) ( 1) nan integer, since the potential is no longer simple harmonic.
. Anharmonic Oscillators Michael Fowler Landau (para 8) considers a simple harmonic oscillator with added small potential energy terms mα + mβ. We ll simpliy slightly by dropping the term, to give an equation
More informationObjectives. By the time the student is finished with this section of the workbook, he/she should be able
FUNCTIONS Quadratic Functions......8 Absolute Value Functions.....48 Translations o Functions..57 Radical Functions...61 Eponential Functions...7 Logarithmic Functions......8 Cubic Functions......91 Piece-Wise
More informationAcoustic forcing of flexural waves and acoustic fields for a thin plate in a fluid
Acoustic orcing o leural waves and acoustic ields or a thin plate in a luid Darryl MCMAHON Maritime Division, Deence Science and Technology Organisation, HMAS Stirling, WA Australia ABSTACT Consistency
More informationBasic mathematics of economic models. 3. Maximization
John Riley 1 January 16 Basic mathematics o economic models 3 Maimization 31 Single variable maimization 1 3 Multi variable maimization 6 33 Concave unctions 9 34 Maimization with non-negativity constraints
More informationThe concept of limit
Roberto s Notes on Dierential Calculus Chapter 1: Limits and continuity Section 1 The concept o limit What you need to know already: All basic concepts about unctions. What you can learn here: What limits
More informationChapter 6 Reliability-based design and code developments
Chapter 6 Reliability-based design and code developments 6. General Reliability technology has become a powerul tool or the design engineer and is widely employed in practice. Structural reliability analysis
More informationAnalog Computing Technique
Analog Computing Technique by obert Paz Chapter Programming Principles and Techniques. Analog Computers and Simulation An analog computer can be used to solve various types o problems. It solves them in
More informationWe would now like to turn our attention to a specific family of functions, the one to one functions.
9.6 Inverse Functions We would now like to turn our attention to a speciic amily o unctions, the one to one unctions. Deinition: One to One unction ( a) (b A unction is called - i, or any a and b in the
More informationFluctuationlessness Theorem and its Application to Boundary Value Problems of ODEs
Fluctuationlessness Theorem and its Application to Boundary Value Problems o ODEs NEJLA ALTAY İstanbul Technical University Inormatics Institute Maslak, 34469, İstanbul TÜRKİYE TURKEY) nejla@be.itu.edu.tr
More informationBayesian Technique for Reducing Uncertainty in Fatigue Failure Model
9IDM- Bayesian Technique or Reducing Uncertainty in Fatigue Failure Model Sriram Pattabhiraman and Nam H. Kim University o Florida, Gainesville, FL, 36 Copyright 8 SAE International ABSTRACT In this paper,
More informationCONVECTIVE HEAT TRANSFER CHARACTERISTICS OF NANOFLUIDS. Convective heat transfer analysis of nanofluid flowing inside a
Chapter 4 CONVECTIVE HEAT TRANSFER CHARACTERISTICS OF NANOFLUIDS Convective heat transer analysis o nanoluid lowing inside a straight tube o circular cross-section under laminar and turbulent conditions
More information9.3 Graphing Functions by Plotting Points, The Domain and Range of Functions
9. Graphing Functions by Plotting Points, The Domain and Range o Functions Now that we have a basic idea o what unctions are and how to deal with them, we would like to start talking about the graph o
More informationRelating axial motion of optical elements to focal shift
Relating aial motion o optical elements to ocal shit Katie Schwertz and J. H. Burge College o Optical Sciences, University o Arizona, Tucson AZ 857, USA katie.schwertz@gmail.com ABSTRACT In this paper,
More informationMath 2412 Activity 1(Due by EOC Sep. 17)
Math 4 Activity (Due by EOC Sep. 7) Determine whether each relation is a unction.(indicate why or why not.) Find the domain and range o each relation.. 4,5, 6,7, 8,8. 5,6, 5,7, 6,6, 6,7 Determine whether
More informationBasic properties of limits
Roberto s Notes on Dierential Calculus Chapter : Limits and continuity Section Basic properties o its What you need to know already: The basic concepts, notation and terminology related to its. What you
More information( x) f = where P and Q are polynomials.
9.8 Graphing Rational Functions Lets begin with a deinition. Deinition: Rational Function A rational unction is a unction o the orm ( ) ( ) ( ) P where P and Q are polynomials. Q An eample o a simple rational
More informationSection 3.4: Concavity and the second Derivative Test. Find any points of inflection of the graph of a function.
Unit 3: Applications o Dierentiation Section 3.4: Concavity and the second Derivative Test Determine intervals on which a unction is concave upward or concave downward. Find any points o inlection o the
More informationGap Excitations and Series Loads in Microstrip Lines: Equivalent Network Characterization with Application to THz Circuits
Gap Ecitations and Series Loads in Microstrip Lines: Equivalent Network Characterization with Application to THz Circuits Andrea Neto and Peter H. Siegel Jet Propulsion Laboratory, Caliornia Institute
More informationRelating axial motion of optical elements to focal shift
Relating aial motion o optical elements to ocal shit Katie Schwertz and J. H. Burge College o Optical Sciences, University o Arizona, Tucson AZ 857, USA katie.schwertz@gmail.com ABSTRACT In this paper,
More informationInternal thermal noise in the LIGO test masses: A direct approach
PHYSICAL EVIEW D VOLUME 57, NUMBE 2 15 JANUAY 1998 Internal thermal noise in the LIGO test masses: A direct approach Yu. Levin Theoretical Astrophysics, Caliornia Institute o Technology, Pasadena, Caliornia
More informationPremixed Flame Response to Acoustic Waves in a Porous-Walled Chamber with Surface Mass Injection. Propulsion Conference and Exhibit
AIAA -369 Premixed Flame Response to Acoustic Waves in a Porous-Walled Chamber with Surace Mass Injection Wen-Wei Chu and Vigor Yang The Pennsylvania State University University Park, PA 68 Propulsion
More informationApplication of Method of Lines on the Bifurcation Analysis of a Horizontal Rijke Tube
th ASPACC July 9, 5 Beijing, China Application o Method o Lines on the Biurcation Analysis o a Horizontal Rijke Tube Xiaochuan Yang, Ali Turan, Adel Nasser, Shenghui Lei School o Mechanical, Aerospace
More informationTLT-5200/5206 COMMUNICATION THEORY, Exercise 3, Fall TLT-5200/5206 COMMUNICATION THEORY, Exercise 3, Fall Problem 1.
TLT-5/56 COMMUNICATION THEORY, Exercise 3, Fall Problem. The "random walk" was modelled as a random sequence [ n] where W[i] are binary i.i.d. random variables with P[W[i] = s] = p (orward step with probability
More informationMicro-canonical ensemble model of particles obeying Bose-Einstein and Fermi-Dirac statistics
Indian Journal o Pure & Applied Physics Vol. 4, October 004, pp. 749-757 Micro-canonical ensemble model o particles obeying Bose-Einstein and Fermi-Dirac statistics Y K Ayodo, K M Khanna & T W Sakwa Department
More informationAnalysis of Foundation Damping from Theoretical Models and Forced Vibration Testing
6 th International Conerence on Earthquae Geotechnical Engineering 1-4 November 2015 Christchurch, New Zealand Analysis o Foundation Damping rom Theoretical Models and Forced Vibration Testing M. J. Givens
More informationAP Calculus Notes: Unit 1 Limits & Continuity. Syllabus Objective: 1.1 The student will calculate limits using the basic limit theorems.
Syllabus Objective:. The student will calculate its using the basic it theorems. LIMITS how the outputs o a unction behave as the inputs approach some value Finding a Limit Notation: The it as approaches
More informationAggregate Growth: R =αn 1/ d f
Aggregate Growth: Mass-ractal aggregates are partly described by the mass-ractal dimension, d, that deines the relationship between size and mass, R =αn 1/ d where α is the lacunarity constant, R is the
More informationIntrinsic Small-Signal Equivalent Circuit of GaAs MESFET s
Intrinsic Small-Signal Equivalent Circuit o GaAs MESFET s M KAMECHE *, M FEHAM M MELIANI, N BENAHMED, S DALI * National Centre o Space Techniques, Algeria Telecom Laboratory, University o Tlemcen, Algeria
More informationScattered Data Approximation of Noisy Data via Iterated Moving Least Squares
Scattered Data Approximation o Noisy Data via Iterated Moving Least Squares Gregory E. Fasshauer and Jack G. Zhang Abstract. In this paper we ocus on two methods or multivariate approximation problems
More informationSyllabus Objective: 2.9 The student will sketch the graph of a polynomial, radical, or rational function.
Precalculus Notes: Unit Polynomial Functions Syllabus Objective:.9 The student will sketch the graph o a polynomial, radical, or rational unction. Polynomial Function: a unction that can be written in
More informationCombustion Instabilities and Oscillations in Solid Propellant Rockets
Combustion Instabilities and Oscillations in Solid Propellant Rockets F.E.C. Culick California Institute of Technology Joint Propulsion Conference 12 July 2006 http://www.culick.caltech.edu 1 Outline Introduction
More informationNumerical Solution of Ordinary Differential Equations in Fluctuationlessness Theorem Perspective
Numerical Solution o Ordinary Dierential Equations in Fluctuationlessness Theorem Perspective NEJLA ALTAY Bahçeşehir University Faculty o Arts and Sciences Beşiktaş, İstanbul TÜRKİYE TURKEY METİN DEMİRALP
More informationRepresentation of Coriolis forces and simulation procedures for moving fluid-conveying pipes
Representation o Coriolis orces and simulation procedures or moving luid-conveying pipes Jörg Gebhardt*, Frank Kassubek** *ABB Corporate Research Center Germany, Department or Automation Device Technologies,
More information8.4 Inverse Functions
Section 8. Inverse Functions 803 8. Inverse Functions As we saw in the last section, in order to solve application problems involving eponential unctions, we will need to be able to solve eponential equations
More information18-660: Numerical Methods for Engineering Design and Optimization
8-66: Numerical Methods or Engineering Design and Optimization Xin Li Department o ECE Carnegie Mellon University Pittsburgh, PA 53 Slide Overview Linear Regression Ordinary least-squares regression Minima
More information. This is the Basic Chain Rule. x dt y dt z dt Chain Rule in this context.
Math 18.0A Gradients, Chain Rule, Implicit Dierentiation, igher Order Derivatives These notes ocus on our things: (a) the application o gradients to ind normal vectors to curves suraces; (b) the generaliation
More informationy2 = 0. Show that u = e2xsin(2y) satisfies Laplace's equation.
Review 1 1) State the largest possible domain o deinition or the unction (, ) = 3 - ) Determine the largest set o points in the -plane on which (, ) = sin-1( - ) deines a continuous unction 3) Find the
More informationAn Alternative Poincaré Section for Steady-State Responses and Bifurcations of a Duffing-Van der Pol Oscillator
An Alternative Poincaré Section or Steady-State Responses and Biurcations o a Duing-Van der Pol Oscillator Jang-Der Jeng, Yuan Kang *, Yeon-Pun Chang Department o Mechanical Engineering, National United
More informationDifferentiation. The main problem of differential calculus deals with finding the slope of the tangent line at a point on a curve.
Dierentiation The main problem o dierential calculus deals with inding the slope o the tangent line at a point on a curve. deinition() : The slope o a curve at a point p is the slope, i it eists, o the
More informationA Simple Explanation of the Sobolev Gradient Method
A Simple Explanation o the Sobolev Gradient Method R. J. Renka July 3, 2006 Abstract We have observed that the term Sobolev gradient is used more oten than it is understood. Also, the term is oten used
More informationTransient Adsorption and Desorption in Micrometer Scale Features
Journal o The Electrochemical Society, 149 8 G461-G473 2002 0013-4651/2002/149 8 /G461/13/$7.00 The Electrochemical Society, Inc. Transient Adsorption and Desorption in Micrometer Scale Features Matthias
More informationOPTIMAL PLACEMENT AND UTILIZATION OF PHASOR MEASUREMENTS FOR STATE ESTIMATION
OPTIMAL PLACEMENT AND UTILIZATION OF PHASOR MEASUREMENTS FOR STATE ESTIMATION Xu Bei, Yeo Jun Yoon and Ali Abur Teas A&M University College Station, Teas, U.S.A. abur@ee.tamu.edu Abstract This paper presents
More informationLecture Outline. Basics of Spatial Filtering Smoothing Spatial Filters. Sharpening Spatial Filters
1 Lecture Outline Basics o Spatial Filtering Smoothing Spatial Filters Averaging ilters Order-Statistics ilters Sharpening Spatial Filters Laplacian ilters High-boost ilters Gradient Masks Combining Spatial
More informationCHAPTER 4 Reactor Statics. Table of Contents
CHAPTER 4 Reactor Statics Prepared by Dr. Benjamin Rouben, & Consulting, Adjunct Proessor, McMaster University & University o Ontario Institute o Technology (UOIT) and Dr. Eleodor Nichita, Associate Proessor,
More informationwhose domain D is a set of n-tuples in is defined. The range of f is the set of all values f x1,..., x n
Grade (MCV4UE) - AP Calculus Etended Page o A unction o n-variales is a real-valued unction... n whose domain D is a set o n-tuples... n in which... n is deined. The range o is the set o all values...
More informationComments on Problems. 3.1 This problem offers some practice in deriving utility functions from indifference curve specifications.
CHAPTER 3 PREFERENCES AND UTILITY These problems provide some practice in eamining utilit unctions b looking at indierence curve maps and at a ew unctional orms. The primar ocus is on illustrating the
More information8. THEOREM If the partial derivatives f x. and f y exist near (a, b) and are continuous at (a, b), then f is differentiable at (a, b).
8. THEOREM I the partial derivatives and eist near (a b) and are continuous at (a b) then is dierentiable at (a b). For a dierentiable unction o two variables z= ( ) we deine the dierentials d and d to
More informationWind-Driven Circulation: Stommel s gyre & Sverdrup s balance
Wind-Driven Circulation: Stommel s gyre & Sverdrup s balance We begin by returning to our system o equations or low o a layer o uniorm density on a rotating earth. du dv h + [ u( H + h)] + [ v( H t y d
More informationOptimal Control of process
VYSOKÁ ŠKOLA BÁŇSKÁ TECHNICKÁ UNIVERZITA OSTRAVA FAKULTA METALURGIE A MATERIÁLOVÉHO INŽENÝRSTVÍ Optimal Control o process Study Support Milan Heger Ostrava 8 Title: Optimal Control o process Code: Author:
More informationPart I: Thin Converging Lens
Laboratory 1 PHY431 Fall 011 Part I: Thin Converging Lens This eperiment is a classic eercise in geometric optics. The goal is to measure the radius o curvature and ocal length o a single converging lens
More informationQuality control of risk measures: backtesting VAR models
De la Pena Q 9/2/06 :57 pm Page 39 Journal o Risk (39 54 Volume 9/Number 2, Winter 2006/07 Quality control o risk measures: backtesting VAR models Victor H. de la Pena* Department o Statistics, Columbia
More informationCHAPTER 1: INTRODUCTION. 1.1 Inverse Theory: What It Is and What It Does
Geosciences 567: CHAPTER (RR/GZ) CHAPTER : INTRODUCTION Inverse Theory: What It Is and What It Does Inverse theory, at least as I choose to deine it, is the ine art o estimating model parameters rom data
More informationFeasibility of a Multi-Pass Thomson Scattering System with Confocal Spherical Mirrors
Plasma and Fusion Research: Letters Volume 5, 044 200) Feasibility o a Multi-Pass Thomson Scattering System with Conocal Spherical Mirrors Junichi HIRATSUKA, Akira EJIRI, Yuichi TAKASE and Takashi YAMAGUCHI
More informationy,z the subscript y, z indicating that the variables y and z are kept constant. The second partial differential with respect to x is written x 2 y,z
8 Partial dierentials I a unction depends on more than one variable, its rate o change with respect to one o the variables can be determined keeping the others ied The dierential is then a partial dierential
More informationChapter 2 Section 3. Partial Derivatives
Chapter Section 3 Partial Derivatives Deinition. Let be a unction o two variables and. The partial derivative o with respect to is the unction, denoted b D1 1 such that its value at an point (,) in the
More informationOn a Closed Formula for the Derivatives of e f(x) and Related Financial Applications
International Mathematical Forum, 4, 9, no. 9, 41-47 On a Closed Formula or the Derivatives o e x) and Related Financial Applications Konstantinos Draais 1 UCD CASL, University College Dublin, Ireland
More informationAnalysis of Friction-Induced Vibration Leading to Eek Noise in a Dry Friction Clutch. Abstract
The 22 International Congress and Exposition on Noise Control Engineering Dearborn, MI, USA. August 19-21, 22 Analysis o Friction-Induced Vibration Leading to Eek Noise in a Dry Friction Clutch P. Wickramarachi
More informationAP Calculus AB Summer Assignment
AP Calculus AB Summer Assignment As Advanced placement students, our irst assignment or the 07-08 school ear is to come to class the ver irst da in top mathematical orm. Calculus is a world o change. While
More information39.1 Gradually Varied Unsteady Flow
39.1 Gradually Varied Unsteady Flow Gradually varied unsteady low occurs when the low variables such as the low depth and velocity do not change rapidly in time and space. Such lows are very common in
More information2. ETA EVALUATIONS USING WEBER FUNCTIONS. Introduction
. ETA EVALUATIONS USING WEBER FUNCTIONS Introduction So ar we have seen some o the methods or providing eta evaluations that appear in the literature and we have seen some o the interesting properties
More information(One Dimension) Problem: for a function f(x), find x 0 such that f(x 0 ) = 0. f(x)
Solving Nonlinear Equations & Optimization One Dimension Problem: or a unction, ind 0 such that 0 = 0. 0 One Root: The Bisection Method This one s guaranteed to converge at least to a singularity, i not
More informationIncreasing and Decreasing Functions and the First Derivative Test. Increasing and Decreasing Functions. Video
SECTION and Decreasing Functions and the First Derivative Test 79 Section and Decreasing Functions and the First Derivative Test Determine intervals on which a unction is increasing or decreasing Appl
More informationOn the Girth of (3,L) Quasi-Cyclic LDPC Codes based on Complete Protographs
On the Girth o (3,L) Quasi-Cyclic LDPC Codes based on Complete Protographs Sudarsan V S Ranganathan, Dariush Divsalar and Richard D Wesel Department o Electrical Engineering, University o Caliornia, Los
More informationNew Functions from Old Functions
.3 New Functions rom Old Functions In this section we start with the basic unctions we discussed in Section. and obtain new unctions b shiting, stretching, and relecting their graphs. We also show how
More informationMODULE 6 LECTURE NOTES 1 REVIEW OF PROBABILITY THEORY. Most water resources decision problems face the risk of uncertainty mainly because of the
MODULE 6 LECTURE NOTES REVIEW OF PROBABILITY THEORY INTRODUCTION Most water resources decision problems ace the risk o uncertainty mainly because o the randomness o the variables that inluence the perormance
More informationINTRODUCTORY MATHEMATICAL ANALYSIS
INTRODUCTORY MATHEMATICAL ANALYSIS For Business, Economics, and the Lie and Social Sciences Chapter 11 Dierentiation 011 Pearson Education, Inc. Chapter 11: Dierentiation Chapter Objectives To compute
More informationPartially fluidized shear granular flows: Continuum theory and molecular dynamics simulations
Partially luidized shear granular lows: Continuum theory and molecular dynamics simulations Dmitri Volson, 1 Lev S. Tsimring, 1 and Igor S. Aranson 2 1 Institute or Nonlinear Science, University o Caliornia,
More informationOBSERVER/KALMAN AND SUBSPACE IDENTIFICATION OF THE UBC BENCHMARK STRUCTURAL MODEL
OBSERVER/KALMAN AND SUBSPACE IDENTIFICATION OF THE UBC BENCHMARK STRUCTURAL MODEL Dionisio Bernal, Burcu Gunes Associate Proessor, Graduate Student Department o Civil and Environmental Engineering, 7 Snell
More informationBOUNDARY LAYER ANALYSIS ALONG A STRETCHING WEDGE SURFACE WITH MAGNETIC FIELD IN A NANOFLUID
Proceedings o the International Conerence on Mechanical Engineering and Reneable Energy 7 (ICMERE7) 8 December, 7, Chittagong, Bangladesh ICMERE7-PI- BOUNDARY LAYER ANALYSIS ALONG A STRETCHING WEDGE SURFACE
More informationPressure derivatives of bulk modulus for materials at extreme compression
Indian Journal o ure & Applied hysics Vol. 5, October, pp. 734-738 ressure derivatives o bulk modulus or materials at extreme compression Singh* & A Dwivedi Department o hysics, Institute o Basic Sciences,
More informationA Method for Assimilating Lagrangian Data into a Shallow-Water-Equation Ocean Model
APRIL 2006 S A L M A N E T A L. 1081 A Method or Assimilating Lagrangian Data into a Shallow-Water-Equation Ocean Model H. SALMAN, L.KUZNETSOV, AND C. K. R. T. JONES Department o Mathematics, University
More informationTHE EFFECT OF THERMAL RADIATION ON THE DYNAMICS OF FLASHOVER IN A COMPARTMENT FIRE
The 6th ASME-JSME Thermal Engineering Joint Conerence March 16-20, 2003 TED-AJ03-127 THE EFFECT OF THERMAL RADIATION ON THE DYNAMICS OF FLASHOVER IN A COMPARTMENT FIRE W. W. YUEN Department o Mechanical
More informationEDGES AND CONTOURS(1)
KOM31 Image Processing in Industrial Sstems Dr Muharrem Mercimek 1 EDGES AND CONTOURS1) KOM31 Image Processing in Industrial Sstems Some o the contents are adopted rom R. C. Gonzalez, R. E. Woods, Digital
More informationLab on Taylor Polynomials. This Lab is accompanied by an Answer Sheet that you are to complete and turn in to your instructor.
Lab on Taylor Polynomials This Lab is accompanied by an Answer Sheet that you are to complete and turn in to your instructor. In this Lab we will approimate complicated unctions by simple unctions. The
More informationStick-slip conditions in the general motion of a planar rigid body
Journal o Mechanical Science and Technology 27 (9) (2013) 2577~2583 www.springerlink.com/content/1738-494 DOI 10.1007/s12206-013-0701-y Stick-slip conditions in the general motion o a planar rigid body
More informationReliability of Axially Loaded Fiber-Reinforced-Polymer Confined Reinforced Concrete Circular Columns
American J. o Engineering and Applied Sciences (1): 31-38, 009 ISSN 1941-700 009 Science Publications Reliability o Axially Loaded Fiber-Reinorced-Polymer Conined Reinorced Concrete Circular Columns Venkatarman
More informationNon-linear Extrapolation of Laminar Flame Properties from Spherically Expanding Flames
Paper # 87LF- Western States Section o the Combustion Institute Caliornia Institute o Technology March -5, 4. Topic: Laminar Flame Non-linear Extrapolation o Laminar Flame Properties rom Spherically Expanding
More informationSIMULATION OF A FLOW THROUGH A RIGID POROUS MEDIUM SUBJECTED TO A CONSTRAINT
SIMULATION O A LOW THROUGH A RIGID POROUS MEDIUM SUBJECTED TO A CONSTRAINT R. M. Saldanha da Gama a, J. J. Pedrosa ilho a, and M. L. Martins-Costa b a Universidade do Estado do Rio de Janeiro Departamento
More information«Develop a better understanding on Partial fractions»
«Develop a better understanding on Partial ractions» ackground inormation: The topic on Partial ractions or decomposing actions is irst introduced in O level dditional Mathematics with its applications
More informationFatigue verification of high loaded bolts of a rocket combustion chamber.
Fatigue veriication o high loaded bolts o a rocket combustion chamber. Marcus Lehmann 1 & Dieter Hummel 1 1 Airbus Deence and Space, Munich Zusammenassung Rocket engines withstand intense thermal and structural
More informationDiffusion in generalized lattice-gas models
PHYSICAL REVIEW B VOLUME 58, NUMBER Diusion in generalized lattice-gas models O. M. Braun Institute o Physics, Ukrainian Academy o Sciences, UA-50 Kiev, Ukraine DECEMBER 998-II C. A. Sholl Department o
More informationThe Deutsch-Jozsa Problem: De-quantization and entanglement
The Deutsch-Jozsa Problem: De-quantization and entanglement Alastair A. Abbott Department o Computer Science University o Auckland, New Zealand May 31, 009 Abstract The Deustch-Jozsa problem is one o the
More informationFRITS VEERMAN, ARJEN DOELMAN
PULSES IN A GIERER-MEINHARDT EQUATION WITH A SLOW NONLINEARITY FRITS VEERMAN, ARJEN DOELMAN Abstract. In this paper, we study in detail the existence and stability o localized pulses in a Gierer- Meinhardt
More informationA Fourier Transform Model in Excel #1
A Fourier Transorm Model in Ecel # -This is a tutorial about the implementation o a Fourier transorm in Ecel. This irst part goes over adjustments in the general Fourier transorm ormula to be applicable
More informationQuadratic Functions. The graph of the function shifts right 3. The graph of the function shifts left 3.
Quadratic Functions The translation o a unction is simpl the shiting o a unction. In this section, or the most part, we will be graphing various unctions b means o shiting the parent unction. We will go
More informationCOMBINED EFFECTS OF RADIATION AND HEAT GENERATION ON MHD NATURAL CONVECTION FLOW ALONG A VERTICAL FLAT PLATE IN PRESENCE OF HEAT CONDUCTION
BRAC University Journal, vol.vi, no., 9, pp 11- COMBINED EFFECTS OF RADIATION AND HEAT GENERATION ON MHD NATURAL CONVECTION FLOW ALONG A VERTICAL FLAT PLATE IN PRESENCE OF HEAT CONDUCTION Mohammad Mokaddes
More informationCHAPTER 8 ANALYSIS OF AVERAGE SQUARED DIFFERENCE SURFACES
CAPTER 8 ANALYSS O AVERAGE SQUARED DERENCE SURACES n Chapters 5, 6, and 7, the Spectral it algorithm was used to estimate both scatterer size and total attenuation rom the backscattered waveorms by minimizing
More informationHYDROMAGNETIC DIVERGENT CHANNEL FLOW OF A VISCO- ELASTIC ELECTRICALLY CONDUCTING FLUID
Rita Choudhury et al. / International Journal o Engineering Science and Technology (IJEST) HYDROAGNETIC DIVERGENT CHANNEL FLOW OF A VISCO- ELASTIC ELECTRICALLY CONDUCTING FLUID RITA CHOUDHURY Department
More informationCHAPTER 5 Reactor Dynamics. Table of Contents
1 CHAPTER 5 Reactor Dynamics prepared by Eleodor Nichita, UOIT and Benjamin Rouben, 1 & 1 Consulting, Adjunct Proessor, McMaster & UOIT Summary: This chapter addresses the time-dependent behaviour o nuclear
More informationA FORMULA FOR THE MEAN CURVATURE OF AN IMPLICIT REGULAR SURFACE
STUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume XLVI, Number 3, September 001 A FORMULA FOR THE MEAN CURVATURE OF AN IMPLICIT REGULAR SURFACE Abstract. In this paper we will ind a ormula or the absolute
More information39th AIAA/ASME/SAE/ASEE Joint Propulsion Conference July 2003 / Huntsville, AL
39th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit 20-23 July 2003, Huntsville, Alabama AIAA 2003-5082 AIAA 2003-5082 Direct Numerical Simulation of Sandwich and Random-Packed Propellant Combustion
More informationANALYSIS OF SPEED-DENSITY TRAFFIC FLOW MODELS ON A MERGE INFLUENCE SECTION IN AN UNINTERRUPTED FACILITY
ANALYSIS OF SPEED-DENSITY TRAFFIC FLOW MODELS ON A MERGE INFLUENCE SECTION IN AN UNINTERRUPTED FACILITY Tcheolwoong DOH Proessor Transportation Engineering University o Hanyang Sa-3dong, Sangrok-gu, Ansan
More informationCOMPUTATIONAL STUDY OF CHEMICALLY REACTING HYPERSONIC FLOW
COMPUTATIONAL STUDY OF CHEMICALLY REACTING HYPERSONIC FLOW Yoshiuru Funahashi Department o Aeronautics and Astronautics, Graduate School o Engineering, The University o Tokyo, Tokyo, JAPAN Keywords: Hypersonic
More informationStrain and Stress Measurements with a Two-Dimensional Detector
Copyright ISSN (C) 97-, JCPDS-International Advances in X-ray Centre Analysis, or Volume Diraction 4 Data 999 5 Strain and Stress Measurements with a Two-Dimensional Detector Baoping Bob He and Kingsley
More informationFeedback Linearization
Feedback Linearization Peter Al Hokayem and Eduardo Gallestey May 14, 2015 1 Introduction Consider a class o single-input-single-output (SISO) nonlinear systems o the orm ẋ = (x) + g(x)u (1) y = h(x) (2)
More informationAH 2700A. Attenuator Pair Ratio for C vs Frequency. Option-E 50 Hz-20 khz Ultra-precision Capacitance/Loss Bridge
0 E ttenuator Pair Ratio or vs requency NEEN-ERLN 700 Option-E 0-0 k Ultra-precision apacitance/loss ridge ttenuator Ratio Pair Uncertainty o in ppm or ll Usable Pairs o Taps 0 0 0. 0. 0. 07/08/0 E E E
More informationNON-SIMILAR SOLUTIONS FOR NATURAL CONVECTION FROM A MOVING VERTICAL PLATE WITH A CONVECTIVE THERMAL BOUNDARY CONDITION
NON-SIMILAR SOLUTIONS FOR NATURAL CONVECTION FROM A MOVING VERTICAL PLATE WITH A CONVECTIVE THERMAL BOUNDARY CONDITION by Asterios Pantokratoras School o Engineering, Democritus University o Thrace, 67100
More informationHYDROELASTIC TAILORING AND OPTIMIZATION OF A COMPOSITE MARINE PROPELLER
HYDROELASTIC TAILORING AND OPTIMIZATION OF A COMPOSITE MARINE PROPELLER José P. Blasques, Christian Berggreen and Poul Andersen Department o Mechanical Engineering, Technical University o Denmark Nils
More information