You need two things. The angular velocity vector $\vec{\omega}$ at some time frame, and the linear velocity vector at some reference point $\vec{v}$.
The axis of rotation is located at a point relative to the reference point
$$ \vec{r} = \frac{\vec{\omega} \times \vec{v}}{ \| \vec{\omega} \|^2} \tag{1}$$
You can apply this at every instance.
See related post1 post2 and post3.
The proof is easy but it requires the use of the vector triple product identity $a \times (b \times c) = b (ac) - c (ab)$. Substitude above $\vec{v} = \vec{r} \times \vec{\omega}$ and carry out simplifications.
Simple example
The center of mass of an object is moving with $\vec{v} = \pmatrix{5 & -1 &0}$ and the rotational vector is $\vec{\omega} = \pmatrix{0 & 0 & 2}$ at some instant.
The point on the axis of rotation closest to the center of mass is located at $$\vec{r} = \frac{\pmatrix{0&0&2} \times \pmatrix{5&-1&0}}{ \| \pmatrix{0&0&2} \|^2 } = \frac{ \pmatrix{2 & 10 & 1} }{2^2} = \pmatrix{0.5 & 2.5 & 5}$$
To verify find the velocity at the center of mass with
$$ \vec{v} = \vec{\omega} \times (-\vec{r} ) = \vec{r} \times \vec{\omega} = \pmatrix{0.5 & 2.5 & 0} \times \pmatrix{0&0&2} = \pmatrix{5&-1&0} \;\checkmark $$
NOTE: That the axis of rotation is any point along the line, which includes $\vec{r}$ and all points parallel to $\vec{\omega}$.