Hatur nuhun pikeun ngadatangan Nature.com.Vérsi browser anu anjeun anggo gaduh dukungan terbatas pikeun CSS. Kanggo pangalaman pangsaéna, kami nyarankeun yén anjeun nganggo browser anu diropéa (atanapi mareuman modeu kasaluyuan dina Internet Explorer).
Percobaan anu dipigawé dina saluran rectangular diblokir ku garis transverse opat rod cylindrical condong.Tekenan dina beungeut rod puseur jeung serelek tekanan sakuliah saluran diukur ku varying sudut inclination rod urang.Tilu rakitan rod diaméterna béda anu diuji.Hasil pangukuran dianalisis ngagunakeun prinsip konservasi moméntum sarta parameter semi-émpiris anu dibangkitkeun revarian diménsi. tekanan dina lokasi kritis sistem ka dimensi karakteristik rod.Prinsip kamerdikaan kapanggih nyekel keur paling angka Euler characterizing tekanan dina lokasi béda, ie lamun tekanan anu dimensionless ngagunakeun proyéksi tina laju inlet normal pikeun rod, susunan nu bebas tina sudut dip.Korelasi semi-émpiris anu dihasilkeun bisa dipaké pikeun Desain hydraulics sarupa.
Seueur alat transfer panas sareng massa diwangun ku sakumpulan modul, saluran atanapi sél anu ngalangkungan cairan dina struktur internal anu langkung atanapi kirang kompleks sapertos rod, panyangga, sisipan, jsb. Langkung-langkung ayeuna, aya minat anu énggal-énggal pikeun ngartos langkung saé ngeunaan mékanisme anu nyambungkeun distribusi tekanan internal sareng gaya dina internal kompleks kana serelek tekanan sakabéh modul. al simulasi, sarta ngaronjatna miniaturization alat.Ulikan eksperimen panganyarna ngeunaan sebaran internal tekanan sarta karugian kaasup saluran roughened ku rupa-rupa tulang rusuk ngawangun 1, sél réaktor éléktrokimia 2, constriction kapilér 3 jeung bahan pigura kisi 4.
Struktur internal paling umum nyaéta rod cylindrical arguably ngaliwatan modul Unit, boh dibuntel atanapi isolated.In exchangers panas, konfigurasi ieu has dina cangkang side.Shell serelek tekanan samping patali jeung desain exchangers panas kayaning generator uap, condensers na evaporators.Dina ulikan panganyarna, Wang et al.5 kapanggih reattachment jeung aliran ko-detachment nagara bagian dina konfigurasi tandem of rods.Liu et al.6 diukur serelek tekanan dina saluran rectangular kalawan diwangun-di ganda U ngawangun bundles tube kalawan sudut inclination béda jeung calibrated model numeris simulating bundles rod kalawan média porous.
Saperti nu diharapkeun, aya sababaraha faktor konfigurasi nu mangaruhan kinerja hidrolik tina bank silinder: tipe susunan (misalna staggered atanapi in-line), dimensi relatif (misalna pitch, diaméterna, panjang), sarta sudut inclination, antara séjén. Sababaraha pangarang fokus kana nyungsi kriteria dimensionless pikeun pituduh desain pikeun nangkep efek gabungan tina parameter geometric.7 diusulkeun model porosity éféktif ngagunakeun panjang sél Unit salaku parameter kontrol, ngagunakeun tandem na arrays staggered jeung angka Reynolds antara 103 jeung 104.Snarski8 diajar kumaha spéktrum kakuatan, ti accelerometers na hydrophones napel silinder dina torowongan cai, beda-beda jeung inclination tina arah aliran et al.Marino.9 diajar sebaran tekanan témbok sabudeureun rod cylindrical di yaw airflow.Mityakov et al.10 plotted widang laju sanggeus silinder yawed maké stereo PIV.Alam dkk.11 ngalaksanakeun ulikan komprehensif ngeunaan silinder tandem, fokus kana épék angka Reynolds jeung ratio geometric on vortex shedding.They éta bisa ngaidentipikasi lima kaayaan, nyaéta ngonci, ngonci intermittent, euweuh ngonci, subharmonic ngonci jeung nagara lapisan geser reattachment.Ulikan numeris panganyarna geus nunjuk kana formasi struktur vortex yawlistricted dina aliran ngaliwatan yawlinder.
Sacara umum, kinerja hidrolik sél unit diperkirakeun gumantung kana konfigurasi jeung géométri tina struktur internal, biasana diitung ku korelasi empiris tina ukuran ékspérimén husus. Dina loba alat nu diwangun ku komponén periodik, pola aliran nu terus-terusan dina unggal sél, sahingga, informasi nu patali jeung sél wawakil bisa dipaké pikeun nganyatakeun sakabéh paripolah hidrolik tina struktur ngaliwatan model multiscale, nu bisa diterapkeun dina darajat conservation umum. conto has nyaeta persamaan ngurangan pikeun hiji pelat orifice 15. Dina kasus husus tina rod condong, naha dina dipasrahkeun atawa kabuka aliran, hiji kriteria metot mindeng dicutat dina literatur jeung dipaké ku désainer nyaéta gedéna hidrolik dominan (misalna, tekanan turun, gaya, vortex shedding frékuénsi, jsb) ) pikeun kontak.) kana komponén aliran perpendikular ieu perpendikular sarta permanén. dinamika aliran didorong utamana ku komponén normal inflow sarta yén pangaruh komponén axial Blok kalawan sumbu silinder nyaeta negligible.Sanajan euweuh konsensus dina literatur dina rentang validitas kriteria ieu, dina loba kasus eta nyadiakeun estimasi mangpaat dina uncertainties ékspérimén has has tina korelasi empiris jeung empiris. sareng dua fase rata-rata drag417.
Dina karya hadir, hasil ulikan ngeunaan tekanan internal tur serelek tekanan dina saluran kalawan garis transverse opat rod cylindrical condong dibere. Ukur tilu rakitan rod kalawan diaméter béda, ngarobah sudut inclination.The Tujuan sakabéh nyaéta pikeun nalungtik mékanisme ku nu distribusi tekanan dina beungeut rod pakait jeung serelek sakabéh tekanan dina dasarna Berkualitas dina channel. pikeun evaluate validitas prinsip kamerdikaan.Ahirna, correlations semi-empiris dimensionless dihasilkeun nu bisa dipaké pikeun ngarancang alat hidrolik sarupa.
Setélan ékspérimén diwangun ku bagian tés rectangular nu narima aliran hawa disadiakeun ku kipas axial. Bagian test ngandung hiji unit diwangun ku dua rod sentral paralel jeung dua satengah rod dipasang dina tembok channel, sakumaha ditémbongkeun dina Gbr. 1e, sakabéh diaméterna sarua. Angka 1a-e nembongkeun géométri lengkep jeung dimensi unggal bagian tina setup.
a bagian Inlet (panjangna dina mm).Jieun b maké Openscad 2021.01, openscad.org.Bagian uji utama (panjangna dina mm).Dijieun kalawan Openscad 2021.01, openscad.org c Pandangan cross-sectional sahiji bagian test utama (panjangna dina mm) .Dijieun maké Openscad 2021.01, kalayan muka openscad.01. 021.01, exploded view tina bagian tés openscad.org e.Created kalawan Openscad 2021.01, openscad.org.
Tilu sét rod diaméter béda ieu tested.Table 1 daptar ciri geometri unggal case.The rod anu dipasang dina protractor a ambéh sudut maranéhna relatif ka arah aliran bisa rupa-rupa antara 90 ° jeung 30 ° (Angka 1b na 3) .All rod dijieunna tina stainless steel sarta aranjeunna dipuseurkeun pikeun ngajaga jarak gap anu sarua antara dua rodrs.
Laju aliran inlet sahiji bagian test diukur ku venturi calibrated, ditémbongkeun saperti dina Gambar 2, sarta diawaskeun maké DP Cell Honeywell SCX. Suhu cairan di outlet sahiji bagian test diukur ku thermometer PT100 tur dikawasa dina 45 ± 1 ° C. kira-kira 4 diaméter hidrolik dipaké antara layar panungtungan sarta rod, sarta panjang stop kontak éta 11 diaméter hidrolik.
Diagram skéma tina tabung Venturi dipaké pikeun ngukur laju aliran inlet (panjangna dina milimeter).Dijieun kalawan Openscad 2021.01, openscad.org.
Ngawas tekanan dina salah sahiji rupa rod puseur ku cara maké 0,5 mm tekanan ketok dina pertengahan pesawat tina test section.The diaméter ketok pakait jeung bentang sudut 5 °;kituna akurasi sudut kira 2 °. Rod diawaskeun bisa diputer ngeunaan sumbu na, ditémbongkeun saperti dina Gambar 3. Bedana antara tekanan permukaan rod jeung tekanan dina lawang ka bagian test diukur ku diferensial DP Cell Honeywell SCX series.Beda tekanan ieu diukur pikeun tiap susunan bar, varying laju aliran, sudut inclination \(azimuth) jeung \(azimuth sudut).
Setélan aliran. Tembok saluran ditémbongkeun dina warna abu. Aliran ngalir ti kénca ka katuhu sarta diblokir ku rod. Catet yén pintonan "A" jejeg sumbu rod. Batang luar semi-dipasang dina tembok channel gurat. A protractor dipaké pikeun ngukur sudut inclination \(\ alpha \).Dijieun kalawan Openscad.org 1,02.
Tujuan tina percobaan nyaéta pikeun ngukur jeung napsirkeun serelek tekanan antara inlets channel sarta tekanan dina beungeut rod puseur, \(\ theta \) jeung \(\ alpha \) pikeun azimuths béda jeung dips. Pikeun nyimpulkeun hasil, tekanan diferensial bakal dinyatakeun dina formulir diménsi salaku angka Euler urang:
dimana \(\rho \) nyaéta dénsitas cairan, \({u}_{i}\) nyaéta laju asupan rata-rata, \({p}_{i}\) nyaéta tekanan inlet, jeung \({p }_{ w}\) nyaéta tekanan dina titik nu tangtu dina témbok rod. Laju inlet dibenerkeun dina tilu rentang béda ditangtukeun ku rentang ti 0, veloc inlet nu ditangtukeun ku rentang bukaan tina klep inlet. Angka Reynolds, \(Re\equiv {u}_{i}H/\nu \) (dimana \(H\) nyaéta jangkungna saluran, jeung \(\nu \) nyaéta viskositas kinematic) antara 40.000 jeung 67.000. Rod Reynolds number (\(Re\ equiv {u}_{i}d/\nu) rentang ti u}_{i}d/\nu. simpangan baku relatif tina sinyal kacatet dina venturi nyaeta 5% rata-rata.
angka 4 nembongkeun korelasi \({Eu}_{w}\) jeung sudut azimuth \(\ theta \), parameterized ku tilu sudut dip, \(\ alpha \) = 30 °, 50 ° jeung 70 ° .The pangukuran dibagi jadi tilu grafik nurutkeun diaméter rod. Ieu bisa ditempo yén dina laju ékspérimén nu teu pasti, dina laju nu teu katangtu. θ nuturkeun trend biasa tekanan témbok sabudeureun perimeter of a obstacle sirkular. Dina sudut aliran-nyanghareup, ie, θ ti 0 nepi ka 90 °, tekanan témbok rod nurun, ngahontal minimum dina 90 °, nu pakait jeung celah antara rod mana laju anu greatest alatan watesan aréa aliran. seragam alatan separation tina lapisan wates pungkur tina témbok rod. Catet yén euweuh parobahan dina sudut tekanan minimum, nu nunjukkeun yén gangguan mungkin tina lapisan geser meungkeut, kayaning épék Coanda, sekundér.
Variasi jumlah Euler témbok sabudeureun rod pikeun sudut inclination béda jeung diaméter rod.Dijieun kalawan Gnuplot 5.4, www.gnuplot.info.
Di handap ieu, urang nganalisis hasil dumasar kana asumsi yén angka Euler bisa diperkirakeun ukur ku parameter geometric, nyaéta babandingan panjang fitur \ (d / g \) jeung \ (d / H \) (dimana \ (H \) nyaeta jangkungna channel urang) jeung inclination \ (\ alpha \).A aturan praktis populér tina thumb nyatakeun yén gaya yawvelocity fluida dina struktur perpendikular nangtukeun gaya yawvelocity cairan. sumbu rod, \({u}_{n}={u}_{i}\mathrm {sin} \alpha \) .Ieu sok disebut prinsip kamerdikaan.Salah sahiji tujuan analisis di handap ieu pikeun nalungtik naha prinsip ieu lumaku pikeun kasus urang, dimana aliran jeung halangan anu dipasrahkeun dina saluran katutup.
Hayu urang nganggap tekanan diukur di hareup beungeut rod panengah, nyaéta θ = 0. Numutkeun persamaan Bernoulli urang, tekanan dina posisi ieu satisfies:
dimana \({u}_{o}\) nyaéta laju fluida deukeut témbok rod dina θ = 0, sarta kami anggap karugian teu bisa balik relatif leutik. Catet yén tekanan dinamis bebas dina istilah énergi kinétik. Lamun \({u}_{o}\) kosong (ie kaayaan ngeuyeumbeu), angka Euler kudu ngahiji. Sanajan kitu, dina Figure nu \ 4 bisa jadi dihijikeun. {Eu}_{w}\) deukeut tapi teu persis sarua jeung nilai ieu, utamana pikeun sudut dip nu leuwih gede. Ieu nunjukkeun yén laju dina beungeut rod teu ngaleungit dina \ (\ theta = 0 \), nu bisa diteken ku deflection ka luhur tina garis ayeuna dijieun ku rod Dengdekkeun. Kusabab aliran ieu dipasrahkeun ka luhur jeung handap dina defleksi bagian test, a decreas sekundér ieu kudu nyieun hiji deflection sekundér jeung handap. Nyanyikeun laju di luhur. Anggap yén gedéna defleksi di luhur nyaéta proyéksi laju inlet dina aci (ie \({u}_{i}\mathrm{cos}\alpha \)), hasil angka Euler anu cocog nyaéta:
Gambar 5 ngabandingkeun persamaan.(3) Ieu nunjukeun kasapukan alus jeung data ékspérimén pakait. Rata simpangan éta 25%, sarta tingkat kapercayaan éta 95%. Catet yén persamaan.(3) Luyu jeung prinsip kamerdikaan. Kitu ogé, Gambar 6 nunjukeun yen angka Euler pakait jeung tekanan dina beungeut pungkur ti rod), (\({0}) tina rod, (\({0}) tina rod, (\({0}) ex, jeung \({0}) ex. {p}_{e}\), Ogé nuturkeun trend sabanding jeung \({\ mathrm {sin}}^{2}\alpha \) .Dina dua kasus, kumaha oge, koefisien gumantung kana diaméter rod, nu lumrah saprak dimungkinkeun dina nangtukeun aréa hindered. Fitur ieu sarupa jeung serelek tekanan tina hiji plat orifice, dimana saluran aliran dicoo sawaréh dina posisi rod diréduksi ku bagian tina rod husus. .Dina hal ieu, tekanan pakait substansi di throttling sarta sawaréh recovers sakumaha eta expands backwards.Considering larangan salaku sumbatan jejeg sumbu rod, serelek tekanan antara hareup jeung pungkur rod bisa ditulis salaku 18:
dimana \({c}_{d}\) nyaéta koefisien seret anu ngajelaskeun pamulihan tekanan parsial antara θ = 90° sareng θ = 180°, sareng \({A}_{m}\) sareng \ ({A}_{f}\) mangrupikeun penampang bébas minimum per unit panjang anu jejeg sumbu rod, sareng hubunganana sareng diameter rod nyaéta {({+A} d diaméterna {\} {\A} katumbiri. )/g\).Angka-angka Euler anu pakait nyaéta:
Tembok Euler angka dina \ (\ theta = 0 \) salaku fungsi dip.kurva Ieu pakait jeung persamaan.(3).Dijieun kalawan Gnuplot 5.4, www.gnuplot.info.
Nomer Tembok Euler robah, dina \(\theta =18{0}^{o}\) (tanda pinuh) jeung kaluar (tanda kosong) kalawan dip.kurva ieu luyu jeung prinsip kamerdikaan, nyaéta \(Eu\propto {\mathrm{sin}}^{2}\alpha \).Dijieun ku Gnuplot 5.4, www.gnuplot.info.gnup.
Gambar 7 nembongkeun gumantungna tina \({Eu}_{0-180}/{\mathrm{sin}}^{2}\alpha \) dina \(d/g\), némbongkeun konsistensi Alus ekstrim.(5).Koéfisién sered anu dimeunangkeun nyaéta \({c}_{d}=1,28\pm 0,02 \) jeung tingkat kayakinan anu sarua antara 0,02 \) kalawan tingkat kapercayaan anu sarua antara 7%. sarta outlet sahiji bagian test nuturkeun trend sarupa, tapi kalawan koefisien béda nu tumut kana akun recovery tekanan dina spasi tukang antara bar jeung outlet saluran.Koefisien sered pakait nyaeta \ ({c}_{d} = 1,00 \ pm 0,05 \) kalawan tingkat kapercayaan 67%.
Koéfisién sered patali jeung \ (d / g \) serelek tekanan fore na aft rod \ (\ kénca ({Eu}_{0-180} \ katuhu) \) jeung total serelek tekanan antara saluran inlet na outlet. Wewengkon abu nyaéta 67% pita kapercayaan pikeun korelasi. Dijieun ku Gnuplot 5.4, www.gnuplot.info.gnuplot.
Tekanan minimum \({p}_{90}\) dina beungeut rod dina θ = 90° merlukeun penanganan husus.Nurutkeun persamaan Bernoulli, sapanjang garis ayeuna ngaliwatan celah antara bar, tekanan di puseur\({p}_{g}\) jeung laju\({u}_{g}\) dina celah antara koin handap jeung faktor handap:
Tekanan \({p}_{g}\) tiasa dipatalikeun sareng tekanan permukaan rod dina θ = 90 ° ku ngahijikeun distribusi tekanan dina sela anu misahkeun rod sentral antara titik tengah sareng témbok (tingali Gambar 8).Kasaimbangan kakuatan masihan 19:
dimana \ (y \) nyaéta koordinat normal kana beungeut rod ti titik puseur celah antara rod sentral, sarta \ (K \) nyaéta curvature tina garis ayeuna dina posisi \ (y \). Pikeun evaluasi analitik tekanan dina beungeut rod, urang nganggap yén \ ({u}_{g}\) geus seragam jeung \ (K \ kénca salaku linearations.A) t témbok rod, curvature ditangtukeun ku bagian elips tina rod dina sudut \ (\ alpha \), ie \ (K \ kénca (g / 2 \ katuhu) = \ kénca (2 / d \ katuhu) {\ mathrm {dosa} } ^ {2} \ alpha \) (tingali Gambar 8). Lajeng, ngeunaan curvature nu universal dina \ koordinat 0. e \(y\) dirumuskeun ku:
Fitur view cross-sectional, hareup (kénca) jeung luhur (handap).Dijieun ku Microsoft Word 2019,
Di sisi séjén, ku konservasi massa, laju rata-rata dina pesawat jejeg aliran di lokasi pangukuran \(\langle {u}_{g}\rangle \) patali jeung laju inlet:
dimana \({A}_{i}\) nyaéta aréa aliran cross-sectional dina inlet saluran jeung \({A}_{g}\) nyaéta aréa aliran cross-sectional di lokasi pangukuran (tingali Gbr. 8) masing-masing ku:
Catet yén \({u}_{g}\) teu sarua jeung \(\langle {u}_{g}\rangle \).Malah, Gambar 9 ngagambarkeun rasio laju \({u}_{g}/\langle {u}_{g}\rangle \), diitung ku persamaan.(10)–(14), plotted nurutkeun rasio \(d/g). polinomial orde dua:
Babandingan laju maksimum\({u}_{g}\) jeung rata-rata\(\langle {u}_{g}\rangle \) tina cross-section puseur saluran\(.\) Kurva padet jeung dashed pakait jeung persamaan.(5) jeung rentang variasi tina koefisien pakait\(\pm 25\%\).
Gambar 10 ngabandingkeun \({Eu}_{90}\) jeung hasil ékspérimén tina persamaan.(16).Mean simpangan relatif 25%, jeung tingkat kapercayaan 95%.
Jumlah Tembok Euler dina \(\theta ={90}^{o}\).Kurva ieu pakait jeung persamaan.(16).Dijieun ku Gnuplot 5.4, www.gnuplot.info.
Gaya net \({f}_{n}\) nu nimpah rod sentral jejeg sumbu na bisa diitung ku ngahijikeun tekanan dina beungeut rod saperti kieu:
dimana koefisien kahiji nyaéta panjang rod dina saluran, sarta integrasi dipigawé antara 0 jeung 2π.
Proyéksi \({f}_{n}\) dina arah aliran cai kedah cocog tekanan antara inlet na outlet saluran, iwal gesekan sajajar jeung rod sarta leutik alatan ngembangkeun lengkep bagian engké.Ku kituna,
angka 11 nembongkeun grafik tina persamaan.(20) némbongkeun perjangjian alus pikeun sakabéh kaayaan ékspérimén. Tapi, aya slight 8% simpangan on katuhu, nu bisa attributed tur dipaké salaku estimasi tina saimbangna moméntum antara saluran inlet na outlet.
Kakuatan saluran kasaimbangan.Garis pakait jeung persamaan.(20).Koefisien korelasi Pearson éta 0.97.Dijieun kalawan Gnuplot 5.4, www.gnuplot.info.
Varying sudut inclination of rod, tekanan dina témbok beungeut rod jeung serelek tekanan dina saluran jeung garis transverse tina opat rod cylindrical condong éta measured.Three rod diaméterna béda rakitan anu tested.In rentang angka Reynolds diuji, antara 2500 jeung 6500, jumlah Euler nu rod cylindrical nu diturutan bebas tina laju aliran séntral, nu dititénan dina laju aliran sentral maksimum. hareup jeung minimum dina celah gurat antara rod, recovering di bagian tukang alatan separation lapisan wates.
Data ékspérimén dianalisis ngagunakeun pertimbangan konservasi moméntum sarta evaluasi semi-empiris pikeun manggihan invarian wilangan dimensionless nu patali angka Euler kana dimensi karakteristik saluran jeung rods.All fitur geometrical blocking pinuh digambarkeun ku rasio antara diameter rod jeung celah antara rod (laterally) jeung jangkungna channel (vertikal).
Prinsip kamerdikaan kapanggih nyekel keur paling angka Euler characterizing tekanan dina lokasi béda, ie lamun tekanan anu dimensionless ngagunakeun proyéksi laju inlet normal rod, susunan nu bebas tina sudut dip.Sajaba ti éta, fitur nu patali jeung massa sarta moméntum aliran The persamaan konservasi anu konsisten tur ngarojong prinsip empiris di luhur. Ngan tekanan permukaan rod dina celah antara rod rada nyimpang tina prinsip ieu.Dimensionless semi-empiris correlations dihasilkeun nu bisa dipaké pikeun ngarancang sarupa devices.The hidrolik pendekatan klasik ieu konsisten jeung panganyarna dilaporkeun aplikasi sarupa, Bernoulli2 hidrolik2, Berequation2, hidrodinamika 2. 24.
Hiji hasil utamana metot batang tina analisis turunna tekanan antara inlet jeung outlet tina bagian test.Dina kateupastian ékspérimén, hasilna koefisien sered sarua persatuan, nu nunjukkeun ayana parameter invarian handap:
Catet ukuran \(\left(d/g+2\right)d/g\) dina pangbagi persamaan.(23) nyaéta gedéna dina tanda kurung dina persamaan.(4), disebutkeun eta bisa diitung kalawan minimum jeung bebas cross-sectional jejeg rod, \({A}_{m}}\) jeung \({A}_{m}}\) jeung \}\{A} s nu tetep dina rentang Reyd nu. ulikan ayeuna (40.000-67.000 pikeun saluran jeung 2500-6500 pikeun rod). Kadé dicatet yén lamun aya béda hawa jero saluran, éta bisa mangaruhan dénsitas cairan. Dina hal ieu, parobahan relatif dina jumlah Euler bisa diperkirakeun ku cara ngalikeun koefisien ékspansi termal ku bédana hawa ekspektasi maksimum.
Ruck, S., Köhler, S., Schlindwein, G., sarta Arbeiter, F. mindahkeun panas sarta pangukuran serelek tekanan dina saluran roughened ku tulang rusuk béda ngawangun dina wall.expert.Heat Transfer 31, 334-354 (2017).
Wu, L., Arenas, L., Graves, J., sarta Walsh, F. Aliran characterization sél: visualisasi aliran, serelek tekanan, sarta angkutan massa dina éléktroda dua diménsi dina channels rectangular.J.Éléktrokimia.Partai Sosialis.167, 043505 (2020).
Liu, S., Dou, X., Zeng, Q. & Liu, J. Parameter konci pangaruh Jamin dina kapilér jeung constricted cross-sections.J.Bensin.science.Britain.196, 107635 (2021).