Due to the fact that your multivariate-normal distribution is independent (the covariance matrix is diagonal), it will be more intuitive to treat each dimension as being its 'own' random variable; i.e., let $(I, J, K) \sim (\log(X), \log(Y), \log(Z))$, where $X, Y, Z$ are all normal distributions with their own mean and variance parameters.
Now, due to the independence, the joint distribution is given by $f(i, j, k) = f_I(i)f_J(j)f_K(k)$, where $f_A(a)$ is the marginal density function of RV $A$. If you want to 'fix', aka condition on knowing $k$, then you want the distribution $I, J|K$. The density for this can be found with Bayes rule:
$$f(i, j|k) = \frac{f(i, j, k)}{f_K(k)}\;;$$
but given we can factorise $f(i, j, k)$ using the above independence assumption we now have
$$f(i, j|k) = \frac{f(i, j, k)}{f_K(k)} = \frac{f_I(i)f_J(j)f_K(k)}{f_K(k)} = f_I(i)f_J(j)\;.$$
So, you can see that the distribution of $I, J|K$ has the same distribution as the product of the two random variables $I$ and $J$. Putting this together, we can deduce that the distribution is still log-normally distributed (indeed, it is the product of two independent log-normal distributions) and they will still be independent.