From b373fc9f50057de52b8e3cd6529b33cf4c48ac45 Mon Sep 17 00:00:00 2001
From: ChaofanChen <chaofan.chen@ufz.de>
Date: Mon, 10 Feb 2020 13:15:57 +0100
Subject: [PATCH] [1P/web] Added benchmark description.
---
.../Analytical_wellbore_heat_transport.zip | Bin 0 -> 1358 bytes
.../heat-transport-bhe/pipe_flow_EBHE.pandoc | 101 ++++++++++++++++++
.../T_out_comparison.png | 3 +
.../absolute_error_fluid_T_30d.png | 3 +
.../pipe_flow_3d_model.png | 3 +
5 files changed, 110 insertions(+)
create mode 100644 web/content/docs/benchmarks/heat-transport-bhe/Analytical_wellbore_heat_transport.zip
create mode 100644 web/content/docs/benchmarks/heat-transport-bhe/pipe_flow_EBHE.pandoc
create mode 100644 web/content/docs/benchmarks/heat-transport-bhe/pipe_flow_EBHE_figures/T_out_comparison.png
create mode 100644 web/content/docs/benchmarks/heat-transport-bhe/pipe_flow_EBHE_figures/absolute_error_fluid_T_30d.png
create mode 100644 web/content/docs/benchmarks/heat-transport-bhe/pipe_flow_EBHE_figures/pipe_flow_3d_model.png
diff --git a/web/content/docs/benchmarks/heat-transport-bhe/Analytical_wellbore_heat_transport.zip b/web/content/docs/benchmarks/heat-transport-bhe/Analytical_wellbore_heat_transport.zip
new file mode 100644
index 0000000000000000000000000000000000000000..c0082f1983d190975c842c85a2976d7a1da7171d
GIT binary patch
literal 1358
zcmWIWW@Zs#-~hr^nO*@5P@uudz@WmQz~Gpdm{VDjnVgstU!IzilaybS8lRDxSQ1}S
zl$cjskY7}ymm3<w%fSB7!!hnD5QFH_3T_5QmKV$n3}B)+G&ldYnZUoi1OFT2dAz%F
z0w>*{w1&A&l%tK&RbWZR#9V`v&P^Q`7C8m^ZI0h>9W?Ffg1cf0+fG)$+xh)*%b%e1
zM`u6S2)lS+l+w!%@^bY#dO}+`NnKP~^fZs5*6LE_sgnI_M>aSg4oQoeaB9+)(t}$A
zYpzR$E%>MCyzJwP!srS7DLeQ4?GswhYP>{5S6|=+<BR}qorq+uP03FCgzl~~I(@U^
z3(Kb~C%r`&CNAB`;&5%Z^dx7gmmgoSioBQ;G^=FE)@RSwhI8JYCH-knKi|uWC0aXm
z8RePpEGj)*6!+=Q(Rt~<)A*;he>=9g=61we*P7m2?biy=K4zW0bbE?ZRq#Z%ohSHq
zD{YprijVu=ps(=%@7&eZY70ds_&bVR_K{T$KNstvbmPXK<@@iaJ4}~7T(x&2k6++5
z*V)o_E>B}N329iaKHJ~9D?umRC_X>qPyNd&IlL!dHJR*5w7qyQEp_g^&sPd&`|NQ{
zio5RP_?7!!))U8E@8WgsMOGzp>emflt$cfKuDERH%8Gw6#<`LkuFvIvw#4IQ(dKs>
zE}7}_wLN5KVNBW_6(Dt^%!T7fWs1b}f-tXwVu8&^Q<d^r&Gc@`WK4~bF<5Qj8adaP
zBk1#^Csze|_PN9zX689l_rvTs|C_bb7F}49$8XJcwBmr)jX1u(#_1LrEq^(FRI@tW
zp3PV&D?Kq-J8D<_>aeGx;!HNdU7c(0vn_qJ(>i76M6L_l4lZ`vXY};w#O3i*KHkX`
zw%uT`?B7&tyA2w#x|+PguYNB-oigu8<Amxd>uO!v3%8Uw-n}}9v*Z8F_#FpLGfj2G
zl{qg5aQHm>&BC{dH)Q^;w{scH_Zr<|NyuzF{FeRuss24bwA#{m*W6b4)_m?~wWwO)
ztMf_s%Xua!$h@iC6&3ORhM}<Tt$$at&bGV^s`oCFtYYs|I-sYf;TGq)a-sEXx07y#
z7mDv*YB!ZD<o?sCqPg_3V59wpTMfSHGCf8MZ&sh2>)^=Sy}7BoZXt`9(8S=ECo{L-
zc~In8u)D{7{>GAk>hqh6i(Gr&3YBwtJUx81l|#=~=bGxhuh%Nw)^j8!Tw25(!+-Pk
z6sP9}2kU-w*iA8r*lJ!oXLf#?lwbVvxiY`jHk?U|mNBWCyiQ}Up&n!U&HG_u_iZc>
zKB@U#(p~)L;~cRFU&r$tiRGs+^EUo%e0{LaihqHM#|LAbxvA;be?7YRvuRyW2hTM#
zjf6c6Uvit0s`oN6SCj^tE$2P&xT%vdYc`*wUoKnL%MZ8gr<U3L>I#d$VE)^#wD>jG
zCH8XHU8ic>>Nl`1^?5LTuarj>tN4@0VN90$&Ijby)TUQBWPMxvF>QH`+mqFY4>R9q
z@cR{U(AMT*9FM)UV(;~*W-LsSlUC01cTCzL$g|42=5zlB!NT@Mdzy~2K7Lg5^Mm{O
zi;~+?>K`s@+_}>C<<4!K%Y`NG_<HK?`?qHK!^Q6><r<qy8kM<Bu1~cJjQDCKZm>;O
z;P8L;0B?4V#Chv~hp;d(fXb=>Z$>5&W<)uLEC(vbU|>ljh((Q3Ex?<V4P+K05EcRH
KV?dWMFaQ88lT-l!
literal 0
HcmV?d00001
diff --git a/web/content/docs/benchmarks/heat-transport-bhe/pipe_flow_EBHE.pandoc b/web/content/docs/benchmarks/heat-transport-bhe/pipe_flow_EBHE.pandoc
new file mode 100644
index 00000000000..d253b4ee351
--- /dev/null
+++ b/web/content/docs/benchmarks/heat-transport-bhe/pipe_flow_EBHE.pandoc
@@ -0,0 +1,101 @@
++++
+author = "Chaofan Chen, Haibing Shao"
+date = "2020-02-24T13:44:00+01:00"
+title = "Wellbore Heat Transport - EBHE"
+weight = 123
+project = "Parabolic/T/BHE_1P/BHE_1P.prj"
+
+[menu]
+ [menu.benchmarks]
+ parent = "heat-transport-bhe"
+
++++
+
+{{< data-link >}}
+
+## Problem description
+
+Ramey (Ramey et al. (1962)) proposed the analytical solution concerning the wellbore heat transmission, which can be used to quantify the fluid temperature change in the wellbore. In order to verify the single pipe flow model in the OGS, the numerical results was compared with the [Ramey's analytical solution](../Analytical_wellbore_heat_transport.zip). The detailed calculation of the Ramey's analytical solution is given below.
+
+## Model Setup
+
+In this benchmark, the length of the wellbore is 30 m as shown in Figure 1 and the cold water is injected into the inlet point of the wellbore with temperature of 20 $^{\circ}$C. The initial temperature of the fluid and grout in the wellbore is 20 $^{\circ}$C, and temperature of the surrounding rock is 55 $^{\circ}$C. The wellbore and pipe diameter are 0.28 m and 0.25826 m, respectively. And the flow rate is 0.0002 $m^3/s$.
+
+{{< img src="../pipe_flow_EBHE_figures/pipe_flow_3d_model.png" width="80">}}
+
+Figure 1: Single pipe flow model
+
+## Ramey's analytical solution
+
+In Ramey's analytical solution (Ramey et al. (1962)), the outlet temperature of the pipe inside the wellbore can be calculated by,
+
+\begin{equation}
+ T_o(t) = T_{s} + (T_i(t) - T_{s})\exp(-\Delta z/X)
+\end{equation}
+
+\noindent where, $q$ is the flow rate of the fluid in the wellbore and coefficient $X$ is determined by,
+
+\begin{equation}
+ X = \frac{q\rho_fc_{p,f}(\lambda_{s}+r_pUf(t))}{2\pi r_pU \lambda_{s}}
+\end{equation}
+
+with dimensionless time $t_D = \frac{\lambda_{s}t}{(\rho_{s}c_{p,s}r_b)}$, the time function $f(t)$ can be calculated by,
+
+\begin{align}
+ &f(t) = [0.4063+0.5\ln(t_D)][1+\frac{0.6}{t_D}], & t_D > 1.5
+ \\
+ &f(t) = 1.1281\sqrt{t_D}(1-0.3\sqrt{t_D}), & t_D \leqslant 1.5
+\end{align}
+
+\noindent and the overall heat transfer coefficient $U$ is written as follows,
+
+\begin{equation}
+ U = [\frac{r_{pi}+t_{pi}}{r_{pi}h}+(r_{pi}+t_{pi})(\frac{\ln{\frac{r_{pi}+t_{pi}}{r_{pi}}}}{\lambda_{pi}}+\frac{\ln{\frac{r_b}{r_{pi}+t_{pi}}}}{\lambda_{grout}})]^{-1}
+\end{equation}
+
+\begin{equation}
+ h = \frac{\lambda_f Nu}{2r_{pi}}
+\end{equation}
+
+The Nusselt number can be determined according to the Gnielinski's equation (Gnielinski et al. (1975)),
+
+\begin{align}
+ & Nu = 4.364, & Re < 2300 \\
+ & Nu = \frac{\frac{f}{8}(Re - 1000)Pr}{1+12.7\sqrt{\frac{f}{8}}(Pr^{\frac{2}{3}}-1)}, & 2300\leqslant Re < 5 \times 10^6
+\end{align}
+
+Pr is the Prandtl number, and the friction factor $f$, is evaluated by Churchill correlation (Churchill et al. (1977)),
+
+\begin{equation}
+ f = \frac{1}{(\frac{1}{[((\frac{8}{Re})^{10}+(\frac{Re}{36500})^{20})]^{1/2}}+[2.21(\ln{\frac{Re}{7}})]^{10})^{1/5}}
+\end{equation}
+
+The Prandtl and Reynolds number can be calculated as follows,
+
+\begin{align}
+ & Pr = \frac{\mu_f c_{p,f}}{\lambda_f}
+ & Re = \frac{\rho_f v d_{pi}}{\mu_f}
+\end{align}
+\noindent where, $\mu_f, \rho_f$ and $\lambda_f$ is the fluid viscosity, density and thermal conductivity.
+
+## Results and discussion
+
+The outlet temperature change over time was compared against analytical solution and presented in Figure 2. After 30 days, the fluid temperature distribution in the wellbore is shown in Figure 3. The maximum relative error between the numerical model and Ramey's analytical solution is less than 0.15 \%.
+
+In numerical model, the outlet temperature at beginning stage is affected by the initial temperature in the pipe inside the wellbore. The initial fluid temperature set in the benchmark means there is water with 20 $^{\circ}$C filled in the wellbore already before injecting water into the wellbore. But in the analytical solution, no initial temperature is set and the temperature keeps equilibrium state at every moment. The impact of initial temperature condition in numerical model is decreasing with increasement of the operational time as shown in Figure 2.
+
+{{< img src="../pipe_flow_EBHE_figures/T_out_comparison.png" width="120">}}
+
+Figure 2: Comparison with analytical solution results
+
+{{< img src="../pipe_flow_EBHE_figures/absolute_error_fluid_T_30d.png" width="200">}}
+
+Figure 3: Distributed temperature of fluid and absolute error.
+
+## References
+
+[1] Ramey Jr, H. J. (1962). Wellbore heat transmission. Journal of petroleum Technology, 14(04), 427-435.
+
+[2] Gnielinski, V. (1975). New equations for heat and mass transfer in the turbulent flow in pipes and channels. NASA STI/recon technical report A, 75, 8-16.
+
+[3] Churchill, S. W. (1977). Comprehensive correlating equations for heat, mass and momentum transfer in fully developed flow in smooth tubes. Industrial & Engineering Chemistry Fundamentals, 16(1), 109-116.
diff --git a/web/content/docs/benchmarks/heat-transport-bhe/pipe_flow_EBHE_figures/T_out_comparison.png b/web/content/docs/benchmarks/heat-transport-bhe/pipe_flow_EBHE_figures/T_out_comparison.png
new file mode 100644
index 00000000000..c6da5181330
--- /dev/null
+++ b/web/content/docs/benchmarks/heat-transport-bhe/pipe_flow_EBHE_figures/T_out_comparison.png
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:07e9a3c4114e57f5aa671423ec59e2dcf07b838cd67daa8d53ce9ce266ded4a4
+size 52500
diff --git a/web/content/docs/benchmarks/heat-transport-bhe/pipe_flow_EBHE_figures/absolute_error_fluid_T_30d.png b/web/content/docs/benchmarks/heat-transport-bhe/pipe_flow_EBHE_figures/absolute_error_fluid_T_30d.png
new file mode 100644
index 00000000000..1a42a0eb177
--- /dev/null
+++ b/web/content/docs/benchmarks/heat-transport-bhe/pipe_flow_EBHE_figures/absolute_error_fluid_T_30d.png
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:d10c3715c55cc99d5f1511c2642badeae60bfe00ccc8378bb8760cecc7f789ed
+size 82189
diff --git a/web/content/docs/benchmarks/heat-transport-bhe/pipe_flow_EBHE_figures/pipe_flow_3d_model.png b/web/content/docs/benchmarks/heat-transport-bhe/pipe_flow_EBHE_figures/pipe_flow_3d_model.png
new file mode 100644
index 00000000000..768c58fcbee
--- /dev/null
+++ b/web/content/docs/benchmarks/heat-transport-bhe/pipe_flow_EBHE_figures/pipe_flow_3d_model.png
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:05ab5d28b42e3f280ec05966371cada8406fafdaa1ea4f6acddf3f7d80fd9060
+size 46647
--
GitLab