The real Stiefel manifold $V\_{n,k}$ of orthogonal $k$-frames in $\mathbb{R}^n$ can be viewed as the reductive homogeneous space $G/H=O(n)/O(n-k)$. If ${\frak{so}}(n)$ is the Lie algebra of $O(n)$, then we have the reductive decomposition
$$
{\frak{so}}(n)={\frak{m}}+{\frak{h}}
$$
where
$$
{\frak{m}}=\left \{
\begin{pmatrix}
 A & B\\
 -B^T & O
\end{pmatrix}: A\>\> \text{is a}\>\> k\times k \>\>\text{skew-symmetric matrix} \right \}
$$
and 
$$
{\frak{h}}=\left \{
\begin{pmatrix}
 O & O\\
 O & C
\end{pmatrix}: C\>\> \text{is a}\>\> (n-k)\times (n-k) \>\>\text{skew-symmetric matrix} \right \}
$$
with ${\frak{so}}(n)$ and ${\frak{h}}$ a reductive pair.


According to Helgason, there is a local diffeomorphism 
$$
(\text{exp}X,h)\mapsto (\text{exp}X)h,\quad \text{where}\>\> X\in {\frak{m}}, h\in H. 
$$
However, there doesn't seem to be anything in the literature where this decomposition is calculated explicitly. I did find the following ``polar decomposition":
$$
V\_{n,k}\times P\_k\rightarrow M\_{n,k},\quad (v,r)\mapsto vr^{1/2}
$$
where $P\_k$ is the set of positive semi-definite symmetric matrices and $M\_{n,k}$ is the set of all $n\times k$ real matrices.


Does this last decomposition have any connection with the polar decompositions of Lie groups? If not, is there an explicit description of the polar decomposition involving Stiefel manifolds somewhere?



Stiefel manifolds and polar decompositions
It is easier for me to examine this geometrically, rather than from the point of view of Lie groups and algebras.


First, the identity matrix $I \in O(n)$ represents the standard orthonormal basis $e\_1, \dots, e\_n$ of $\mathbb{R}^n$, and its coset $I\cdot O(n-k)$ represents the $k$-plane spanned by the first $k$ basis vectors $e\_1, \dots, e\_k$.


The local splitting you're looking for is equivalent to extending the standard orthonormal basis to a local map from $V\_{n,i}$ to $O(n)$, viewed as the space of orthonormal frames in $\mathbb{R}^n$, such that the span of the first $k$ vectors in each frame is equal to the corresponding $k$-plane in $V\_{n,k}$. Helgason is basically pointing out that one way to do this is to use the connection naturally induced by the bi-invariant Riemannian metric on $O(n)$ to parallel transport the initial frame along geodesic rays starting at the identity matrix.


This local section can be described as a solution to a linear second first order system of ODE's along each geodesic (essentially the Jacobi equations along each geodesic), and in principle these ODE's could be integrated to give some kind of formula for this section. Unfortunately, I've never computed this and don't know what it looks like. And it's not clear to me how to transfer this description to the setting of Lie groups and algebras. My wild guess is that you can write the system of ODE's in terms of the Lie algebra, but you won't be able to integrate them explicitly.


One possibly easier way to do this is to using the dual description using invariant differential forms as presented in:


Griffiths, P.
On Cartan's method of Lie groups and moving frames as applied to uniqueness and existence questions in differential geometry.
Duke Math. J. 41 (1974), 775–814.


