Mathematical Modeling of Self - Oscillations in the Combustion Chamber of Liquid Rocket Engine with Variable Latency Combustion
|
|
- Gavin Small
- 5 years ago
- Views:
Transcription
1 Physics Joural Vol., No. 3, 05, pp Mathematical Modelig of Self - Oscillatios i the Combustio Chamber of Liquid Rocket Egie with Variable Latecy Combustio B. I. Basok, V. V. otsuleko * Departmet of Thermophysical Fudametals of Eergy-Savig Techologies, Istitute of Egieerig Thermal Physics of NAS of Ukraie, Kiev, Ukraie Abstract Self-oscillatios ad certai of their regularities determied by solutio of a system of differetial equatios with variable delay argumet equatios that is used i cosiderig combustio istability i combustio chambers of liquid-propellat rocket egies are modeled mathematically. Periodic solutios of the system of equatios of ostatioary motio of a medium i a liquid-propellat rocket egie were obtaied, with the aid of which the possibility of lowerig the amplitude of the logitudial self-oscillatios of vibratio combustio or their complete removal has bee substatiated. Aalytically determied critical time delay combustio, above which a statioary combustio becomes ustable ad self-excited oscillatios. Keywords Vibratio Combustio, Istability, Self - Oscillatio, Limitig Cycle, the Time Delay of Combustio Received: September 6, 05 / Accepted: December, 05 / Published olie: December 9, 05 The Authors. Published by America Istitute of Sciece. This Ope Access article is uder the CC BY-NC licese. Itroductio The ostatioarity of the process of burig of a fuel i differet devices is geerated by the chemical-kietic, diffusio-thermal, ad covective pheomea []. Furthermore, accordig to the first low of thermo- dyamics, the heat released i combustio of the fuel i the flow is coverted to the iteral eergy ad compoets of the total head whose value is depedet o the flow rate. The ecessary coditio of excitatio of self-oscillatios i the case i questio is the presece of the ascedig brach o the head characteristic [, 3], which is formed by the heat-tohead coversio. I [4], it has bee established that, i burig of the fuel i solid-fuel egies, the pressure depedece of the rate of formatio of the gas ca be such that the system will lose stability [5]. I liquid-propellat rocket egies (LPRE, the coditio of itrachamber istability lies i the formatio of the ascedig brach of the p = F ( depedece of the pressure i the combustio chamber o the flow rate. rowth i the head or the pressure p i the combustio chamber with flow rate features promietly amog the causes of excitatio of vibratioal-combustio self-oscillatios; it has ot bee cosidered i theoretical descriptio of the pheomeo of istability of combustio i liquid-propellat rocket egies. The reaso is that the coditios of formatio of the ascedig brach o the p = F ( depedece of the head characteristic of the combustio chamber o the flow rate remaied ukow. Thus, e.g., i [5, 6], the above depedece of the head characteristic has bee represeted as beig mootoically decreasig. I [5], the stability of the statioary regime with such a depedece has bee substatiated. Therefore, combustio istability i liquid-propellat rocket egies * Correspodig author address: basok@ittf.kiev.ua (B. I. Basok, gosul@ukr.et (V. V. otsuleko
2 344 B. I. Basok ad V. V. otsuleko: Mathematical Modelig of Self - Oscillatios i the Combustio Chamber of Liquid Rocket Egie with Variable Latecy Combustio was maily determied by the mechaism of pheomeological time delay of the combustio of the fuel [5 7] ad by the maifestatio of differet feedbacks which are iteral ad act uiversally.. Formulatio of the problem I this paper, we cosider the problem of determiig the critical time delay of combustio ad the costructio of the boudaries of the stability regio whe the delay is variable, depedig o the pressure iside the combustio chamber. I Fig. shows a diagram of the combustio chamber rocket egie. 3. System of Equatios of Itrachamber Combustio Istability The equatio of motio for the sectios of the flow of a fuel 0 0 ad of combustio products 3 3 (Fig. i the form of d mw = p h h p S, takig ito accout that the f T fr mass of the combustio is m = ρls, will be writte as [8, 9] where La,c.ch L a,c.ch d = F ( p, ( = l S acoustic mass of the combustio chamber; mass flow rate i the combustio chamber; f p = H pressure characteristic of parallel coectio cetrifugal pumps oxidat ad fuel; F = p h h pressure head characteristic f T fr of fuel flow i a combustio chamber [8-9]; ht ( pressure losses because of heat supply (thermal resistace; hfr ( viscous losses alog the legth of the combustio chamber. The mass equatio i the combustio chamber is represeted, accordig to [8], i the form ( τ ( dm = t p, where = f + Cox, f flow rate of a fuel, C ox flow rate of oxidizer, τ ( p the time delay of combustio of fuel, mass flow rate through the ozzle. Sice dp d ρ is equal to c, where c velocity of soud i the flow, it may be reduced to the form dp Ca,c.ch = t τ ( p φ ( p, ( where φ ( p = characteristic of the ozzle of a liquidpropellat rocket egie, Ca,c.ch = V c acoustic flexibility of the combustio chamber, V the volume of the combustio chamber. Experimetal studies [7] showed that the depedece of τ ( p is mootoically decreasig. Figure shows the experimetal poits ad their approximatio usig the followig formula ( p bexp τ = ap p, where a = 0., b = Calculatio of Hydraulic Characteristics i the LPRE Combustio Chamber The viscous losses occurrig over the combustio chamber legth i the smoke gas motio segmet were calculated usig the Darcy - Weisbach formula l w hfr ( = λ ρ, (3 d where λ is the coefficiet of hydraulic losses. Heat iput to the fuel flow with a mass rate of flow = ρws i the chael of the LPRE combustio chamber with a costat ormal cross-sectio area (Fig. leads to a decrease i the desity ρ of the flow of combustio products ad a icrease i its velocity w. Therefore, the heat resistace to such flow [0] is determied by the local hydraulic resistace i the regio of heat iput. I [0], we obtaied a fuctio for determiig the eergy loss i the combustio chamber ht ( at polytropic heat iput to the fuel flow. Let us defie the fuctio ht ( at isobaric heat iput to the fuel flow i the LPRE combustio chamber. We write the eergy equatio for sectios ad (Fig., ad from this equatio we determie the thermal resistace, which is the local oe, appearig i the regio of heat supply [0] p w p w q + u u h ρ + + = + + +, (4 T ρ where ht is the eergy loss i the flow due to the heat iput. Sice at isobaric combustio heat iput to the ideal gas flow q = сp( T T, the chage i the iteral eergy of the u u = с ( T T, ad cp cv = R, the gas v q u = RT ( T. Further, i view of the relatio q u = p ρ p ρ from Eq. (4 it follows that
3 Physics Joural Vol., No. 3, 05, pp h = w w, or i pressure uits, assumig h T = ρ h, we have T T h w w = ρ. w T Makig use of the cotiuity equatio ρw = ρw for the gas flow i the combustio chamber with a costat ormal cross-sectio area ad of the relatio betwee the isobaric process parameters T T = ρ ρ, we obtai [] shows the graphics data, whe the pressure characteristic F ( of the combustio chamber is mootoically decreasig. 5. Calculatig the Boudaries of the Stable Combustio The parameters of statioary combustio mode are determied from the system of equatios ( ( believig i it T T = ρs T h. (5 Determie the formula for the hydraulic characteristics of the jet ozzle. The maximum mass flow rate of gas i the critical sectio of ozzle or where φ = + Smi k k = + p Smik k ( p k + k p w, = φ, (6 k + ( k p krt ; R gas costat of combustio products; k adiabatic expoet; T absolute temperature of the flowig gas. Foud above depedecies (3, (5, (6 makes the system of equatios ( ( fully defied. Itegratig this system of equatios allows to determie the limit cycles ad selfoscillatios i the combustio chamber rocket egie. I Fig. 3 shows the limit cycles ad their correspodig forms of self-oscillatios i a saddle-ode discharge characteristics F ( of the combustio chamber. Respectively i Fig. 4 d dp = 0, = 0. (7 = p = p Usig coditios (7, we obtai p F ( = φ ( p. Also Usig the Taylor decompositio d ( t τ = ( t τ + O ( τ, = ad the system of equatios ( - (, accurate to quatities of order O ( τ, ca be writte i the followig form L a,c.ch d = F ( p, ( p ( φ dp τ C = F p p, (8 a,c.ch La Further, it is more coveiet to switch to dimesioless variables: x =, y p p =. (9 p I the ew variables (9 the system of equatios (8 ca be writte i the followig form: dx La,c.ch = F ( + x p ( y +, dy F ( + x p ( y + pc a,c.ch = ( x + τ ( p ( y + φ ( p ( y +. La,c.ch (0 Thus, the ature of stability of statioary combustio mode is reduced to the study of stability of zero equilibrium of the dyamic system (0. Usig Lyapuov's first method of stability aalysis let us cosider the Jacobia matrix of the system (0 is calculated i the zero positio of equilibrium F p La La,c.ch J =. ( τ ( p τ ( p F ( φ ( p pc a,c.ch La,c.ch Ca,c.ch La,c.ch
4 346 B. I. Basok ad V. V. otsuleko: Mathematical Modelig of Self - Oscillatios i the Combustio Chamber of Liquid Rocket Egie with Variable Latecy Combustio df where ( + x F ( =, φ ( p dx x= 0 ( + dφ p p y =. dy To determie the critical time delay of combustio is ecessary to pre-compute the roots of the characteristic equatio ( λ y= 0 det J E = 0. ( Calculatig the determiat of (, we obtai ± i tr J 4det J tr J λ, =, where tr ( J trace, ad det ( J the determiat of the Jacobia matrix J. Moreover, accordig to (, obtai the followig represetatios for these characteristics: τ ( p F tr ( J = + φ ( p, La,c.ch Ca,c.ch La,c.ch det ( J = ( F ( φ ( p. L C a,c.ch a,c.ch Thus, the critical delay time of combustio is determied from the coditio: {,} Re λ = 0 ( J ( J tr = 0, det > 0. F τ + = 0, ( p kr φ La,c.ch C a,c.ch L a,c.ch F ( φ ( p <, depedece for the critical time delay from combustiochamber pressure p τ 3 3 L a,c.ch F0 p F0 p kr ( p = 3kC F a,c.ch p kf kf. (4 I Fig. 5 shows a plot of (4 to the combustio chamber legth l = 0.5 m ad diameter d = 0. m, whe as the oxidizig aget used liquid oxyge ad fuel hydroge. 6. Coclusios By umerically itegratig the system of equatios of the ostatioary motio of a medium i a LPRE, i trasitio of it to a degeerate form, for the first time the possibility of cotrollig the amplitude of oscillatios of vibratio combustio at differet pressure head characteristic of a combustio chamber has bee show. The possibility of complete suppressio of self-oscillatios ad attaimet of a statioary regime of combustio i a LPRE has bee established. Aalytically derived the ratio that determies the critical time delay of combustio of gaseous fuel that must be exceeded i the combustio chamber of rocket egies statioary combustio mode becomes ustable ad excited vibratio combustio. Also for the combustio chamber with a mootoically decreasig pressure characteristic was obtaied depedece of the critical delay time of combustio. where fially, we obtai that ( p C F ( τ = L φ. (3 kr a,c.ch a,c.ch Thus, the equilibrium of the dyamic system (0 that determies the parameters of statioary combustio mode becomes ustable if true the followig iequality τ p τ kr >. < Note also that the coditio ( ( p F φ is executed automatically whe pressure characteristic F ( of the combustio chamber is a mootoically decreasig fuctio. For further aalysis of the obtaied formula (3 we well approximate the pressure characteristic of the combustio chamber by a polyomial of third degree, equatig 3 F F k. The from (3 we get the followig 0 F Fig.. Positio of the sectios i the LPRE combustio chamber that are used for costructio of the equatios of motio of combustio products. Fig.. The depedece of the time delay from combustio chamber pressure.
5 Physics Joural Vol., No. 3, 05, pp Fig. 3. Limit cycles ad self - oscillatios of vibratio combustio whe the characteristic c = 3.5 kg sec ; d = 4. kg sec F is a saddle-ode: a = kg sec ; b =.5 kg sec ; Fig. 4. Limit cycles ad self - oscillatios of vibratio combustio whe the characteristic =4 kg sec ; c = 5 kg sec. F is mootoically decreasig: a = kg sec ; b Refereces [] Ya. B. Zel dovich,. I. Bareblatt, V. B. Librovich, ad. M. Makhviladze, Mathematical Theory of Combustio ad Explosio [i Russia], Nauka, Moscow (980. [] V. V. otsuleko, Special modes of the Rijke pheomeo, Izh.-Fiz. Zh., 6, No. 9, (005. [3] V. V. otsuleko, Mathematical modelig of the decrease i the amplitudes of vibratios of vibratioal combustio i large idustrial uits, Mat. Modelir., 7, No., 6 4 (005. Fig. 5. The depedece of the critical time delay of combustio τkr ( p. [4] Ya. B. Zel dovich, O. I. Leipuskii, ad V. B. Librovich, Theory of Nostatioary Combustio of upowder [i Russia], Nauka, Moscow (975.
6 348 B. I. Basok ad V. V. otsuleko: Mathematical Modelig of Self - Oscillatios i the Combustio Chamber of Liquid Rocket Egie with Variable Latecy Combustio [5] M. S. Natazo, Istability of Combustio [i Russia], Mashiostroeie, Moscow (986. [6] K. I. Artamoov, Thermoacoustic Stability [i Russia], Mashiostroeie, Moscow (98. [7] L. Crocco ad Si-I Cheg, Theory of Combustio Istability i Liquid-Propellat Rocket Egies [Russia traslatio], IL, Moscow (958. [8] V. V. otsuleko ad V. N. otsuleko, Mathematical modelig of the self-oscillatios of vibratioal combustio i liquid-propellat rocket egies due to combustio heat release, Mat. Modelyuv. (Deprodzerzhisk os. Tekh. Uiv., No., (5, 6 4 (006. [9] otsuleko V. V. ad otsuleko V. N. O the Problem of Cotrol of the Amplitude of Self-Oscillatios of a Sigig Flame // Joural of Egieerig Physics ad Thermophysics, Volume 87, Issue, (04. [0] B. I. Basok ad V.V. otsuleko, Negative thermal resistace i oe-dimesioal steady iviscid flow of a perfect gas, Joural of the Moscow physical-techical Istitute, Vol 6, No. 4(4, (04. [] Basok B. I. ad otsuleko V.V. Calculatig the Parameters of Self - Oscillatios i the Vertical Combustio Chamber of the Blast - Furace Air Heater durig Ustable Combustio // Thermal Egieerig, Vol. 6, No., (05.
Damped Vibration of a Non-prismatic Beam with a Rotational Spring
Vibratios i Physical Systems Vol.6 (0) Damped Vibratio of a No-prismatic Beam with a Rotatioal Sprig Wojciech SOCHACK stitute of Mechaics ad Fudametals of Machiery Desig Uiversity of Techology, Czestochowa,
More informationPOWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS
Joural of Applied Mathematics ad Computatioal Mechaics 4 3(3) 3-8 POWER SERIES SOLUION OF FIRS ORDER MARIX DIFFERENIAL EQUAIONS Staisław Kukla Izabela Zamorska Istitute of Mathematics Czestochowa Uiversity
More informationDETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.
ICSV4 Cairs Australia 9- July 7 DTRMINATION OF MCHANICAL PROPRTIS OF A NON- UNIFORM BAM USING TH MASURMNT OF TH XCITD LONGITUDINAL LASTIC VIBRATIONS Pavel Aokhi ad Vladimir Gordo Departmet of the mathematics
More informationComparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series
Applied Mathematical Scieces, Vol. 7, 03, o. 6, 3-337 HIKARI Ltd, www.m-hikari.com http://d.doi.org/0.988/ams.03.3430 Compariso Study of Series Approimatio ad Covergece betwee Chebyshev ad Legedre Series
More informationKinetics of Complex Reactions
Kietics of Complex Reactios by Flick Colema Departmet of Chemistry Wellesley College Wellesley MA 28 wcolema@wellesley.edu Copyright Flick Colema 996. All rights reserved. You are welcome to use this documet
More informationBoundary layer problem on conveyor belt. Gabriella Bognár University of Miskolc 3515 Miskolc-Egyetemváros, Hungary
Boudary layer problem o coveyor belt Gabriella Bogár Uiversity of Miskolc 355 Miskolc-Egyetemváros, Hugary e-mail: matvbg@ui-miskolc.hu Abstract: A techologically importat source of the boudary layer pheomeo
More information2C09 Design for seismic and climate changes
2C09 Desig for seismic ad climate chages Lecture 02: Dyamic respose of sigle-degree-of-freedom systems I Daiel Grecea, Politehica Uiversity of Timisoara 10/03/2014 Europea Erasmus Mudus Master Course Sustaiable
More informationNumerical Methods in Fourier Series Applications
Numerical Methods i Fourier Series Applicatios Recall that the basic relatios i usig the Trigoometric Fourier Series represetatio were give by f ( x) a o ( a x cos b x si ) () where the Fourier coefficiets
More informationFree Surface Hydrodynamics
Water Sciece ad Egieerig Free Surface Hydrodyamics y A part of Module : Hydraulics ad Hydrology Water Sciece ad Egieerig Dr. Shreedhar Maskey Seior Lecturer UNESCO-IHE Istitute for Water Educatio S. Maskey
More informationStreamfunction-Vorticity Formulation
Streamfuctio-Vorticity Formulatio A. Salih Departmet of Aerospace Egieerig Idia Istitute of Space Sciece ad Techology, Thiruvaathapuram March 2013 The streamfuctio-vorticity formulatio was amog the first
More information1. Hydrogen Atom: 3p State
7633A QUANTUM MECHANICS I - solutio set - autum. Hydroge Atom: 3p State Let us assume that a hydroge atom is i a 3p state. Show that the radial part of its wave fuctio is r u 3(r) = 4 8 6 e r 3 r(6 r).
More informationEXPERIMENT OF SIMPLE VIBRATION
EXPERIMENT OF SIMPLE VIBRATION. PURPOSE The purpose of the experimet is to show free vibratio ad damped vibratio o a system havig oe degree of freedom ad to ivestigate the relatioship betwee the basic
More informationTaylor polynomial solution of difference equation with constant coefficients via time scales calculus
TMSCI 3, o 3, 129-135 (2015) 129 ew Treds i Mathematical Scieces http://wwwtmscicom Taylor polyomial solutio of differece equatio with costat coefficiets via time scales calculus Veysel Fuat Hatipoglu
More informationAN OPEN-PLUS-CLOSED-LOOP APPROACH TO SYNCHRONIZATION OF CHAOTIC AND HYPERCHAOTIC MAPS
http://www.paper.edu.c Iteratioal Joural of Bifurcatio ad Chaos, Vol. 1, No. 5 () 119 15 c World Scietific Publishig Compay AN OPEN-PLUS-CLOSED-LOOP APPROACH TO SYNCHRONIZATION OF CHAOTIC AND HYPERCHAOTIC
More informationMETHOD OF FUNDAMENTAL SOLUTIONS FOR HELMHOLTZ EIGENVALUE PROBLEMS IN ELLIPTICAL DOMAINS
Please cite this article as: Staisław Kula, Method of fudametal solutios for Helmholtz eigevalue problems i elliptical domais, Scietific Research of the Istitute of Mathematics ad Computer Sciece, 009,
More information62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +
62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of
More informationDECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS. Park Road, Islamabad, Pakistan
Mathematical ad Computatioal Applicatios, Vol. 9, No. 3, pp. 30-40, 04 DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS Muhammad Aslam Noor, Khalida Iayat Noor ad Asif Waheed
More informationA numerical Technique Finite Volume Method for Solving Diffusion 2D Problem
The Iteratioal Joural Of Egieerig d Sciece (IJES) Volume 4 Issue 10 Pages PP -35-41 2015 ISSN (e): 2319 1813 ISSN (p): 2319 1805 umerical Techique Fiite Volume Method for Solvig Diffusio 2D Problem 1 Mohammed
More informationSimilarity between quantum mechanics and thermodynamics: Entropy, temperature, and Carnot cycle
Similarity betwee quatum mechaics ad thermodyamics: Etropy, temperature, ad Carot cycle Sumiyoshi Abe 1,,3 ad Shiji Okuyama 1 1 Departmet of Physical Egieerig, Mie Uiversity, Mie 514-8507, Japa Istitut
More informationResearch Article Nonautonomous Discrete Neuron Model with Multiple Periodic and Eventually Periodic Solutions
Discrete Dyamics i Nature ad Society Volume 21, Article ID 147282, 6 pages http://dx.doi.org/1.11/21/147282 Research Article Noautoomous Discrete Neuro Model with Multiple Periodic ad Evetually Periodic
More informationDynamical stability and parametrical vibrations of the laminated plates with complex shape
(3) 75 88 Dyamical stability ad parametrical vibratios of the lamiated plates with comple shape Abstract The problem of oliear vibratios ad stability aalysis for the symmetric lamiated plates with comple
More informationChapter 4. Fourier Series
Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,
More informationTime-Domain Representations of LTI Systems
2.1 Itroductio Objectives: 1. Impulse resposes of LTI systems 2. Liear costat-coefficiets differetial or differece equatios of LTI systems 3. Bloc diagram represetatios of LTI systems 4. State-variable
More informationRay Optics Theory and Mode Theory. Dr. Mohammad Faisal Dept. of EEE, BUET
Ray Optics Theory ad Mode Theory Dr. Mohammad Faisal Dept. of, BUT Optical Fiber WG For light to be trasmitted through fiber core, i.e., for total iteral reflectio i medium, > Ray Theory Trasmissio Ray
More information1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations
. Liearizatio of a oliear system give i the form of a system of ordiary differetial equatios We ow show how to determie a liear model which approximates the behavior of a time-ivariat oliear system i a
More informationExercises and Problems
HW Chapter 4: Oe-Dimesioal Quatum Mechaics Coceptual Questios 4.. Five. 4.4.. is idepedet of. a b c mu ( E). a b m( ev 5 ev) c m(6 ev ev) Exercises ad Problems 4.. Model: Model the electro as a particle
More informationANALYSIS OF DAMPING EFFECT ON BEAM VIBRATION
Molecular ad Quatum Acoustics vol. 7, (6) 79 ANALYSIS OF DAMPING EFFECT ON BEAM VIBRATION Jerzy FILIPIAK 1, Lech SOLARZ, Korad ZUBKO 1 Istitute of Electroic ad Cotrol Systems, Techical Uiversity of Czestochowa,
More informationThe Random Walk For Dummies
The Radom Walk For Dummies Richard A Mote Abstract We look at the priciples goverig the oe-dimesioal discrete radom walk First we review five basic cocepts of probability theory The we cosider the Beroulli
More information1 Adiabatic and diabatic representations
1 Adiabatic ad diabatic represetatios 1.1 Bor-Oppeheimer approximatio The time-idepedet Schrödiger equatio for both electroic ad uclear degrees of freedom is Ĥ Ψ(r, R) = E Ψ(r, R), (1) where the full molecular
More informationModified Decomposition Method by Adomian and. Rach for Solving Nonlinear Volterra Integro- Differential Equations
Noliear Aalysis ad Differetial Equatios, Vol. 5, 27, o. 4, 57-7 HIKARI Ltd, www.m-hikari.com https://doi.org/.2988/ade.27.62 Modified Decompositio Method by Adomia ad Rach for Solvig Noliear Volterra Itegro-
More information[ ] sin ( ) ( ) = 2 2 ( ) ( ) ( ) ˆ Mechanical Spectroscopy II
Solid State Pheomea Vol. 89 (003) pp 343-348 (003) Tras Tech Publicatios, Switzerlad doi:0.408/www.scietific.et/ssp.89.343 A New Impulse Mechaical Spectrometer to Study the Dyamic Mechaical Properties
More information(c) Write, but do not evaluate, an integral expression for the volume of the solid generated when R is
Calculus BC Fial Review Name: Revised 7 EXAM Date: Tuesday, May 9 Remiders:. Put ew batteries i your calculator. Make sure your calculator is i RADIAN mode.. Get a good ight s sleep. Eat breakfast. Brig:
More informationA Negative Result. We consider the resolvent problem for the scalar Oseen equation
O Osee Resolvet Estimates: A Negative Result Paul Deurig Werer Varhor 2 Uiversité Lille 2 Uiversität Kassel Laboratoire de Mathématiques BP 699, 62228 Calais cédex Frace paul.deurig@lmpa.uiv-littoral.fr
More informationProblem 1. Problem Engineering Dynamics Problem Set 9--Solution. Find the equation of motion for the system shown with respect to:
2.003 Egieerig Dyamics Problem Set 9--Solutio Problem 1 Fid the equatio of motio for the system show with respect to: a) Zero sprig force positio. Draw the appropriate free body diagram. b) Static equilibrium
More informationLecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods
Itroductory lecture otes o Partial Differetial Equatios - c Athoy Peirce. Not to be copied, used, or revised without explicit writte permissio from the copyright ower. 1 Lecture 8: Solvig the Heat, Laplace
More informationTHE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE ABSTRACT
Europea Joural of Egieerig ad Techology Vol. 3 No., 5 ISSN 56-586 THE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE Atif Nazir, Tahir Mahmood ad
More informationChemical Kinetics CHAPTER 14. Chemistry: The Molecular Nature of Matter, 6 th edition By Jesperson, Brady, & Hyslop. CHAPTER 14 Chemical Kinetics
Chemical Kietics CHAPTER 14 Chemistry: The Molecular Nature of Matter, 6 th editio By Jesperso, Brady, & Hyslop CHAPTER 14 Chemical Kietics Learig Objectives: Factors Affectig Reactio Rate: o Cocetratio
More informationTrue Nature of Potential Energy of a Hydrogen Atom
True Nature of Potetial Eergy of a Hydroge Atom Koshu Suto Key words: Bohr Radius, Potetial Eergy, Rest Mass Eergy, Classical Electro Radius PACS codes: 365Sq, 365-w, 33+p Abstract I cosiderig the potetial
More informationFinite Difference Derivations for Spreadsheet Modeling John C. Walton Modified: November 15, 2007 jcw
Fiite Differece Derivatios for Spreadsheet Modelig Joh C. Walto Modified: November 15, 2007 jcw Figure 1. Suset with 11 swas o Little Platte Lake, Michiga. Page 1 Modificatio Date: November 15, 2007 Review
More informationSimilarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall
Iteratioal Mathematical Forum, Vol. 9, 04, o. 3, 465-475 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/imf.04.48 Similarity Solutios to Usteady Pseudoplastic Flow Near a Movig Wall W. Robi Egieerig
More informationAnalytic Theory of Probabilities
Aalytic Theory of Probabilities PS Laplace Book II Chapter II, 4 pp 94 03 4 A lottery beig composed of umbered tickets of which r exit at each drawig, oe requires the probability that after i drawigs all
More informationTHE KALMAN FILTER RAUL ROJAS
THE KALMAN FILTER RAUL ROJAS Abstract. This paper provides a getle itroductio to the Kalma filter, a umerical method that ca be used for sesor fusio or for calculatio of trajectories. First, we cosider
More informationTMA4245 Statistics. Corrected 30 May and 4 June Norwegian University of Science and Technology Department of Mathematical Sciences.
Norwegia Uiversity of Sciece ad Techology Departmet of Mathematical Scieces Corrected 3 May ad 4 Jue Solutios TMA445 Statistics Saturday 6 May 9: 3: Problem Sow desity a The probability is.9.5 6x x dx
More informationThermodynamical analysis for a variable generalized Chaplygin gas
Available olie at www.worldscietificews.com WS 66 (7) 49-6 EISS 9-9 hermodyamical aalysis for a variable geeralized haplygi gas Mauel Malaver Departmet of asic Scieces, Maritime iversity of the aribbea,
More informationControl Charts for Mean for Non-Normally Correlated Data
Joural of Moder Applied Statistical Methods Volume 16 Issue 1 Article 5 5-1-017 Cotrol Charts for Mea for No-Normally Correlated Data J. R. Sigh Vikram Uiversity, Ujjai, Idia Ab Latif Dar School of Studies
More informationCastiel, Supernatural, Season 6, Episode 18
13 Differetial Equatios the aswer to your questio ca best be epressed as a series of partial differetial equatios... Castiel, Superatural, Seaso 6, Episode 18 A differetial equatio is a mathematical equatio
More informationNUMERICAL METHOD FOR SINGULARLY PERTURBED DELAY PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS
THERMAL SCIENCE, Year 07, Vol., No. 4, pp. 595-599 595 NUMERICAL METHOD FOR SINGULARLY PERTURBED DELAY PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS by Yula WANG *, Da TIAN, ad Zhiyua LI Departmet of Mathematics,
More information1. pn junction under bias 2. I-Vcharacteristics
Lecture 10 The p Juctio (II) 1 Cotets 1. p juctio uder bias 2. I-Vcharacteristics 2 Key questios Why does the p juctio diode exhibit curret rectificatio? Why does the juctio curret i forward bias icrease
More informationChE 471 Lecture 10 Fall 2005 SAFE OPERATION OF TUBULAR (PFR) ADIABATIC REACTORS
SAFE OPERATION OF TUBULAR (PFR) ADIABATIC REACTORS I a exothermic reactio the temperature will cotiue to rise as oe moves alog a plug flow reactor util all of the limitig reactat is exhausted. Schematically
More informationIntroduction to Astrophysics Tutorial 2: Polytropic Models
Itroductio to Astrophysics Tutorial : Polytropic Models Iair Arcavi 1 Summary of the Equatios of Stellar Structure We have arrived at a set of dieretial equatios which ca be used to describe the structure
More informationThe axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.
5.4 Applicatio of Perturbatio Methods to the Dispersio Model for Tubular Reactors The axial dispersio model for tubular reactors at steady state ca be described by the followig equatios: d c Pe dz z =
More informationA STUDY ON MHD BOUNDARY LAYER FLOW OVER A NONLINEAR STRETCHING SHEET USING IMPLICIT FINITE DIFFERENCE METHOD
IRET: Iteratioal oural of Research i Egieerig ad Techology eissn: 39-63 pissn: 3-7308 A STUDY ON MHD BOUNDARY LAYER FLOW OVER A NONLINEAR STRETCHING SHEET USING IMPLICIT FINITE DIFFERENCE METHOD Satish
More informationFluid Physics 8.292J/12.330J % (1)
Fluid Physics 89J/133J Problem Set 5 Solutios 1 Cosider the flow of a Euler fluid i the x directio give by for y > d U = U y 1 d for y d U + y 1 d for y < This flow does ot vary i x or i z Determie the
More informationA NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS
Jural Karya Asli Loreka Ahli Matematik Vol. No. (010) page 6-9. Jural Karya Asli Loreka Ahli Matematik A NEW CLASS OF -STEP RATIONAL MULTISTEP METHODS 1 Nazeeruddi Yaacob Teh Yua Yig Norma Alias 1 Departmet
More informationNumerical Simulation of Thermomechanical Problems in Applied Mechanics: Application to Solidification Problem
Leoardo Joural of Scieces ISSN 1583-0233 Issue 9, July-December 2006 p. 25-32 Numerical Simulatio of Thermomechaical Problems i Applied Mechaics: Applicatio to Solidificatio Problem Vicet Obiajulu OGWUAGWU
More informationStopping oscillations of a simple harmonic oscillator using an impulse force
It. J. Adv. Appl. Math. ad Mech. 5() (207) 6 (ISSN: 2347-2529) IJAAMM Joural homepage: www.ijaamm.com Iteratioal Joural of Advaces i Applied Mathematics ad Mechaics Stoppig oscillatios of a simple harmoic
More informationWHAT IS THE PROBABILITY FUNCTION FOR LARGE TSUNAMI WAVES? ABSTRACT
WHAT IS THE PROBABILITY FUNCTION FOR LARGE TSUNAMI WAVES? Harold G. Loomis Hoolulu, HI ABSTRACT Most coastal locatios have few if ay records of tsuami wave heights obtaied over various time periods. Still
More information17 Phonons and conduction electrons in solids (Hiroshi Matsuoka)
7 Phoos ad coductio electros i solids Hiroshi Matsuoa I this chapter we will discuss a miimal microscopic model for phoos i a solid ad a miimal microscopic model for coductio electros i a simple metal.
More informationA CONFINEMENT MODEL OF HIGH STRENGTH CONCRETE
3 th World Coferece o Earthquake Egieerig Vacouver, B.C., Caada August -6, 24 Paper No. 873 A CONFINEMENT MODEL OF HIGH STRENGTH CONCRETE Nobutaka NAKAZAWA, Kazuhiko KAWASHIMA 2, Gakuho WATANABE 3, Ju-ichi
More information2.004 Dynamics and Control II Spring 2008
MIT OpeCourseWare http://ocw.mit.edu 2.004 Dyamics ad Cotrol II Sprig 2008 For iformatio about citig these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Massachusetts Istitute of Techology
More informationNumerical Method for Blasius Equation on an infinite Interval
Numerical Method for Blasius Equatio o a ifiite Iterval Alexader I. Zadori Omsk departmet of Sobolev Mathematics Istitute of Siberia Brach of Russia Academy of Scieces, Russia zadori@iitam.omsk.et.ru 1
More informationLOSS-MINIMIZATION CONTROL OF SCALAR- CONTROLLED INDUCTION MOTOR DRIVES
LOSS-MINIMIZATION CONTROL OF SCALAR- CONTROLLED INDUCTION MOTOR DRIVES Hussei Sarha, Rateb Al-Issa, ad Qazem Jaber Departmet of Mechatroics Egieerig, Faculty of Egieerig Techology Al-Balqa Applied Uiversity,
More informationR is a scalar defined as follows:
Math 8. Notes o Dot Product, Cross Product, Plaes, Area, ad Volumes This lecture focuses primarily o the dot product ad its may applicatios, especially i the measuremet of agles ad scalar projectio ad
More informationModeling and Estimation of a Bivariate Pareto Distribution using the Principle of Maximum Entropy
Sri Laka Joural of Applied Statistics, Vol (5-3) Modelig ad Estimatio of a Bivariate Pareto Distributio usig the Priciple of Maximum Etropy Jagathath Krisha K.M. * Ecoomics Research Divisio, CSIR-Cetral
More information3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,
3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [
More informationAppendix: The Laplace Transform
Appedix: The Laplace Trasform The Laplace trasform is a powerful method that ca be used to solve differetial equatio, ad other mathematical problems. Its stregth lies i the fact that it allows the trasformatio
More informationFREE VIBRATION RESPONSE OF A SYSTEM WITH COULOMB DAMPING
Mechaical Vibratios FREE VIBRATION RESPONSE OF A SYSTEM WITH COULOMB DAMPING A commo dampig mechaism occurrig i machies is caused by slidig frictio or dry frictio ad is called Coulomb dampig. Coulomb dampig
More informationNumerical Solution of the First-Order Hyperbolic Partial Differential Equation with Point-Wise Advance
Iteratioal oural of Sciece ad Research (ISR) ISSN (Olie): 39-74 Ide Copericus Value (3): 4 Impact Factor (3): 4438 Numerical Solutio of the First-Order Hyperbolic Partial Differetial Equatio with Poit-Wise
More informationExact Solutions for a Class of Nonlinear Singular Two-Point Boundary Value Problems: The Decomposition Method
Exact Solutios for a Class of Noliear Sigular Two-Poit Boudary Value Problems: The Decompositio Method Abd Elhalim Ebaid Departmet of Mathematics, Faculty of Sciece, Tabuk Uiversity, P O Box 741, Tabuki
More informationFIR Filter Design: Part II
EEL335: Discrete-Time Sigals ad Systems. Itroductio I this set of otes, we cosider how we might go about desigig FIR filters with arbitrary frequecy resposes, through compositio of multiple sigle-peak
More informationProof of Fermat s Last Theorem by Algebra Identities and Linear Algebra
Proof of Fermat s Last Theorem by Algebra Idetities ad Liear Algebra Javad Babaee Ragai Youg Researchers ad Elite Club, Qaemshahr Brach, Islamic Azad Uiversity, Qaemshahr, Ira Departmet of Civil Egieerig,
More informationTEACHER CERTIFICATION STUDY GUIDE
COMPETENCY 1. ALGEBRA SKILL 1.1 1.1a. ALGEBRAIC STRUCTURES Kow why the real ad complex umbers are each a field, ad that particular rigs are ot fields (e.g., itegers, polyomial rigs, matrix rigs) Algebra
More informationPhysics Supplement to my class. Kinetic Theory
Physics Supplemet to my class Leaers should ote that I have used symbols for geometrical figures ad abbreviatios through out the documet. Kietic Theory 1 Most Probable, Mea ad RMS Speed of Gas Molecules
More informationJournal of Computational Physics 149, (1999) Article ID jcph , available online at
Joural of Computatioal Physics 149, 418 422 (1999) Article ID jcph.1998.6131, available olie at http://www.idealibrary.com o NOTE Defiig Wave Amplitude i Characteristic Boudary Coditios Key Words: Euler
More informationThe Space Redundant Robotic Manipulator Chaotic Motion Dynamics Control Algorithm
Sesors & rasducers, Vol. 75, Issue 7, July 24, pp. 27-3 Sesors & rasducers 24 by IFSA Publishig, S. L. http://www.sesorsportal.com he Space Redudat Robotic Maipulator Chaotic Motio Dyamics Cotrol Algorithm
More informationNew Version of the Rayleigh Schrödinger Perturbation Theory: Examples
New Versio of the Rayleigh Schrödiger Perturbatio Theory: Examples MILOŠ KALHOUS, 1 L. SKÁLA, 1 J. ZAMASTIL, 1 J. ČÍŽEK 2 1 Charles Uiversity, Faculty of Mathematics Physics, Ke Karlovu 3, 12116 Prague
More informationECONOMIC OPERATION OF POWER SYSTEMS
ECOOMC OEATO OF OWE SYSTEMS TOUCTO Oe of the earliest applicatios of o-lie cetralized cotrol was to provide a cetral facility, to operate ecoomically, several geeratig plats supplyig the loads of the system.
More informationAlgorithm of Superposition of Boolean Functions Given with Truth Vectors
IJCSI Iteratioal Joural of Computer Sciece Issues, Vol 9, Issue 4, No, July ISSN (Olie: 694-84 wwwijcsiorg 9 Algorithm of Superpositio of Boolea Fuctios Give with Truth Vectors Aatoly Plotikov, Aleader
More informationSubject: Differential Equations & Mathematical Modeling-III
Power Series Solutios of Differetial Equatios about Sigular poits Subject: Differetial Equatios & Mathematical Modelig-III Lesso: Power series solutios of differetial equatios about Sigular poits Lesso
More informationCS322: Network Analysis. Problem Set 2 - Fall 2009
Due October 9 009 i class CS3: Network Aalysis Problem Set - Fall 009 If you have ay questios regardig the problems set, sed a email to the course assistats: simlac@staford.edu ad peleato@staford.edu.
More informationInelastic spherical caps with defects
Ielastic spherical caps with defects JAAN LLLP RNST TUNGL Istitute of athematics Uiversity of Tartu Liivi str Tartu 549 STONIA aalellep@utee Tartu College Talli Uiversity of Techology Puiestee 78 58 Tartu
More informationMechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter
Time Respose & Frequecy Respose d -Order Dyamic System -Pole, Low-Pass, Active Filter R 4 R 7 C 5 e i R 1 C R 3 - + R 6 - + e out Assigmet: Perform a Complete Dyamic System Ivestigatio of the Two-Pole,
More informationA PROCEDURE TO MODIFY THE FREQUENCY AND ENVELOPE CHARACTERISTICS OF EMPIRICAL GREEN'S FUNCTION. Lin LU 1 SUMMARY
A POCEDUE TO MODIFY THE FEQUENCY AND ENVELOPE CHAACTEISTICS OF EMPIICAL GEEN'S FUNCTION Li LU SUMMAY Semi-empirical method, which divides the fault plae of large earthquake ito mets ad uses small groud
More informationResearch Article Health Monitoring for a Structure Using Its Nonstationary Vibration
Advaces i Acoustics ad Vibratio Volume 2, Article ID 69652, 5 pages doi:.55/2/69652 Research Article Health Moitorig for a Structure Usig Its Nostatioary Vibratio Yoshimutsu Hirata, Mikio Tohyama, Mitsuo
More informationQuantum Annealing for Heisenberg Spin Chains
LA-UR # - Quatum Aealig for Heiseberg Spi Chais G.P. Berma, V.N. Gorshkov,, ad V.I.Tsifriovich Theoretical Divisio, Los Alamos Natioal Laboratory, Los Alamos, NM Istitute of Physics, Natioal Academy of
More informationThe Wave Function and Quantum Reality
The Wave Fuctio ad Quatum Reality Sha Gao Uit for History ad Philosophy of Sciece & Cetre for Time, SOPHI Uiversity of Sydey, Sydey, NSW 006, Australia Abstract. We ivestigate the meaig of the wave fuctio
More informationDISTRIBUTION LAW Okunev I.V.
1 DISTRIBUTION LAW Okuev I.V. Distributio law belogs to a umber of the most complicated theoretical laws of mathematics. But it is also a very importat practical law. Nothig ca help uderstad complicated
More informationUNIFORM FLOW. U x. U t
UNIFORM FLOW if : 1) there are o appreciable variatios i the chael geometry (width, slope, roughess/grai size), for a certai legth of a river reach ) flow discharge does ot vary the, UNIFORM FLOW coditios
More information577. Estimation of surface roughness using high frequency vibrations
577. Estimatio of surface roughess usig high frequecy vibratios V. Augutis, M. Sauoris, Kauas Uiversity of Techology Electroics ad Measuremets Systems Departmet Studetu str. 5-443, LT-5368 Kauas, Lithuaia
More informationStability Analysis and Bifurcation Control of Hysteresis Current Controlled Ćuk Converter Using Filippov s Method
Stability Aalysis ad Bifurcatio otrol of Hysteresis urret otrolled Ću overter Usig Filippov s Method I. Daho*, D. Giaouris * (Member IE, B. Zahawi*(Member IE, V. Picer*(Member IE ad S. Baerjee** *School
More informationCUMULATIVE DAMAGE ESTIMATION USING WAVELET TRANSFORM OF STRUCTURAL RESPONSE
CUMULATIVE DAMAGE ESTIMATION USING WAVELET TRANSFORM OF STRUCTURAL RESPONSE Ryutaro SEGAWA 1, Shizuo YAMAMOTO, Akira SONE 3 Ad Arata MASUDA 4 SUMMARY Durig a strog earthquake, the respose of a structure
More informationMeasurement uncertainty of the sound absorption
Measuremet ucertaity of the soud absorptio coefficiet Aa Izewska Buildig Research Istitute, Filtrowa Str., 00-6 Warsaw, Polad a.izewska@itb.pl 6887 The stadard ISO/IEC 705:005 o the competece of testig
More informationHolistic Approach to the Periodic System of Elements
Holistic Approach to the Periodic System of Elemets N.N.Truov * D.I.Medeleyev Istitute for Metrology Russia, St.Peterburg. 190005 Moskovsky pr. 19 (Dated: February 20, 2009) Abstract: For studyig the objectivity
More informationSPEC/4/PHYSI/SPM/ENG/TZ0/XX PHYSICS PAPER 1 SPECIMEN PAPER. 45 minutes INSTRUCTIONS TO CANDIDATES
SPEC/4/PHYSI/SPM/ENG/TZ0/XX PHYSICS STANDARD LEVEL PAPER 1 SPECIMEN PAPER 45 miutes INSTRUCTIONS TO CANDIDATES Do ot ope this examiatio paper util istructed to do so. Aswer all the questios. For each questio,
More informationA Simplified Binet Formula for k-generalized Fibonacci Numbers
A Simplified Biet Formula for k-geeralized Fiboacci Numbers Gregory P. B. Dresde Departmet of Mathematics Washigto ad Lee Uiversity Lexigto, VA 440 dresdeg@wlu.edu Zhaohui Du Shaghai, Chia zhao.hui.du@gmail.com
More informationUniform Strict Practical Stability Criteria for Impulsive Functional Differential Equations
Global Joural of Sciece Frotier Research Mathematics ad Decisio Scieces Volume 3 Issue Versio 0 Year 03 Type : Double Blid Peer Reviewed Iteratioal Research Joural Publisher: Global Jourals Ic (USA Olie
More informationFinally, we show how to determine the moments of an impulse response based on the example of the dispersion model.
5.3 Determiatio of Momets Fially, we show how to determie the momets of a impulse respose based o the example of the dispersio model. For the dispersio model we have that E θ (θ ) curve is give by eq (4).
More informationChimica Inorganica 3
himica Iorgaica Irreducible Represetatios ad haracter Tables Rather tha usig geometrical operatios, it is ofte much more coveiet to employ a ew set of group elemets which are matrices ad to make the rule
More informationTR/46 OCTOBER THE ZEROS OF PARTIAL SUMS OF A MACLAURIN EXPANSION A. TALBOT
TR/46 OCTOBER 974 THE ZEROS OF PARTIAL SUMS OF A MACLAURIN EXPANSION by A. TALBOT .. Itroductio. A problem i approximatio theory o which I have recetly worked [] required for its solutio a proof that the
More informationMATH 1080: Calculus of One Variable II Fall 2017 Textbook: Single Variable Calculus: Early Transcendentals, 7e, by James Stewart.
MATH 1080: Calculus of Oe Variable II Fall 2017 Textbook: Sigle Variable Calculus: Early Trascedetals, 7e, by James Stewart Uit 3 Skill Set Importat: Studets should expect test questios that require a
More information