De hevel Rik Schepens 0772841
Rob Wu 0787817
23 maart 2012
Begeleider: Arris Tijsseling
Modelleren A Vakcode: 2WH01
Inhoudsopgave Samenvatting
1
1 Inleiding
1
2 Theorie
2
3 Model
3
4 Resultaten en conclusie
4
5 Discussie
4
A Berekeningen A.1 Verwerking formules . . . . . . . . . . . . . . . . . . . . . . . A.2 Een aantal uitkomsten . . . . . . . . . . . . . . . . . . . . . .
5 5 7
B Programma
8
Referenties
11
Samenvatting Een hevel is een (niet noodzakelijk continue dalende) buis die wordt gebruikt om een vloeistof van een reservoir naar een lager gelegen reservoir. Onder invloed van de luchtdruk vloeit het water omhoog, waarna het zakt vanwege de zwaartekracht. De stroomsnelheid is afhankelijk van de druk, de gravitatiekracht, het hoogteverschil en de gebruikte vloeistof. Deze kan worden gevonden met behulp van de wet van Bernoulli, en de vergelijking van Darcy-Weisbach, beide afkomstig uit de vloeistofdynamica. In de praktijk is de weerstand van de hevel te verwaarlozen. Wanneer de verhouding tussen de lengte en het kwadraat van de diameter groter is dan 1000, is de wrijving niet meer te verwaarlozen. De maximale hoogte waarover geheveld kan worden is 10 meter. De maximale lengte is afhankelijk van de diameter van de buis. Zolang de verhoulengte ding diameter 2 lager is als 200.000, is de hevel nog werkzaam, met een redelijke stroomsnelheid.
1
Inleiding
Een hevel is een buis of slang die wordt gebruikt om te hevelen. Hevelen is het overbrengen van vloeistof van het ene reservoir naar een ander lager reservoir. Het bijzondere hieraan is dat dit niet hoeft te gebeuren via een continu dalende slang/buis. Dankzij natuurkundige principes kan deze buis eerst een stuk omhoog lopen en daarna pas naar beneden. Dit is handig voor een groot aantal toepassingen in het dagelijks leven. Van het legen van een zwembad tot het stelen van benzine uit de tank van een auto. Dit is een verslag van het onderzoek naar de eigenschappen (waaronder de stroomsnelheid en invloed van weerstand) en fysieke grenzen van dit principe. Hiervoor worden er eerst twee belangrijke formules uit de vloeistofdynamica gepresenteerd, die in het model verwerkt worden. Dit model is ge¨ımplementeerd in een programma, waarmee de gezochte resultaten worden berekend.
1
Figuur 1: Een hevel
2
Theorie
Wet van Bernoulli Volgens de Wet van Bernoulli een relatie tussen de stroomsnelheid van de vloeistof, de valversnelling, het hoogteverschil, de druk en de dichtheid. Deze wet wordt beschreven door de volgende formule: v2 P + gh + = constant 2 ρ
(1)
Waarbij, voor een gegeven punt geldt: v is de stroomsnelheid van de vloeistof in het punt (m/s) g is de valversnelling (m · s−2 ) h is het hoogteverschil (m) P is de druk in het punt (Pa = kg · m · s−2 ) ρ is de dichtheid van de vloeistof (kg · m−3 )
Darcy-Weisbach vergelijking Om de invloed van wrijving te berekenen, wordt de Darcy-Weisbach vergelijking gebruikt. Het resultaat van deze formule is een indicatie van het energieverlies door wrijving. Dit verlies is in de formule van Bernoulli terug te zien in de vorm van een verlaagde druk: hf = f ·
L V2 · , D 2g
(2)
Hierbij geldt voor een gegeven punt: hf is de stijghoogte ten gevolge van wrijving (m) f is de Darcy-wrijvingsco¨effici¨ent (dimensieloos), gegeven door 64 f = Re , waarbij Re het getal van Reynolds is, gegeven door Re = ρvD µ ρ is de dichtheid van de vloeistof (kg · m−3 ) µ is de dynamische viscositeit van de vloeistof (Pa· s) L is de lengte van de buis (m) D is de hydraulische diameter van de buis (m). In het geval van een buis is dit precies gelijk aan de binnenste diameter. V is de gemiddelde stroomsnelheid van de vloeistof (m · s−1 ) (Merk op, omdat de hevel een gesloten circuit is, is deze gelijk aan de snelheid v) 2
g is de valversnelling (m · s−2 ) Met de stijghoogte kan vervolgens het drukverlies worden berekend: ∆p = ρ · g · hf
(3)
∆p is het drukverlies (kg · m · s−2 ) ρ is de dichtheid van de vloeistof (kg · m−3 ) g is de valversnelling (m · s−2 ) hf is de stijghoogte (m)
3
Model
De hevel is verwerkt in een eenvoudig model (zie het programma A.2). Aan de hand van een beschrijving van de hevel (zie fig 3) berekent dit programma de gevraagde onbekenden. Hiervoor zijn eerst de formules van Bernoulli (vergelijking 1) en Darcy-Weisbach (vergelijking 2) omgeschreven, om op handige wijze de waarden van onder andere de maximale hoogte te vinden. De uitwerking hiervan is te vinden in bijlage (A). Zo kan dit model voor elke beschreven hevel de mondingssnelheid uitrekenen. Deze wordt gevonden door de volgende berekening uit te voeren: v=
p 2 · g · hc
(4)
Echter omdat de hevel een aaneengesloten systeem is, is deze snelheid begrensd. Het stuk dat het water omhoog moet lopen kan namelijk niet onbeperkt snel. De maximale snelheid wordt vervolgens gegeven door: s P − g · hb (5) vmax = 2 ρ
3
Figuur 2: Illustratie van het hevelmechanisme (Bron: [?] In het voorgaande is de invloed van wrijving niet vermeld, aangezien deze in de praktijk vaak verwaarloosbaar is. Wanneer er echter wordt gekeken naar buitengewone gevallen, zoals hevelen over een grote afstand door een zeer dunne buis, dan is de invloed van de wrijving niet meer verwaarloosbaar. Als de DL2 verhouding 1000 of hoger is, kan de wrijving niet meer worden verwaarloosd. In dit geval moet het drukverlies ∆p worden berekend (formule 2 en 3), zodat de waarde van de druk P in formule 1 kan worden aangepast.
4
Resultaten en conclusie
Na verwerking van de gegevens (zie A) volgt dat de maximale hoogte waarover water kan worden geheveld ongeveer 10 meter is. Ook blijkt dat de maximale lengte van de buis heel groot kan zijn, zolang de diameter ten minste groot genoeg is. Anders gezegd, zolang de DL2 -verhouding hoogstens 200.000 is. Wanneer de waarde van deze verhouding hoger is, is de snelheid z´o laag, dat het niet meer rendabel is.(zie A.2
5
Discussie
In het model is alleen rekening gehouden met de mogelijkheden tot opschalen van de hevel. Wanneer er onderzocht moet worden hoe de hevel op kleine schaal zich gedraagt, moet ook nog rekening worden gehouden met de eventuele capillaire werking.
4
A A.1
Berekeningen Verwerking formules
Hieronder is de toepassing van de wet van Bernoulli en de Darcy-Weisbach vergelijking worden verder uitgewerkt. In deze uitwerkingen worden de volgende variabelen gebruikt: v is de stroomsnelheid van de vloeistof in het punt (m/s) g is de valversnelling (m · s−2 ) y is het hoogteverschil (m) P is de druk in het punt (Pa = kg · m · s−2 ) ρ is de dichtheid van de vloeistof (kg · m−3 ) hf is de stijghoogte ten gevolge van wrijving (m) f is de Darcy-wrijvingsco¨effici¨ent (dimensieloos), gegeven door 64 , waarbij Re het getal van Reynolds is, gegeven door Re = ρvD f = Re µ ρ is de dichtheid van de vloeistof (kg · m−3 ) µ is de dynamische viscositeit van de vloeistof (Pa· s) L is de lengte van de buis (m) D is de hydraulische diameter van de buis (m) V is de gemiddelde stroomsnelheid van de vloeistof (m · s−1 ) g is de valversnelling (m · s−2 ) µ is de dynamische viscositeit (Pa · s)
Figuur 3: Betekenis hb en hc [1] 5
Pas de vergelijking van Bernoulli toe op het bovenste reservoir. De oppervlakte van het water daalt in feite continu, maar er wordt aangenomen dat deze gelijk blijft. Dit is in de meeste gevallen voldoende. Het model kan wel worden uitgebreid om ook het verloop van het waterniveau in bovenste reservoir mee te nemen. Verder wordt er vanuit gegaan dat de druk aan de oppervlakte (punt B) en bij de uitgang (punt c) gelijk zijn aan de atmosferische druk, en dat er wordt geheveld met vloeibaar water. Dit levert: Patm 02 + g(0) + = constant 2 ρ
(6)
Pas nu de vergelijking van Bernoulli toe op punt A, vA2 PA −g·d+ = constant 2 ρ
(7)
Pas nu de vergelijking van Bernoulli toe op punt B, vB2 PB + g · hB + = constant 2 ρ
(8)
Pas nu de vergelijking van Bernoulli toe op punt C, Patm vC2 − g · hC + = constant 2 ρ
(9)
Omdat de hevel ´e´en systeem is, moet deze constante overal gelijk zijn. Het combineren van vergelijkingen 6 en 9 levert de volgende vergelijking: 02 Patm v2 Patm +g·0+ = c − g · hc + 2 ρ 2 ρ
(10)
Dit levert vervolgens: vc =
p 2 · g · hc
(11)
Omdat de snelheid in de hevel overal even groot is levert dit: vc =
p 2 · g · hc
(12)
Maar aangezien in het hele systeem de snelheid gelijk is, wordt deze snelheid begrensd. In het punt B moet de druk positief zijn. Om deze reden worden de vergelijkingen 6 en 8 aan elkaar gelijk gesteld. Dit levert: vB2 PB 02 Patm + g · hB + = + g(0) + 2 ρ 2 ρ 6
(13)
Neem PB = 0. Dan geldt voor vmax : s Patm − g · hb vmax = 2 ρ
(14)
Door 13 op te lossen naar hb levert de volgende vergelijking voor hb : hb =
Patm v2 − ρ·g 2g
(15)
De hoogst mogelijke waarde voor hb wordt dus bereikt als v = 0, dus: hb,max =
A.2
Een aantal uitkomsten
7
Patm ρ·g
(16)
B 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44
Programma
/∗ ∗ ∗ ∗ @author s 1 1 1 1 8 6 ∗ Rik Schepens ∗/ import j a v a . u t i l . Scanner ; p u b l i c c l a s s Hevel { d o u b l e v , hc , hb , g , p , rho , l , d , mu, f , dp ; Scanner s c = new Scanner ( System . i n ) ; // Get t h e d e s c r i p t i o n by t h e u s e r , i n t h i s c a s e i t i s a l r e a d y g i v e n t h a t you u s e water on Earth . Changes can be made e a s i l y void getVariables ( ) { System . out . p r i n t l n ( ”What i s t h e a c c e l e r a t i o n due t o g r a v i t y ? (m/ s ˆ 2 ) ” ) ; // g = s c . nextDouble ( ) ; g = 9.81; System . out . p r i n t l n ( ”What i s t h e e l e v a t i o n d i f f e r e n c e between t h e s u r f a c e o f t h e top r e s e r v o i r and t h e s u r f a c e o f t h e bottom r e s e r v o i r ? (m) ” ) ; hc = s c . nextDouble ( ) ; System . out . p r i n t l n ( ”What i s t h e maximum s i p h o n e l e v a t i o n compared t o t h e top r e s e r v o i r ? (m) ” ) ; hb = s c . nextDouble ( ) ; System . out . p r i n t l n ( ”What i s t h e p r e s s u r e a t t h e c h o s e n p o i n t ? ( Pa ) ” ) ; //p = s c . nextDouble ( ) ; p = 101300; System . out . p r i n t l n ( ”What i s t h e d e n s i t y o f t h e used f l u i d ? ( kg /mˆ 3 ) ” ) ; // rho = s c . nextDouble ( ) ; rho = 1 0 0 0 ; System . out . p r i n t l n ( ”What i s t h e l e n g t h o f t h e p i p e ? (m) ” ) ; l = s c . nextDouble ( ) ; System . out . p r i n t l n ( ”What i s t h e h y d r a u l i c d i a m e t e r o f t h e p i p e ? (m) ” ) ; d = s c . nextDouble ( ) ; System . out . p r i n t l n ( ”What i s t h e dynamic v i s c o s i t y o f t h e f l u i d ? ( Pa∗ s ) ” ) ; //mu = s c . nextDouble ( ) ; mu = 0 . 0 0 0 8 9 ; } double calculateV ( ) { d o u b l e t = Math . s q r t ( 2 ∗ g ∗ hc ) ; i f ( t > calculateVmax ( ) ) { r e t u r n calculateVmax ( ) ; } else { return t ;
8
45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93
} } d o u b l e calculateVmax ( ) { r e t u r n Math . s q r t ( p / rho − g ∗ hb ) ; } double calculateHb ( ) { r e t u r n ( p / ( rho ∗ g ) − ( v ∗ v ) / ( 2 ∗ g ) ) ; } d o u b l e calculateHbMax ( ) { r e t u r n ( p / ( rho ∗ g ) ) ; } double c a l c u l a t e R e y no l d s ( ) { r e t u r n rho ∗ v ∗ d / mu; } double c a l c u l a t e F ( ) { r e t u r n 64 / c a l c u l a t e R e y n o l d s ( ) ; } double calculateDp ( ) { r e t u r n f ∗ l / d ∗ rho ∗ v ∗ v / 2 ; } // Use t h i s t o c a l c u l a t e what you want t o know . This i s an example we used a l o t . v o i d run ( ) { getVariables () ;
v = calculateV () ; System . out . p r i n t l n ( ” Without f r i c t i o n : ” ) ; System . out . p r i n t l n ( ”vmax = ” + calculateVmax ( ) ) ; System . out . p r i n t l n ( ”hbmax = ” + calculateHbMax ( ) ) ; f = calculateF () ; dp = c a l c u l a t e D p ( ) ; p = p − dp ; v = calculateV () ; System . out . p r i n t l n ( ”With f r i c t i o n : ” ) ; System . out . p r i n t l n ( ”v = ” + c a l c u l a t e V ( ) ) ; System . out . p r i n t l n ( ”vmax = ” + calculateVmax ( ) ) ; System . out . p r i n t l n ( ”hbmax = ” + calculateHbMax ( ) ) ; } p u b l i c s t a t i c v o i d main ( S t r i n g [ ] a r g s ) { new Hevel ( ) . run ( ) ; }
9
94
}
10
Referenties [1] Siphon. Wikipedia. URL: http://en.wikipedia.org/wiki/Siphon [2] Bernoulli’s equation. Wikipedia. URL: http://en.wikipedia.org/ wiki/Bernoulli’s_equation [3] Darcy–Weisbach equation. Wikipedia. URL: http://en.wikipedia. org/wiki/Darcy-Weisbach_equation [4] Darcy friction factor. Wikipedia. URL: http://en.wikipedia.org/ wiki/Darcy_friction_factor_formulae#Laminar_flow [5] Reynolds number. Wikipedia. URL: http://en.wikipedia.org/wiki/ Reynolds_number#Flow_in_pipe
11
