Physics · Boek 4 · Bachelor Year 2

Universitaire natuurkunde — jaar 2

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.
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\vect v = v_x(y)\,\vect e_x oefent het fluïdum boven een vlak y=y = const op het fluïdum eronder (door een oppervlakte-element  ⁣dS\dd S) de tangentiële kracht

 ⁣dF=η ⁣dvx ⁣dy ⁣dS  ex,\dd\vect F = \eta\,\frac{\dd v_x}{\dd y}\,\dd S\;\vect e_x ,

uit, waarmee het de tragere laag vooruittrekt en zelf wordt teruggetrokken. De schuifspanning is τ=η ⁣dvx/ ⁣dy\tau = \eta\,\dd v_x/\dd y; de coëfficiënt η\eta (eenheid Pas\mathrm{Pa}\,\mathrm{s}) is de dynamische viscositeit van het fluïdum, en ν=η/ρ\nu = \eta/\rho (m2/s\mathrm{m}^{2}/\mathrm{s}) zijn kinematische viscositeit. Een fluïdum waarvoor η\eta niet van de afschuifsnelheid afhangt, heet newtoniaans: water, lucht, oliën, alcohol; bloed, verf, ketchup en polymeersmelten zijn het niet. Ordegrootten bij 20C20{}^{\circ}\mathrm{C}: lucht η=1.8×105Pas\eta = 1.8 \times 10^{-5}\,\mathrm{Pa}\,\mathrm{s} (ν=1.5×105m2/s\nu = 1.5 \times 10^{-5}\,\mathrm{m}^{2}/\mathrm{s}), water 1.0×103Pas1.0 \times 10^{-3}\,\mathrm{Pa}\,\mathrm{s} (ν=1.0×106m2/s\nu = 1.0 \times 10^{-6}\,\mathrm{m}^{2}/\mathrm{s}), olijfolie 0.08Pas0.08\,\mathrm{Pa}\,\mathrm{s}, glycerol 1.5Pas1.5\,\mathrm{Pa}\,\mathrm{s}, honing 10Pas\sim10\,\mathrm{Pa}\,\mathrm{s}. Vloeistoffen worden dunner bij verwarming, gassen dikker.

Propositie 4.2 (Viskeuze krachtdichtheid; hechtvoorwaarde)

In de vlakke afschuifstroming bedraagt de netto viskeuze kracht op een plakje fluïdum tussen yy en y+ ⁣dyy + \dd y, per volume-eenheid, η ⁣d2vx/ ⁣dy2\eta\,\dd^2v_x/\dd y^2: viscositeit werkt als een volumekrachtdichtheid ηΔv\eta\,\Delta\vect 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\eta\,v_x'(y + \dd y)\dd S getrokken en van onder met ηvx(y) ⁣dS-\eta\,v_x'(y)\dd S: netto ηvx ⁣dy ⁣dS\eta v_x''\,\dd y\,\dd S. 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.

Stelling 4.3 (Vergelijking van Navier–Stokes)

Een onsamendrukbaar newtoniaans fluïdum voldoet aan

ρ(vt+(vgrad)v)=gradP+ρg+ηΔv,divv=0,\rho\Bigl(\frac{\partial\vect v}{\partial t} + (\vect v\cdot\operatorname{\vect{grad}}) \vect v\Bigr) = -\operatorname{\vect{grad}}P + \rho\vect g + \eta\,\Delta\vect v , \qquad \operatorname{div}\vect v = 0 ,

met de hechtvoorwaarde op elke wand. De vergelijking van Euler is het geval η=0\eta = 0.

Bewijs. Op dit niveau zonder bewijs aangenomen.

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.
Links: vlakke stroming van Couette tussen een vaste plaat en een plaat die met UU 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.1mm0.1\,\mathrm{mm} (η=0.1Pas\eta = 0.1\,\mathrm{Pa}\,\mathrm{s}) tussen een plaat en een tafel, met de plaat glijdend op 1m/s1\,\mathrm{m}/\mathrm{s}: τ=ηU/h=1kPa\tau = \eta U/h = 1\,\mathrm{kPa}, een kracht van 100N100\,\mathrm{N} 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, 10Pas10\,\mathrm{Pa}\,\mathrm{s}, 1mm1\,\mathrm{mm} dik, aflopend met 1mm/s1\,\mathrm{mm}/\mathrm{s}: τ=10Pa\tau = 10\,\mathrm{Pa}, 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 UU en lengte LL is het reynoldsgetal

Re=ρULη=ULν\mathrm{Re} = \frac{\rho UL}{\eta} = \frac{UL}{\nu}

de verhouding van de ordegrootten van de traagheidsterm ρ(vgrad)vρU2/L\rho(\vect v\cdot\operatorname{\vect{grad}})\vect v \sim \rho U^2/L en de viskeuze term ηΔvηU/L2\eta\Delta\vect v \sim \eta U/L^2 in de vergelijking van Navier–Stokes. Bij Re1\mathrm{Re} \ll 1 is de traagheid verwaarloosbaar en kruipt de stroming (viskeus, omkeerbaar, glad); bij Re1\mathrm{Re} \gg 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 10310^3 (ongeveer 20002000 tot 30003000 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µm1\,\text{µ}\mathrm{m}, 30µm/s30\,\text{µ}\mathrm{m}/\mathrm{s}, water): Re=3×105\mathrm{Re} = 3 \times 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: 10310^{-3}; een vallende stofkorrel: 10210^{-2}; een goudvis: 10310^3; een zwemmer: 10610^6; een auto op 30m/s30\,\mathrm{m}/\mathrm{s} (4m4\,\mathrm{m}): 8×1068 \times 10^6; een verkeersvliegtuig: 10810^8; de Golfstroom: 101110^{11}. Twee stromingen met dezelfde meetkunde en hetzelfde reynoldsgetal zijn gelijkvormig: daarom kan een scheepsmodel van 1m1\,\mathrm{m} in een sleeptank worden beproefd, of een vliegtuig in een windtunnel.

4.3 Laminaire stromingen: Couette, Poiseuille, Stokes

Propositie 4.7 (Stroming van Poiseuille in een buis)

In een horizontale cilindrische buis met straal RR en lengte LL, gedreven door het drukverschil ΔP=P1P2\Delta P = P_1 - P_2 tussen haar uiteinden, is de stationaire laminaire stroming

v(r)=ΔP4ηL(R2r2),DV=πR48ηLΔP:v(r) = \frac{\Delta P}{4\eta L}\,(R^2 - r^2) , \qquad D_V = \frac{\pi R^4}{8\eta L}\,\Delta P :

een parabolisch profiel, maximaal op de as, met een gemiddelde snelheid v=vmax/2\langle v\rangle = v_{\max}/2. Het debiet is evenredig met de drukval (wet van Poiseuille), met de hydraulische weerstand Rh=ΔP/DV=8ηL/πR4R_h = \Delta P/D_V = 8\eta L/\pi R^4 — en met de vierde macht van de straal. De schuifspanning aan de wand is τw=ΔPR/2L=4ηv/R\tau_w = \Delta P\,R/2L = 4\eta\langle v \rangle/R. Geldig zolang Re=2Rv/ν2000\mathrm{Re} = 2R\langle v\rangle/\nu \lesssim 2000.

Bewijs. Stationair, v=v(r)ez\vect v = v(r)\vect e_z (zodat de convectieve term wegvalt: vgrad\vect v \cdot\operatorname{\vect{grad}} werkt langs zz, waarlangs niets verandert). Navier–Stokes langs zz: 0=zP+ηΔv0 = -\partial_zP + \eta\Delta v, en radiaal rP=0\partial_rP = 0: PP hangt alleen van zz af, Δv\Delta v alleen van rr, zodat beide leden gelijk zijn aan een constante, zP=ΔP/L\partial_zP = -\Delta P/L. In cilindercoördinaten is Δv=1r ⁣d ⁣dr(r ⁣dv ⁣dr)\Delta v = \frac1r\frac{\dd}{\dd r}\bigl(r\frac{\dd v}{\dd r}\bigr) (hier gesteld; het is de balans van de viskeuze krachten op een cilinderschil, 2πrLηv2\pi rL\,\eta v' in rr en in r+ ⁣drr + \dd r). Twee keer integreren met vv eindig op de as en v(R)=0v(R) = 0 geeft de parabool; DV=0Rv2πr ⁣drD_V = \int_0^Rv\,2\pi r\,\dd r. De wandspanning is ηv(R)-\eta v'(R); even goed houdt de drukkracht πR2ΔP\pi R^2\Delta P de wandwrijving 2πRLτw2\pi RL\tau_w in evenwicht.

Links: de stroming van Poiseuille — het parabolische profiel in een buis, gedreven door de drukval P_1 - P_2. 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.
Links: de stroming van Poiseuille — het parabolische profiel in een buis, gedreven door de drukval P1P2P_1 - P_2. 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.4mm0.4\,\mathrm{mm} en lengte 3cm3\,\mathrm{cm}, water, duimdruk 20kPa20\,\mathrm{kPa}: DV=π(2×104)4×2×104/(8×103×0.03)=0.42mL/sD_V = \pi(2 \times 10^{-4})^4 \times 2 \times 10^4/(8 \times 10^{-3} \times 0.03) = 0.42\,\mathrm{mL}/\mathrm{s}; een twee keer zo brede naald geeft zestien keer meer. Een slagader die door een plaque tot 70%70\% van haar straal is vernauwd, biedt (1/0.7)4=4.2(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 RR die met snelheid vv door een fluïdum beweegt bij Re=2ρvR/η1\mathrm{Re} = 2\rho vR/\eta \ll 1, ondervindt de weerstand

F=6πηRv,F = 6\pi\eta Rv ,

tegengesteld aan zijn snelheid. Een bol met dichtheid ρs\rho_s die valt in een fluïdum met dichtheid ρ\rho bereikt de eindsnelheid v=2R2(ρsρ)g9ηv_\infty = \dfrac{2R^2(\rho_s - \rho)g}{9\eta} — 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, 43πR3(ρsρ)g\tfrac43\pi R^3(\rho_s - \rho)g, is gelijk aan de weerstand; controleer achteraf altijd Re1\mathrm{Re} \ll 1.

Voorbeeld 4.10 (Mist, regen en stof)

Een mistdruppel met diameter 10µm10\,\text{µ}\mathrm{m} valt in lucht met 2×(5×106)2×104/(9×1.8×105)=3mm/s2 \times (5 \times 10^{-6})^2 \times 10^4/(9 \times 1.8 \times 10^{-5}) = 3\,\mathrm{mm}/\mathrm{s} (Re=2×103\mathrm{Re} = 2 \times 10^{-3}): de kleinste opwaartse luchtstroom houdt haar in de lucht, en daarom blijven wolken drijven. Een regendruppel van 1mm1\,\mathrm{mm} zou 30m/s30\,\mathrm{m}/\mathrm{s} geven en Re2000\mathrm{Re} \approx 2000: Stokes geldt dan niet meer, en de kwadratische weerstand van de volgende paragraaf geeft de waargenomen 4m/s4\,\mathrm{m}/\mathrm{s}. Een stofkorrel van 10µm10\,\text{µ}\mathrm{m} (ρs=2500kg/m3\rho_s = 2500\,\mathrm{kg}/\mathrm{m}^{3}) bezinkt met 8mm/s8\,\mathrm{mm}/\mathrm{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 UU tot nul daalt (hechting). Langs een vlakke plaat groeit haar dikte als

δ(x)νxU=xRex,Rex=Uxν,\delta(x) \sim \sqrt{\frac{\nu x}{U}} = \frac{x}{\sqrt{\mathrm{Re}_x}} , \qquad \mathrm{Re}_x = \frac{Ux}{\nu} ,

met xx de afstand tot de voorrand; de laag zelf wordt turbulent voorbij Rex5×105\mathrm{Re}_x \sim 5 \times 10^5. 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/Ux/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\nu\,\partial^2v/\partial y^2 is een diffusie van impuls met diffusiecoëfficiënt ν\nu, die de afstand νt\sqrt{\nu t} bereikt in de tijd tt (Hoofdstuk 24). Vandaar δνx/U\delta \sim \sqrt{\nu 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 SS dat met vv door een fluïdum beweegt, wordt geschreven als

F=12ρv2SCx,F = \tfrac12\rho v^2\,S\,C_x ,

waarin de weerstandscoëfficiënt CxC_x van de vorm en van Re\mathrm{Re} afhangt. Voor een bol is Cx=24/ReC_x = 24/\mathrm{Re} bij Re1\mathrm{Re} \ll 1 (Stokes, met Re\mathrm{Re} op de diameter gebaseerd), daarna ongeveer 0.40.40.50.5, vrijwel constant, van Re103\mathrm{Re} \approx 10^3 tot 2×1052 \times 10^5, waar zij plotseling tot ongeveer 0.10.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 Cx0.05C_x \approx 0.05; een auto 0.30.3; een fietser 0.90.9; een vlakke plaat dwars op de stroming 1.21.2. Het vermogen om met vv te bewegen is Fvv3Fv \propto v^3.

Bewijs. Op dit niveau zonder bewijs aangenomen.

Links: de weerstandscoëfficiënt van een bol tegen het reynoldsgetal (logaritmische schalen): de wet van Stokes bij lage Re, een plateau bij 0.45, en de weerstandscrisis nabij Re = 2 × 105. Rechts: de grenslaag die langs een vlakke plaat dikker wordt; daarbuiten is de stroming die van een volmaakt fluïdum. Links: de weerstandscoëfficiënt van een bol tegen het reynoldsgetal (logaritmische schalen): de wet van Stokes bij lage Re, een plateau bij 0.45, en de weerstandscrisis nabij Re = 2 × 105. Rechts: de grenslaag die langs een vlakke plaat dikker wordt; daarbuiten is de stroming die van een volmaakt fluïdum.
Links: de weerstandscoëfficiënt van een bol tegen het reynoldsgetal (logaritmische schalen): de wet van Stokes bij lage Re\mathrm{Re}, een plateau bij 0.450.45, en de weerstandscrisis nabij Re=2×105\mathrm{Re} = 2 \times 10^5. Rechts: de grenslaag die langs een vlakke plaat dikker wordt; daarbuiten is de stroming die van een volmaakt fluïdum.

Voorbeeld 4.13 (Een fietser, een auto)

Een fietser (SCx=0.35m2SC_x = 0.35\,\mathrm{m}^{2}) op 10m/s10\,\mathrm{m}/\mathrm{s}: F=0.5×1.2×100×0.35=21NF = 0.5 \times 1.2 \times 100 \times 0.35 = 21\,\mathrm{N}, P=210WP = 210\,\mathrm{W} — een inspanning die vol te houden is; op 15m/s15\,\mathrm{m}/\mathrm{s} (54km/h54\,\mathrm{km}/\mathrm{h}) zou zij 710W710\,\mathrm{W} vergen, en daarom rijden sprinters in het wiel. Een auto (SCx=0.7m2SC_x = 0.7\,\mathrm{m}^{2}) op 130km/h130\,\mathrm{km}/\mathrm{h}: 550N550\,\mathrm{N}, 20kW20\,\mathrm{kW} 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 γ\gamma (N/m\mathrm{N}/\mathrm{m}; water 0.072N/m0.072\,\mathrm{N}/\mathrm{m}), omdat de moleculen aan het oppervlak de helft van hun buren missen. De druk binnen een druppel met straal RR overtreft de druk erbuiten met de laplacedruk 2γ/R2\gamma/R (1.4kPa1.4\,\mathrm{kPa} in een druppel van 0.1mm0.1\,\mathrm{mm}, 14bar14\,\mathrm{bar} in een van 10nm10\,\mathrm{nm}); de oppervlaktespanning wint het van de zwaartekracht beneden de capillaire lengte γ/ρg=2.7mm\sqrt{\gamma/\rho g} = 2.7\,\mathrm{mm} voor water — de maat van dauwdruppels, van de meniscus in een glas, van de stijging in een capillair buisje (h=2γ/ρgrh = 2\gamma/\rho 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\mathrm{Re}. Re1\mathrm{Re} \ll 1: Stokes en kruipstroming, weerstand v\propto v, drukvallen \propto debiet. Re\mathrm{Re} tot  ⁣2000\sim\!2000 in een buis: laminair, Poiseuille. Daarboven: turbulent — Poiseuille onderschat het verlies dan sterk; gebruik de empirische wet ΔP=λ(L/D)12ρv2\Delta P = \lambda(L/D)\,\tfrac12\rho\langle v\rangle^2 met λ0.3/Re1/4\lambda \approx 0.3/\mathrm{Re}^{1/4} voor gladde buizen, en F=12ρv2SCxF = \tfrac12\rho v^2SC_x voor lichamen. Ga in elk geval achteraf na dat de Re\mathrm{Re} die uit het antwoord volgt, overeenstemt met het veronderstelde regime.

4.5 Oefeningen

Oefening 4.1

Een plaat van 0.50m20.50\,\mathrm{m}^{2} glijdt met 2.0m/s2.0\,\mathrm{m}/\mathrm{s} over een oliefilm van 0.20mm0.20\,\mathrm{mm} (η=0.25Pas\eta = 0.25\,\mathrm{Pa}\,\mathrm{s}). Schuifspanning, kracht op de plaat, gedissipeerd vermogen; vergelijk met droge wrijving (f=0.3f = 0.3) als de plaat 10kg10\,\mathrm{kg} weegt.

Oplossing

Oplossing van Oefening 4.1.

τ=ηU/h=0.25×2/2×104=2.5kPa\tau = \eta U/h = 0.25 \times 2/2 \times 10^{-4} = 2.5\,\mathrm{kPa}; F=1.25kNF = 1.25\,\mathrm{kN}; P=Fv=2.5kWP = Fv = 2.5\,\mathrm{kW}. Droge wrijving: 0.3×98=29N0.3 \times 98 = 29\,\mathrm{N} — hier verzet de dunne, sterk afgeschoven film zich veel meer; de viskeuze weerstand groeit als U/hU/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µm2\,\text{µ}\mathrm{m} bij 30µm/s30\,\text{µ}\mathrm{m}/\mathrm{s} in water; (b) een zwemmer (1.8m1.8\,\mathrm{m}, 1.5m/s1.5\,\mathrm{m}/\mathrm{s}); (c) een auto (4m4\,\mathrm{m}, 30m/s30\,\mathrm{m}/\mathrm{s}) in lucht; (d) bloed in de aorta (2.5cm2.5\,\mathrm{cm}, 0.33m/s0.33\,\mathrm{m}/\mathrm{s}, ν=3.3×106m2/s\nu = 3.3 \times 10^{-6}\,\mathrm{m}^{2}/\mathrm{s}); (e) olie (η=0.1Pas\eta = 0.1\,\mathrm{Pa}\,\mathrm{s}, ρ=900kg/m3\rho = 900\,\mathrm{kg}/\mathrm{m}^{3}) bij 0.5m/s0.5\,\mathrm{m}/\mathrm{s} in een buis van 1cm1\,\mathrm{cm}. Geef het regime van elk.

Oplossing

Oplossing van Oefening 4.2.

(a) 2×106×3×105/106=6×1052 \times 10^{-6} \times 3 \times 10^{-5}/10^{-6} = 6 \times 10^{-5}, kruipend. (b) 1.8×1.5/106=2.7×1061.8 \times 1.5/10^{-6} = 2.7 \times 10^6, turbulent. (c) 4×30/1.5×105=8×1064 \times 30/1.5 \times 10^{-5} = 8 \times 10^6, turbulent. (d) 0.025×0.33/3.3×106=25000.025 \times 0.33/3.3 \times 10^{-6} = 2500, net laminair. (e) 900×0.5×0.01/0.1=45900 \times 0.5 \times 0.01/0.1 = 45, laminair.

Oefening 4.3

Water wordt met 20kPa20\,\mathrm{kPa} door een naald van 0.40mm0.40\,\mathrm{mm} (inwendige diameter) en 3.0cm3.0\,\mathrm{cm} geperst. Debiet; tijd om 5mL5\,\mathrm{mL} in te spuiten; maximale snelheid en reynoldsgetal; schuifspanning aan de wand. En met een naald van 0.80mm0.80\,\mathrm{mm}?

Oplossing

Oplossing van Oefening 4.3.

DV=πR4ΔP/8ηL=π(2×104)4×2×104/(8×103×0.03)=0.42mL/sD_V = \pi R^4\Delta P/8\eta L = \pi(2 \times 10^{-4})^4 \times 2 \times 10^4/(8 \times 10^{-3} \times 0.03) = 0.42\,\mathrm{mL}/\mathrm{s}; 12s12\,\mathrm{s} voor 5mL5\,\mathrm{mL}. vmax=ΔPR2/4ηL=6.7m/sv_{\max} = \Delta PR^2/ 4\eta L = 6.7\,\mathrm{m}/\mathrm{s}, v=3.3m/s\langle v\rangle = 3.3\,\mathrm{m}/\mathrm{s}, Re=2Rv/ν=1300\mathrm{Re} = 2R\langle v \rangle/\nu = 1300: laminair. τw=ΔPR/2L=67Pa\tau_w = \Delta PR/2L = 67\,\mathrm{Pa}. Bij 0.8mm0.8\,\mathrm{mm}: zestien keer het debiet, 6.7mL/s6.7\,\mathrm{mL}/\mathrm{s} — maar Re104\mathrm{Re} \approx 10^4: de stroming wordt turbulent en de winst is kleiner.

Oefening 4.4

Eindsnelheden in lucht (η=1.8×105Pas\eta = 1.8 \times 10^{-5}\,\mathrm{Pa}\,\mathrm{s}) van (a) een stuifmeelkorrel van 20µm20\,\text{µ}\mathrm{m} (ρs=1100kg/m3\rho_s = 1100\,\mathrm{kg}/\mathrm{m}^{3}), (b) een rookdeeltje van 2µm2\,\text{µ}\mathrm{m} (ρs=2000kg/m3\rho_s = 2000\,\mathrm{kg}/\mathrm{m}^{3}), (c) een motregendruppel van 0.2mm0.2\,\mathrm{mm}. Controleer telkens Re\mathrm{Re}; tijd om 1m1\,\mathrm{m} te vallen.

Oplossing

Oplossing van Oefening 4.4.

v=2R2Δρg/9ηv_\infty = 2R^2\Delta\rho\,g/9\eta. (a) 1.3cm/s1.3\,\mathrm{cm}/\mathrm{s}, Re=0.02\mathrm{Re} = 0.02, 75s75\,\mathrm{s} per meter. (b) 0.24mm/s0.24\,\mathrm{mm}/\mathrm{s}, Re=3×105\mathrm{Re} = 3 \times 10^{-5}, 70min70\,\mathrm{min} per meter. (c) 1.2m/s1.2\,\mathrm{m}/\mathrm{s} maar Re=16\mathrm{Re} = 16: Stokes overschat (de werkelijke snelheid is ongeveer 0.7m/s0.7\,\mathrm{m}/\mathrm{s}), ongeveer een seconde per meter.

Oefening 4.5 ★★

Een waterleiding met diameter 30cm30\,\mathrm{cm} voert 0.10m3/s0.10\,\mathrm{m}^{3}/\mathrm{s} over 1.0km1.0\,\mathrm{km}. (a) Gemiddelde snelheid en reynoldsgetal. (b) Drukval als de wet van Poiseuille zou gelden. (c) Werkelijke val met λ=0.3/Re1/4\lambda = 0.3/\mathrm{Re}^{1/4}; verhouding tot (b); pompvermogen. (d) Hetzelfde debiet in twee evenwijdige leidingen van 21cm21\,\mathrm{cm} (dezelfde totale doorsnede): val en vermogen.

Oplossing

Oplossing van Oefening 4.5.

(a) v=0.1/0.0707=1.4m/s\langle v\rangle = 0.1/0.0707 = 1.4\,\mathrm{m}/\mathrm{s}, Re=4.2×105\mathrm{Re} = 4.2 \times 10^5. (b) 8ηLDV/πR4=500Pa8\eta LD_V/\pi R^4 = 500\,\mathrm{Pa}. (c) λ=0.3/25.5=0.0118\lambda = 0.3/25.5 = 0.0118, ΔP=0.0118×3333×1000=39kPa\Delta P = 0.0118 \times 3333 \times 1000 = 39\,\mathrm{kPa}: 8080 keer meer; P=ΔPDV=3.9kWP = \Delta PD_V = 3.9\,\mathrm{kW}. (d) Dezelfde snelheid, Re=3×105\mathrm{Re} = 3 \times 10^5, λ=0.0128\lambda = 0.0128, ΔP=0.0128×4717×1000=60kPa\Delta P = 0.0128 \times 4717 \times 1000 = 60\,\mathrm{kPa}, 6kW6\,\mathrm{kW}: meer wand per eenheid debiet, meer verlies.

Oefening 4.6 ★★

Cellen met diameter 10µm10\,\text{µ}\mathrm{m} en dichtheid 1050kg/m31050\,\mathrm{kg}/\mathrm{m}^{3} in water. (a) Bezinksnelheid onder de zwaartekracht; tijd om 5cm5\,\mathrm{cm} te zakken. (b) In een centrifuge bij 1000g1000\,g. (c) Minimaal toerental voor 1000g1000\,g op 10cm10\,\mathrm{cm} straal. (d) Waarom kunnen de grootste deeltjes het eerst worden afgescheiden?

Oplossing

Oplossing van Oefening 4.6.

(a) v=2×25×1012×50×9.81/(9×103)=2.7µm/sv = 2 \times 25 \times 10^{-12} \times 50 \times 9.81/(9 \times 10^{-3}) = 2.7\,\text{µ}\mathrm{m}/\mathrm{s}; 5cm5\,\mathrm{cm} in 1.8×1041.8 \times 10^4 s, vijf uur. (b) 2.7mm/s2.7\,\mathrm{mm}/\mathrm{s}, 18s18\,\mathrm{s}. (c) Ω=1000g/r=313rad/s=3000rpm\Omega = \sqrt{1000g/r} = 313\,\mathrm{rad}/\mathrm{s} = 3000\,\mathrm{rpm}. (d) vR2v \propto R^2: grote deeltjes slaan het eerst neer; opeenvolgende centrifugegangen met stijgend toerental scheiden de maten.

Oefening 4.7 ★★

Een fietser (SCx=0.35m2SC_x = 0.35\,\mathrm{m}^{2}, rolweerstand 4N4\,\mathrm{N}) rijdt 12m/s12\,\mathrm{m}/\mathrm{s}. (a) Luchtweerstand en totaal vermogen. (b) Bij 15m/s15\,\mathrm{m}/\mathrm{s}. (c) Tegen 5m/s5\,\mathrm{m}/\mathrm{s} tegenwind bij 12m/s12\,\mathrm{m}/\mathrm{s} grondsnelheid. (d) In het wiel rijden verlaagt CxC_x met 30%30\%: bespaard vermogen. (e) Een helling van 5%5\% af zonder trappen (massa 80kg80\,\mathrm{kg}): eindsnelheid.

Oplossing

Oplossing van Oefening 4.7.

(a) F=0.6×144×0.35=30NF = 0.6 \times 144 \times 0.35 = 30\,\mathrm{N}, +4N+4\,\mathrm{N}: P=34×12=410WP = 34 \times 12 = 410\,\mathrm{W}. (b) 47+4=51N47 + 4 = 51\,\mathrm{N}, 770W770\,\mathrm{W}. (c) Relatieve snelheid 17m/s17\,\mathrm{m}/\mathrm{s}: 61+4=65N61 + 4 = 65\,\mathrm{N}, maal de grondsnelheid: 780W780\,\mathrm{W}. (d) Luchtweerstand 21N21\,\mathrm{N}: 300W300\,\mathrm{W}, een besparing van 110W110\,\mathrm{W}. (e) mgsinα=39N=4+0.21v2mg\sin\alpha = 39\,\mathrm{N} = 4 + 0.21v^2: v=13m/sv = 13\,\mathrm{m}/\mathrm{s}, 47km/h47\,\mathrm{km}/\mathrm{h}.

Oefening 4.8 ★★

Kwadratische weerstand. (a) Eindsnelheid van een golfbal (46g46\,\mathrm{g}, 43mm43\,\mathrm{mm}, Cx=0.25C_x = 0.25) en van een parachutist (80kg80\,\mathrm{kg}, SCx=0.7m2SC_x = 0.7\,\mathrm{m}^{2}, en daarna 25m225\,\mathrm{m}^{2} met open parachute). (b) Bewegingsvergelijking voor een val uit rust met kwadratische weerstand; toon aan dat v=vtanh(gt/v)v = v_\infty\tanh(gt/v_\infty). (c) Tijd en afstand waarin de parachutist 90%90\% van vv_\infty bereikt.

Oplossing

Oplossing van Oefening 4.8.

(a) Golfbal: 12ρCxS=2.2×104kg/m\tfrac12\rho C_xS = 2.2 \times 10^{-4}\,\mathrm{kg}/\mathrm{m}, v=0.45/2.2×104=45m/sv_\infty = \sqrt{0.45/2.2 \times 10^{-4}} = 45\,\mathrm{m}/\mathrm{s}; parachutist 785/0.42=43m/s\sqrt{785/0.42} = 43\,\mathrm{m}/\mathrm{s}; met parachute 785/15=7.2m/s\sqrt{785/15} = 7.2\,\mathrm{m}/\mathrm{s}. (b) mv˙=mgkv2m\dot v = mg - kv^2, d.w.z. v˙=g(1v2/v2)\dot v = g(1 - v^2/v_\infty^2), met als oplossing vanuit rust vtanh(gt/v)v_\infty\tanh(gt/v_\infty). (c) tanhx=0.9\tanh x = 0.9 bij x=1.47x = 1.47: t=1.47×43/9.81=6.5st = 1.47 \times 43/9.81 = 6.5\,\mathrm{s}; z=(v2/g)lncoshx=190×0.83=160mz = (v_\infty^2/g)\ln\cosh x = 190 \times 0.83 = 160\,\mathrm{m}.

Oefening 4.9 ★★

Viscosimeter van Couette. Een cilinder met straal 5.0cm5.0\,\mathrm{cm} en hoogte 10cm10\,\mathrm{cm} draait met 10rad/s10\,\mathrm{rad}/\mathrm{s} binnen een vaste cilinder; de spleet van 1.0mm1.0\,\mathrm{mm} is gevuld met de te onderzoeken vloeistof; het benodigde koppel is 2.0×102Nm2.0 \times 10^{-2}\,\mathrm{N}\,\mathrm{m}. (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/103=500s1\Omega R/e = 0.5/10^{-3} = 500\,\mathrm{s}^{-1}; Γ=τ2πRhR\Gamma = \tau\cdot2\pi Rh \cdot R: τ=0.02/(2π×2.5×103×0.1)=12.7Pa\tau = 0.02/(2\pi \times 2.5 \times 10^{-3} \times 0.1) = 12.7\,\mathrm{Pa}. (b) η=12.7/500=2.5×102Pas\eta = 12.7/500 = 2.5 \times 10^{-2}\,\mathrm{Pa}\,\mathrm{s}. (c) Re=ρUe/η1000×0.5×103/0.025=20\mathrm{Re} = \rho Ue/\eta \approx 1000 \times 0.5 \times 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\eta = 3.5\,\mathrm{mPa}\,\mathrm{s}) stroomt met 5.0L/min5.0\,\mathrm{L}/\mathrm{min}. (a) Hydraulische weerstand van de aorta (straal 1.25cm1.25\,\mathrm{cm}, lengte 40cm40\,\mathrm{cm}) en drukval erlangs. (b) De 101010^{10} haarvaten (straal 4µm4\,\text{µ}\mathrm{m}, lengte 1mm1\,\mathrm{mm}) staan parallel: weerstand van één, van alle samen, en de drukval over het haarvatbed. (c) De arteriolen (10810^8, straal 10µm10\,\text{µ}\mathrm{m}, lengte 2mm2\,\mathrm{mm}): dezelfde vragen; waar zit de weerstand van de bloedsomloop? (d) De totale val is 13kPa13\,\mathrm{kPa}: vermogen van de linkerkamer van het hart, en welke fractie van de 100W100\,\mathrm{W} van het lichaam. (e) Een geneesmiddel verwijdt de arteriolen met 10%10\%: gevolg voor de totale weerstand.

Oplossing

Oplossing van Oefening 4.10.

DV=8.3×105m3/sD_V = 8.3 \times 10^{-5}\,\mathrm{m}^{3}/\mathrm{s}. (a) Rh=8ηL/πR4=8×3.5×103×0.4/(π×2.4×108)=1.5×105Pas/m3R_h = 8\eta L/\pi R^4 = 8 \times 3.5 \times 10^{-3} \times 0.4/(\pi \times 2.4 \times 10^{-8}) = 1.5 \times 10^{5}\,\mathrm{Pa}\,\mathrm{s}/\mathrm{m}^{3}, ΔP=12Pa\Delta P = 12\,\mathrm{Pa}. (b) Eén haarvat 3.5×1016Pas/m33.5 \times 10^{16}\,\mathrm{Pa}\,\mathrm{s}/\mathrm{m}^{3}; 101010^{10} parallel: 3.5×106Pas/m33.5 \times 10^{6}\,\mathrm{Pa}\,\mathrm{s}/\mathrm{m}^{3}; ΔP=290Pa\Delta P = 290\,\mathrm{Pa} (2mmHg2\,\mathrm{mmHg}). (c) Eén arteriole 8×3.5×103×2×103/(π×1020)=1.8×1015Pas/m38 \times 3.5 \times 10^{-3} \times 2 \times 10^{-3}/(\pi \times 10^{-20}) = 1.8 \times 10^{15}\,\mathrm{Pa}\,\mathrm{s}/\mathrm{m}^{3}; het geheel: 1.8×107Pas/m31.8 \times 10^{7}\,\mathrm{Pa}\,\mathrm{s}/\mathrm{m}^{3}; ΔP=1.5kPa\Delta P = 1.5\,\mathrm{kPa}: de arteriolen overheersen de twee andere (de drukval zit in de kleine vaten, niet in de grote). (d) P=13000×8.3×105=1.1WP = 13000 \times 8.3 \times 10^{-5} = 1.1\,\mathrm{W}, ongeveer 1%1\%. (e) (1/1.1)4=0.68(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 LL en breedte bb beweegt met de rand vooruit met UU. In het laminaire regime is de schuifspanning aan de wand τw=0.332ρU2/Rex\tau_w = 0.332\,\rho U^2/\sqrt{\mathrm{Re}_x} (Blasius). (a) Toon aan dat de wrijvingskracht op één vlak F=0.664ρU2bL/ReLF = 0.664\,\rho U^2bL/\sqrt{\mathrm{Re}_L} is. (b) Dikte van de grenslaag aan de achterrand van een plaat van 1m1\,\mathrm{m} in water bij 0.3m/s0.3\,\mathrm{m}/\mathrm{s}, en de kracht op een plaat van 1m21\,\mathrm{m}^{2} (beide vlakken). (c) Bij 10m/s10\,\mathrm{m}/\mathrm{s}: ReL\mathrm{Re}_L, en waarom de formule dan niet meer geldt. (d) De romp van een schip (100m100\,\mathrm{m}, 10m/s10\,\mathrm{m}/\mathrm{s}): schat de dikte van de grenslaag aan de achtersteven met de turbulente wet δ0.37x/Rex1/5\delta \approx 0.37x/\mathrm{Re}_x^{1/5}.

Oplossing

Oplossing van Oefening 4.11.

(a) F=b0Lτw ⁣dx=0.332ρU2bν/U0L ⁣dx/x=0.664ρU2bνL/U=0.664ρU2bL/ReLF = b\int_0^L\tau_w\,\dd x = 0.332\rho U^2b\sqrt{\nu/U}\int_0^L\dd x/\sqrt x = 0.664\rho U^2b\sqrt{\nu L/U} = 0.664\rho U^2bL/\sqrt{\mathrm{Re}_L}. (b) ReL=3×105\mathrm{Re}_L = 3 \times 10^5 (laminair); δ5νL/U=9mm\delta \approx 5\sqrt{\nu L/U} = 9\,\mathrm{mm}; F=0.664×1000×0.09×1/548=0.11NF = 0.664 \times 1000 \times 0.09 \times 1/548 = 0.11\,\mathrm{N} per vlak, 0.22N0.22\,\mathrm{N}. (c) ReL=107\mathrm{Re}_L = 10^7: de laag is voorbij de eerste centimeters turbulent en de laminaire formule faalt (de werkelijke weerstand is enkele malen groter). (d) Rex=109\mathrm{Re}_x = 10^9: δ0.37×100/63=0.6m\delta \approx 0.37 \times 100/63 = 0.6\,\mathrm{m}.

Oefening 4.12 ★★★

Een afstromende film. Een vloeistoffilm met dikte hh 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\vect v = v(y)\vect e_z (yy de afstand tot de wand, zz omlaag) en de randvoorwaarden. (b) Toon aan dat v(y)=ρg2η(2hyy2)v(y) = \dfrac{\rho g}{2\eta}(2hy - y^2) en bereken het debiet per eenheid breedte q=ρgh3/3ηq = \rho gh^3/3\eta. (c) Een verse laag verf (η=0.5Pas\eta = 0.5\,\mathrm{Pa}\,\mathrm{s}, ρ=1200kg/m3\rho = 1200\,\mathrm{kg}/\mathrm{m}^{3}) van 0.2mm0.2\,\mathrm{mm} dik: snelheid aan het oppervlak en debiet; hoe ver beweegt het oppervlak in 10min10\,\mathrm{min}? Waarom zakt verf uit als zij te dik wordt aangebracht? (d) Honing (η=10Pas\eta = 10\,\mathrm{Pa}\,\mathrm{s}) van 2mm2\,\mathrm{mm} dik op een verticale lepel: snelheid aan het oppervlak.

Oplossing

Oplossing van Oefening 4.12.

(a) 0=ρg+ηv(y)0 = \rho g + \eta v''(y), v(0)=0v(0) = 0 (hechting), v(h)=0v'(h) = 0 (geen spanning aan het vrije oppervlak). (b) Twee keer integreren: v=(ρg/2η)(2hyy2)v = (\rho g/2\eta)(2hy - y^2); q=0hv ⁣dy=ρgh3/3ηq = \int_0^hv\,\dd y = \rho gh^3/3\eta. (c) vs=ρgh2/2η=1200×9.81×4×108/1=0.47mm/sv_s = \rho gh^2/2\eta = 1200 \times 9.81 \times 4 \times 10^{-8}/1 = 0.47\,\mathrm{mm}/\mathrm{s}; q=6.3×108m2/sq = 6.3 \times 10^{-8}\,\mathrm{m}^{2}/\mathrm{s}; 28cm28\,\mathrm{cm} 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 vsh2v_s \propto h^2. (d) 1400×9.81×4×106/20=2.7mm/s1400 \times 9.81 \times 4 \times 10^{-6}/20 = 2.7\,\mathrm{mm}/\mathrm{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).
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\rho = 1060\,\mathrm{kg}/\mathrm{m}^{3}, η=3.5×103Pas\eta = 3.5 \times 10^{-3}\,\mathrm{Pa}\,\mathrm{s} (als newtoniaans behandeld); hartminuutvolume DV=5.0L/minD_V = 5.0\,\mathrm{L}/\mathrm{min}.

  1. Aorta, straal 1.25cm1.25\,\mathrm{cm}: gemiddelde snelheid, reynoldsgetal, regime.
  2. 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 1Pa1\,\mathrm{Pa}.)
  3. Een haarvat met straal 4.0µm4.0\,\text{µ}\mathrm{m} voert bloed met 0.30mm/s0.30\,\mathrm{mm}/\mathrm{s}: reynoldsgetal; drukval over zijn lengte van 1.0mm1.0\,\mathrm{mm}, in pascal en in millimeters kwik (1mmHg1\,\mathrm{mmHg} = 133Pa133\,\mathrm{Pa}).
  4. Tijd waarin een rode bloedcel het haarvat doorloopt; waarom is dat de juiste orde voor gasuitwisseling?
  5. Arteriolen (straal 10µm10\,\text{µ}\mathrm{m}, lengte 2.0mm2.0\,\mathrm{mm}, 10810^8 parallel): weerstand van het geheel en de drukval erover bij het volle hartminuutvolume.
  6. De totale val van slagader naar ader is 13kPa13\,\mathrm{kPa}: mechanisch vermogen dat de linkerkamer levert. Vergelijk met een lichaam van 100W100\,\mathrm{W}.
  7. Een stenose verkleint de straal van een kransslagader over een korte lengte met 40%40\%: factor waarmee haar weerstand stijgt; waarom de hartspier stroomafwaarts in rust nog wel, maar bij inspanning niet meer voldoende doorbloed kan zijn.
  8. Bij zware inspanning bereikt het hartminuutvolume 25L/min25\,\mathrm{L}/\mathrm{min}: reynoldsgetal in de aorta, en regime.
  9. 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.2mD = 1.2\,\mathrm{m}, lengte L=1300kmL = 1300\,\mathrm{km}, debiet DV=3.0m3/sD_V = 3.0\,\mathrm{m}^{3}/\mathrm{s}; olie ρ=850kg/m3\rho = 850\,\mathrm{kg}/\mathrm{m}^{3}, η=2.0×102Pas\eta = 2.0 \times 10^{-2}\,\mathrm{Pa}\,\mathrm{s} als zij warm wordt gehouden. Turbulente wrijving: ΔP=λ(L/D)12ρv2\Delta P = \lambda (L/D)\tfrac12\rho\langle v\rangle^2, λ=0.316/Re1/4\lambda = 0.316/\mathrm{Re}^{1/4}.

  1. Gemiddelde snelheid en reynoldsgetal; regime.
  2. Drukval per kilometer en over de hele leiding.
  3. De buis verdraagt 80bar80\,\mathrm{bar}: minimaal aantal pompstations langs de leiding.
  4. Totaal pompvermogen; energie per kubieke meter geleverde olie, vergeleken met de energie-inhoud van de olie (36GJ/m336\,\mathrm{GJ}/\mathrm{m}^{3}).
  5. 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).
  6. In de winter zou de onverwarmde olie η=2.0Pas\eta = 2.0\,\mathrm{Pa}\,\mathrm{s} 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).
  7. Doorlooptijd van de olie van het ene uiteinde naar het andere.
  8. Schuifspanning aan de wand in het warme, turbulente geval, uit de drukval; totale wrijvingskracht op de buiswand.
  9. Hetzelfde debiet in een buis van 1.0m1.0\,\mathrm{m}: met welke factor verandert de turbulente drukval?
  10. Een zandkorrel van 0.10mm0.10\,\mathrm{mm} (ρs=2650kg/m3\rho_s = 2650\,\mathrm{kg}/\mathrm{m}^{3}) 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\langle v \rangle: bezinkt het zand in de draaiende pijpleiding? En in een stilgelegde?

Deel III — Druppels. Lucht: ρ=1.2kg/m3\rho = 1.2\,\mathrm{kg}/\mathrm{m}^{3}, η=1.8×105Pas\eta = 1.8 \times 10^{-5}\,\mathrm{Pa}\,\mathrm{s}; water ρw=1000kg/m3\rho_w = 1000\,\mathrm{kg}/\mathrm{m}^{3}, γ=0.072N/m\gamma = 0.072\,\mathrm{N}/\mathrm{m}.

  1. Wolkendruppel, diameter 10µm10\,\text{µ}\mathrm{m}: eindsnelheid volgens Stokes, reynoldsgetal, tijd om 1km1\,\mathrm{km} te vallen; waarom de wolk blijft hangen.
  2. Motregen, 0.20mm0.20\,\mathrm{mm}: snelheid volgens Stokes en reynoldsgetal; is Stokes nog betrouwbaar?
  3. Regendruppel, 2.0mm2.0\,\mathrm{mm}: toon aan dat Stokes faalt, en bepaal de eindsnelheid met Cx=0.5C_x = 0.5. Tijd om van 1km1\,\mathrm{km} te vallen.
  4. Stijgt of daalt de wolkendruppel in een opwaartse stroom van 1cm/s1\,\mathrm{cm}/\mathrm{s}? Welke opwaartse stroom zou de motregendruppel dragen?
  5. Laplacedruk binnen de drie druppels; vergelijk telkens met de dynamische druk 12ρv2\tfrac12\rho v_\infty^2, en verklaar waarom grote regendruppels afplatten en boven ongeveer 5mm5\,\mathrm{mm} uiteenvallen.
  6. 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.9cm2S = 4.9\,\mathrm{cm}^{2}, v=83/4.9=0.17m/s\langle v\rangle = 83/4.9 = 0.17\,\mathrm{m}/\mathrm{s}; Re=1060×0.17×0.025/3.5×103=1300\mathrm{Re} = 1060 \times 0.17 \times 0.025/3.5 \times 10^{-3} = 1300: laminair (de pulserende pieken halen 5000\sim 5000).

2. vmax=0.34m/sv_{\max} = 0.34\,\mathrm{m}/\mathrm{s}; τw=4ηv/R=0.19Pa\tau_w = 4\eta\langle v\rangle/R = 0.19\,\mathrm{Pa}.

3. Re=1060×3×104×8×106/3.5×103=7×104\mathrm{Re} = 1060 \times 3 \times 10^{-4} \times 8 \times 10^{-6}/3.5 \times 10^{-3} = 7 \times 10^{-4}; ΔP=8ηLv/R2=8×3.5×103×103×3×104/1.6×1011=525Pa=3.9mmHg\Delta P = 8\eta L\langle v\rangle/R^2 = 8 \times 3.5 \times 10^{-3} \times 10^{-3} \times 3 \times 10^{-4}/1.6 \times 10^{-11} = 525\,\mathrm{Pa} = 3.9\,\mathrm{mmHg}.

4. 1/0.3=3.3s1/0.3 = 3.3\,\mathrm{s}: lang genoeg opdat zuurstof de enkele micrometers naar het weefsel kan diffunderen (Hoofdstuk 24).

5. Eén arteriole 8ηL/πR4=1.8×1015Pas/m38\eta L/\pi R^4 = 1.8 \times 10^{15}\,\mathrm{Pa}\,\mathrm{s}/\mathrm{m}^{3}; 10810^8 parallel: 1.8×107Pas/m31.8 \times 10^{7}\,\mathrm{Pa}\,\mathrm{s}/\mathrm{m}^{3}; ΔP=1.8×107×8.3×105=1.5kPa=11mmHg\Delta P = 1.8 \times 10^7 \times 8.3 \times 10^{-5} = 1.5\,\mathrm{kPa} = 11\,\mathrm{mmHg}.

6. P=ΔPDV=13000×8.3×105=1.1WP = \Delta P\,D_V = 13000 \times 8.3 \times 10^{-5} = 1.1\,\mathrm{W}, ongeveer 1%1\% van het vermogen van het lichaam.

7. (1/0.6)4=7.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: Re6500\mathrm{Re} \approx 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.

10. v=3/1.13=2.65m/s\langle v\rangle = 3/1.13 = 2.65\,\mathrm{m}/\mathrm{s}; Re=850×2.65×1.2/0.02=1.35×105\mathrm{Re} = 850 \times 2.65 \times 1.2/0.02 = 1.35 \times 10^5: turbulent.

11. λ=0.316/19.2=0.0165\lambda = 0.316/19.2 = 0.0165; ΔP/L=(0.0165/1.2)×12×850×7.0=41Pa/m=41kPa/km\Delta P/L = (0.0165/1.2) \times \tfrac12 \times 850 \times 7.0 = 41\,\mathrm{Pa}/\mathrm{m} = 41\,\mathrm{kPa}/\mathrm{km}; 534bar534\,\mathrm{bar} over de leiding.

12. 534/80=6.7534/80 = 6.7: minstens zeven stations.

13. P=ΔPDV=5.34×107×3=160MWP = \Delta P\,D_V = 5.34 \times 10^7 \times 3 = 160\,\mathrm{MW}; 53MJ/m353\,\mathrm{MJ}/\mathrm{m}^{3}, 0.15%0.15\% van de energie van de olie.

14. ΔP=32ηLv/D2=32×0.02×1.3×106×2.65/1.44=15bar\Delta P = 32\eta L\langle v\rangle/D^2 = 32 \times 0.02 \times 1.3 \times 10^6 \times 2.65/1.44 = 15\,\mathrm{bar}: 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\mathrm{Re} = 1350: laminair; Poiseuille: ΔP=32×2×1.3×106×2.65/1.44=1500bar\Delta P = 32 \times 2 \times 1.3 \times 10^6 \times 2.65/1.44 = 1500\,\mathrm{bar} — onmogelijk; de olie moet warm worden gehouden, vandaar de isolatie en de verwarming.

16. 1.3×106/2.65=4.9×105s=5.7dagen1.3 \times 10^6/2.65 = 4.9 \times 10^{5}\,\mathrm{s} = 5.7\,\mathrm{dagen}.

17. τw=(ΔP/L)D/4=41×0.3=12Pa\tau_w = (\Delta P/L)D/4 = 41 \times 0.3 = 12\,\mathrm{Pa}; F=τwπDL=6×107NF = \tau_w\pi DL = 6 \times 10^{7}\,\mathrm{N}.

18. ΔPλv2/DD4.75\Delta P \propto \lambda\langle v\rangle^2/D \propto D^{-4.75}: 1.24.75=2.41.2^{4.75} = 2.4.

19. v=2R2Δρg/9η=2×2.5×109×1800×9.81/0.18=0.49mm/sv = 2R^2\Delta\rho\,g/9\eta = 2 \times 2.5 \times 10^{-9} \times 1800 \times 9.81/ 0.18 = 0.49\,\mathrm{mm}/\mathrm{s}; Re=2×103\mathrm{Re} = 2 \times 10^{-3}; 41min41\,\mathrm{min} om 1.2m1.2\,\mathrm{m} 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×1011×1000×9.81/1.62×104=3.0mm/sv = 2 \times 2.5 \times 10^{-11} \times 1000 \times 9.81/1.62 \times 10^{-4} = 3.0\,\mathrm{mm}/\mathrm{s}; Re=2×103\mathrm{Re} = 2 \times 10^{-3}; 1km1\,\mathrm{km} in 3.3×1053.3 \times 10^5 s, vier dagen: elke opwaartse stroom boven enkele millimeters per seconde houdt de druppel vast.

21. 1.2m/s1.2\,\mathrm{m}/\mathrm{s}, Re=16\mathrm{Re} = 16: Stokes overschat met ongeveer een factor twee (in werkelijkheid 0.7m/s0.7\,\mathrm{m}/\mathrm{s}).

22. Stokes zou 120m/s120\,\mathrm{m}/\mathrm{s} en Re104\mathrm{Re} \sim 10^4 geven: absurd. Kwadratische weerstand, 43πR3ρwg=12ρv2πR2Cx\tfrac43\pi R^3\rho_wg = \tfrac12\rho v^2\pi R^2C_x: v=8Rρwg/3ρCx=6.6m/sv = \sqrt{8R \rho_wg/3\rho C_x} = 6.6\,\mathrm{m}/\mathrm{s} (Re=900\mathrm{Re} = 900, consistent); 150s150\,\mathrm{s} vanaf 1km1\,\mathrm{km}.

23. 1cm/s1\,\mathrm{cm}/\mathrm{s} >> 3mm/s3\,\mathrm{mm}/\mathrm{s}: de druppel stijgt; de motregendruppel heeft ongeveer 1m/s1\,\mathrm{m}/\mathrm{s} nodig, wat convectieve wolken leveren.

24. 2γ/R2\gamma/R: 29kPa29\,\mathrm{kPa}, 1.4kPa1.4\,\mathrm{kPa}, 144Pa144\,\mathrm{Pa}; dynamische drukken 12ρv2\tfrac12\rho v^2: 5×1065 \times 10^{-6}, 0.30.3 en 26Pa26\,\mathrm{Pa}. De verhouding bereikt 0.20.2 bij 2mm2\,\mathrm{mm}; bij 5mm5\,\mathrm{mm} (v9m/sv \approx 9\,\mathrm{m}/\mathrm{s}, 12ρv2=49Pa\tfrac12\rho v^2 = 49\,\mathrm{Pa} tegen 58Pa58\,\mathrm{Pa}) plat de aerodynamische druk de druppel af en breekt hem: geen regendruppel wordt groter dan ongeveer 5mm5\,\mathrm{mm}.

25. Aorta 10310^3 (laminair, dicht bij de overgang); haarvat 10310^{-3} (Poiseuille); arteriole 10210^{-2} (Poiseuille); warme pijpleiding 10510^5 (turbulent, λ(Re)\lambda(\mathrm{Re})); koude pijpleiding 10310^3 (Poiseuille); zand in olie 10310^{-3} (Stokes); wolkendruppel 10310^{-3} (Stokes); regendruppel 10310^3 (kwadratische weerstand, Cx0.5C_x \approx 0.5).

Begrippen gedefinieerd in dit hoofdstuk

Bekijk alle 393 begrippen in de begrippenlijst