Numerical Modelling of the Effects of Biofouling on Ship Resistance and Propulsion Characteristics

Autor: Farkas, Andrea
Přispěvatelé: Parunov, Joško, Kožuh, Stjepan, Bauer, Branko, Duić, Neven, Jokić, Andrej, Landek, Darko, Lisjak, Dragutin, Lulić, Zoran, Majetić, Dubravko, Matijević, Božidar, Runje, Biserka, Sorić, Jurica, Terze, Zdravko (ur.)., Degiuli, Nastia
Jazyk: angličtina
Rok vydání: 2018
Předmět:
Popis: The effect of biofouling on the hydrodynamic characteristics of ship resistance and propulsion in calm water is very important from both an economic and environmental point of view. Hydrodynamic performance of a ship is disrupted because of the presence of biofouling organisms, which results in increased fuel consumption, ship speed reduction and increased emission of harmful gases. Presently, there is no comprehensive procedure, which could reliably predict the effect of biofouling on the ship hydrodynamic characteristics. Consequently, International Towing Tank Conference (ITTC) has advised scientists to present new formulae or methods based on the experimental data to determine the effect of biofouling on the ship resistance and propulsion characteristics. Since biofouling depends on many parameters and it is very difficult to predict how long will antifouling coatings prevent fouling of a ship, the proposed research is focused on the effects of predetermined surface conditions on the ship hydrodynamic characteristics. Biofouling can be classified into the soft, hard and composite fouling. In this thesis the effects of biofilm and hard fouling on ship resistance, propeller open water and ship self-propulsion characteristics are investigated. Within the proposed research commercial software package is used. The mathematical model is based on the averaged continuity equation and Reynolds Averaged Navier-Stokes (RANS) equations. Governing equations are discretised utilizing Finite Volume Method (FVM). After the analysis of several turbulence models for the closure of set of equations and their influence on the obtained ship hydrodynamic characteristics has been performed, k − Shear Stress Transport turbulence model is selected as a compromise between accuracy and computational time. Volume of Fluid method is utilized for tracking and locating the free surface. The effects of biofouling are modelled through the implementation of roughness function model within the wall function of Computational Fluid Dynamics (CFD) solver. The research is based on the wall similarity hypothesis which claims that roughness effects are limited to inner layer of turbulent boundary layer. The validity of numerical procedure is examined through verification and validation of the obtained results. The validation of the obtained numerical results for the smooth surface condition is carried out by comparison with the extrapolated towing tank results and other numerical studies available in the literature. The validation of numerical drag characterization study is performed by comparison of the obtained numerical results and experimental ones. Also, the obtained numerical results in terms of the increase in frictional resistance for a flat plate having the same length as a ship are compared with the ones obtained using the Granville similarity law scaling method. Thereafter, the applicability of the proposed CFD approach is demonstrated on the example of three full-scale merchant ships. Also, the newly proposed performance prediction method for fouled surfaces is presented, which can account for fouling effects on the ship performance. The applicability of this method is demonstrated for fouling conditions with lower fouling rates. Thus, a robust and rapid assessment of the effects of biofouling on the ship hydrodynamic characteristics in calm water is enabled. Uvod Obraštanje je nakupljanje mikroorganizama, biljaka, algi ili životinja na površinama koje su u kontaktu s vodom te predstavlja rastući problem s ekonomskog, ali i ekološkog stajališta. Iako je obraštanje prisutno u mnogim različitim područjima, pomorska industrija je najviše pogođena prisutnošću i rastom organizama. Prisutnost raznih organizama na trupu broda uzrokuje povećanje hrapavosti trupa broda. Zbog toga se otpor trenja povećava, što uzrokuje povećanje potrošnje goriva, smanjene brzine broda i povećanje opterećenja motora. Pokazano je kako su svi troškovi vezani uz razvoj, nabavu i primjenu, kao i relativno skupo tehničko ili plansko rješenje za antivegetativnu zaštitu brodova, ekonomski opravdani ukoliko i u maloj mjeri poboljšavaju stanje obraštenog trupa. Još jedan problem vezan uz obraštanje trupa broda je rasprostiranje biljnih i životinjskih vrsta svjetskim morima i posljedični poremećaj i promjena postojećih bio sustava. Energetska učinkovitost broda, koja se može razmatrati kroz potrošnju goriva po satu, brzinu i nosivost broda, ovisi o otporu broda, stupnju djelovanja propulzije, interakciji trupa broda i brodskog vijka, specifičnoj potrošnji goriva motora itd. Iako je u procesu osnivanja broda jedan od glavnih zadataka postizanje određene brzine tijekom pokusne plovidbe na određenom gazu uz određenu potrošnju goriva, jednako je tako važno točno predvidjeti brzinu broda i potrošnju goriva u uvjetima službe, što u velikoj mjeri ovisi o obraštanju. Važnost određivanja utjecaja obraštanja na hidrodinamičke značajke broda u mirnoj vodi jednako je važna kao i povećanje energetske učinkovitosti postojećih brodova kroz primjenu novih antivegetativnih premaza, budući da može ukazati na važnost čišćenja trupa broda i brodskog vijka. Naime, iako postoje razne mjere za povećanje energetske učinkovitosti brodova, brodovlasnici ili operateri broda oklijevaju s njihovom primjenom zbog nedostataka pouzdanih informacija vezanih uz potencijalne uštede prilikom primjene pojedine mjere. Točan proračun povrata investicije za brodovlasnike ili operatere broda je nužan prilikom odlučivanja o implementaciji pojedine mjere. Optimizacija čišćenja trupa broda i brodskog vijka jedna je od mjera za poboljšavanje energetske učinkovitosti brodova nad kojom brodovlasnik ili operater broda ima veliki utjecaj. Točna procjena utjecaja obraštanja na hidrodinamičke značajke broda, a time i utjecaja čišćenja na iste omogućuje procjenu ekonomskog aspekta ove mjere, što omogućuje točniju procjenu povrata investicije. Trenutno, ne postoji sveobuhvatna procedura pomoću koje bi se mogao pouzdano predvidjeti utjecaj obraštanja na hidrodinamičke značajke broda u službi. Iz tog razloga je International Towing Tank Conference (ITTC) izdao preporuke istraživačima da predlože nove izraze ili metode temeljene na eksperimentalnim podacima za predviđanje utjecaja obraštanja. Metodologija Glavni cilj ovog doktorskog rada je razvoj numeričkog postupka temeljenog na računalnoj dinamici fluida (RDF), koji može simulirati utjecaj obraštanja na strujanje oko oplakane površine obraslog trupa broda i brodskog vijka. Primjenom ovoga modela moguće je odrediti utjecaj obraštanja na hidrodinamičke značajke otpora i propulzije u mirnoj vodi. Model se temelji na hipotezi o sličnosti zida (eng. wall similarity hypothesis), odnosno na hipotezi da je utjecaj hrapavosti ograničen na unutarnji sloj turbulentnog graničnog sloja. Tako je ponašanje pojedine hrapave površine moguće opisati funkcijom hrapavosti, koja ustvari predstavlja pomak profila srednje brzine u logaritamskom području prema dolje. Kako ne postoji funkcija hrapavosti koja bi mogla opisati svaki pojedini tip hrapavosti, potrebno je izraditi studiju karakterizacije otpora (eng. drag characterization study) određenih tipova hrapavosti i na taj način odrediti tip funkcije i skalu hrapavosti, koja se koristi u definiciji Reynoldsovog broja hrapavosti. Za potrebe razvoja modela temeljenog na RDF-u, koji može simulirati utjecaj obraštanja, potrebno je odrediti model funkcije hrapavosti, odnosno ovisnost funkcije hrapavosti o Reynoldsovom broju hrapavosti. U ovome radu, temeljem studija karakterizacije otpora dostupnih u literaturi za dva tipa obraštanja, tj. biofilm i tvrdi obraštaj, predloženi su modeli funkcije hrapavosti za ova dva tipa obraštanja. Ovi modeli implementirani su u zidnu funkciju RDF rješavača te je time omogućeno određivanje utjecaja ovih tipova obraštanja na hidrodinamičke značajke proizvoljnog tijela. Budući da obraštanje ovisi o mnogim parametrima te da još uvijek nije moguće predvidjeti koliko dugo će antivegetativni premazi sprječavati obraštanje, ovo istraživanje je provedeno za unaprijed određena stanja obraštanja biofilmom i tvrdim obraštanjem. Razvijeni RDF model temelji se na Reynoldsovim osrednjenim Navier-Stokesovim jednadžbama (eng. Reynolds Averaged Navier-Stokes, RANS), koje su diskretizirane primjenom metode konačnih volumena. Provedena je detaljna analiza utjecaja odabira modela turbulencije na dobivene hidrodinamičke značajke broda te je pokazano kako k − Shear Stress Transport model turbulencije predstavlja najbolji kompromis između točnosti i proračunskog vremena za zatvaranje sustava RANS jednadžbi i osrednjene jednadžbe kontinuiteta. Za praćenje i određivanje položaja slobodne površine primijenjena je metoda udjela fluida u volumenu. Numeričke simulacije pokusa slobodne vožnje vijka provedene su primjenom metode višestrukih referentnih koordinatnih sustava, dok je utjecaj brodskog vijka u numeričkim simulacijama pokusa vlastitog pogona modeliran pomoću metode virtualnog diska. Ukupni otpor broda određen je numeričkim simulacijama temeljenim na viskoznom strujanju oko trupa broda, koje uključuju utjecaj slobodne površine (eng. Free Surface Simulation, FSS). Viskozni otpor određen je primjenom numeričkih simulacija s udvojenim modelom (eng. Double Body Simulation, DBS), koje ne uzimaju u obzir utjecaj slobodne površine. U DBS simulacijama simulira se strujanje oko duboko uronjenog udvojenog modela broda ili broda te je tako dobiveni ukupni otpor jednak viskoznom otporu. U DBS simulacijama otpor trenja određuje se integracijom tangencijalnih naprezanja po oplakanoj površini, dok se viskozni otpor tlaka određuje integracijom tlaka po oplakanoj površini. Otpor valova određuje se oduzimanjem ukupnog otpora dobivenog u FSS i viskoznog otpora u DBS simulacijama. Tako je omogućeno istraživanje utjecaja različitih stanja površine na pojedine komponente ukupnog otpora. Važan parametar za ukupni stupanj djelovanja propulzije predstavlja sustrujanje, koje se može podijeliti na efektivno i nominalno sustrujanje. Utjecaj stanja površine trupa broda na nominalno sustrujanje do danas je gotovo u potpunosti neistražen, zbog nemogućnosti mjerenja nominalnog sustrujanja na brodu u naravi. U okviru predloženog istraživanja nominalno sustrujanje broda u naravi određeno je prema proceduri opisanoj u ITTC preporukama, a utjecaj obraštanja na nominalno sustrujanje primjenom predloženih modela funkcija hrapavosti i određene skale hrapavosti te usporedbom s nominalnim sustrujanjem glatkog trupa broda. Utjecaj obraštanja na hidrodinamičke značajke slobodne vožnje vijka određen je primjenom predloženih modela funkcija hrapavosti i određene skale hrapavosti u numeričkim simulacijama pokusa slobodne vožnje vijka te usporedbom s hidrodinamičkim značajkama slobodne vožnje glatkog vijka. Utjecaj obraštanja trupa broda i brodskog vijka na propulzijske značajke istražen je numeričkim simulacijama pokusa vlastitog pogona obraslog i glatkog broda te analizom rezultata, koristeći rezultate pokusa otpora i rezultate pokusa slobodne vožnje brodskog vijka. Konačno, primjenom numeričkog postupka moguće je odrediti utjecaj obraštanja na radnu točku brodskog vijka, definiranu brzinom, snagom predanom vijku te brzinom vrtnje vijka. Premda procjena utjecaja obraštanja na hidrodinamičke značajke otpora i propulzije primjenom razvijenog numeričkog postupka omogućuje detaljan uvid u strujanje oko obraslog broda, primjena ovog postupka može biti složena za neiskusne korisnike RDF-a. Također, kako se RDF analiza sastoji od pripreme numeričkih simulacija, proračuna i analize dobivenih rezultata vrijeme trajanja proračuna je relativno dugo. U okviru optimizacije rasporeda čišćenja trupa broda i brodskog vijka bilo bi od velike koristi predložiti metodu, koja bi omogućila brzu procjenu utjecaja obraštanja na hidrodinamičke značajke otpora i propulzije. Štoviše, ITTC je izdao preporuke istraživačima da modificiraju trenutne poluempirijske izraze za koeficijent otpora i uzgona profila, a koji bi uzeo u obzir utjecaj obraštanja. U ovom radu, predložena je nova metoda za procjenu utjecaja obraštanja na hidrodinamičke značajke otpora i propulzije broda. Svrha ove metode je brza procjena utjecaja obraštanja na hidrodinamičke značajke broda. Štoviše, ova metoda omogućuje preciznije predviđanje negativnih učinaka obraštanja od onih dobivenih primjenom Granvilleove metode zakona sličnosti, kojom je moguće odrediti samo porast ukupnog otpora obraslog broda. Konačno, ovom metodom je moguće odrediti zaseban utjecaj obraštanja trupa broda i brodskoga vijka na ukupne hidrodinamičke značajke broda. Verifikacija i validacija Za potvrdu točnosti i vjerodostojnosti numeričkog postupka, potrebno je provesti opsežan postupak verifikacije i validacije dobivenih rezultata. Postupak verifikacije i procjene numeričke nesigurnosti u određivanju hidrodinamičkih značajki proveden je metodom indeksa konvergencije mreže (eng. Grid Convergence Index, GCI). Validacija dobivenih numeričkih rezultata za glatki brod i brodski vijak u naravi provedena je usporedbom s ekstrapoliranim rezultatima ispitivanja u bazenu, kao i s rezultatima dostupnih numeričkih istraživanja. Postupak validacije numeričkog postupka proveden je usporedbom dobivenih numeričkih rezultata s eksperimentalnim rezultatima studija karakterizacije otpora dostupnim u literaturi. Nadalje, dobiveni porasti koeficijenta otpora trenja primjenom predloženog numeričkog postupka uspoređeni su s porastima koeficijenta otpora trenja dobivenih primjenom Granvilleove metode zakona sličnosti. Promjene značajki otpora i propulzije broda dobivene primjenom novo predložene metode za procjenu utjecaja obraštanja uspoređene su s dobivenim rezultatima predloženog numeričkog postupka. Znanstveni doprinos Ovaj rad je rezultirao predloženim numeričkim postupkom modeliranja obraštanja, kao i novo predloženom metodom predviđanja značajki, koji se mogu primijeniti za procjenu utjecaja obraštanja na hidrodinamičke značajke otpora i propulzije broda u mirnoj vodi. To uključuje određivanje utjecaja obraštanja trupa broda i brodskog vijka na ukupni otpor, otpor trenja, otpor viskoznog tlaka i otpor valova, kao i na strujanje oko trupa broda i brodskog vijka. Nadalje, omogućeno je određivanje utjecaja obraštanja na značajke slobodne vožnje brodskog vijka, koeficijent upijanja, nominalno i efektivno sustrujanje, koeficijent napredovanja za točku vlastitoga pogona, stupanj djelovanja u slobodnoj vožnji, koeficijent prijelaza te kvazi propulzivni koeficijent. Uz to omogućeno je određivanje radne točke brodskog vijka, kao i energetske učinkovitosti broda za različita stanja obraštanja trupa broda i brodskog vijka. Procjena ovih utjecaja može se koristiti u alatima za optimizaciju rasporeda čišćenja trupa broda i brodskog vijka, gdje je predviđanje optimalnog vremena čišćenja od presudne važnosti. Dobiveni rezultati ukazuju na važnost čišćenja trupa broda i brodskog vijka. Glavni znanstveni doprinosi sastoje se od predloženih modela funkcija hrapavosti za određivanje otpora trenja površina prekrivenih biofilmom te njihova implementacija u programski paket za RDF, određivanje utjecaja biofilma na značajke otpora i nominalnog sustrujanja broda, značajke slobodne vožnje brodskog vijka i propulzijske značajke broda. Daljnji znanstveni doprinos je implementacija modela funkcije hrapavosti za određivanje otpora trenja površina prekrivenih tvrdim obraštanjem, određivanje utjecaja tvrdog obraštanja na značajke otpora i nominalnog sustrujanja broda, značajke slobodne vožnje brodskog vijka i propulzijske značajke broda. Konačno, znanstveni doprinos predstavlja novo predložena metoda za brzu procjenu utjecaja obraštanja na značajke otpora i propulzije broda.
Databáze: OpenAIRE
Externí odkaz: