afvoer
Moderator: Rhiannon
- Berichten: 2.474
Re: afvoer
Als vergelijkingen kom ik op
$$g\frac{V(t)}{A}-\frac{32 \mu L v(t) }{d^2}=\frac{v^2(t)}{2}$$
$$\frac{dV(t)}{dt}=-\frac{v(t) \pi d^2}{4}$$
met \(V(0)=AH_1\).
Als je mijn eerste vergelijking in de tweede substitueert, dan krijg je een DV in v(t). Oplossen naar v(t) en terugsubstitueren in de eerste levert V(t) op en dan uit
$$V(t) = AH_1 - 1.5$$
de tijd t berekenen.
Relatief lange berekening en in mijn numerieke uitkomsten bij jouw oefeningen zitten toch altijd fouten .
$$g\frac{V(t)}{A}-\frac{32 \mu L v(t) }{d^2}=\frac{v^2(t)}{2}$$
$$\frac{dV(t)}{dt}=-\frac{v(t) \pi d^2}{4}$$
met \(V(0)=AH_1\).
Als je mijn eerste vergelijking in de tweede substitueert, dan krijg je een DV in v(t). Oplossen naar v(t) en terugsubstitueren in de eerste levert V(t) op en dan uit
$$V(t) = AH_1 - 1.5$$
de tijd t berekenen.
Relatief lange berekening en in mijn numerieke uitkomsten bij jouw oefeningen zitten toch altijd fouten .
- Berichten: 4.620
- Berichten: 2.474
Re: afvoer
In vorige post van mij moet het zijn
$$g\frac{V(t)}{A}-\frac{32 \mu L v(t) }{\rho d^2}=\frac{v^2(t)}{2}$$
$$g\frac{V(t)}{A}-\frac{32 \mu L v(t) }{\rho d^2}=\frac{v^2(t)}{2}$$
- Berichten: 2.474
Re: afvoer
Ik heb het eens numeriek geprobeerd ipv analytisch.
12862s kan dat? Is wel lang.
Code: Selecteer alles
import matplotlib.pyplot as plt
import numpy as np
from scipy.integrate import solve_ivp
from scipy.optimize import fsolve
from scipy.interpolate import UnivariateSpline
rho=1100;
g=9.81;
A=4;
L=20;
d=0.025;
H1=2;
mu=0.08;
def f(v):
return g*H1 - (32*mu*v*L / (rho*d**2)) - v**2/2
v0 = fsolve(f,1);
print('Beginsnelheid vloeistof = ', v0)
F = lambda t, v: -(np.pi*(d**2)*g/(4*A)) / (32*mu*L / (rho*v*(d**2)) + 1) ;
sol = solve_ivp(F, [0, 50000], v0, t_eval=np.linspace(0, 50000, 1000) )
time_array=sol.t
solution_array=sol.y
plt.xlabel('t')
plt.ylabel('v(t)')
plt.plot(time_array, solution_array[0,:])
plt.show()
volume_array = A/g*( solution_array**2 / 2 + 32*mu*solution_array*L / (rho*d**2))
plt.xlabel('t')
plt.ylabel('V(t)')
plt.plot(time_array, volume_array[0,:])
plt.show()
# calculate the difference array
difference_array = np.absolute(volume_array-6.5)
# find the index of minimum element from the array
index = difference_array.argmin()
print("Tijd is : ", time_array[index])
- Berichten: 4.620
Re: afvoer
@wnvl1 Dat ziet er goed uit
Bij turbulente stroming zal het een stuk sneller zijn verwacht ik.- Berichten: 2.474
Re: afvoer
Voor het turbulente geval kom ik op 1151s.
Code: Selecteer alles
import matplotlib.pyplot as plt
import numpy as np
from scipy.integrate import solve_ivp
from scipy.optimize import fsolve
from scipy.interpolate import UnivariateSpline
rho=1100;
g=9.81;
A=4;
L=20;
d=0.025;
H1=2;
mu=0.08;
f=0.005
F = lambda t, V: -(np.pi*(d**2))/4 * (g*V/A / (f*L/(2*d)+0.5))**0.5 ;
sol = solve_ivp(F, [0, 50000], [A*H1], t_eval=np.linspace(0, 50000, 1000) )
time_array=sol.t
solution_array=sol.y
plt.xlabel('t')
plt.ylabel('V(t)')
plt.plot(time_array, solution_array[0,:])
plt.show()
# calculate the difference array
difference_array = np.absolute(solution_array-6.5)
# find the index of minimum element from the array
index = difference_array.argmin()
print("Tijd is : ", time_array[index])
- Berichten: 10.584
Re: afvoer
wnvl1 gebruikt de Darcy frictiefactor, ukster de Fanning. Die verschillen een factor 4.
Cetero censeo Senseo non esse bibendum