= Test your code using the following values: p = 2000 kg/m, To = 2500 K,...
Fantastic news! We've Found the answer you've been seeking!
Question:
Transcribed Image Text:
= Test your code using the following values: p = 2000 kg/m, To = 2500 K, k = 0.002 m/s/atm n = 0.5, y: 1.2, M = 10, R = 0.1 m, R = 0.2 m, L = 2 m, A* = R. Give a plot of stagnation pressure over time. It should show an increase, i.e., progressive thrust. (Al- though we are not plotting thrust, the thrust is proportional to the chamber pressure.) Also include two plots zooming in to show the start and stop transients. 0.5 write Matlab code to simulate the firing of a solid rocket motor with a single cylindrical grain with a circular perforation. In such a motor, the radius of the perforation (the hole in the propellant grain) increases, so that the wetted perimeter, the mass flow, the chamber pressure, and the thrust all increase with time. This is called a progressive burn. Your code should use the Matlab function ode 45 to integrate the equations of change for two state variables: the perforation radius r(t) and the chamber pressure Po(t). These equa- tions are dr dt where k is the burn rate constant and n is the pressure exponent, and dPo dt = mg = kpn, RTO V Here R is the gas constant, given by R = R/M, where R = 8314 J/(kmol K) is the uni- versal gas constant and M is the average molecular weight of the exhaust gas; To is the flame temperature (this will be given; it depends on the solid chemistry); V = rL is the chamber volume, with L the length of the motor; mg and m, are the rates of gas generation and gas leaving the nozzle, respectively. The gas generation rate is m = A* = -(g mn). dr dt -2rLp, where p is the propellant density. The rate of gas leaving the nozzle is (assuming ideal gas flow) 2 VRTO (1+1) where A* is the sonic (throat) area and y is the specific heat ratio. 2(7-1) Your function should start out like function [t, P0] = solidSim ( rho, TO, k, n, gamma, Mbar, Ri, Rf, L, Astar ) which will compute two arrays: the time and the chamber pressure. (Here Ri is the initial value of the perforation radius r, and Rf is the burnout value, which is the same as the inner radius of the motor wall.) Your code should set the burn rate to zero when the radius reaches the burnout value, but continue the integration after burnout until the chamber pressure de- creases to atmospheric pressure (1 atm = 1.01 10 Pa). = Test your code using the following values: p = 2000 kg/m, To = 2500 K, k = 0.002 m/s/atm n = 0.5, y: 1.2, M = 10, R = 0.1 m, R = 0.2 m, L = 2 m, A* = R. Give a plot of stagnation pressure over time. It should show an increase, i.e., progressive thrust. (Al- though we are not plotting thrust, the thrust is proportional to the chamber pressure.) Also include two plots zooming in to show the start and stop transients. 0.5 write Matlab code to simulate the firing of a solid rocket motor with a single cylindrical grain with a circular perforation. In such a motor, the radius of the perforation (the hole in the propellant grain) increases, so that the wetted perimeter, the mass flow, the chamber pressure, and the thrust all increase with time. This is called a progressive burn. Your code should use the Matlab function ode 45 to integrate the equations of change for two state variables: the perforation radius r(t) and the chamber pressure Po(t). These equa- tions are dr dt where k is the burn rate constant and n is the pressure exponent, and dPo dt = mg = kpn, RTO V Here R is the gas constant, given by R = R/M, where R = 8314 J/(kmol K) is the uni- versal gas constant and M is the average molecular weight of the exhaust gas; To is the flame temperature (this will be given; it depends on the solid chemistry); V = rL is the chamber volume, with L the length of the motor; mg and m, are the rates of gas generation and gas leaving the nozzle, respectively. The gas generation rate is m = A* = -(g mn). dr dt -2rLp, where p is the propellant density. The rate of gas leaving the nozzle is (assuming ideal gas flow) 2 VRTO (1+1) where A* is the sonic (throat) area and y is the specific heat ratio. 2(7-1) Your function should start out like function [t, P0] = solidSim ( rho, TO, k, n, gamma, Mbar, Ri, Rf, L, Astar ) which will compute two arrays: the time and the chamber pressure. (Here Ri is the initial value of the perforation radius r, and Rf is the burnout value, which is the same as the inner radius of the motor wall.) Your code should set the burn rate to zero when the radius reaches the burnout value, but continue the integration after burnout until the chamber pressure de- creases to atmospheric pressure (1 atm = 1.01 10 Pa).
Expert Answer:
Related Book For
Posted Date:
Students also viewed these mechanical engineering questions
-
QUIZ... Let D be a poset and let f : D D be a monotone function. (i) Give the definition of the least pre-fixed point, fix (f), of f. Show that fix (f) is a fixed point of f. [5 marks] (ii) Show that...
-
CANMNMM January of this year. (a) Each item will be held in a record. Describe all the data structures that must refer to these records to implement the required functionality. Describe all the...
-
You have recently been hired as a consultant for a personal financial planning firm. One of your first projects is creating a retirement plan for a couple, Tom and Helena Keeley. They have just...
-
Sound Audio manufactures and sells audio equipment for automobiles. Engineers notified management in December 2018 of a circuit flaw in an amplifier that poses a potential fire hazard. An intense...
-
An investor wishes to buy euros spot (at $1.3908) and sell euros forward for 180 days (at $1.3996). Is there any premium on 180-day Euro or is there any discount on 180-day forward? Show your workings
-
Water flows from a large tank as shown in Fig. P3.90. Atmospheric pressure is \(14.5 \mathrm{psia}\), and the vapor pressure is 1.60 psia. If viscous effects are neglected, at what height, \(h\),...
-
Virginia Partners, Ltd. (Virginia Partners), a limited partnership organized under the laws of Florida, conducted business in Kentucky but failed to register as a foreign limited partnership, as...
-
Write an essay that demonstrates your research on your technology topic and addresses the following: 1. Explanation of the technology and why it was developed. 2. Explain how the technology is being...
-
Bob's utility function is given by U(x,y) = x/2,1/2, His income is $100, the price of x is $20, and the price of y is $1. If the price of x decreases to $10, how much money would have to be given to...
-
What is wrong with the histogram in Figure 2.35? 0.3 0.2 0.1 0 0 1 2 3 4 5 Figure 2.35
-
Researchers used a random sample of size 100 to test the null hypothesis that the average difference between the household income of immigrant women and their adult daughters is equal to 0. Their...
-
What is wrong with the histogram in Figure 2.36? 0.4 0.3 0.2 0.1 0.0 -0.1 0 1 2 3 4 5 6 7 8 Figure 2.36
-
What is wrong with Figure 2.34, showing the distribution of family income? How would you define the middle class? Under $5,000 $5,000-$9,999 $10,000-$14,999 $15,000-$24,999 $25,000-$34,999...
-
To demonstrate the central limit theorem, draw 100 samples of size 5 from a random-number table and calculate the sample mean for each of the 100 samples. Construct a frequency distribution of sample...
-
A cylindrical garbage can of depth 3ft and radius 1 ft fills with rainwater up to a depth of 2ft. How much work would be done in pumping the water up to the top edge of the can? (water weighs 62.4...
-
Diamond Walker sells homemade knit scarves for $25 each at local craft shows. Her contribution margin ratio is 60%. Currently, the craft show entrance fees cost Diamond $1,500 per year. The craft...
-
a. List additional QTAs that could have been taken in addition to the ones listed for each CTA in the table. b. Apply the seven types of actions of CTA to Example 13-2 and write an action that should...
-
The production of ethylene glycol from ethylene chlorohydrin and sodium bicarbonate CH2OHCH2Cl + NaHCO3 (CH2OH)2 + NaCl + CO2 is carried out in a semibatch reactor. A 1.5-molar solution of ethylene...
-
The liquid-phase irreversible reaction A B + C is carried out in a CSTR. To learn the rate law, the volumetric flow rate, 0 , (hence = V/ 0 ) is varied and the effluent concentrations of species A...
-
Derive the differential equation governing the motion of the one degree-of-freedom system by applying the appropriate form(s) of Newton's laws to the appropriate free-body diagrams. Use the...
-
Determine the differential equations governing the motion of the system by using the equivalent systems method. Use the generalized coordinates shown in Figures P2.52. k FIGURE P 2.52 E x 2k www
-
Determine the differential equations governing the motion of the system by using the equivalent systems method. Use the generalized coordinates shown in Figures P2.53. 2r k ww T m FIGURE P 2.53
Study smarter with the SolutionInn App