Mediante una descomposición ortogonal adecuada (POD) para reducir la dimensionalidad de los vectores de coeficientes de la solución de elementos finitos (FE) desconocidos en el método FE mixto Crank-Nicolson (CN) (CNMFE) para ecuaciones de Stokes no estacionarias bidimensionales (2D) en relación con las funciones de corriente de vorticidad, se construye un método reducido de CNMFE recursivo (RDR-CNMFE). En este caso, el método RDR-CNMFE tiene las mismas funciones de base de FE y precisión que el método CNMFE. La existencia, estabilidad y errores de las soluciones de RDR-CNMFE se analizan mediante análisis matricial, lo que resulta en una teoría de análisis muy simple. Se realizan algunos intentos numéricos para verificar la validez del método RDR-CNMFE. El método RDR-CNMFE tiene una precisión temporal de segundo orden y pocos desconocidos para poder acortar el tiempo de CPU y retardar la acumulación de errores en el proceso de cálculo de simulación, y mejorar la precisión de cálculo en tiempo real.