7

Based on this question about an impossible barbell planet, how hot would the bar be if it extended from halfway in the Earth's mantle to an altitude of 2r of Earth's radius? The diameter of the bar itself is 20% of the radius of Earth. The rod is made of mild steel (we are handwaving the structural integrity of molten steel and ignoring that a structure this high would instantly collapse.) Also, the rod is permanently affixed at halfway through the mantle. The orientation and distance between Earth's core and the bar are fixed points.

Barbell Planet, Broo

We are going to ignore orbital mechanics, structure strength of megastructures, and structural strength of the Earth's crust surrounding this megastructure.

Given the Earth's core temperature and the conductivity of mild steel, how high would a person have to go on the rod in order to touch the rod with their bare hands and not get burned?

Bonus Question(s): What would happen to the weather near the base of the rod where it intersects with earth? If you want to talk about how fast Earth's core would completely cool off, that's cool too.

Extra Bonus Question(s): Account for orbital mechanics, interaction with the magnetosphere, and compression heating of the rod against the earth's core and friction with the crust.

Remember, this is a question. Equations, official government sources, and journal references are most appreciated. A minimum of handwavium, please.

Green
  • 52,750
  • 10
  • 130
  • 260
  • 1
    This question would be more suitable for physics.SE. – user6760 Jul 29 '15 at 04:32
  • @user6760 I think it's fine here. – HDE 226868 Jul 29 '15 at 14:23
  • 1
    I suppose you're ignoring the melting point of steel and the fact that the core can't rotate with a rod sticking out of it up through the crust? – Samuel Jul 29 '15 at 17:13
  • @Samuel, the melting point of steel is ignored along the length of the rod and let's suspend the rod halfway through the mantle. The thermal dynamics for a question like this are complicated enough, so no need to further complicate the question with interactions with Earth's magnetosphere. Below the surface, I expect the rod will be the same temperature as the surrounding mantle so I don't think it needs to be stuck to the core. – Green Jul 29 '15 at 18:54
  • 1
    @Green Yeah, I started to calculate the heat transfer and temperature in terms of distance from the core already. It is complex, the rod is too thick to be considered isothermal for a given cross-section and the shift from rock, to soil+water, to air, to vacuum (and back!) is a lot to deal with. Hard-science is the perfect description for this and I have abandoned my answer (no time this week, not that much!). – Samuel Jul 29 '15 at 19:05
  • @Samuel, I figured that either you or HDE would be able to answer this question. Take your time! :) This question is much more a fun puzzle than anything I desperately need. And considering that to solve this in a real life circumstance would require some extensive CFD model building and compute time, the roughest of answers will suffice. – Green Jul 29 '15 at 19:16
  • The rod would not stop the planet from rotating if it was in one of the poles. – Loren Pechtel Jul 30 '15 at 03:59
  • @LorenPechtel I suppose you could but I remember seeing a child's toy where a short wooden dowel is inserted into wooden sphere. The rod is then spun and dropped on a hard surface. In a short period the rod and sphere are inverted and instead of spinning with the sphere side down, it's rod side down. I don't understand/know the physics behind why this happens but I want to avoid as much complexity as possible. – Green Jul 30 '15 at 15:58
  • Unless the planets are also made out of unobtanium, you'd end up with one big planet with a skewer in the middle, as they would just slide over the bar. – Serban Tanasa Jul 30 '15 at 19:47
  • @SerbanTanasa, thank you for pointing that out. I've fixed Earth and the bar in a constant reference frame. – Green Jul 30 '15 at 19:52
  • 3
    "We're going to hand wave all these things. Oh yeah, and hard science, a minimum of handwavium, please." – AndyD273 Aug 03 '15 at 17:35

1 Answers1

9

As a first approximation, we treat the rod as a one-dimensional rod of length $L$. On each end is a heat source (the mantles of each of the two "Earths". We can begin to model the system as following the one-dimensional heat equation: $$\frac{\partial T}{\partial t}=k\frac{\partial^2T}{\partial x^2}\tag{1}$$ where $T$ represents temperature, $x$ is the distance from one end of the rod, $t$ is time and $k$ is a constant, the thermal diffusivity. Consider a case with the following boundary conditions: $$T(0,t)=T_1,\quad T(L,t)=T_2,\quad T(x,0)=f(x)$$ The first condition is because both heat sources are identical. We can use the method of separation of variables to arrive at our solution: $$T(x,t)=T_1+\frac{T_2-T_1}{L}x+\sum_{n=1}^{\infty}B_n\sin\left(\frac{n\pi x}{L}\right)e^{-k\left(\frac{n\pi}{L}\right)^2t}$$ where $$B_n=\frac{2}{L}\int_0^L(f(x)-u_E(x))\sin\left(\frac{n\pi x}{L}\right)dx,\quad u_E(x)=T_1+\frac{T_2-T_1}{L}x$$ In our case, $T_1=T_2$ and $f(x)=T_0$ is uniform. Then, if $\Delta T=T_0-T_1$, $$B_n=\frac{2\Delta T(1-\cos(\pi n))}{\pi n}$$ and $$T(x,t)=T_1+2\Delta T\sum_{n=1}^{\infty}\left(\frac{1-\cos(\pi n)}{\pi n}\right)\sin\left(\frac{n\pi x}{L}\right)e^{-k\left(\frac{n\pi}{L}\right)^2t}$$ Wikipedia cites steel as having roughly $k=2\times10^{-5}\text{ m}^2\text{ s}^{-1}$. Let's say $L=4R_{\oplus}$, $T_1=T_2=6000\text{ K}$ (according to estimates), and $T_0=300\text{ K}$.

It turns out that the bridge between the cores heats up really slowly. I computed the first 100 terms of $T(x,t)$ at a variety of times, and plotted them.

Temperature plot

Second temperature plot

The strong oscillations at the ends are just examples of the Gibbs phenomenon, and don't have physical significance.

It becomes dangerous to walk on the rod after timescales longer than about 100 million to 1 billion years, from what I can tell. The center should remain habitable.

There are a couple of things we haven't considered:

  • Radiative cooling, which I think will be important. This also might mean that the ends of the rods - at least, the sections immediately protruding from the surfaces of the planets - could appear quite bright. After all, 6000 K is approximately the temperature of the surface of the Sun!
  • The fact that the rod is a cylinder, not one-dimensional. I don't think this is a major factor when it comes to the relevant timescales, though, especially if the rod is thin.
  • The ends of the rod are embedded inside planets, and heat will diffuse into the outer layers of the planets. There should also be a non-zero gradient at the cores, at $t=0$.
  • Atmospheres have an impact on temperature circulation; the (cooler) sections of the atmosphere around the two planets should take in heat from the (hotter) sections around the rod.

That said, I think this simple 1-D model can give us an order-of-magnitude estimate.

HDE 226868
  • 101,188
  • 25
  • 303
  • 542
  • I'm not physicist or mechanical engineer but I think your approach is sound. – Green Jul 30 '15 at 20:51
  • 1
    On the one hand, you're ignoring radiative cooling once the rod leaves the crust. On the other hand, you're misapplying the geothermal gradient, since heat will diffuse away from the rod into the crust as soon as the rod temperature gets above the gradient temp. On the gripping hand, atmospheric temperatures aren't enough: you also need to factor in both the changing density and the added cooling caused by the massive convective currents which will accompany the increased temperatures. – WhatRoughBeast Jul 31 '15 at 03:25
  • 2
    It's definitely heading in the right direction though :) – Tim B Jul 31 '15 at 08:50
  • I already mentioned in the chat, but for those who don't go there, you are missing a few points (so far), like day-night variations, climate influence and seasonal variations. All those will make some differences. But, and maybe more importantly, the rod will influence its environment. Intuitively, it might start to cool the Earth, and affect the climate. At first close by and with butterfly effect on the whole surface of the Earth. Those will affect any values you actually put in your model. So as Tim has it, yes it a good starting point. The question might be a bit demanding, though... – clem steredenn Jul 31 '15 at 09:46
  • When I started to calculate this I found this document incredibly helpful. You can treat the rod like a pin fin in a fluid (except in space you stop subtracting the Q lost to convection out the sides of the rod). – Samuel Jul 31 '15 at 16:33
  • @WhatRoughBeast This is a dashed-off approximation. I was away for two days, and I certainly don't intend to leave it as is. – HDE 226868 Aug 02 '15 at 14:10
  • @TimB Thanks. See my comment to WhatRoughBeast. – HDE 226868 Aug 02 '15 at 14:10
  • @bilbo_pingouin See my comment to WhatRoughBeast. – HDE 226868 Aug 02 '15 at 14:10