UNIVERSIDAD PRIVADA DEL VALLE

FACULTAD DE TECNOLOGIA
INGENIERIA CIVIL
CAMPUS LA PAZ

FISICA II

INFORME N.º 7
Viscosidad de Stokes

Grupo “C”

Estudiante: Facundo Ismael Soto

Guerrero

Docente: Lic. Valero Choque

La Paz 08 de mayo de 2026

Gestión II –2026
1. OBJETIVOS
La presente práctica de laboratorio tiene los siguientes objetivos:
• Determinar experimentalmente el coeficiente de viscosidad dinámico (η) de la glicerina
mediante el método de caída de esferas de Stokes.
• Verificar la proporcionalidad entre la velocidad terminal de caída y el cuadrado del radio de
las esferas, establecida por la Ley de Stokes.
• Aplicar métodos estadísticos de regresión lineal por mínimos cuadrados y t de Student para
cuantificar la incertidumbre experimental.
• Comparar el valor experimental de la viscosidad con el valor de referencia reportado en
tablas para la glicerina.

2. MARCO TEÓRICO
La viscosidad es una propiedad intrínseca de los fluidos que mide su resistencia al flujo o a la
deformación por esfuerzo cortante. Un fluido de alta viscosidad fluye lentamente, mientras que uno
de baja viscosidad fluye con facilidad. Esta propiedad es fundamental en múltiples aplicaciones de
la ingeniería, como el diseño de tuberías, la lubricación de maquinaria y la caracterización de
polímeros.

2.1 Ley de Stokes

Cuando una esfera sólida cae a través de un fluido viscoso en régimen laminar (número de
Reynolds Re << 1), la fuerza de arrastre viscoso que experimenta, conocida como fuerza de
Stokes, está dada por:
F_arrastre = 6π η r v
donde η es el coeficiente de viscosidad dinámico del fluido (Pa·s), r es el radio de la esfera (m) y v
es la velocidad de la esfera respecto al fluido (m/s).

2.2 Velocidad terminal

Una esfera que cae en un fluido viscoso experimenta tres fuerzas: su peso W, el empuje
hidrostático E (principio de Arquímedes) y la fuerza de arrastre viscoso F_d. Cuando la aceleración
neta es cero, la esfera alcanza la velocidad terminal v_t. El equilibrio dinámico de fuerzas impone:
W = E + F_d
ρ_s · (4/3)πr³ · g = ρ_f · (4/3)πr³ · g + 6π η r v_t
Despejando el coeficiente de viscosidad η, se obtiene la ecuación de trabajo del experimento:
η = [2r²(ρ_s - ρ_f) · g] / (9 · v_t) ... (1)
donde ρ_s es la densidad de la esfera (kg/m³), ρ_f es la densidad del fluido (kg/m³), g es la
aceleración gravitacional (9,81 m/s² para La Paz) y v_t es la velocidad terminal medida
experimentalmente.

2.3 Linealización para mínimos cuadrados

De la ecuación (1) se puede expresar la velocidad terminal como función del cuadrado del radio:
v_t = [2(ρ_s - ρ_f) · g / (9η)] · r² ... (2)
Esta relación es lineal de la forma v_t = m·r² + b, donde la pendiente m = 2(ρ_s - ρ_f)g / (9η)
permite obtener η una vez determinada m por regresión lineal.

2.4 Densidad de la esfera

La densidad de cada esfera metálica se calcula a partir de su masa medida con balanza y su
volumen geométrico determinado con el tornillo micrométrico:
ρ_s = m / [(4/3)πr³] ... (3)

3. MATERIALES Y PROCEDIMIENTO
3.1 Materiales y equipos

N° Instrumento / Material Resolución / Incertidumbre

Apreciación instrumental δ
1 Cronómetro digital 0,01 s δt = 0,01 s
2 Regla graduada (escala mm) 1,0 mm δL = 0,5 mm
3 Balanza analítica 0,01 g δm = 0,005 g
4 Tornillo micrométrico 0,001 mm δd = 0,001 mm
5 Probeta graduada con glicerina (≥ 30 cm) — —
6 5 esferas metálicas de diferentes — —
diámetros
7 Pinzas, hilo guía, cronómetro — —
3.2 Procedimiento experimental
1. Se preparó la probeta llenándola con glicerina pura hasta una altura mínima de 30 cm,
marcando con cinta adhesiva las posiciones de inicio y fin de la zona de medición (L =
200,0 mm), asegurando que la esfera ya hubiera alcanzado la velocidad terminal antes de
cruzar la marca de inicio.
2. Con el tornillo micrométrico se midió el diámetro de cada esfera en cuatro posiciones
distintas (dos posiciones angulares perpendiculares). Se registraron los valores d₁, d₂, d₃ y
d₄ para calcular el diámetro promedio.
3. Con la balanza analítica se midió la masa de cada esfera m_i.
4. Se soltó cada esfera desde el reposo por encima de la marca superior. Con el cronómetro
se midió el tiempo de caída t entre las dos marcas. Se realizaron 5 repeticiones por esfera
para reducir el error aleatorio.
5. Se calculó la velocidad terminal v_t = L / t_promedio y posteriormente el coeficiente de
viscosidad η para cada esfera usando la ecuación (1).
6. Se graficó v_t en función de r² y se aplicó regresión lineal por mínimos cuadrados para
obtener la pendiente m y calcular η de forma independiente.

4. DATOS EXPERIMENTALES
Distancia de caída medida con regla graduada: L = 200,0 mm = 0,2000 m (δL = 0,5 mm)
Densidad de la glicerina (referencia a 20 °C): ρ_f = 1 260,0 kg/m³
Aceleración gravitacional local (La Paz, Bolivia): g = 9,81 m/s²

4.1 Diámetros de las esferas (tornillo micrométrico, δd = 0,001 mm)

Esfera d₁ (mm) d₂ (mm) d₃ (mm) d₄ (mm) d̄ (mm) r̄ (mm)
1 3.954 3.952 3.956 3.953 3.954 1.977
2 5.012 5.015 5.010 5.013 5.012 2.506
3 6.351 6.348 6.352 6.350 6.350 3.175
4 7.894 7.891 7.896 7.893 7.894 3.947
5 9.125 9.128 9.122 9.126 9.125 4.563

4.2 Masas de las esferas (balanza analítica, δm = 0,005 g)

Esfera m (g) m (kg) ρ_s (kg/m³) Material aprox.
1 0.25 0.00025 7725.3 Acero / hierro
2 0.52 0.00052 7885.7 Acero / hierro
3 1.05 0.00105 7831.0 Acero / hierro
4 2.01 0.00201 7805.3 Acero / hierro
5 3.10 0.00310 7791.6 Acero / hierro

4.3 Tiempos de caída (cronómetro, δt = 0,01 s, L = 200,0 mm)

Esfera t₁ (s) t₂ (s) t₃ (s) t₄ (s) t₅ (s) t̄ (s) v_t
(m/s)
1 4.82 4.79 4.85 4.81 4.83 4.820 0.0415
2 2.38 2.41 2.35 2.39 2.37 2.380 0.0840
3 1.19 1.22 1.18 1.21 1.20 1.200 0.1667
4 0.62 0.60 0.63 0.61 0.62 0.616 0.3247
5 0.36 0.38 0.35 0.37 0.36 0.364 0.5495
5. CÁLCULOS Y ANÁLISIS DE ERRORES
5.1 Viscosidad η por esfera (ecuación de Stokes)
Aplicando la ecuación (1) para cada esfera: η_i = 2·r²·(ρ_s − ρ_f)·g / (9·v_t)
Ejemplo de cálculo detallado — Esfera 1:
η₁ = 2 × (1,977×10⁻³)² × (7725.3 − 1260,0) × 9,81 / (9 × 0.0415)
η₁ = 2 × 3,908×10⁻⁶ × 6465.3 × 9,81 / 0.3734
η₁ = 1.3275 Pa·s

Esfera r (×10⁻³ m) ρ_s (kg/m³) Δρ (kg/m³) v_t (m/s) η_i (Pa·s) η_i
(cP)
1 1.977 7725.3 6465.3 0.0415 1.3275 1327.5
2 2.506 7885.7 6625.7 0.0840 1.0797 1079.7
3 3.175 7831.0 6571.0 0.1667 0.8665 866.5
4 3.947 7805.3 6545.3 0.3247 0.6846 684.6
5 4.563 7791.6 6531.6 0.5495 0.5395 539.5

5.2 Error estadístico por t de Student (gl = 4, 95 %)

Con n = 5 mediciones independientes de η (una por esfera), se calcula la desviación estándar
muestral s y el error de la media mediante el estadístico t de Student con gl = n − 1 = 4:
s_η = √[Σ(η_i − η̄)² / (n−1)]
s_η = 0.3132 Pa·s
Δη = t(gl=4, 95%) · s_η / √n = 2,776 × s_η / √5
Δη = 2,776 × 0.3132 / √5 = 0.3889 Pa·s
η = (0.8995 ± 0.3889) Pa·s

5.3 Regresión lineal — Mínimos cuadrados (v_t vs r²)

Se linealiza la Ley de Stokes graficando v_t en función de r². La tabla de sumas auxiliares es:

i r² (m²) v_t (m/s) (r²)² (m⁴) (v_t)² (m²/s²) r²·v_t (m³/s)

1 3.9080e-6 0.0415 1.5273e-11 0.00172 1.6216e-7
2 6.2813e-6 0.0840 3.9455e-11 0.00706 5.2784e-7
3 1.0081e-5 0.1667 1.0164e-10 0.02778 1.6802e-6
4 1.5577e-5 0.3247 2.4264e-10 0.10541 5.0574e-6
5 2.0818e-5 0.5495 4.3337e-10 0.30190 1.1438e-5
Σ 5.6665e-5 1.1663 8.3237e-10 0.44387 1.8866e-5

Cálculo de la pendiente m y el intercepto b:

m = (n·Σx·y − Σx·Σy) / (n·Σx² − (Σx)²) = 2.9697e+4 m·s⁻¹/m²
b = (Σy − m·Σx) / n = -0.1033 m/s
r = (n·Σxy − Σx·Σy) / √[(n·Σx² − (Σx)²)(n·Σy² − (Σy)²)] = 0.98805
Incertidumbre en pendiente e intercepto (usando desviación estándar de residuos S_r):
S_r = √[Σ(y_i − ŷ_i)² / (n−2)] = 3.6889e-2 m/s
Δm = S_r · √(n/Δ) = 2.6749e+3 m·s⁻¹/m²
Δb = S_r · √(Σx²/Δ) = 3.4513e-2 m/s

Determinación de η a partir de la pendiente (usando ρ_s promedio = 7807.8 kg/m³):

m = 2(ρ_s − ρ_f)g / (9η) → η = 2(ρ_s − ρ_f)g / (9m)
η = 2 × (7807.8 − 1260,0) × 9,81 / (9 × 2.9697e+4)
η (pendiente) = 0.4807 Pa·s

5.4 Propagación de errores instrumentales

La incertidumbre instrumental en la velocidad terminal v_t = L/t se propaga por cuadratura:
δv_t / v_t = √[(δL/L)² + (δt/t)²]
Ejemplo para la esfera 1 (t̄ = 4.820 s, L = 0,2000 m):
δv_t = v_t · √[(0,0005/0,2000)² + (0,01/4.820)²] = 0.00013 m/s

La propagación para η a partir de r y v_t:

δη / η = √[(2·δr/r)² + (δv_t/v_t)²]
δη₁ = η₁ · √[(2×5×10⁻⁷/1.9769e-3)² + (0.00013/0.0415)²] = 0.0044 Pa·s

6. RESULTADOS
Parámetro Valor obtenido
η (t de Student, n=5) (0.8995 ± 0.3889) Pa·s = (899.5 ± 388.9) cP
η (mínimos cuadrados / pendiente) 0.4807 Pa·s = 480.7 cP
Coeficiente de correlación r 0.98805
η referencia glicerina 20°C ≈ 1,490 Pa·s = 1490 cP
Error relativo (t de Student) |η_exp − η_ref| / η_ref × 100% = 39.6 %

7. GRÁFICAS
7.1 Velocidad terminal v_t en función de r² (Ley de Stokes linealizada)
La gráfica muestra la relación lineal entre la velocidad terminal v_t y el cuadrado del radio r²,
consistente con la Ley de Stokes. La recta de regresión se ajusta con una correlación r = 0.9880,
indicando una dependencia lineal significativa.
 Figura 1. v_t vs r² | y = 2.970e+4·x + (-0.1033) | r = 0.9880

7.2 Viscosidad calculada por esfera

La siguiente gráfica muestra los valores de η calculados individualmente para cada esfera
mediante la ecuación de Stokes. La línea punteada roja representa el valor promedio.

Figura 2. Viscosidad η calculada por esfera | η̄ = 0.8995 Pa·s

8. ANÁLISIS DE RESULTADOS
Los resultados experimentales muestran que la velocidad terminal de caída de las esferas aumenta
con el radio, en concordancia con la predicción de la Ley de Stokes. El coeficiente de correlación r
= 0.98805 en la gráfica v_t vs r² confirma que la relación lineal se mantiene con buen grado de
ajuste, lo que valida la hipótesis de flujo laminar alrededor de las esferas.
El coeficiente de viscosidad promedio obtenido mediante t de Student es η = (0.8995 ± 0.3889)
Pa·s, equivalente a (899.5 ± 388.9) cP. Comparado con el valor de referencia para la glicerina a 20
°C (≈ 1 490 cP), se obtiene un error relativo del 39.6 %.
La discrepancia entre el valor experimental y el de referencia puede atribuirse a: (1) la temperatura
del fluido durante la práctica, ya que la viscosidad de la glicerina es extremadamente sensible a la
temperatura (varía en un orden de magnitud entre 20 °C y 30 °C); (2) posibles impurezas en la
glicerina que modifican su densidad y viscosidad; (3) efectos de pared de la probeta que aumentan
el arrastre efectivo cuando la relación d_esfera/d_probeta no es despreciable; (4) la dificultad de
asegurar que la esfera haya alcanzado exactamente la velocidad terminal antes de cruzar la marca
de inicio.
Los valores de η obtenidos muestran una tendencia decreciente con el radio de la esfera (η₁ > η₂
> … > η₅), lo que sugiere que las esferas más grandes podrían no estar completamente en
régimen de Stokes (Re << 1). Para las esferas de mayor diámetro, el número de Reynolds puede
ser suficientemente alto para que aparezcan efectos no lineales en la fuerza de arrastre,
subestimando la resistencia real y por tanto el η calculado.

9. CONCLUSIONES
• Se determinó experimentalmente el coeficiente de viscosidad dinámico de la glicerina
mediante el método de Stokes, obteniendo η = (0.8995 ± 0.3889) Pa·s, con un error relativo
del 39.6 % respecto al valor de referencia.
• La gráfica v_t vs r² exhibe un comportamiento lineal con r = 0.9880, verificando
experimentalmente la proporcionalidad establecida por la Ley de Stokes entre la velocidad
terminal y el cuadrado del radio.
• La pendiente de la recta de regresión permitió obtener una segunda estimación de la
viscosidad (η = 0.4807 Pa·s), independiente del promedio directo, proporcionando un
método de verificación cruzada.
• Las principales fuentes de incertidumbre identificadas son la temperatura de la glicerina, los
efectos de pared de la probeta y la posible salida del régimen laminar para las esferas de
mayor diámetro.
• Los métodos de análisis de errores empleados —t de Student para el promedio de η y
mínimos cuadrados para la pendiente— son adecuados para cuantificar la incertidumbre
experimental en este tipo de práctica.

10. CUESTIONARIO
1. ¿Qué condiciones deben cumplirse para que la Ley de Stokes sea válida?
La Ley de Stokes es válida bajo las siguientes condiciones: (a) Flujo laminar en torno a la esfera,
es decir, número de Reynolds Re = ρ_f·v·d/η << 1 (generalmente Re < 0,5). (b) El fluido debe ser
incompresible, newtoniano y de extensión infinita (o al menos mucho mayor que la esfera). (c) La
esfera debe ser rígida, lisa y esférica. (d) La velocidad debe ser la terminal, es decir, la aceleración
neta debe ser cero. En la práctica, se verifica aproximadamente el cumplimiento de Re << 1
calculando el número de Reynolds con los datos obtenidos.
2. ¿Por qué la viscosidad de la glicerina depende fuertemente de la temperatura?
La viscosidad de los fluidos líquidos disminuye al aumentar la temperatura porque la energía
térmica que adquieren las moléculas les permite superar más fácilmente las fuerzas
intermoleculares de cohesión que se oponen al deslizamiento relativo entre capas. En la glicerina,
que es un poliol de alta cohesión molecular (forma puentes de hidrógeno extensos), esta
dependencia es especialmente pronunciada: su viscosidad puede pasar de ≈ 12 000 cP a 0 °C a
solo ≈ 200 cP a 40 °C. Por esta razón, es imprescindible registrar y reportar la temperatura del
fluido durante el experimento.

3. ¿Cómo afecta el efecto de pared de la probeta a los resultados?

Cuando el diámetro de la esfera (d) no es despreciable respecto al diámetro interno de la probeta
(D), la velocidad terminal medida es menor que la velocidad en un fluido infinito debido al aumento
del gradiente de velocidad cerca de las paredes. Este efecto se corrige mediante el factor de
Faxén: v_t_corregida = v_t_medida / (1 − 2,104·(d/D) + 2,09·(d/D)³ − …). En la práctica, para
minimizar el error, se recomienda que d/D < 0,1. Si no se corrige este efecto, el coeficiente de
viscosidad calculado resulta sobreestimado.

4. ¿Qué es el número de Reynolds y cómo se verifica el régimen laminar en este

experimento?
El número de Reynolds es un parámetro adimensional que representa la razón entre las fuerzas
inerciales y las fuerzas viscosas en un fluido: Re = ρ_f·v·d/η. Para la esfera 1: Re₁ = 1260 × 0.0415
× 0.00395 / 0.8995 ≈ 0.230. Un Re << 1 confirma flujo laminar (régimen de Stokes). Para las
esferas de mayor diámetro y velocidad, Re puede alcanzar valores mayores, lo que justifica la
desviación observada en los valores de η.

5. ¿Cuál es la diferencia entre viscosidad dinámica y viscosidad cinemática?

La viscosidad dinámica η (Pa·s = kg/(m·s)) mide la resistencia interna al flujo de un fluido bajo la
acción de un esfuerzo cortante, según la Ley de Newton de la viscosidad: τ = η·(dv/dy). La
viscosidad cinemática ν (m²/s) relaciona la viscosidad dinámica con la densidad del fluido: ν =
η/ρ_f. La viscosidad cinemática es la que aparece en el número de Reynolds y es relevante cuando
el flujo está gobernado tanto por la inercia como por la viscosidad. En la Ley de Stokes se trabaja
directamente con η (dinámica), ya que la fuerza de arrastre es proporcional a η y no a ν.