Yesterday I got (but I haven't tested it numerically) that $$\sum_{n=1}^\infty \sum_{k=1}^\infty\frac{G_n}{k}\frac{\Gamma(k+\frac{3}{2})\Gamma(n+\frac{1}{2})}{(n+k+1)!}=\frac{\pi(1-8\log 2)}{8},$$ where $G_n$ denotes the sequence of Gregory coefficients or reciprocal logarithmic numbers, see the Wikipedia Gregory coefficients. Using the generating function for this sequence and the result of the Problem 4383 from Crux Mathematicorum ([1]), I can then write the closed-form of a series of the type

$$\sum_{n=1}^\infty \sum_{k=1}^\infty(\text{a function of }G_n\text{, or }|G_n|,\text{ and }k)\cdot(\text{particular values of the gamma function}).$$

Question. Are known series of previous type and how are evaluated? If it is in the literature refer it, and I try to search and read those statements from the literature. In other case, what work can be done to deduce similar formulas to get the closed-form for double series of the previous type? Many thanks.

Here I add my example of series, how I got the result for the series of the first paragraph.

Example. All issues of convergence are satisfied by uniform convergence of the series in the unit interval. Multiplying the Maclaurin series of the Gregory coefficients by $\sqrt{\frac{x}{1-x}}\log x$ and integrating the equation over the unit interval, one gets invoking the Problem 4383, due to Michel Bataille, that $$\sum_{n=1}^\infty G_n\left(\int_0^1(\log x)\sqrt{\frac{x}{1-x}}x^n dx\right)=\frac{\pi(-1+8\log 2)}{8}.\tag{1}$$ We change the varialbe of the integral $x=1-z$ and use the Taylor series of the logarithm to state that $$\sum_{n=1}^\infty \sum_{k=1}^\infty\frac{G_n}{k}\left(\int_0^1 z^k (1-z)^n \sqrt{\frac{1-z}{z}}dz\right)=\frac{\pi(1-8\log 2)}{8},\tag{2}$$ and finally I've used Wolfram Alpha online calculator to get a closed-form of the integral in terms of particular values of the gamma function.

References:

[1] Problem 4383, Crux Mathematicorum, Volume 45, Number 6, July 2019.