$\begingroup$

This is not an answer, but an experiment for moderately large $n$, making the limit of $\approx 1.5$ unlikely. $$S_n=\sup\left\{\sum_{k=1}^{n}\frac{k}{\sum_{j=1}^{k}a_j} : a_1,\ldots,a_n>0,\ \sum_{k=1}^{n}\frac{1}{a_k}=1\right\}$$ The method of Lagrange works here. Writing $$\mathcal{L}(a_1,\ldots,a_n,\lambda)=\sum_{k=1}^{n}\frac{k}{\sum_{j=1}^{k}a_j}-\lambda\sum_{k=1}^{n}\frac{1}{a_k},\\\frac{\partial\mathcal{L}}{\partial a_m}=\frac{\lambda}{a_m^2}-\sum_{k=m}^{n}\frac{k}{\left(\sum_{j=1}^{k}a_j\right)^2}\quad(1\leqslant m\leqslant n),$$ we obtain the following system for $(a_1,\ldots,a_n,\lambda)$: $$\sum_{k=1}^n\frac{1}{a_k}=1;\quad\frac{k}{\left(\sum_{j=1}^{k}a_j\right)^2}=\frac{\lambda}{a_k^2}-\frac{\lambda}{a_{k+1}^2}\quad(1\leqslant k\leqslant n)$$ where we set $a_{n+1}:=\infty$ so that the very last term vanishes when $k=n$.

For fixed $\lambda$, the ratios $a_{k+1}/a_1$ can be computed from the second equation taken at $1\leqslant k\leqslant n-1$; then, the first equation determines $a_1$, and the last equation at $k=n$, finally, can be viewed as a single equation for $\lambda$. Observe also that if $(a_1,\ldots,a_n,\lambda)$ is the solution, then summation by parts yields simply $$S_n=\sum_{k=1}^{n}\left(\frac{\lambda}{a_k^2}-\frac{\lambda}{a_{k+1}^2}\right)\sum_{j=1}^{k}a_j=\lambda\sum_{k=1}^{n}\frac{1}{a_k}=\lambda.$$

This handles $\color{red}{n\sim 10^6}$ easily. For powers of $2$, running

foo(n,t)={my(b=1,r=0,s=0);for(k=1,n,if(b<0,return(-1));my(a=sqrt(b));r+=a;s+=1/a;b-=k/t/s^2);return(b)}; for(n=1,20,printf("%2u... ",n);print(n*log(2)*(2-solve(t=1.1,2.0,foo(2^n,t)))))

in PARI/GP (under different precision settings for sure), I get

1... 0.60650378298995214574[+++] 2... 1.02433225912743257353[+++] 3... 1.27241880529843009127[+++] 4... 1.39985861561679724754[+++] 5... 1.45256362011495626110[+++] 6... 1.46164482213152111129[+++] 7... 1.44578727489838766722[+++] 8... 1.41589954712815847210[+++] 9... 1.37839977433573904591[+++] 10... 1.33712471597920952819[+++] 11... 1.29439637045726263776[+++] 12... 1.25162537387072306196[+++] 13... 1.20966180786244720779[+++] 14... 1.16900506057647903574[+++] 15... 1.12993258591926390364[+++] 16... 1.09258058237638236418[+++] 17... 1.05699540545119896406[+++] 18... 1.02316675997623301901[+++] 19... 0.99104933243682303143[+++] 20... 0.96057697237757525984[+++]

Thus, the expression $L_n=(2-S_n)\ln n$ under the limit becomes $<1$. In fact $$L_{430605}=1.00000009015459070849\!+\\L_{430606}=0.99999998376961512393\!-$$