Non-transferred arc hydrogen plasma torch simulations in OpenFOAM

Abstract

Dette prosjektet utførte en litteratur studie og numeriske simuleringsforsøk på termisk hydrogen plasma. Hovedmålet var å modellere hydrogen plasma i en ikke-overført lysbue plasma fakkel. Hovedmotivasjonen bak er å gå fra karbondioksid basert metalloksidreduksjon til renere prosesser som bruker hydrogenplasma. Plasma reaktorer er dyre å bygge, og oppførselen til plasma er dårlig forstått siden det er vanskelig å hente ut eksperimentell data fra plasma fakler. Simulering i OpenFOAM er derfor et perfekt verktøy for å genere data som kan hjelpe til med forståelsen av hvordan plasma oppføre seg, og hvordan en designer bedre reaktorer. En modell og løser prosedyre i OpenFOAM ble bygget ved hjelp av dataene og modellen fra litteraturstudiet. Modellen var et stasjoner, stabilt laminært system med en aksesymmetrisk geometri. Denne aksesymmetrien var sentrert på symmetriaksen til katoden. Styreligninger for masse, moment og energi ble implementert sammen med en modell for elektromagnetiske felter. Disse modellen muliggjorde koblingen mellom den elektriske strømmen og bevaringsligningen til plasmaet via Joule oppvarming og Lorentzkraften. I tillegg ble thermofysiske polymer implementert for gjenskaper den dynamiske responsen plasmaet har til varme fra den elektriske strømmen. Simuleringen av hydrogen og argon plasma ble gjennomført på modellen og en enkel parameter studie ble gjennomført. Argon resultatene ble sammenlignet med tidligere arbeid på samme modell som i dette prosjektet. Parametertestene på hydrogen og argon ble sammenlignet med tidligere arbeid i feltet. Resultatene viser at temperaturfeltet ligner på tidligere arbeid, men temperaturen var lavere enn forventet. Denne lave temperaturen kan skyldes at polynomene implementert ikke var nøyaktige nok. Parametertestingen viser at noe av responsen til plasmaet er som forventet, men lysbuen oppførte seg ikke slik litteraturen tilsier. Dette resultatet viser at den elektromagnetiske modellen og implementasjonen av Lorentzkraft trenger mer arbeid.

This project conducted a literature study and numerical simulations on thermal hydrogen plasma. The main goal was to model hydrogen plasma in a non-transferred plasma torch. The primary motivation behind the simulation is the move from carbon dioxide-based reduction processes of metal oxides to cleaner reduction processes with hydrogen plasma. Plasma reactors are expensive to build, and the inner workings of the torch are not entirely understood as it is challenging to extract experimental data. Hence, modeling it in software such as OpenFOAM can create data that might help understand hydrogen plasma behavior and design better plasma reactors. A case and solver procedure were constructed in OpenFOAM using data and models from the literature study. The case was a steady-state laminar system with an axisymmetric geometry. The axisymmetry was centered on the symmetry axis of the cathode in the non-transferred plasma torch. Governing equations for mass, momentum, and energy were implemented, along with an electromagnetic model. These models enabled the coupling of an electric current with the conservation equations of the plasma through the source terms of Joule heating and Lorentz force. In addition, polynomials of thermophysical properties were implemented to capture the dynamic response of the plasma to the heat from the electric current.

Hydrogen and argon plasma was simulated using the case setup, and a simple parameter study was performed. The argon results were compared to earlier work, which uses the same geometry and parameters. In addition, the results of the hydrogen and argon parameter study were also compared to the literature of parameter studies on hydrogen and argon. The results show that the temperature field resembles the earlier work, but the temperature was lower than expected. The low temperature could be attributed to the polynomial regression, which is not accurate enough to capture the plasma properties. Furthermore, the response to the parameters such as inflow velocity and current resembles litterateur results, but the arc behaves drastically different from what is expected. This result shows that the electromagnetic model and Lorentz force implementation needs more attention. In addition, the hydrogen results are questionable as they did not meet the requirement for convergence.