{"id":3192,"date":"2014-04-22T19:02:11","date_gmt":"2014-04-22T23:02:11","guid":{"rendered":"http:\/\/pages.vassar.edu\/magnes\/?p=3192"},"modified":"2014-04-29T18:41:44","modified_gmt":"2014-04-29T22:41:44","slug":"preliminary-results-modeling-e-and-b-fields-of-a-toroidal-conductor-to-explore-tokamak-plasmas","status":"publish","type":"post","link":"https:\/\/pages.vassar.edu\/magnes\/2014\/04\/22\/preliminary-results-modeling-e-and-b-fields-of-a-toroidal-conductor-to-explore-tokamak-plasmas\/","title":{"rendered":"Preliminary Results: Modeling E and B Fields of a Toroidal Conductor to Explore Tokamak Plasmas"},"content":{"rendered":"<p>Thus far, I have determined values for the quantities of interest (except for $\\mu$ of the plasma) and I have derived an expression for\u00a0<strong>B\u00a0<\/strong>and\u00a0<strong>H<\/strong>. \u00a0I will base the parameters of my model off of the Joint European Torus (as per Albanese, et al., 03; Crisanti, et al., 03). It appears that there will be no\u00a0<strong>D\u00a0<\/strong>field, since the plasma is a conductor and would not polarize, but simply experience a current.<\/p>\n<ul>\n<li>$q(r) = \\frac{r B_{\\phi}(r)}{R B_{\\theta}(r)} = 3$<\/li>\n<li>$r = .95~m$<\/li>\n<li>$R = 2.85~m$<\/li>\n<li>$I = 6.0~MA$<\/li>\n<\/ul>\n<p>I am going to assume that N, the total number of turns, is proportional to the safety factor: $N = \\alpha q(r)$ where $\\alpha$ will probably be related to the circumference of the torus (which can be found from the radii).<\/p>\n<p>The Biot-Savart Law gives:<\/p>\n<p>$d\\textbf{B} = \\frac{\\mu_{0}\u00a0\\textbf{I} \\times \\textbf{r} dl}{4 \\pi r^{3}}$.<\/p>\n<p>If we assume that the torus is in the x-y plane, centered at the origin, and consider a test point at $\\textbf{r}$ and source charge elements at various $\\textbf{r&#8217;}$. \u00a0Then $\\scriptr =\u00a0\\textbf{r} &#8211; \\textbf{r&#8217;}$. \u00a0Consider a test point on the x-z plane.<\/p>\n<p>Assuming that the winding is tight enough, we can say that $ \\textbf{I} = I_{s} \\hat{s} +\u00a0I_{z} \\hat{z}$ which, in cylindrical coordinates, is $(I_{s}cos(\\phi), I_{s}sin(\\phi), I_{z})$.<\/p>\n<p>In evaluating the Biot-Savart Law, we have the term $\\textbf{I} \\times \\textbf{r}$. \u00a0Due to the symmetry of the tours, the x and z components cancel out, leaving just a y component. \u00a0If the test point is in the x-z plane, this means that the magnetic field will be circumferential, or in the positive $\\hat{\\phi}$ direction. \u00a0The magnetic field can then easily be found using Ampere&#8217;s Law:<\/p>\n<p>$\\textbf{B(r)} = \\frac{\\mu_{0} N I}{2 \\pi s} \\hat{\\phi}$<\/p>\n<p>for all points within the minor radius (inside the torus) and<\/p>\n<p>$\\textbf{B(r)} = 0$<\/p>\n<p>for all points outside the torus.<\/p>\n<p>From here, it is also easy to find the Auxiliary field, from $\\textbf{H} = (1\/\\mu) \\textbf{B}$, so<\/p>\n<p>$\\textbf{H} = \\frac{1}{\\mu}\u00a0\\frac{\\mu_{0} N I}{2 \\pi s} \\hat{\\phi}$<\/p>\n<p>inside the torus and<\/p>\n<p>$\\textbf{H} = 0$<\/p>\n<p>outside the torus, since there is no magnetic field.<\/p>\n<p>I have also begun work on my Mathematica model. \u00a0I began by trying to figure out how to get the vector field to look right, and have been working with a simplified version of my\u00a0<strong>B\u00a0<\/strong>equation, $\\textbf{B}\u00a0= k\/r$. \u00a0I have found that the transform to spherical coordinates is easiest to work with (since toroidal geometries do not fit easily into either spherical or cylindrical), which is why my expression uses r instead of s. \u00a0This has given the following vector field:<\/p>\n<div id=\"attachment_3193\" style=\"width: 310px\" class=\"wp-caption aligncenter\"><a href=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/B_of-_torus-_-2014_-04_-22.jpeg\"><img loading=\"lazy\" decoding=\"async\" aria-describedby=\"caption-attachment-3193\" class=\"size-medium wp-image-3193\" alt=\"The (simplified) circumferential magnetic field due to a toroidal solenoid (without taking into account the fact that the field is zero outside of the loop of the solenoid).\" src=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/B_of-_torus-_-2014_-04_-22-300x295.jpeg\" width=\"300\" height=\"295\" srcset=\"https:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/B_of-_torus-_-2014_-04_-22-300x295.jpeg 300w, https:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/B_of-_torus-_-2014_-04_-22.jpeg 360w\" sizes=\"auto, (max-width: 300px) 100vw, 300px\" \/><\/a><p id=\"caption-attachment-3193\" class=\"wp-caption-text\">The (simplified) circumferential magnetic field due to a toroidal solenoid (without taking into account the fact that the field is zero outside of the loop of the solenoid).<\/p><\/div>\n<p>The vector field is the shape that I need it to be, but I&#8217;m having some trouble getting the field to be zero outside of the solenoid. \u00a0I tried to use the piecewise function of Mathematica, but keep getting an error. \u00a0I think I may need to define two fields (maybe 3) such that they cancel out in the regions where the field should be zero.<\/p>\n<p>The next step will be to input the numbers into my expression, and work out the Mathematica model. \u00a0Then I will be able to easily vary\u00a0<em>q(r)<\/em> and observe the results.<\/p>\n<p>Link to Mathematica notebook:\u00a0<a href=\"https:\/\/drive.google.com\/file\/d\/0B-C9MvBAfmyIQS12dTZJOGh5bm8\/edit?usp=sharing\">https:\/\/drive.google.com\/file\/d\/0B-C9MvBAfmyIQS12dTZJOGh5bm8\/edit?usp=sharing<\/a><\/p>\n<p><strong>References<\/strong><\/p>\n<p>Ra\ufb00aele Albanese, G Calabr`o, M Mattei, and F Villone. Plasma response models<br \/>\nfor current, shape and position control in jet. Fusion engineering and design,<br \/>\n66:715\u2013718, 2003.<\/p>\n<p>F Crisanti, R Albanese, G Ambrosino, M Ariola, J Lister, M Mattei, F Milani,<br \/>\nA Pironti, F Sartori, and F Villone. Upgrade of the present jet shape and vertical<br \/>\nstability controller. Fusion engineering and design, 66:803\u2013807, 2003.<\/p>\n<p>Griffiths, D. J., &amp; Reed College. (1999).\u00a0<i>Introduction to electrodynamics<\/i>\u00a0(Vol. 3). Upper Saddle River, NJ: prentice Hall.<\/p>\n","protected":false},"excerpt":{"rendered":"<p>Thus far, I have determined values for the quantities of interest (except for $\\mu$ of the plasma) and I have derived an expression for\u00a0B\u00a0and\u00a0H. \u00a0I will base the parameters of my model off of the Joint European Torus (as per Albanese, et al., 03; Crisanti, et al., 03). It appears that there will be no\u00a0D\u00a0field, [&hellip;]<\/p>\n","protected":false},"author":1773,"featured_media":0,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[4101,54056,54190],"tags":[],"class_list":["post-3192","post","type-post","status-publish","format-standard","hentry","category-advanced-em","category-derek","category-spring-2014"],"_links":{"self":[{"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/posts\/3192","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/users\/1773"}],"replies":[{"embeddable":true,"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/comments?post=3192"}],"version-history":[{"count":23,"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/posts\/3192\/revisions"}],"predecessor-version":[{"id":3572,"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/posts\/3192\/revisions\/3572"}],"wp:attachment":[{"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/media?parent=3192"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/categories?post=3192"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/tags?post=3192"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}