afvoer

Moderator: Rhiannon

Reageer
Gebruikersavatar
Berichten: 4.534

afvoer

open vat
open vat.png
open vat.png (1.97 KiB) 1211 keer bekeken
ρ=1100kg/m3
g=9.81m/s2
A=4m2
L=20m
d=25mm
Beginhoogte H1=2m

Wat is de benodigde tijd voor afvoer van 1,5m3 door de pijp bij:
1. Laminaire stroming μ=80mPas
2. Turbulente stroming f=0,005

Gebruikersavatar
Berichten: 2.312

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 :( .

Gebruikersavatar
Berichten: 4.534

Re: afvoer

wnvl1 schreef: zo 18 feb 2024, 15:58 Relatief lange berekeningen
Valt reuze mee...
Hagen- Poiseuiille substitueren in de Continuiteitsvergelijking
Scheiden van variabelen
Integreren

Gebruikersavatar
Berichten: 2.312

Re: afvoer

We zullen proberen.

Gebruikersavatar
Berichten: 2.312

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}$$

Gebruikersavatar
Berichten: 2.312

Re: afvoer

Ik heb het eens numeriek geprobeerd ipv analytisch.

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])
Uitlaatvolume.png
Uitlaatvolume.png (9.57 KiB) 945 keer bekeken
12862s kan dat? Is wel lang.

Gebruikersavatar
Berichten: 4.534

Re: afvoer

@wnvl1 Dat ziet er goed uit :)
1.png
2.png
2.png (6.89 KiB) 920 keer bekeken
3.png
3.png (10.43 KiB) 920 keer bekeken
Bij turbulente stroming zal het een stuk sneller zijn verwacht ik.

Gebruikersavatar
Berichten: 2.312

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])
Uitlaatvolumeturbulent.png
Uitlaatvolumeturbulent.png (10.25 KiB) 887 keer bekeken

Gebruikersavatar
Berichten: 4.534

Re: afvoer

Met deze formule kom ik op praktisch een factor 2 hoger.
turbulente stroming.png
turbulente stroming.png (3.15 KiB) 836 keer bekeken

Gebruikersavatar
Berichten: 2.312

Re: afvoer

Ik zal eerstdaags eens naar mijn fout zoeken.

Gebruikersavatar
Berichten: 10.560

Re: afvoer

wnvl1 gebruikt de Darcy frictiefactor, ukster de Fanning. Die verschillen een factor 4.
Cetero censeo Senseo non esse bibendum

Gebruikersavatar
Berichten: 2.312

Re: afvoer

Goed gezien, dan hebben we allebei gelijk.

Reageer