{"id":865,"date":"2021-09-03T17:25:47","date_gmt":"2021-09-03T17:25:47","guid":{"rendered":"https:\/\/courses.lumenlearning.com\/calculus3\/?post_type=chapter&#038;p=865"},"modified":"2022-10-29T01:50:24","modified_gmt":"2022-10-29T01:50:24","slug":"partial-differential-equations","status":"publish","type":"chapter","link":"https:\/\/courses.lumenlearning.com\/calculus3\/chapter\/partial-differential-equations\/","title":{"raw":"Partial Differential Equations","rendered":"Partial Differential Equations"},"content":{"raw":"<div class=\"textbox learning-objectives\">\r\n<h3>Learning Outcomes<\/h3>\r\n<ul>\r\n \t<li>Explain the meaning of a partial differential equation and give an example.<\/li>\r\n<\/ul>\r\n<\/div>\r\nIn <a href=\"https:\/\/courses.lumenlearning.com\/calculus2\/chapter\/introduction-4-2\/\" target=\"_blank\" rel=\"noopener\" data-book-uuid=\"1d39a348-071f-4537-85b6-c98912458c3c\" data-page-slug=\"4-introduction\">Introduction to Differential Equations<\/a>, we studied differential equations in which the unknown function had one independent variable. A <strong><span id=\"b1ffd73f-b9d7-41a1-8546-c25c366b2e90_term170\" data-type=\"term\">partial differential equation<\/span><\/strong> is an equation that involves an unknown function of more than one independent variable and one or more of its partial derivatives. Examples of partial differential equations are\r\n<p style=\"text-align: center;\">[latex]\\large{u_t=c^{2}\\,(u_{xx}+u_{yy})}[\/latex]<\/p>\r\n(heat equation in two dimensions)\r\n<p style=\"text-align: center;\">[latex]\\large{u_{tt}=c^{2}\\,(u_{xx}+u_{yy})}[\/latex]<\/p>\r\n(wave equation in two dimensions)\r\n<p style=\"text-align: center;\">[latex]\\large{u_{xx}+u_{yy}=0}[\/latex]<\/p>\r\n(Laplace's equation in two dimensions)\r\n\r\nIn the first two equations, the unknown function [latex]u[\/latex] has three independent variables[latex]-t,\\ x[\/latex], and [latex]y-[\/latex]and [latex]c[\/latex] is an arbitrary constant. The independent variables [latex]x[\/latex] and [latex]y[\/latex] are considered to be spatial variables, and the variable [latex]t[\/latex] represents time. In Laplace's equation, the unknown function [latex]u[\/latex] has two independent variables [latex]x[\/latex] and [latex]y[\/latex].\r\n<div class=\"textbox exercises\">\r\n<h3>Example: A Solution to the Wave Equation<\/h3>\r\nVerify that\r\n<p style=\"text-align: center;\">[latex]\\large{u\\,(x,\\ y,\\ t)=5\\sin{(3\\pi x)}\\sin{(4\\pi y)}\\cos{(10\\pi t)}}[\/latex]<\/p>\r\nis a solution to the wave equation\r\n<p style=\"text-align: center;\">[latex]\\large{u_{tt}=4\\,(u_{xx}+u_{yy})}.[\/latex]<\/p>\r\n[reveal-answer q=\"fs-id1577773930324\"]Show Solution[\/reveal-answer]\r\n[hidden-answer a=\"fs-id1577773930324\"]\r\n<div style=\"text-align: left;\">First, we calculate [latex]u_{tt},\\ u_{xx}[\/latex], and [latex]u_{yy}[\/latex]:<\/div>\r\n<div style=\"text-align: center;\">[latex]\\begin{array}{ccc}\\hfill {u_{tt}} &amp; =\\hfill &amp; {\\frac{\\partial}{\\partial t}\\left[\\frac{\\partial u}{\\partial t}\\right]} \\hfill \\\\ \\hfill &amp; =\\hfill &amp; {\\frac{\\partial}{\\partial t}[5\\sin{(3\\pi x)}\\sin{(4\\pi y)}(-10\\pi\\sin{(10\\pi t)})]}\\hfill \\\\ \\hfill &amp; =\\hfill &amp; {\\frac{\\partial}{\\partial t}[-50\\pi\\sin{(3\\pi x)}\\sin{(4\\pi y)}\\sin{(10\\pi t)}]}\\hfill \\\\ \\hfill &amp; =\\hfill &amp; {-500\\pi^{2}\\sin{(3\\pi x)}\\sin{(4\\pi y)}\\cos{(10\\pi t)}}\\hfill \\\\ \\hfill {u_{xx}} &amp; =\\hfill &amp; {\\frac{\\partial}{\\partial x}\\left[\\frac{\\partial u}{\\partial x}\\right]} \\hfill \\\\ \\hfill &amp; =\\hfill &amp; {\\frac{\\partial}{\\partial x}[15\\pi\\cos{(3\\pi x)}\\sin{(4\\pi y)}\\cos{(10\\pi t)}]}\\hfill \\\\ \\hfill &amp; =\\hfill &amp; {-45\\pi^{2}\\sin{(3\\pi x)}\\sin{(4\\pi y)}\\cos{(10\\pi t)}}\\hfill \\\\ \\hfill {u_{yy}} &amp; =\\hfill &amp; {\\frac{\\partial}{\\partial y}\\left[\\frac{\\partial u}{\\partial y}\\right]}\\hfill \\\\ \\hfill &amp; =\\hfill &amp; {\\frac{\\partial}{\\partial y}[5\\sin{(3\\pi x)}(4\\pi\\cos{(4\\pi y)})\\cos{(10\\pi t)}]}\\hfill \\\\ \\hfill &amp; =\\hfill &amp; {\\frac{\\partial}{\\partial y}[20\\pi\\sin{(3\\pi x)}\\cos{(4\\pi y)}\\cos{(10\\pi t)}]}\\hfill \\\\ \\hfill &amp; =\\hfill &amp; {-80\\pi^{2}\\sin{(3\\pi x)}\\sin{(4\\pi y)}\\cos{(10\\pi t)}.}\\hfill \\\\ \\hfill \\end{array}[\/latex]<\/div>\r\n<div>Next, we substitute each of these into the right-hand side of our wave equation and simplify:<\/div>\r\n<div><\/div>\r\n<div>\r\n\r\n[latex] \\text{ }\\\\ \\begin{alignat}{2} 4(u_{xx}+u_{yy}) &amp;= 4(-45\\pi^2\\sin{(3\\pi x)}\\sin{(4\\pi y)}\\cos{(10\\pi t)} + -80\\pi^2\\sin{(3\\pi x)}\\sin{(4\\pi y)}\\cos{(10\\pi t)}\\\\\r\n\r\n\\hspace{2cm}&amp;=4(-125\\pi^2\\sin{(3\\pi x)}\\sin{(4\\pi y)}\\cos{(10\\pi t)} &amp;\\quad\\\\\r\n\r\n&amp;=-500\\pi^2\\sin{(3\\pi x)}\\sin{(4\\pi y)}\\cos{(10\\pi t)} \\\\\r\n\r\n&amp;=u_{tt}.\\\\\r\n\r\n\\end{alignat}[\/latex]\r\n\r\nThis verifies the solution.\r\n\r\n<\/div>\r\n[\/hidden-answer]\r\n\r\n<\/div>\r\n<div class=\"textbox key-takeaways\">\r\n<h3>try it<\/h3>\r\nVerify that [latex]u\\,(x,\\ y,\\ t)=2\\,\\sin\\,(\\frac{x}{3})\\,\\sin\\,(\\frac{y}{4})\\,e^{-25t\/16}[\/latex] is a solution to the heat equation\r\n<p style=\"text-align: center;\">[latex]u_{t}=9(u_{xx}+u_{yy})[\/latex].<\/p>\r\n[reveal-answer q=\"fs-id1367573930124\"]Show Solution[\/reveal-answer]\r\n[hidden-answer a=\"fs-id1367573930124\"]\r\n<div style=\"text-align: center;\"><\/div>\r\n[latex]u\\,(x,\\ y,\\ t)=2\\,\\sin\\,(\\frac{x}{3})\\,\\sin\\,(\\frac{y}{4})\\,e^{-25t\/16}[\/latex] is a solution to the heat equation. See worked example video below to verify the solution.\r\n[\/hidden-answer]\r\n\r\n<\/div>\r\n\r\n[caption]Watch the following video to see the worked solution to the above Try It[\/caption]\r\n\r\n<center><iframe src=\"\/\/plugin.3playmedia.com\/show?mf=8186156&amp;p3sdk_version=1.10.1&amp;p=20361&amp;pt=375&amp;video_id=pMoH9Key-B8&amp;video_target=tpm-plugin-k0fmirvu-pMoH9Key-B8\" width=\"800px\" height=\"450px\" frameborder=\"0\" marginwidth=\"0px\" marginheight=\"0px\"><\/iframe><\/center><center>You can view the <a href=\"https:\/\/course-building.s3.us-west-2.amazonaws.com\/Calculus+3\/Calc+3+transcripts\/CP4.18_transcript.html\">transcript for \u201cCP 4.18\u201d here (opens in new window).<\/a><\/center>Since the solution to the two-dimensional heat equation is a function of three variables, it is not easy to create a visual representation of the solution. We can graph the solution for fixed values of\u00a0[latex]t[\/latex], which amounts to snapshots of the heat distributions at fixed times. These snapshots show how the heat is distributed over a two-dimensional surface as time progresses. The graph of the preceding solution at time [latex]t=0[\/latex] appears in the following figure. As time progresses, the extremes level out, approaching zero as [latex]t[\/latex] approaches infinity.\r\n\r\n[caption id=\"attachment_985\" align=\"aligncenter\" width=\"456\"]<img class=\"size-full wp-image-985\" src=\"https:\/\/s3-us-west-2.amazonaws.com\/courses-images\/wp-content\/uploads\/sites\/5667\/2021\/09\/28230317\/4-3-3.jpeg\" alt=\"A complicated curve in xyz space with many sinusoidally alternating local maxima and minima.\" width=\"456\" height=\"379\" \/> Figure 1.[\/caption]\r\n\r\nIf we consider the heat equation in one dimension, then it is possible to graph the solution over time. The heat equation in one dimension becomes\r\n<p style=\"text-align: center;\">[latex]u_{t}=c^{2}u_{xx}[\/latex],<\/p>\r\nwhere [latex]c^{2}[\/latex] represents the thermal diffusivity of the material in question. A solution of this differential equation can be written in the form\r\n<p style=\"text-align: center;\">[latex]u_{m}\\,(x,\\ t)=e^{-\\pi^{2}m^{2}c^{2}t}\\sin{(m\\pi{x})}[\/latex]<\/p>\r\nwhere [latex]m[\/latex] is any positive integer. A graph of this solution using [latex]m=1[\/latex] appears in\u00a0Figure 4, where the initial temperature distribution over a wire of length [latex]1[\/latex] is given by [latex]u\\,(x,\\ 0)=\\sin{\\pi{x}}[\/latex]. Notice that as time progresses, the wire cools off. This is seen because, from left to right, the highest temperature (which occurs in the middle of the wire) decreases and changes color from red to blue.\r\n\r\n[caption id=\"attachment_986\" align=\"aligncenter\" width=\"707\"]<img class=\"size-full wp-image-986\" src=\"https:\/\/s3-us-west-2.amazonaws.com\/courses-images\/wp-content\/uploads\/sites\/5667\/2021\/09\/28230533\/4-3-4.jpeg\" alt=\"A curve in xtu space with a local maximum at (0.5, 0, 12). From this maximum, the values decrease with increasing t and for any value of x.\" width=\"707\" height=\"648\" \/> Figure 2.\u00a0Graph of a solution of the heat equation in one dimension over time.[\/caption]\r\n\r\n<div class=\"textbox tryit\">\r\n<h3 style=\"text-align: center;\" data-type=\"title\">Lord Kelvin and the Age of Earth<\/h3>\r\n\r\n<hr \/>\r\n\r\n[caption id=\"attachment_987\" align=\"aligncenter\" width=\"731\"]<img class=\"size-full wp-image-987\" src=\"https:\/\/s3-us-west-2.amazonaws.com\/courses-images\/wp-content\/uploads\/sites\/5667\/2021\/09\/28230657\/4-3-5.jpeg\" alt=\"This figure consists of two figures marked a and b. Figure a show Lord Kelvin, dressed well and with a beard. Figure b shows an image of the planet Earth taken from space.\" width=\"731\" height=\"429\" \/> Figure 3.\u00a0(a) William Thomson (Lord Kelvin), 1824-1907, was a British physicist and electrical engineer; (b) Kelvin used the heat diffusion equation to estimate the age of Earth (credit: modification of work by NASA).[\/caption]\r\n\r\nDuring the late 1800s, the scientists of the new field of geology were coming to the conclusion that Earth must be \u201cmillions and millions\u201d of years old. At about the same time, Charles Darwin had published his treatise on evolution. Darwin\u2019s view was that evolution needed many millions of years to take place, and he made a bold claim that the Weald chalk fields, where important fossils were found, were the result of [latex]300[\/latex] million years of erosion.\r\n\r\nAt that time, eminent physicist <span id=\"term174\" class=\"no-emphasis\" data-type=\"term\">William Thomson (Lord Kelvin)<\/span> used an important partial differential equation, known as the <em data-effect=\"italics\">heat diffusion equation<\/em>, to estimate the age of Earth by determining how long it would take Earth to cool from molten rock to what we had at that time. His conclusion was a range of [latex]20[\/latex] to [latex]400[\/latex] million years, but most likely about [latex]50[\/latex] million years. For many decades, the proclamations of this irrefutable icon of science did not sit well with geologists or with Darwin.\r\n<div class=\"textbox\">\r\n<h3 style=\"text-align: center;\" data-type=\"title\">Interactive<\/h3>\r\n\r\n<hr \/>\r\n\r\nRead Kelvin\u2019s <a href=\"http:\/\/www.openstax.org\/l\/20_KelEarthAge\" target=\"_blank\" rel=\"noopener nofollow\">paper<\/a> on estimating the age of the Earth.\r\n\r\n<\/div>\r\n<p id=\"fs-id1167794213455\" class=\" \">Kelvin made reasonable assumptions based on what was known in his time, but he also made several assumptions that turned out to be wrong. One incorrect assumption was that Earth is solid and that the cooling was therefore via conduction only, hence justifying the use of the diffusion equation. But the most serious error was a forgivable one\u2014omission of the fact that Earth contains radioactive elements that continually supply heat beneath Earth\u2019s mantle. The discovery of radioactivity came near the end of Kelvin\u2019s life and he acknowledged that his calculation would have to be modified.<\/p>\r\n<p id=\"fs-id1167794213464\" class=\" \">Kelvin used the simple one-dimensional model applied only to Earth\u2019s outer shell, and derived the age from graphs and the roughly known temperature gradient near Earth\u2019s surface. Let\u2019s take a look at a more appropriate version of the diffusion equation in radial coordinates, which has the form<\/p>\r\n\r\n<div style=\"text-align: center;\">[latex]\\large{\\frac{\\partial{T}}{\\partial t}=K\\left[\\frac{\\partial^{2}T}{\\partial^{2}r}+\\frac{2}{r}\\frac{\\partial T}{\\partial r}\\right]}[\/latex].<\/div>\r\n&nbsp;\r\n\r\nHere, [latex]T\\,(r,\\ t)[\/latex] is temperature as a function of [latex]r[\/latex] (measured from the center of Earth) and time [latex]t[\/latex]. [latex]K[\/latex] is the heat conductivity\u00a0conductivity\u2014for molten rock, in this case. The standard method of solving such a partial differential equation is by separation of variables, where we express the solution as the product of functions containing each variable separately. In this case, we would write the temperature as\r\n<div style=\"text-align: center;\">[latex]\\large{T\\,(r,\\ t)=R\\,(r)\\,f\\,(t)}[\/latex].<\/div>\r\n&nbsp;\r\n<ol>\r\n \t<li>Substitute this form into our first equation from this section and, noting that [latex]f(t)[\/latex] is constant with respect to distance [latex](r)[\/latex] and [latex]R(r)[\/latex] is constant with respect to time [latex](t)[\/latex], show that\r\n<div style=\"text-align: center;\">[latex]\\large{\\frac{1}{f}\\frac{\\partial f}{\\partial t}=\\frac{K}{R}\\left[\\frac{\\partial^{2}R}{\\partial r^{2}}+\\frac{2}{r}\\frac{\\partial R}{\\partial r}\\right]}[\/latex].<\/div>\r\n&nbsp;<\/li>\r\n \t<li>This equation represents the separation of variables we want. The left-hand side is only a function of [latex]t[\/latex] and the right-hand side is only a function of [latex]r[\/latex], and they must be equal for all values of [latex]r[\/latex] and [latex]t[\/latex]. Therefore, they both must be equal to a constant. Let's call that constant [latex]-\\lambda^{2}[\/latex]. (The convenience of this choice is seen on substitution.) So, we have\r\n<div style=\"text-align: center;\">[latex]\\large{\\frac{1}{f}\\frac{\\partial f}{\\partial t}=-\\lambda^{2}}[\/latex] and [latex]\\large{\\frac{K}{R}\\left[\\frac{\\partial^{2}R}{\\partial r^{2}}+\\frac{2}{r}\\frac{\\partial R}{\\partial r}\\right]=-\\lambda^{2}}[\/latex].<\/div>\r\n&nbsp;\r\n\r\nNow, we can verify through direct substitution for each equation that the solutions are [latex]f\\,(t)=Ae^{-\\lambda^{2}t}[\/latex] and [latex]R\\,(r)=B\\left(\\frac{\\sin{\\alpha r}}{r}\\right)+C\\left(\\frac{\\cos{\\alpha r}}{r}\\right)[\/latex], where [latex]\\alpha=\\lambda \/\\sqrt{K}[\/latex]. Note that [latex]f\\,(t)=Ae^{+\\lambda n^{2}t}[\/latex] is also a valid solution, so we could have chosen [latex]+\\lambda^{2}[\/latex] for our constant. Can you see why it would not be valid for this case as time increases?<\/li>\r\n \t<li>Let's now apply boundary conditions.\r\n<ol>\r\n \t<li>The temperature must be finite at the center of Earth, [latex]r=0[\/latex]. Which of the two constants, [latex]B[\/latex] or [latex]C[\/latex], must therefore be zero to keep [latex]R[\/latex] finite at [latex]r=0[\/latex]? (Recall that [latex]\\sin{(\\alpha r)}\/r\\to a=[\/latex] as [latex]r\\to 0[\/latex], but [latex]\\cos{(\\alpha r)}\/r[\/latex] behaves very differently.)<\/li>\r\n \t<li>Kelvin argued that when magma reaches Earth\u2019s surface, it cools very rapidly. A person can often touch the surface within weeks of the flow. Therefore, the surface reached a moderate temperature very early and remained nearly constant at a surface temperature [latex]T_s[\/latex]. For simplicity, let's set [latex]T=0[\/latex] at [latex]r=R_E[\/latex] and find [latex]\\alpha[\/latex] such that this is the temperature there for all time [latex]t[\/latex]. (Kelvin took the value to be [latex]300\\text{ K}\\approx 80^{\\circ}F[\/latex]. We can add this [latex]300\\text{ K}[\/latex] constant to our solution later.) For this to be true, the sine argument must be zero at [latex]r=R_E[\/latex]. Note that [latex]\\alpha[\/latex] represents a valid solution (each with its own value for [latex]A[\/latex]). The total or general solution is the sum of all these solutions.<\/li>\r\n \t<li>At [latex]t=0[\/latex], we assume that all of Earth was at an initial hot temperature [latex]T_0[\/latex] (Kelvin took this to be about [latex]7000\\text{ K}[\/latex].) The application of this boundary condition involves the more advanced application of Fourier coefficients. As noted in part 2. each value of [latex]\\alpha_n[\/latex] represents a valid solution, and the general solution is a sum of all these solutions. This results in a series solution:\r\n<div style=\"text-align: center;\">[latex]\\large{T(r,t)=\\left(\\frac{T_0 R_E}{\\pi}\\right)\\displaystyle\\sum_{n}^{}\\frac{(-1)^{n-1}}{n}e^{-\\lambda n^2t}\\frac{\\sin{(\\alpha_nr)}}{r}}[\/latex], where [latex]\\large{\\alpha_n=n\\pi \/R_E}[\/latex].<\/div>\r\n&nbsp;<\/li>\r\n<\/ol>\r\n<\/li>\r\n<\/ol>\r\nNote how the values of [latex]\\alpha_n[\/latex] come from the boundary condition applied in part 2. The term [latex]\\frac{-1^{n-1}}{n}[\/latex] is the constant [latex]A_n[\/latex] for each term in the series, determined from applying the Fourier method. Letting [latex]\\beta=\\frac{\\pi}{R_E}[\/latex], examine the first few terms of this solution shown here and note how [latex]\\lambda^{2}[\/latex] in the exponential causes the higher terms to decrease quickly as time progresses:\r\n<div>\r\n\r\n[latex]\\hspace{4cm}T(r,t)=\\frac{T_0 R_E}{\\pi r}\\Bigg(\\begin{alignat}{2} e^{-K\\beta^{2}t}(\\sin{\\beta r})&amp;- \\frac{1}{2}e^{-4K\\beta^{2}t}(\\sin{2\\beta r})+\\frac{1}{3}e^{-9K\\beta^{2}t}(\\sin{3\\beta r})\\\\\r\n\r\n&amp;-\\frac{1}{4}e^{-16K\\beta^{2}t}(\\sin{4\\beta r})+\\frac{1}{5}e^{-25K\\beta^{2}t}(\\sin{5\\beta r})\\ldots&amp;\\quad\\\\\r\n\r\n\\end{alignat}\\Bigg).[\/latex]\r\n\r\n<\/div>\r\nNear time [latex]t=0[\/latex], many terms of the solution are needed for accuracy. Inserting values for the conductivity [latex]K[\/latex] and [latex]\\beta=\\pi \/R_E[\/latex] for time approaching merely thousands of years, only the first few terms make a significant contribution. Kelvin only needed to look at the solution near Earth's surface\u00a0(Figure 6) and, after a long time, determine what time best yielded the estimated temperature gradient known during his era ([latex]1^{\\circ}\\text{F}[\/latex] increase per [latex]50\\text{ ft}[\/latex]). He simply chose a range of times with a gradient close to this value. In\u00a0Figure 6, the solutions are plotted and scaled, with the [latex]300-\\text{K}[\/latex] surface temperature added. Note that the center of Earth would be relatively cool. At the time, it was thought Earth must be solid.\r\n\r\n[caption id=\"attachment_989\" align=\"aligncenter\" width=\"899\"]<img class=\"size-full wp-image-989\" src=\"https:\/\/s3-us-west-2.amazonaws.com\/courses-images\/wp-content\/uploads\/sites\/5667\/2021\/09\/28230809\/4-3-6.jpeg\" alt=\"This figure consists of two figures labeled a and b. Figure a shows three curves labeled 20, 50, and 200 million years on a chart showing fraction of the earth\u2019s radius vs. temperature (K). The highest curve is the 20 million one, then the 50 million one, and then the 200 million one, with all of them starting with a mildly decreasing slope until the slope decreases more steeply around x = 0.2 and then they all intersect at roughly (1, 315). Figure b shows a close up near (1, 315) with the x axis marked 4.0 miles below Earth\u2019s surface; the curves all appear linear in this close up, with the slopes increasing as the value of the curve does.\" width=\"899\" height=\"535\" \/> Figure 4. Temperature versus radial distance from the center of Earth. (a) Kelvin\u2019s results, plotted to scale. (b) A close-up of the results at a depth of [latex]\\small{4.0\\text{ mi}}[\/latex] below Earth\u2019s surface.[\/caption]<strong>Epilogue: <\/strong>On May 20, 1904, physicist <span id=\"term175\" class=\"no-emphasis\" data-type=\"term\">Ernest Rutherford<\/span> spoke at the Royal Institution to announce a revised calculation that included the contribution of radioactivity as a source of Earth\u2019s heat. In Rutherford\u2019s own words:\r\n<p id=\"fs-id1167793432220\" class=\" \">\u201cI came into the room, which was half-dark, and presently spotted Lord Kelvin in the audience, and realised that I was in for trouble at the last part of my speech dealing with the age of the Earth, where my views conflicted with his. To my relief, Kelvin fell fast asleep, but as I came to the important point, I saw the old bird sit up, open an eye and cock a baleful glance at me.<\/p>\r\n<p id=\"fs-id1167793432227\" class=\" \">Then a sudden inspiration came, and I said Lord Kelvin had limited the age of the Earth, <em data-effect=\"italics\">provided no new source<\/em> [<em data-effect=\"italics\">of heat<\/em>] <em data-effect=\"italics\">was discovered<\/em>. That prophetic utterance referred to what we are now considering tonight, radium! Behold! The old boy beamed upon me.\u201d<\/p>\r\n<p id=\"fs-id1167793432244\" class=\" \">Rutherford calculated an age for Earth of about [latex]500[\/latex] million years. Today's accepted value of Earth's age is about [latex]4.6[\/latex] billion years.<\/p>\r\n\r\n<\/div>\r\n&nbsp;","rendered":"<div class=\"textbox learning-objectives\">\n<h3>Learning Outcomes<\/h3>\n<ul>\n<li>Explain the meaning of a partial differential equation and give an example.<\/li>\n<\/ul>\n<\/div>\n<p>In <a href=\"https:\/\/courses.lumenlearning.com\/calculus2\/chapter\/introduction-4-2\/\" target=\"_blank\" rel=\"noopener\" data-book-uuid=\"1d39a348-071f-4537-85b6-c98912458c3c\" data-page-slug=\"4-introduction\">Introduction to Differential Equations<\/a>, we studied differential equations in which the unknown function had one independent variable. A <strong><span id=\"b1ffd73f-b9d7-41a1-8546-c25c366b2e90_term170\" data-type=\"term\">partial differential equation<\/span><\/strong> is an equation that involves an unknown function of more than one independent variable and one or more of its partial derivatives. Examples of partial differential equations are<\/p>\n<p style=\"text-align: center;\">[latex]\\large{u_t=c^{2}\\,(u_{xx}+u_{yy})}[\/latex]<\/p>\n<p>(heat equation in two dimensions)<\/p>\n<p style=\"text-align: center;\">[latex]\\large{u_{tt}=c^{2}\\,(u_{xx}+u_{yy})}[\/latex]<\/p>\n<p>(wave equation in two dimensions)<\/p>\n<p style=\"text-align: center;\">[latex]\\large{u_{xx}+u_{yy}=0}[\/latex]<\/p>\n<p>(Laplace&#8217;s equation in two dimensions)<\/p>\n<p>In the first two equations, the unknown function [latex]u[\/latex] has three independent variables[latex]-t,\\ x[\/latex], and [latex]y-[\/latex]and [latex]c[\/latex] is an arbitrary constant. The independent variables [latex]x[\/latex] and [latex]y[\/latex] are considered to be spatial variables, and the variable [latex]t[\/latex] represents time. In Laplace&#8217;s equation, the unknown function [latex]u[\/latex] has two independent variables [latex]x[\/latex] and [latex]y[\/latex].<\/p>\n<div class=\"textbox exercises\">\n<h3>Example: A Solution to the Wave Equation<\/h3>\n<p>Verify that<\/p>\n<p style=\"text-align: center;\">[latex]\\large{u\\,(x,\\ y,\\ t)=5\\sin{(3\\pi x)}\\sin{(4\\pi y)}\\cos{(10\\pi t)}}[\/latex]<\/p>\n<p>is a solution to the wave equation<\/p>\n<p style=\"text-align: center;\">[latex]\\large{u_{tt}=4\\,(u_{xx}+u_{yy})}.[\/latex]<\/p>\n<div class=\"qa-wrapper\" style=\"display: block\"><span class=\"show-answer collapsed\" style=\"cursor: pointer\" data-target=\"qfs-id1577773930324\">Show Solution<\/span><\/p>\n<div id=\"qfs-id1577773930324\" class=\"hidden-answer\" style=\"display: none\">\n<div style=\"text-align: left;\">First, we calculate [latex]u_{tt},\\ u_{xx}[\/latex], and [latex]u_{yy}[\/latex]:<\/div>\n<div style=\"text-align: center;\">[latex]\\begin{array}{ccc}\\hfill {u_{tt}} & =\\hfill & {\\frac{\\partial}{\\partial t}\\left[\\frac{\\partial u}{\\partial t}\\right]} \\hfill \\\\ \\hfill & =\\hfill & {\\frac{\\partial}{\\partial t}[5\\sin{(3\\pi x)}\\sin{(4\\pi y)}(-10\\pi\\sin{(10\\pi t)})]}\\hfill \\\\ \\hfill & =\\hfill & {\\frac{\\partial}{\\partial t}[-50\\pi\\sin{(3\\pi x)}\\sin{(4\\pi y)}\\sin{(10\\pi t)}]}\\hfill \\\\ \\hfill & =\\hfill & {-500\\pi^{2}\\sin{(3\\pi x)}\\sin{(4\\pi y)}\\cos{(10\\pi t)}}\\hfill \\\\ \\hfill {u_{xx}} & =\\hfill & {\\frac{\\partial}{\\partial x}\\left[\\frac{\\partial u}{\\partial x}\\right]} \\hfill \\\\ \\hfill & =\\hfill & {\\frac{\\partial}{\\partial x}[15\\pi\\cos{(3\\pi x)}\\sin{(4\\pi y)}\\cos{(10\\pi t)}]}\\hfill \\\\ \\hfill & =\\hfill & {-45\\pi^{2}\\sin{(3\\pi x)}\\sin{(4\\pi y)}\\cos{(10\\pi t)}}\\hfill \\\\ \\hfill {u_{yy}} & =\\hfill & {\\frac{\\partial}{\\partial y}\\left[\\frac{\\partial u}{\\partial y}\\right]}\\hfill \\\\ \\hfill & =\\hfill & {\\frac{\\partial}{\\partial y}[5\\sin{(3\\pi x)}(4\\pi\\cos{(4\\pi y)})\\cos{(10\\pi t)}]}\\hfill \\\\ \\hfill & =\\hfill & {\\frac{\\partial}{\\partial y}[20\\pi\\sin{(3\\pi x)}\\cos{(4\\pi y)}\\cos{(10\\pi t)}]}\\hfill \\\\ \\hfill & =\\hfill & {-80\\pi^{2}\\sin{(3\\pi x)}\\sin{(4\\pi y)}\\cos{(10\\pi t)}.}\\hfill \\\\ \\hfill \\end{array}[\/latex]<\/div>\n<div>Next, we substitute each of these into the right-hand side of our wave equation and simplify:<\/div>\n<div><\/div>\n<div>\n<p>[latex]\\text{ }\\\\ \\begin{alignat}{2} 4(u_{xx}+u_{yy}) &= 4(-45\\pi^2\\sin{(3\\pi x)}\\sin{(4\\pi y)}\\cos{(10\\pi t)} + -80\\pi^2\\sin{(3\\pi x)}\\sin{(4\\pi y)}\\cos{(10\\pi t)}\\\\    \\hspace{2cm}&=4(-125\\pi^2\\sin{(3\\pi x)}\\sin{(4\\pi y)}\\cos{(10\\pi t)} &\\quad\\\\    &=-500\\pi^2\\sin{(3\\pi x)}\\sin{(4\\pi y)}\\cos{(10\\pi t)} \\\\    &=u_{tt}.\\\\    \\end{alignat}[\/latex]<\/p>\n<p>This verifies the solution.<\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"textbox key-takeaways\">\n<h3>try it<\/h3>\n<p>Verify that [latex]u\\,(x,\\ y,\\ t)=2\\,\\sin\\,(\\frac{x}{3})\\,\\sin\\,(\\frac{y}{4})\\,e^{-25t\/16}[\/latex] is a solution to the heat equation<\/p>\n<p style=\"text-align: center;\">[latex]u_{t}=9(u_{xx}+u_{yy})[\/latex].<\/p>\n<div class=\"qa-wrapper\" style=\"display: block\"><span class=\"show-answer collapsed\" style=\"cursor: pointer\" data-target=\"qfs-id1367573930124\">Show Solution<\/span><\/p>\n<div id=\"qfs-id1367573930124\" class=\"hidden-answer\" style=\"display: none\">\n<div style=\"text-align: center;\"><\/div>\n<p>[latex]u\\,(x,\\ y,\\ t)=2\\,\\sin\\,(\\frac{x}{3})\\,\\sin\\,(\\frac{y}{4})\\,e^{-25t\/16}[\/latex] is a solution to the heat equation. See worked example video below to verify the solution.\n<\/p><\/div>\n<\/div>\n<\/div>\n<p>Watch the following video to see the worked solution to the above Try It<\/p>\n<div style=\"text-align: center;\"><iframe loading=\"lazy\" src=\"\/\/plugin.3playmedia.com\/show?mf=8186156&amp;p3sdk_version=1.10.1&amp;p=20361&amp;pt=375&amp;video_id=pMoH9Key-B8&amp;video_target=tpm-plugin-k0fmirvu-pMoH9Key-B8\" width=\"800px\" height=\"450px\" frameborder=\"0\" marginwidth=\"0px\" marginheight=\"0px\"><\/iframe><\/div>\n<div style=\"text-align: center;\">You can view the <a href=\"https:\/\/course-building.s3.us-west-2.amazonaws.com\/Calculus+3\/Calc+3+transcripts\/CP4.18_transcript.html\">transcript for \u201cCP 4.18\u201d here (opens in new window).<\/a><\/div>\n<p>Since the solution to the two-dimensional heat equation is a function of three variables, it is not easy to create a visual representation of the solution. We can graph the solution for fixed values of\u00a0[latex]t[\/latex], which amounts to snapshots of the heat distributions at fixed times. These snapshots show how the heat is distributed over a two-dimensional surface as time progresses. The graph of the preceding solution at time [latex]t=0[\/latex] appears in the following figure. As time progresses, the extremes level out, approaching zero as [latex]t[\/latex] approaches infinity.<\/p>\n<div id=\"attachment_985\" style=\"width: 466px\" class=\"wp-caption aligncenter\"><img loading=\"lazy\" decoding=\"async\" aria-describedby=\"caption-attachment-985\" class=\"size-full wp-image-985\" src=\"https:\/\/s3-us-west-2.amazonaws.com\/courses-images\/wp-content\/uploads\/sites\/5667\/2021\/09\/28230317\/4-3-3.jpeg\" alt=\"A complicated curve in xyz space with many sinusoidally alternating local maxima and minima.\" width=\"456\" height=\"379\" \/><\/p>\n<p id=\"caption-attachment-985\" class=\"wp-caption-text\">Figure 1.<\/p>\n<\/div>\n<p>If we consider the heat equation in one dimension, then it is possible to graph the solution over time. The heat equation in one dimension becomes<\/p>\n<p style=\"text-align: center;\">[latex]u_{t}=c^{2}u_{xx}[\/latex],<\/p>\n<p>where [latex]c^{2}[\/latex] represents the thermal diffusivity of the material in question. A solution of this differential equation can be written in the form<\/p>\n<p style=\"text-align: center;\">[latex]u_{m}\\,(x,\\ t)=e^{-\\pi^{2}m^{2}c^{2}t}\\sin{(m\\pi{x})}[\/latex]<\/p>\n<p>where [latex]m[\/latex] is any positive integer. A graph of this solution using [latex]m=1[\/latex] appears in\u00a0Figure 4, where the initial temperature distribution over a wire of length [latex]1[\/latex] is given by [latex]u\\,(x,\\ 0)=\\sin{\\pi{x}}[\/latex]. Notice that as time progresses, the wire cools off. This is seen because, from left to right, the highest temperature (which occurs in the middle of the wire) decreases and changes color from red to blue.<\/p>\n<div id=\"attachment_986\" style=\"width: 717px\" class=\"wp-caption aligncenter\"><img loading=\"lazy\" decoding=\"async\" aria-describedby=\"caption-attachment-986\" class=\"size-full wp-image-986\" src=\"https:\/\/s3-us-west-2.amazonaws.com\/courses-images\/wp-content\/uploads\/sites\/5667\/2021\/09\/28230533\/4-3-4.jpeg\" alt=\"A curve in xtu space with a local maximum at (0.5, 0, 12). From this maximum, the values decrease with increasing t and for any value of x.\" width=\"707\" height=\"648\" \/><\/p>\n<p id=\"caption-attachment-986\" class=\"wp-caption-text\">Figure 2.\u00a0Graph of a solution of the heat equation in one dimension over time.<\/p>\n<\/div>\n<div class=\"textbox tryit\">\n<h3 style=\"text-align: center;\" data-type=\"title\">Lord Kelvin and the Age of Earth<\/h3>\n<hr \/>\n<div id=\"attachment_987\" style=\"width: 741px\" class=\"wp-caption aligncenter\"><img loading=\"lazy\" decoding=\"async\" aria-describedby=\"caption-attachment-987\" class=\"size-full wp-image-987\" src=\"https:\/\/s3-us-west-2.amazonaws.com\/courses-images\/wp-content\/uploads\/sites\/5667\/2021\/09\/28230657\/4-3-5.jpeg\" alt=\"This figure consists of two figures marked a and b. Figure a show Lord Kelvin, dressed well and with a beard. Figure b shows an image of the planet Earth taken from space.\" width=\"731\" height=\"429\" \/><\/p>\n<p id=\"caption-attachment-987\" class=\"wp-caption-text\">Figure 3.\u00a0(a) William Thomson (Lord Kelvin), 1824-1907, was a British physicist and electrical engineer; (b) Kelvin used the heat diffusion equation to estimate the age of Earth (credit: modification of work by NASA).<\/p>\n<\/div>\n<p>During the late 1800s, the scientists of the new field of geology were coming to the conclusion that Earth must be \u201cmillions and millions\u201d of years old. At about the same time, Charles Darwin had published his treatise on evolution. Darwin\u2019s view was that evolution needed many millions of years to take place, and he made a bold claim that the Weald chalk fields, where important fossils were found, were the result of [latex]300[\/latex] million years of erosion.<\/p>\n<p>At that time, eminent physicist <span id=\"term174\" class=\"no-emphasis\" data-type=\"term\">William Thomson (Lord Kelvin)<\/span> used an important partial differential equation, known as the <em data-effect=\"italics\">heat diffusion equation<\/em>, to estimate the age of Earth by determining how long it would take Earth to cool from molten rock to what we had at that time. His conclusion was a range of [latex]20[\/latex] to [latex]400[\/latex] million years, but most likely about [latex]50[\/latex] million years. For many decades, the proclamations of this irrefutable icon of science did not sit well with geologists or with Darwin.<\/p>\n<div class=\"textbox\">\n<h3 style=\"text-align: center;\" data-type=\"title\">Interactive<\/h3>\n<hr \/>\n<p>Read Kelvin\u2019s <a href=\"http:\/\/www.openstax.org\/l\/20_KelEarthAge\" target=\"_blank\" rel=\"noopener nofollow\">paper<\/a> on estimating the age of the Earth.<\/p>\n<\/div>\n<p id=\"fs-id1167794213455\" class=\"\">Kelvin made reasonable assumptions based on what was known in his time, but he also made several assumptions that turned out to be wrong. One incorrect assumption was that Earth is solid and that the cooling was therefore via conduction only, hence justifying the use of the diffusion equation. But the most serious error was a forgivable one\u2014omission of the fact that Earth contains radioactive elements that continually supply heat beneath Earth\u2019s mantle. The discovery of radioactivity came near the end of Kelvin\u2019s life and he acknowledged that his calculation would have to be modified.<\/p>\n<p id=\"fs-id1167794213464\" class=\"\">Kelvin used the simple one-dimensional model applied only to Earth\u2019s outer shell, and derived the age from graphs and the roughly known temperature gradient near Earth\u2019s surface. Let\u2019s take a look at a more appropriate version of the diffusion equation in radial coordinates, which has the form<\/p>\n<div style=\"text-align: center;\">[latex]\\large{\\frac{\\partial{T}}{\\partial t}=K\\left[\\frac{\\partial^{2}T}{\\partial^{2}r}+\\frac{2}{r}\\frac{\\partial T}{\\partial r}\\right]}[\/latex].<\/div>\n<p>&nbsp;<\/p>\n<p>Here, [latex]T\\,(r,\\ t)[\/latex] is temperature as a function of [latex]r[\/latex] (measured from the center of Earth) and time [latex]t[\/latex]. [latex]K[\/latex] is the heat conductivity\u00a0conductivity\u2014for molten rock, in this case. The standard method of solving such a partial differential equation is by separation of variables, where we express the solution as the product of functions containing each variable separately. In this case, we would write the temperature as<\/p>\n<div style=\"text-align: center;\">[latex]\\large{T\\,(r,\\ t)=R\\,(r)\\,f\\,(t)}[\/latex].<\/div>\n<p>&nbsp;<\/p>\n<ol>\n<li>Substitute this form into our first equation from this section and, noting that [latex]f(t)[\/latex] is constant with respect to distance [latex](r)[\/latex] and [latex]R(r)[\/latex] is constant with respect to time [latex](t)[\/latex], show that\n<div style=\"text-align: center;\">[latex]\\large{\\frac{1}{f}\\frac{\\partial f}{\\partial t}=\\frac{K}{R}\\left[\\frac{\\partial^{2}R}{\\partial r^{2}}+\\frac{2}{r}\\frac{\\partial R}{\\partial r}\\right]}[\/latex].<\/div>\n<p>&nbsp;<\/li>\n<li>This equation represents the separation of variables we want. The left-hand side is only a function of [latex]t[\/latex] and the right-hand side is only a function of [latex]r[\/latex], and they must be equal for all values of [latex]r[\/latex] and [latex]t[\/latex]. Therefore, they both must be equal to a constant. Let&#8217;s call that constant [latex]-\\lambda^{2}[\/latex]. (The convenience of this choice is seen on substitution.) So, we have\n<div style=\"text-align: center;\">[latex]\\large{\\frac{1}{f}\\frac{\\partial f}{\\partial t}=-\\lambda^{2}}[\/latex] and [latex]\\large{\\frac{K}{R}\\left[\\frac{\\partial^{2}R}{\\partial r^{2}}+\\frac{2}{r}\\frac{\\partial R}{\\partial r}\\right]=-\\lambda^{2}}[\/latex].<\/div>\n<p>&nbsp;<\/p>\n<p>Now, we can verify through direct substitution for each equation that the solutions are [latex]f\\,(t)=Ae^{-\\lambda^{2}t}[\/latex] and [latex]R\\,(r)=B\\left(\\frac{\\sin{\\alpha r}}{r}\\right)+C\\left(\\frac{\\cos{\\alpha r}}{r}\\right)[\/latex], where [latex]\\alpha=\\lambda \/\\sqrt{K}[\/latex]. Note that [latex]f\\,(t)=Ae^{+\\lambda n^{2}t}[\/latex] is also a valid solution, so we could have chosen [latex]+\\lambda^{2}[\/latex] for our constant. Can you see why it would not be valid for this case as time increases?<\/li>\n<li>Let&#8217;s now apply boundary conditions.\n<ol>\n<li>The temperature must be finite at the center of Earth, [latex]r=0[\/latex]. Which of the two constants, [latex]B[\/latex] or [latex]C[\/latex], must therefore be zero to keep [latex]R[\/latex] finite at [latex]r=0[\/latex]? (Recall that [latex]\\sin{(\\alpha r)}\/r\\to a=[\/latex] as [latex]r\\to 0[\/latex], but [latex]\\cos{(\\alpha r)}\/r[\/latex] behaves very differently.)<\/li>\n<li>Kelvin argued that when magma reaches Earth\u2019s surface, it cools very rapidly. A person can often touch the surface within weeks of the flow. Therefore, the surface reached a moderate temperature very early and remained nearly constant at a surface temperature [latex]T_s[\/latex]. For simplicity, let&#8217;s set [latex]T=0[\/latex] at [latex]r=R_E[\/latex] and find [latex]\\alpha[\/latex] such that this is the temperature there for all time [latex]t[\/latex]. (Kelvin took the value to be [latex]300\\text{ K}\\approx 80^{\\circ}F[\/latex]. We can add this [latex]300\\text{ K}[\/latex] constant to our solution later.) For this to be true, the sine argument must be zero at [latex]r=R_E[\/latex]. Note that [latex]\\alpha[\/latex] represents a valid solution (each with its own value for [latex]A[\/latex]). The total or general solution is the sum of all these solutions.<\/li>\n<li>At [latex]t=0[\/latex], we assume that all of Earth was at an initial hot temperature [latex]T_0[\/latex] (Kelvin took this to be about [latex]7000\\text{ K}[\/latex].) The application of this boundary condition involves the more advanced application of Fourier coefficients. As noted in part 2. each value of [latex]\\alpha_n[\/latex] represents a valid solution, and the general solution is a sum of all these solutions. This results in a series solution:\n<div style=\"text-align: center;\">[latex]\\large{T(r,t)=\\left(\\frac{T_0 R_E}{\\pi}\\right)\\displaystyle\\sum_{n}^{}\\frac{(-1)^{n-1}}{n}e^{-\\lambda n^2t}\\frac{\\sin{(\\alpha_nr)}}{r}}[\/latex], where [latex]\\large{\\alpha_n=n\\pi \/R_E}[\/latex].<\/div>\n<p>&nbsp;<\/li>\n<\/ol>\n<\/li>\n<\/ol>\n<p>Note how the values of [latex]\\alpha_n[\/latex] come from the boundary condition applied in part 2. The term [latex]\\frac{-1^{n-1}}{n}[\/latex] is the constant [latex]A_n[\/latex] for each term in the series, determined from applying the Fourier method. Letting [latex]\\beta=\\frac{\\pi}{R_E}[\/latex], examine the first few terms of this solution shown here and note how [latex]\\lambda^{2}[\/latex] in the exponential causes the higher terms to decrease quickly as time progresses:<\/p>\n<div>\n<p>[latex]\\hspace{4cm}T(r,t)=\\frac{T_0 R_E}{\\pi r}\\Bigg(\\begin{alignat}{2} e^{-K\\beta^{2}t}(\\sin{\\beta r})&- \\frac{1}{2}e^{-4K\\beta^{2}t}(\\sin{2\\beta r})+\\frac{1}{3}e^{-9K\\beta^{2}t}(\\sin{3\\beta r})\\\\    &-\\frac{1}{4}e^{-16K\\beta^{2}t}(\\sin{4\\beta r})+\\frac{1}{5}e^{-25K\\beta^{2}t}(\\sin{5\\beta r})\\ldots&\\quad\\\\    \\end{alignat}\\Bigg).[\/latex]<\/p>\n<\/div>\n<p>Near time [latex]t=0[\/latex], many terms of the solution are needed for accuracy. Inserting values for the conductivity [latex]K[\/latex] and [latex]\\beta=\\pi \/R_E[\/latex] for time approaching merely thousands of years, only the first few terms make a significant contribution. Kelvin only needed to look at the solution near Earth&#8217;s surface\u00a0(Figure 6) and, after a long time, determine what time best yielded the estimated temperature gradient known during his era ([latex]1^{\\circ}\\text{F}[\/latex] increase per [latex]50\\text{ ft}[\/latex]). He simply chose a range of times with a gradient close to this value. In\u00a0Figure 6, the solutions are plotted and scaled, with the [latex]300-\\text{K}[\/latex] surface temperature added. Note that the center of Earth would be relatively cool. At the time, it was thought Earth must be solid.<\/p>\n<div id=\"attachment_989\" style=\"width: 909px\" class=\"wp-caption aligncenter\"><img loading=\"lazy\" decoding=\"async\" aria-describedby=\"caption-attachment-989\" class=\"size-full wp-image-989\" src=\"https:\/\/s3-us-west-2.amazonaws.com\/courses-images\/wp-content\/uploads\/sites\/5667\/2021\/09\/28230809\/4-3-6.jpeg\" alt=\"This figure consists of two figures labeled a and b. Figure a shows three curves labeled 20, 50, and 200 million years on a chart showing fraction of the earth\u2019s radius vs. temperature (K). The highest curve is the 20 million one, then the 50 million one, and then the 200 million one, with all of them starting with a mildly decreasing slope until the slope decreases more steeply around x = 0.2 and then they all intersect at roughly (1, 315). Figure b shows a close up near (1, 315) with the x axis marked 4.0 miles below Earth\u2019s surface; the curves all appear linear in this close up, with the slopes increasing as the value of the curve does.\" width=\"899\" height=\"535\" \/><\/p>\n<p id=\"caption-attachment-989\" class=\"wp-caption-text\">Figure 4. Temperature versus radial distance from the center of Earth. (a) Kelvin\u2019s results, plotted to scale. (b) A close-up of the results at a depth of [latex]\\small{4.0\\text{ mi}}[\/latex] below Earth\u2019s surface.<\/p>\n<\/div>\n<p><strong>Epilogue: <\/strong>On May 20, 1904, physicist <span id=\"term175\" class=\"no-emphasis\" data-type=\"term\">Ernest Rutherford<\/span> spoke at the Royal Institution to announce a revised calculation that included the contribution of radioactivity as a source of Earth\u2019s heat. In Rutherford\u2019s own words:<\/p>\n<p id=\"fs-id1167793432220\" class=\"\">\u201cI came into the room, which was half-dark, and presently spotted Lord Kelvin in the audience, and realised that I was in for trouble at the last part of my speech dealing with the age of the Earth, where my views conflicted with his. To my relief, Kelvin fell fast asleep, but as I came to the important point, I saw the old bird sit up, open an eye and cock a baleful glance at me.<\/p>\n<p id=\"fs-id1167793432227\" class=\"\">Then a sudden inspiration came, and I said Lord Kelvin had limited the age of the Earth, <em data-effect=\"italics\">provided no new source<\/em> [<em data-effect=\"italics\">of heat<\/em>] <em data-effect=\"italics\">was discovered<\/em>. That prophetic utterance referred to what we are now considering tonight, radium! Behold! The old boy beamed upon me.\u201d<\/p>\n<p id=\"fs-id1167793432244\" class=\"\">Rutherford calculated an age for Earth of about [latex]500[\/latex] million years. Today&#8217;s accepted value of Earth&#8217;s age is about [latex]4.6[\/latex] billion years.<\/p>\n<\/div>\n<p>&nbsp;<\/p>\n\n\t\t\t <section class=\"citations-section\" role=\"contentinfo\">\n\t\t\t <h3>Candela Citations<\/h3>\n\t\t\t\t\t <div>\n\t\t\t\t\t\t <div id=\"citation-list-865\">\n\t\t\t\t\t\t\t <div class=\"licensing\"><div class=\"license-attribution-dropdown-subheading\">CC licensed content, Original<\/div><ul class=\"citation-list\"><li>CP 4.18. <strong>Authored by<\/strong>: Ryan Melton. <strong>License<\/strong>: <em><a target=\"_blank\" rel=\"license\" href=\"https:\/\/creativecommons.org\/licenses\/by\/4.0\/\">CC BY: Attribution<\/a><\/em><\/li><\/ul><div class=\"license-attribution-dropdown-subheading\">CC licensed content, Shared previously<\/div><ul class=\"citation-list\"><li>Calculus Volume 3. <strong>Authored by<\/strong>: Gilbert Strang, Edwin (Jed) Herman. <strong>Provided by<\/strong>: OpenStax. <strong>Located at<\/strong>: <a target=\"_blank\" href=\"https:\/\/openstax.org\/books\/calculus-volume-3\/pages\/1-introduction\">https:\/\/openstax.org\/books\/calculus-volume-3\/pages\/1-introduction<\/a>. <strong>License<\/strong>: <em><a target=\"_blank\" rel=\"license\" href=\"https:\/\/creativecommons.org\/licenses\/by-nc-sa\/4.0\/\">CC BY-NC-SA: Attribution-NonCommercial-ShareAlike<\/a><\/em>. <strong>License Terms<\/strong>: Access for free at https:\/\/openstax.org\/books\/calculus-volume-3\/pages\/1-introduction<\/li><\/ul><\/div>\n\t\t\t\t\t\t <\/div>\n\t\t\t\t\t <\/div>\n\t\t\t <\/section>","protected":false},"author":349141,"menu_order":14,"template":"","meta":{"_candela_citation":"[{\"type\":\"cc\",\"description\":\"Calculus Volume 3\",\"author\":\"Gilbert Strang, Edwin (Jed) Herman\",\"organization\":\"OpenStax\",\"url\":\"https:\/\/openstax.org\/books\/calculus-volume-3\/pages\/1-introduction\",\"project\":\"\",\"license\":\"cc-by-nc-sa\",\"license_terms\":\"Access for free at https:\/\/openstax.org\/books\/calculus-volume-3\/pages\/1-introduction\"},{\"type\":\"original\",\"description\":\"CP 4.18\",\"author\":\"Ryan Melton\",\"organization\":\"\",\"url\":\"\",\"project\":\"\",\"license\":\"cc-by\",\"license_terms\":\"\"}]","CANDELA_OUTCOMES_GUID":"","pb_show_title":"on","pb_short_title":"","pb_subtitle":"","pb_authors":[],"pb_section_license":""},"chapter-type":[],"contributor":[],"license":[],"class_list":["post-865","chapter","type-chapter","status-publish","hentry"],"part":22,"_links":{"self":[{"href":"https:\/\/courses.lumenlearning.com\/calculus3\/wp-json\/pressbooks\/v2\/chapters\/865","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/courses.lumenlearning.com\/calculus3\/wp-json\/pressbooks\/v2\/chapters"}],"about":[{"href":"https:\/\/courses.lumenlearning.com\/calculus3\/wp-json\/wp\/v2\/types\/chapter"}],"author":[{"embeddable":true,"href":"https:\/\/courses.lumenlearning.com\/calculus3\/wp-json\/wp\/v2\/users\/349141"}],"version-history":[{"count":141,"href":"https:\/\/courses.lumenlearning.com\/calculus3\/wp-json\/pressbooks\/v2\/chapters\/865\/revisions"}],"predecessor-version":[{"id":5905,"href":"https:\/\/courses.lumenlearning.com\/calculus3\/wp-json\/pressbooks\/v2\/chapters\/865\/revisions\/5905"}],"part":[{"href":"https:\/\/courses.lumenlearning.com\/calculus3\/wp-json\/pressbooks\/v2\/parts\/22"}],"metadata":[{"href":"https:\/\/courses.lumenlearning.com\/calculus3\/wp-json\/pressbooks\/v2\/chapters\/865\/metadata\/"}],"wp:attachment":[{"href":"https:\/\/courses.lumenlearning.com\/calculus3\/wp-json\/wp\/v2\/media?parent=865"}],"wp:term":[{"taxonomy":"chapter-type","embeddable":true,"href":"https:\/\/courses.lumenlearning.com\/calculus3\/wp-json\/pressbooks\/v2\/chapter-type?post=865"},{"taxonomy":"contributor","embeddable":true,"href":"https:\/\/courses.lumenlearning.com\/calculus3\/wp-json\/wp\/v2\/contributor?post=865"},{"taxonomy":"license","embeddable":true,"href":"https:\/\/courses.lumenlearning.com\/calculus3\/wp-json\/wp\/v2\/license?post=865"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}