En esta tesis estudiamos la solución exacta de un modelo de cosmología viscosa con métrica Friedman-Lemaître-Robertson-Wolker para describir la dinámica de la expansión tardía del Universo. A través del método de factorización, encontramos una solución analítica a la ecuación diferencial de Hubble, y restringimos los parámetros libres del modelo utilizando un método de estadística bayesiana que se basa en cadenas de Markov por Monte Carlo y toma como referencia las ultimas observaciones del parámetro de Hubble. Fue encontrado un valor para la constante de Hubble $H_0 = 67{.}51^{+0{.}41}_{-0{.}40}$, y un valor para el parámetro de desaceleración medido al día de hoy $q_0=-0{.}31^{+0{.}02}_{-0{.}02}$. Es mostrado que la solución encontrada conduce a una expansión acelerada del universo y presenta una transición suave de una fase desacelerada a una acelerada. Los valores numéricos obtenidos son consistentes con los reportados por la colaboración de Planck (Planck 2018 results. VI. Cosmological parameters).