Děkujeme za návštěvu webu Nature.com. Verze prohlížeče, kterou používáte, má omezenou podporu pro CSS. Pro dosažení nejlepšího zážitku doporučujeme používat aktualizovaný prohlížeč (nebo vypnout režim kompatibility v prohlížeči Internet Explorer). Mezitím budeme web zobrazovat bez stylů a JavaScriptu, abychom zajistili jeho nepřetržitou podporu.
Experimenty byly provedeny v obdélníkovém kanálu blokovaném příčnými liniemi čtyř nakloněných válcových tyčí. Tlak na středním povrchu tyče a pokles tlaku v kanálu byly měřeny změnou úhlu sklonu tyče. Byly testovány tři sestavy tyčí s různým průměrem. Výsledky měření jsou analyzovány pomocí principu zachování hybnosti a semiempirických úvah. Je generováno několik invariantních sad bezrozměrných parametrů, které vztahují tlak v kritických místech systému k charakteristickým rozměrům tyče. Bylo zjištěno, že princip nezávislosti platí pro většinu Eulerových čísel charakterizujících tlak v různých místech, tj. pokud je tlak bezrozměrný při použití projekce vstupní rychlosti kolmé k tyči, je sada nezávislá na úhlu sklonu. Výsledná semiempirická korelace může být použita pro návrh podobných hydraulických systémů.
Mnoho zařízení pro přenos tepla a hmoty se skládá ze sady modulů, kanálů nebo buněk, kterými procházejí tekutiny ve více či méně složitých vnitřních strukturách, jako jsou tyče, nárazníky, vložky atd. V poslední době se obnovil zájem o lepší pochopení mechanismů spojujících vnitřní rozložení tlaku a síly působící na složité vnitřní prvky s celkovým poklesem tlaku v modulu. Tento zájem byl mimo jiné podpořen inovacemi v materiálové vědě, rozšířením výpočetních možností pro numerické simulace a rostoucí miniaturizací zařízení. Nedávné experimentální studie vnitřního rozložení tlaku a ztrát zahrnují kanály zdrsněné různě tvarovanými žebry 1, elektrochemické reakční články 2, kapilární zúžení 3 a materiály mřížkových rámů 4.
Nejběžnějšími vnitřními strukturami jsou pravděpodobně válcové tyče procházející jednotkovými moduly, buď svazky, nebo izolované. U výměníků tepla je tato konfigurace typická na straně pláště. Pokles tlaku na straně pláště souvisí s konstrukcí výměníků tepla, jako jsou parogenerátory, kondenzátory a výparníky. V nedávné studii Wang a kol. 5 zjistili stavy proudění s opětovným připojením a společným odpojením v tandemové konfiguraci tyčí. Liu a kol. 6 měřili pokles tlaku v obdélníkových kanálech s vestavěnými dvojitými svazky trubek ve tvaru U s různými úhly sklonu a kalibrovali numerický model simulující svazky tyčí s porézním médiem.
Jak se očekávalo, existuje řada konfiguračních faktorů, které ovlivňují hydraulický výkon válcové řady: typ uspořádání (např. střídavé nebo řadové), relativní rozměry (např. rozteč, průměr, délka) a úhel sklonu, mimo jiné. Několik autorů se zaměřilo na nalezení bezrozměrných kritérií, která by řídila návrhy a zachytila ​​kombinované účinky geometrických parametrů. V nedávné experimentální studii Kim a kol. 7 navrhli efektivní model pórovitosti s využitím délky jednotkové buňky jako řídicího parametru, s použitím tandemových a střídavých polí a Reynoldsových čísel mezi 103 a 104. Snarski 8 studoval, jak se výkonové spektrum z akcelerometrů a hydrofonů připojených k válci ve vodním tunelu mění se sklonem směru proudění. Marino a kol. 9 studovali rozložení tlaku na stěně kolem válcové tyče v proudění vzduchu s otočením. Mityakov a kol. 10 vynesli rychlostní pole po otočení válce pomocí stereo PIV. Alam a kol. 11 provedli komplexní studii tandemových válců se zaměřením na vliv Reynoldsova čísla a geometrického poměru na uvolňování vírů. Byli schopni identifikovat pět stavů, a to blokování, přerušované blokování, žádné blokování, subharmonické blokování a stav opětovného připojení smykové vrstvy. Nedávné numerické studie poukázaly na tvorbu vírových struktur v proudění skrz válce s omezeným stáčením.
Obecně se očekává, že hydraulický výkon jednotkové buňky bude záviset na konfiguraci a geometrii vnitřní struktury, obvykle kvantifikované empirickými korelacemi specifických experimentálních měření. V mnoha zařízeních složených z periodických složek se proudění v každé buňce opakuje, a proto lze informace týkající se reprezentativních buněk použít k vyjádření celkového hydraulického chování struktury pomocí víceměřítkových modelů. V těchto symetrických případech lze často snížit stupeň specifičnosti, s nímž se uplatňují obecné principy zachování. Typickým příkladem je rovnice výtoku pro clonu 15. Ve zvláštním případě šikmých tyčí, ať už v uzavřeném nebo otevřeném proudění, je zajímavým kritériem často uváděným v literatuře a používaným konstruktéry dominantní hydraulická veličina (např. pokles tlaku, síla, frekvence uvolnění vírů atd.), se kterou se má složka proudění dotýkat kolmá k ose válce. Toto se často označuje jako princip nezávislosti a předpokládá, že dynamika proudění je primárně řízena normálovou složkou přítoku a že vliv axiální složky zarovnané s osou válce je zanedbatelný. Ačkoli v literatuře neexistuje shoda ohledně rozsahu platnosti tohoto kritéria, v mnoha případech poskytuje užitečné odhady v rámci experimentálních podmínek. nejistoty typické pro empirické korelace. Mezi nedávné studie o platnosti nezávislého principu patří vibrace vyvolané víry16 a jednofázový a dvoufázový průměrovaný odpor417.
V předkládané práci jsou prezentovány výsledky studie vnitřního tlaku a tlakového spádu v kanálu s příčnou linií čtyř nakloněných válcových tyčí. Byly změřeny tři sestavy tyčí s různými průměry, přičemž se mění úhel sklonu. Celkovým cílem je zkoumat mechanismus, kterým je rozložení tlaku na povrchu tyče spojeno s celkovým tlakovým spádem v kanálu. Experimentální data jsou analyzována s použitím Bernoulliho rovnice a principu zachování hybnosti za účelem vyhodnocení platnosti principu nezávislosti. Nakonec jsou generovány bezrozměrné semiempirické korelace, které lze použít k návrhu podobných hydraulických zařízení.
Experimentální uspořádání se skládalo z obdélníkové testovací sekce, do které byl proud vzduchu zajišťován axiálním ventilátorem. Testovací sekce obsahovala jednotku sestávající ze dvou rovnoběžných centrálních tyčí a dvou polovičních tyčí zapuštěných do stěn kanálu, jak je znázorněno na obr. 1e, všechny o stejném průměru. Obrázky 1a–e znázorňují podrobnou geometrii a rozměry každé části experimentálního uspořádání. Obrázek 3 znázorňuje uspořádání procesu.
a Vstupní řez (délka v mm). Vytvořeno pomocí Openscad 2021.01, openscad.org. Hlavní testovací řez (délka v mm). Vytvořeno pomocí Openscad 2021.01, openscad.org c Průřez hlavního testovacího řezu (délka v mm). Vytvořeno pomocí Openscad 2021.01, openscad.org d Exportovaný řez (délka v mm). Vytvořeno pomocí Openscad 2021.01, rozložený pohled na testovací řez z openscad.org e. Vytvořeno pomocí Openscad 2021.01, openscad.org.
Byly testovány tři sady tyčí různých průměrů. Tabulka 1 uvádí geometrické vlastnosti každého případu. Tyče jsou upevněny na úhloměru tak, aby se jejich úhel vzhledem ke směru proudění mohl měnit mezi 90° a 30° (obrázky 1b a 3). Všechny tyče jsou vyrobeny z nerezové oceli a jsou vycentrovány tak, aby mezi nimi byla zachována stejná vzdálenost mezer. Vzájemná poloha tyčí je fixována dvěma distančními podložkami umístěnými vně testované sekce.
Vstupní průtok testovací sekce byl měřen kalibrovanou Venturiho trubicí, jak je znázorněno na obrázku 2, a monitorován pomocí DP Cell Honeywell SCX. Teplota kapaliny na výstupu ze testovací sekce byla měřena teploměrem PT100 a udržována na 45±1 °C. Aby bylo zajištěno rovinné rozložení rychlosti a snížena úroveň turbulence na vstupu do kanálu, je vstupní proud vody protlačován třemi kovovými síty. Mezi posledním sítem a tyčí byla použita usazovací vzdálenost přibližně 4 hydraulických průměrů a délka výstupu byla 11 hydraulických průměrů.
Schéma zapojení Venturiho trubice používané k měření rychlosti proudění na vstupu (délka v milimetrech). Vytvořeno pomocí Openscad 2021.01, openscad.org.
Sledujte tlak na jedné z ploch středové tyče pomocí tlakového odběru o průměru 0,5 mm ve střední rovině testované sekce. Průměr odběru odpovídá úhlovému rozpětí 5°, proto je úhlová přesnost přibližně 2°. Sledovanou tyč lze otáčet kolem své osy, jak je znázorněno na obrázku 3. Rozdíl mezi tlakem na povrchu tyče a tlakem na vstupu do testované sekce se měří diferenciálním DP Cell řady Honeywell SCX. Tento tlakový rozdíl se měří pro každé uspořádání tyčí, přičemž se mění rychlost proudění, úhel sklonu (α) a azimutální úhel (θ).
Nastavení proudění. Stěny kanálu jsou zobrazeny šedě. Proud proudí zleva doprava a je blokován tyčí. Všimněte si, že pohled „A“ je kolmý k ose tyče. Vnější tyče jsou částečně zapuštěny do bočních stěn kanálu. K měření úhlu sklonu (alfa) se používá úhloměr. Vytvořeno pomocí Openscad 2021.01, openscad.org.
Účelem experimentu je změřit a interpretovat tlakový spád mezi vstupy kanálu a tlak na povrchu středové tyče, \(\theta\) a \(\alfa\) pro různé azimuty a sklony. Pro shrnutí výsledků bude tlakový spád vyjádřen v bezrozměrné formě jako Eulerovo číslo:
kde ρ je hustota kapaliny, u_{i} je střední vstupní rychlost, p_{i} je vstupní tlak a p_{w} je tlak v daném bodě na stěně tyče. Vstupní rychlost je pevně stanovena ve třech různých rozsazích určených otevřením vstupního ventilu. Výsledné rychlosti se pohybují od 6 do 10 m/s, což odpovídá Reynoldsovu číslu kanálu, Re\equiv u_{i}H/nu (kde H je výška kanálu a kinematická viskozita) mezi 40 000 a 67 000. Reynoldsovo číslo tyče (Re\equiv u_{i}d/nu) se pohybuje od 2 500 do 6 500. Intenzita turbulence odhadnutá relativní směrodatnou odchylkou signálů zaznamenaných ve Venturiho trubici je v průměru 5 %.
Obrázek 4 ukazuje korelaci \({Eu}_{w}\) s azimutálním úhlem \(\theta \), parametrizovaným třemi úhly sklonu, \(\alpha \) = 30°, 50° a 70°. Měření jsou rozdělena do tří grafů podle průměru tyče. Je vidět, že v rámci experimentální nejistoty jsou získaná Eulerova čísla nezávislá na průtoku. Obecná závislost na θ sleduje obvyklý trend tlaku na stěně kolem obvodu kruhové překážky. Při úhlech směřujících do proudění, tj. θ od 0 do 90°, tlak na stěně tyče klesá a dosahuje minima při 90°, což odpovídá mezeře mezi tyčemi, kde je rychlost největší v důsledku omezení plochy proudění. Následně dochází k obnovení tlaku θ od 90° do 100°, po kterém tlak zůstává rovnoměrný v důsledku oddělení zadní mezní vrstvy stěny tyče. Je třeba poznamenat, že nedochází ke změně úhlu minimálního tlaku, což naznačuje možné poruchy ze sousedních smykových vrstev, jako například Coandovy efekty, jsou sekundární.
Změna Eulerova čísla stěny kolem tyče pro různé úhly sklonu a průměry tyčí. Vytvořeno pomocí Gnuplotu 5.4, www.gnuplot.info.
V následující části analyzujeme výsledky na základě předpokladu, že Eulerova čísla lze odhadnout pouze pomocí geometrických parametrů, tj. poměrů délek prvku d/g a d/H (kde H je výška kanálu) a sklonu alfa. Oblíbené praktické pravidlo říká, že strukturální síla tekutiny na stáčivou tyč je určena projekcí vstupní rychlosti kolmé k ose tyče, u_{n} = u_{i} sin alfa. Tomu se někdy říká princip nezávislosti. Jedním z cílů následující analýzy je zkoumat, zda se tento princip vztahuje na náš případ, kdy jsou proudění a překážky omezeny na uzavřené kanály.
Uvažujme tlak měřený na přední straně mezilehlé plochy tyče, tj. θ = 0. Podle Bernoulliho rovnice tlak v této poloze (p)_{o} splňuje:
kde \({u}_{o}\) je rychlost tekutiny v blízkosti stěny tyče při θ = 0 a předpokládáme relativně malé nevratné ztráty. Všimněte si, že dynamický tlak je nezávislý na členu kinetické energie. Pokud je \({u}_{o}\) prázdné (tj. stagnující stav), Eulerova čísla by měla být sjednocena. Na obrázku 4 je však vidět, že při \(\theta = 0\) je výsledná \({Eu}_{w}\) blízká této hodnotě, ale ne přesně rovna, zejména pro větší úhly sklonu. To naznačuje, že rychlost na povrchu tyče při \(\theta = 0\) nemizí, což může být potlačeno vychýlením proudových čar směrem nahoru, vytvořeným nakloněním tyče. Protože proudění je omezeno na horní a dolní část testované sekce, mělo by toto vychýlení vytvořit sekundární recirkulaci, čímž se zvýší axiální rychlost dole a sníží rychlost nahoře. Za předpokladu, že velikost výše uvedeného vychýlení je projekcí vstupní rychlosti na hřídel (tj. \({u}_{i}\mathrm{cos}\alpha \)), odpovídající výsledek Eulerova čísla je:
Obrázek 5 porovnává rovnice.(3) Ukazuje dobrou shodu s odpovídajícími experimentálními daty. Průměrná odchylka byla 25 % a míra spolehlivosti 95 %. Všimněte si, že rovnice.(3) V souladu s principem nezávislosti. Podobně obrázek 6 ukazuje, že Eulerovo číslo odpovídá tlaku na zadní ploše tyče, \({p}_{180}\), a na výstupu z testovacího segmentu, \({p}_{e}\), také sleduje trend úměrný \({\mathrm{sin}}^{2}\alpha \). V obou případech však koeficient závisí na průměru tyče, což je rozumné, protože ten určuje oblast omezení. Tato vlastnost je podobná poklesu tlaku v cloně, kde je průtokový kanál částečně zmenšen na určitých místech. V této testovací části hraje roli otvoru mezera mezi tyčemi. V tomto případě tlak v místě škrcení podstatně klesá a částečně se obnovuje, když se rozšiřuje zpět. Omezení uvažujeme jako blokádu kolmou k tyči ose, tlaková ztráta mezi přední a zadní částí tyče může být zapsána jako 18:
kde \({c}_{d}\) je koeficient odporu vysvětlující zotavení parciálního tlaku mezi θ = 90° a θ = 180° a \({A}_{m}\) a \({A}_{f}\) je minimální volný průřez na jednotku délky kolmý k ose tyče a jeho vztah k průměru tyče je \({A}_{f}/{A}_{m}=\ ​​(g+d\vpravo)/g\). Odpovídající Eulerova čísla jsou:
Eulerovo číslo zdi v bodě \(\theta = 0\) jako funkce dipu. Tato křivka odpovídá rovnici (3). Vytvořeno pomocí Gnuplot 5.4, www.gnuplot.info.
Eulerovo číslo se mění s dipem v \(\theta =18{0}^{o}\) (plné znaménko) a exitu (prázdné znaménko). Tyto křivky odpovídají principu nezávislosti, tj. \(Eu\propto {\mathrm{sin}}^{2}\alpha \). Vytvořeno pomocí Gnuplotu 5.4, www.gnuplot.info.
Obrázek 7 ukazuje závislost \({Eu}_{0-180}/{\mathrm{sin}}^{2}\alpha \) na \(d/g\), což ukazuje extrémně dobrou konzistenci.(5). Získaný koeficient odporu je \({c}_{d}=1,28\pm 0,02\) s hladinou spolehlivosti 67 %. Stejný graf také ukazuje, že celkový pokles tlaku mezi vstupem a výstupem z testované sekce sleduje podobný trend, ale s různými koeficienty, které zohledňují zotavení tlaku v zadním prostoru mezi příčkou a výstupem kanálu. Odpovídající koeficient odporu je \({c}_{d}=1,00\pm 0,05\) s hladinou spolehlivosti 67 %.
Součinitel odporu souvisí s poklesem tlaku \(d/g\) před a za tyčí \(\left({Eu}_{0-180}\right)\) a s celkovým poklesem tlaku mezi vstupem a výstupem kanálu. Šedá oblast představuje 67% pásmo spolehlivosti pro korelaci. Vytvořeno pomocí Gnuplot 5.4, www.gnuplot.info.
Minimální tlak \({p}_{90}\) na povrchu tyče při θ = 90° vyžaduje speciální zacházení. Podle Bernoulliho rovnice, podél proudové čáry skrz mezeru mezi tyčemi, tlak ve středu \({p}_{g}\) a rychlost \({u}_{g}\) v mezeře mezi tyčemi (která se shoduje se středem kanálu) souvisí s následujícími faktory:
Tlak \({p}_{g}\) lze vztáhnout k tlaku na povrchu tyče při θ = 90° integrací rozložení tlaku v mezeře oddělující centrální tyč mezi středem a stěnou (viz obrázek 8). Rovnováha sil dává 19:
kde y je souřadnice kolmá k povrchu tyče ze středu mezery mezi centrálními tyčemi a K je zakřivení proudové čáry v poloze y. Pro analytické vyhodnocení tlaku na povrchu tyče předpokládáme, že u_{g} je rovnoměrné a K(y) je lineární. Tyto předpoklady byly ověřeny numerickými výpočty. Na stěně tyče je zakřivení určeno elipsovitým průřezem tyče pod úhlem alfa, tj. K(g/2) = (2/d)^{2} alfa (viz obrázek 8). Pak s ohledem na zakřivení proudnice mizící v bodě y=0 v důsledku symetrie je zakřivení v univerzální souřadnici y dáno vztahem:
Průřezový pohled na prvek, zepředu (vlevo) a shora (dole). Vytvořeno v programu Microsoft Word 2019.
Na druhou stranu, z důvodu zachování hmotnosti, je průměrná rychlost v rovině kolmé k proudění v místě měření (\langle {u}_{g}\rangle \) vztažena k vstupní rychlosti:
kde \({A}_{i}\) je průřezová plocha proudění na vstupu do kanálu a \({A}_{g}\) je průřezová plocha proudění v místě měření (viz obr. 8) podle následujícího vzorce:
Všimněte si, že \({u}_{g}\) se nerovná \(\langle {u}_{g}\rangle \). Obrázek 9 ve skutečnosti znázorňuje poměr rychlostí \({u}_{g}/\langle {u}_{g}\rangle \), vypočítaný pomocí rovnice (10)–(14), vynesený do grafu podle poměru \(d/g\). Navzdory určité diskrétnosti lze identifikovat trend, který je aproximován polynomem druhého řádu:
Poměr maximální\({u}_{g}\) a průměrné\(\langle {u}_{g}\rangle \) rychlosti průřezu středu kanálu\(.\). Plné a přerušované křivky odpovídají rovnicím.(5) a rozsahu variace odpovídajících koeficientů\(\pm 25\%\). Vytvořeno pomocí Gnuplot 5.4, www.gnuplot.info.
Obrázek 10 porovnává \({Eu}_{90}\) s experimentálními výsledky rovnice.(16). Průměrná relativní odchylka byla 25 % a míra spolehlivosti 95 %.
Wallovo Eulerovo číslo v bodě \(\theta ={90}^{o}\). Tato křivka odpovídá rovnici. (16). Vytvořeno pomocí Gnuplot 5.4, www.gnuplot.info.
Výslednou sílu \({f}_{n}\) působící na centrální tyč kolmo k její ose lze vypočítat integrací tlaku na povrchu tyče takto:
kde první koeficient je délka tyče v kanálu a integrace se provádí mezi 0 a 2π.
Projekce \({f}_{n}\) ve směru proudění vody by měla odpovídat tlaku mezi vstupem a výstupem z kanálu, pokud tření není rovnoběžné s tyčí a menší kvůli neúplnému rozvinutí pozdější části. Tok hybnosti je proto nevyvážený.
Obrázek 11 ukazuje graf rovnic.(20) vykazoval dobrou shodu pro všechny experimentální podmínky. Vpravo je však mírná 8% odchylka, kterou lze připsat a použít jako odhad nerovnováhy hybnosti mezi vstupem a výstupem kanálu.
Rovnováha výkonu kanálů. Čára odpovídá rovnici (20). Pearsonův korelační koeficient byl 0,97. Vytvořeno pomocí Gnuplotu 5.4, www.gnuplot.info.
Změnou úhlu sklonu tyče byl měřen tlak na stěně povrchu tyče a pokles tlaku v kanálu s příčnými liniemi čtyř nakloněných válcových tyčí. Byly testovány tři sestavy tyčí s různým průměrem. V testovaném rozsahu Reynoldsových čísel mezi 2500 a 6500 je Eulerovo číslo nezávislé na průtoku. Tlak na centrálním povrchu tyče sleduje obvyklý trend pozorovaný u válců, přičemž je maximální vpředu a minimální v boční mezeře mezi tyčemi a v zadní části se zotavuje v důsledku oddělení mezní vrstvy.
Experimentální data jsou analyzována s využitím zákonu zachování hybnosti a semiempirických výpočtů za účelem nalezení invariantních bezrozměrných čísel, která vztahují Eulerova čísla k charakteristickým rozměrům kanálů a tyčí. Všechny geometrické vlastnosti blokování jsou plně reprezentovány poměrem mezi průměrem tyče a mezerou mezi tyčemi (laterálně) a výškou kanálu (vertikálně).
Bylo zjištěno, že princip nezávislosti platí pro většinu Eulerových čísel charakterizujících tlak v různých místech, tj. pokud je tlak bezrozměrný pomocí projekce vstupní rychlosti kolmé k tyči, je množina nezávislá na úhlu sklonu. Kromě toho tento rys souvisí s hmotností a hybností proudění. Rovnice zachování jsou konzistentní a podporují výše uvedený empirický princip. Pouze tlak na povrchu tyče v mezeře mezi tyčemi se od tohoto principu mírně odchyluje. Jsou generovány bezrozměrné semiempirické korelace, které lze použít k návrhu podobných hydraulických zařízení. Tento klasický přístup je v souladu s nedávno publikovanými podobnými aplikacemi Bernoulliho rovnice v hydraulice a hemodynamice20,21,22,23,24.
Obzvláště zajímavý výsledek vyplývá z analýzy tlakového spádu mezi vstupem a výstupem zkušební sekce. V rámci experimentální nejistoty je výsledný koeficient odporu roven jedné, což naznačuje existenci následujících invariantních parametrů:
Všimněte si velikosti \(\left(d/g+2\right)d/g\) ve jmenovateli rovnice.(23) je velikost v závorkách v rovnici.(4), jinak ji lze vypočítat s minimálním a volným průřezem kolmým k tyči, \({A}_{m}\) a \({A}_{f}\). To naznačuje, že Reynoldsova čísla se předpokládají v rozsahu aktuální studie (40 000–67 000 pro kanály a 2 500–6 500 pro tyče). Je důležité poznamenat, že pokud je uvnitř kanálu teplotní rozdíl, může to ovlivnit hustotu kapaliny. V tomto případě lze relativní změnu Eulerova čísla odhadnout vynásobením koeficientu tepelné roztažnosti maximálním očekávaným teplotním rozdílem.
Ruck, S., Köhler, S., Schlindwein, G. a Arbeiter, F. Měření přenosu tepla a tlakové ztráty v kanálu zdrsněném stěnou s různě tvarovanými žebry. expert. Heat Transfer 31, 334–354 (2017).
Wu, L., Arenas, L., Graves, J. a Walsh, F. Charakterizace průtokové cely: vizualizace proudění, pokles tlaku a transport hmoty ve dvourozměrných elektrodách v obdélníkových kanálech. J. Electrochemistry. Socialist Party. 167, 043505 (2020).
Liu, S., Dou, X., Zeng, Q. a Liu, J. Klíčové parametry Jaminova jevu v kapilárách se zúženými průřezy. J. Gasoline.science.Britain.196, 107635 (2021).