Optimum Fan Pressure Ratio for Bypass Engines with Separate or Mixed Exhaust Streams

Size: px
Start display at page:

Download "Optimum Fan Pressure Ratio for Bypass Engines with Separate or Mixed Exhaust Streams"

Transcription

1 JOURNAL OF PROPULSION AND POWER Vol. 17, No. 5, September October 2001 Optimum Fn Pressure Rtio for Bypss Engines with Seprte or Mixed Exhust Strems Abhijit Guh University of Bristol, Bristol, Englnd BS8 1TR, United Kingdom The optimum fn pressure rtio is determined both numericlly nd nlyticlly for seprte-strem s well s mixed-strem bypss engines. The results re pplicble to civil nd militry engines. The optimum fn pressure rtio is shown to be predominntly function of the speci c thrust nd wek function of the bypss rtio. Two simple, explicit equtions, one for ech type of engine, hve been derived tht specify the optimum fn pressure rtio. The predictions compre very well ginst numericl optimiztion performed by specilist computer pckge employing itertive nd dvnced serch techniques nd rel gs properties. The nlyticl reltions ccelerte the optimiztion process nd offer physicl insight. It hs been shown tht the optimum fn pressure rtio chieves the conditionv jc /V jh = KE in seprte-strem engine nd the conditionv 163 /V 63 ¼ KE in mixedstrem engine. The condition V 163 /V 63 ¼ KE pplies even under situtions when signi cnt deprtures from the normlly ssumed condition p 016 /p 06 ¼ 1 cn occur. Nomenclture B = bypss rtio c p = speci c het of ir t constnt pressure c pg = speci c het of combustion products t constnt pressure E k = kinetic energy dded by the core engine F N = net thrust OF N = speci c thrust M = Mch number of the ircrft Pm c = mss ow rte through the bypss duct (cold ow) Pm f = mss ow rte of fuel Pm h = mss ow rte through the core engine (hot ow) p = sttic pressure p 0 = totl pressure R = speci c gs constnt of ir T = sttic temperture T = temperture of mbient ir T 0 = totl temperture V = velocity V = forwrd speed of the ircrft V j = jet speed in mixed strem engine V j c = fully expnded jet speed of the cold bypss strem V j h = fully expnded jet speed of the hot core strem V m = men jet speed [given by Eq. (8)] = isentropic index of ir g = isentropic index of combustion products f = isentropic ef ciency of the fn or low pressure compressor KE = ef ciency of energy trnsfer between the core nd bypss ow LPT = isentropic ef ciency of the low pressure turbine NB = isentropic ef ciency of the bypss nozzle Subscript op = optimum vlue Introduction IN this pper simple, explicit reltions hve been derived from fundmentl principles for determining the optimum fn (com- Received 30 My 2000; revision received 20 Februry 2001; ccepted for publiction 29 My Copyright c 2001 by Abhijit Guh. Published by the Americn Institute of Aeronutics nd Astronutics, Inc., with permission. Lecturer, Aerospce Engineering Deprtment, University Wlk; A.Guh@bristol.c.uk pressor) pressurertio in bypss engines. Seprte nlysishs been presented for the mixed-strem nd seprte-strem engines becuse the ow physics is different in the two cses. In both cses, the present theory gives results in good greement with numericl optimiztion clcultions. Determintion of the optimum fn pressure rtio (FPR op ) is importntbecuse, givenll other vriblessuch s the overll pressure rtio (OPR), bypss rtio B, nd turbine entry temperture T 04 the optimum vlue of FPR simultneously ensures mximum speci c thrust nd minimum speci c fuel consumption. Thus, the FPR op chieves the two, usully contrdictory, gols of jet engine design. For civil engines, prticulrlyfor medium to long rnge ircrfts, the FPR op cn be clculted under cruise conditions. For militry engines, optimiztionmy be mission speci c. However, Millhouse et l. 1 hve shown tht, lthough the most fuel-ef cient vlues of such prmeters s bypss rtio nd OPR depend on the integrted effects of the vrious ight ltitudes, Mch numbers, nd thrust required throughout the mission, it is possible to locte heuristiclly the most ef cient FPR to complement other engine prmeters independent of the mission. Thus, the reltions derived in this pper cn be used in the design of both civil nd militry engines. The design wisdom is tht only single-stge fn is used for civil ircrft engine with lrge bypss rtios, nd so the mximum FPR is bout 1.8. Engines with low bypss rtios, for exmple, militry engines, normlly use low-pressure (LP) compressor hving 3 5 stges. In this cse, the mximum possible FPR cn be much higher. The mixing of the two strems in mixed-strem turbofn engine offers performnce gin (lower speci c fuel consumption nd higher speci c thrust), even though mixing incurs loss in totl pressure becuse the enthlpiesof the two strems re redistributed. Moreover, only one rehet system nd only one vrible nozzle re needed. A long ncelle cn contin the fn/compressor noise more effectively.for the sme speci c thrust, the mixed-strem engine hs lower optimum FPR thn seprte-strem engine; thus, the LP turbine my hve one fewer stge if both engines hve the sme bypss rtio. The outer cold strem lso protects the jet pipe fromseveretemperturesresultingfrom fterburning.low-speci c- thrust, high-bypss-rtioengines, on the other hnd, do not employ fterburning,nd mixing of the two stremswould necessitte long ncelle incresing weight, cost, ncelle drg, nd interference effects. These engines, therefore, use seprte-stremcon gurtions. When common nozzle is, however, used for high-bypss-rtioengines,the smll xilextentof the mixingzone is usullyinsuf cient to chieve complete mixing of the two strems. Clcultion nd nlysis of the performnce of turbofn engines re discussed in detils in severl excellent texts. 2 4 In comprison,

2 1118 GUHA the issue of optimized design hs received less ttention, nd nlyticl solutions did not exist. The rel store house of expertise nd knowledge obviously exists with the eroengine mnufcturers who hve ctively pursued the optimiztion of gs turbine performnce for the pst 60 yers. Ruf es 5 shows tht there hs been 50% improvement in the speci c fuel consumption (SFC) of the bre engine, s the engines evolved from the turbojet of 1958 to the high-bypsstrent enginesof moderntime. Improvementsin propulsive ef ciency,componentef ciencynd cycle ef ciencycontribute pproximtely one-third ech to this gin in performnce. The publictions mde by experts from industry (e.g., Jckson, 6 Bennett, 7 Wilde, 8 nd Birch 9 ) contin welth of experience, informtion, nd clcultion results. However, the detils of clcultion method re not known to the reders. The industry would typiclly hve its own, sophistictedcomputer pckges to which the public does not hve ccess. Moreover, the publishedresults often would involve n dopted fmily of engine designs. It is then not lwys obvious to the generl reder how to generlize the results or how to proceed on clen-sheetnlysis. In this pper nlyticlreltionshve been derivedfor clculting FPR op. The nlyticl results not only ccelerte the optimiztion process but lso offer vluble physicl insight. It is shown tht the FPR op is predominntly function of the speci c thrust nd only wekly depends on the bypss rtio. The FPR op monotoniclly increses with incresing speci c thrust t two differentlevels corresponding to the seprte-strem nd mixed-strem turbofn engines. The theory hs been veri ed by compring its predictions with clcultionsperformed by the computer progrm GsTurb, 10 which is generl purpose commercil pckge for the clcultion of design nd off-design performnce of vrious types of gs turbine engines, turbojet, turboshft, seprte- or mixed-strem turbofn, nd geredturbofn,ll with single-or twin-spoolcon gurtions.it uses rel gs propertieswith dissocitionnd hs severl fetures including optimiztion, itertion, trnsient, nd Monte Crlo simultion. Ech numericl point subsequently illustrted hs been clculted by GsTurb by lengthy optimiztion process tht uses dvnced serch techniques to nd the optimum FPR tht minimizes the SFC while, t ech tril, the primry vribles re iterted to give the prescribed speci c thrust. This numericl optimiztion is used not only to verify the theoreticl prediction, but lso to discover wht exctly hppens to the mny vribles when the FPR op is chieved. Bsed on the principles discussed in the present pper, new methodologyfor the thermodynmicoptimiztionof bypssengines (turbofn or dvncedpropulsors) hs been developed by Guh 11 in which the optimum combintion of ll vribles is determined concurrently. The process strts with estblishing n optimum speci c thrust for the engine bsed on n economic nlysis (instlltion constrints, noise regultions,etc., lso need to be considered). The tsk of the optimiztion process is then to nd the combintion of optimum vlues of OPR, B, FPR, nd T 04 concurrently tht minimizes SFC t the xed speci c thrust.this procedureis quite different from the usul prmetric studies of engine performnce, 2;12 in which single prmeteris vried ech time, while keeping ll other prmeters xed, therefore t their nonoptimum vlues. Moreover the usul single-vribleprmetric studies my involve lrge excursion in the vlue of spe cic thrust uncceptble for speci c mission. Guh 11 hs lso discussed t length the determintion of optiml jet velocity, with the derivtion of new nlyticl expression tht performs well ginst numericl optimiztion results. Following industril prctice, the SFC nd the speci c thrust, re expressed here respectively in pounds mss per hour per pounds force nd pounds force per pounds mss per second where 1 lbf/lbm/s D 9:81 m/s nd 1 lbm/hr/lbf D 28: kg/s/n. Seprte-Strem Bypss Engines The schemtic rrngement nd ow structure in seprtestrem bypss engine is shown in Fig. 1. The description of n equivlent turbojet engine is lso shown in Fig. 1 s reference point in the nlysis. The core gs genertor produces, for pr- Turbojet Seprte-strem bypss engine Fig. 1 Schemtic description of seprte-strem turbofn engine nd n equivlent turbojet engine used for nlysis; energy is the sme t point in both engines. ticulr vlue of fuel consumption rte, jet with speed V j;ref in the bsence of ny bypss ow. The kinetic energy dded by the gs genertor is then, neglecting the fuel mss ow rte which is usully very smll proportion of the ir mss ow rte in jet engine, equl to E k D 1 Pm 2 h V 2 j;ref (1) In the bypss engine, the whole of this kinetic energy is, of course, not present in the core jet. The LP turbine extrcts portion of this energy nd turns the fn, which, in turn, dds kinetic energy to the bypss ow. The ef ciency KE of this energy trnsfer between the core nd bypss ow depends on the component ef ciencies of the LP turbine, the fn, nd the bypss nozzle nd lso on the ow lossesin the ductsndmechniclef ciencyof the shfts.assuming the mechnicl ef ciency is close to 1, KE ¼ LPT f NB (2) At the current level of technologies, 4;9 LPT ¼ 0:9, f ¼ 0:9, nd NB ¼ 1. Eqution (2), therefore, suggests tht KE ¼ 0:8. In bypss engine, the kinetic energy given by Eq. (1) is relted to the sum of the kinetic energies of the two strems by E k D 1 Pm 2 h V 2 j;ref D 1 Pm 2 h C 1 2 Pm h.b= KE / j c V 2 j h In Eq. (3) we hve expressedthe mss ow rte of the bypssstrem, Pm c, in terms of Pm h nd B. Similrly the net thrust F N produced by the bypss engine is the sum of the thrusts of the two strems nd is given by (3) F N D Pm h.v j h V / C B Pm h.v jc V / (4) The pressure thrust hs been tken into ccount becuse V jh nd V j c re the fully expnded jet velocities of the two strems. The objectiveof the optimiztionprocess is to determinethe rtio V jc =V j h, which would mximize the net thrust F N while keeping the fuel consumption xed. The condition for mximum net thrust is given F jh D 0 (5) At constnt therml ef ciency, xed fuel consumption mens E k in Eq. (3) is constnt. k jh D 0 (6)

3 GUHA 1119 For constnt bypss rtio B, Eqs. (3 6) cn be combined to give.v j c =V j h / op D KE (7) While evlutingthe prtilderivtivesin Eqs. (5) nd (6), B nd Pm h re ssumed constnt.togetherthey specify constnt inlet ir ow. Becuse the speci c thrustis the rtio of net thrustnd mss ow rte of ir t the inlet, Eq. (5) is lso the condition for mximum speci c thrust. SFC is the rtio of mss ow rte of fuel nd net thrust. Eqution (6) implies xed mss ow rte of fuel; Eq. (5) is, thus, lso the condition for minimum SFC. Becuse Eq. (7) gives both mximum speci c thrust nd minimum SFC, this condition would be used in the nlyticl formultion of the optimiztion process tht follows. A designer cn stisfy Eq. (7) in rel engine by choosing the FPR op, s will be shown. The men jet speed V m cn be expressed, using Eq. (7), s Fig. 3 Comprison of present theory with numericl optimiztion results for FPR op in seprte-strem bypss engines. (Altitude = 11 km, M = 0:82, isentropic ef ciency of compressors nd turbines = 0:9). V m D [1=.1 C B/].B C 1= KE /V jc (8) nd it follows tht OF N D F N =. Pm h C Pm c / D V m V (9) Suppose the rise in totl temperture in the bypss ow for hypotheticl, isentropic compression by the fn, for the sme pressure rtio tht exists cross rel fn, is 1T 0;isen. If the fn nd the bypss nozzle were isentropic, then the sttic temperture of the fully expnded jet would be equl to the mbient temperture. Appliction of the stedy ow energy eqution then gives 1T 0;isen D.1=2c p / jc (10) Equtions (8 10) cn be combined with the isentropic pressure temperture reltion. One then obtins, fter some lgebric mnipultion, the nl expression for the FPR op :. 1/.FPR/ op D 1 C 2 C. 1/M 2 ".1 C B/ 2 OF p N C M.B C 1= KE / 2 RT ¼ 2 M 2 # (11) All quntities in Eq. (11) re nondimensionl numbers. Figure 2 shows grphiclly the prediction of Eq. (11). One interesting feture of Eq. (11) is tht if KE D 1, tht is, in the cse of idel energy trnsfer from the core strem to the bypss strem, the bypss rtio B drops out of the eqution. FPR op would then hve depended only on OF N. At ny rte, Eq. (11) predicts tht, for xed OF N, the dependenceof FPR op on B is wek, especilly t highervluesof B (B > 4). Eqution(11) lso producesthe expected limiting result: lim B! 1 OF N! 0.FPR/ op D 1 For chosen speci c thrust OF N nd bypss rtio B, Eq. (11) immeditely speci es the FPR op. The llowble rnge of FPR is rther restricted. The design wisdom is tht only single-stge fn is used for civil ircrft with lrge-bypss rtios, the mximum Fig. 2 Vrition of FPR op with speci c thrust in seprte-strem bypss engines: prediction of Eq. (11); M = 0:82, T = 216:65 K, KE = 0:81, nd = 1:4. Fig. 4 Dt for vrious current civil turbofn engines under cruise conditions (typiclly 35,000 ft, M = ). FPR is, therefore, restricted to 1.8. A limit of minimum FPR my lso rise due to fn instbilities t prtlod, off-design conditions. According to Rüd nd Lichtfuss, 13 below n FPR of 1.4, the engine will require vrible geometry either vi vrible pitch fn or vrible re bypss nozzle to provide erodynmic fn stbility under prtlod conditions. Therefore, knowing FPR op t n erly stge of the design is n dvntge. If its vlue does not lie within the desired limit, the choice of OF N nd B cn be ltered (of course, nother option would be not to dhere strictly to the optimum vlue of FPR). Eqution (11) hs been tested ginst optimiztion performed by the computer progrm GsTurb. 10 Figure 3 shows the comprison. Ech numericl point hs been clculted by GsTurb by lengthy optimiztion process tht uses dvnced serch techniques to nd the optimum FPR tht minimizes the SFC while, t ech tril, the primry vribles re iterted to give the prescribed speci c thrust. Figure 3 shows tht Eq. (11) performs well ginst sophisticted numericl optimiztion nd cn, therefore, be used for the design of rel engines. To pprecite the current design stndrd, Fig. 4 plots the SFC nd speci c thrust of some current civil turbofn engines of vrious mnufcturers. 14 Speci c thrust dt in Fig. 4 re pproximte becuse lthough net thrust is known t cruise, the mss ow rte is estimted from given vlues t se level sttic conditions by dynmic scling, tht is, ssuming the sme nondimensionloperting point. The speci c thrust dt indictesover which rnge ny theory needs to be pplicble in order to be relevnt for current nd future designs. Figure 4 shows tht for most existing engines, the cruise speci c thrust lies in the bnd lbf/lbm/s. Over the pst 40 yers of civil engine design, the speci c thrust hs reduced signi cntly (producingpprecibleimprovementin propulsiveef ciency) while the bypss rtio hs incresed from 1 2 in the 1960s to 7 9 in the 1990s. 9 The forecst, 9;13 is tht the design driver for future engines would be towrd even lower speci c thrust. A surge in fuel price nd/or the introduction of more stringent noise regultion my necessitte such designs. Future gered turbofn or dvnced ducted propulsors would reduce the speci c thrust signi cntly. Figure 3 thus shows tht Eq. (11) is pplicble to both current nd future bypss engines. Mixed-Strem Turbofn Engines Before developing the theory for FPR op for mixed-strem turbofn engines, the computer pckge GsTurb ws used to determine

4 1120 GUHA the sme numericlly. Figure 5 shows the results of this numericl optimiztion. Clcultions re shown t two bypss rtios. Mixedstrem engines used for militry purposes employ much higher speci c thrust thn tht used for seprte-strem civil pplictions. These engines usully hve smller bypss rtio. Therefore, mximum tolerble turbine entry tempertures cn produce lrge speci c thrust. Of course, such high vlues of speci c thrust would men low propulsive ef ciency giving higher fuel consumption. The reduction in engine size nd weight is, however, more crucil for such pplictions. Figure 6 shows the stndrd sttion numbering for mixed- ow turbofn engine. In the following nlysis these numbers re used s subscripts to denote ow vribles t vrious loctions. The following representtive vlues give resonble pproximtion to the propertiesof ir nd combustionproducts: R D 287 J/kg/K, D 1:4, c p D R=. 1/, g D 1:33, c pg D g R=. g 1/. A simple, nlyticl procedurefor clcultingthe propertiesof combustion products s function of temperture, fuel ir rtio, nd fuel composition is given by Guh. 15 The nlysis cn be formulted more esily if the optimum FPR ndthe speci c thrustre expressedin terms of the primry vribles OPR, B, nd T 04. An energy blnce gives (with T 05 D T 06 nd T 013 D T 016 ).c pg =c p /.T 04 T 06 / D.T 03 T 02 / C B.T 016 T 02 / (12) where the vrious totl tempertures cn be clculted from T 02 =T D 1 C 0:5. 1/M 2 (13) T 03 T 02 D.1= c/t 02 OPR 1 (14) T 016 T 02 D.1= f /T 02 FPR 1 (15) T 04 T 06 D t T 04 1.OPR/FPR/.1 g /= g (16) where c is the isentropic ef ciency of the compression between points 2 nd 3 nd t is the isentropic ef ciency of the expnsion between points 4 nd 6. While writing Eq. (16) it is ssumed tht p 016 =p 06 D 1 when optimum fn pressure rtio is chieved. Strictly speking, this condition is not lwys fully stis ed. However, the vlue of the rtio p 016 =p 06 rises under off-designconditions:if the vlue of this rtio used t the design point is greter thn one, serious mixing losses my result t off-designopertions.therefore, it is prudent to choose the vlue of p 016 =p 06 between 0.95 nd 1 t the design point even if it is slightly suboptimlt tht prticulropertingpoint. Numericlclcultions of Millhouse et l. 1 show tht if the condition p 016 =p 06 ¼ 1 is used t the design point, the computed FPR then gives the optimum vlue independent of speci c mission of the ircrft. Equtions (12 16) cn be solved to determine the optimum FPR. This will, however, require n itertive solution. For smll-bypss rtios, s is usul for mixed-strem militry engines, the second term in the right-hnd side of Eq. (12) is smll compred to the rst term. The index for FPR in Eq. (15) cn then be pproximted by substituting g for. With this pproximtion,the equtionscn be combined to give n explicit reltion for the FPR op : µ FPR. g 1/= g BT02 op C.c pg=c p / t T 04 c pg D f OPR. g 1/= g t T 04 c p 1 T 02 OPR c 1 C BT 02 f (17) If V j is the jet speed nd V is the ircrft speed, the speci c thrust is clculted from where OF N D V j V (18) V j D p 2c pg.t 064 T j / (19) T 064 D T 06 C BT C B (20) T j D T 064.FPR op f M /.1 g/= g (21) f M D [1 C 0:5. 1/M 2 ] =. 1/ (22) Eqution (20) gives the totl temperture of the mixed strem nd Eq. (21) gives the sttic temperture of the jet. While writing the precedingequtions, the losses in the fn, compressorsnd turbines hve been ccounted for, but other losses, for exmple, tht in the vrious ducts, hve been neglected. Figure 7 shows the comprison of the present theory nd the numericl optimiztion results using GsTurb. Ech numericl point hs been clculted by GsTurb to nd the optimum FPR tht minimizes the SFC, while, t ech tril, the primry vriblesre iterted to give the prescribedspeci c thrust. The nozzle re ws lso optimized t ech point so tht it produced fully expndedjet. The convergenceof the numericloptimiztionmy be slow nd sometimes, dependingon the strtingpoint, the numericlserch techniquemy not nd the optimum or my nd locl rther thn the globl optimum. The nlyticlequtions, on the other hnd, redily provide the nswer. Figure 7 shows tht the present theory gives ccurte results: The itertive nlyticl solution [Eqs. (12 16)] grees lmost Fig. 5 FPR op for mixedstrem turbofn engines determined by numericl optimiztion (GsTurb); for ll clcultions, OPR = 17:5, M = 0:82, isentropic component ef ciencies = 0:9 nd mix = 1. Fig. 7 Comprison of the present theory with numericl optimiztion results for FPR op in mixed-strem turbofn engines; for ll clcultions, OPR = 17:5, B = 0:5, M = 0:82, isentropic component ef ciencies = 0:9 nd mix = 1. Fig. 6 Stndrd nomenclture for vrious loctions in mixed-strem turbofn engine (only hlf of the engine is shown).

5 GUHA 1121 Tble 1 Dependence of FPR op on OPR for mixed-strem engines Tble 2 Comprison of theory nd numericl results for optimum V jc /V jh OF N, lbf/lbm/s OPR T 04, K B op FPR op identiclly with GsTurb optimiztion; the explicit eqution (17) lso gives quite cceptble nswer until the speci c thrust becomes very high. The present theory cn be extended to include the effects of incompletemixing by introducing prmeter mix. It is then ssumed tht three seprte strems hving totl tempertures T 016, T 06, nd T 064 expnd through the sme pressure rtio. The totl thrust is the sum of the thrust produced by these individul strems: p.1 /= V j16 D 2c p.t 016 T j16 /; T j16 D T 016.FPR op f M / p.1 g /= g V j6 D 2c pg.t 06 T j6 /; T j 6 D T 06.FPR op f M / OF N D mix V j C.1 mix /.V j6 C BV j16 /=.1 C B/ V (23) Seprte Versus Mixed Strem Engines Figure 5 cn be compred with Figs. 2 nd 3. Three observtions cn be mde. 1) The FPR op for mixed-stremengines rises monotonicllywith the speci c thrust, s it does for seprte strem engines. 2) The FPR op is predominntly function of the speci c thrust nd only wekly depends on the bypss rtio. The dependence on bypss rtio is weker in mixed-stremenginesthn tht in seprtestrems engines. 3) At prticulr vlue of speci c thrust, the FPR op for mixedstrem engines is lower thn tht in the seprte-strems engines. In seprte-strem engine, ltering the vlues of OPR does not lter the vlue of FPR op, so long s the speci c thrust nd bypss rtio re kept xed. This is predicted by Eq. (11) nd borne out by numericl optimiztion. In mixed-strem engine, even t xed speci c thrust, vrious vlues of OPR nd M chnge the vlue of FPR op. Either the nlyticl reltions (12 16) or the numericl optimiztions show this. Tble 1 shows the results of optimiztion performed by GsTurb (both B nd FPR re optimized t xed speci c thrust t two levels of OPR, minimising SFC). Effects of FPR op on the Conditions of Other Vribles The behviorof other ow vribles re now exmined when the FPR op is chieved. Seprte discussion is provided for the mixedstrem nd seprte-stremengines becuse differentphysicl principles re involved. Seprte-Strem Bypss Engines The conditions to which the FPR op correspond hve been determined nlyticlly in this cse. Eqution (7) shows tht there is prticulr rtio of the fully expnded jet velocities of the cold nd hot strem [.V j c =V j h / op D KE ] tht simultneously chieves minimum SFC nd mximum speci c thrust. The designer cn mke the two jet velocities tke this vlue by choosing the FPR op. The computtionl progrm GsTurb ws used to nd the rtio V j c =V j h tht occurs in n engine whose FPR hs been numericlly optimized for minimizing SFC. Tble 2 shows these vlues (for OPR D 30, T 04 D 1200 K, LPT D 0:9, f D 0:9, nd NB D 1), ginst the prediction of Eq. (7). Numericl optimiztions were performed t severl other OPR nd T 04 : V j c =V j h ly between0.77nd 0.82 (for the ssumedcomponent ef ciencies). The numericl clcultionsshow tht the derived nlyticl reltion [Eq. (7)] is pproximtely vlid. Mixed-Strem Bypss Engines The rtio of the totl pressure of the two strems t the beginning of the mixer zone, p 016 =p 06, is n importnt vrible nd the.v j c =V jh / op Bypss Numericl optimiztion.v jc =V j h / op D KE rtio with GsTurb [ KE ¼ LPT f NB ] Tble 3 Vlidity of Eq. (24) for vrious component ef ciencies f LPT. p 016 =p 06 / op.v 163 =V 63 / op KE D f LPT Clculted by numericl optimiztion with GsTurb. Fig. 8 Effects of design mixer Mch number on the two chrcteristic rtios when FPR op is chieved; for ll clcultions, M = 0:82, B = 0:5, ÃF N = 60 lb/lb/s, nd component isentropic ef ciencies = 0:9. FPR controls this. It is generlly recognized, 1;10 tht when FPR op is chieved, the rtio p 016 =p 06 is close to one. The physicl reson is tht, when the totl pressures of the two strems re nerly equl, the loss (entropy genertion) due to mixing would be smll. The computer progrm GsTurb ws used to determine the vlue of FPR (by rndom dptive serch) tht minimizes the SFC t severl speci c thrust levels tht were xed by iterting OPR nd T 04. It ws found tht optimum p 016 =p 06 did not remin xed but vried. Optimum p 016 =p 06 incresed with incresing speci c thrust nd bypss rtios. The optimum vlue of p 016 =p 06 ws found to vry prticulrly strongly with design mixer Mch number M 64, s shown in Fig. 8. However, the present study discovered new reltion between the velocities of the two strems t the mixer inlet under optimum conditions:.v 163 =V 63 / op ¼ KE (24) where KE D f LPT. At the FPR op, the reltion V 163 =V 63 ¼ KE remined nerly vlid t ll vlues of speci c thrust, bypss rtio, nd design mixer Mch number. The similrity of this condition with Eq. (7) derivedfor seprte-strembypssenginesis noticeble.the vrition in the velocity rtio is shown in Fig. 8 (with KE D 0:81). To ssess the ccurcy of Eq. (24), lrge number of numericl experiments were conducted with GsTurb in which the isentropic ef ciencies of vrious components were chnged systemticlly. It ws found tht indeed the ef ciencies of the high-pressure (HP) compressor, HPC, nd the HP turbine, HPT, did not in uence the vlue of.v 163 =V 63 / op, s Eq. (24) suggests. Eqution (24) ws lso stis ed when vrious vlues of the LP compressor nd LP turbine ef ciencies were tested. A few of these numericl optimiztion results re given in Tble 3 to demonstrte this point. (For ll clcultions shown in Tble 3, OF N D 60 1b/1b/s, B D 1, OPR D 17:5, M 64 D 0:3, M D 0:82, nd mix D 1. Similr results re obtined t other vlues of these vribles.) Conclusions The optimum fn (compressor) pressure rtio is determined both numericlly nd nlyticlly for seprte-strem s well s mixedstrem bypss engines. The FPR op is shown to be predominntly function of the speci c thrust nd wek function of the bypss

6 1122 GUHA rtio.for mixed-stremengines,the dependenceof FPR op on bypss rtio is very wek, nd FPR op lso depends on OPR. At the sme speci c thrustndbypssrtio,the FPR op for mixed-stremengine is lower thn tht of seprte-stremengine. Two simple, explicitequtionshvebeenderivedfrom fundmentl principles.eqution(11) gives the FPR op for seprte-stremengines,wheres Eq. (17) [or the equtionset (12 16)] givesthe FPR op for mixed-strem engines. The ccurcy of the nlyticl formuls hs been estblished through extensive veri ction by numericl optimiztion results of the commercil computer pckge GsTurb (Figs. 3 nd 7). The nlyticl results ccelerte the optimiztion process nd offer physicl insight. It hs been shown tht the FPR op chieves the condition V j c =V j h D KE in seprte-strem engine nd the condition V 163 =V 63 ¼ KE in mixed-strem engine. The condition V 163 =V 63 ¼ KE pplies even under situtions when signi - cnt deprtures from the normlly ssumed condition p 016 =p 06 ¼ 1 occur. Acknowledgment The uthor is grteful to the referees for their helpful comments. References 1 Millhouse,P. T., Krmer, S. C., King, P. I., nd Mykytk,E. F., Identifying Optiml Fn Compressor Pressure Rtios for themixed-strem Turbofn Engine, Journl of Propulsion nd Power, Vol. 16, No. 1, 2000, pp Hill, P. G., nd Peterson, C. R., Mechnics nd Thermodynmics of Propulsion, 2nd ed., Addison Wesley Longmn, Reding, MA, 1992, Chp Cohen, H., Roger, G. F. C., nd Srvnmuttoo, H. I. H., Gs Turbine Theory, 4th ed., Longmn, Hrlow, Englnd, U.K., 1996, Chp Cumpsty, N. A., Jet Propulsion, Cmbridge Univ. Press, Cmbridge, Englnd, U.K., Ruf es, P. C., The Future of Aircrft Propulsion, Journl of Mechnicl Engineering Science, Vol. 214, No. C1, 2000, pp Jckson, A. J. B., Some Future Trends in Aero Engine Design for Subsonic Trnsport, Journl of Engineering for Power, Vol. 98, April 1976, pp Bennett, H. W., Aero Engine Development for the Future, Proceedings of the Inst. of Mechnicl Engineers, Series A, Vol. 197, July 1983, pp Wilde, G. L., Future Lrge Civil Turbofns nd Powerplnts, Aeronuticl Journl, Vol. 82, July 1978, pp Birch, N. T., 2020Vision: the Prospects for Lrge Civil Aircrft Propulsion, Aeronuticl Journl, Vol. 104, No. 1038, 2000, pp Kurzke, J., Mnul for GsTurb 7.0 for Windows, MTU, München, Germny, Guh, A., Optimiztion of Aero Gs Turbine Engines, Aeronuticl Journl, Vol. 105, No. 1049, 2001, pp Mttingly, J. D., Elements of Gs Turbine Propulsion, McGrw Hill, New York, 1996, pp Rüd, K., nd Lichtfuss, H. J., Trends in Aero-Engines Development, Aspects of Engine Airfrme Integrtion for Trnsport Aircrft, Proceedings of DLR Workshop, DLR-Mitteilung, Brunschweig, Germny, Mrch Rolls Royce Aero Dt, TS 1491, Issue 9, Derby, Englnd, U.K., Guh, A., An Ef cient Generic Method for Clculting the Properties of Combustion Products, Journl of Power nd Energy, Vol. 215, No. A3, 2001, pp

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations ME 3600 Control ystems Chrcteristics of Open-Loop nd Closed-Loop ystems Importnt Control ystem Chrcteristics o ensitivity of system response to prmetric vritions cn be reduced o rnsient nd stedy-stte responses

More information

New Expansion and Infinite Series

New Expansion and Infinite Series Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University

More information

New data structures to reduce data size and search time

New data structures to reduce data size and search time New dt structures to reduce dt size nd serch time Tsuneo Kuwbr Deprtment of Informtion Sciences, Fculty of Science, Kngw University, Hirtsuk-shi, Jpn FIT2018 1D-1, No2, pp1-4 Copyright (c)2018 by The Institute

More information

Entropy ISSN

Entropy ISSN Entropy 006, 8[], 50-6 50 Entropy ISSN 099-4300 www.mdpi.org/entropy/ ENTROPY GENERATION IN PRESSURE GRADIENT ASSISTED COUETTE FLOW WITH DIFFERENT THERMAL BOUNDARY CONDITIONS Abdul Aziz Deprtment of Mechnicl

More information

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn

More information

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies Stte spce systems nlysis (continued) Stbility A. Definitions A system is sid to be Asymptoticlly Stble (AS) when it stisfies ut () = 0, t > 0 lim xt () 0. t A system is AS if nd only if the impulse response

More information

In the next lecture... Tutorial on ideal cycles and component performance.

In the next lecture... Tutorial on ideal cycles and component performance. In the next lecture... utoril on idel cycles nd component performnce. rof. Bhskr Roy, rof. A M rdeep, Deprtment of Aerospce, II Bomby roblem # Lect-9 A Bryton cycle opertes with regenertor of 75% effectiveness.

More information

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by. NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with

More information

Tests for the Ratio of Two Poisson Rates

Tests for the Ratio of Two Poisson Rates Chpter 437 Tests for the Rtio of Two Poisson Rtes Introduction The Poisson probbility lw gives the probbility distribution of the number of events occurring in specified intervl of time or spce. The Poisson

More information

Review of Calculus, cont d

Review of Calculus, cont d Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some

More information

Synoptic Meteorology I: Finite Differences September Partial Derivatives (or, Why Do We Care About Finite Differences?

Synoptic Meteorology I: Finite Differences September Partial Derivatives (or, Why Do We Care About Finite Differences? Synoptic Meteorology I: Finite Differences 16-18 September 2014 Prtil Derivtives (or, Why Do We Cre About Finite Differences?) With the exception of the idel gs lw, the equtions tht govern the evolution

More information

CBE 291b - Computation And Optimization For Engineers

CBE 291b - Computation And Optimization For Engineers The University of Western Ontrio Fculty of Engineering Science Deprtment of Chemicl nd Biochemicl Engineering CBE 9b - Computtion And Optimiztion For Engineers Mtlb Project Introduction Prof. A. Jutn Jn

More information

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus Unit #9 : Definite Integrl Properties; Fundmentl Theorem of Clculus Gols: Identify properties of definite integrls Define odd nd even functions, nd reltionship to integrl vlues Introduce the Fundmentl

More information

CHAPTER 4a. ROOTS OF EQUATIONS

CHAPTER 4a. ROOTS OF EQUATIONS CHAPTER 4. ROOTS OF EQUATIONS A. J. Clrk School o Engineering Deprtment o Civil nd Environmentl Engineering by Dr. Ibrhim A. Asskk Spring 00 ENCE 03 - Computtion Methods in Civil Engineering II Deprtment

More information

Numerical integration

Numerical integration 2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter

More information

13: Diffusion in 2 Energy Groups

13: Diffusion in 2 Energy Groups 3: Diffusion in Energy Groups B. Rouben McMster University Course EP 4D3/6D3 Nucler Rector Anlysis (Rector Physics) 5 Sept.-Dec. 5 September Contents We study the diffusion eqution in two energy groups

More information

Math 8 Winter 2015 Applications of Integration

Math 8 Winter 2015 Applications of Integration Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl

More information

Ordinary differential equations

Ordinary differential equations Ordinry differentil equtions Introduction to Synthetic Biology E Nvrro A Montgud P Fernndez de Cordob JF Urchueguí Overview Introduction-Modelling Bsic concepts to understnd n ODE. Description nd properties

More information

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS. THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem

More information

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite Unit #8 : The Integrl Gols: Determine how to clculte the re described by function. Define the definite integrl. Eplore the reltionship between the definite integrl nd re. Eplore wys to estimte the definite

More information

Duality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below.

Duality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below. Dulity #. Second itertion for HW problem Recll our LP emple problem we hve been working on, in equlity form, is given below.,,,, 8 m F which, when written in slightly different form, is 8 F Recll tht we

More information

CHM Physical Chemistry I Chapter 1 - Supplementary Material

CHM Physical Chemistry I Chapter 1 - Supplementary Material CHM 3410 - Physicl Chemistry I Chpter 1 - Supplementry Mteril For review of some bsic concepts in mth, see Atkins "Mthemticl Bckground 1 (pp 59-6), nd "Mthemticl Bckground " (pp 109-111). 1. Derivtion

More information

Math 1B, lecture 4: Error bounds for numerical methods

Math 1B, lecture 4: Error bounds for numerical methods Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the

More information

INTRODUCTION. The three general approaches to the solution of kinetics problems are:

INTRODUCTION. The three general approaches to the solution of kinetics problems are: INTRODUCTION According to Newton s lw, prticle will ccelerte when it is subjected to unblnced forces. Kinetics is the study of the reltions between unblnced forces nd the resulting chnges in motion. The

More information

Review of Gaussian Quadrature method

Review of Gaussian Quadrature method Review of Gussin Qudrture method Nsser M. Asi Spring 006 compiled on Sundy Decemer 1, 017 t 09:1 PM 1 The prolem To find numericl vlue for the integrl of rel vlued function of rel vrile over specific rnge

More information

An approximation to the arithmetic-geometric mean. G.J.O. Jameson, Math. Gazette 98 (2014), 85 95

An approximation to the arithmetic-geometric mean. G.J.O. Jameson, Math. Gazette 98 (2014), 85 95 An pproximtion to the rithmetic-geometric men G.J.O. Jmeson, Mth. Gzette 98 (4), 85 95 Given positive numbers > b, consider the itertion given by =, b = b nd n+ = ( n + b n ), b n+ = ( n b n ) /. At ech

More information

The steps of the hypothesis test

The steps of the hypothesis test ttisticl Methods I (EXT 7005) Pge 78 Mosquito species Time of dy A B C Mid morning 0.0088 5.4900 5.5000 Mid Afternoon.3400 0.0300 0.8700 Dusk 0.600 5.400 3.000 The Chi squre test sttistic is the sum of

More information

The Moving Center of Mass of a Leaking Bob

The Moving Center of Mass of a Leaking Bob The Moving Center of Mss of Leking Bob rxiv:1002.956v1 [physics.pop-ph] 21 Feb 2010 P. Arun Deprtment of Electronics, S.G.T.B. Khls College University of Delhi, Delhi 110 007, Indi. Februry 2, 2010 Abstrct

More information

SUMMER KNOWHOW STUDY AND LEARNING CENTRE

SUMMER KNOWHOW STUDY AND LEARNING CENTRE SUMMER KNOWHOW STUDY AND LEARNING CENTRE Indices & Logrithms 2 Contents Indices.2 Frctionl Indices.4 Logrithms 6 Exponentil equtions. Simplifying Surds 13 Opertions on Surds..16 Scientific Nottion..18

More information

Part I: Basic Concepts of Thermodynamics

Part I: Basic Concepts of Thermodynamics Prt I: Bsic Concepts o Thermodynmics Lecture 4: Kinetic Theory o Gses Kinetic Theory or rel gses 4-1 Kinetic Theory or rel gses Recll tht or rel gses: (i The volume occupied by the molecules under ordinry

More information

Week 10: Line Integrals

Week 10: Line Integrals Week 10: Line Integrls Introduction In this finl week we return to prmetrised curves nd consider integrtion long such curves. We lredy sw this in Week 2 when we integrted long curve to find its length.

More information

4.4 Areas, Integrals and Antiderivatives

4.4 Areas, Integrals and Antiderivatives . res, integrls nd ntiderivtives 333. Ares, Integrls nd Antiderivtives This section explores properties of functions defined s res nd exmines some connections mong res, integrls nd ntiderivtives. In order

More information

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1 3 9. SEQUENCES AND SERIES 63. Representtion of functions s power series Consider power series x 2 + x 4 x 6 + x 8 + = ( ) n x 2n It is geometric series with q = x 2 nd therefore it converges for ll q =

More information

The First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).

The First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a). The Fundmentl Theorems of Clculus Mth 4, Section 0, Spring 009 We now know enough bout definite integrls to give precise formultions of the Fundmentl Theorems of Clculus. We will lso look t some bsic emples

More information

Improper Integrals, and Differential Equations

Improper Integrals, and Differential Equations Improper Integrls, nd Differentil Equtions October 22, 204 5.3 Improper Integrls Previously, we discussed how integrls correspond to res. More specificlly, we sid tht for function f(x), the region creted

More information

Psychrometric Applications

Psychrometric Applications Psychrometric Applictions The reminder of this presenttion centers on systems involving moist ir. A condensed wter phse my lso be present in such systems. The term moist irrefers to mixture of dry ir nd

More information

4 The dynamical FRW universe

4 The dynamical FRW universe 4 The dynmicl FRW universe 4.1 The Einstein equtions Einstein s equtions G µν = T µν (7) relte the expnsion rte (t) to energy distribution in the universe. On the left hnd side is the Einstein tensor which

More information

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner

More information

Calculus of Variations

Calculus of Variations Clculus of Vritions Com S 477/577 Notes) Yn-Bin Ji Dec 4, 2017 1 Introduction A functionl ssigns rel number to ech function or curve) in some clss. One might sy tht functionl is function of nother function

More information

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION XX IMEKO World Congress Metrology for Green Growth September 9,, Busn, Republic of Kore THERMAL EXPANSION COEFFICIENT OF WATER FOR OLUMETRIC CALIBRATION Nieves Medin Hed of Mss Division, CEM, Spin, mnmedin@mityc.es

More information

Numerical Integration

Numerical Integration Chpter 5 Numericl Integrtion Numericl integrtion is the study of how the numericl vlue of n integrl cn be found. Methods of function pproximtion discussed in Chpter??, i.e., function pproximtion vi the

More information

ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS

ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS F. Tkeo 1 nd M. Sk 1 Hchinohe Ntionl College of Technology, Hchinohe, Jpn; Tohoku University, Sendi, Jpn Abstrct:

More information

CHAPTER 20: Second Law of Thermodynamics

CHAPTER 20: Second Law of Thermodynamics CHAER 0: Second Lw of hermodynmics Responses to Questions 3. kg of liquid iron will hve greter entropy, since it is less ordered thn solid iron nd its molecules hve more therml motion. In ddition, het

More information

Non-Linear & Logistic Regression

Non-Linear & Logistic Regression Non-Liner & Logistic Regression If the sttistics re boring, then you've got the wrong numbers. Edwrd R. Tufte (Sttistics Professor, Yle University) Regression Anlyses When do we use these? PART 1: find

More information

Riemann Sums and Riemann Integrals

Riemann Sums and Riemann Integrals Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 2013 Outline 1 Riemnn Sums 2 Riemnn Integrls 3 Properties

More information

Simulated Performance of Packed Bed Solar Energy Storage System having Storage Material Elements of Large Size - Part I

Simulated Performance of Packed Bed Solar Energy Storage System having Storage Material Elements of Large Size - Part I The Open Fuels & Energy Science Journl, 2008, 1, 91-96 91 Open Access Simulted Performnce of Pcked Bed Solr Energy Storge System hving Storge Mteril Elements of Lrge Size - Prt I Rnjit Singh *,1, R.P.

More information

The Predom module. Predom calculates and plots isothermal 1-, 2- and 3-metal predominance area diagrams. Predom accesses only compound databases.

The Predom module. Predom calculates and plots isothermal 1-, 2- and 3-metal predominance area diagrams. Predom accesses only compound databases. Section 1 Section 2 The module clcultes nd plots isotherml 1-, 2- nd 3-metl predominnce re digrms. ccesses only compound dtbses. Tble of Contents Tble of Contents Opening the module Section 3 Stoichiometric

More information

A5682: Introduction to Cosmology Course Notes. 4. Cosmic Dynamics: The Friedmann Equation. = GM s

A5682: Introduction to Cosmology Course Notes. 4. Cosmic Dynamics: The Friedmann Equation. = GM s 4. Cosmic Dynmics: The Friedmnn Eqution Reding: Chpter 4 Newtonin Derivtion of the Friedmnn Eqution Consider n isolted sphere of rdius R s nd mss M s, in uniform, isotropic expnsion (Hubble flow). The

More information

CS 188 Introduction to Artificial Intelligence Fall 2018 Note 7

CS 188 Introduction to Artificial Intelligence Fall 2018 Note 7 CS 188 Introduction to Artificil Intelligence Fll 2018 Note 7 These lecture notes re hevily bsed on notes originlly written by Nikhil Shrm. Decision Networks In the third note, we lerned bout gme trees

More information

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007 A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus

More information

Acceptance Sampling by Attributes

Acceptance Sampling by Attributes Introduction Acceptnce Smpling by Attributes Acceptnce smpling is concerned with inspection nd decision mking regrding products. Three spects of smpling re importnt: o Involves rndom smpling of n entire

More information

Predict Global Earth Temperature using Linier Regression

Predict Global Earth Temperature using Linier Regression Predict Globl Erth Temperture using Linier Regression Edwin Swndi Sijbt (23516012) Progrm Studi Mgister Informtik Sekolh Teknik Elektro dn Informtik ITB Jl. Gnesh 10 Bndung 40132, Indonesi 23516012@std.stei.itb.c.id

More information

Lecture 1. Functional series. Pointwise and uniform convergence.

Lecture 1. Functional series. Pointwise and uniform convergence. 1 Introduction. Lecture 1. Functionl series. Pointwise nd uniform convergence. In this course we study mongst other things Fourier series. The Fourier series for periodic function f(x) with period 2π is

More information

Terminal Velocity and Raindrop Growth

Terminal Velocity and Raindrop Growth Terminl Velocity nd Rindrop Growth Terminl velocity for rindrop represents blnce in which weight mss times grvity is equl to drg force. F 3 π3 ρ L g in which is drop rdius, g is grvittionl ccelertion,

More information

Riemann Sums and Riemann Integrals

Riemann Sums and Riemann Integrals Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 203 Outline Riemnn Sums Riemnn Integrls Properties Abstrct

More information

Applications of Bernoulli s theorem. Lecture - 7

Applications of Bernoulli s theorem. Lecture - 7 Applictions of Bernoulli s theorem Lecture - 7 Prcticl Applictions of Bernoulli s Theorem The Bernoulli eqution cn be pplied to gret mny situtions not just the pipe flow we hve been considering up to now.

More information

Monte Carlo method in solving numerical integration and differential equation

Monte Carlo method in solving numerical integration and differential equation Monte Crlo method in solving numericl integrtion nd differentil eqution Ye Jin Chemistry Deprtment Duke University yj66@duke.edu Abstrct: Monte Crlo method is commonly used in rel physics problem. The

More information

On the Linear Stability of Compound Capillary Jets

On the Linear Stability of Compound Capillary Jets ILASS Americs, th Annul Conference on Liquid Atomiztion nd Spry Systems, Chicgo, IL, My 7 On the Liner Stbility of Compound Cpillry Jets Mksud (Mx) Ismilov, Stephen D Heister School of Aeronutics nd Astronutics,

More information

Probabilistic Investigation of Sensitivities of Advanced Test- Analysis Model Correlation Methods

Probabilistic Investigation of Sensitivities of Advanced Test- Analysis Model Correlation Methods Probbilistic Investigtion of Sensitivities of Advnced Test- Anlysis Model Correltion Methods Liz Bergmn, Mtthew S. Allen, nd Dniel C. Kmmer Dept. of Engineering Physics University of Wisconsin-Mdison Rndll

More information

( dg. ) 2 dt. + dt. dt j + dh. + dt. r(t) dt. Comparing this equation with the one listed above for the length of see that

( dg. ) 2 dt. + dt. dt j + dh. + dt. r(t) dt. Comparing this equation with the one listed above for the length of see that Arc Length of Curves in Three Dimensionl Spce If the vector function r(t) f(t) i + g(t) j + h(t) k trces out the curve C s t vries, we cn mesure distnces long C using formul nerly identicl to one tht we

More information

Section 6.1 INTRO to LAPLACE TRANSFORMS

Section 6.1 INTRO to LAPLACE TRANSFORMS Section 6. INTRO to LAPLACE TRANSFORMS Key terms: Improper Integrl; diverge, converge A A f(t)dt lim f(t)dt Piecewise Continuous Function; jump discontinuity Function of Exponentil Order Lplce Trnsform

More information

CALCULUS WITHOUT LIMITS

CALCULUS WITHOUT LIMITS CALCULUS WITHOUT LIMITS The current stndrd for the clculus curriculum is, in my opinion, filure in mny spects. We try to present it with the modern stndrd of mthemticl rigor nd comprehensiveness but of

More information

I1 = I2 I1 = I2 + I3 I1 + I2 = I3 + I4 I 3

I1 = I2 I1 = I2 + I3 I1 + I2 = I3 + I4 I 3 2 The Prllel Circuit Electric Circuits: Figure 2- elow show ttery nd multiple resistors rrnged in prllel. Ech resistor receives portion of the current from the ttery sed on its resistnce. The split is

More information

1.9 C 2 inner variations

1.9 C 2 inner variations 46 CHAPTER 1. INDIRECT METHODS 1.9 C 2 inner vritions So fr, we hve restricted ttention to liner vritions. These re vritions of the form vx; ǫ = ux + ǫφx where φ is in some liner perturbtion clss P, for

More information

Chapter 0. What is the Lebesgue integral about?

Chapter 0. What is the Lebesgue integral about? Chpter 0. Wht is the Lebesgue integrl bout? The pln is to hve tutoril sheet ech week, most often on Fridy, (to be done during the clss) where you will try to get used to the ides introduced in the previous

More information

Lecture 19: Continuous Least Squares Approximation

Lecture 19: Continuous Least Squares Approximation Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for

More information

The Regulated and Riemann Integrals

The Regulated and Riemann Integrals Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue

More information

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0) 1 Tylor polynomils In Section 3.5, we discussed how to pproximte function f(x) round point in terms of its first derivtive f (x) evluted t, tht is using the liner pproximtion f() + f ()(x ). We clled this

More information

1. a) Describe the principle characteristics and uses of the following types of PV cell: Single Crystal Silicon Poly Crystal Silicon

1. a) Describe the principle characteristics and uses of the following types of PV cell: Single Crystal Silicon Poly Crystal Silicon 2001 1. ) Describe the principle chrcteristics nd uses of the following types of PV cell: Single Crystl Silicon Poly Crystl Silicon Amorphous Silicon CIS/CIGS Gllium Arsenide Multijunction (12 mrks) b)

More information

Purpose of the experiment

Purpose of the experiment Newton s Lws II PES 6 Advnced Physics Lb I Purpose of the experiment Exmine two cses using Newton s Lws. Sttic ( = 0) Dynmic ( 0) fyi fyi Did you know tht the longest recorded flight of chicken is thirteen

More information

Lecture 6: Singular Integrals, Open Quadrature rules, and Gauss Quadrature

Lecture 6: Singular Integrals, Open Quadrature rules, and Gauss Quadrature Lecture notes on Vritionl nd Approximte Methods in Applied Mthemtics - A Peirce UBC Lecture 6: Singulr Integrls, Open Qudrture rules, nd Guss Qudrture (Compiled 6 August 7) In this lecture we discuss the

More information

Lecture 13 - Linking E, ϕ, and ρ

Lecture 13 - Linking E, ϕ, and ρ Lecture 13 - Linking E, ϕ, nd ρ A Puzzle... Inner-Surfce Chrge Density A positive point chrge q is locted off-center inside neutrl conducting sphericl shell. We know from Guss s lw tht the totl chrge on

More information

Driving Cycle Construction of City Road for Hybrid Bus Based on Markov Process Deng Pan1, a, Fengchun Sun1,b*, Hongwen He1, c, Jiankun Peng1, d

Driving Cycle Construction of City Road for Hybrid Bus Based on Markov Process Deng Pan1, a, Fengchun Sun1,b*, Hongwen He1, c, Jiankun Peng1, d Interntionl Industril Informtics nd Computer Engineering Conference (IIICEC 15) Driving Cycle Construction of City Rod for Hybrid Bus Bsed on Mrkov Process Deng Pn1,, Fengchun Sun1,b*, Hongwen He1, c,

More information

First Law of Thermodynamics. Control Mass (Closed System) Conservation of Mass. Conservation of Energy

First Law of Thermodynamics. Control Mass (Closed System) Conservation of Mass. Conservation of Energy First w of hermodynmics Reding Problems 3-3-7 3-0, 3-5, 3-05 5-5- 5-8, 5-5, 5-9, 5-37, 5-0, 5-, 5-63, 5-7, 5-8, 5-09 6-6-5 6-, 6-5, 6-60, 6-80, 6-9, 6-, 6-68, 6-73 Control Mss (Closed System) In this section

More information

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz

More information

221B Lecture Notes WKB Method

221B Lecture Notes WKB Method Clssicl Limit B Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using

More information

Studies on Nuclear Fuel Rod Thermal Performance

Studies on Nuclear Fuel Rod Thermal Performance Avilble online t www.sciencedirect.com Energy Procedi 1 (1) 1 17 Studies on Nucler Fuel od herml Performnce Eskndri, M.1; Bvndi, A ; Mihndoost, A3* 1 Deprtment of Physics, Islmic Azd University, Shirz

More information

1 Online Learning and Regret Minimization

1 Online Learning and Regret Minimization 2.997 Decision-Mking in Lrge-Scle Systems My 10 MIT, Spring 2004 Hndout #29 Lecture Note 24 1 Online Lerning nd Regret Minimiztion In this lecture, we consider the problem of sequentil decision mking in

More information

Review of basic calculus

Review of basic calculus Review of bsic clculus This brief review reclls some of the most importnt concepts, definitions, nd theorems from bsic clculus. It is not intended to tech bsic clculus from scrtch. If ny of the items below

More information

Numerical Integration

Numerical Integration Chpter 1 Numericl Integrtion Numericl differentition methods compute pproximtions to the derivtive of function from known vlues of the function. Numericl integrtion uses the sme informtion to compute numericl

More information

ODE: Existence and Uniqueness of a Solution

ODE: Existence and Uniqueness of a Solution Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =

More information

Numerical Analysis: Trapezoidal and Simpson s Rule

Numerical Analysis: Trapezoidal and Simpson s Rule nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =

More information

Euler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), )

Euler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), ) Euler, Iochimescu nd the trpezium rule G.J.O. Jmeson (Mth. Gzette 96 (0), 36 4) The following results were estblished in recent Gzette rticle [, Theorems, 3, 4]. Given > 0 nd 0 < s

More information

Consequently, the temperature must be the same at each point in the cross section at x. Let:

Consequently, the temperature must be the same at each point in the cross section at x. Let: HW 2 Comments: L1-3. Derive the het eqution for n inhomogeneous rod where the therml coefficients used in the derivtion of the het eqution for homogeneous rod now become functions of position x in the

More information

Lecture 14: Quadrature

Lecture 14: Quadrature Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl

More information

Introduction to Finite Element Method

Introduction to Finite Element Method Introduction to Finite Element Method Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pn.pl/ tzielins/ Tble of Contents 1 Introduction 1 1.1 Motivtion nd generl concepts.............

More information

Classical Mechanics. From Molecular to Con/nuum Physics I WS 11/12 Emiliano Ippoli/ October, 2011

Classical Mechanics. From Molecular to Con/nuum Physics I WS 11/12 Emiliano Ippoli/ October, 2011 Clssicl Mechnics From Moleculr to Con/nuum Physics I WS 11/12 Emilino Ippoli/ October, 2011 Wednesdy, October 12, 2011 Review Mthemtics... Physics Bsic thermodynmics Temperture, idel gs, kinetic gs theory,

More information

Physics 201 Lab 3: Measurement of Earth s local gravitational field I Data Acquisition and Preliminary Analysis Dr. Timothy C. Black Summer I, 2018

Physics 201 Lab 3: Measurement of Earth s local gravitational field I Data Acquisition and Preliminary Analysis Dr. Timothy C. Black Summer I, 2018 Physics 201 Lb 3: Mesurement of Erth s locl grvittionl field I Dt Acquisition nd Preliminry Anlysis Dr. Timothy C. Blck Summer I, 2018 Theoreticl Discussion Grvity is one of the four known fundmentl forces.

More information

Operations with Polynomials

Operations with Polynomials 38 Chpter P Prerequisites P.4 Opertions with Polynomils Wht you should lern: How to identify the leding coefficients nd degrees of polynomils How to dd nd subtrct polynomils How to multiply polynomils

More information

8 Laplace s Method and Local Limit Theorems

8 Laplace s Method and Local Limit Theorems 8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved

More information

Estimation of Binomial Distribution in the Light of Future Data

Estimation of Binomial Distribution in the Light of Future Data British Journl of Mthemtics & Computer Science 102: 1-7, 2015, Article no.bjmcs.19191 ISSN: 2231-0851 SCIENCEDOMAIN interntionl www.sciencedomin.org Estimtion of Binomil Distribution in the Light of Future

More information

Conservation Law. Chapter Goal. 5.2 Theory

Conservation Law. Chapter Goal. 5.2 Theory Chpter 5 Conservtion Lw 5.1 Gol Our long term gol is to understnd how mny mthemticl models re derived. We study how certin quntity chnges with time in given region (sptil domin). We first derive the very

More information

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying Vitli covers 1 Definition. A Vitli cover of set E R is set V of closed intervls with positive length so tht, for every δ > 0 nd every x E, there is some I V with λ(i ) < δ nd x I. 2 Lemm (Vitli covering)

More information

Chapters 4 & 5 Integrals & Applications

Chapters 4 & 5 Integrals & Applications Contents Chpters 4 & 5 Integrls & Applictions Motivtion to Chpters 4 & 5 2 Chpter 4 3 Ares nd Distnces 3. VIDEO - Ares Under Functions............................................ 3.2 VIDEO - Applictions

More information

Markscheme May 2016 Mathematics Standard level Paper 1

Markscheme May 2016 Mathematics Standard level Paper 1 M6/5/MATME/SP/ENG/TZ/XX/M Mrkscheme My 06 Mthemtics Stndrd level Pper 7 pges M6/5/MATME/SP/ENG/TZ/XX/M This mrkscheme is the property of the Interntionl Bcclurete nd must not be reproduced or distributed

More information

Lecture 20: Numerical Integration III

Lecture 20: Numerical Integration III cs4: introduction to numericl nlysis /8/0 Lecture 0: Numericl Integrtion III Instructor: Professor Amos Ron Scribes: Mrk Cowlishw, Yunpeng Li, Nthnel Fillmore For the lst few lectures we hve discussed

More information

Minimum Energy State of Plasmas with an Internal Transport Barrier

Minimum Energy State of Plasmas with an Internal Transport Barrier Minimum Energy Stte of Plsms with n Internl Trnsport Brrier T. Tmno ), I. Ktnum ), Y. Skmoto ) ) Formerly, Plsm Reserch Center, University of Tsukub, Tsukub, Ibrki, Jpn ) Plsm Reserch Center, University

More information

Rel Gses 1. Gses (N, CO ) which don t obey gs lws or gs eqution P=RT t ll pressure nd tempertures re clled rel gses.. Rel gses obey gs lws t extremely low pressure nd high temperture. Rel gses devited

More information

ARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac

ARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac REVIEW OF ALGEBRA Here we review the bsic rules nd procedures of lgebr tht you need to know in order to be successful in clculus. ARITHMETIC OPERATIONS The rel numbers hve the following properties: b b

More information

INTRODUCTION TO INTEGRATION

INTRODUCTION TO INTEGRATION INTRODUCTION TO INTEGRATION 5.1 Ares nd Distnces Assume f(x) 0 on the intervl [, b]. Let A be the re under the grph of f(x). b We will obtin n pproximtion of A in the following three steps. STEP 1: Divide

More information