Consider the Hessian matrix of $D(v)$:
$$
H=
\begin{bmatrix}
\frac{\partial^2 D(v)}{\partial v_1^2} & \frac{\partial^2 D(v)}{\partial v_1 \partial v_2}& \dots & \frac{\partial^2 D(v)}{\partial v_1 \partial v_k} \\
\frac{\partial^2 D(v)}{\partial v_2 \partial v_1}&\frac{\partial^2 D(v)}{\partial v_2^2} & \dots & \frac{\partial^2 D(v)}{\partial v_2 \partial v_k} \\
\vdots & \vdots & \ddots & \vdots\\
\frac{\partial^2 D(v)}{\partial v_k \partial v_1} & \frac{\partial^2 D(v)}{\partial v_k \partial v_2} & \dots & \frac{\partial^2 D(v)}{\partial v_k^2}
\end{bmatrix}
$$
Observe that the diagonal of the matrix consists of positive elements of the form:
$$\frac{ \partial^2 D(v)}{\partial v_i^2} = \frac{2 \left(k-1 \right) \prod^k_{j\neq i}v_i^2 \left( \sum_{j \in (1,k]} \prod_{p\neq j} v_p \right) }{\left( \sum^k_{j=1} \prod_{i\neq j}^k v_i \right)^3 }.$$
Off-diagonal elements are represented by negative mixed partial derivatives of the form:
$$ \frac{ \partial^2 D(v)}{\partial v_i \partial v_j} =-\frac{2 \left(k-1 \right) v_i^2 v_j^2 \prod_{p=1}^k v_i }{\left( \sum^k_{j=1} \prod_{i\neq j}^k v_i \right)^3 }. $$
Now, the first leading principal minor is:
$$ \Delta_1 = \frac{2 \left(k-1 \right) \prod_{j\neq 1}^k v_j^2 \left( \sum_{j \in (1,k]} \prod_{p\neq j} v_p \right) }{\left( \sum^k_{j=1} \prod_{i\neq j}^k v_i \right)^3 } > 0. $$
The second leading principal is:
$$ \Delta_2 =\Delta_1 \frac{2(k-1)\prod_{j\neq 2}^k v_j^2 \left( \sum_{j \in (2,k]} \prod_{p\neq j} v_p \right)}{ \left( \sum^k_{j=1} \prod_{i\neq j}^k v_j \right)^2 \left( \sum_{j \in (1,k]} \prod_{p\neq j} v_p \right) } >0 .$$
Assume that $m<k$ leading principal minor $\Delta_m>0$. Now I verify that $\Delta_{m+1}$ is also positive:
\begin{equation}
\label{m plus one}
\Delta_{m+1} = \Delta_m \frac{2(k-1)\prod_{j\neq m}^k v_j^2 \left( \overbrace{\sum_{j \in (m+1,k]} \prod_{p\neq j} v_p }^{n_1}\right)}{ \left( \sum^k_{j=1} \prod_{i\neq j}^k v_j \right)^2 \left( \sum_{j \in (m,k]} \prod_{p\neq j} v_p \right) }>0,
\end{equation}
which is always positive too.
Therefore, the first $m<k$ leading principal minors are positive by induction.
Also observe that term in $\Delta_{m+1}$ term $ n_1=0$ if $m+1=k$ and the last leading principal minor is $\Delta_k=0$.
Given that matrix, $H$ enjoys leading implies all (LIA) properties and all leading principal minors are positive, it immediately follows that $H$ is positive semi-definite and $D(v)$ is convex.
FunctionConvexityfunction on mathematica.stackexchange. user64494 suggested that it checks the non-negativity of Laplacian instead of checking the definition of the convexity directly. I am not familiar with a concept, but I am definitely going to read about it. – Oleh Dec 18 '21 at 19:31