I ran the numbers for a Carnot engine running between poles and equator and estimate diameter of the tunnel after accounting for pumping costs.
Calculation is done for an engine generating 1 Terawatt power, running at a temperature difference of 90 degrees, operating between 5 and 95 degrees. We assume 20 pc of the energy generated is used for pumping.
I will write down all the formulas so that it can be checked.
For the calculation I will assume a
friction factor of .01. (The actual friction factor is slightly lower, which we will calculate at the end. We are overestimating friction by ~30pc, to leave some headroom)
$$ f = .01$$
The head loss per unit length is given by the formula
$$\frac{\Delta p}{L} = \frac{f \rho v^2}{2 d} $$
We assume the engine works at
efficiency at peak power.( ## T_h = 95 ## and ## T_c = 5 ##)
$$ \eta = 1- \sqrt{\frac{T_c}{T_h}} = .131$$
The potential to do work per Kg of water is ( c = specific heat of water )
$$ w = c \eta \Delta T = 49300 J/Kg$$
A fraction of the energy ##p_f = .2 ## is used for pumping. We assume pumping efficiency is 75 pc. In one cycle work can be extracted twice from the same mass, one at the pole when cooling, and another at the equator when heating, we will assume it is identical. (## h_{head} = \frac{p_f w}{g} ##, L = 10000 Km} )
$$ \frac{\Delta p}{L} = .75 \frac{p g h_{head} }{L} = .739 pascal/m $$
mass flow rate M to generate 1 Terawatt of power is ## M = \frac{10^{12}}{w (1-p_f)} = 2.5 \times 10^7 kg/s## and Volume flow rate ##Q = \frac{M}{\rho} = 25\times 10^3 m^3/s##
The formula $$\frac{\Delta p}{L} = \frac{f \rho v^2}{2 d} = \frac{ 8 f \rho Q^2 }{\pi^2 d^5}$$ can be inverted to get the diameter
$$d = \sqrt[5]{\frac{ 8 f Q^2 \rho }{ \pi^2 (\Delta P/L) } } = 93 m$$
velocity of flow is ## v= \frac{Q}{\pi d^2} = .92 m/s ##
We use
Colebrook-White equation iteratively to estimate the friction factor for these flow parameters, with surface roughness ##e = .25 mm ## and reynolds number ##Re = 3.9 \times 10^8##. We get ##f = 0.0067 ## In principle the calculation can iterated over, and if that is done we get ##d = 86m##. But some head room may be useful.
I will repeat the same for other cases of interest for different kinds of cases and make general comments on a separate post. But is the general method OK for the problem as stated?