En este trabajo se presenta la solución numérica-analítica para el fujo de un ferrofuido en el espacio anular entre  dos cilindros concéntricos e inducido por un campo magnético rotativo. A diferencia de soluciones analíticas  presentadas en la literatura, este análisis tiene en cuenta el efecto de los términos ω×M y en la ecuación  de magnetización los cuales son comúnmente despreciados con el fn de desacoplar el problema magnético del  hidrodinámico y entonces poder obtener una solución analítica. No obstante, se encontró que su efecto es despreciable  bajo la suposición de que el vector de magnetización del ferrofuido es proporcional al vector del campo magnético. Los resultados numéricos muestran un buen ajuste con la solución asintótica reportadas por Chaves et al. 2010.  También se revisa el rango de aplicación de este análisis y otros reportados en la literatura bajo la suposición de  proporcionalidad de los vectores de magnetización y campo magnético.    

Este artículo describe una estrategia de discretización de la ecuación de difusión en mallas no estructuradas aplicando el método de los volúmenes finitos con las variables calculadas en el centroide de cada volumen. Esta aproximación está basada en el trabajo de Date [1] que usa una técnica iterativa conocida como corrección diferida para solucionar el cálculo del flujo difusivo en mallas no ortogonales. Se comprobó que para ángulos internos del elemento menores a 50°, el método propuesto por Date no converge y entonces se propone una nueva forma de calcular el gradiente que favorece la convergencia del problema. Se muestra un estudio de la convergencia donde se demuestra la alta efectividad del método propuesto. A partir de la solución de un problema típico, basado en la solución de la ecuación de Poisson, se realizó la comparación de los resultados obtenidos con los de la solución analítica, donde se observó una alta correspondencia de los resultados sin comprometer el tiempo de cálculo. Finalmente, se demostró la flexibilidad de la aproximación implementada al realizar simulaciones sobre mallas estructuradas y no estructuradas, usando elementos en forma de cuadriláteros y triángulos en 2D, y cubos curvilíneos y tetraedros en 3D.