22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 9: The High Beta Tokamak

Size: px
Start display at page:

Download "22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 9: The High Beta Tokamak"

Transcription

1 .65, MHD Theory of Fusion Sysems Prof. Freidberg Lecure 9: The High e Tokmk Summry of he Properies of n Ohmic Tokmk. Advnges:. good euilibrium (smll shif) b. good sbiliy ( ) c. good confinemen ( τ nr ) d. good ohmic heing ( T T kev ). Disdvnges: e i. low β p β : β, β = = % b. exernl heing is reuired: joule heing is no deue. Exernl heing is expensive, bu lso rises β c. igh spec rio is reuired o rise β : his is echnologiclly difficul d. pulsed operion is reuired unless curren drive works efficienly 3. The high β okmk resolves he problem of low β. I llows okmk o opere higher β 5 %. This leds o more economic devices b. How do we chieve higher β? We pply ddiionl uxiliry heing, keeping fixed. This rises p relive o High β Tokmk Expnsion μ.. We gin ssume lrge spec rio: R R. Sbiliy is produced by lrge oroidl field: 3. Thus, s in he ohmic okmk r R, implying h.65, MHD Theory of Fusion Sysems Lecure 9 Prof. Freidberg Pge of

2 4. Rdil pressure blnce, however is produced by he oroidl field = + δ dimgneic ' ' p + + ( r ) = μ μ r δ + p + μ μ r ' ' ( r ) = neglec, confines only β 5. To improve β over h chievble in he ohmic okmk we need p δ μ μ 6. Or in erms of β δ β 7. How lrge cn β nd δ ge? 8. The limiing condiion is deermined by oroidl force blnce which is sill ccomplished by combinion of I nd v 9. Incresing β increses he oroidl shif. The lrges possible β occurs when he shif becomes of order uniy Δ. Recll h in n ohmic okmk Δ R..65, MHD Theory of Fusion Sysems Lecure 9 Prof. Freidberg Pge of

3 . Le us esime he shif using he smll shif relion ' ' b dr r r ' μ dp = y y dy r dy μp neglec since μ p. Therefore Δ R R μp μp β β ' p. For, hen Δ when β 3. This suggess he following ordering for he high β okmk δ β,, Δ Comprison of Expnsions Ohmic Tokmk High e Tokmk β μ p δ (pr) δ β μ p p β μ p (di) βp μ p Expnsion of Grd-Shfrnov Euion. Since Δ, oroidl force blnce nd rdil pressure blnce ener ogeher in zeroh order.. Good news: we need only he zeroh order euions. No firs order correcions re reuired. 3. d news: The zero oher euions re sill nonliner pril differenil euions..65, MHD Theory of Fusion Sysems Lecure 9 Prof. Freidberg Pge 3 of

4 4. Expnsion:. = r, +..., rr b. p = p +..., μ p c. F R = F ( ) F = R p μ + new free funcion d. This uomiclly produces pinch pressure blnce p + μ cons 5. Subsiue he expnsion ino he Grd-Shfrnov euion dp d F = μ ( R + r cos) + cos sin = T R r r T cos = sin R r r R 3 neglec dp d F T = μ ( R + r cos ) d F dp + μrp μrr cos +... = d dp = R μp + + μ p μrr cos = R d ( ) R dp μrr cos μp.65, MHD Theory of Fusion Sysems Lecure 9 Prof. Freidberg Pge 4 of

5 6. Therefore, o leding order he Grd-Shfrnov euion reduces o d = dp ( ) R μrr cos 7. Noe h μ RJ, so h on circulr plsm flux surfce μrj = R d verge over We see h d flux surfce. is proporionl o he verge oroidl curren wihin given 8. Even hough he euion is simpler, i is sill nonliner PDE. 9. In generl, i mus be solved numericlly.. The difficuly rises becuse he shifs re finie nd cnno be reed perurbively.. We shll deermine generl feures of high β okmk by exmining specil cse. Specil Cse. Choose dp μr = C C = cons d R = A A= cons. This implies p C (ssume (boundry)=) nd J A.65, MHD Theory of Fusion Sysems Lecure 9 Prof. Freidberg Pge 5 of

6 3. Soluion: wih hese choices he Grd-Shfrnov euion becomes r + = A + Cr cos r r r r 4. oundry condiions: We ssume circulr plsm of rdius r = (, ) = (normlizion of flux funcion is rbirry) ( r, ) regulr for r The circulr ssumpion is mde for simpliciy nd cn be generlized o oher cross secions. 5. Soluion: (We need only cos n erms becuse of up-down symmery.) 3 r r pr = A + C cos 4 8 n n k4 cos hom 3 n n r n= = k + k ln r + k r cos + + r + b r cos n no regulr no regulr no regulr = krcos + k 3 only erms which re regulr nd reuired o blnce on he boundry r = pr For ll n, i follows h = n k ( ) 6. Choose k nd o mke 3 A C 3 ( r ) ( r ) ( r r), =, = + cos Then C = = Ar 3r cos R r R + 4 R R 3 C C C μ p = = A r + r r cos 8.65, MHD Theory of Fusion Sysems Lecure 9 Prof. Freidberg Pge 6 of

7 Physicl Inerpreion Le us express A nd C in erms of more physicl uniies: β, I or euivlenly β, * I. * is prmeer reled o kink sbiliy nd he surfce MHD sfey fcor * A p μ RI π = μ RI. I = p dl = (, ) d μ ON A CIRCLE π C A πa = A d = R + cos π = R R 3. μr πa A = = * π R μ 4. μ μ β = p = pr dr d π C + μ C π 3 = r dr d A r r r cos 8 R π μ 4 π 4 5. Thus = AC = 8R 4R A = * C * C 4R = * β.65, MHD Theory of Fusion Sysems Lecure 9 Prof. Freidberg Pge 7 of

8 6. Generl euilibrium relion for high β circulr okmk (minor degression). β p = pr dr d μ I 8 π 8π π = β μ I μ π = μi β * β, = R b. The generl euilibrium relion is given by β = βp * β* 7. Subsiue A nd C bck ino he soluions. Define ν = = βp, ρ = r * 3 = ρ + ν ρ ρ cos ν = + 3 cos ρ ρ * r μ p ( ρ ) ν = sin * = β ρ + νρ cos μ RJ = + ν cos * ( ρ ) cos - cos = ρ β ρ + νρ.65, MHD Theory of Fusion Sysems Lecure 9 Prof. Freidberg Pge 8 of

9 Properies of High e Tokmk Euilibri. Skeched below re ypicl midplne profiles (Z=) showing rdil pressure blnce.. Shown here re p nd J profiles long he midplne for differen ν. Incresing ν implies higher β..65, MHD Theory of Fusion Sysems Lecure 9 Prof. Freidberg Pge 9 of

10 . Observe he incresed shif of he mgneic xis s β increses. b. Observe he buildup of curren on he ouside of he orus o produce oroidl force blnce higher β. c. Observe J reversing on he inside of he orus when ν >. 3. The flux surfces re round, bu re no circles excep for he boundry 3 ρ + ν ρ ρ cos =cons 4. Le us clcule he mgneic xis shif Δ by finding he vlue of r where p r r. A = =. y symmery his occurs when =( orπ)..65, MHD Theory of Fusion Sysems Lecure 9 Prof. Freidberg Pge of

11 3 ρ + ν ρ ρ b. Se Δ ρ =, = = ρ c. This yields Δ + ν 3 = Δ d. The vlue of Δ is given by Δ ν = ( ν ) + +3 Δ e. For he HT ν Δ ν Ohmic Tν 5. Find he shpe of he flux surfces ner he mgneic xis.. Le x = ρ cos y = ρ sin b. Then x + y + ν x Δ c. Expnd x = + δ x y = δy δx, δy nd define x =Δ d. Subsiue x + x δ x + δ x + δy + νx + νδ x = cons = x + ν x + δ x x + x + x ν ν ν or + δ y + νx + δ x + νx + ν x definiion of x.65, MHD Theory of Fusion Sysems Lecure 9 Prof. Freidberg Pge of

12 ( νx )( x) ( νx )( y) + 3 δ + + δ = cons e. This is he euion of n ellipse elonged flux surfces, sushed ner he ouside f. The elongion κ is defined by b + 3ν x + ν x κ = = ( 3ν ) = + ( ν ) ν ( ν ) κ for ν nd κ = 3 for ν =..65, MHD Theory of Fusion Sysems Lecure 9 Prof. Freidberg Pge of

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 10: The High Beta Tokamak Con d and the High Flux Conserving Tokamak.

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 10: The High Beta Tokamak Con d and the High Flux Conserving Tokamak. .615, MHD Theory of Fusion Sysems Prof. Freidberg Lecure 1: The High Be Tokmk Con d nd he High Flux Conserving Tokmk Proeries of he High Tokmk 1. Evlue he MHD sfey fcor: ψ B * ( ) 1 3 ρ 1+ ν ρ ρ cosθ *

More information

September 20 Homework Solutions

September 20 Homework Solutions College of Engineering nd Compuer Science Mechnicl Engineering Deprmen Mechnicl Engineering A Seminr in Engineering Anlysis Fll 7 Number 66 Insrucor: Lrry Creo Sepember Homework Soluions Find he specrum

More information

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 7: The First Order Grad Shafranov Equation. dp 1 dp

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 7: The First Order Grad Shafranov Equation. dp 1 dp First Order Eqution.65, MHD Theory of Fusion Systems Prof. Freiderg Lecture 7: The First Order Grd Shfrnov Eqution The first order Grd Shfrnov eqution is given y d p d dp d + μr + R B B = μr r cos + cos

More information

e t dt e t dt = lim e t dt T (1 e T ) = 1

e t dt e t dt = lim e t dt T (1 e T ) = 1 Improper Inegrls There re wo ypes of improper inegrls - hose wih infinie limis of inegrion, nd hose wih inegrnds h pproch some poin wihin he limis of inegrion. Firs we will consider inegrls wih infinie

More information

5.1-The Initial-Value Problems For Ordinary Differential Equations

5.1-The Initial-Value Problems For Ordinary Differential Equations 5.-The Iniil-Vlue Problems For Ordinry Differenil Equions Consider solving iniil-vlue problems for ordinry differenil equions: (*) y f, y, b, y. If we know he generl soluion y of he ordinry differenil

More information

MTH 146 Class 11 Notes

MTH 146 Class 11 Notes 8.- Are of Surfce of Revoluion MTH 6 Clss Noes Suppose we wish o revolve curve C round n is nd find he surfce re of he resuling solid. Suppose f( ) is nonnegive funcion wih coninuous firs derivive on he

More information

Magnetostatics Bar Magnet. Magnetostatics Oersted s Experiment

Magnetostatics Bar Magnet. Magnetostatics Oersted s Experiment Mgneosics Br Mgne As fr bck s 4500 yers go, he Chinese discovered h cerin ypes of iron ore could rc ech oher nd cerin mels. Iron filings "mp" of br mgne s field Crefully suspended slivers of his mel were

More information

3. Renewal Limit Theorems

3. Renewal Limit Theorems Virul Lborories > 14. Renewl Processes > 1 2 3 3. Renewl Limi Theorems In he inroducion o renewl processes, we noed h he rrivl ime process nd he couning process re inverses, in sens The rrivl ime process

More information

4.8 Improper Integrals

4.8 Improper Integrals 4.8 Improper Inegrls Well you ve mde i hrough ll he inegrion echniques. Congrs! Unforunely for us, we sill need o cover one more inegrl. They re clled Improper Inegrls. A his poin, we ve only del wih inegrls

More information

The solution is often represented as a vector: 2xI + 4X2 + 2X3 + 4X4 + 2X5 = 4 2xI + 4X2 + 3X3 + 3X4 + 3X5 = 4. 3xI + 6X2 + 6X3 + 3X4 + 6X5 = 6.

The solution is often represented as a vector: 2xI + 4X2 + 2X3 + 4X4 + 2X5 = 4 2xI + 4X2 + 3X3 + 3X4 + 3X5 = 4. 3xI + 6X2 + 6X3 + 3X4 + 6X5 = 6. [~ o o :- o o ill] i 1. Mrices, Vecors, nd Guss-Jordn Eliminion 1 x y = = - z= The soluion is ofen represened s vecor: n his exmple, he process of eliminion works very smoohly. We cn elimine ll enries

More information

ENGR 1990 Engineering Mathematics The Integral of a Function as a Function

ENGR 1990 Engineering Mathematics The Integral of a Function as a Function ENGR 1990 Engineering Mhemics The Inegrl of Funcion s Funcion Previously, we lerned how o esime he inegrl of funcion f( ) over some inervl y dding he res of finie se of rpezoids h represen he re under

More information

0 for t < 0 1 for t > 0

0 for t < 0 1 for t > 0 8.0 Sep nd del funcions Auhor: Jeremy Orloff The uni Sep Funcion We define he uni sep funcion by u() = 0 for < 0 for > 0 I is clled he uni sep funcion becuse i kes uni sep = 0. I is someimes clled he Heviside

More information

ECE Microwave Engineering

ECE Microwave Engineering EE 537-635 Microwve Engineering Adped from noes y Prof. Jeffery T. Willims Fll 8 Prof. Dvid R. Jcson Dep. of EE Noes Wveguiding Srucures Pr 7: Trnsverse Equivlen Newor (N) Wveguide Trnsmission Line Model

More information

T-Match: Matching Techniques For Driving Yagi-Uda Antennas: T-Match. 2a s. Z in. (Sections 9.5 & 9.7 of Balanis)

T-Match: Matching Techniques For Driving Yagi-Uda Antennas: T-Match. 2a s. Z in. (Sections 9.5 & 9.7 of Balanis) 3/0/018 _mch.doc Pge 1 of 6 T-Mch: Mching Techniques For Driving Ygi-Ud Anenns: T-Mch (Secions 9.5 & 9.7 of Blnis) l s l / l / in The T-Mch is shun-mching echnique h cn be used o feed he driven elemen

More information

S Radio transmission and network access Exercise 1-2

S Radio transmission and network access Exercise 1-2 S-7.330 Rdio rnsmission nd nework ccess Exercise 1 - P1 In four-symbol digil sysem wih eqully probble symbols he pulses in he figure re used in rnsmission over AWGN-chnnel. s () s () s () s () 1 3 4 )

More information

A Kalman filtering simulation

A Kalman filtering simulation A Klmn filering simulion The performnce of Klmn filering hs been esed on he bsis of wo differen dynmicl models, ssuming eiher moion wih consn elociy or wih consn ccelerion. The former is epeced o beer

More information

Average & instantaneous velocity and acceleration Motion with constant acceleration

Average & instantaneous velocity and acceleration Motion with constant acceleration Physics 7: Lecure Reminders Discussion nd Lb secions sr meeing ne week Fill ou Pink dd/drop form if you need o swich o differen secion h is FULL. Do i TODAY. Homework Ch. : 5, 7,, 3,, nd 6 Ch.: 6,, 3 Submission

More information

REAL ANALYSIS I HOMEWORK 3. Chapter 1

REAL ANALYSIS I HOMEWORK 3. Chapter 1 REAL ANALYSIS I HOMEWORK 3 CİHAN BAHRAN The quesions re from Sein nd Shkrchi s e. Chper 1 18. Prove he following sserion: Every mesurble funcion is he limi.e. of sequence of coninuous funcions. We firs

More information

ECE Microwave Engineering. Fall Prof. David R. Jackson Dept. of ECE. Notes 10. Waveguides Part 7: Transverse Equivalent Network (TEN)

ECE Microwave Engineering. Fall Prof. David R. Jackson Dept. of ECE. Notes 10. Waveguides Part 7: Transverse Equivalent Network (TEN) EE 537-635 Microwve Engineering Fll 7 Prof. Dvid R. Jcson Dep. of EE Noes Wveguides Pr 7: Trnsverse Equivlen Newor (N) Wveguide Trnsmission Line Model Our gol is o come up wih rnsmission line model for

More information

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 22

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 22 allooning mode equaion π δ W = ψ ψ μ d W.615, MHD Theory o Fusion Sysems Pro. Freidberg Lecure μ R W Jd k k k k k J χ dψ 1 X d ( ψ ) = χ ( n + ) ( κn nκ ) Mercier Crierion X 1. In using he quasimode reresenaion

More information

Solutions for Nonlinear Partial Differential Equations By Tan-Cot Method

Solutions for Nonlinear Partial Differential Equations By Tan-Cot Method IOSR Journl of Mhemics (IOSR-JM) e-issn: 78-578. Volume 5, Issue 3 (Jn. - Feb. 13), PP 6-11 Soluions for Nonliner Pril Differenil Equions By Tn-Co Mehod Mhmood Jwd Abdul Rsool Abu Al-Sheer Al -Rfidin Universiy

More information

Procedia Computer Science

Procedia Computer Science Procedi Compuer Science 00 (0) 000 000 Procedi Compuer Science www.elsevier.com/loce/procedi The Third Informion Sysems Inernionl Conference The Exisence of Polynomil Soluion of he Nonliner Dynmicl Sysems

More information

1.0 Electrical Systems

1.0 Electrical Systems . Elecricl Sysems The ypes of dynmicl sysems we will e sudying cn e modeled in erms of lgeric equions, differenil equions, or inegrl equions. We will egin y looking fmilir mhemicl models of idel resisors,

More information

INTEGRALS. Exercise 1. Let f : [a, b] R be bounded, and let P and Q be partitions of [a, b]. Prove that if P Q then U(P ) U(Q) and L(P ) L(Q).

INTEGRALS. Exercise 1. Let f : [a, b] R be bounded, and let P and Q be partitions of [a, b]. Prove that if P Q then U(P ) U(Q) and L(P ) L(Q). INTEGRALS JOHN QUIGG Eercise. Le f : [, b] R be bounded, nd le P nd Q be priions of [, b]. Prove h if P Q hen U(P ) U(Q) nd L(P ) L(Q). Soluion: Le P = {,..., n }. Since Q is obined from P by dding finiely

More information

MA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions

MA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions MA 14 Calculus IV (Spring 016) Secion Homework Assignmen 1 Soluions 1 Boyce and DiPrima, p 40, Problem 10 (c) Soluion: In sandard form he given firs-order linear ODE is: An inegraing facor is given by

More information

MATH 124 AND 125 FINAL EXAM REVIEW PACKET (Revised spring 2008)

MATH 124 AND 125 FINAL EXAM REVIEW PACKET (Revised spring 2008) MATH 14 AND 15 FINAL EXAM REVIEW PACKET (Revised spring 8) The following quesions cn be used s review for Mh 14/ 15 These quesions re no cul smples of quesions h will pper on he finl em, bu hey will provide

More information

M r. d 2. R t a M. Structural Mechanics Section. Exam CT5141 Theory of Elasticity Friday 31 October 2003, 9:00 12:00 hours. Problem 1 (3 points)

M r. d 2. R t a M. Structural Mechanics Section. Exam CT5141 Theory of Elasticity Friday 31 October 2003, 9:00 12:00 hours. Problem 1 (3 points) Delf Universiy of Technology Fculy of Civil Engineering nd Geosciences Srucurl echnics Secion Wrie your nme nd sudy numer he op righ-hnd of your work. Exm CT5 Theory of Elsiciy Fridy Ocoer 00, 9:00 :00

More information

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion

More information

Minimum Squared Error

Minimum Squared Error Minimum Squred Error LDF: Minimum Squred-Error Procedures Ide: conver o esier nd eer undersood prolem Percepron y i > for ll smples y i solve sysem of liner inequliies MSE procedure y i = i for ll smples

More information

Minimum Squared Error

Minimum Squared Error Minimum Squred Error LDF: Minimum Squred-Error Procedures Ide: conver o esier nd eer undersood prolem Percepron y i > 0 for ll smples y i solve sysem of liner inequliies MSE procedure y i i for ll smples

More information

Some Inequalities variations on a common theme Lecture I, UL 2007

Some Inequalities variations on a common theme Lecture I, UL 2007 Some Inequliies vriions on common heme Lecure I, UL 2007 Finbrr Hollnd, Deprmen of Mhemics, Universiy College Cork, fhollnd@uccie; July 2, 2007 Three Problems Problem Assume i, b i, c i, i =, 2, 3 re rel

More information

Properties of Logarithms. Solving Exponential and Logarithmic Equations. Properties of Logarithms. Properties of Logarithms. ( x)

Properties of Logarithms. Solving Exponential and Logarithmic Equations. Properties of Logarithms. Properties of Logarithms. ( x) Properies of Logrihms Solving Eponenil nd Logrihmic Equions Properies of Logrihms Produc Rule ( ) log mn = log m + log n ( ) log = log + log Properies of Logrihms Quoien Rule log m = logm logn n log7 =

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

Jackson 2.7 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

Jackson 2.7 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell Jckson.7 Homework Problem Solution Dr. Christopher S. Bird University of Msschusetts Lowell PROBLEM: Consider potentil problem in the hlf-spce defined by, with Dirichlet boundry conditions on the plne

More information

Physics 2A HW #3 Solutions

Physics 2A HW #3 Solutions Chper 3 Focus on Conceps: 3, 4, 6, 9 Problems: 9, 9, 3, 41, 66, 7, 75, 77 Phsics A HW #3 Soluions Focus On Conceps 3-3 (c) The ccelerion due o grvi is he sme for boh blls, despie he fc h he hve differen

More information

IX.1.1 The Laplace Transform Definition 700. IX.1.2 Properties 701. IX.1.3 Examples 702. IX.1.4 Solution of IVP for ODEs 704

IX.1.1 The Laplace Transform Definition 700. IX.1.2 Properties 701. IX.1.3 Examples 702. IX.1.4 Solution of IVP for ODEs 704 Chper IX The Inegrl Trnform Mehod IX. The plce Trnform November 4, 7 699 IX. THE APACE TRANSFORM IX.. The plce Trnform Definiion 7 IX.. Properie 7 IX..3 Emple 7 IX..4 Soluion of IVP for ODE 74 IX..5 Soluion

More information

Contraction Mapping Principle Approach to Differential Equations

Contraction Mapping Principle Approach to Differential Equations epl Journl of Science echnology 0 (009) 49-53 Conrcion pping Principle pproch o Differenil Equions Bishnu P. Dhungn Deprmen of hemics, hendr Rn Cmpus ribhuvn Universiy, Khmu epl bsrc Using n eension of

More information

( ) ( ) ( ) ( ) ( ) ( y )

( ) ( ) ( ) ( ) ( ) ( y ) 8. Lengh of Plne Curve The mos fmous heorem in ll of mhemics is he Pyhgoren Theorem. I s formulion s he disnce formul is used o find he lenghs of line segmens in he coordine plne. In his secion you ll

More information

Probability, Estimators, and Stationarity

Probability, Estimators, and Stationarity Chper Probbiliy, Esimors, nd Sionriy Consider signl genered by dynmicl process, R, R. Considering s funcion of ime, we re opering in he ime domin. A fundmenl wy o chrcerize he dynmics using he ime domin

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

6. Gas dynamics. Ideal gases Speed of infinitesimal disturbances in still gas

6. Gas dynamics. Ideal gases Speed of infinitesimal disturbances in still gas 6. Gs dynmics Dr. Gergely Krisóf De. of Fluid echnics, BE Februry, 009. Seed of infiniesiml disurbnces in sill gs dv d, dv d, Coninuiy: ( dv)( ) dv omenum r r heorem: ( ( dv) ) d 3443 4 q m dv d dv llievi

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

More information

A new model for limit order book dynamics

A new model for limit order book dynamics Anewmodelforlimiorderbookdynmics JeffreyR.Russell UniversiyofChicgo,GrdueSchoolofBusiness TejinKim UniversiyofChicgo,DeprmenofSisics Absrc:Thispperproposesnewmodelforlimiorderbookdynmics.Thelimiorderbookconsiss

More information

FM Applications of Integration 1.Centroid of Area

FM Applications of Integration 1.Centroid of Area FM Applicions of Inegrion.Cenroid of Are The cenroid of ody is is geomeric cenre. For n ojec mde of uniform meril, he cenroid coincides wih he poin which he ody cn e suppored in perfecly lnced se ie, is

More information

ME 391 Mechanical Engineering Analysis

ME 391 Mechanical Engineering Analysis Fall 04 ME 39 Mechanical Engineering Analsis Eam # Soluions Direcions: Open noes (including course web posings). No books, compuers, or phones. An calculaor is fair game. Problem Deermine he posiion of

More information

EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR A SECOND-ORDER ITERATIVE BOUNDARY-VALUE PROBLEM

EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR A SECOND-ORDER ITERATIVE BOUNDARY-VALUE PROBLEM Elecronic Journl of Differenil Equions, Vol. 208 (208), No. 50, pp. 6. ISSN: 072-669. URL: hp://ejde.mh.xse.edu or hp://ejde.mh.un.edu EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR A SECOND-ORDER ITERATIVE

More information

Reduction of Cross-Field Plasma Transport in Tokamaks due to Power Input Redistribution and Sheared Flow Profile Modification

Reduction of Cross-Field Plasma Transport in Tokamaks due to Power Input Redistribution and Sheared Flow Profile Modification TH/P8-6 Reducion of Cross-Field Plsm Trnspor in Tokmks due o Power Inpu Redisribuion nd Shered Flow Profile Modificion V.P. Psukhov, N.V. Chudin RRC «Kurchov Insiue», 8 Moscow, Russin Federion e-mil conc

More information

Chapter Direct Method of Interpolation

Chapter Direct Method of Interpolation Chper 5. Direc Mehod of Inerpolion Afer reding his chper, you should be ble o:. pply he direc mehod of inerpolion,. sole problems using he direc mehod of inerpolion, nd. use he direc mehod inerpolns o

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

Lecture 4 Notes (Little s Theorem)

Lecture 4 Notes (Little s Theorem) Lecure 4 Noes (Lile s Theorem) This lecure concerns one of he mos imporan (and simples) heorems in Queuing Theory, Lile s Theorem. More informaion can be found in he course book, Bersekas & Gallagher,

More information

..,..,.,

..,..,., 57.95. «..» 7, 9,,. 3 DOI:.459/mmph7..,..,., E-mil: yshr_ze@mil.ru -,,. -, -.. -. - - ( ). -., -. ( - ). - - -., - -., - -, -., -. -., - - -, -., -. : ; ; - ;., -,., - -, []., -, [].,, - [3, 4]. -. 3 (

More information

Chapter 7: Solving Trig Equations

Chapter 7: Solving Trig Equations Haberman MTH Secion I: The Trigonomeric Funcions Chaper 7: Solving Trig Equaions Le s sar by solving a couple of equaions ha involve he sine funcion EXAMPLE a: Solve he equaion sin( ) The inverse funcions

More information

5.7 Improper Integrals

5.7 Improper Integrals 458 pplictions of definite integrls 5.7 Improper Integrls In Section 5.4, we computed the work required to lift pylod of mss m from the surfce of moon of mss nd rdius R to height H bove the surfce of the

More information

Thermal neutron self-shielding factor in foils: a universal curve

Thermal neutron self-shielding factor in foils: a universal curve Proceedings of he Inernionl Conference on Reserch Recor Uilizion, Sfey, Decommissioning, Fuel nd Wse Mngemen (Snigo, Chile, -4 Nov.3) Pper IAEA-CN-/, IAEA Proceedings Series, Vienn, 5 Therml neuron self-shielding

More information

THE 2-BODY PROBLEM. FIGURE 1. A pair of ellipses sharing a common focus. (c,b) c+a ROBERT J. VANDERBEI

THE 2-BODY PROBLEM. FIGURE 1. A pair of ellipses sharing a common focus. (c,b) c+a ROBERT J. VANDERBEI THE 2-BODY PROBLEM ROBERT J. VANDERBEI ABSTRACT. In his shor noe, we show ha a pair of ellipses wih a common focus is a soluion o he 2-body problem. INTRODUCTION. Solving he 2-body problem from scrach

More information

Chapter 2: Evaluative Feedback

Chapter 2: Evaluative Feedback Chper 2: Evluive Feedbck Evluing cions vs. insrucing by giving correc cions Pure evluive feedbck depends olly on he cion ken. Pure insrucive feedbck depends no ll on he cion ken. Supervised lerning is

More information

IX.1.1 The Laplace Transform Definition 700. IX.1.2 Properties 701. IX.1.3 Examples 702. IX.1.4 Solution of IVP for ODEs 704

IX.1.1 The Laplace Transform Definition 700. IX.1.2 Properties 701. IX.1.3 Examples 702. IX.1.4 Solution of IVP for ODEs 704 Chper IX The Inegrl Trnform Mehod IX. The plce Trnform November 6, 8 699 IX. THE APACE TRANSFORM IX.. The plce Trnform Definiion 7 IX.. Properie 7 IX..3 Emple 7 IX..4 Soluion of IVP for ODE 74 IX..5 Soluion

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

More information

EE 315 Notes. Gürdal Arslan CLASS 1. (Sections ) What is a signal?

EE 315 Notes. Gürdal Arslan CLASS 1. (Sections ) What is a signal? EE 35 Noes Gürdal Arslan CLASS (Secions.-.2) Wha is a signal? In his class, a signal is some funcion of ime and i represens how some physical quaniy changes over some window of ime. Examples: velociy of

More information

Solutions to Problems from Chapter 2

Solutions to Problems from Chapter 2 Soluions o Problems rom Chper Problem. The signls u() :5sgn(), u () :5sgn(), nd u h () :5sgn() re ploed respecively in Figures.,b,c. Noe h u h () :5sgn() :5; 8 including, bu u () :5sgn() is undeined..5

More information

Chapter 2. Motion along a straight line. 9/9/2015 Physics 218

Chapter 2. Motion along a straight line. 9/9/2015 Physics 218 Chper Moion long srigh line 9/9/05 Physics 8 Gols for Chper How o describe srigh line moion in erms of displcemen nd erge elociy. The mening of insnneous elociy nd speed. Aerge elociy/insnneous elociy

More information

Y 0.4Y 0.45Y Y to a proper ARMA specification.

Y 0.4Y 0.45Y Y to a proper ARMA specification. HG Jan 04 ECON 50 Exercises II - 0 Feb 04 (wih answers Exercise. Read secion 8 in lecure noes 3 (LN3 on he common facor problem in ARMA-processes. Consider he following process Y 0.4Y 0.45Y 0.5 ( where

More information

P441 Analytical Mechanics - I. Coupled Oscillators. c Alex R. Dzierba

P441 Analytical Mechanics - I. Coupled Oscillators. c Alex R. Dzierba Lecure 3 Mondy - Deceber 5, 005 Wrien or ls upded: Deceber 3, 005 P44 Anlyicl Mechnics - I oupled Oscillors c Alex R. Dzierb oupled oscillors - rix echnique In Figure we show n exple of wo coupled oscillors,

More information

1 jordan.mcd Eigenvalue-eigenvector approach to solving first order ODEs. -- Jordan normal (canonical) form. Instructor: Nam Sun Wang

1 jordan.mcd Eigenvalue-eigenvector approach to solving first order ODEs. -- Jordan normal (canonical) form. Instructor: Nam Sun Wang jordnmcd Eigenvlue-eigenvecor pproch o solving firs order ODEs -- ordn norml (cnonicl) form Insrucor: Nm Sun Wng Consider he following se of coupled firs order ODEs d d x x 5 x x d d x d d x x x 5 x x

More information

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous

More information

Space Curves. Recall the parametric equations of a curve in xy-plane and compare them with parametric equations of a curve in space.

Space Curves. Recall the parametric equations of a curve in xy-plane and compare them with parametric equations of a curve in space. Clculus 3 Li Vs Spce Curves Recll the prmetric equtions of curve in xy-plne nd compre them with prmetric equtions of curve in spce. Prmetric curve in plne x = x(t) y = y(t) Prmetric curve in spce x = x(t)

More information

Laplace Transforms. Examples. Is this equation differential? y 2 2y + 1 = 0, y 2 2y + 1 = 0, (y ) 2 2y + 1 = cos x,

Laplace Transforms. Examples. Is this equation differential? y 2 2y + 1 = 0, y 2 2y + 1 = 0, (y ) 2 2y + 1 = cos x, Laplace Transforms Definiion. An ordinary differenial equaion is an equaion ha conains one or several derivaives of an unknown funcion which we call y and which we wan o deermine from he equaion. The equaion

More information

Transforms II - Wavelets Preliminary version please report errors, typos, and suggestions for improvements

Transforms II - Wavelets Preliminary version please report errors, typos, and suggestions for improvements EECS 3 Digil Signl Processing Universiy of Cliforni, Berkeley: Fll 007 Gspr November 4, 007 Trnsforms II - Wveles Preliminry version plese repor errors, ypos, nd suggesions for improvemens We follow n

More information

Physics 1402: Lecture 22 Today s Agenda

Physics 1402: Lecture 22 Today s Agenda Physics 142: ecure 22 Today s Agenda Announcemens: R - RV - R circuis Homework 6: due nex Wednesday Inducion / A curren Inducion Self-Inducance, R ircuis X X X X X X X X X long solenoid Energy and energy

More information

PHYSICS 1210 Exam 1 University of Wyoming 14 February points

PHYSICS 1210 Exam 1 University of Wyoming 14 February points PHYSICS 1210 Em 1 Uniersiy of Wyoming 14 Februry 2013 150 poins This es is open-noe nd closed-book. Clculors re permied bu compuers re no. No collborion, consulion, or communicion wih oher people (oher

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

F This leads to an unstable mode which is not observable at the output thus cannot be controlled by feeding back.

F This leads to an unstable mode which is not observable at the output thus cannot be controlled by feeding back. Lecure 8 Las ime: Semi-free configuraion design This is equivalen o: Noe ns, ener he sysem a he same place. is fixed. We design C (and perhaps B. We mus sabilize if i is given as unsable. Cs ( H( s = +

More information

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Swiss Federl Insiue of Pge 1 The Finie Elemen Mehod for he Anlysis of Non-Liner nd Dynmic Sysems Prof. Dr. Michel Hvbro Fber Dr. Nebojs Mojsilovic Swiss Federl Insiue of ETH Zurich, Swizerlnd Mehod of

More information

Phys 110. Answers to even numbered problems on Midterm Map

Phys 110. Answers to even numbered problems on Midterm Map Phys Answers o een numbered problems on Miderm Mp. REASONING The word per indices rio, so.35 mm per dy mens.35 mm/d, which is o be epressed s re in f/cenury. These unis differ from he gien unis in boh

More information

Mathematics 805 Final Examination Answers

Mathematics 805 Final Examination Answers . 5 poins Se he Weiersrss M-es. Mhemics 85 Finl Eminion Answers Answer: Suppose h A R, nd f n : A R. Suppose furher h f n M n for ll A, nd h Mn converges. Then f n converges uniformly on A.. 5 poins Se

More information

An integral having either an infinite limit of integration or an unbounded integrand is called improper. Here are two examples.

An integral having either an infinite limit of integration or an unbounded integrand is called improper. Here are two examples. Improper Inegrls To his poin we hve only considered inegrls f(x) wih he is of inegrion nd b finie nd he inegrnd f(x) bounded (nd in fc coninuous excep possibly for finiely mny jump disconinuiies) An inegrl

More information

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem) Week 1 Lecure Problems, 5 Wha if somehing oscillaes wih no obvious spring? Wha is ω? (problem se problem) Sar wih Try and ge o SHM form E. Full beer can in lake, oscillaing F = m & = ge rearrange: F =

More information

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux Gues Lecures for Dr. MacFarlane s EE3350 Par Deux Michael Plane Mon., 08-30-2010 Wrie name in corner. Poin ou his is a review, so I will go faser. Remind hem o go lisen o online lecure abou geing an A

More information

Electromagnetism Answers to Problem Set 10 Spring 2006

Electromagnetism Answers to Problem Set 10 Spring 2006 Electromgnetism 76 Answers to Problem Set 1 Spring 6 1. Jckson Prob. 5.15: Shielded Bifilr Circuit: Two wires crrying oppositely directed currents re surrounded by cylindricl shell of inner rdius, outer

More information

Power of Random Processes 1/40

Power of Random Processes 1/40 Power of Random Processes 40 Power of a Random Process Recall : For deerminisic signals insananeous power is For a random signal, is a random variable for each ime. hus here is no single # o associae wih

More information

ENGI 9420 Engineering Analysis Assignment 2 Solutions

ENGI 9420 Engineering Analysis Assignment 2 Solutions ENGI 940 Engineering Analysis Assignmen Soluions 0 Fall [Second order ODEs, Laplace ransforms; Secions.0-.09]. Use Laplace ransforms o solve he iniial value problem [0] dy y, y( 0) 4 d + [This was Quesion

More information

Lecture #8 Redfield theory of NMR relaxation

Lecture #8 Redfield theory of NMR relaxation Lecure #8 Redfield heory of NMR relaxaion Topics The ineracion frame of reference Perurbaion heory The Maser Equaion Handous and Reading assignmens van de Ven, Chapers 6.2. Kowalewski, Chaper 4. Abragam

More information

Problems on transformer main dimensions and windings

Problems on transformer main dimensions and windings Probles_Trn_winding Probles on rnsforer in diensions nd windings. Deerine he in diensions of he core nd window for 500 ka, /400, 50Hz, Single phse core ype, oil iersed, self cooled rnsforer. Assue: Flux

More information

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires

More information

DEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING EXAMINATIONS 2008

DEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING EXAMINATIONS 2008 [E5] IMPERIAL COLLEGE LONDON DEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING EXAMINATIONS 008 EEE/ISE PART II MEng BEng and ACGI SIGNALS AND LINEAR SYSTEMS Time allowed: :00 hours There are FOUR quesions

More information

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow 1 KEY Mah 4 Miderm I Fall 8 secions 1 and Insrucor: Sco Glasgow Please do NOT wrie on his eam. No credi will be given for such work. Raher wrie in a blue book, or on our own paper, preferabl engineering

More information

14 Autoregressive Moving Average Models

14 Autoregressive Moving Average Models 14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class

More information

Economics 8105 Macroeconomic Theory Recitation 6

Economics 8105 Macroeconomic Theory Recitation 6 Economics 8105 Macroeconomic Theory Reciaion 6 Conor Ryan Ocober 11h, 2016 Ouline: Opimal Taxaion wih Governmen Invesmen 1 Governmen Expendiure in Producion In hese noes we will examine a model in which

More information

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities: Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial

More information

Motion. Part 2: Constant Acceleration. Acceleration. October Lab Physics. Ms. Levine 1. Acceleration. Acceleration. Units for Acceleration.

Motion. Part 2: Constant Acceleration. Acceleration. October Lab Physics. Ms. Levine 1. Acceleration. Acceleration. Units for Acceleration. Moion Accelerion Pr : Consn Accelerion Accelerion Accelerion Accelerion is he re of chnge of velociy. = v - vo = Δv Δ ccelerion = = v - vo chnge of velociy elpsed ime Accelerion is vecor, lhough in one-dimensionl

More information

1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx.

1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx. . Use Simpson s rule wih n 4 o esimae an () +. Soluion: Since we are using 4 seps, 4 Thus we have [ ( ) f() + 4f + f() + 4f 3 [ + 4 4 6 5 + + 4 4 3 + ] 5 [ + 6 6 5 + + 6 3 + ]. 5. Our funcion is f() +.

More information

Signals and Systems Profs. Byron Yu and Pulkit Grover Fall Midterm 1 Solutions

Signals and Systems Profs. Byron Yu and Pulkit Grover Fall Midterm 1 Solutions 8-90 Signals and Sysems Profs. Byron Yu and Pulki Grover Fall 07 Miderm Soluions Name: Andrew ID: Problem Score Max 0 8 4 6 5 0 6 0 7 8 9 0 6 Toal 00 Miderm Soluions. (0 poins) Deermine wheher he following

More information

14.4. Lengths of curves and surfaces of revolution. Introduction. Prerequisites. Learning Outcomes

14.4. Lengths of curves and surfaces of revolution. Introduction. Prerequisites. Learning Outcomes Lengths of curves nd surfces of revolution 4.4 Introduction Integrtion cn be used to find the length of curve nd the re of the surfce generted when curve is rotted round n xis. In this section we stte

More information

For example, the comb filter generated from. ( ) has a transfer function. e ) has L notches at ω = (2k+1)π/L and L peaks at ω = 2π k/l,

For example, the comb filter generated from. ( ) has a transfer function. e ) has L notches at ω = (2k+1)π/L and L peaks at ω = 2π k/l, Comb Filers The simple filers discussed so far are characeried eiher by a single passband and/or a single sopband There are applicaions where filers wih muliple passbands and sopbands are required The

More information

Solutions from Chapter 9.1 and 9.2

Solutions from Chapter 9.1 and 9.2 Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is

More information

Vanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law

Vanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law Vanishing Viscosiy Mehod. There are anoher insrucive and perhaps more naural disconinuous soluions of he conservaion law (1 u +(q(u x 0, he so called vanishing viscosiy mehod. This mehod consiss in viewing

More information

Composite Reinforcement of Cylindrical Pressure Vessels

Composite Reinforcement of Cylindrical Pressure Vessels Comosie einorcemen o Cylindricl Pressure Vessels Cylindricl Pressure Vessels Cylindricl ressure vessels re in idesred use or vriey o licions SCBA nd SCUBA nks Prone nks Comressed Nurl Gs (CNG) nd ydrogen

More information

Green s Functions and Comparison Theorems for Differential Equations on Measure Chains

Green s Functions and Comparison Theorems for Differential Equations on Measure Chains Green s Funcions nd Comprison Theorems for Differenil Equions on Mesure Chins Lynn Erbe nd Alln Peerson Deprmen of Mhemics nd Sisics, Universiy of Nebrsk-Lincoln Lincoln,NE 68588-0323 lerbe@@mh.unl.edu

More information

1. Introduction. 1 b b

1. Introduction. 1 b b Journl of Mhemicl Inequliies Volume, Number 3 (007), 45 436 SOME IMPROVEMENTS OF GRÜSS TYPE INEQUALITY N. ELEZOVIĆ, LJ. MARANGUNIĆ AND J. PEČARIĆ (communiced b A. Čižmešij) Absrc. In his pper some inequliies

More information

Phys. 506 Electricity and Magnetism Winter 2004 Prof. G. Raithel Problem Set 2 Total 40 Points. 1. Problem Points

Phys. 506 Electricity and Magnetism Winter 2004 Prof. G. Raithel Problem Set 2 Total 40 Points. 1. Problem Points Phys. 56 Electricity nd Mgnetism Winter 4 Prof. G. ithel Problem Set Totl 4 Points 1. Problem 8.6 1 Points : TM mnp : ω mnp = 1 µɛ x mn + p π y with y = L where m, p =, 1,.. nd n = 1,,.. nd x mn is the

More information