Estuve haciendo una pequeña simulación y la fuerza de rozamiento obtenida por la expresión de Stokes sale como 4 órdenes de magnitud menos, por lo que prácticamente no afecta a la aceleración de bajada. Posiblemente el error esté en las unidades de algún dato.
Aquí mi simulación:
%Datos líquido
de=1000;
vis=1;
%Datos esfera
r=10;
densidadcobre=8500;
m=
(densidadcobre*
(4/
3)*
pi*
(r.^
3));
%Datos iniciales: altura inicial, velocidad inicial, aceleración y altura
%del líquido
yo=15;
vyo=0;
g=-9.8;
yent=5;
%%Ecuaciones caida libre
tent=
(-vyo-
sqrt(vyo^
2-
2*g*
(yo-yent
)))/g;
t=0:0.001:tent;
vy = vyo + g*t;
y=yo+vyo*t+(1/2)*g*t.^2;
while y1>0
al=(g*m+b)/m;
v1= vyo + al*
(t
(length(t
))+
0.01);
y1=yo+vyo*
(t
(length(t
))+
0.01)+
(1/
2)*al*
(t
(length(t
))+
0.01)^
2;
y=[y,y1];
vy=[vy,v1];
end
%%Ecuaciones Fluidos
ylabel('Velocidad en Y - Vy (m/s)')
Eso está muy poco optimizado y es un primer ejemplo solo para el caso en el que se tira desde una altitud superior a la entrada en el líquido.