Equilibrium Reconstruction with Grad-Shafranov Equation
|
|
- Spencer Green
- 5 years ago
- Views:
Transcription
1 Equilibium Reconstuction with Gad-Shafanov Equation Qian Peng (qp215) Apil 14, Gad-Shafanov equation [1]Fo axial symmetic case, the tooidal B eld is B to = F φ. Whee 2π F (, θ) is the total poloidal cuent though the cicle with adius at(, θ). The total B eld is given by B = F φ φ ψ p notice that hee 2π ψ p is the polodial ux, but dened slightly dieent fom that in Plasma II. ψ p is at magnetic axis and is the integation of counte-clockwise B eld fom magnetic axis to the suface. Then we have thus µ j = B = d F d ψ ψ p φ (R p R ( 1 R = d F d ψ p ψ p φ ψp φ j toodial = ψ p µ R ) Z ) ψ p φ It should be pointed out that cuently is dened in cylindical coodinates while othe functions ae dened in ux coodinates. Fo MHD model: µ j B = µ p d F ( F d ψ ψp )( φ) 2 p = d(µ p) d ψ p ψp = 2 d(µ p) d ψ p = 2 P T F d F d ψ p 1
2 Gad-Shafonov equation is ψp = 2 P ( ψ p ) T ( ψ p ) (1) and we have µ j toodial = P + T 1.1 Geen's function Seach fo Geens function G(, z,, z ) such that[2],z G(, z,, z ) = δ(, z z ) (2) G,z In fact the poloidal component of B eld: φ ψ p = ( ψ p φ) and ψ p is thus the φ component of magnetic potential A φ coesponding to the B eld. The physical meaning of the RHS of (2) is a ing cuent (1/µ ) locating at (, z ) and its magnectic potential is A( ) = µ 4π = ( φ) 4π j( ) d3 π π cos(θ) cos(θ) dθ thus π G(, z,, z ) = 4π cos(θ) dθ (3) π cos(θ) 1.2 Solution of GSh equation We only want to solve GSh equation within a domain: ψ p = f(, z) (, z) Ω (4) we can extend Ω to the whole space to include the extenal coils and have ψ p (, z) = f(, z )G(, z,, z )d dz (5) +Σ i I i G(, z, i, z i) whee I i denote extenal coils. This is not vey numeically ecient because if Ω is a N by N space and f(, z) is non-zeo almost evey whee in the 2
3 domain then this would equie O(N 4 ) opeations to obtain the value of ψ p at evey point. Anothe way to solve this is to st constuct the bounday value ψ p (, z) (,z) Ω using (5) (Which will take O(N 3 ) opeations) and then discetize (4) and solve it numeically (Which is expected to take O(N 3 ) opeations). [4] To constuct ψ p on the bounday, we still need to calculate Geen's function O(N 3 ) times befoehand and stoe them in the memoy. We can impove it futhe moe.[3] u = ( (1 = ( 1 u) ) z )u using Geen's identity G(, z,, z ) ψ p (, z) ψ p (, z) G(, z,, z )ddz = G(, z,, z )f(, z) ψ p (, z)δ(, z z )ddz = G(, z,, z )f(, z)ddz ψ p (, z ) since = = we have Ω Ω Ω ψ p (, z ) = G(, z,, z ) ψ p (, z) ψ p (, z) G(, z,, z )ddz G(, z,, z ) ( 1 ψ p (, z)) ψ p (, z) ( 1 G(, z,, z ))ddz G(, z,, z ) + Ω Ω n ψ p (, z) ψ p (, z) n G(, z,, z ) ds,z G(, z,, z )f(, z)ddz ψ p (, z) n G(, z,, z ) G(, z,, z ) n ψ p (, z) ds,z (6) Notice that the st tem of RHS of (6) is what we need in (5), we can obtain the value of this tem by the following pocedue: Solve the equation numeically by discetizing the opeato ψ p = f(, z) (, z) Ω ψ p Ω = 3
4 whee the domain Ω is lage than Ω. Then calculate ψ n p on the bounday and since ψ p is zeo on the bounday, we have G(, z,, z )f(, z)ddz = G(, z,, z )f(, z)ddz Ω Ω = ψ p (, z G(, z,, z ) ) + Ω n ψ p (, z) ds,z and we only need to know the value of G(, z,, z )with (, z) Ω, (, z ) Ω. Substitute this into (5) we have ψ p (, z) = ψ G(, z,, z) p (, z) + Ω n ψ p (, z ) ds,z (7) +ΣI i G(, z, i i, z i) Use (7) to calculate ψ p (, z) on the bounday and then solve it numeically. 1.3 Non-lineaity of GSh equation. In the pevious section we ae assuming the RHS of GSh equation is f(, z) while in eality the RHS of (1) is f( ψ p ) which depends on the function of p( ψ p ) and F ( ψ p ). Even assuming that p and F is given. f is a function of (, z) in a sense that we aleady know the function ψ p (, z) and f(, z) = f( ψ p (, z)). But ψ p (, z) is the function that we ae solving fo! One way to teat this poblem is to solve this equation iteatively. Stat with an initial guess of ψ p, (, z) and use this to calculate f n+1 (, z) = f( ψ p,n (, z)) and then solve the equation to obtain ψ p,n+1, iteate until ψp conveges. In fact, the plasma bounday can also change duing each iteation. Plasma bounday is detemined by a set of limites located at ( s, z s ) in this case and the limite with lowest ψ p,n,min = ψ p,n ( s, z s ) detemines the bounday of plasma. The function of p and F ae actually functions of nomalized poloidal ux p( ψ p (, z)/ ψ p,min (, z)), F ( ψp (, z)/ ψ p,min (, z)) because they ae dened on the plasma egion. 1.4 GSh equation on ux coodinates [1]Solving GSh equation on (, z) coodinates gives simplicity in discetizing the opeato, but it it may not be the most ecient way. Use geneal coodinates (a, θ) to descibe the (, z) plane, whee a can be teated as a suface label and θ the poloidal angle. Let D(, z) D = D(a, θ) 4
5 then whee f = ( 1 f) = 1 D a (DG a,a a f) + 1 D a (DG a,θ θ f) + 1 D θ (DG θ,a a f) + 1 D θ (DG θ,θ θ f) G a,a = x a x G a,θ = G θ,a = x a a x θ The GSh equation is thus D ψ p = DP D T G θ,θ = x θ x θ = V P LT whee V D, L = D If (a, θ) is aleady the ux coodinate (a, θ) that ψ p = ψ p (a ), p = p(a ), F = F (a ) then we shall have a (DG a,a aveaging ove θ we have ψ p) + θ (DG θ,a ψ p) = V P LT (8) (( DG a,a ) ψ p ) = T L P V, whee (f) = 1 2π Dene the total coss-section cuent J(a ) = 1 2π = 1 2π = Compae it with (9) we have: and since we have Which is an 1D equation. a 2π 2π dθ dθ a a V P + L T da µ j tooidal D(, z) D(a, θ) da DP + DT da (( DG a,a ) ψ p ) a = J(a ) B θ (a = ) = ψ p a= = fdθ (9) ( DG a,a ) ψ p = J (1) 5
6 1.5 Petubed GSh equation [1]When solving the GSh equation, we don't eally have the ux coodinate (a, θ) until we get the tue solution. Howeve, we can use the coodinate (a, θ) geneated fom each appoximation of ψp and assume that (a, θ) is close to (a, θ) in the sense that a = a ξ(a, θ) then (, z) which is initially a function of (a, θ) can be expessed as = (a, θ) = (a + ξ, θ) (a, θ) + aξ z = z(a, θ) = z(a + ξ, θ) z(a, θ) + z aξ the suface of (a, θ) has to satisfy that ψ p (a, θ) ψ p + ψ(a, θ) = const whee ψ p is the main ode appoximation and can be obtained by solving the equation (1), in the coodinate (a, θ) but taken as (a, θ). Thus: ψ pξ + ψ = const Notice that ξ have the pat ξ and oscillatoy pat ξ thus ξ = ψ ψ p ξ = ψ p (a, θ) ψ p ψ ψ p chose ψ p (a, θ) = ψ p then ξ = ψ ψ p Now, GSh equation can be lineaized: ψp = 2 P (a ) T (a ) with P (a ) = P (a) P ξ + δp T (a ) = T (a) T ξ + δt (δp and δt denote the vaiations of the given function pole which will be used late.) such that L ψ p = V P LT + V P ξ + LT ξ V δp LδT (11) Fo the pevious algoithms, at each step, when ψ p n n+1 ψ p, only the change to the LHS of the above equation is consideed. If one can also calculated the elated ξ due to the change in ψ p, and use both of them in the above equation, a bette conveging ate should be achieved. 6
7 2 Equilibium econstuction Solving GSh equation is one majo pat of econstucting the equilibium pole. Yet, one also needs to have the pole of P ( ψ p ) and T ( ψ p ) to solve the equation. These poles ae the input paametes fo GSh equation solve and ae often detemined by the measuement fom pobes. We stat with an initial guess of the poles with paametes P ( ψ p ) = P (α 1, α 2,..., α n, ψ p ) T ( ψ p ) = T (β 1, β 2,..., β n, ψ p ) and solve fo ψ p (, z). Then we have P (α 1, α 2,..., α n,, z) and T (β 1, β 2,..., β n,, z), fom which we can get the expect eading fom the ith pobe E i (α 1, α 2,..., α n, β 1, β 2,..., β n ) Since we also have the actually eadings, the paametes can be adjusted to bette t the eading. Update the paametes and solve the GSh equation again until it conveges. [4] One can see that the convegence ate does not tend to be vey fast, because the tting pocedue does not take into account that ψ p (, z) would be changed when the paametes change. A bette way to do it might be using δp and δt in the equation (11) and obtain an estimated esponse of ψp. Taking that into account when doing the tting pocesue may esult in a bette convegence ate. 2.1 Eects of measuing techniques on econstuction. It is a commen sense that the moe accuate you measuements ae, the bette econstuction you will get. But, it would be bette if one can tell quantitatively how does the accuacy of the pobe eadings aects the cedibility of one's econstuction esult. In L.E. Zakhaov's pape[5], he gives the elations quantitatively. The idea is to view eos as petubations and see how lage the magnitude of the petubation could be as long as they ae undetectable by the device. When we change the P and T pole by δp and δt, it will esult in a change of δ ψ p and eventually a change of signals we ead, δs. Fo smal petubation, we can lineaize the elationship and have δs = AδX whee δs is a N dimensional vecto of all the signals fom pobes and δx ae a M dimensional discete vesion of the change of δp and δt. Since all pobes have cetain degees of eo allowance, we use δɛ to denote that. Nomalize each ow of A(i, :) with δɛ i and we get Ã. Thus, δ S = ÃδX 7
8 and each enty of δ S has to have a magnitude less than 1 so that the petubation is within the eo ange of pobes. One can single value decompse à = UW V T and fo δx k = γv k we have δ S k = γw k U k. Set cetain constains fo the nom of δs and we'll get the eo ange δx fo the econstuction. Refeences [1] `Theoy of petubed equilibia fo solving the GadShafanov equation.' Zakhaov, L E Pletze, A (1999). [2] `Numeical Solution of Gad-Shafanov Equation fo the Distibution of Magnetic Flux in Nuclea Fusion Devices.' Aydin, Selcuk Han, Teze-Sezgin Munevve. (28). [3] `Computation of ideal MHD equilibia.' Lackne K, (1976). [4] `Numeical detemination of axisymmetic tooidal magnetohydodynamic equilibia.' Johnson, J L Dalhed, H E Geene, J M Gimm, R C Hsieh, Y Y Jadin, S C Manickam, J Okabayashi, M Stoe, R G Todd, A M M Voss, D E Weime, K E. (1979). [5] `Reconstuction of the p and q poles in ITER fom extenal and intenal measuements.' Zakhaov. L, Foley. E, Levinton. F, Yuh. H. (28). 8
12th WSEAS Int. Conf. on APPLIED MATHEMATICS, Cairo, Egypt, December 29-31,
th WSEAS Int. Conf. on APPLIED MATHEMATICS, Caio, Egypt, Decembe 9-3, 7 5 Magnetostatic Field calculations associated with thick Solenoids in the Pesence of Ion using a Powe Seies expansion and the Complete
More informationAbsorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere
Applied Mathematics, 06, 7, 709-70 Published Online Apil 06 in SciRes. http://www.scip.og/jounal/am http://dx.doi.og/0.46/am.06.77065 Absoption Rate into a Small Sphee fo a Diffusing Paticle Confined in
More informationConservative Averaging Method and its Application for One Heat Conduction Problem
Poceedings of the 4th WSEAS Int. Conf. on HEAT TRANSFER THERMAL ENGINEERING and ENVIRONMENT Elounda Geece August - 6 (pp6-) Consevative Aveaging Method and its Application fo One Heat Conduction Poblem
More informationJ. Electrical Systems 1-3 (2005): Regular paper
K. Saii D. Rahem S. Saii A Miaoui Regula pape Coupled Analytical-Finite Element Methods fo Linea Electomagnetic Actuato Analysis JES Jounal of Electical Systems In this pape, a linea electomagnetic actuato
More informationReview: Electrostatics and Magnetostatics
Review: Electostatics and Magnetostatics In the static egime, electomagnetic quantities do not vay as a function of time. We have two main cases: ELECTROSTATICS The electic chages do not change postion
More informationBoundary Integral Equation Approach Based on. Polynomial Expansion of Current Distribution to. Reconstruct Current Density Profile in Tokamak Plasma
Title Bounday integal equation appoach based on a polyn cuent density pofile in tokamak plasmas Autho(s)Itagaki, Masafumi; Yamaguchi, Satoki; Fukunaga, Taka CitationNuclea Fusion, 45(3): 153-162 Issue
More informationOn the integration of the equations of hydrodynamics
Uebe die Integation de hydodynamischen Gleichungen J f eine u angew Math 56 (859) -0 On the integation of the equations of hydodynamics (By A Clebsch at Calsuhe) Tanslated by D H Delphenich In a pevious
More informationMAGNETIC FIELD AROUND TWO SEPARATED MAGNETIZING COILS
The 8 th Intenational Confeence of the Slovenian Society fo Non-Destuctive Testing»pplication of Contempoay Non-Destuctive Testing in Engineeing«Septembe 1-3, 5, Potoož, Slovenia, pp. 17-1 MGNETIC FIELD
More information11) A thin, uniform rod of mass M is supported by two vertical strings, as shown below.
Fall 2007 Qualifie Pat II 12 minute questions 11) A thin, unifom od of mass M is suppoted by two vetical stings, as shown below. Find the tension in the emaining sting immediately afte one of the stings
More information1D2G - Numerical solution of the neutron diffusion equation
DG - Numeical solution of the neuton diffusion equation Y. Danon Daft: /6/09 Oveview A simple numeical solution of the neuton diffusion equation in one dimension and two enegy goups was implemented. Both
More informationPOISSON S EQUATION 2 V 0
POISSON S EQUATION We have seen how to solve the equation but geneally we have V V4k We now look at a vey geneal way of attacking this poblem though Geen s Functions. It tuns out that this poblem has applications
More informationEEO 401 Digital Signal Processing Prof. Mark Fowler
EEO 41 Digital Signal Pocessing Pof. Mak Fowle Note Set #31 Linea Phase FIR Design Optimum Equiipple (Paks-McClellan) Reading: Sect. 1.2.4 1.2.6 of Poakis & Manolakis 1/2 Motivation The window method and
More informationConjugate Gradient Methods. Michael Bader. Summer term 2012
Gadient Methods Outlines Pat I: Quadatic Foms and Steepest Descent Pat II: Gadients Pat III: Summe tem 2012 Pat I: Quadatic Foms and Steepest Descent Outlines Pat I: Quadatic Foms and Steepest Descent
More informationFE FORMULATIONS FOR PLASTICITY
G These slides ae designed based on the book: Finite Elements in Plasticity Theoy and Pactice, D.R.J. Owen and E. Hinton, 970, Pineidge Pess Ltd., Swansea, UK. Couse Content: A INTRODUCTION AND OVERVIEW
More informationI. CONSTRUCTION OF THE GREEN S FUNCTION
I. CONSTRUCTION OF THE GREEN S FUNCTION The Helmohltz equation in 4 dimensions is 4 + k G 4 x, x = δ 4 x x. In this equation, G is the Geen s function and 4 efes to the dimensionality. In the vey end,
More information( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx.
9. LAGRANGIAN OF THE ELECTROMAGNETIC FIELD In the pevious section the Lagangian and Hamiltonian of an ensemble of point paticles was developed. This appoach is based on a qt. This discete fomulation can
More informationLecture 8 - Gauss s Law
Lectue 8 - Gauss s Law A Puzzle... Example Calculate the potential enegy, pe ion, fo an infinite 1D ionic cystal with sepaation a; that is, a ow of equally spaced chages of magnitude e and altenating sign.
More informationMathematical Model of Magnetometric Resistivity. Sounding for a Conductive Host. with a Bulge Overburden
Applied Mathematical Sciences, Vol. 7, 13, no. 7, 335-348 Mathematical Model of Magnetometic Resistivity Sounding fo a Conductive Host with a Bulge Ovebuden Teeasak Chaladgan Depatment of Mathematics Faculty
More informationA 1. EN2210: Continuum Mechanics. Homework 7: Fluid Mechanics Solutions
EN10: Continuum Mechanics Homewok 7: Fluid Mechanics Solutions School of Engineeing Bown Univesity 1. An ideal fluid with mass density ρ flows with velocity v 0 though a cylindical tube with cosssectional
More informationAs is natural, our Aerospace Structures will be described in a Euclidean three-dimensional space R 3.
Appendix A Vecto Algeba As is natual, ou Aeospace Stuctues will be descibed in a Euclidean thee-dimensional space R 3. A.1 Vectos A vecto is used to epesent quantities that have both magnitude and diection.
More informationApplication of homotopy perturbation method to the Navier-Stokes equations in cylindrical coordinates
Computational Ecology and Softwae 5 5(): 9-5 Aticle Application of homotopy petubation method to the Navie-Stokes equations in cylindical coodinates H. A. Wahab Anwa Jamal Saia Bhatti Muhammad Naeem Muhammad
More informationCoarse Mesh Radiation Transport Code COMET Radiation Therapy Application*
Coase Mesh Radiation Tanspot Code COMT Radiation Theapy Application* Fazad Rahnema Nuclea & Radiological ngineeing and Medical Physics Pogams Geogia Institute of Technology Computational Medical Physics
More informationSOLVING THE VISCOUS COMPOSITE CYLINDER PROBLEM BY SOKOLOV S METHOD
1 SOLVING THE VISCOUS COMPOSITE CYLINDER PROBLEM BY SOKOLOV S METHOD NGO By CO. H. TRAN, PHONG. T. Faculty of Mathematics & Infomatics, Univesity of Natual Sciences VNU-HCM Abstact : The pape pesents some
More informationAn Exact Solution of Navier Stokes Equation
An Exact Solution of Navie Stokes Equation A. Salih Depatment of Aeospace Engineeing Indian Institute of Space Science and Technology, Thiuvananthapuam, Keala, India. July 20 The pincipal difficulty in
More informationTUTORIAL 9. Static magnetic field
TUTOIAL 9 Static magnetic field Vecto magnetic potential Null Identity % & %$ A # Fist postulation # " B such that: Vecto magnetic potential Vecto Poisson s equation The solution is: " Substitute it into
More informationChapter 3 Optical Systems with Annular Pupils
Chapte 3 Optical Systems with Annula Pupils 3 INTRODUCTION In this chapte, we discuss the imaging popeties of a system with an annula pupil in a manne simila to those fo a system with a cicula pupil The
More informationQuantum Mechanics II
Quantum Mechanics II Pof. Bois Altshule Apil 25, 2 Lectue 25 We have been dicussing the analytic popeties of the S-matix element. Remembe the adial wave function was u kl () = R kl () e ik iπl/2 S l (k)e
More informationPROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.
POBLM S # SOLUIONS by obet A. DiStasio J. Q. he Bon-Oppenheime appoximation is the standad way of appoximating the gound state of a molecula system. Wite down the conditions that detemine the tonic and
More information15 Solving the Laplace equation by Fourier method
5 Solving the Laplace equation by Fouie method I aleady intoduced two o thee dimensional heat equation, when I deived it, ecall that it taes the fom u t = α 2 u + F, (5.) whee u: [0, ) D R, D R is the
More informationMulti-Objective Optimization Algorithms for Finite Element Model Updating
Multi-Objective Optimization Algoithms fo Finite Element Model Updating E. Ntotsios, C. Papadimitiou Univesity of Thessaly Geece Outline STRUCTURAL IDENTIFICATION USING MEASURED MODAL DATA Weighted Modal
More informationMagnetohydrodynamics (MHD) I
Magnetohydodynamics (MHD) I Yong-Su Na National Fusion Reseach Cente POSTECH, Koea, 8-10 May, 006 Contents 1. Review Confinement & Single Paticle Motion. Plasmas as Fluids Fluid Equations 3. MHD Equations
More informationSupplementary Figure 1. Circular parallel lamellae grain size as a function of annealing time at 250 C. Error bars represent the 2σ uncertainty in
Supplementay Figue 1. Cicula paallel lamellae gain size as a function of annealing time at 50 C. Eo bas epesent the σ uncetainty in the measued adii based on image pixilation and analysis uncetainty contibutions
More information1 Spherical multipole moments
Jackson notes 9 Spheical multipole moments Suppose we have a chage distibution ρ (x) wheeallofthechageiscontained within a spheical egion of adius R, as shown in the diagam. Then thee is no chage in the
More information3-Axis Vector Magnet: Construction and Characterisation of Split Coils at RT. Semester Project Petar Jurcevic
3-Axis Vecto Magnet: Constuction and Chaacteisation of Split Coils at RT Semeste Poject Peta Jucevic Outline Field Calculation and Simulation Constuction Details Field Calculations Chaacteization at RT
More informationEM Boundary Value Problems
EM Bounday Value Poblems 10/ 9 11/ By Ilekta chistidi & Lee, Seung-Hyun A. Geneal Desciption : Maxwell Equations & Loentz Foce We want to find the equations of motion of chaged paticles. The way to do
More informationTHE LAPLACE EQUATION. The Laplace (or potential) equation is the equation. u = 0. = 2 x 2. x y 2 in R 2
THE LAPLACE EQUATION The Laplace (o potential) equation is the equation whee is the Laplace opeato = 2 x 2 u = 0. in R = 2 x 2 + 2 y 2 in R 2 = 2 x 2 + 2 y 2 + 2 z 2 in R 3 The solutions u of the Laplace
More informationAP-C WEP. h. Students should be able to recognize and solve problems that call for application both of conservation of energy and Newton s Laws.
AP-C WEP 1. Wok a. Calculate the wok done by a specified constant foce on an object that undegoes a specified displacement. b. Relate the wok done by a foce to the aea unde a gaph of foce as a function
More informationMath 2263 Solutions for Spring 2003 Final Exam
Math 6 Solutions fo Sping Final Exam ) A staightfowad appoach to finding the tangent plane to a suface at a point ( x, y, z ) would be to expess the cuve as an explicit function z = f ( x, y ), calculate
More informationGoodness-of-fit for composite hypotheses.
Section 11 Goodness-of-fit fo composite hypotheses. Example. Let us conside a Matlab example. Let us geneate 50 obsevations fom N(1, 2): X=nomnd(1,2,50,1); Then, unning a chi-squaed goodness-of-fit test
More informationPhysics 2A Chapter 10 - Moment of Inertia Fall 2018
Physics Chapte 0 - oment of netia Fall 08 The moment of inetia of a otating object is a measue of its otational inetia in the same way that the mass of an object is a measue of its inetia fo linea motion.
More informationworking pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50
woking pages fo Paul Richads class notes; do not copy o ciculate without pemission fom PGR 2004/11/3 10:50 CHAPTER7 Solid angle, 3D integals, Gauss s Theoem, and a Delta Function We define the solid angle,
More informationHopefully Helpful Hints for Gauss s Law
Hopefully Helpful Hints fo Gauss s Law As befoe, thee ae things you need to know about Gauss s Law. In no paticula ode, they ae: a.) In the context of Gauss s Law, at a diffeential level, the electic flux
More informationLight Time Delay and Apparent Position
Light Time Delay and ppaent Position nalytical Gaphics, Inc. www.agi.com info@agi.com 610.981.8000 800.220.4785 Contents Intoduction... 3 Computing Light Time Delay... 3 Tansmission fom to... 4 Reception
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Seies UG Examination 2015 16 FLUID DYNAMICS WITH ADVANCED TOPICS MTH-MD59 Time allowed: 3 Hous Attempt QUESTIONS 1 and 2, and THREE othe questions.
More information22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 4: Toroidal Equilibrium and Radial Pressure Balance
.615, MHD Theoy of Fusion Systems Pof. Feidbeg Lectue 4: Tooidal Equilibium and Radial Pessue Balance Basic Poblem of Tooidal Equilibium 1. Radial pessue balance. Tooidal foce balance Radial Pessue Balance
More information8 Separation of Variables in Other Coordinate Systems
8 Sepaation of Vaiables in Othe Coodinate Systems Fo the method of sepaation of vaiables to succeed you need to be able to expess the poblem at hand in a coodinate system in which the physical boundaies
More informationA dual-reciprocity boundary element method for axisymmetric thermoelastodynamic deformations in functionally graded solids
APCOM & ISCM 11-14 th Decembe, 013, Singapoe A dual-ecipocity bounday element method fo axisymmetic themoelastodynamic defomations in functionally gaded solids *W. T. Ang and B. I. Yun Division of Engineeing
More informationr cos, and y r sin with the origin of coordinate system located at
Lectue 3-3 Kinematics of Rotation Duing ou peious lectues we hae consideed diffeent examples of motion in one and seeal dimensions. But in each case the moing object was consideed as a paticle-like object,
More informationHigh precision computer simulation of cyclotrons KARAMYSHEVA T., AMIRKHANOV I. MALININ V., POPOV D.
High pecision compute simulation of cyclotons KARAMYSHEVA T., AMIRKHANOV I. MALININ V., POPOV D. Abstact Effective and accuate compute simulations ae highly impotant in acceleatos design and poduction.
More informationElectron Density Distribution in HSX
Electon Density Distibution in HSX C. Deng, D.L. Bowe Electical Engineeing Depatment Univesity of Califonia, Los Angeles J. Canik, S.P. Gehadt, D.T. Andeson, P. Pobet, F.S.B. Andeson The HSX Plasma Laboatoy
More informationSchool of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007
School of Electical and Compute Engineeing, Conell Univesity ECE 303: Electomagnetic Fields and Waves Fall 007 Homewok 8 Due on Oct. 19, 007 by 5:00 PM Reading Assignments: i) Review the lectue notes.
More informationTo Feel a Force Chapter 7 Static equilibrium - torque and friction
To eel a oce Chapte 7 Chapte 7: Static fiction, toque and static equilibium A. Review of foce vectos Between the eath and a small mass, gavitational foces of equal magnitude and opposite diection act on
More information-Δ u = λ u. u(x,y) = u 1. (x) u 2. (y) u(r,θ) = R(r) Θ(θ) Δu = 2 u + 2 u. r = x 2 + y 2. tan(θ) = y/x. r cos(θ) = cos(θ) r.
The Laplace opeato in pola coodinates We now conside the Laplace opeato with Diichlet bounday conditions on a cicula egion Ω {(x,y) x + y A }. Ou goal is to compute eigenvalues and eigenfunctions of the
More informationNewton s Laws, Kepler s Laws, and Planetary Orbits
Newton s Laws, Keple s Laws, and Planetay Obits PROBLEM SET 4 DUE TUESDAY AT START OF LECTURE 28 Septembe 2017 ASTRONOMY 111 FALL 2017 1 Newton s & Keple s laws and planetay obits Unifom cicula motion
More information= e2. = 2e2. = 3e2. V = Ze2. where Z is the atomic numnber. Thus, we take as the Hamiltonian for a hydrogenic. H = p2 r. (19.4)
Chapte 9 Hydogen Atom I What is H int? That depends on the physical system and the accuacy with which it is descibed. A natual stating point is the fom H int = p + V, (9.) µ which descibes a two-paticle
More information2 Governing Equations
2 Govening Equations This chapte develops the govening equations of motion fo a homogeneous isotopic elastic solid, using the linea thee-dimensional theoy of elasticity in cylindical coodinates. At fist,
More informationBut for simplicity, we ll define significant as the time it takes a star to lose all memory of its original trajectory, i.e.,
Stella elaxation Time [Chandasekha 1960, Pinciples of Stella Dynamics, Chap II] [Ostike & Davidson 1968, Ap.J., 151, 679] Do stas eve collide? Ae inteactions between stas (as opposed to the geneal system
More information4/18/2005. Statistical Learning Theory
Statistical Leaning Theoy Statistical Leaning Theoy A model of supevised leaning consists of: a Envionment - Supplying a vecto x with a fixed but unknown pdf F x (x b Teache. It povides a desied esponse
More informationMATH 220: SECOND ORDER CONSTANT COEFFICIENT PDE. We consider second order constant coefficient scalar linear PDEs on R n. These have the form
MATH 220: SECOND ORDER CONSTANT COEFFICIENT PDE ANDRAS VASY We conside second ode constant coefficient scala linea PDEs on R n. These have the fom Lu = f L = a ij xi xj + b i xi + c i whee a ij b i and
More informationContact impedance of grounded and capacitive electrodes
Abstact Contact impedance of gounded and capacitive electodes Andeas Hödt Institut fü Geophysik und extateestische Physik, TU Baunschweig The contact impedance of electodes detemines how much cuent can
More information1 Fundamental Solutions to the Wave Equation
1 Fundamental Solutions to the Wave Equation Physical insight in the sound geneation mechanism can be gained by consideing simple analytical solutions to the wave equation. One example is to conside acoustic
More informationThe geometric construction of Ewald sphere and Bragg condition:
The geometic constuction of Ewald sphee and Bagg condition: The constuction of Ewald sphee must be done such that the Bagg condition is satisfied. This can be done as follows: i) Daw a wave vecto k in
More informationDesigning a Sine-Coil for Measurement of Plasma Displacements in IR-T1 Tokamak
Designing a Sine-Coil fo Measuement of Plasma Displacements in IR-T Tokamak Pejman Khoshid, M. Razavi, M. Ghoanneviss, M. Molaii, A. TalebiTahe, R. Avin, S. Mohammadi and A. NikMohammadi Dept. of Physics,
More informationCBE Transport Phenomena I Final Exam. December 19, 2013
CBE 30355 Tanspot Phenomena I Final Exam Decembe 9, 203 Closed Books and Notes Poblem. (20 points) Scaling analysis of bounday laye flows. A popula method fo measuing instantaneous wall shea stesses in
More informationPhysics 2212 GH Quiz #2 Solutions Spring 2016
Physics 2212 GH Quiz #2 Solutions Sping 216 I. 17 points) Thee point chages, each caying a chage Q = +6. nc, ae placed on an equilateal tiangle of side length = 3. mm. An additional point chage, caying
More informationScientific Computing II
Scientific Computing II Conjugate Gadient Methods Michael Bade Summe 2014 Conjugate Gadient Methods, Summe 2014 1 Families of Iteative Solves elaxation methods: Jacobi-, Gauss-Seidel-Relaxation,... Ove-Relaxation-Methods
More information2. Plane Elasticity Problems
S0 Solid Mechanics Fall 009. Plane lasticity Poblems Main Refeence: Theoy of lasticity by S.P. Timoshenko and J.N. Goodie McGaw-Hill New Yok. Chaptes 3..1 The plane-stess poblem A thin sheet of an isotopic
More informationElectromagnetic scattering. Graduate Course Electrical Engineering (Communications) 1 st Semester, Sharif University of Technology
Electomagnetic scatteing Gaduate Couse Electical Engineeing (Communications) 1 st Semeste, 1390-1391 Shaif Univesity of Technology Geneal infomation Infomation about the instucto: Instucto: Behzad Rejaei
More informationLecture 5 Solving Problems using Green s Theorem. 1. Show how Green s theorem can be used to solve general electrostatic problems 2.
Lectue 5 Solving Poblems using Geen s Theoem Today s topics. Show how Geen s theoem can be used to solve geneal electostatic poblems. Dielectics A well known application of Geen s theoem. Last time we
More informationQualifying Examination Electricity and Magnetism Solutions January 12, 2006
1 Qualifying Examination Electicity and Magnetism Solutions Januay 12, 2006 PROBLEM EA. a. Fist, we conside a unit length of cylinde to find the elationship between the total chage pe unit length λ and
More informationInternet Appendix for A Bayesian Approach to Real Options: The Case of Distinguishing Between Temporary and Permanent Shocks
Intenet Appendix fo A Bayesian Appoach to Real Options: The Case of Distinguishing Between Tempoay and Pemanent Shocks Steven R. Genadie Gaduate School of Business, Stanfod Univesity Andey Malenko Gaduate
More informationDymore User s Manual Two- and three dimensional dynamic inflow models
Dymoe Use s Manual Two- and thee dimensional dynamic inflow models Contents 1 Two-dimensional finite-state genealized dynamic wake theoy 1 Thee-dimensional finite-state genealized dynamic wake theoy 1
More informationRotational Motion. Lecture 6. Chapter 4. Physics I. Course website:
Lectue 6 Chapte 4 Physics I Rotational Motion Couse website: http://faculty.uml.edu/andiy_danylov/teaching/physicsi Today we ae going to discuss: Chapte 4: Unifom Cicula Motion: Section 4.4 Nonunifom Cicula
More informationReconstruction of the Eddy Current Distribution on the Vacuum Vessel in a Reversed Field Pinch Device based on the External Magnetic Sensor Signals
Reconstuction of the Eddy Cuent Distibution on the Vacuum Vessel in a Revesed Field Pinch Device based on the Extenal Magnetic Senso Signals Masafumi ITAGAKI, Akio SANPEI ), Sadao MASAMUNE ) and Kiyomasa
More informationEddy Currents and Magnetic Calibrations in LDX using a Copper Plasma. D.P. Boyle, PPPL M.E. Mauel, D.T. Garnier, Columbia J.
Eddy Cuents and Magneti Calibations in LDX using a Coppe Plasma D.P. Boyle PPPL M.E. Mauel D.T. Ganie Columbia J. Kesne MIT PSFC Coppe Plasma Oveview LDX Magnetis Goals Calibate magneti diagnostis positions
More informationInseting this into the left hand side of the equation of motion above gives the most commonly used algoithm in classical molecula dynamics simulations
Chem465 in 2000 Univesity of Washington Lectue notes Hannes Jonsson Classical dynamics When we ae dealing with heavy atoms and high enough enegy o tempeatue, it is often suciently accuate to neglect quantum
More informationIn the previous section we considered problems where the
5.4 Hydodynamically Fully Developed and Themally Developing Lamina Flow In the pevious section we consideed poblems whee the velocity and tempeatue pofile wee fully developed, so that the heat tansfe coefficient
More informationASTR415: Problem Set #6
ASTR45: Poblem Set #6 Cuan D. Muhlbege Univesity of Mayland (Dated: May 7, 27) Using existing implementations of the leapfog and Runge-Kutta methods fo solving coupled odinay diffeential equations, seveal
More informationGA A22703 POLOIDAL GRADIENTS, DIVERTOR DETACHMENT AND MARFES
GA A22703 POLOIDAL GRADIENTS, DIVERTOR DETACHMENT by NOVEMBER 1997 DISCLAIMER This epot was pepaed as an account of wok sponsoed by an agency of the United States Govenment. Neithe the United States Govenment
More informationOn a quantity that is analogous to potential and a theorem that relates to it
Su une quantité analogue au potential et su un théoème y elatif C R Acad Sci 7 (87) 34-39 On a quantity that is analogous to potential and a theoem that elates to it By R CLAUSIUS Tanslated by D H Delphenich
More informationLecture 23. Representation of the Dirac delta function in other coordinate systems
Lectue 23 Repesentation of the Diac delta function in othe coodinate systems In a geneal sense, one can wite, ( ) = (x x ) (y y ) (z z ) = (u u ) (v v ) (w w ) J Whee J epesents the Jacobian of the tansfomation.
More informationCh 30 - Sources of Magnetic Field! The Biot-Savart Law! = k m. r 2. Example 1! Example 2!
Ch 30 - Souces of Magnetic Field 1.) Example 1 Detemine the magnitude and diection of the magnetic field at the point O in the diagam. (Cuent flows fom top to bottom, adius of cuvatue.) Fo staight segments,
More information( )( )( ) ( ) + ( ) ( ) ( )
3.7. Moel: The magnetic fiel is that of a moving chage paticle. Please efe to Figue Ex3.7. Solve: Using the iot-savat law, 7 19 7 ( ) + ( ) qvsinθ 1 T m/a 1.6 1 C. 1 m/s sin135 1. 1 m 1. 1 m 15 = = = 1.13
More information= 4 3 π( m) 3 (5480 kg m 3 ) = kg.
CHAPTER 11 THE GRAVITATIONAL FIELD Newton s Law of Gavitation m 1 m A foce of attaction occus between two masses given by Newton s Law of Gavitation Inetial mass and gavitational mass Gavitational potential
More informationNumerical Methods for Two-Fluid Dispersive Fast MHD Phenomena
Numeical Methods fo Two-Fluid Dispesive Fast MHD Phenomena B. Sinivasan 1,, A. Hakim 2, and U. Shumlak 3 1 Applied Mathematics and Plasma Physics, Los Alamos National Laboatoy, Los Alamos, NM 87545 2 Tech-X
More information0606 ADDITIONAL MATHEMATICS
UNIVERSITY OF CAMBRIDGE INTERNATIONAL EXAMINATIONS Intenational Geneal Cetificate of Seconday Education MARK SCHEME fo the Octobe/Novembe 011 question pape fo the guidance of teaches 0606 ADDITIONAL MATHEMATICS
More informationDiffusion and Transport. 10. Friction and the Langevin Equation. Langevin Equation. f d. f ext. f () t f () t. Then Newton s second law is ma f f f t.
Diffusion and Tanspot 10. Fiction and the Langevin Equation Now let s elate the phenomena of ownian motion and diffusion to the concept of fiction, i.e., the esistance to movement that the paticle in the
More informationUniversity of Illinois at Chicago Department of Physics. Electricity & Magnetism Qualifying Examination
E&M poblems Univesity of Illinois at Chicago Depatment of Physics Electicity & Magnetism Qualifying Examination Januay 3, 6 9. am : pm Full cedit can be achieved fom completely coect answes to 4 questions.
More informationtransformation Earth V-curve (meridian) λ Conical projection. u,v curves on the datum surface projected as U,V curves on the projection surface
. CONICAL PROJECTIONS In elementay texts on map pojections, the pojection sufaces ae often descibed as developable sufaces, such as the cylinde (cylindical pojections) and the cone (conical pojections),
More informationECE 3318 Applied Electricity and Magnetism. Spring Prof. David R. Jackson ECE Dept. Notes 13
ECE 338 Applied Electicity and Magnetism ping 07 Pof. David R. Jackson ECE Dept. Notes 3 Divegence The Physical Concept Find the flux going outwad though a sphee of adius. x ρ v0 z a y ψ = D nˆ d = D ˆ
More informationMONTE CARLO STUDY OF PARTICLE TRANSPORT PROBLEM IN AIR POLLUTION. R. J. Papancheva, T. V. Gurov, I. T. Dimov
Pliska Stud. Math. Bulga. 14 (23), 17 116 STUDIA MATHEMATICA BULGARICA MOTE CARLO STUDY OF PARTICLE TRASPORT PROBLEM I AIR POLLUTIO R. J. Papancheva, T. V. Guov, I. T. Dimov Abstact. The actual tanspot
More informationA Most Useful Device of Studying Electrode Processes: The Rotating Disk Electrode
A Most Useful Device of Studying Electode Pocesses: The Rotating Disk Electode the theoetical basis Soma Vesztegom Laboatoy of Electochemisty & Electoanalytical Chemisty Eötvös Loánd Univesity of Budapest
More informationNumerical Integration
MCEN 473/573 Chapte 0 Numeical Integation Fall, 2006 Textbook, 0.4 and 0.5 Isopaametic Fomula Numeical Integation [] e [ ] T k = h B [ D][ B] e B Jdsdt In pactice, the element stiffness is calculated numeically.
More informationDOING PHYSICS WITH MATLAB COMPUTATIONAL OPTICS
DOING PHYIC WITH MTLB COMPUTTIONL OPTIC FOUNDTION OF CLR DIFFRCTION THEORY Ian Coope chool of Physics, Univesity of ydney ian.coope@sydney.edu.au DOWNLOD DIRECTORY FOR MTLB CRIPT View document: Numeical
More informationPearson s Chi-Square Test Modifications for Comparison of Unweighted and Weighted Histograms and Two Weighted Histograms
Peason s Chi-Squae Test Modifications fo Compaison of Unweighted and Weighted Histogams and Two Weighted Histogams Univesity of Akueyi, Bogi, v/noduslód, IS-6 Akueyi, Iceland E-mail: nikolai@unak.is Two
More informationStanford University CS259Q: Quantum Computing Handout 8 Luca Trevisan October 18, 2012
Stanfod Univesity CS59Q: Quantum Computing Handout 8 Luca Tevisan Octobe 8, 0 Lectue 8 In which we use the quantum Fouie tansfom to solve the peiod-finding poblem. The Peiod Finding Poblem Let f : {0,...,
More informationNuclear size corrections to the energy levels of single-electron atoms
Nuclea size coections to the enegy levels of single-electon atoms Babak Nadii Nii a eseach Institute fo Astonomy and Astophysics of Maagha (IAAM IAN P. O. Box: 554-44. Abstact A study is made of nuclea
More informationA Bijective Approach to the Permutational Power of a Priority Queue
A Bijective Appoach to the Pemutational Powe of a Pioity Queue Ia M. Gessel Kuang-Yeh Wang Depatment of Mathematics Bandeis Univesity Waltham, MA 02254-9110 Abstact A pioity queue tansfoms an input pemutation
More informationThree dimensional flow analysis in Axial Flow Compressors
1 Thee dimensional flow analysis in Axial Flow Compessos 2 The ealie assumption on blade flow theoies that the flow inside the axial flow compesso annulus is two dimensional means that adial movement of
More informationInformation Retrieval Advanced IR models. Luca Bondi
Advanced IR models Luca Bondi Advanced IR models 2 (LSI) Pobabilistic Latent Semantic Analysis (plsa) Vecto Space Model 3 Stating point: Vecto Space Model Documents and queies epesented as vectos in the
More information