Universitaire natuurkunde — jaar 2 · Bachelor Year 2
4Viskeuze stromingen
Kantel een pot honing en hij loopt in een traag, dik lint uit; kantel een pot water en het is weg. Trek een lepel door beide: de honing verzet zich, en blijft zich bij elke snelheid verzetten; het water merkt er nauwelijks iets van tot u snel beweegt, en dan gaat het kolken. Het verschil is de viscositeit — de inwendige wrijving van een fluïdum — en de strijd tussen viscositeit en traagheid, samengevat in één getal, dat van Reynolds, beslist of een stroming in ordelijke lagen kruipt of in turbulentie omslaat. Dit hoofdstuk voegt de viskeuze kracht toe aan de vergelijking van Euler, lost de twee laminaire stromingen op die elke ingenieur gebruikt, en legt uit waarom de viscositeit zich bij een hoog reynoldsgetal in een dunne laag aan de wand verstopt en toch elke auto, elk schip en elk vliegtuig het grootste deel van zijn brandstof kost.
Honing loopt in een dun, gestaag lint uit dat bij de landing oprolt: een viscositeit tienduizend keer die van water maakt traagheid onbelangrijk — een kruipstroming.
4.1 Viscositeit en de vergelijking van Navier–Stokes
Definitie 4.1(Newtoniaanse viscositeit)
In een vlakke afschuifstroming v=vx(y)ex oefent het fluïdum boven een vlak y= const op het fluïdum eronder (door een oppervlakte-element dS) de tangentiële kracht
dF=ηdydvxdSex,
uit, waarmee het de tragere laag vooruittrekt en zelf wordt teruggetrokken. De schuifspanning is τ=ηdvx/dy; de coëfficiënt η (eenheid Pas) is de dynamische viscositeit van het fluïdum, en ν=η/ρ (m2/s) zijn kinematische viscositeit. Een fluïdum waarvoor η niet van de afschuifsnelheid afhangt, heet newtoniaans: water, lucht, oliën, alcohol; bloed, verf, ketchup en polymeersmelten zijn het niet. Ordegrootten bij 20∘C: lucht η=1.8×10−5Pas (ν=1.5×10−5m2/s), water 1.0×10−3Pas (ν=1.0×10−6m2/s), olijfolie 0.08Pas, glycerol 1.5Pas, honing ∼10Pas. Vloeistoffen worden dunner bij verwarming, gassen dikker.
In de vlakke afschuifstroming bedraagt de netto viskeuze kracht op een plakje fluïdum tussen y en y+dy, per volume-eenheid, ηd2vx/dy2: viscositeit werkt als een volumekrachtdichtheid ηΔv (met de laplaciaan component voor component genomen), een resultaat dat voor elke onsamendrukbare stroming geldt. Aan een vaste wand kleeft het fluïdum: zijn snelheid is die van de wand (de hechtvoorwaarde), en daarom blijft er stof op een ventilatorblad zitten en voelt een wand een viskeuze weerstand.
Bewijs. Het plakje wordt van boven met ηvx′(y+dy)dS getrokken en van onder met −ηvx′(y)dS: netto ηvx′′dydS. De algemene vorm wordt aangenomen (zij volgt uit het schrijven van de spanning als een symmetrische lineaire functie van de snelheidsgradiënten). De hechtvoorwaarde is experimenteel. ∎
Links: vlakke stroming van Couette tussen een vaste plaat en een plaat die met U beweegt — een lineair profiel, met dezelfde schuifspanning op elke laag. Rechts: de viskeuze krachten op een plakje fluïdum; zij heffen elkaar op tenzij het snelheidsprofiel gekromd is.
Voorbeeld 4.4(Schuifspanning op een tafel)
Een oliefilm van 0.1mm (η=0.1Pas) tussen een plaat en een tafel, met de plaat glijdend op 1m/s: τ=ηU/h=1kPa, een kracht van 100N op een tiende vierkante meter — het beginsel van het gesmeerde lager, waar de viskeuze weerstand de honderd keer grotere droge wrijving vervangt. Honing op een lepel, 10Pas, 1mm dik, aflopend met 1mm/s: τ=10Pa, precies het gewicht van een millimeter honing per oppervlakte-eenheid — en daarom loopt zij juist met die snelheid af.
4.2 Het reynoldsgetal
Definitie 4.5(Reynoldsgetal)
Voor een stroming met karakteristieke snelheid U en lengte L is het reynoldsgetal
Re=ηρUL=νUL
de verhouding van de ordegrootten van de traagheidsterm ρ(v⋅grad)v∼ρU2/L en de viskeuze term ηΔv∼ηU/L2 in de vergelijking van Navier–Stokes. Bij Re≪1 is de traagheid verwaarloosbaar en kruipt de stroming (viskeus, omkeerbaar, glad); bij Re≫1 is de viscositeit verwaarloosbaar behalve nabij wanden en is de stroming die van een volmaakt fluïdum — totdat zij, voorbij een drempel van orde 103 (ongeveer 2000 tot 3000 in een buis, op de diameter gebaseerd), turbulent wordt: niet-stationair, chaotisch, vol wervels op elke schaal. Daartussenin heet een stroming laminair: stationair en in lagen.
Voorbeeld 4.6(Reynoldsgetallen van alledaagse stromingen)
Een bacterie (1µm, 30µm/s, water): Re=3×10−5 — zij leeft in een wereld zonder traagheid, waar het stilzetten van haar zweepstaart haar binnen een atoomdiameter tot stilstand brengt. Bloed in een haarvat: 10−3; een vallende stofkorrel: 10−2; een goudvis: 103; een zwemmer: 106; een auto op 30m/s (4m): 8×106; een verkeersvliegtuig: 108; de Golfstroom: 1011. Twee stromingen met dezelfde meetkunde en hetzelfde reynoldsgetal zijn gelijkvormig: daarom kan een scheepsmodel van 1m in een sleeptank worden beproefd, of een vliegtuig in een windtunnel.
Propositie 4.7(Stroming van Poiseuille in een buis)
In een horizontale cilindrische buis met straal R en lengte L, gedreven door het drukverschil ΔP=P1−P2 tussen haar uiteinden, is de stationaire laminaire stroming
v(r)=4ηLΔP(R2−r2),DV=8ηLπR4ΔP:
een parabolisch profiel, maximaal op de as, met een gemiddelde snelheid ⟨v⟩=vmax/2. Het debiet is evenredig met de drukval (wet van Poiseuille), met de hydraulische weerstandRh=ΔP/DV=8ηL/πR4 — en met de vierde macht van de straal. De schuifspanning aan de wand is τw=ΔPR/2L=4η⟨v⟩/R. Geldig zolang Re=2R⟨v⟩/ν≲2000.
Bewijs. Stationair, v=v(r)ez (zodat de convectieve term wegvalt: v⋅grad werkt langs z, waarlangs niets verandert). Navier–Stokes langs z: 0=−∂zP+ηΔv, en radiaal ∂rP=0: P hangt alleen van z af, Δv alleen van r, zodat beide leden gelijk zijn aan een constante, ∂zP=−ΔP/L. In cilindercoördinaten is Δv=r1drd(rdrdv) (hier gesteld; het is de balans van de viskeuze krachten op een cilinderschil, 2πrLηv′ in r en in r+dr). Twee keer integreren met v eindig op de as en v(R)=0 geeft de parabool; DV=∫0Rv2πrdr. De wandspanning is −ηv′(R); even goed houdt de drukkracht πR2ΔP de wandwrijving 2πRLτw in evenwicht. ∎
Links: de stroming van Poiseuille — het parabolische profiel in een buis, gedreven door de drukval P1−P2. Rechts: de krachtenbalans op een cilinderschil fluïdum: het drukverschil op haar eindvlakken wordt in evenwicht gehouden door de viskeuze spanningen op haar binnen- en buitenoppervlak.
Voorbeeld 4.8(De vierde macht)
Een injectienaald met inwendige diameter 0.4mm en lengte 3cm, water, duimdruk 20kPa: DV=π(2×10−4)4×2×104/(8×10−3×0.03)=0.42mL/s; een twee keer zo brede naald geeft zestien keer meer. Een slagader die door een plaque tot 70% van haar straal is vernauwd, biedt (1/0.7)4=4.2 keer zo veel weerstand: het lichaam compenseert dat door de druk stroomopwaarts te verhogen, en dat is een deel van de reden waarom aderverkalking en hoge bloeddruk samengaan.
Propositie 4.9(Weerstand van Stokes)
Een bol met straal R die met snelheid v door een fluïdum beweegt bij Re=2ρvR/η≪1, ondervindt de weerstand
F=6πηRv,
tegengesteld aan zijn snelheid. Een bol met dichtheid ρs die valt in een fluïdum met dichtheid ρ bereikt de eindsnelheidv∞=9η2R2(ρs−ρ)g — de wet van de bezinking, van centrifuges en van de oliedruppels van Millikan.
Bewijs. De weerstand wordt aangenomen (zij is de exacte oplossing van de vergelijkingen van de kruipstroming om een bol, voor twee derde viskeuze wrijving en voor een derde druk). Eindsnelheid: gewicht min opwaartse kracht, 34πR3(ρs−ρ)g, is gelijk aan de weerstand; controleer achteraf altijd Re≪1. ∎
Voorbeeld 4.10(Mist, regen en stof)
Een mistdruppel met diameter 10µm valt in lucht met 2×(5×10−6)2×104/(9×1.8×10−5)=3mm/s (Re=2×10−3): de kleinste opwaartse luchtstroom houdt haar in de lucht, en daarom blijven wolken drijven. Een regendruppel van 1mm zou 30m/s geven en Re≈2000: Stokes geldt dan niet meer, en de kwadratische weerstand van de volgende paragraaf geeft de waargenomen 4m/s. Een stofkorrel van 10µm (ρs=2500kg/m3) bezinkt met 8mm/s, een meter in twee minuten — fijner stof blijft uren hangen.
4.4 Hoog reynoldsgetal: grenslagen en weerstand
Propositie 4.11(De grenslaag)
Bij een hoog reynoldsgetal is de stroming om een lichaam overal die van een volmaakt fluïdum, behalve in een dunne grenslaag langs de wand, waarin de snelheid van haar buitenwaarde U tot nul daalt (hechting). Langs een vlakke plaat groeit haar dikte als
δ(x)∼Uνx=Rexx,Rex=νUx,
met x de afstand tot de voorrand; de laag zelf wordt turbulent voorbij Rex∼5×105. De viskeuze wrijving op de wand, en het loslaten van de laag achter een stomp lichaam (dat een turbulent zog met lage druk achterlaat), zijn de twee bronnen van weerstand.
Bewijs. Ordegrootte: een fluïdumdeeltje doet er de tijd x/U over om langs de plaat te reizen, en ondertussen dringt de invloed van de wand het fluïdum binnen — de viskeuze term ν∂2v/∂y2 is een diffusie van impuls met diffusiecoëfficiënt ν, die de afstand νt bereikt in de tijd t (Hoofdstuk 24). Vandaar δ∼νx/U. De volledige theorie (van Prandtl) wordt aangenomen. ∎
Propositie 4.12(Weerstand bij een hoog reynoldsgetal)
De weerstand op een lichaam met frontaal oppervlak S dat met v door een fluïdum beweegt, wordt geschreven als
F=21ρv2SCx,
waarin de weerstandscoëfficiëntCx van de vorm en van Re afhangt. Voor een bol is Cx=24/Re bij Re≪1 (Stokes, met Re op de diameter gebaseerd), daarna ongeveer 0.4–0.5, vrijwel constant, van Re≈103 tot 2×105, waar zij plotseling tot ongeveer 0.1 valt (de weerstandscrisis: de grenslaag wordt turbulent, blijft langer aan het oppervlak kleven en het zog krimpt — de putjes van een golfbal lokken dat vroeger uit). Gestroomlijnde lichamen hebben Cx≈0.05; een auto 0.3; een fietser 0.9; een vlakke plaat dwars op de stroming 1.2. Het vermogen om met v te bewegen is Fv∝v3.
Een fietser (SCx=0.35m2) op 10m/s: F=0.5×1.2×100×0.35=21N, P=210W — een inspanning die vol te houden is; op 15m/s (54km/h) zou zij 710W vergen, en daarom rijden sprinters in het wiel. Een auto (SCx=0.7m2) op 130km/h: 550N, 20kW tegen de lucht alleen — het brandstofverbruik op de snelweg is grotendeels luchtweerstand.
Opmerking 4.14(Oppervlaktespanning)
Aan de grenzen van een vloeistof leeft een tweede kracht: haar vrije oppervlak gedraagt zich als een gespannen vlies met spanning γ (N/m; water 0.072N/m), omdat de moleculen aan het oppervlak de helft van hun buren missen. De druk binnen een druppel met straal R overtreft de druk erbuiten met de laplacedruk2γ/R (1.4kPa in een druppel van 0.1mm, 14bar in een van 10nm); de oppervlaktespanning wint het van de zwaartekracht beneden de capillaire lengteγ/ρg=2.7mm voor water — de maat van dauwdruppels, van de meniscus in een glas, van de stijging in een capillair buisje (h=2γ/ρgr, Jurin). Daarom lopen kleine insecten over water, daarom kleven de haren van een penseel samen als hij nat is, en daarom breekt een vloeistofstraal in druppels uiteen.
Methode 4.15(Welk regime?)
Bereken vóór elke berekening aan een viskeuze stroming eerst Re. Re≪1: Stokes en kruipstroming, weerstand ∝v, drukvallen ∝ debiet. Re tot ∼2000 in een buis: laminair, Poiseuille. Daarboven: turbulent — Poiseuille onderschat het verlies dan sterk; gebruik de empirische wet ΔP=λ(L/D)21ρ⟨v⟩2 met λ≈0.3/Re1/4 voor gladde buizen, en F=21ρv2SCx voor lichamen. Ga in elk geval achteraf na dat de Re die uit het antwoord volgt, overeenstemt met het veronderstelde regime.
4.5 Oefeningen
Oefening 4.1★
Een plaat van 0.50m2 glijdt met 2.0m/s over een oliefilm van 0.20mm (η=0.25Pas). Schuifspanning, kracht op de plaat, gedissipeerd vermogen; vergelijk met droge wrijving (f=0.3) als de plaat 10kg weegt.
Oplossing
Oplossing van Oefening 4.1.
τ=ηU/h=0.25×2/2×10−4=2.5kPa; F=1.25kN; P=Fv=2.5kW. Droge wrijving: 0.3×98=29N — hier verzet de dunne, sterk afgeschoven film zich veel meer; de viskeuze weerstand groeit als U/h, zodat smering loont bij lagere snelheden en dikkere films (en zij loopt nooit vast).
Oefening 4.2★
Reynoldsgetallen van: (a) een bacterie van 2µm bij 30µm/s in water; (b) een zwemmer (1.8m, 1.5m/s); (c) een auto (4m, 30m/s) in lucht; (d) bloed in de aorta (2.5cm, 0.33m/s, ν=3.3×10−6m2/s); (e) olie (η=0.1Pas, ρ=900kg/m3) bij 0.5m/s in een buis van 1cm. Geef het regime van elk.
Water wordt met 20kPa door een naald van 0.40mm (inwendige diameter) en 3.0cm geperst. Debiet; tijd om 5mL in te spuiten; maximale snelheid en reynoldsgetal; schuifspanning aan de wand. En met een naald van 0.80mm?
Oplossing
Oplossing van Oefening 4.3.
DV=πR4ΔP/8ηL=π(2×10−4)4×2×104/(8×10−3×0.03)=0.42mL/s; 12s voor 5mL. vmax=ΔPR2/4ηL=6.7m/s, ⟨v⟩=3.3m/s, Re=2R⟨v⟩/ν=1300: laminair. τw=ΔPR/2L=67Pa. Bij 0.8mm: zestien keer het debiet, 6.7mL/s — maar Re≈104: de stroming wordt turbulent en de winst is kleiner.
Oefening 4.4★
Eindsnelheden in lucht (η=1.8×10−5Pas) van (a) een stuifmeelkorrel van 20µm (ρs=1100kg/m3), (b) een rookdeeltje van 2µm (ρs=2000kg/m3), (c) een motregendruppel van 0.2mm. Controleer telkens Re; tijd om 1m te vallen.
Oplossing
Oplossing van Oefening 4.4.
v∞=2R2Δρg/9η. (a) 1.3cm/s, Re=0.02, 75s per meter. (b) 0.24mm/s, Re=3×10−5, 70min per meter. (c) 1.2m/s maar Re=16: Stokes overschat (de werkelijke snelheid is ongeveer 0.7m/s), ongeveer een seconde per meter.
Oefening 4.5★★
Een waterleiding met diameter 30cm voert 0.10m3/s over 1.0km. (a) Gemiddelde snelheid en reynoldsgetal. (b) Drukval als de wet van Poiseuille zou gelden. (c) Werkelijke val met λ=0.3/Re1/4; verhouding tot (b); pompvermogen. (d) Hetzelfde debiet in twee evenwijdige leidingen van 21cm (dezelfde totale doorsnede): val en vermogen.
Oplossing
Oplossing van Oefening 4.5.
(a) ⟨v⟩=0.1/0.0707=1.4m/s, Re=4.2×105. (b) 8ηLDV/πR4=500Pa. (c) λ=0.3/25.5=0.0118, ΔP=0.0118×3333×1000=39kPa: 80 keer meer; P=ΔPDV=3.9kW. (d) Dezelfde snelheid, Re=3×105, λ=0.0128, ΔP=0.0128×4717×1000=60kPa, 6kW: meer wand per eenheid debiet, meer verlies.
Oefening 4.6★★
Cellen met diameter 10µm en dichtheid 1050kg/m3 in water. (a) Bezinksnelheid onder de zwaartekracht; tijd om 5cm te zakken. (b) In een centrifuge bij 1000g. (c) Minimaal toerental voor 1000g op 10cm straal. (d) Waarom kunnen de grootste deeltjes het eerst worden afgescheiden?
Oplossing
Oplossing van Oefening 4.6.
(a) v=2×25×10−12×50×9.81/(9×10−3)=2.7µm/s; 5cm in 1.8×104 s, vijf uur. (b) 2.7mm/s, 18s. (c) Ω=1000g/r=313rad/s=3000rpm. (d) v∝R2: grote deeltjes slaan het eerst neer; opeenvolgende centrifugegangen met stijgend toerental scheiden de maten.
Oefening 4.7★★
Een fietser (SCx=0.35m2, rolweerstand 4N) rijdt 12m/s. (a) Luchtweerstand en totaal vermogen. (b) Bij 15m/s. (c) Tegen 5m/s tegenwind bij 12m/s grondsnelheid. (d) In het wiel rijden verlaagt Cx met 30%: bespaard vermogen. (e) Een helling van 5% af zonder trappen (massa 80kg): eindsnelheid.
Oplossing
Oplossing van Oefening 4.7.
(a) F=0.6×144×0.35=30N, +4N: P=34×12=410W. (b) 47+4=51N, 770W. (c) Relatieve snelheid 17m/s: 61+4=65N, maal de grondsnelheid: 780W. (d) Luchtweerstand 21N: 300W, een besparing van 110W. (e) mgsinα=39N=4+0.21v2: v=13m/s, 47km/h.
Oefening 4.8★★
Kwadratische weerstand. (a) Eindsnelheid van een golfbal (46g, 43mm, Cx=0.25) en van een parachutist (80kg, SCx=0.7m2, en daarna 25m2 met open parachute). (b) Bewegingsvergelijking voor een val uit rust met kwadratische weerstand; toon aan dat v=v∞tanh(gt/v∞). (c) Tijd en afstand waarin de parachutist 90% van v∞ bereikt.
Oplossing
Oplossing van Oefening 4.8.
(a) Golfbal: 21ρCxS=2.2×10−4kg/m, v∞=0.45/2.2×10−4=45m/s; parachutist 785/0.42=43m/s; met parachute 785/15=7.2m/s. (b) mv˙=mg−kv2, d.w.z. v˙=g(1−v2/v∞2), met als oplossing vanuit rust v∞tanh(gt/v∞). (c) tanhx=0.9 bij x=1.47: t=1.47×43/9.81=6.5s; z=(v∞2/g)lncoshx=190×0.83=160m.
Oefening 4.9★★
Viscosimeter van Couette. Een cilinder met straal 5.0cm en hoogte 10cm draait met 10rad/s binnen een vaste cilinder; de spleet van 1.0mm is gevuld met de te onderzoeken vloeistof; het benodigde koppel is 2.0×10−2Nm. (a) Behandel de spleet als vlakke stroming van Couette: afschuifsnelheid en schuifspanning. (b) Viscositeit. (c) Reynoldsgetal van de stroming in de spleet; is zij laminair? (d) Waarom wordt de spleet dun gemaakt?
Oplossing
Oplossing van Oefening 4.9.
(a) Afschuifsnelheid ΩR/e=0.5/10−3=500s−1; Γ=τ⋅2πRh⋅R: τ=0.02/(2π×2.5×10−3×0.1)=12.7Pa. (b) η=12.7/500=2.5×10−2Pas. (c) Re=ρUe/η≈1000×0.5×10−3/0.025=20: laminair. (d) Een dunne spleet maakt de afschuifsnelheid uniform (de vlakke benadering) en houdt de stroming laminair.
Oefening 4.10★★★
Een vaatnetwerk. Bloed (η=3.5mPas) stroomt met 5.0L/min. (a) Hydraulische weerstand van de aorta (straal 1.25cm, lengte 40cm) en drukval erlangs. (b) De 1010 haarvaten (straal 4µm, lengte 1mm) staan parallel: weerstand van één, van alle samen, en de drukval over het haarvatbed. (c) De arteriolen (108, straal 10µm, lengte 2mm): dezelfde vragen; waar zit de weerstand van de bloedsomloop? (d) De totale val is 13kPa: vermogen van de linkerkamer van het hart, en welke fractie van de 100W van het lichaam. (e) Een geneesmiddel verwijdt de arteriolen met 10%: gevolg voor de totale weerstand.
Oplossing
Oplossing van Oefening 4.10.
DV=8.3×10−5m3/s. (a) Rh=8ηL/πR4=8×3.5×10−3×0.4/(π×2.4×10−8)=1.5×105Pas/m3, ΔP=12Pa. (b) Eén haarvat 3.5×1016Pas/m3; 1010 parallel: 3.5×106Pas/m3; ΔP=290Pa (2mmHg). (c) Eén arteriole 8×3.5×10−3×2×10−3/(π×10−20)=1.8×1015Pas/m3; het geheel: 1.8×107Pas/m3; ΔP=1.5kPa: de arteriolen overheersen de twee andere (de drukval zit in de kleine vaten, niet in de grote). (d) P=13000×8.3×10−5=1.1W, ongeveer 1%. (e) (1/1.1)4=0.68: een derde van de arteriolaire weerstand weg — zo verlagen vaatverwijders de bloeddruk.
Oefening 4.11★★★
Grenslaag en wrijvingsweerstand. Een vlakke plaat met lengte L en breedte b beweegt met de rand vooruit met U. In het laminaire regime is de schuifspanning aan de wand τw=0.332ρU2/Rex (Blasius). (a) Toon aan dat de wrijvingskracht op één vlak F=0.664ρU2bL/ReL is. (b) Dikte van de grenslaag aan de achterrand van een plaat van 1m in water bij 0.3m/s, en de kracht op een plaat van 1m2 (beide vlakken). (c) Bij 10m/s: ReL, en waarom de formule dan niet meer geldt. (d) De romp van een schip (100m, 10m/s): schat de dikte van de grenslaag aan de achtersteven met de turbulente wet δ≈0.37x/Rex1/5.
Oplossing
Oplossing van Oefening 4.11.
(a) F=b∫0Lτwdx=0.332ρU2bν/U∫0Ldx/x=0.664ρU2bνL/U=0.664ρU2bL/ReL. (b) ReL=3×105 (laminair); δ≈5νL/U=9mm; F=0.664×1000×0.09×1/548=0.11N per vlak, 0.22N. (c) ReL=107: de laag is voorbij de eerste centimeters turbulent en de laminaire formule faalt (de werkelijke weerstand is enkele malen groter). (d) Rex=109: δ≈0.37×100/63=0.6m.
Oefening 4.12★★★
Een afstromende film. Een vloeistoffilm met dikte h stroomt stationair langs een verticale wand omlaag onder de zwaartekracht; de lucht oefent geen spanning uit op haar vrije oppervlak. (a) Schrijf Navier–Stokes op voor v=v(y)ez (y de afstand tot de wand, z omlaag) en de randvoorwaarden. (b) Toon aan dat v(y)=2ηρg(2hy−y2) en bereken het debiet per eenheid breedte q=ρgh3/3η. (c) Een verse laag verf (η=0.5Pas, ρ=1200kg/m3) van 0.2mm dik: snelheid aan het oppervlak en debiet; hoe ver beweegt het oppervlak in 10min? Waarom zakt verf uit als zij te dik wordt aangebracht? (d) Honing (η=10Pas) van 2mm dik op een verticale lepel: snelheid aan het oppervlak.
Oplossing
Oplossing van Oefening 4.12.
(a) 0=ρg+ηv′′(y), v(0)=0 (hechting), v′(h)=0 (geen spanning aan het vrije oppervlak). (b) Twee keer integreren: v=(ρg/2η)(2hy−y2); q=∫0hvdy=ρgh3/3η. (c) vs=ρgh2/2η=1200×9.81×4×10−8/1=0.47mm/s; q=6.3×10−8m2/s; 28cm in tien minuten — een newtoniaanse verf van die dikte zou eraf lopen; echte verven worden dunner bij afschuiving en zijn thixotroop, en zakken toch uit als zij te dik zijn, omdat vs∝h2. (d) 1400×9.81×4×10−6/20=2.7mm/s.
Een wervelstraat van von Kármán in de wolken stroomafwaarts van een eiland, gezien vanuit de ruimte: het zog van een cilinder bij een reynoldsgetal van enkele honderden, getekend op een schaal van honderd kilometer (NASA).
4.6 Vraagstuk: Bloed, ruwe olie en regen
Probleem 4.1
Weekendvraagstuk — dezelfde viskeuze wetten in een slagader, in een pijpleiding van duizend kilometer, en in een vallende druppel
Deel I — Bloed. Bloed: ρ=1060kg/m3, η=3.5×10−3Pas (als newtoniaans behandeld); hartminuutvolume DV=5.0L/min.
Als de stroming in de aorta die van Poiseuille was, wat zouden dan de maximale snelheid en de schuifspanning aan de wand zijn? (Het endotheel voelt spanningen van orde 1Pa.)
Een haarvat met straal 4.0µm voert bloed met 0.30mm/s: reynoldsgetal; drukval over zijn lengte van 1.0mm, in pascal en in millimeters kwik (1mmHg = 133Pa).
Tijd waarin een rode bloedcel het haarvat doorloopt; waarom is dat de juiste orde voor gasuitwisseling?
Arteriolen (straal 10µm, lengte 2.0mm, 108 parallel): weerstand van het geheel en de drukval erover bij het volle hartminuutvolume.
De totale val van slagader naar ader is 13kPa: mechanisch vermogen dat de linkerkamer levert. Vergelijk met een lichaam van 100W.
Een stenose verkleint de straal van een kransslagader over een korte lengte met 40%: factor waarmee haar weerstand stijgt; waarom de hartspier stroomafwaarts in rust nog wel, maar bij inspanning niet meer voldoende doorbloed kan zijn.
Bij zware inspanning bereikt het hartminuutvolume 25L/min: reynoldsgetal in de aorta, en regime.
Bloed is niet newtoniaans: zijn schijnbare viscositeit daalt bij een hoge afschuifsnelheid en stijgt bij een lage. Welk van de bovenstaande vaten wordt daardoor het sterkst beïnvloed, en in welke richting?
Deel II — De pijpleiding. Een pijpleiding voor ruwe olie: diameter D=1.2m, lengte L=1300km, debiet DV=3.0m3/s; olie ρ=850kg/m3, η=2.0×10−2Pas als zij warm wordt gehouden. Turbulente wrijving: ΔP=λ(L/D)21ρ⟨v⟩2, λ=0.316/Re1/4.
De buis verdraagt 80bar: minimaal aantal pompstations langs de leiding.
Totaal pompvermogen; energie per kubieke meter geleverde olie, vergeleken met de energie-inhoud van de olie (36GJ/m3).
Welke val zou de wet van Poiseuille voorspellen? Geef commentaar op de kosten van turbulentie (en op waarom in sommige pijpleidingen polymeren worden ingespoten die de turbulentie uitstellen).
In de winter zou de onverwarmde olie η=2.0Pas bereiken: nieuw reynoldsgetal en regime; drukval over de leiding met de passende wet. Conclusie voor het ontwerp (de echte leiding is geïsoleerd en de olie wordt warm gehouden).
Doorlooptijd van de olie van het ene uiteinde naar het andere.
Hetzelfde debiet in een buis van 1.0m: met welke factor verandert de turbulente drukval?
Een zandkorrel van 0.10mm (ρs=2650kg/m3) wordt in de warme olie meegevoerd: bezinksnelheid volgens Stokes, reynoldsgetal van de korrel, tijd om door de buis te bezinken. Vergelijk met de turbulente snelheidsschommelingen, enkele procenten van ⟨v⟩: bezinkt het zand in de draaiende pijpleiding? En in een stilgelegde?
Deel III — Druppels. Lucht: ρ=1.2kg/m3, η=1.8×10−5Pas; water ρw=1000kg/m3, γ=0.072N/m.
Wolkendruppel, diameter 10µm: eindsnelheid volgens Stokes, reynoldsgetal, tijd om 1km te vallen; waarom de wolk blijft hangen.
Motregen, 0.20mm: snelheid volgens Stokes en reynoldsgetal; is Stokes nog betrouwbaar?
Regendruppel, 2.0mm: toon aan dat Stokes faalt, en bepaal de eindsnelheid met Cx=0.5. Tijd om van 1km te vallen.
Stijgt of daalt de wolkendruppel in een opwaartse stroom van 1cm/s? Welke opwaartse stroom zou de motregendruppel dragen?
Laplacedruk binnen de drie druppels; vergelijk telkens met de dynamische druk21ρv∞2, en verklaar waarom grote regendruppels afplatten en boven ongeveer 5mm uiteenvallen.
Geef in één tabel voor de acht stromingen van dit vraagstuk het reynoldsgetal en de wet die het beheerst.
Oplossing
Oplossing van Probleem 4.1.
1.S=4.9cm2, ⟨v⟩=83/4.9=0.17m/s; Re=1060×0.17×0.025/3.5×10−3=1300: laminair (de pulserende pieken halen ∼5000).
6.P=ΔPDV=13000×8.3×10−5=1.1W, ongeveer 1% van het vermogen van het lichaam.
7.(1/0.6)4=7.7. In rust verwijden de vaten stroomafwaarts zich en compenseren dat; bij inspanning stijgt de vraag vier- tot vijfvoudig en kan de vernauwde slagader die niet leveren: angina pectoris.
8. Vijf keer de snelheid: Re≈6500, turbulent — de geruisen die een stethoscoop tijdens inspanning hoort.
9. In de trage stroming met lage afschuiving in kleine aders en haarvaten klonteren de rode bloedcellen samen en stijgt de schijnbare viscositeit (in de snelle aorta juist andersom); alle schattingen van de drukval in kleine vaten zijn dus ondergrenzen.
11.λ=0.316/19.2=0.0165; ΔP/L=(0.0165/1.2)×21×850×7.0=41Pa/m=41kPa/km; 534bar over de leiding.
12.534/80=6.7: minstens zeven stations.
13.P=ΔPDV=5.34×107×3=160MW; 53MJ/m3, 0.15% van de energie van de olie.
14.ΔP=32ηL⟨v⟩/D2=32×0.02×1.3×106×2.65/1.44=15bar: vijfendertig keer minder. Turbulentie is duur; enkele delen per miljoen aan lange polymeren stellen haar uit en verlagen de val met tientallen procenten.
15.Re=1350: laminair; Poiseuille: ΔP=32×2×1.3×106×2.65/1.44=1500bar — onmogelijk; de olie moet warm worden gehouden, vandaar de isolatie en de verwarming.
16.1.3×106/2.65=4.9×105s=5.7dagen.
17.τw=(ΔP/L)D/4=41×0.3=12Pa; F=τwπDL=6×107N.
18.ΔP∝λ⟨v⟩2/D∝D−4.75: 1.24.75=2.4.
19.v=2R2Δρg/9η=2×2.5×10−9×1800×9.81/0.18=0.49mm/s; Re=2×10−3; 41min om 1.2m over te steken. Turbulente schommelingen van enkele centimeters per seconde houden de korrel zwevend zolang de olie loopt; in een stilgelegde leiding bezinkt zij binnen het uur.
20.v=2×2.5×10−11×1000×9.81/1.62×10−4=3.0mm/s; Re=2×10−3; 1km in 3.3×105 s, vier dagen: elke opwaartse stroom boven enkele millimeters per seconde houdt de druppel vast.
21.1.2m/s, Re=16: Stokes overschat met ongeveer een factor twee (in werkelijkheid 0.7m/s).
22. Stokes zou 120m/s en Re∼104 geven: absurd. Kwadratische weerstand, 34πR3ρwg=21ρv2πR2Cx: v=8Rρwg/3ρCx=6.6m/s (Re=900, consistent); 150s vanaf 1km.
23.1cm/s>3mm/s: de druppel stijgt; de motregendruppel heeft ongeveer 1m/s nodig, wat convectieve wolken leveren.
24.2γ/R: 29kPa, 1.4kPa, 144Pa; dynamische drukken 21ρv2: 5×10−6, 0.3 en 26Pa. De verhouding bereikt 0.2 bij 2mm; bij 5mm (v≈9m/s, 21ρv2=49Pa tegen 58Pa) plat de aerodynamische druk de druppel af en breekt hem: geen regendruppel wordt groter dan ongeveer 5mm.