Systems of Linear First Order Ordinary Differential Equations Example Problems
|
|
- Imogene Hunter
- 5 years ago
- Views:
Transcription
1 Systes of Linear First Order Ordinary Differential Equations Eaple Probles David Keffer Departent of Cheial Engineering University of Tennessee Knoville, TN 79 Last Updated: Septeber 4, Eaple. Transient Behavior of a Bath Reator Eaple. Transient Behavior of a Continuous Stirred-Tan Reator 6 Eaple. Transient Vibrational Behavior of Carbon Dioide
2 D. Keffer, ChE 55, Dept. of Cheial Engineering, University of Tennessee, Septeber, Eaple. Transient Behavior of a Bath Reator Consider that you have a three-oponent reative iture in a bath reator, all undergoing reversible reations, as pitured below: In this piture, the s are onentrations of the three speies and the s are rate onstants. n eaple of this syste is the ineti equilibriu between para-, eta-, and ortho-ylene. Now suppose we want to now what the onentration is as a funtion of tie. We an write the ass balanes for eah oponent. There are no in and out ters (the reator is a bath reator. There is only the auulation ter and the reation ters. lso, assue eah reation is first order in onentration. d d d (. We an gather lie ters and rearrange the right hand side: d d d ( + ( ( (. and we hange this syste of equations into atri & vetor for:
3 D. Keffer, ChE 55, Dept. of Cheial Engineering, University of Tennessee, Septeber, d X (. where ( + X ( + ( + (.4 The atri X is singular. It has a deterinant of zero and a ran of. Therefore, it has one zero value eigenvalue. Nevertheless, the atri has three distint real eigenvalues, λ, λ, λ. Corresponding to eah of these eigenvalues is a real, distint eigenvetor, w,,w,, w,. The eigenvetors of X are nown as the eigenrows of X: w r,,w r,, wr,. The solution is then: ( t ( wr,i o ( w w w,i ep[ λi( t to ] r,i,i n i (.5 In order to obtain nuerial values for this proble, we would have to first obtain the eigenvalues, eigenvetors, and eigenrows. This an be done nuerially, using routines disussed in the setion on nuerial ethods for solving systes of linear algebrai equations. We provide a Matlab ode to solve this proble on the following page.
4 D. Keffer, ChE 55, Dept. of Cheial Engineering, University of Tennessee, Septeber, Matlab ode to solve Eaple. % % Solution to the Bath Reator Proble % % David Keffer % Departent of Cheial Engineering % University of Tennessee % Knoville Tennessee % Septeber, % % input.5;.5;.;.5;.;.5; [ (--,, ;, (--, ;,, (--]; yo [./.;./.;./.]; to.; tf.; n a(size(; % eigenvalues and eigenvetors [wol,labda] eig(; wolinv inv(wol; % set up disretized solution grid npoints ; (tf-to/npoints; for i ::npoints+ tp(i (i-* +to; % alulate solution eplabda zeros(n,n; yp zeros(n,npoints; for i ::npoints+ for j :n eplabda(j,j ep(labda(j,j*(tp(i - to ; yp(:,i wol*eplabda*wolinv*yo; %plot for j ::n if (j plot (tp,yp(j,:,'g-', label( 't', ylabel ( 'y' elseif (j plot (tp,yp(j,:,'r-', label( 't', ylabel ( 'y' elseif (j plot (tp,yp(j,:,'b-', label( 't', ylabel ( 'y' hold on hold off
5 D. Keffer, ChE 55, Dept. of Cheial Engineering, University of Tennessee, Septeber, Plot of Solution to Eaple : Transient Behavior in a Bath Reator ole fration tie (s Coents: t all ties, the su of the ole frations is unity. The initial ondition is y o [/ / /]. The steady state equilibriu appears to be about y [ ]. 4
6 D. Keffer, ChE 55, Dept. of Cheial Engineering, University of Tennessee, Septeber,.9.8 ole fration tie (s Coents: t all ties, the su of the ole frations is unity. The initial ondition is y o [ ]. The steady state equilibriu appears to be about y [ ]. This is the sae steady state oposition as was obtained for the other initial ondition, plotted on the previous page. The equilibriu oposition is indepent of the initial oposition, in this ase. This is generally true, eept where ultiple steady states eist. 5
7 D. Keffer, ChE 55, Dept. of Cheial Engineering, University of Tennessee, Septeber, Eaple. Cheial Reation Equilibria with addition/deletion of oponents b( We an add onstant flowrates to our syste by adding a onstant ter to the ODEs. d d d ( + ( ( F + F + F (4. The flowrates an either be feed or effluent, deping upon the sign. To eep the proble siple, let's enfore onservation of volue by stipulating that i F i (4. and we hange this syste of equations into atri & vetor for: d (4. X + F We now the solution to the nonhoogeneous equation is ( t W ep Λ( t [ ] ep[ Λ( t ] W u(t + u(t o o o nh (4.4 The only question is the for of u. u [ Λ( ] W b( d ep (. [ ] Sine, for this proble, b is a onstant, all the funtionality lies in ep Λ( u ( t ep Λ( [ ] dw b t (4.5 6
8 D. Keffer, ChE 55, Dept. of Cheial Engineering, University of Tennessee, Septeber, u ( t Λ ep Λ( t [ ] W b (4.6 Plugging ba into the solution we find nh [ ] ( t W ep Λ( t ep Λ [ ( ] Λ to W o + Λ ep Λ( t o ep [ ] W b [ (] Λ t W b (4.7 opliation arises beause one of the eigenvalues is zero. Thus there is no -depeny in the eponential beause there is no eponential. Then we would have a ter of the for u t t (4.8 ( t ep[ ] dw b dw b [] t W b tw b so that we would have u ( t UW b (4.9 ep λ U [ λ( t ] ep [ λ( t ] λ for the ase where we had three eigenvalues, the first two of whih are non-zero. t (4. We provide a Matlab ode to solve this proble on the following page. 7
9 D. Keffer, ChE 55, Dept. of Cheial Engineering, University of Tennessee, Septeber, Matlab ode to solve Eaple. % % Solution to the Continuous Stirred Tan Reator Proble % % David Keffer % Departent of Cheial Engineering % University of Tennessee % Knoville Tennessee % Septeber, % tol.e-4; % input.5;.5;.;.5;.;.5; [ (--,, ;, (--, ;,, (--]; n a(size(; yo [./.;./.;./.]; %yo [.;.;.]; % onstant ters in ODE b [-.; -.;.]; to.; tf.; % eigenvalues and eigenvetors [wol,labda] eig(; wolinv inv(wol; % alulate ep (-labda*to atri eplabdano zeros(n,n; for j ::di eplabdano(j,j ep(-labda(j,j*to; % alulate u(to uo zeros(n,n; for j ::n if (abs(labda(j,j > tol uo(j,j -eplabdano(j,j/labda(j,j; else uo(j,j to; % set up disretized solution grid nintervals ; npoints nintervals + ; (tf-to/nintervals; for i ::npoints tp(i (i-* +to; % alulate solution yp zeros(n,npoints; for i ::npoints eplabdan zeros(n,n; eplabda zeros(n,n; for j ::di eplabdan(j,j ep(-labda(j,j*tp(i; 8
10 D. Keffer, ChE 55, Dept. of Cheial Engineering, University of Tennessee, Septeber, eplabda(j,j ep(labda(j,j*tp(i; % alulate u(t u zeros(n,n; for j ::n if (abs(labda(j,j > tol u(j,j -eplabdan(j,j/labda(j,j; else u(j,j tp(i; ter eplabdano*wolinv*yo; ter (u-uo*wolinv*b; yp(:,i wol*eplabda*(ter + ter; %plot for j ::n if (j plot (tp,yp(j,:,'g-', label( 't', ylabel ( 'y' elseif (j plot (tp,yp(j,:,'r-', label( 't', ylabel ( 'y' elseif (j plot (tp,yp(j,:,'b-', label( 't', ylabel ( 'y' hold on hold off label('tie (s'; ylabel('ole fration'; leg(' ',' ', ' '; 9
11 D. Keffer, ChE 55, Dept. of Cheial Engineering, University of Tennessee, Septeber, ole fration tie (s Coents: t all ties, the su of the ole frations is unity. The initial ondition is y o [/ / /]. The steady state equilibriu appears to be about y [..6.65]. The flowrates for eah oponent are b[-...].
12 D. Keffer, ChE 55, Dept. of Cheial Engineering, University of Tennessee, Septeber, ole fration tie (s Coents: t all ties, the su of the ole frations is unity. The initial ondition is y o [/ / ]. The steady state equilibriu appears to be about y [..6.65]. The flowrates for eah oponent are b[-...]. This is the sae steady state oposition as was obtained for the other initial ondition, plotted on the previous page. The equilibriu oposition is indepent of the initial oposition, in this ase. This is generally true, eept where ultiple steady states eist.
13 D. Keffer, ChE 55, Dept. of Cheial Engineering, University of Tennessee, Septeber, Eaple. Transient Vibrational Behavior of Carbon Dioide Consider that we want to investigate the vibrational properties of arbon dioide, CO. Our odel of the oleule loos lie this: spring ( spring ( O C O,,, We odel the interation between oleules as Hooian springs. For a Hooian spring, the potential energy, U, is U ( (5. and the fore, F, is F ( (5. where is the spring onstant (units of g/s, is the equilibriu displaeent, and is the atual displaeent. When both s of the spring are obile we an write, for the spring that onnets ass and : U (5. (( and we an write, for the spring that onnets ass and : U (5.4 ((
14 D. Keffer, ChE 55, Dept. of Cheial Engineering, University of Tennessee, Septeber, The fores are then: (( U F (( (( U F + (5.5 (( U F With only a slight sleight of hand, we redefine our variables to be: + (5.6 This eliinates the equilibriu bond distanes fro the alulation. So, with this definition, (5.7 at equilibriu. This does not affet our equations of otion beause the derivatives of our variables before and after the transforation of the variables are the sae. Our fores beoe: ( U F ( ( U F + (5.8 ( U F We an write Newton s equations of otion for the three oleules: ( ( ( ( O C O F a F a F a + (5.9 We an generalized this to a non-syetri linear tri-atoi oleule by writing:
15 D. Keffer, ChE 55, Dept. of Cheial Engineering, University of Tennessee, Septeber, a a a F F F ( ( + ( ( (5. Knowing that the aeleration is the seond derivative of the position, we an rewrite the above equations in atri for as (first divide both side of all of the equations by the asses d (5. where (5. Solving this syste of seond order linear differential equations yields the integrated equations of otion for arbon dioide. We an onvert the three seond order ODEs into si first order ODEs as shown below. Tae ato nuber one. whih we write as ( a F (5. d ( (5.4 Now we an rewrite that seond order ODE as two first order ODEs. d 4 (5.5 4
16 D. Keffer, ChE 55, Dept. of Cheial Engineering, University of Tennessee, Septeber, 5 ( 4 d (5.6 We an write analogous equation for the other two atos, to oe with ODEs for variables 5 and 6. Of ourse, variables,, and are still the positions of the atos. Variables 4, 5, and 6 are now the veloities of the atos. ll of these equations are hoogeneous. Our atri beoes a 66 atri: (5.7 where the solution vetor,, is now defined as: & & & (5.8 Now, learly, we an solve this proble eatly as we solve any syste of hoogeneous linear first order ODEs. The solution will be given as ( ( [ ] o o h W t t ep W t Λ (.4 We provide a Matlab ode to solve this proble on the following page. The only quantitative differene between this ode and the ode for eaple, is that this proble has iaginary eigenvalues, whih eventually yield real solutions via the Euler Identities. Near the of the ode, we eliinate any residual round-off errors of the iaginary oponent of the solution.
17 D. Keffer, ChE 55, Dept. of Cheial Engineering, University of Tennessee, Septeber, Matlab ode to solve Eaple. % % Solution to vibrational odes of CO Proble % % David Keffer % Departent of Cheial Engineering % University of Tennessee % Knoville Tennessee % Septeber, % lear all; lose all; % tol.e-4; % input ass( 6; ass( ; ass( 6; 48; 48; [ -/ass( /ass( /ass( (--/ass( /ass( /ass( -/ass( ]; n a(size(; yo [.;.; -.5;.;.;.]; % onstant ters in ODE b zeros(n,; to.; tf.; % eigenvalues and eigenvetors [wol,labda] eig(; wolinv inv(wol; % alulate ep (-labda*to atri eplabdano zeros(n,n; for j ::n eplabdano(j,j ep(-labda(j,j*to; 6
18 D. Keffer, ChE 55, Dept. of Cheial Engineering, University of Tennessee, Septeber, % alulate u(to uo zeros(n,n; for j ::n if (abs(labda(j,j > tol uo(j,j -eplabdano(j,j/labda(j,j; else uo(j,j to; % set up disretized solution grid nintervals ; npoints nintervals + ; (tf-to/nintervals; for i ::npoints tp(i (i-* +to; % alulate solution yp zeros(n,npoints; for i ::npoints eplabdan zeros(n,n; eplabda zeros(n,n; for j ::n eplabdan(j,j ep(-labda(j,j*tp(i; eplabda(j,j ep(labda(j,j*tp(i; % alulate u(t u zeros(n,n; for j ::n if (abs(labda(j,j > tol u(j,j -eplabdan(j,j/labda(j,j; else u(j,j tp(i; ter eplabdano*wolinv*yo; ter (u-uo*wolinv*b; yp(:,i wol*eplabda*(ter + ter; % reove residual iaginary oponents due to round-off for i ::npoints for j ::n 7
19 D. Keffer, ChE 55, Dept. of Cheial Engineering, University of Tennessee, Septeber, if (iag(yp(j,i < tol yp(j,i real( yp(j,i ; %plot for j ::n if (j plot (tp,yp(j,:,'g-' elseif (j plot (tp,yp(j,:,'r-' elseif (j plot (tp,yp(j,:,'b-' elseif (j4 plot (tp,yp(j,:,'g:' elseif (j5 plot (tp,yp(j,:,'r:' elseif (j6 plot (tp,yp(j,:,'b:' hold on hold off label('tie '; ylabel('position & veloity '; leg('pos O','pos C', 'pos O','vel O','vel C', 'vel O'; 8
20 D. Keffer, ChE 55, Dept. of Cheial Engineering, University of Tennessee, Septeber,.5.5 position & veloity pos O pos C pos O vel O vel C vel O tie Coents: The initial ondition is y o [.5 ]. Clearly periodi behavior is observed. lso observe visually that the veloities do indeed orrespond to the tie derivatives of the positions. There is no net enter of ass otion, sine in the initial ondition there was no enter of ass otion, and onservation of oentu is a onsequene of Newton s equations of otion. 9
21 D. Keffer, ChE 55, Dept. of Cheial Engineering, University of Tennessee, Septeber, 5 4 pos O pos C pos O vel O vel C vel O position & veloity tie Coents: The initial ondition is y o [ ]. Clearly periodi behavior is observed. lso observe visually that the veloities do indeed orrespond to the tie derivatives of the positions. There is a net enter of ass otion, sine in the initial ondition there was a enter of ass otion, and onservation of oentu is a onsequene of Newton s equations of otion.
Included in this hand-out are five examples of problems requiring the solution of a system of linear algebraic equations.
he Lecture Notes, Dept. of heical Engineering, Univ. of TN, Knoville - D. Keffer, 5/9/98 (updated /) Eaple pplications of systes of linear equations Included in this hand-out are five eaples of probles
More informationReference. R. K. Herz,
Identifiation of CVD kinetis by the ethod of Koiyaa, et al. Coparison to 1D odel (2012) filenae: CVD_Koiyaa_1D_odel Koiyaa, et al. (1999) disussed ethods to identify the iportant steps in a CVD reation
More informationLecture 24: Spinodal Decomposition: Part 3: kinetics of the
Leture 4: Spinodal Deoposition: Part 3: kinetis of the oposition flutuation Today s topis Diffusion kinetis of spinodal deoposition in ters of the onentration (oposition) flutuation as a funtion of tie:
More informationWorked Solutions to Problems
rd International Cheistry Olypiad Preparatory Probles Wored Solutions to Probles. Water A. Phase diagra a. he three phases of water oeist in equilibriu at a unique teperature and pressure (alled the triple
More informationAnalytical Study of Stability of Systems of ODEs
D. Keffer, ChE 55,Universit of Tennessee, Ma, 999 Analtial Stud of Stabilit of Sstems of ODEs David Keffer Department of Chemial Engineering Universit of Tennessee, Knoxville date begun: September 999
More informationGreen s Function for Potential Field Extrapolation
Green s Funtion for Potential Field Extrapolation. Soe Preliinaries on the Potential Magneti Field By definition, a potential agneti field is one for whih the eletri urrent density vanishes. That is, J
More informationx(t) y(t) c c F(t) F(t) EN40: Dynamics and Vibrations Homework 6: Forced Vibrations Due Friday April 5, 2018
EN40: Dynais and Vibrations Hoewor 6: Fored Vibrations Due Friday April 5, 2018 Shool of Engineering Brown University 1. The vibration isolation syste shown in the figure has =20g, = 19.8 N / = 1.259 Ns
More informationDerivation of Non-Einsteinian Relativistic Equations from Momentum Conservation Law
Asian Journal of Applied Siene and Engineering, Volue, No 1/13 ISSN 35-915X(p); 37-9584(e) Derivation of Non-Einsteinian Relativisti Equations fro Moentu Conservation Law M.O.G. Talukder Varendra University,
More information1. Which two values of temperature are equivalent to the nearest degree when measured on the Kelvin and on the
. Whih two values of teperature are equivalent to the nearest degree when easured on the Kelvin and on the Celsius sales of teperature? Kelvin sale Celsius sale A. 40 33 B. 273 00 C. 33 40 D. 373 0 2.
More informationCHAPTER 3 PROBLEMS. δ = where A is the cross-sectional area, and E is the modulus of elasticity.
CHPTER 3 PROLEMS d SPRING-MSS-DMPER PPLICTIONS Proble 3.1 The buoy shown in figure P3.1 has a irular ross-setion with diaeter d and has length L. Most of the weight of the buoy, w, is onentrated in the
More informationModeling of vibration systems
Modeling of vibration systes Atual syste Mae design deision Choose physial paraeters, hange or augent syste if neessary Physial odeling Mae siple approiations based on engineering judgeent Physial odel
More informationForce and dynamics with a spring, analytic approach
Force and dynaics with a spring, analytic approach It ay strie you as strange that the first force we will discuss will be that of a spring. It is not one of the four Universal forces and we don t use
More informationLecture 23: Spinodal Decomposition: Part 2: regarding free energy. change and interdiffusion coefficient inside the spinodal
Leture 3: Spinodal eoposition: Part : regarding free energy hange and interdiffusion oeffiient inside the spinodal Today s topis ontinue to understand the basi kinetis of spinodal deoposition. Within the
More informationInternational Journal of Thermodynamics, Vol. 18, No. 1, P (2015). Sergey G.
International Journal of Therodynais Vol. 8 No. P. 3-4 (5). http://dx.doi.org/.554/ijot.5343 Four-diensional equation of otion for visous opressible and harged fluid with regard to the aeleration field
More information(Newton s 2 nd Law for linear motion)
PHYSICS 3 Final Exaination ( Deeber Tie liit 3 hours Answer all 6 questions You and an assistant are holding the (opposite ends of a long plank when oops! the butterfingered assistant drops his end If
More informationKinematics of Elastic Neutron Scattering
.05 Reator Physis - Part Fourteen Kineatis of Elasti Neutron Sattering. Multi-Group Theory: The next ethod that we will study for reator analysis and design is ulti-group theory. This approah entails dividing
More informationME357 Problem Set The wheel is a thin homogeneous disk that rolls without slip. sin. The wall moves with a specified motion x t. sin..
ME357 Proble Set 3 Derive the equation(s) of otion for the systes shown using Newton s Method. For ultiple degree of freedo systes put you answer in atri for. Unless otherwise speified the degrees of freedo
More informationDynamics of Structures. Giacomo Boffi. Definitions. Dynamics of Structures. Giacomo Boffi. Introduction. Characteristics of a Dynamical Problem
An to Dipartiento di Ingegneria Civile e Abientale, Politenio di Milano Part I Marh 1, 014 Definitions Definitions Let s start with soe definitions Dynais the branh of ehanis onerned with the effets of
More informationChapter 3. Problem Solutions
Capter. Proble Solutions. A poton and a partile ave te sae wavelengt. Can anyting be said about ow teir linear oenta opare? About ow te poton's energy opares wit te partile's total energy? About ow te
More informationDoppler Effect (Text 1.3)
Doppler Effet (et 1.3) Consider a light soure as a soure sending out a tik eery 1/ν and these tiks are traeling forward with speed. tik tik tik tik Doppler Effet (et 1.3) Case 1. Obserer oing transersely.
More informationThe Seesaw Mechanism
The Seesaw ehanis By obert. Klauber www.quantufieldtheory.info 1 Bakground It ay see unusual to have suh low values for asses of neutrinos, when all other partiles like eletrons, quarks, et are uh heavier,
More informationSimple Harmonic Motion of Spring
Nae P Physics Date iple Haronic Motion and prings Hooean pring W x U ( x iple Haronic Motion of pring. What are the two criteria for siple haronic otion? - Only restoring forces cause siple haronic otion.
More informationAdvanced Computational Fluid Dynamics AA215A Lecture 4
Advaned Computational Fluid Dynamis AA5A Leture 4 Antony Jameson Winter Quarter,, Stanford, CA Abstrat Leture 4 overs analysis of the equations of gas dynamis Contents Analysis of the equations of gas
More informationThe Gravitation As An Electric Effect
The Gravitation As An Eletri Effet Hans-Jörg Hoheker Donaustr 30519 Hannover e-ail: johoer@yahoode Web-Site: http://wwwhohekereu Abstrat: The eletri fores are iensely great in oparison with the gravitational
More informationLecture 4 Normal Modes
Lecture 4 Noral Modes Coupled driven oscillators Double pendulu The daped driven pendulu = g/l +k y+fcost y = y gy/l k y d dt + d dt + g + k l k k d dt + d dt + g + k l y = F 0 Re eit y =Re X Y eit CF
More informationand ζ in 1.1)? 1.2 What is the value of the magnification factor M for system A, (with force frequency ω = ωn
EN40: Dynais and Vibrations Hoework 6: Fored Vibrations, Rigid Body Kineatis Due Friday April 7, 017 Shool of Engineering Brown University 1. Syste A in the figure is ritially daped. The aplitude of the
More informationThe nuclear fusion reaction rate based on relativistic equilibrium velocity distribution
he nulear fusion reation rate based on relativisti equilibriu veloity distribution Jian-Miin Liu* Departent of Physis, Nanjing University Nanjing, he People's Republi of China *On leave. E-ail address:
More informationMAC Calculus II Summer All you need to know on partial fractions and more
MC -75-Calulus II Summer 00 ll you need to know on partial frations and more What are partial frations? following forms:.... where, α are onstants. Partial frations are frations of one of the + α, ( +
More informationOptimization of the CBSMAP Queueing Model
July 3-5 23 London UK Optiization of the CBSMAP Queueing Model Kondrashova EV Kashtanov VA Abstrat The present paper is devoted to the researh of ontrolled queueing odels at ontrol of CBSMAP-flow Controlled
More informationPhysics 221B: Solution to HW # 6. 1) Born-Oppenheimer for Coupled Harmonic Oscillators
Physics B: Solution to HW # 6 ) Born-Oppenheier for Coupled Haronic Oscillators This proble is eant to convince you of the validity of the Born-Oppenheier BO) Approxiation through a toy odel of coupled
More informationSimple Harmonic Motion
Siple Haronic Motion Physics Enhanceent Prograe for Gifted Students The Hong Kong Acadey for Gifted Education and Departent of Physics, HKBU Departent of Physics Siple haronic otion In echanical physics,
More informationElectromagnetic Waves
Eletroagneti Waves Physis 6C Eletroagneti (EM) waves an be produed by atoi transitions (ore on this later), or by an alternating urrent in a wire. As the harges in the wire osillate bak and forth, the
More informationCongruences and Modular Arithmetic
Congruenes and Modular Aritheti 6-17-2016 a is ongruent to b od n eans that n a b. Notation: a = b (od n). Congruene od n is an equivalene relation. Hene, ongruenes have any of the sae properties as ordinary
More informationPAST YEAR EXAM PAPER SOLUTION SEMESTER II, ACADEMIC YEAR PAP261 Introduction to Lasers (Solution by: Phann Sophearin)
SPMS Physis Club PAST YEAR EXAM PAPER SOLUTION SEMESTER II, ACADEMIC YEAR 009-0 PAP6 Introdution to Lasers (Solution by: Phann Sophearin) Question (a) - Three basi eleents required: avity, gain, and pup.
More informationMath 220A - Fall 2002 Homework 8 Solutions
Math A - Fall Homework 8 Solutions 1. Consider u tt u = x R 3, t > u(x, ) = φ(x) u t (x, ) = ψ(x). Suppose φ, ψ are supported in the annular region a < x < b. (a) Find the time T 1 > suh that u(x, t) is
More informationEFFECTIVE MODAL MASS & MODAL PARTICIPATION FACTORS Revision I
EFFECTIVE MODA MASS & MODA PARTICIPATION FACTORS Revision I B To Irvine Eail: to@vibrationdata.co Deceber, 5 Introduction The effective odal ass provides a ethod for judging the significance of a vibration
More informationTUTORIAL 1 SIMPLE HARMONIC MOTION. Instructor: Kazumi Tolich
TUTORIAL 1 SIMPLE HARMONIC MOTION Instructor: Kazui Tolich About tutorials 2 Tutorials are conceptual exercises that should be worked on in groups. Each slide will consist of a series of questions that
More informationWhere as discussed previously we interpret solutions to this partial differential equation in the weak sense: b
Consider the pure initial value problem for a homogeneous system of onservation laws with no soure terms in one spae dimension: Where as disussed previously we interpret solutions to this partial differential
More informationSHM stuff the story continues
SHM stuff the story continues Siple haronic Motion && + ω solution A cos t ( ω + α ) Siple haronic Motion + viscous daping b & + ω & + Viscous daping force A e b t Viscous daped aplitude Viscous daped
More informationNumerical Studies of Counterflow Turbulence
Nonae anusript No. will be inserted by the editor Nuerial Studies of Counterflow Turbulene Veloity Distribution of Vorties Hiroyuki Adahi Makoto Tsubota Reeived: date Aepted: date Abstrat We perfored the
More informationOscillations: Review (Chapter 12)
Oscillations: Review (Chapter 1) Oscillations: otions that are periodic in tie (i.e. repetitive) o Swinging object (pendulu) o Vibrating object (spring, guitar string, etc.) o Part of ediu (i.e. string,
More informationOverview. Review Multidimensional PDEs. Review Diffusion Solutions. Review Separation of Variables. More diffusion solutions February 2, 2009
More diffusion solutions February, 9 More Diffusion Equation Solutions arry Caretto Mechanical Engineering 5B Seinar in Engineering Analysis February, 9 Overview Review last class Separation of variables
More informationENGI 3424 Engineering Mathematics Problem Set 1 Solutions (Sections 1.1 and 1.2)
ENGI 344 Engineering Matheatics Proble Set 1 Solutions (Sections 1.1 and 1.) 1. Find the general solution of the ordinary differential equation y 0 This ODE is not linear (due to the product y ). However,
More informationVelocity Addition in Space/Time David Barwacz 4/23/
Veloity Addition in Spae/Time 003 David arwaz 4/3/003 daveb@triton.net http://members.triton.net/daveb Abstrat Using the spae/time geometry developed in the previous paper ( Non-orthogonal Spae- Time geometry,
More informationTHE ESSENCE OF QUANTUM MECHANICS
THE ESSENCE OF QUANTUM MECHANICS Capter belongs to te "Teory of Spae" written by Dariusz Stanisław Sobolewski. Http: www.tsengines.o ttp: www.teoryofspae.info E-ail: info@tsengines.o All rigts resered.
More informationPY241 Solutions Set 9 (Dated: November 7, 2002)
PY241 Solutions Set 9 (Dated: Noveber 7, 2002) 9-9 At what displaceent of an object undergoing siple haronic otion is the agnitude greatest for the... (a) velocity? The velocity is greatest at x = 0, the
More informationIntro Wb3303 Kineto-statics and Dynamics
Wb333 Kineto-statis and Dynais 8 ineto-statis and Dynais...8. 8. Introdution...8. 8. Virtual wor and equilibriu equations (statis)...8.3 8.. General priniple...8.3 8.. Driving fores...8.5 8..3 Support
More informationRELATIVISTIC GRAVITY AND THE ORIGIN OF INERTIA AND INERTIAL MASS
RELTIVISTIC GRVITY ND THE ORIGIN OF INERTI ND INERTIL MSS K Tsarouhas To ite this version: K Tsarouhas. RELTIVISTIC GRVITY ND THE ORIGIN OF INERTI ND INERTIL MSS. 07. HL Id: hal-047498 https://hal.arhives-ouvertes.fr/hal-047498
More informationExperiment 03: Work and Energy
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physis Department Physis 8.01 Purpose of the Experiment: Experiment 03: Work and Energy In this experiment you allow a art to roll down an inlined ramp and run into
More informationV. Interacting Particles
V. Interating Partiles V.A The Cumulant Expansion The examples studied in the previous setion involve non-interating partiles. It is preisely the lak of interations that renders these problems exatly solvable.
More informationPHYS 1443 Section 003 Lecture #21 Wednesday, Nov. 19, 2003 Dr. Mystery Lecturer
PHYS 443 Section 003 Lecture # Wednesday, Nov. 9, 003 Dr. Mystery Lecturer. Fluid Dyanics : Flow rate and Continuity Equation. Bernoulli s Equation 3. Siple Haronic Motion 4. Siple Bloc-Spring Syste 5.
More informationHORIZONTAL MOTION WITH RESISTANCE
DOING PHYSICS WITH MATLAB MECHANICS HORIZONTAL MOTION WITH RESISTANCE Ian Cooper School of Physics, Uniersity of Sydney ian.cooper@sydney.edu.au DOWNLOAD DIRECTORY FOR MATLAB SCRIPTS ec_fr_b. This script
More informationMODULE 3: MASS TRANSFER COEFFICIENTS
MOULE 3: MASS TRASFER COEFFICIETS LECTURE O. 6 3.5. Penetration theory Most of the industrial proesses of mass transfer is unsteady state proess. In suh ases, the ontat time between phases is too short
More information1. Simplify this circuit to find the total power absorbed by all resistors. R1 R3 47
Hoework 1. iplify this iruit to find the total power asored y all resistors. R1 R 6 R4 0 47 R 47 Ans: Use series and parallel oinations to siplify iruit. R1R= 47 47 =94Ω and this resistane is in parallel
More information23.1 Tuning controllers, in the large view Quoting from Section 16.7:
Lesson 23. Tuning a real ontroller - modeling, proess identifiation, fine tuning 23.0 Context We have learned to view proesses as dynami systems, taking are to identify their input, intermediate, and output
More informationECE-320 Linear Control Systems. Winter 2013, Exam 1. No calculators or computers allowed, you may leave your answers as fractions.
ECE-320 Linear Control Systems Winter 2013, Exam 1 No alulators or omputers allowed, you may leave your answers as frations. All problems are worth 3 points unless noted otherwise. Total /100 1 Problems
More informationChapter 1: Basics of Vibrations for Simple Mechanical Systems
Chapter 1: Basics of Vibrations for Siple Mechanical Systes Introduction: The fundaentals of Sound and Vibrations are part of the broader field of echanics, with strong connections to classical echanics,
More informationChameleon mechanism. Lecture 2
Chaeleon ehanis Leture Cosi aeleration Many independent data sets indiate that the expansion of the Universe is aelerating Siilar to preise tests of GR? Dark energy v Dark gravity Standard odel based on
More informationFurther refutation of the de Broglie Einstein theory in the case of general Compton scattering
Further refutation of the de Broglie Einstein theory 7 Journal of Foundations of Physis and Cheistry, 0, vol () 7 37 Further refutation of the de Broglie Einstein theory in the ase of general Copton sattering
More informationPhysics 2107 Oscillations using Springs Experiment 2
PY07 Oscillations using Springs Experient Physics 07 Oscillations using Springs Experient Prelab Read the following bacground/setup and ensure you are failiar with the concepts and theory required for
More informationShear Force and Bending Moment
Shear Fore and Bending oent Shear Fore: is the algebrai su of the vertial fores ating to the left or right of a ut setion along the span of the bea Bending oent: is the algebrai su of the oent of the fores
More informationDynamical behaviour of a compound elastic pendulum
MATEC Web of Conferenes 45, 3 (8) NCTAM 7 https://doi.org/.5/ateonf/8453 Dynaial behaviour of a opound elasti pendulu Svetoslav Nikolov,,* and Daniela Zaharieva Institute of Mehanis-BAS, Bioehanis Departent,
More informationAnnouncements Review: Relativistic mechanics Room: G1B30 From last class: total energy Example: Deuterium fusion Example: Deuterium fusion
Announeents Review: Relativisti ehanis Reading for Monday: Chapters 1 &! Relativisti oentu: dr p propper =γ u HW 4 due Wed. Do it before the ea! a 1 in 4 days. It overs Chapters 1 &. Roo: G1B3 (net to
More informationSubject: Introduction to Component Matching and Off-Design Operation % % ( (1) R T % (
16.50 Leture 0 Subjet: Introdution to Component Mathing and Off-Design Operation At this point it is well to reflet on whih of the many parameters we have introdued (like M, τ, τ t, ϑ t, f, et.) are free
More information' ' , and z ' components ( u u u'
Mesosale Meteorology: Gravity Waves 3 April 07 Introdution Here, we priarily onsider internal gravity waves, or waves that propagate in a density-stratified fluid (noinally, a stably-stratified fluid,
More informationA new Lagrangian of the simple harmonic oscillator 1 revisited
A new Lagrangian of the siple haronic oscillator 1 revisited Faisal Ain Yassein Abdelohssin Sudan Institute for Natural Sciences, P.O.BOX 3045, Khartou, Sudan Abstract A better and syetric new Lagrangian
More informationChapter 2: Solution of First order ODE
0 Chapter : Solution of irst order ODE Se. Separable Equations The differential equation of the form that is is alled separable if f = h g; In order to solve it perform the following steps: Rewrite the
More informationCOUPLED OSCILLATORS. Two identical pendulums
COUPED OSCIATORS A real physical object can be rearded as a lare nuber of siple oscillators coupled toether (atos and olecules in solids. The question is: how does the couplin affect the behavior of each
More information22.54 Neutron Interactions and Applications (Spring 2004) Chapter 6 (2/24/04) Energy Transfer Kernel F(E E')
22.54 Neutron Interations and Appliations (Spring 2004) Chapter 6 (2/24/04) Energy Transfer Kernel F(E E') Referenes -- J. R. Lamarsh, Introdution to Nulear Reator Theory (Addison-Wesley, Reading, 1966),
More informationOptimal sliding mode control of the pendubot
International Researh Journal of Coputer Siene and Inforation Systes (IRJCSIS Vol. ( pp. 45-5, April, Available online http://www.interesjournals.org/irjcsis Copyright International Researh Journals Full
More informationThe Lagrangian Method vs. other methods (COMPARATIVE EXAMPLE)
The Lagrangian ethod vs. other ethods () This aterial written by Jozef HANC, jozef.hanc@tuke.sk Technical University, Kosice, Slovakia For Edwin Taylor s website http://www.eftaylor.co/ 6 January 003 The
More informationPhysically Based Modeling CS Notes Spring 1997 Particle Collision and Contact
Physically Based Modeling CS 15-863 Notes Spring 1997 Particle Collision and Contact 1 Collisions with Springs Suppose we wanted to ipleent a particle siulator with a floor : a solid horizontal plane which
More information9 HOOKE S LAW AND SIMPLE HARMONIC MOTION
Experient 9 HOOKE S LAW AND SIMPLE HARMONIC MOTION Objectives 1. Verify Hoo s law,. Measure the force constant of a spring, and 3. Measure the period of oscillation of a spring-ass syste and copare it
More informationFrame with 6 DOFs. v m. determining stiffness, k k = F / water tower deflected water tower dynamic response model
CE 533, Fall 2014 Undaped SDOF Oscillator 1 / 6 What is a Single Degree of Freedo Oscillator? The siplest representation of the dynaic response of a civil engineering structure is the single degree of
More informationDeveloping Excel Macros for Solving Heat Diffusion Problems
Session 50 Developing Exel Maros for Solving Heat Diffusion Problems N. N. Sarker and M. A. Ketkar Department of Engineering Tehnology Prairie View A&M University Prairie View, TX 77446 Abstrat This paper
More informationCHAPTER 15: Vibratory Motion
CHAPTER 15: Vibratory Motion courtesy of Richard White courtesy of Richard White 2.) 1.) Two glaring observations can be ade fro the graphic on the previous slide: 1.) The PROJECTION of a point on a circle
More informationMass Transfer (Stoffaustausch) Fall 2012
Mass Transfer (Stoffaustaush) Fall Examination 9. Januar Name: Legi-Nr.: Edition Diffusion by E. L. Cussler: none nd rd Test Duration: minutes The following materials are not permitted at your table and
More informationMassachusetts Institute of Technology Quantum Mechanics I (8.04) Spring 2005 Solutions to Problem Set 4
Massachusetts Institute of Technology Quantu Mechanics I (8.04) Spring 2005 Solutions to Proble Set 4 By Kit Matan 1. X-ray production. (5 points) Calculate the short-wavelength liit for X-rays produced
More informationMAIN TOPICS iensional i l Analysis Bukingha Pi Theore eterination of Pi Ters Coents about iensional Analysis Coon iensionless Groups in Fluid Mehanis
FUNAMENTALS OF FLUI MECHANICS Chapter 7 iensional Analysis Modeling, and Siilitude MAIN TOPICS iensional i l Analysis Bukingha Pi Theore eterination of Pi Ters Coents about iensional Analysis Coon iensionless
More informationPhysics (Theory) There are 30 questions in total. Question Nos. 1 to 8 are very short answer type questions and carry one mark each.
Physis (Theory) Tie allowed: 3 hours] [Maxiu arks:7 General Instrutions: (i) ll uestions are opulsory. (ii) (iii) (iii) (iv) (v) There are 3 uestions in total. Question Nos. to 8 are very short answer
More informationA NETWORK SIMPLEX ALGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM
NETWORK SIMPLEX LGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM Cen Çalışan, Utah Valley University, 800 W. University Parway, Orem, UT 84058, 801-863-6487, en.alisan@uvu.edu BSTRCT The minimum
More informationTWO WAYS TO DISTINGUISH ONE INERTIAL FRAME FROM ANOTHER
TWO WAYS TO DISTINGUISH ONE INERTIAL FRAME FROM ANOTHER (No general ausality without superluminal veloities) by Dr. Tamas Lajtner Correspondene via web site: www.lajtnemahine.om ABSTRACT...2 1. SPACETIME
More informationFinal Exam: know your section, bring your ID!
Chapter 15: Equilibrium Part 1 Read: BLB 15.1 3 HW: BLB 15:13,14, 21 Supplemental 15:1 4 Know: Chemial Equilibrium Catalysts Equilibrium Constant Equilibrium onstant expression Homogeneous/Heterogeneous
More informationDefinition of Work, The basics
Physics 07 Lecture 16 Lecture 16 Chapter 11 (Work) v Eploy conservative and non-conservative forces v Relate force to potential energy v Use the concept of power (i.e., energy per tie) Chapter 1 v Define
More informationUnit 14 Harmonic Motion. Your Comments
Today s Concepts: Periodic Motion Siple - Mass on spring Daped Forced Resonance Siple - Pendulu Unit 1, Slide 1 Your Coents Please go through the three equations for siple haronic otion and phase angle
More informationMECHANICS OF MATERIALS
00 The Graw-Hill Copanies, n. All rights reserved. Third E CHAPTER Pure ECHANCS OF ATERALS Ferdinand P. Beer E. Russell Johnston, Jr. John T. DeWolf Bending Leture Notes: J. Walt Oler Teas Teh Universit
More informationChapter 15 Equilibrium. Reversible Reactions & Equilibrium. Reversible Reactions & Equilibrium. Reversible Reactions & Equilibrium 2/3/2014
Amount of reatant/produt //01 quilibrium in Chemial Reations Lets look bak at our hypothetial reation from the kinetis hapter. A + B C Chapter 15 quilibrium [A] Why doesn t the onentration of A ever go
More informationCDS 101/110: Lecture 2.1 System Modeling. Model-Based Analysis of Feedback Systems
CDS 11/11: Leture 21 System Modeling Rihard M Murray 5 Otober 215 Goals:! Define a model and its use in answering questions about a system! Introdue the onepts of state, dynamis, inputs and outputs! Review
More informationDynamics of Bass Reflex Loudspeaker Systems (3)
MCAP5E Dynaics of Bass Refle Loudspeaker Systes (3) Deriving Equations of Motion of Bass Refle Speaker Systes 3. Deriving Equations of Motion Shigeru Suzuki Released in May 5th, 8 in Japanese Released
More informationMath 151 Introduction to Eigenvectors
Math 151 Introdution to Eigenvetors The motivating example we used to desrie matrixes was landsape hange and vegetation suession. We hose the simple example of Bare Soil (B), eing replaed y Grasses (G)
More informationChapter 5, Conceptual Questions
Chapter 5, Conceptual Questions 5.1. Two forces are present, tension T in the cable and gravitational force 5.. F G as seen in the figure. Four forces act on the block: the push of the spring F, sp gravitational
More informationModern Control Systems (ECEG-4601) Instructor: Andinet Negash. Chapter 1 Lecture 3: State Space, II
Modern Control Systes (ECEG-46) Instructor: Andinet Negash Chapter Lecture 3: State Space, II Eaples Eaple 5: control o liquid levels: in cheical plants, it is oten necessary to aintain the levels o liquids.
More informationName Solutions to Test 1 September 23, 2016
Name Solutions to Test 1 September 3, 016 This test onsists of three parts. Please note that in parts II and III, you an skip one question of those offered. Possibly useful formulas: F qequb x xvt E Evpx
More informationDistillation. The Continuous Column. Learning Outcomes. Recap - VLE for Meth H 2 O. Gavin Duffy School of Electrical Engineering DIT Kevin Street
Distillation The Continuous Colun Gavin Duffy School of Electrical Engineering DIT Kevin Street Learning Outcoes After this lecture you should be able to.. Describe how continuous distillation works List
More informationA Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair
Proceedings of the 6th SEAS International Conference on Siulation, Modelling and Optiization, Lisbon, Portugal, Septeber -4, 006 0 A Siplified Analytical Approach for Efficiency Evaluation of the eaving
More informationAnalytical Analysis and Numerical Prediction of Seven DOF Human Vibratory Model for the Various Cars Driving Posture Swami Mahesh, Kosbe Pradnya
ISSN: - ISO 9:8 Certified International Journal of Engineering and Innovative Tehnology (IJEIT) Volue, Issue, May Analytial Analysis and Nuerial Predition of Seven DOF Huan Vibratory Model for the Various
More informationDetermination of the reaction order
5/7/07 A quote of the wee (or amel of the wee): Apply yourself. Get all the eduation you an, but then... do something. Don't just stand there, mae it happen. Lee Iaoa Physial Chemistry GTM/5 reation order
More information27 Oscillations: Introduction, Mass on a Spring
Chapter 7 Oscillations: Introduction, Mass on a Spring 7 Oscillations: Introduction, Mass on a Spring If a siple haronic oscillation proble does not involve the tie, you should probably be using conservation
More informationDETERMINATION OF TWO LAYER EARTH STRUCTURE PARAMETERS
DETERMINATION OF TWO LAYER EARTH STRUCTURE PARAMETERS Ioannis F. GONOS, Vassiliki T. KONTARGYRI, Ioannis A. STATHOPULOS, National Tehnial University of Athens, Shool of Eletrial and Coputer Engineering,
More informationThe Laws of Acceleration
The Laws of Aeleration The Relationships between Time, Veloity, and Rate of Aeleration Copyright 2001 Joseph A. Rybzyk Abstrat Presented is a theory in fundamental theoretial physis that establishes the
More information