{"id":3608,"date":"2014-04-30T03:32:07","date_gmt":"2014-04-30T07:32:07","guid":{"rendered":"http:\/\/pages.vassar.edu\/magnes\/?p=3608"},"modified":"2014-04-30T05:00:17","modified_gmt":"2014-04-30T09:00:17","slug":"final-data-rlc-circuits","status":"publish","type":"post","link":"https:\/\/pages.vassar.edu\/magnes\/2014\/04\/30\/final-data-rlc-circuits\/","title":{"rendered":"Final Data- RLC Circuits"},"content":{"rendered":"<p align=\"JUSTIFY\"><strong>Overview<\/strong><\/p>\n<p align=\"JUSTIFY\">After talking with Professor Magnes after my preliminary results, I decided to first look at an RLC circuit with no voltage source, which is much simpler to model. Doing this allowed me to better understand the mechanism behind the circuitry and gain more experience with Mathematica as well.<\/p>\n<p align=\"JUSTIFY\">After doing so, I returned to my progress with an RLC circuit with a voltage source and determined that the odd structural variations with the long-term structure were due to computational limitations of Mathematica (again, thanks to Professor Magnes for the suggestion). Since I had only developed the approximate solution, I went on to graphically solve for the full solution and see what the damping terms did to the solution.<\/p>\n<p align=\"JUSTIFY\">Note: I included the findings from my preliminary data in this final data post (so it may seem redundant for a portion in the middle) since it was used in my final results as well. I figured it helps to have everything in one place for the reader.<\/p>\n<p style=\"text-align: center\" align=\"JUSTIFY\"><span style=\"text-decoration: underline\"><strong>RLC Circuit with no voltage source<\/strong><\/span><\/p>\n<p style=\"text-align: left\" align=\"JUSTIFY\"><strong>Circuit Diagram and Mathematics<\/strong><\/p>\n<p style=\"text-align: left\" align=\"JUSTIFY\">The first circuit we will examine is an RLC circuit with no voltage source. The initial potential difference is stored on the capacitor, but can not be recharged due to the lack of a voltage source. The three main components of this circuit are a resistor (R), an inductor (L), and a capacitor (C).<\/p>\n<p align=\"JUSTIFY\"><a href=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/1.png\"><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-3637\" alt=\"1\" src=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/1.png\" width=\"242\" height=\"199\" \/><\/a>[Figure 1]<\/p>\n<p align=\"JUSTIFY\">The differential equation that governs this circuit is<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-dff481c7d798be95ebe92a3dcd18af4e_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#76;&#73;&#39;&#39;&#40;&#116;&#41;&#43;&#82;&#73;&#39;&#40;&#116;&#41;&#43;&#92;&#102;&#114;&#97;&#99;&#123;&#49;&#125;&#123;&#67;&#125;&#73;&#40;&#116;&#41;&#61;&#32;&#48;\" title=\"Rendered by QuickLaTeX.com\" height=\"22\" width=\"217\" style=\"vertical-align: -6px;\"\/> [1]<\/p>\n<p>where L is the inductance, R is the resistance, C is the capacitance, I is the current through the circuit (which is the same at all places at one time since this is a series circuit), and t is the time.<\/p>\n<p><strong>Solving the Differential Equation<\/strong><\/p>\n<p>We can set arbitrary values for R, L, and C since they depend on the specific pieces of equipment themselves:<\/p>\n<ul>\n<li>Inductance= 1 Henry<\/li>\n<li>Resistance= 10 ohms<\/li>\n<li>Capacitance= .1 farads<\/li>\n<\/ul>\n<p>In Mathematica, the command to solve [1] looks like<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-edf43b2ac32669349aa6c8d04ba9e32e_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#92;&#116;&#101;&#120;&#116;&#123;&#68;&#83;&#111;&#108;&#118;&#101;&#125;&#92;&#108;&#101;&#102;&#116;&#91;&#92;&#102;&#114;&#97;&#99;&#123;&#92;&#116;&#101;&#120;&#116;&#123;&#67;&#117;&#114;&#114;&#101;&#110;&#116;&#125;&#40;&#116;&#41;&#125;&#123;&#92;&#116;&#101;&#120;&#116;&#123;&#67;&#97;&#112;&#97;&#99;&#105;&#116;&#97;&#110;&#99;&#101;&#125;&#125;&#43;&#92;&#116;&#101;&#120;&#116;&#123;&#73;&#110;&#100;&#117;&#99;&#116;&#97;&#110;&#99;&#101;&#125;&#32;&#92;&#116;&#101;&#120;&#116;&#123;&#67;&#117;&#114;&#114;&#101;&#110;&#116;&#125;&#39;&#39;&#40;&#116;&#41;&#43;&#92;&#116;&#101;&#120;&#116;&#123;&#82;&#101;&#115;&#105;&#115;&#116;&#97;&#110;&#99;&#101;&#125;&#32;&#92;&#116;&#101;&#120;&#116;&#123;&#67;&#117;&#114;&#114;&#101;&#110;&#116;&#125;&#39;&#40;&#116;&#41;&#61;&#48;&#44;&#92;&#116;&#101;&#120;&#116;&#123;&#67;&#117;&#114;&#114;&#101;&#110;&#116;&#125;&#40;&#116;&#41;&#44;&#116;&#92;&#114;&#105;&#103;&#104;&#116;&#93;\" title=\"Rendered by QuickLaTeX.com\" height=\"32\" width=\"667\" style=\"vertical-align: -11px;\"\/> [2]<\/p>\n<p>which when solved, yields:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-13f25fdefef2d4e488eccdb4ff9cc983_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#73;&#40;&#116;&#41;&#61;&#32;&#99;&#95;&#49;&#32;&#101;&#94;&#123;&#45;&#56;&#46;&#56;&#55;&#50;&#57;&#56;&#32;&#116;&#125;&#43;&#99;&#95;&#50;&#32;&#101;&#94;&#123;&#45;&#49;&#46;&#49;&#50;&#55;&#48;&#50;&#32;&#116;&#125;\" title=\"Rendered by QuickLaTeX.com\" height=\"19\" width=\"245\" style=\"vertical-align: -4px;\"\/> [3]<\/p>\n<p><strong>Completing the Full Solution via Graphical Solving<\/strong><\/p>\n<p>The two constants depend on the initial conditions of the circuit, which must satisfy Ohm&#8217;s Law of <img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-0d913679a42337de0fe847ab1e3c2d62_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#73;&#61;&#92;&#102;&#114;&#97;&#99;&#123;&#86;&#125;&#123;&#82;&#125;\" title=\"Rendered by QuickLaTeX.com\" height=\"22\" width=\"47\" style=\"vertical-align: -6px;\"\/>. So, the initial current in the circuit must be 0.1 amps. Using a graphical guess-and-check method in Mathematica, I determined the value of the two constants to be approximately 0.05 each, resulting in a solution that is given by<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-374d28dbecf617b5a36f9d8b03c8d888_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#73;&#40;&#116;&#41;&#61;&#32;&#48;&#46;&#48;&#53;&#32;&#101;&#94;&#123;&#45;&#56;&#46;&#56;&#55;&#50;&#57;&#56;&#32;&#116;&#125;&#43;&#48;&#46;&#48;&#53;&#32;&#101;&#94;&#123;&#45;&#49;&#46;&#49;&#50;&#55;&#48;&#50;&#32;&#116;&#125;\" title=\"Rendered by QuickLaTeX.com\" height=\"19\" width=\"278\" style=\"vertical-align: -4px;\"\/> [4]<\/p>\n<p>The plot of this equation gives us the current as a function of time, and looks like<\/p>\n<p><a href=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata11.jpg\"><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-3652\" alt=\"RLCCircuitsworking_finaldata1\" src=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata11.jpg\" width=\"396\" height=\"269\" srcset=\"https:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata11.jpg 396w, https:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata11-300x203.jpg 300w\" sizes=\"auto, (max-width: 396px) 100vw, 396px\" \/><\/a>[Figure 2]<\/p>\n<p>This graph tells us that although the RLC circuit starts out with a current of 0.1 amps, it quickly dies down due to the lack of a voltage source and is not very interesting after a few seconds. Refer to the upcoming conclusion blog post for a bit more analysis.<\/p>\n<p>Now, lets examine the RLC circuit with a voltage source (much more interesting!)<\/p>\n<p style=\"text-align: center\" align=\"JUSTIFY\"><span style=\"text-decoration: underline\"><strong>RLC Circuit with AC voltage source<\/strong><\/span><\/p>\n<p style=\"text-align: left\" align=\"JUSTIFY\"><strong>Circuit Diagram and Mathematics<\/strong><\/p>\n<p align=\"JUSTIFY\">Now we can model an RLC circuit that consists of a capacitor, voltage source, resistor, and an inductor.\u00a0To easily visualize this, I have constructed a basic circuit diagram (Figure 3).<\/p>\n<p align=\"JUSTIFY\"><a href=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/circuitwithnumbers11.png\"><img loading=\"lazy\" decoding=\"async\" alt=\"circuitwithnumbers1\" src=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/circuitwithnumbers11.png\" width=\"242\" height=\"199\" \/><\/a>(Figure 3)<\/p>\n<p align=\"JUSTIFY\">The differential equation that governs this RLC circuit is given by<\/p>\n<p align=\"JUSTIFY\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-eac0cedfb7457178f885cdcf95c0db3d_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#76;&#73;&#39;&#39;&#40;&#116;&#41;&#43;&#82;&#73;&#39;&#40;&#116;&#41;&#43;&#92;&#102;&#114;&#97;&#99;&#123;&#49;&#125;&#123;&#67;&#125;&#73;&#40;&#116;&#41;&#61;&#32;&#69;&#95;&#123;&#48;&#125;&#32;&#92;&#111;&#109;&#101;&#103;&#97;&#32;&#99;&#111;&#115;&#40;&#92;&#111;&#109;&#101;&#103;&#97;&#32;&#116;&#41;\" title=\"Rendered by QuickLaTeX.com\" height=\"22\" width=\"296\" style=\"vertical-align: -6px;\"\/>\u00a0[5]<\/p>\n<p align=\"JUSTIFY\">where L is the inductance, R is the resistance, C is the capacitance, I is the current, <img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-11ccf66e35a6ad701edfee2663acb536_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#92;&#111;&#109;&#101;&#103;&#97;\" title=\"Rendered by QuickLaTeX.com\" height=\"8\" width=\"11\" style=\"vertical-align: 0px;\"\/>\u00a0is the\u00a0resonant\u00a0angular frequency, <img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-2cf56434dc031c28acb706e4308b1df8_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#32;&#69;&#95;&#123;&#48;&#125;\" title=\"Rendered by QuickLaTeX.com\" height=\"15\" width=\"20\" style=\"vertical-align: -3px;\"\/> is\u00a0an initial value dependent on the voltage source,\u00a0and t is the time. The inductance\u00a0(L), resistance\u00a0(R),\u00a0and\u00a0capacitance\u00a0(C)are all determined by their respective circuit components,\u00a0which are easily changed out\u00a0for different ones.\u00a0<img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-11ccf66e35a6ad701edfee2663acb536_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#92;&#111;&#109;&#101;&#103;&#97;\" title=\"Rendered by QuickLaTeX.com\" height=\"8\" width=\"11\" style=\"vertical-align: 0px;\"\/> depends\u00a0on the inductance and capacitance, and its value is given by<\/p>\n<p align=\"JUSTIFY\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-003c7233f95d3889807c2089378f7804_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#92;&#111;&#109;&#101;&#103;&#97;&#32;&#61;&#92;&#102;&#114;&#97;&#99;&#123;&#49;&#125;&#123;&#92;&#115;&#113;&#114;&#116;&#123;&#76;&#67;&#125;&#125;\" title=\"Rendered by QuickLaTeX.com\" height=\"27\" width=\"70\" style=\"vertical-align: -11px;\"\/> [6]<\/p>\n<p align=\"JUSTIFY\">What we have left in terms of variables are current (I) and time (t). For our purposes, time is the independent variable, which leaves the current in the circuit at any time to be the dependent variable. Now that we have established this, we can solve our differential equation [5] for current as a function of time (<img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-3fdb038c5615c4490dae465213f7e0dd_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#73;&#40;&#116;&#41;\" title=\"Rendered by QuickLaTeX.com\" height=\"18\" width=\"28\" style=\"vertical-align: -4px;\"\/>).<\/p>\n<p align=\"JUSTIFY\"><strong>Solving the Differential Equation<\/strong><\/p>\n<p>In order to solve [5], I will use the equation solving powers of Mathematica. In particular, the DSolve function is used.\u00a0\u00a0When plugged into Mathematica, the input line looks like:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-4fd6cd2ab6f0ba47fdfaab15b71e9d9e_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#92;&#116;&#101;&#120;&#116;&#123;&#68;&#83;&#111;&#108;&#118;&#101;&#125;&#92;&#108;&#101;&#102;&#116;&#91;&#92;&#102;&#114;&#97;&#99;&#123;&#92;&#116;&#101;&#120;&#116;&#123;&#67;&#117;&#114;&#114;&#101;&#110;&#116;&#125;&#40;&#116;&#41;&#125;&#123;&#92;&#116;&#101;&#120;&#116;&#123;&#67;&#97;&#112;&#97;&#99;&#105;&#116;&#97;&#110;&#99;&#101;&#125;&#125;&#43;&#92;&#116;&#101;&#120;&#116;&#123;&#73;&#110;&#100;&#117;&#99;&#116;&#97;&#110;&#99;&#101;&#125;&#32;&#92;&#116;&#101;&#120;&#116;&#123;&#67;&#117;&#114;&#114;&#101;&#110;&#116;&#125;&#39;&#39;&#40;&#116;&#41;&#43;&#92;&#116;&#101;&#120;&#116;&#123;&#82;&#101;&#115;&#105;&#115;&#116;&#97;&#110;&#99;&#101;&#125;&#32;&#92;&#116;&#101;&#120;&#116;&#123;&#67;&#117;&#114;&#114;&#101;&#110;&#116;&#125;&#39;&#40;&#116;&#41;&#61;&#69;&#95;&#123;&#48;&#125;&#32;&#92;&#111;&#109;&#101;&#103;&#97;&#32;&#92;&#99;&#111;&#115;&#32;&#40;&#92;&#111;&#109;&#101;&#103;&#97;&#32;&#116;&#41;&#44;&#92;&#116;&#101;&#120;&#116;&#123;&#67;&#117;&#114;&#114;&#101;&#110;&#116;&#125;&#40;&#116;&#41;&#44;&#116;&#92;&#114;&#105;&#103;&#104;&#116;&#93;\" title=\"Rendered by QuickLaTeX.com\" height=\"32\" width=\"749\" style=\"vertical-align: -11px;\"\/> [7]<\/p>\n<p>We must set values for resistance, inductance, capacitance,\u00a0and\u00a0<img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-2cf56434dc031c28acb706e4308b1df8_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#32;&#69;&#95;&#123;&#48;&#125;\" title=\"Rendered by QuickLaTeX.com\" height=\"15\" width=\"20\" style=\"vertical-align: -3px;\"\/>. For this solution, I chose the same values as the first circuit :<\/p>\n<ul>\n<li>Inductance= 1 Henry<\/li>\n<li>Resistance= 10 ohms<\/li>\n<li>Capacitance= .1 farads<\/li>\n<li><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-2cf56434dc031c28acb706e4308b1df8_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#32;&#69;&#95;&#123;&#48;&#125;\" title=\"Rendered by QuickLaTeX.com\" height=\"15\" width=\"20\" style=\"vertical-align: -3px;\"\/>= 1 Volt<\/li>\n<\/ul>\n<p>When the command is run with the above values set, Mathematica generates the general solution, which is:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-e92f0ada61a2748426a72df614104599_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#73;&#40;&#116;&#41;&#61;&#99;&#95;&#49;&#32;&#101;&#94;&#123;&#45;&#56;&#46;&#56;&#55;&#50;&#57;&#56;&#32;&#116;&#125;&#43;&#99;&#95;&#50;&#32;&#101;&#94;&#123;&#45;&#49;&#46;&#49;&#50;&#55;&#48;&#50;&#32;&#116;&#125;&#43;&#48;&#46;&#48;&#52;&#48;&#56;&#50;&#52;&#40;&#45;&#67;&#111;&#115;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#43;&#101;&#94;&#123;&#54;&#46;&#54;&#54;&#49;&#42;&#49;&#48;&#94;&#45;&#49;&#54;&#116;&#125;&#67;&#111;&#115;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#45;&#48;&#46;&#51;&#53;&#54;&#52;&#83;&#105;&#110;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#43;&#50;&#46;&#56;&#48;&#54;&#101;&#94;&#123;&#54;&#46;&#54;&#54;&#49;&#42;&#49;&#48;&#94;&#45;&#49;&#54;&#116;&#125;&#83;&#105;&#110;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#41;\" title=\"Rendered by QuickLaTeX.com\" height=\"41\" width=\"589\" style=\"vertical-align: -4px;\"\/> [8]<\/p>\n<p>First, lets make sure that this solution is reasonable. The known form of the solution to the equation for such a circuit is given by:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-f1e55e04d02c3795b6e87f169c4164d1_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#73;&#40;&#116;&#41;&#61;&#99;&#95;&#49;&#32;&#101;&#94;&#123;&#114;&#95;&#123;&#49;&#125;&#116;&#125;&#43;&#99;&#95;&#50;&#32;&#101;&#94;&#123;&#114;&#95;&#123;&#50;&#125;&#116;&#125;&#43;&#65;&#115;&#105;&#110;&#40;&#92;&#111;&#109;&#101;&#103;&#97;&#32;&#116;&#45;&#92;&#118;&#97;&#114;&#112;&#104;&#105;&#32;&#41;\" title=\"Rendered by QuickLaTeX.com\" height=\"19\" width=\"284\" style=\"vertical-align: -4px;\"\/> [9]<\/p>\n<p>The form of our equation [8] seems to match up pretty nicely with this known solution [9]. Now, lets examine [8] in more detail.<\/p>\n<p>Our solution [8] looks a little daunting, but we can break it down into its components. The first two terms have constants that depend on the initial current running through the system, i.e. when <img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-1eecbd67175c41e11d7a3bc6a02ffbf8_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#73;&#40;&#116;&#41;&#61;&#73;&#40;&#48;&#41;\" title=\"Rendered by QuickLaTeX.com\" height=\"18\" width=\"84\" style=\"vertical-align: -4px;\"\/>. However, we can approximate the equation by only looking at the last term, which is a superposition of sin and cosine functions along with a couple exponential terms thrown in. The reason that we can apply this approximation is that the first two terms &#8220;damp out&#8221; quickly, and thus have exponentially less of an effect as time progresses. The reader may still be curious to see how much of an effect these terms have on the current flowing through the system: to that I say don&#8217;t worry! After my initial approximation, I will solve for these constants and compare the two cases.<\/p>\n<p><b>Approximate<\/b><strong>\u00a0Solution and Graphs<\/strong><\/p>\n<p>For now, we will deal with the approximate solution, which is:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-ee1e7c19febbeffa364281e06d9c85a5_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#73;&#40;&#116;&#41;&#61;&#48;&#46;&#48;&#52;&#48;&#56;&#50;&#52;&#40;&#45;&#67;&#111;&#115;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#43;&#101;&#94;&#123;&#54;&#46;&#54;&#54;&#49;&#42;&#49;&#48;&#94;&#45;&#49;&#54;&#116;&#125;&#67;&#111;&#115;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#45;&#48;&#46;&#51;&#53;&#54;&#52;&#83;&#105;&#110;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#43;&#50;&#46;&#56;&#48;&#54;&#101;&#94;&#123;&#54;&#46;&#54;&#54;&#49;&#42;&#49;&#48;&#94;&#45;&#49;&#54;&#116;&#125;&#83;&#105;&#110;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#41;\" title=\"Rendered by QuickLaTeX.com\" height=\"41\" width=\"582\" style=\"vertical-align: -4px;\"\/> [10]<\/p>\n<p>Now that we have this equation, we can use Mathematica to give us an idea of what the changing current looks like over a longer period of time. Figure 2 shows the current in the circuit as a function of time (<img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-3fdb038c5615c4490dae465213f7e0dd_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#73;&#40;&#116;&#41;\" title=\"Rendered by QuickLaTeX.com\" height=\"18\" width=\"28\" style=\"vertical-align: -4px;\"\/>) plotted over a 30-second interval.<\/p>\n<p><a href=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/30secondplot.jpg\"><img loading=\"lazy\" decoding=\"async\" alt=\"30secondplot\" src=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/30secondplot.jpg\" width=\"542\" height=\"363\" \/><\/a>(Figure 4)<\/p>\n<p>This looks like something we would expect! A sinusoidal wave that describes the changing current over time, just like our equations dictate. But, this is a short time-scale. What happens when we look at the changing current over a longer period of time? In Figure 3, I show the changing current over a longer time period (3000 seconds\/50 minutes). In Figure 4, I repeat the process, but with the timescale being 5 years.<\/p>\n<p><a href=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/50minuteplot.jpg\"><img loading=\"lazy\" decoding=\"async\" alt=\"50minuteplot\" src=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/50minuteplot.jpg\" width=\"526\" height=\"337\" \/><\/a>(Figure 5)<\/p>\n<p><a href=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/5yearplot1.jpg\"><img loading=\"lazy\" decoding=\"async\" alt=\"5yearplot\" src=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/5yearplot1.jpg\" width=\"532\" height=\"318\" \/><\/a>(Figure 6)<\/p>\n<p>In both Figures 5 and 6, we see some very interesting long-term variations in the current. Why does this occur?<\/p>\n<p><b>Computational Limits<\/b><\/p>\n<p>It\u00a0is clear that Mathematica is displaying some odd variations that do not seem to have a set pattern. To investigate this and determine its cause, lets look closer at a range of time where the current seems to be irregular in the previous two figures. Specifically, lets look at the time scale between 300 and 400 seconds because it seems to have very strange behavior.<\/p>\n<p><a href=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata1.jpg\"><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-3616\" alt=\"RLCCircuitsworking_finaldata1\" src=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata1.jpg\" width=\"396\" height=\"250\" srcset=\"https:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata1.jpg 396w, https:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata1-300x189.jpg 300w\" sizes=\"auto, (max-width: 396px) 100vw, 396px\" \/><\/a>[Figure 7]<\/p>\n<p>Looking at the graph above, we see that it looks completely normal! There are no strange variations and the current continues the predicted sinusoidal trend. This leads me to the conclusion that Mathematica has limitations that caused the irregular variations that we saw before, and that they do not actually exist in the data\/circuit.<\/p>\n<p><strong>Completing the Full Solution via Graphical Solving<\/strong><\/p>\n<p>The next and final step in completing our model is accounting for the two initial terms that we left out. Lets take another look at the full equation.<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-e92f0ada61a2748426a72df614104599_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#73;&#40;&#116;&#41;&#61;&#99;&#95;&#49;&#32;&#101;&#94;&#123;&#45;&#56;&#46;&#56;&#55;&#50;&#57;&#56;&#32;&#116;&#125;&#43;&#99;&#95;&#50;&#32;&#101;&#94;&#123;&#45;&#49;&#46;&#49;&#50;&#55;&#48;&#50;&#32;&#116;&#125;&#43;&#48;&#46;&#48;&#52;&#48;&#56;&#50;&#52;&#40;&#45;&#67;&#111;&#115;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#43;&#101;&#94;&#123;&#54;&#46;&#54;&#54;&#49;&#42;&#49;&#48;&#94;&#45;&#49;&#54;&#116;&#125;&#67;&#111;&#115;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#45;&#48;&#46;&#51;&#53;&#54;&#52;&#83;&#105;&#110;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#43;&#50;&#46;&#56;&#48;&#54;&#101;&#94;&#123;&#54;&#46;&#54;&#54;&#49;&#42;&#49;&#48;&#94;&#45;&#49;&#54;&#116;&#125;&#83;&#105;&#110;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#41;\" title=\"Rendered by QuickLaTeX.com\" height=\"41\" width=\"589\" style=\"vertical-align: -4px;\"\/> [11]<\/p>\n<p>We can see that they are exponential functions with a separate constant for each term. Solving for these constants quantitatively can prove to be difficult, so I will attempt to solve for them via a guess-and-check method using plots in Mathematica. I need to match the initial condition so that Ohm&#8217;s law is satisfied.\u00a0So,\u00a0<img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-0d913679a42337de0fe847ab1e3c2d62_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#73;&#61;&#92;&#102;&#114;&#97;&#99;&#123;&#86;&#125;&#123;&#82;&#125;\" title=\"Rendered by QuickLaTeX.com\" height=\"22\" width=\"47\" style=\"vertical-align: -6px;\"\/> must hold in the initial state. This means that the initial current must be 0.1 amps. Lets try a few values for\u00a0<img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-86d57aecc373976151de585b5b2bec5a_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#67;&#91;&#49;&#93;\" title=\"Rendered by QuickLaTeX.com\" height=\"18\" width=\"31\" style=\"vertical-align: -5px;\"\/> and\u00a0<img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-12b2dbf736e307ece734ae42b09665aa_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#67;&#91;&#50;&#93;\" title=\"Rendered by QuickLaTeX.com\" height=\"18\" width=\"31\" style=\"vertical-align: -5px;\"\/> and figure out graphically what they should be.<\/p>\n<p>First, lets try values of 0.5 for each, resulting in the following equation:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-d760534d4dd55ff8bd60f98ea22a7b2b_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#73;&#40;&#116;&#41;&#61;&#48;&#46;&#53;&#101;&#94;&#123;&#45;&#56;&#46;&#56;&#55;&#50;&#57;&#56;&#32;&#116;&#125;&#43;&#48;&#46;&#53;&#101;&#94;&#123;&#45;&#49;&#46;&#49;&#50;&#55;&#48;&#50;&#32;&#116;&#125;&#43;&#48;&#46;&#48;&#52;&#48;&#56;&#50;&#52;&#40;&#45;&#67;&#111;&#115;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#43;&#101;&#94;&#123;&#54;&#46;&#54;&#54;&#49;&#42;&#49;&#48;&#94;&#45;&#49;&#54;&#116;&#125;&#67;&#111;&#115;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#45;&#48;&#46;&#51;&#53;&#54;&#52;&#83;&#105;&#110;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#43;&#50;&#46;&#56;&#48;&#54;&#101;&#94;&#123;&#54;&#46;&#54;&#54;&#49;&#42;&#49;&#48;&#94;&#45;&#49;&#54;&#116;&#125;&#83;&#105;&#110;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#41;\" title=\"Rendered by QuickLaTeX.com\" height=\"41\" width=\"589\" style=\"vertical-align: -4px;\"\/> [12]<\/p>\n<p>The graph resulting from this equation is:<\/p>\n<p><a href=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata6.jpg\"><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-3623\" alt=\"RLCCircuitsworking_finaldata6\" src=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata6.jpg\" width=\"396\" height=\"255\" srcset=\"https:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata6.jpg 396w, https:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata6-300x193.jpg 300w\" sizes=\"auto, (max-width: 396px) 100vw, 396px\" \/><\/a>[Figure 8]<\/p>\n<p>We can see that the initial current that this combination of constants is too high! It gives 1.0 amps, when we are looking for 0.1 amps.<\/p>\n<p>Lets try slightly higher values for each, as well as having different values for the two. The equation that results is:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-225e385a4b28a26cc571b90eb307f7b7_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#73;&#40;&#116;&#41;&#61;&#50;&#101;&#94;&#123;&#45;&#56;&#46;&#56;&#55;&#50;&#57;&#56;&#32;&#116;&#125;&#43;&#51;&#101;&#94;&#123;&#45;&#49;&#46;&#49;&#50;&#55;&#48;&#50;&#32;&#116;&#125;&#43;&#48;&#46;&#48;&#52;&#48;&#56;&#50;&#52;&#40;&#45;&#67;&#111;&#115;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#43;&#101;&#94;&#123;&#54;&#46;&#54;&#54;&#49;&#42;&#49;&#48;&#94;&#45;&#49;&#54;&#116;&#125;&#67;&#111;&#115;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#45;&#48;&#46;&#51;&#53;&#54;&#52;&#83;&#105;&#110;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#43;&#50;&#46;&#56;&#48;&#54;&#101;&#94;&#123;&#54;&#46;&#54;&#54;&#49;&#42;&#49;&#48;&#94;&#45;&#49;&#54;&#116;&#125;&#83;&#105;&#110;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#41;\" title=\"Rendered by QuickLaTeX.com\" height=\"41\" width=\"589\" style=\"vertical-align: -4px;\"\/> [13]<\/p>\n<p>Which produces the following graph:<\/p>\n<p><a href=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata7.jpg\"><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-3625\" alt=\"RLCCircuitsworking_finaldata7\" src=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata7.jpg\" width=\"396\" height=\"255\" srcset=\"https:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata7.jpg 396w, https:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata7-300x193.jpg 300w\" sizes=\"auto, (max-width: 396px) 100vw, 396px\" \/><\/a>[Figure 9]<\/p>\n<p>Looks like we have gone the wrong way with our guess! Lets try lowering both the constants to arrive at the correct graphical solution (in reality, this took me many more iterations of guessing and checking):<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/pages.vassar.edu\/magnes\/wp-content\/ql-cache\/quicklatex.com-0d29e1920900a12c78e1f9855406b3f1_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#73;&#40;&#116;&#41;&#61;&#48;&#46;&#48;&#53;&#101;&#94;&#123;&#45;&#56;&#46;&#56;&#55;&#50;&#57;&#56;&#32;&#116;&#125;&#43;&#48;&#46;&#48;&#53;&#101;&#94;&#123;&#45;&#49;&#46;&#49;&#50;&#55;&#48;&#50;&#32;&#116;&#125;&#43;&#48;&#46;&#48;&#52;&#48;&#56;&#50;&#52;&#40;&#45;&#67;&#111;&#115;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#43;&#101;&#94;&#123;&#54;&#46;&#54;&#54;&#49;&#42;&#49;&#48;&#94;&#45;&#49;&#54;&#116;&#125;&#67;&#111;&#115;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#45;&#48;&#46;&#51;&#53;&#54;&#52;&#83;&#105;&#110;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#43;&#50;&#46;&#56;&#48;&#54;&#101;&#94;&#123;&#54;&#46;&#54;&#54;&#49;&#42;&#49;&#48;&#94;&#45;&#49;&#54;&#116;&#125;&#83;&#105;&#110;&#40;&#51;&#46;&#49;&#54;&#50;&#116;&#41;&#41;\" title=\"Rendered by QuickLaTeX.com\" height=\"41\" width=\"589\" style=\"vertical-align: -4px;\"\/> [14]<\/p>\n<p>Which results in the following graph:<\/p>\n<p><a href=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata8.jpg\"><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-3626\" alt=\"RLCCircuitsworking_finaldata8\" src=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata8.jpg\" width=\"396\" height=\"255\" srcset=\"https:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata8.jpg 396w, https:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata8-300x193.jpg 300w\" sizes=\"auto, (max-width: 396px) 100vw, 396px\" \/><\/a>[Figure 10]<\/p>\n<p><strong>The Graphical Solution<\/strong><\/p>\n<p>Now that we have arrived at the correct graphical solution, lets change the scaling of the graph to better match that of [Figure 2], which was the initial approximation that did not take into account the two damping exponential terms.<\/p>\n<p><a href=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata9.jpg\"><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-3630\" alt=\"RLCCircuitsworking_finaldata9\" src=\"http:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata9.jpg\" width=\"396\" height=\"238\" srcset=\"https:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata9.jpg 396w, https:\/\/pages.vassar.edu\/magnes\/files\/2014\/04\/RLCCircuitsworking_finaldata9-300x180.jpg 300w\" sizes=\"auto, (max-width: 396px) 100vw, 396px\" \/><\/a>[Figure 11]<\/p>\n<p>&nbsp;<\/p>\n<p>My conclusion from this final data will follow in my next blog post.<\/p>\n<p>&nbsp;<\/p>\n<p><strong>Resources<\/strong><\/p>\n<p>1. Circuit Diagram Maker-\u00a0http:\/\/www.circuit-diagram.org\/<\/p>\n<p>2.\u00a0The RLC Circuit- University of British Columbia- \u00a0<a href=\"http:\/\/www.math.ubc.ca\/~feldman\/m121\/RLC.pdf\">http:\/\/www.math.ubc.ca\/~feldman\/m121\/RLC.pdf<\/a><\/p>\n<p>3.\u00a0Mathematica Cookbook by Sal Mangano<\/p>\n<p>4. Electronic Circuit Analysis for Scientists by James A. McCray and Thomas A. Cahill<\/p>\n<p>5. Dynamical Systems with Applications using Mathematica by Stephen Lynch<\/p>\n<p>6.\u00a0Class Notes- Mathematics 228 (Methods of Applied Mathematics) taught by Matthew Miller<\/p>\n<p>&nbsp;<\/p>\n<p><strong>Mathematica Notebook<\/strong><\/p>\n<p><a title=\"Mathematica Notebook\" href=\"https:\/\/drive.google.com\/file\/d\/0B3mtB6CQNnpjbXROcHl0VURzYlE\/edit?usp=sharing\" target=\"_blank\">https:\/\/drive.google.com\/file\/d\/0B3mtB6CQNnpjbXROcHl0VURzYlE\/edit?usp=sharing<\/a><\/p>\n","protected":false},"excerpt":{"rendered":"<p>Overview After talking with Professor Magnes after my preliminary results, I decided to first look at an RLC circuit with no voltage source, which is much simpler to model. Doing this allowed me to better understand the mechanism behind the circuitry and gain more experience with Mathematica as well. After doing so, I returned to [&hellip;]<\/p>\n","protected":false},"author":1826,"featured_media":0,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[4101,54058,188,54190],"tags":[],"class_list":["post-3608","post","type-post","status-publish","format-standard","hentry","category-advanced-em","category-gagandeep","category-mathematica","category-spring-2014"],"_links":{"self":[{"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/posts\/3608","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\/1826"}],"replies":[{"embeddable":true,"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/comments?post=3608"}],"version-history":[{"count":8,"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/posts\/3608\/revisions"}],"predecessor-version":[{"id":3718,"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/posts\/3608\/revisions\/3718"}],"wp:attachment":[{"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/media?parent=3608"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/categories?post=3608"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/tags?post=3608"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}