Copyright 1979, by the author(s). All rights reserved.

Size: px
Start display at page:

Download "Copyright 1979, by the author(s). All rights reserved."

Transcription

1 Cpyright 1979, by the authr(s). All rights reserved. Permissin t make digital r hard cpies f all r part f this wrk fr persnal r classrm use is granted withut fee prvided that cpies are nt made r distributed fr prfit r cmmercial advantage and that cpies bear this ntice and the full citatin n the first page. T cpy therwise, t republish, t pst n servers r t redistribute t lists, requires prir specific permissin.

2 A FLUID MODEL FOR MULTIPLE-MIRROR AXIAL PLASMA FLOW PART I PRELIMINARIES by R. Bravenec Memrandum N. UCB/ERL M79/49 11 June 1979 ELECTRONICS RESEARCH LABORATORY Cllege f Engineering University f Califrnia, Berkeley 94720

3 A FLUID MODEL FOR MULTIPLE-MIRROR AXIAL PLASMA FLOW PART I PRELIMINARIES* R. Bravenec Department f Electrical Engineering and Cmputer Science Electrnics Research Labratry University f Califrnia Berkeley CA 9^720 ABSTRACT The fluid equatins are slved fr a straight magnetic field and a single mirrr under the assumptin f a steady-state, nedimensinal, isthermal flw. The latter assumptin is fund t be a pr ne during supersnic, shcked flw; therwise, it is accept able. Flw thrugh a magnetic mirrr is studied as a functin f the relative viscsity f the plasma and cmpared t fluid flw thrugh a V cnverging-diverging nzzle. Finally, applicatins t multiplemirrr flw are discussed. 'Spnsred by Department f Energy Grant N. DE-AS03-76F0003zf-PA# DE-ATOE-76ET

4 A FLUID MODEL FOR MULTIPLE-MIRROR AXIAL PLASMA FLOW PART I PRELIMINARIES I. INTRODUCTION This reprt lays the fundatin fr studying multiplemirrr axial plasma flw via a fluid mdel. The fluid equatins are slved fr a unifrm magnetic field and fr a cnverging-diverging field (a magnetic mirrr). The latter cnfiguratin is f interest because a multiple-mirrr device is a successin f single mirrrs. Our interest in a unifrm field is in studying the basic character istics f the equatins and their slutins. Frequently we will cmpare the plasma flw alng the mag netic field t actual fluid flw thrugh pipes f varying crss-sectin. The similarity is clear after making the fllwing assumptins: the plasma-vacuum interface is sharp and the plasma is f lw 3 > s that the magnetic field is cnsidered a knwn functin f z. Thrugh cnservatin f magnetic flux the plasma crss-sectin A is als knwn.1 Thus we may cmpare the plasma flw t that f a fluid flwing thrugh a pipe f crss-sectin A(z). Thrughut this reprt the flw is assumed steady, nedimensinal and isthermal. A steady-state results frm an arbitrary plasma surce at the system rigin and a sink at the exit. Radial mtin *Fr high 3, A(z) is nt s simply related t the vacuum magnetic field but must' be slved fr via pressure balance, as in Taylr and Wessn [1]

5 is neglected by making a "lng, thin apprximatin" with B large enugh t prevent any significant radial diffusin. Finally, isthermy riginates frm a large electrn thermal cnductivity which maintains Te(z)=T, a cnstant. This is cmmunicated t the ins via inelectrn cllisins n the time scale x. ^ x. *x» /m./m t.., where the le le i e ii superscript e dentes an ene/igy exchange time, and where t. and i e t.. are the in-electrn and in-in bcjouttojilyiq times, respectively. (See Trubnikv [2].) Thus, fr isthermal ins we require that x. «x, le s where xg is the time fr a fluid element t traverse the scale length f inhmgeneity f the fluid I. All ther electrn effects s are ignred in the analysis which fllws. With this assumptin the similarity between the plasma flw and that f a viscus fluid is even mre prnunced. In fact, the equatin f mtin used here differs frm the Navier-Stkes Equatin nly in the frm f the viscsity term. In a straight system they are identical. Sectin II" presents the fluid equatins and the derivatin f a secnd rder ODE fr the fluid velcity. This equatin is cast in dimensin less frm t select ut the imprtant parameters f the prblem and becmes the basis frmst f this reprt. Sectin III presents the analytical slutins f this equatin fr a unifrm magnet ic field. In Sec. IV an analytical slutin neglecting viscsity is btained fr a magnetic mirrr; hwever, this slutin des nt give transnic slutins (transitins frm subsnic t supersnic flw, r vice versa). The transnic case is discussed qualitatively, drawing heavily frm Shapir [3]. Finally, the mirrr-flw prblem, including - 3 -

6 viscsity, is slved numerically and the results cmpared t the inviscid case. Sectin V ummarizes the results and discusses the applica tins t multiple-mirrr flw. - h -

7 II. THE FLUID EQUATIONS The fllwing are taken directly frm Dawsn and Uman [k] ignring the electrstatic ptential and setting 3/3t = 0. ION EQUATION OF MOTION nm,v --&bti>** K'isf) +<-T«>i: <' > CONTINUITY EQUATION ^(nav) =0 (lb) PERPENDICULAR ION TEMPERATURE dti Ti-Tj 2vm /Hv/X2 T±f V a - +! L dz t.. \dz/ A dz + m i r7 [f<«>+t(t»e-t>] i ei "i 33f^[nA". dt{tl+2tj.,)] <lc> PARALLEL ION TEMPERATURE J dtl 2(Tl - Tjj ) 2v.m. dz x /dv\ \dl) ' 2T» i dv d7 1 m r«+-l--s. r^(te-t') - 5 -

8 + _la rna i(ti+2ti,] (ld) 3nA dz L i dz " J- J FLUX CONSERVATION B A = A^ -2- (ie) In the abve, n is the in density (assumed equal t the electrn density), A is the plasma crss-sectin, v is the z-cmpnent f the flw velcity (same fr ins and electrns), Tn ' and j» the electrn and in" temperatures (in energy units) i respectively paralare lei and perpendicular t the magnetic field, m. and m are the in and electrn masses, v. is the in bulk viscsity, x.. and x. are i 11 ei the in-in and electrn-in cllisin times, and k. is the in thermal i 2 cnductivity. The reader is referred t Dawsn and Uman fr a thrugh discussin f these equatins. Althugh we ultimately set T = Tn =T, Eqs. (1c) and (1d) are included t evaluate the last term f Eq. (1a). We cannt simply set T = Tn here since this term is fund t be f the same rder as the bulk viscsity term. Setting dt1 dz dtji' dz in Eqs. (1c) and (id), subtracting them frm each ther and then setting T - Tn - T thereafter (except in the difference T - Tn ), we btain 2 Thrughut this reprt all quantities subscripted with a zer are evalu ated at z = 0, i.e., A =A(0), B = B(0), etc

9 X1 I-' 1 T A d /V \ /0v T - Tn ~ x.. T ( J (2) 1 II 3 11 v dz \ /\ / Thus, the temperature difference is prprtinal t the small quantity x.., as is the in bulk viscsity v. = x.. T (3) 1 m. 11 v-" ( a f rder unity). Substituting Eqs. (2) and (3) int Eq. (1a), we have ydva _zl d(nt) } a _d_ dz m.n dz m.na dz \ dz/, 3,... m.v dz dz \ A / 1 1 x ' 1 (T..TnA^)+lIlil " d l \ ( Next we define the fllwing dimensinless quantities "af- vi7" a5f <«sf- M5^= v (5) c s Here i^ is an apprpriate reference length in the prblem, i.e., the cell length f a multiple mirrr, and v = /T/m. is the sund speed s 1 r (als the in thermal speed). Substituting int Eq. (4) and using the 3/2 relatin x..=x/vs«t /n, where X is the in-in mean-free-path, we btain v±l=z±l*r +?lj>_±j_ I dv\.1^ _L J_ da d /v2\ tc, d? M2 n d5 c Mq na d? TdJ 3 *c Mq nv dc dc IT7 (6a) The cntinuity equatin becmes simply nav = 1 (6b) - 7 -

10 Using Eq. (6b) t eliminate n, Eq. (6a) becmes dv_j_j_ d(av) = ^J_ dc M2 av dc I M M_ c av d / dv\ a_ da d_ /v d? \adj 3d;d^ay (7) Equatin (7) is a secnd-rder ODE fr the nrmalized flw velcity v in terms f the knwn area functin a(?) and the tw parameters X / and M, hereafter referred t as the "flw equatin". Its c M slutin is cmpletely determined given the tw bundary cnditins v(0) and dv/d?(0) (r v(0) and v(l/& )), where L is the system c length). The first f each set is trivial since by definitin v(0) =1 (The bundary cnditin v(0) makes its appearance thrugh the para meter M.) Recall that the terms n the left-hand side represent the in inertia and scalar pressure, respectively. The term n the right represents the cmbined effects f pressure anistrpy and will here after be referred t as simply the "viscsity term." Fr v^dv/de; 2 2 ^d v/d? ^0(1), the relative rder f each term is easily seen t be inertia : pressure : viscsity = M ' : 1 : -t M & c Thus fr M,X /I «1 the pressure term is dminant. Hwever, fr c r ' v=1/m s that M=1, the rdering becmes 1 : 1 : (X / )(l/m ) s c that all terms may be cmparable. Fr future reference it is useful t rewrite the cnditin fr unifrm in temperature in terms f ur dimensinless variables. First it can be shwn that x = (l/m) U /% ) (z /x)x.. Since s seen A«1/n this can be written x = (l/m) U A ) U /l ) (x../n). Hence, S C O S C - 8 -

11 the cnditin fr isthermal ins becmes x A u O c M. i ^ /n\ Finally, we shuld nte that Eq. (4) differs frm that used by Miller [5], wh derives the frce equatin directly frm Braginskii [6]. The differences, which as f yet have nt been reslved, lie in the explicit frm f the viscsity term. Since this term is usually small, we will nt verly cncern urselves with this discrep ancy. J - 9 -

12 III. FLOW ALONG A UNIFORM MAGNETIC FIELD Setting a = 1 in Eq. (7) we have X 2 dv _]_ j_ dv av d v, * d? " M2 vd? I M.2 {3) M " cd? Ntice that the temperature anistrpy term has drpped ut, leaving us with the Navier-Stkes Equatin. Dividing by v and integrating nce, there results -' i(h - rr( -) M c where *;s < > Slving fr dv/d?, we btain. ir (, i (1.,\ +v. (10) dc A a I ' M \v /J This equatin can be integrated by first separating variables. Fr v«?*0 vdv M 2 /m ^ 1 X0 A.. 1 Cd?,,-x Mv - I M +TT -r~av,lv+-7- = r 2- (11) \ M j M X. a c /

13 The results f the integratin depend n whether M <1 r M >1. Fr v'l K -'): «_ _L M (M +1) ' X 2a..,..2 the expressins are particularly simple3 a *c 1-M2 1/M' - v 1. \ 1/M - 1 '1 + u' - v> H In ' M < 1 (12a) C = a *cm2-1 f] + y» - v^ * 1 v- 1/M' jr M ln\ Vl-1/M2 2 O J, M > 1 (12b) where ^ 5 c am M (12c) There are tw distinct prfiles given by each f Eqs. (12), depending n the sign f V. Fr v^<0, Eq. (12a) gives the same prfile as (12b) except shifted in c. Chsing v' >0 in Eq. (12b) leads t an entirely unphysical prfile in which the velcity increases withut bund Thus, we are left with the tw cases M <1, v' > 0 and M >1, v1 <0 3 This inequality is impsed nly t btain the simple frms f Eqs. (12), and nly restricts ur placement f the rigin relative t the prfiles

14 As an example f the frmer, cnsider the chices M =0.5, v^o.01, a=1 fr which Eq. (12a) gives the prfiles shwn in Fig. 1(a) Here X /& is set t ne s that the units f c are mean-free-paths. c ^ r Superimpsed n the curve M = M v is that f n, btained frm the cn tinuity relatin nv = 1. The mst bvius feature f the curves is that M-*0 and n^00 at a pint c. t the left f the rigin, where mm c. = r In mm M (1-M2) y 1-M2 Althugh we have, impsed an upper limit n v" t btain Eqs. (12), slutins certainly exist fr v1 large. Fr example, shifting the curves f Fig. 1(a) t the right by an amunt less than \c,. rep resents such a slutin. The resulting flw is then characterized by an acceleratin and pressure drp within a few mean-free-paths int the system. Whatever the surce f the initial velcity gradient, the impr tant pint here is that bulk viscsity des nt lead t an increasingly accelerating flw, as is the case during fluid flw thrugh pipes with wall frictin [3]. Instead, it simply smths ut the initial gradients, causing the prfiles t asymptte t steady values. Nte that fr v1 =0 M arbitrary, Eq. (9) has the trivial slutin n = v=1, which may be easily seen frm Eq. (10). We have just shwn that the equatins have a preference fr such a unifrm slutin; thus we cnclude that the imp sitin f v* >0 des nt represent an easily realizable physical system. This cnclusin will be applied in the next sectin n flw thrugh a mag netic mi rrr

15 An example f the case M >1, v1 <0 is presented in Fig. 1(b) fr M =2.0, v»=-0.01, a=1. As required by basic fluid thery [3], the velcity decreases and the density increases during vis cus supersnic flw. These prfiles represent a shck frnt, resulting in a transitin t subsnic flw ( M = 0.5 ) ver a distance f abut fur mean-free-paths. Referring back t Eq. (12b), we see that in general the tw asympttes f v(c) are v - 1 and v. =1/M, such that max mm 2 the shck strength is v /v. =M. In a frame with X /% «1, Eq. max m in c (12b) predicts an almst discntinuus prfile. This is expected since the shck thickness depends n the magnitude f the viscsity term which in turn is prprtinal t X /I. Furthermre, n either side f the c shck dv/dc is very small. Thus, fr a fixed value f M, varying v1 within the range *, (M2-1)2 0 < v' < C 1 M (M + 1) has the effect f shifting the shck alng c (Letting \>' -+Q drives the shck tward +«.) This behavir will reappear in the fllwing sectin n mirrr flw. Befre cntinuing, we shuld realize that the previus results are nt t be interpreted t literally. First f all we have n justifi catin in accepting slutins which vary significantly ver a few meanfree-paths. This is a direct cnsequence f the limitatins f the fluid thery. Thus, we shuld be wary f the "bundary-layer effect" predicted by Eq. (12a) r the details f the shck structure f Eq. (12b). As t

16 the latter, hwever, we are justified in accepting the values attained by the prfiles dwnstream f the shck. This is because the equatins cnserve such quantities as mass and mmentum flux [3]. That is, the shck strength predicted by Eq. (12b) is crrect under the assumptin f isthermy. Hwever, as mentined in the Intrductin, this assumptin is valid nly fr T (z)=t, a cnstant, and x. «x. Referring t e le s 3 Eq. (8) with I ^X, M^ 1 and n^ 1, the latter cnditin is badly vilated and isthermy is a pr assumptin. We will nt let this deter us, hwever, since ur ultimate intentin is nt t study shck frnts but subsnic, shck-free flw thrugn a multiple-mirrr, fr which is thermy is a fair assumptin. If desired we may imagine the temperature maintained unifrm acrss the shck by sme external means. 4?

17 IV. FLOW THROUGH A MAGNETIC MIRROR Fr this sectin and the remainder f this reprt we use the area functin (C) = ^f+^fcs 2w5 (13) where Rs B /B., the mirrr rati. This functin is sketched fr nia j\ m i n R= 2 in Fig. 2. Substituting Eq. (13) int Eq. (7) results in an expres sin that can be slved easily nly in the limit X /l =0. In this c case, Eq. (7) becmes dv _ J_ d In av Slving fr dv/d?, dv 1 da d? 3dC(Mv)2-1 d4) This can be integrated t yield 2 v 2 - in av» 1 (15a) M Using Eq. (6b) we can write this in terms f n a 2'n2(l-^lnn) =1 (15b)

18 Eqs. (15) are sketched in Fig. 3 fr R=2 and the tw chices (a) Mq =0.25 and (b) Mq =2.0. As expected the subsnic slutin is char acterized by a density minimum and a velcity maximum in the mirrr thrat, and the supersnic slutin by just the ppsite. Ntice further that fr a(c) peridic with perid unity the slutins are similarly peridic. This means that there is n verall density (pressure) drp acrss a successin f such mirrrs. We will return t this pint in the final sectin. Referring t Eq. (14), the slutin becmes singular at v=1/m (M=1), preventing a transnic slutin. We nw discuss qualitatively the family f slutins, bth subsnic and transnic, fr an inviscid fluid flwing thrugh a cnverg ing-diverging nzzle. We imagine being given the inlet pressure p and investigating the flw as a functin f the exit pressure p. Refer ences are repeatedly made t Fig. 4, taken frm Shapir, p. 140, which is fr adlabatic flw f a perfect gas. Initially, entirely subsnic flw results when p1 is reduced belw p, as in Fig. 3(a). This is shwn by case (a) f Fig. 4. Lwering p1 further, state (b) is even tually reached, where the velcity in the thrat has reached M=1 (snic flw), but the flw dwnstream is still subsnic. The flw is nw chked, meaning that n further reductin in p affects the flw upstream f the thrat r increases the fluid flux. Between cases (b) and (f) there exist n values f p. that result in smth shck-free flw. Case (f) crrespnds t a unique value f p1 fr which shck-free supersnic flw des exist in the diverging sectin. Cases (c) and (d) represent shck frnts invlving abrupt transitins frm supersnic t subsnic flw. As P1 is reduced frm that f case (b) the shck is seen' t prgressively

19 mve dwn the nzzle. The flw appraches that f case (d), in which the shck ccurs at the exit. Except fr case (f), further reductin f the exit pressure results in blique shcks utside the nzzle, as represented by cases (e) and (g). These effects are tw-dimensinal and cannt be treated in ur simple thery. The flw characterized in Fig. 4 is that f an adiabatic, inviscid fluid. We wish nw t investigate the effects f isthermy and viscsity. Fr X /% «1 we expect Fig. 4 t be qualitatively duplicated since the viscsity term f Eq. (7) is then small. The slight amunt f viscsity simply gives structure t the shck frnts (Fig. 1(b)) s that n singularities ccur in the flw equatin. Fur thermre we expect that ur impsitin f cnstant temperature nly affects the magnitude f variatin f v and n, but nt the general behavir characterized in Fig. 4. The flw equatin (7) is slved numerically using an inte grating rutine devised by Gear and Hindmarsh [7]. The prgram requires a system f first-rder ODE's necessitating the definitin f a new va r iab 1e Eq. (7) then becmes w >L = Ida. [-LL*±/2\v^] 1 c 1 I"/,. t 1 \, 1 idal dc adcvl3aadc V+3a j VJ +a X~ iv LlV"FTVJ V mt I di"j (16b)

20 In the fllwing, a is set t 0.96 [6]. The rutine further requires all bundary cnditins t be given at C= 0, namely v(0) and v'(0). (As mentined earlier, we have defined v(0) =1 s that v(0) appears thrugh the parameter M^.) Thus, given a value fr X /, the slu " c tin t the system f Eqs. (16) is cmpletely determined by the set M and v' h v' (0). T generate curves similar t thse f Fig. 4, X /% is c set t 0.01, s that the fluid is nearly inviscid. Frm the findings f the last sectin, v1 is set t zer since the plasma initially sees 1-7 a straight system da/d^(0)=0. Tw cases f entirely subsnic flw are presented in Fig. 5 fr M =0.25 and M =0.3. We immediately nte that the small amunt f viscsity shifts the velcity maximum and density minimum int the divergent sectin f the mirrr. As M is in creased t duplicate case (b) f Fig. 4, a value is reached fr which snic flw first ccurs at a pint slightly beynd the mirrr thrat. This situ atin is represented by the case labeled M = As M is increased further, the flw becmes supersnic ver a shrt distance with a subse quent smth transitin back t subsnic flw. This behavir is inter preted as a weak shck develping in the divergent sectin f the mirrr. As Mq is increased still further, appraching snic flw in the thrat, the shck frnt grws in amplitude, steepens, and mves dwnstream. Curves fr which M = 1 in the thrat are presented in Fig. 5, crrespnding k The value f M in the mirrr thrat is simply estimated frm the curves. Hence, there is a small uncertainty invlved in the actual value

21 t M = An interesting phenmenn ccurs if M is chsen s just slightly larger s that M is barely greater than unity in the thrat the slutin fr v diverges, indicating an unphysical slu tin. An example is the case labeled M = Thus, the numerical slutins preserve the phenmenn f chking n physical steady-state slutins fr M larger than a critical value. Furthermre, it appears that chking ccurs fr M=1 in the thrat even during isthermal vis cus flw. In rder t prceed further and generate curves analgus t cases (b) thrugh (f) f Fig. 4, we begin the integratin at C= 0.5 (in the mirrr thrat) where we may accurately impse the cnditin M=1. This leaves nly the quantity v'(.5) with which t determine the flw. As argued in the inviscid fluid case, the flw upstream f the thrat is unaffected by reducing the exit pressure. Thus, in the fllwing cases the prfiles upstream f the thrat are thse fr critical flw. It shuld be recalled that the relevant physical parameter that is being varied here is the exit pressure (density) via the bundary cnditin v'(.5). Thus, nly thse values f v'(.5) giving n(1.0)<n. (1.0) are cr 11 valid. This cnditin translates t v'(.5)>v'. (.5). It is fund en t that there exists a maximum value f v'(.5) abve which the slutin diverges. This valve, which we will dente by v' (.5), crrespnds t ' max ' v the instance f smth, shck-free supersnic flw, as in case (f) f Fig. 4. Furthermre vmax(-5) is fund t be nly slightly larger than v' :«.(*5), leaving us nly the restricted range v1. (.5)<v'(,5)<v' (.5) crit crit - - max This case is hereafter referred t as "critical flw". In the fllwing discussin its slutins n and v are subscripted with "crit"

22 with which t wrk. In rder t generate a shck frnt and mve it dwn stream, v'(.5) must be chsen very near v1 (.5). T accmplish this, max r ' the prgram was mdified t search fr that value f v'(.5) which yields a specified value f n(1.0), the physical bundary cnditin. The results fr three chices f n(1.0) are presented in Fig. 5, labeled n(l.o) =0.85, 0.65, and 0.2. As evident frm the curves, the differ ences between their values f v1 (.5) are unnticable. The case n(1.0) =0.20 crrespnds t all but the "tail" f the shck having mved ut the system. Thus, we have seen that the bundary cnditin n dv/dc is nly imprtant in generating the flw prfiles dwnstream f the thrat in a chked system. Otherwise M has by far the dminant rle. We have als seen that the effect f v'(.5) is the same as in ur previus study f shcks in a straight field that is, it shifts the lcatin f the shck. Finally, we shuld mentin that isthermy is a fair assumptin nly fr the case labeled M =0.25. That this is s can be seen using 3 Eq. (8) with X /A =0.01, a / a (frm the curve fr M ), M*0.5 and n^o^ (wrst cases). Fr larger values f M, the scale lengths fr density and velcity variatins are t shrt fr in-electrn clli sins t maintain T. a cnstant. 1 T cmplete ur study f viscus flw thrugh a single mirrr, we nw cnsider the flw f a mre viscus plasma. T d this we increase vj the value f X / in Eqs. (7) and (16b). Hwever, we are limited in the fluid thery t X «Z, where as befre l is a typical scale s s length fr variatins in n, v r A. Analgus t Eq. (8) this cnditin can be written in terms f ur dimensin less variables as

23 Fr M>j/nTTrnT, this cnditin is less stringent than that f Eq. (8). T begin, we define I in terms f A(z) fr which I /% «\. Fr s s c n=1, Eq. (17) then demands that X /l «. Chsing X /I =0.1 c c and vq =0, we present tw cases f purely subsnic flw in Fig. 6 as the curves labeled M =0.15 and M =0.25. We make the fllwing cm 3 parisns t Fig. 5: (1) the pints f velcity maxima and density minima are displaced farther dwnstream f the thrat. (2) Fr the same value f Mq (fr example, Mq =0.25 ) the pressure (density) drp acrss the system is greater. As Mq is increased, M=1 first ccurs at 5=0.675, pre sented by the curves labeled M = Increasing M further results in a regin f supersnic flw fllwed by a very slw transitin back t subsnic flw, such as the case labeled M = As M is in creased still further, the regin f supersnic flw increases until fr MQ =0.331 the flw dwnstream f the thrat is entirely supersnic. Thus, we see that shck-like behavir vanishes in such a viscus plasma. As Mq is increased further, the psitin f the velcity maximum and density minimum mves ut the exit. The curves fr M* 1 in the mirrr thrat & are given by thse labeled M = Unlike the case X /% =0.01, c n catastrphic behavir ccurs fr M>1 in the thrat. Referring back t Fig. 5, we ntice that fr M<1 in the thrat the velcity prfiles exhibit inflectin pints near the thrat. Hwever, fr M>1 jn tne Fr M=1 in the thrat it can be shwn using Eq. (16b) that this psitin is the inflectin pint

24 thrat, n such pint exists and the prfiles diverge. In the present case, hwever, the inflectin pint mves relatively slwly tward the exit as M is made greater than unity in the thrat. Tw examples are the curves labeled M =0.36 and M =0.4. In this sense the numerical slutins n lnger demnstrate an bvius chking effect. We instead resrt t intuitin and a cursry study f the flw equatin in asserting that chking still ccurs fr M=1 in the thrat, even with substan tial viscsity, and that the prfiles with M>1 here are unphysical. Since shck frnts d nt ccur there is n need t initialize the inte gratins in the thrat, where we specify M=1 and vary v'(.5). This is because we are able t generate all exist densities (pressures) dwn t that analgus t case (f) f Fig. 4 via the parameter M0 alne. Althugh the curves are much smther than thse f Fig. 5, isthermy is fund t be a pr assumptin fr all the cases shwn in Fig. 6. Even fr M =0.15 where %/% =0.6, M =0.3 and n-0.9, s c the inequality (8) is badly vilated. This is because X /I is a fac c tr f ten larger than that f Fig. 5, resulting in relatively few inelectrn cllisins. Finally we check if Eq. (17) is satisfied with &s/&c defined in terms f n and v, and taking int accunt the substantial decay f n alng? (which lengthens the mean-free-path). Because W& «s large where n is small, Eq. (17) is still fund t hld (althugh barely)

25 V. CONCLUSIONS Cmparing Figs. 4 and 5, we see that the flw f a nearly inviscid plasma thrugh a magnetic mirrr is very similar t fluid flw thrugh a cnverging-diverging nzzle. Hwever, the flw f a mderately viscus plasma is quite different, in that shck frnts d nt develp. The fact that isthermy was a pr assumptin des nt hinder us in applying the results t subsnic multiple-mirrr flw, in which shcks are excluded and the flw is usually slwer. T see hw this is dne, cnsider a "half-system" with z = 0 crrespnding t the midpint f the multiple-mirrr and z = L an exit. In the limit X /% =0, Eqs. c ^ (15) give the slutin fr a(&) peridic With a(c) given by Eq. (13), Fig. 3(a) represents an example f the prfiles at any particular mirrr f the device. As mentined at that time, there is n change in the ver all levels f density and velcity frm mirrr t mirrr. Hwever, the inclusin f viscsity results in an verall decrease in density (pressure) and an acceleratin f the flw frm mirrr t mirrr, as evident frm the entirely subsnic prfiles f Figs. 5 and 6. This effect is mre dramatic with greater viscsity. Fr small pressure differences between the ends f ur half-system, the flw will be as described. As the pres sure difference is increased, snic flw and subsequent chking will first ccur in the exit mirrr. In physical systems in which the exit pres V? sure is lw, a chked flw is a reasnable assumptin, and M=1 in the exit mirrr thrat becmes a "natural" bundary cnditin. Since we are nt cncerned with the flw beynd the exit mirrr thrat, shck^ frnts d nt ccur in the multiple-mirrr. Furthermre, the flw is everywhere subsnic, adding t ur justificatin in assuming isthermy. These bservatins will be made mre quantitative in a frthcming reprt

26 REFERENCES [1] Taylr, J. B., Wessn, J. A., Nuclear Fusin 5. (1965) 159. [2] Trubnikv, B. A., "Particle Interactins in a Fully Inized Plasma", in Reviews f Plasma Physics, Vl. 1, Cnsultants Bureau, New Yrk, [3] Shapir, A. H., The Dynamics and Thermdynamics f Cm pressible Fluid Flw, Vl. 1, The Rnald Press, New Yrk, [4] Dawsn, J. M., Uman, M. F., Nuclear Fusin 5. (1965) 242. [5] Miller, G., Ls Alams Scientific Labratry Reprt LA MS, Dec [6] Braginskii, S. I., "Transprt Prcesses in a Plasma", in Reviews f Plasma Physics, Vl. 1, Cnsultants Bureau, New Yrk, [7] Hindmarsh, A. C, "Gear-Ordinary Differential Equatin System Slver", UCID Rev. 3, Lawrence Livermre Labr atry, Dec

27 - 3 i M 2.0 M V H I k min (a) 0.0 M (b) Fig. 1 Viscus flw alng a unifrm magnetic field ( c = z/x ) (a) Subsnic flw ( M =0.5 ) (b) Supersnic flw ( M =2.0 )

28 0.0 F'9' 2 The nrmalized plasma crss-sectin vs =z/a

29 0.0 Fig. 3 Inviscid flw thrugh a magnetic mirrr ( X /I =0 ) (a) Subsnic flw ( M =0.25 ) (b) Supersnic flw ( M =2.0 )

30 (cnstant) Lcus f states dwnstream f nrmal shck Distance alng nzzle VJ Fig. 4 Fluid flw thrugh a cnverging-diverging nzzle

31 M<1 M>1 2.0 Y M>1 M > M<1 0.0 Fig. 5 Viscus flw thrugh a magnetic mirrr ( X /I =0.01 ) c (a) nrmalized density vs =z/l (b) Mach number vs z, = z/l

32 M< M> y M>l M Y M<l 0.0 Fig. 6 Viscus flw thrugh a magnetic mirrr ( X_/A_ =0.1 ) (a) nrmalized density vs c=z/% c (b) Mach number vs = z/jl c

Pressure And Entropy Variations Across The Weak Shock Wave Due To Viscosity Effects

Pressure And Entropy Variations Across The Weak Shock Wave Due To Viscosity Effects Pressure And Entrpy Variatins Acrss The Weak Shck Wave Due T Viscsity Effects OSTAFA A. A. AHOUD Department f athematics Faculty f Science Benha University 13518 Benha EGYPT Abstract:-The nnlinear differential

More information

Examiner: Dr. Mohamed Elsharnoby Time: 180 min. Attempt all the following questions Solve the following five questions, and assume any missing data

Examiner: Dr. Mohamed Elsharnoby Time: 180 min. Attempt all the following questions Solve the following five questions, and assume any missing data Benha University Cllege f Engineering at Banha Department f Mechanical Eng. First Year Mechanical Subject : Fluid Mechanics M111 Date:4/5/016 Questins Fr Final Crrective Examinatin Examiner: Dr. Mhamed

More information

Aircraft Performance - Drag

Aircraft Performance - Drag Aircraft Perfrmance - Drag Classificatin f Drag Ntes: Drag Frce and Drag Cefficient Drag is the enemy f flight and its cst. One f the primary functins f aerdynamicists and aircraft designers is t reduce

More information

Study Group Report: Plate-fin Heat Exchangers: AEA Technology

Study Group Report: Plate-fin Heat Exchangers: AEA Technology Study Grup Reprt: Plate-fin Heat Exchangers: AEA Technlgy The prblem under study cncerned the apparent discrepancy between a series f experiments using a plate fin heat exchanger and the classical thery

More information

(1.1) V which contains charges. If a charge density ρ, is defined as the limit of the ratio of the charge contained. 0, and if a force density f

(1.1) V which contains charges. If a charge density ρ, is defined as the limit of the ratio of the charge contained. 0, and if a force density f 1.0 Review f Electrmagnetic Field Thery Selected aspects f electrmagnetic thery are reviewed in this sectin, with emphasis n cncepts which are useful in understanding magnet design. Detailed, rigrus treatments

More information

ELECTRON CYCLOTRON HEATING OF AN ANISOTROPIC PLASMA. December 4, PLP No. 322

ELECTRON CYCLOTRON HEATING OF AN ANISOTROPIC PLASMA. December 4, PLP No. 322 ELECTRON CYCLOTRON HEATING OF AN ANISOTROPIC PLASMA by J. C. SPROTT December 4, 1969 PLP N. 3 These PLP Reprts are infrmal and preliminary and as such may cntain errrs nt yet eliminated. They are fr private

More information

A Few Basic Facts About Isothermal Mass Transfer in a Binary Mixture

A Few Basic Facts About Isothermal Mass Transfer in a Binary Mixture Few asic Facts but Isthermal Mass Transfer in a inary Miture David Keffer Department f Chemical Engineering University f Tennessee first begun: pril 22, 2004 last updated: January 13, 2006 dkeffer@utk.edu

More information

Compressibility Effects

Compressibility Effects Definitin f Cmpressibility All real substances are cmpressible t sme greater r lesser extent; that is, when yu squeeze r press n them, their density will change The amunt by which a substance can be cmpressed

More information

Flipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System

Flipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System Flipping Physics Lecture Ntes: Simple Harmnic Mtin Intrductin via a Hrizntal Mass-Spring System A Hrizntal Mass-Spring System is where a mass is attached t a spring, riented hrizntally, and then placed

More information

Flipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System

Flipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System Flipping Physics Lecture Ntes: Simple Harmnic Mtin Intrductin via a Hrizntal Mass-Spring System A Hrizntal Mass-Spring System is where a mass is attached t a spring, riented hrizntally, and then placed

More information

7.0 Heat Transfer in an External Laminar Boundary Layer

7.0 Heat Transfer in an External Laminar Boundary Layer 7.0 Heat ransfer in an Eternal Laminar Bundary Layer 7. Intrductin In this chapter, we will assume: ) hat the fluid prperties are cnstant and unaffected by temperature variatins. ) he thermal & mmentum

More information

CHAPTER 3 INEQUALITIES. Copyright -The Institute of Chartered Accountants of India

CHAPTER 3 INEQUALITIES. Copyright -The Institute of Chartered Accountants of India CHAPTER 3 INEQUALITIES Cpyright -The Institute f Chartered Accuntants f India INEQUALITIES LEARNING OBJECTIVES One f the widely used decisin making prblems, nwadays, is t decide n the ptimal mix f scarce

More information

, which yields. where z1. and z2

, which yields. where z1. and z2 The Gaussian r Nrmal PDF, Page 1 The Gaussian r Nrmal Prbability Density Functin Authr: Jhn M Cimbala, Penn State University Latest revisin: 11 September 13 The Gaussian r Nrmal Prbability Density Functin

More information

Lecture 5: Equilibrium and Oscillations

Lecture 5: Equilibrium and Oscillations Lecture 5: Equilibrium and Oscillatins Energy and Mtin Last time, we fund that fr a system with energy cnserved, v = ± E U m ( ) ( ) One result we see immediately is that there is n slutin fr velcity if

More information

Thermodynamics Partial Outline of Topics

Thermodynamics Partial Outline of Topics Thermdynamics Partial Outline f Tpics I. The secnd law f thermdynamics addresses the issue f spntaneity and invlves a functin called entrpy (S): If a prcess is spntaneus, then Suniverse > 0 (2 nd Law!)

More information

February 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA

February 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA February 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA Mental Experiment regarding 1D randm walk Cnsider a cntainer f gas in thermal

More information

SPH3U1 Lesson 06 Kinematics

SPH3U1 Lesson 06 Kinematics PROJECTILE MOTION LEARNING GOALS Students will: Describe the mtin f an bject thrwn at arbitrary angles thrugh the air. Describe the hrizntal and vertical mtins f a prjectile. Slve prjectile mtin prblems.

More information

ChE 471: LECTURE 4 Fall 2003

ChE 471: LECTURE 4 Fall 2003 ChE 47: LECTURE 4 Fall 003 IDEL RECTORS One f the key gals f chemical reactin engineering is t quantify the relatinship between prductin rate, reactr size, reactin kinetics and selected perating cnditins.

More information

ENGI 4430 Parametric Vector Functions Page 2-01

ENGI 4430 Parametric Vector Functions Page 2-01 ENGI 4430 Parametric Vectr Functins Page -01. Parametric Vectr Functins (cntinued) Any nn-zer vectr r can be decmpsed int its magnitude r and its directin: r rrˆ, where r r 0 Tangent Vectr: dx dy dz dr

More information

Module 4: General Formulation of Electric Circuit Theory

Module 4: General Formulation of Electric Circuit Theory Mdule 4: General Frmulatin f Electric Circuit Thery 4. General Frmulatin f Electric Circuit Thery All electrmagnetic phenmena are described at a fundamental level by Maxwell's equatins and the assciated

More information

Physics 2010 Motion with Constant Acceleration Experiment 1

Physics 2010 Motion with Constant Acceleration Experiment 1 . Physics 00 Mtin with Cnstant Acceleratin Experiment In this lab, we will study the mtin f a glider as it accelerates dwnhill n a tilted air track. The glider is supprted ver the air track by a cushin

More information

Chapter 2 GAUSS LAW Recommended Problems:

Chapter 2 GAUSS LAW Recommended Problems: Chapter GAUSS LAW Recmmended Prblems: 1,4,5,6,7,9,11,13,15,18,19,1,7,9,31,35,37,39,41,43,45,47,49,51,55,57,61,6,69. LCTRIC FLUX lectric flux is a measure f the number f electric filed lines penetrating

More information

Electric Current and Resistance

Electric Current and Resistance Electric Current and Resistance Electric Current Electric current is the rate f flw f charge thrugh sme regin f space The SI unit f current is the ampere (A) 1 A = 1 C / s The symbl fr electric current

More information

ENGINEERING COUNCIL CERTIFICATE LEVEL THERMODYNAMIC, FLUID AND PROCESS ENGINEERING C106 TUTORIAL 5 THE VISCOUS NATURE OF FLUIDS

ENGINEERING COUNCIL CERTIFICATE LEVEL THERMODYNAMIC, FLUID AND PROCESS ENGINEERING C106 TUTORIAL 5 THE VISCOUS NATURE OF FLUIDS ENGINEERING COUNCIL CERTIFICATE LEVEL THERMODYNAMIC, FLUID AND PROCESS ENGINEERING C106 TUTORIAL 5 THE VISCOUS NATURE OF FLUIDS On cmpletin f this tutrial yu shuld be able t d the fllwing. Define viscsity

More information

Kinetics of Particles. Chapter 3

Kinetics of Particles. Chapter 3 Kinetics f Particles Chapter 3 1 Kinetics f Particles It is the study f the relatins existing between the frces acting n bdy, the mass f the bdy, and the mtin f the bdy. It is the study f the relatin between

More information

Lecture 7: Damped and Driven Oscillations

Lecture 7: Damped and Driven Oscillations Lecture 7: Damped and Driven Oscillatins Last time, we fund fr underdamped scillatrs: βt x t = e A1 + A csω1t + i A1 A sinω1t A 1 and A are cmplex numbers, but ur answer must be real Implies that A 1 and

More information

AP Physics Kinematic Wrap Up

AP Physics Kinematic Wrap Up AP Physics Kinematic Wrap Up S what d yu need t knw abut this mtin in tw-dimensin stuff t get a gd scre n the ld AP Physics Test? First ff, here are the equatins that yu ll have t wrk with: v v at x x

More information

Surface and Contact Stress

Surface and Contact Stress Surface and Cntact Stress The cncept f the frce is fundamental t mechanics and many imprtant prblems can be cast in terms f frces nly, fr example the prblems cnsidered in Chapter. Hwever, mre sphisticated

More information

Differentiation Applications 1: Related Rates

Differentiation Applications 1: Related Rates Differentiatin Applicatins 1: Related Rates 151 Differentiatin Applicatins 1: Related Rates Mdel 1: Sliding Ladder 10 ladder y 10 ladder 10 ladder A 10 ft ladder is leaning against a wall when the bttm

More information

General Chemistry II, Unit I: Study Guide (part I)

General Chemistry II, Unit I: Study Guide (part I) 1 General Chemistry II, Unit I: Study Guide (part I) CDS Chapter 14: Physical Prperties f Gases Observatin 1: Pressure- Vlume Measurements n Gases The spring f air is measured as pressure, defined as the

More information

NGSS High School Physics Domain Model

NGSS High School Physics Domain Model NGSS High Schl Physics Dmain Mdel Mtin and Stability: Frces and Interactins HS-PS2-1: Students will be able t analyze data t supprt the claim that Newtn s secnd law f mtin describes the mathematical relatinship

More information

Revision: August 19, E Main Suite D Pullman, WA (509) Voice and Fax

Revision: August 19, E Main Suite D Pullman, WA (509) Voice and Fax .7.4: Direct frequency dmain circuit analysis Revisin: August 9, 00 5 E Main Suite D Pullman, WA 9963 (509) 334 6306 ice and Fax Overview n chapter.7., we determined the steadystate respnse f electrical

More information

Rigid Body Dynamics (continued)

Rigid Body Dynamics (continued) Last time: Rigid dy Dynamics (cntinued) Discussin f pint mass, rigid bdy as useful abstractins f reality Many-particle apprach t rigid bdy mdeling: Newtn s Secnd Law, Euler s Law Cntinuus bdy apprach t

More information

Unit code: H/ QCF level: 5 Credit value: 15 OUTCOME 3 - STATIC AND DYNAMIC FLUID SYSTEMS TUTORIAL 3 - VISCOSITY

Unit code: H/ QCF level: 5 Credit value: 15 OUTCOME 3 - STATIC AND DYNAMIC FLUID SYSTEMS TUTORIAL 3 - VISCOSITY Unit 43: Plant and Prcess Principles Unit cde: H/601 44 QCF level: 5 Credit value: 15 OUTCOME 3 - STATIC AND DYNAMIC FLUID SYSTEMS TUTORIAL 3 - VISCOSITY 3 Understand static and namic fluid systems with

More information

Advanced Heat and Mass Transfer by Amir Faghri, Yuwen Zhang, and John R. Howell

Advanced Heat and Mass Transfer by Amir Faghri, Yuwen Zhang, and John R. Howell 6.5 Natural Cnvectin in Enclsures Enclsures are finite spaces bunded by walls and filled with fluid. Natural cnvectin in enclsures, als knwn as internal cnvectin, takes place in rms and buildings, furnaces,

More information

MODULE 1. e x + c. [You can t separate a demominator, but you can divide a single denominator into each numerator term] a + b a(a + b)+1 = a + b

MODULE 1. e x + c. [You can t separate a demominator, but you can divide a single denominator into each numerator term] a + b a(a + b)+1 = a + b . REVIEW OF SOME BASIC ALGEBRA MODULE () Slving Equatins Yu shuld be able t slve fr x: a + b = c a d + e x + c and get x = e(ba +) b(c a) d(ba +) c Cmmn mistakes and strategies:. a b + c a b + a c, but

More information

Interference is when two (or more) sets of waves meet and combine to produce a new pattern.

Interference is when two (or more) sets of waves meet and combine to produce a new pattern. Interference Interference is when tw (r mre) sets f waves meet and cmbine t prduce a new pattern. This pattern can vary depending n the riginal wave directin, wavelength, amplitude, etc. The tw mst extreme

More information

We need to do review for some thermodynamic relations: Equation of state P=ρ R T, h= u + pv = u + RT dh= du +R dt Cp dt=cv dt + R dt.

We need to do review for some thermodynamic relations: Equation of state P=ρ R T, h= u + pv = u + RT dh= du +R dt Cp dt=cv dt + R dt. Cmpressible Flw Cmpressible flw is the study f fluids flwing at speeds cmparable t the lcal speed f sund. his ccurs when fluid speeds are abut 30% r mre f the lcal acustic velcity. hen, the fluid density

More information

On Boussinesq's problem

On Boussinesq's problem Internatinal Jurnal f Engineering Science 39 (2001) 317±322 www.elsevier.cm/lcate/ijengsci On Bussinesq's prblem A.P.S. Selvadurai * Department f Civil Engineering and Applied Mechanics, McGill University,

More information

Lyapunov Stability Stability of Equilibrium Points

Lyapunov Stability Stability of Equilibrium Points Lyapunv Stability Stability f Equilibrium Pints 1. Stability f Equilibrium Pints - Definitins In this sectin we cnsider n-th rder nnlinear time varying cntinuus time (C) systems f the frm x = f ( t, x),

More information

Computational modeling techniques

Computational modeling techniques Cmputatinal mdeling techniques Lecture 4: Mdel checing fr ODE mdels In Petre Department f IT, Åb Aademi http://www.users.ab.fi/ipetre/cmpmd/ Cntent Stichimetric matrix Calculating the mass cnservatin relatins

More information

22.54 Neutron Interactions and Applications (Spring 2004) Chapter 11 (3/11/04) Neutron Diffusion

22.54 Neutron Interactions and Applications (Spring 2004) Chapter 11 (3/11/04) Neutron Diffusion .54 Neutrn Interactins and Applicatins (Spring 004) Chapter (3//04) Neutrn Diffusin References -- J. R. Lamarsh, Intrductin t Nuclear Reactr Thery (Addisn-Wesley, Reading, 966) T study neutrn diffusin

More information

Kinematic transformation of mechanical behavior Neville Hogan

Kinematic transformation of mechanical behavior Neville Hogan inematic transfrmatin f mechanical behavir Neville Hgan Generalized crdinates are fundamental If we assume that a linkage may accurately be described as a cllectin f linked rigid bdies, their generalized

More information

Supporting information

Supporting information Electrnic Supplementary Material (ESI) fr Physical Chemistry Chemical Physics This jurnal is The wner Scieties 01 ydrgen perxide electrchemistry n platinum: twards understanding the xygen reductin reactin

More information

Lead/Lag Compensator Frequency Domain Properties and Design Methods

Lead/Lag Compensator Frequency Domain Properties and Design Methods Lectures 6 and 7 Lead/Lag Cmpensatr Frequency Dmain Prperties and Design Methds Definitin Cnsider the cmpensatr (ie cntrller Fr, it is called a lag cmpensatr s K Fr s, it is called a lead cmpensatr Ntatin

More information

bulk velocity through orifice,

bulk velocity through orifice, 150A Review Sessin Other Frictin Lsses Bernulli hf accunts fr all types f drag: is drag due t skin frictin is drag due t fittings (tabulated fractin f the velcity head) is drag due t units (a given r calculated

More information

MXNSSTRY OF SUPPLY AERONAUTICAL RESEARCH COUNCIL REPORTS AND MEMORANDA. of 7 a Pitot ~ntake the Air Stream. X. PLACE and R.

MXNSSTRY OF SUPPLY AERONAUTICAL RESEARCH COUNCIL REPORTS AND MEMORANDA. of 7 a Pitot ~ntake the Air Stream. X. PLACE and R. P~. & M. N. 2621 A.~.G. Technical Reprt t.i ~ ;2 I!Z:, ' MXNSSTRY OF SUPPLY AERONAUTICAL RESEARCH COUNCIL REPORTS AND MEMORANDA :!' ~ i " ',' ~i~ i''~ 'i j i[;_ " ~', c,-, '....~ "-:~ 7 ):, ~Fhe, Efiqciency

More information

Example 1. A robot has a mass of 60 kg. How much does that robot weigh sitting on the earth at sea level? Given: m. Find: Relationships: W

Example 1. A robot has a mass of 60 kg. How much does that robot weigh sitting on the earth at sea level? Given: m. Find: Relationships: W Eample 1 rbt has a mass f 60 kg. Hw much des that rbt weigh sitting n the earth at sea level? Given: m Rbt = 60 kg ind: Rbt Relatinships: Slutin: Rbt =589 N = mg, g = 9.81 m/s Rbt = mrbt g = 60 9. 81 =

More information

Chapter 4. Unsteady State Conduction

Chapter 4. Unsteady State Conduction Chapter 4 Unsteady State Cnductin Chapter 5 Steady State Cnductin Chee 318 1 4-1 Intrductin ransient Cnductin Many heat transfer prblems are time dependent Changes in perating cnditins in a system cause

More information

Yeu-Sheng Paul Shiue, Ph.D 薛宇盛 Professor and Chair Mechanical Engineering Department Christian Brothers University 650 East Parkway South Memphis, TN

Yeu-Sheng Paul Shiue, Ph.D 薛宇盛 Professor and Chair Mechanical Engineering Department Christian Brothers University 650 East Parkway South Memphis, TN Yeu-Sheng Paul Shiue, Ph.D 薛宇盛 Prfessr and Chair Mechanical Engineering Department Christian Brthers University 650 East Parkway Suth Memphis, TN 38104 Office: (901) 321-3424 Rm: N-110 Fax : (901) 321-3402

More information

COASTAL ENGINEERING Chapter 2

COASTAL ENGINEERING Chapter 2 CASTAL ENGINEERING Chapter 2 GENERALIZED WAVE DIFFRACTIN DIAGRAMS J. W. Jhnsn Assciate Prfessr f Mechanical Engineering University f Califrnia Berkeley, Califrnia INTRDUCTIN Wave diffractin is the phenmenn

More information

Computational modeling techniques

Computational modeling techniques Cmputatinal mdeling techniques Lecture 2: Mdeling change. In Petre Department f IT, Åb Akademi http://users.ab.fi/ipetre/cmpmd/ Cntent f the lecture Basic paradigm f mdeling change Examples Linear dynamical

More information

Q1. A) 48 m/s B) 17 m/s C) 22 m/s D) 66 m/s E) 53 m/s. Ans: = 84.0 Q2.

Q1. A) 48 m/s B) 17 m/s C) 22 m/s D) 66 m/s E) 53 m/s. Ans: = 84.0 Q2. Phys10 Final-133 Zer Versin Crdinatr: A.A.Naqvi Wednesday, August 13, 014 Page: 1 Q1. A string, f length 0.75 m and fixed at bth ends, is vibrating in its fundamental mde. The maximum transverse speed

More information

Suggested reading: Lackmann (2011), Sections

Suggested reading: Lackmann (2011), Sections QG Thery and Applicatins: Apprximatins and Equatins Atms 5110 Synptic Dynamic Meterlgy I Instructr: Jim Steenburgh jim.steenburgh@utah.edu 801-581-8727 Suite 480/Office 488 INSCC Suggested reading: Lackmann

More information

Materials Engineering 272-C Fall 2001, Lecture 7 & 8 Fundamentals of Diffusion

Materials Engineering 272-C Fall 2001, Lecture 7 & 8 Fundamentals of Diffusion Materials Engineering 272-C Fall 2001, Lecture 7 & 8 Fundamentals f Diffusin Diffusin: Transprt in a slid, liquid, r gas driven by a cncentratin gradient (r, in the case f mass transprt, a chemical ptential

More information

Introduction: A Generalized approach for computing the trajectories associated with the Newtonian N Body Problem

Introduction: A Generalized approach for computing the trajectories associated with the Newtonian N Body Problem A Generalized apprach fr cmputing the trajectries assciated with the Newtnian N Bdy Prblem AbuBar Mehmd, Syed Umer Abbas Shah and Ghulam Shabbir Faculty f Engineering Sciences, GIK Institute f Engineering

More information

Thermodynamics and Equilibrium

Thermodynamics and Equilibrium Thermdynamics and Equilibrium Thermdynamics Thermdynamics is the study f the relatinship between heat and ther frms f energy in a chemical r physical prcess. We intrduced the thermdynamic prperty f enthalpy,

More information

Fall 2013 Physics 172 Recitation 3 Momentum and Springs

Fall 2013 Physics 172 Recitation 3 Momentum and Springs Fall 03 Physics 7 Recitatin 3 Mmentum and Springs Purpse: The purpse f this recitatin is t give yu experience wrking with mmentum and the mmentum update frmula. Readings: Chapter.3-.5 Learning Objectives:.3.

More information

37 Maxwell s Equations

37 Maxwell s Equations 37 Maxwell s quatins In this chapter, the plan is t summarize much f what we knw abut electricity and magnetism in a manner similar t the way in which James Clerk Maxwell summarized what was knwn abut

More information

Physics 212. Lecture 12. Today's Concept: Magnetic Force on moving charges. Physics 212 Lecture 12, Slide 1

Physics 212. Lecture 12. Today's Concept: Magnetic Force on moving charges. Physics 212 Lecture 12, Slide 1 Physics 1 Lecture 1 Tday's Cncept: Magnetic Frce n mving charges F qv Physics 1 Lecture 1, Slide 1 Music Wh is the Artist? A) The Meters ) The Neville rthers C) Trmbne Shrty D) Michael Franti E) Radiatrs

More information

Phys102 Final-061 Zero Version Coordinator: Nasser Wednesday, January 24, 2007 Page: 1

Phys102 Final-061 Zero Version Coordinator: Nasser Wednesday, January 24, 2007 Page: 1 Crdinatr: Nasser Wednesday, January 4, 007 Page: 1 Q1. Tw transmitters, S 1 and S shwn in the figure, emit identical sund waves f wavelength λ. The transmitters are separated by a distance λ /. Cnsider

More information

Dispersion Ref Feynman Vol-I, Ch-31

Dispersion Ref Feynman Vol-I, Ch-31 Dispersin Ref Feynman Vl-I, Ch-31 n () = 1 + q N q /m 2 2 2 0 i ( b/m) We have learned that the index f refractin is nt just a simple number, but a quantity that varies with the frequency f the light.

More information

Physics 2B Chapter 23 Notes - Faraday s Law & Inductors Spring 2018

Physics 2B Chapter 23 Notes - Faraday s Law & Inductors Spring 2018 Michael Faraday lived in the Lndn area frm 1791 t 1867. He was 29 years ld when Hand Oersted, in 1820, accidentally discvered that electric current creates magnetic field. Thrugh empirical bservatin and

More information

making triangle (ie same reference angle) ). This is a standard form that will allow us all to have the X= y=

making triangle (ie same reference angle) ). This is a standard form that will allow us all to have the X= y= Intrductin t Vectrs I 21 Intrductin t Vectrs I 22 I. Determine the hrizntal and vertical cmpnents f the resultant vectr by cunting n the grid. X= y= J. Draw a mangle with hrizntal and vertical cmpnents

More information

EHed of Curvature on the Temperature Profiles

EHed of Curvature on the Temperature Profiles PROC. OF THE OKLA. ACAD. OF SCI. FOR 1967 EHed f Curvature n the Temperature Prfiles in Cnduding Spines J. E. FRANCIS add R. V. KASER, University f Oklahma, Nrman and GORDON SCOFIELD, University f Missuri,

More information

I. Analytical Potential and Field of a Uniform Rod. V E d. The definition of electric potential difference is

I. Analytical Potential and Field of a Uniform Rod. V E d. The definition of electric potential difference is Length L>>a,b,c Phys 232 Lab 4 Ch 17 Electric Ptential Difference Materials: whitebards & pens, cmputers with VPythn, pwer supply & cables, multimeter, crkbard, thumbtacks, individual prbes and jined prbes,

More information

Phys101 Final Code: 1 Term: 132 Wednesday, May 21, 2014 Page: 1

Phys101 Final Code: 1 Term: 132 Wednesday, May 21, 2014 Page: 1 Phys101 Final Cde: 1 Term: 1 Wednesday, May 1, 014 Page: 1 Q1. A car accelerates at.0 m/s alng a straight rad. It passes tw marks that are 0 m apart at times t = 4.0 s and t = 5.0 s. Find the car s velcity

More information

Chapter 3 Kinematics in Two Dimensions; Vectors

Chapter 3 Kinematics in Two Dimensions; Vectors Chapter 3 Kinematics in Tw Dimensins; Vectrs Vectrs and Scalars Additin f Vectrs Graphical Methds (One and Tw- Dimensin) Multiplicatin f a Vectr b a Scalar Subtractin f Vectrs Graphical Methds Adding Vectrs

More information

Chapter 16. Capacitance. Capacitance, cont. Parallel-Plate Capacitor, Example 1/20/2011. Electric Energy and Capacitance

Chapter 16. Capacitance. Capacitance, cont. Parallel-Plate Capacitor, Example 1/20/2011. Electric Energy and Capacitance summary C = ε A / d = πε L / ln( b / a ) ab C = 4πε 4πε a b a b >> a Chapter 16 Electric Energy and Capacitance Capacitance Q=CV Parallel plates, caxial cables, Earth Series and parallel 1 1 1 = + +..

More information

FIELD QUALITY IN ACCELERATOR MAGNETS

FIELD QUALITY IN ACCELERATOR MAGNETS FIELD QUALITY IN ACCELERATOR MAGNETS S. Russenschuck CERN, 1211 Geneva 23, Switzerland Abstract The field quality in the supercnducting magnets is expressed in terms f the cefficients f the Furier series

More information

This section is primarily focused on tools to aid us in finding roots/zeros/ -intercepts of polynomials. Essentially, our focus turns to solving.

This section is primarily focused on tools to aid us in finding roots/zeros/ -intercepts of polynomials. Essentially, our focus turns to solving. Sectin 3.2: Many f yu WILL need t watch the crrespnding vides fr this sectin n MyOpenMath! This sectin is primarily fcused n tls t aid us in finding rts/zers/ -intercepts f plynmials. Essentially, ur fcus

More information

Kinetic Model Completeness

Kinetic Model Completeness 5.68J/10.652J Spring 2003 Lecture Ntes Tuesday April 15, 2003 Kinetic Mdel Cmpleteness We say a chemical kinetic mdel is cmplete fr a particular reactin cnditin when it cntains all the species and reactins

More information

CHAPTER 8b Static Equilibrium Units

CHAPTER 8b Static Equilibrium Units CHAPTER 8b Static Equilibrium Units The Cnditins fr Equilibrium Slving Statics Prblems Stability and Balance Elasticity; Stress and Strain The Cnditins fr Equilibrium An bject with frces acting n it, but

More information

Solution to HW14 Fall-2002

Solution to HW14 Fall-2002 Slutin t HW14 Fall-2002 CJ5 10.CQ.003. REASONING AND SOLUTION Figures 10.11 and 10.14 shw the velcity and the acceleratin, respectively, the shadw a ball that underges unirm circular mtin. The shadw underges

More information

Chapter VII Electrodynamics

Chapter VII Electrodynamics Chapter VII Electrdynamics Recmmended prblems: 7.1, 7., 7.4, 7.5, 7.7, 7.8, 7.10, 7.11, 7.1, 7.13, 7.15, 7.17, 7.18, 7.0, 7.1, 7., 7.5, 7.6, 7.7, 7.9, 7.31, 7.38, 7.40, 7.45, 7.50.. Ohm s Law T make a

More information

Plan o o. I(t) Divide problem into sub-problems Modify schematic and coordinate system (if needed) Write general equations

Plan o o. I(t) Divide problem into sub-problems Modify schematic and coordinate system (if needed) Write general equations STAPLE Physics 201 Name Final Exam May 14, 2013 This is a clsed bk examinatin but during the exam yu may refer t a 5 x7 nte card with wrds f wisdm yu have written n it. There is extra scratch paper available.

More information

Chapter 9 Compressible Flow 667

Chapter 9 Compressible Flow 667 Chapter 9 Cmpreible Flw 667 9.57 Air flw frm a tank thrugh a nzzle int the tandard atmphere, a in Fig. P9.57. A nrmal hck tand in the exit f the nzzle, a hwn. Etimate (a) the tank preure; and (b) the ma

More information

Relationship Between Amplifier Settling Time and Pole-Zero Placements for Second-Order Systems *

Relationship Between Amplifier Settling Time and Pole-Zero Placements for Second-Order Systems * Relatinship Between Amplifier Settling Time and Ple-Zer Placements fr Secnd-Order Systems * Mark E. Schlarmann and Randall L. Geiger Iwa State University Electrical and Cmputer Engineering Department Ames,

More information

Lecture 13: Electrochemical Equilibria

Lecture 13: Electrochemical Equilibria 3.012 Fundamentals f Materials Science Fall 2005 Lecture 13: 10.21.05 Electrchemical Equilibria Tday: LAST TIME...2 An example calculatin...3 THE ELECTROCHEMICAL POTENTIAL...4 Electrstatic energy cntributins

More information

Chapter 23 Electromagnetic Waves Lecture 14

Chapter 23 Electromagnetic Waves Lecture 14 Chapter 23 Electrmagnetic Waves Lecture 14 23.1 The Discvery f Electrmagnetic Waves 23.2 Prperties f Electrmagnetic Waves 23.3 Electrmagnetic Waves Carry Energy and Mmentum 23.4 Types f Electrmagnetic

More information

ON THE EFFECTIVENESS OF POROSITY ON UNSTEADY COUETTE FLOW AND HEAT TRANSFER BETWEEN PARALLEL POROUS PLATES WITH EXPONENTIAL DECAYING PRESSURE GRADIENT

ON THE EFFECTIVENESS OF POROSITY ON UNSTEADY COUETTE FLOW AND HEAT TRANSFER BETWEEN PARALLEL POROUS PLATES WITH EXPONENTIAL DECAYING PRESSURE GRADIENT 17 Kragujevac J. Sci. 8 (006) 17-4. ON THE EFFECTIVENESS OF POROSITY ON UNSTEADY COUETTE FLOW AND HEAT TRANSFER BETWEEN PARALLEL POROUS PLATES WITH EXPONENTIAL DECAYING PRESSURE GRADIENT Hazem Ali Attia

More information

The Sputtering Problem James A Glackin, James V. Matheson

The Sputtering Problem James A Glackin, James V. Matheson The Sputtering Prblem James A Glackin, James V. Mathesn I prpse t cnsider first the varius elements f the subject, next its varius parts r sectins, and finally the whle in its internal structure. In ther

More information

Assume that the water in the nozzle is accelerated at a rate such that the frictional effect can be neglected.

Assume that the water in the nozzle is accelerated at a rate such that the frictional effect can be neglected. 1 HW #3: Cnservatin f Linear Mmentum, Cnservatin f Energy, Cnservatin f Angular Mmentum and Turbmachines, Bernulli s Equatin, Dimensinal Analysis, and Pipe Flws Prblem 1. Cnservatins f Mass and Linear

More information

Department of Economics, University of California, Davis Ecn 200C Micro Theory Professor Giacomo Bonanno. Insurance Markets

Department of Economics, University of California, Davis Ecn 200C Micro Theory Professor Giacomo Bonanno. Insurance Markets Department f Ecnmics, University f alifrnia, Davis Ecn 200 Micr Thery Prfessr Giacm Bnann Insurance Markets nsider an individual wh has an initial wealth f. ith sme prbability p he faces a lss f x (0

More information

OF SIMPLY SUPPORTED PLYWOOD PLATES UNDER COMBINED EDGEWISE BENDING AND COMPRESSION

OF SIMPLY SUPPORTED PLYWOOD PLATES UNDER COMBINED EDGEWISE BENDING AND COMPRESSION U. S. FOREST SERVICE RESEARCH PAPER FPL 50 DECEMBER U. S. DEPARTMENT OF AGRICULTURE FOREST SERVICE FOREST PRODUCTS LABORATORY OF SIMPLY SUPPORTED PLYWOOD PLATES UNDER COMBINED EDGEWISE BENDING AND COMPRESSION

More information

Numerical Simulation of the Thermal Resposne Test Within the Comsol Multiphysics Environment

Numerical Simulation of the Thermal Resposne Test Within the Comsol Multiphysics Environment Presented at the COMSOL Cnference 2008 Hannver University f Parma Department f Industrial Engineering Numerical Simulatin f the Thermal Respsne Test Within the Cmsl Multiphysics Envirnment Authr : C. Crradi,

More information

Aerodynamic Separability in Tip Speed Ratio and Separability in Wind Speed- a Comparison

Aerodynamic Separability in Tip Speed Ratio and Separability in Wind Speed- a Comparison Jurnal f Physics: Cnference Series OPEN ACCESS Aerdynamic Separability in Tip Speed Rati and Separability in Wind Speed- a Cmparisn T cite this article: M L Gala Sants et al 14 J. Phys.: Cnf. Ser. 555

More information

ENSC Discrete Time Systems. Project Outline. Semester

ENSC Discrete Time Systems. Project Outline. Semester ENSC 49 - iscrete Time Systems Prject Outline Semester 006-1. Objectives The gal f the prject is t design a channel fading simulatr. Upn successful cmpletin f the prject, yu will reinfrce yur understanding

More information

Phy 213: General Physics III 6/14/2007 Chapter 28 Worksheet 1

Phy 213: General Physics III 6/14/2007 Chapter 28 Worksheet 1 Ph 13: General Phsics III 6/14/007 Chapter 8 Wrksheet 1 Magnetic Fields & Frce 1. A pint charge, q= 510 C and m=110-3 m kg, travels with a velcit f: v = 30 ˆ s i then enters a magnetic field: = 110 T ˆj.

More information

L a) Calculate the maximum allowable midspan deflection (w o ) critical under which the beam will slide off its support.

L a) Calculate the maximum allowable midspan deflection (w o ) critical under which the beam will slide off its support. ecture 6 Mderately arge Deflectin Thery f Beams Prblem 6-1: Part A: The department f Highways and Public Wrks f the state f Califrnia is in the prcess f imprving the design f bridge verpasses t meet earthquake

More information

Sections 15.1 to 15.12, 16.1 and 16.2 of the textbook (Robbins-Miller) cover the materials required for this topic.

Sections 15.1 to 15.12, 16.1 and 16.2 of the textbook (Robbins-Miller) cover the materials required for this topic. Tpic : AC Fundamentals, Sinusidal Wavefrm, and Phasrs Sectins 5. t 5., 6. and 6. f the textbk (Rbbins-Miller) cver the materials required fr this tpic.. Wavefrms in electrical systems are current r vltage

More information

Computational modeling techniques

Computational modeling techniques Cmputatinal mdeling techniques Lecture 11: Mdeling with systems f ODEs In Petre Department f IT, Ab Akademi http://www.users.ab.fi/ipetre/cmpmd/ Mdeling with differential equatins Mdeling strategy Fcus

More information

CHAPTER 6 -- ENERGY. Approach #2: Using the component of mg along the line of d:

CHAPTER 6 -- ENERGY. Approach #2: Using the component of mg along the line of d: Slutins--Ch. 6 (Energy) CHAPTER 6 -- ENERGY 6.) The f.b.d. shwn t the right has been prvided t identify all the frces acting n the bdy as it mves up the incline. a.) T determine the wrk dne by gravity

More information

Problem 1 Known: Dimensions and materials of the composition wall, 10 studs each with 2.5m high

Problem 1 Known: Dimensions and materials of the composition wall, 10 studs each with 2.5m high Prblem Knwn: Dimensins and materials f the cmpsitin wall, 0 studs each with.5m high Unknwn:. Thermal resistance assciate with wall when surfaces nrmal t the directin f heat flw are isthermal. Thermal resistance

More information

Information for Physics 1201 Midterm I Wednesday, February 20

Information for Physics 1201 Midterm I Wednesday, February 20 My lecture slides are psted at http://www.physics.hi-state.edu/~humanic/ Infrmatin fr Physics 1201 Midterm I Wednesday, February 20 1) Frmat: 10 multiple chice questins (each wrth 5 pints) and tw shw-wrk

More information

Modeling the Nonlinear Rheological Behavior of Materials with a Hyper-Exponential Type Function

Modeling the Nonlinear Rheological Behavior of Materials with a Hyper-Exponential Type Function www.ccsenet.rg/mer Mechanical Engineering Research Vl. 1, N. 1; December 011 Mdeling the Nnlinear Rhelgical Behavir f Materials with a Hyper-Expnential Type Functin Marc Delphin Mnsia Département de Physique,

More information

Lab 1 The Scientific Method

Lab 1 The Scientific Method INTRODUCTION The fllwing labratry exercise is designed t give yu, the student, an pprtunity t explre unknwn systems, r universes, and hypthesize pssible rules which may gvern the behavir within them. Scientific

More information

THE FLUXOID QUANTUM AND ELECTROGRAVITATIONAL DYNAMICS. Chapter 8. This work extends chapter 6 titled, "Field Mass Generation and Control", while

THE FLUXOID QUANTUM AND ELECTROGRAVITATIONAL DYNAMICS. Chapter 8. This work extends chapter 6 titled, Field Mass Generation and Control, while 133 THE FLUXOID QUANTUM AND ELECTROGRAVITATIONAL DYNAMICS Chapter 8 This wrk extends chapter 6 titled, "Field Mass Generatin and Cntrl", while als develping a new cnceptual apprach t mass-field vehicle

More information

1. Transformer A transformer is used to obtain the approximate output voltage of the power supply. The output of the transformer is still AC.

1. Transformer A transformer is used to obtain the approximate output voltage of the power supply. The output of the transformer is still AC. PHYSIS 536 Experiment 4: D Pwer Supply I. Intrductin The prcess f changing A t D is investigated in this experiment. An integrated circuit regulatr makes it easy t cnstruct a high-perfrmance vltage surce

More information

AP Statistics Notes Unit Two: The Normal Distributions

AP Statistics Notes Unit Two: The Normal Distributions AP Statistics Ntes Unit Tw: The Nrmal Distributins Syllabus Objectives: 1.5 The student will summarize distributins f data measuring the psitin using quartiles, percentiles, and standardized scres (z-scres).

More information