NOTE: End date changed to 01/19/2025 per NSSC information (Ed., 3/10/24).

PeleLM is a state-of-the-art reacting flow code tailored with a unique integration scheme that efficiently couples high-speed processes with the slower evolution of the flow structures using the adaptive mesh refinement (AMR) strategy. This project is to improve the ion chemistry (particularly employing reduced chemical mechanisms available in the published literature) to generalize the electric field configurations in PeleLM, and to then simulate zero-g coflow flames in the E-FIELD Flames dataset. Comparisons include V-I (voltage-ion current) curves and flame shapes (deduced from CH* images), varying with fuel types, compositions and flow rates. The data includes step function changes in voltage that can be used to evaluate the time response of the flame. It is important to also consider the impact of image and ion current collection timescales, which will be critical for model validation.

The work takes advantage of the unique E-FIELD Flames data, measured in microgravity, to validate the proposed PeleLM model augmentations, to permit a detailed investigation of the ion dynamics under the influence of electric fields, and to evaluate the subsequent interactions with non-charged species, as well as to assess the potential of electric field forcing to affect thermal transport and sooting in diffusion flames. Once the simulation is validated against the zero-g data, it will be possible to further extend the simulations to flames in earth gravity, where potential applications for improved heat release and emission reduction can be explored. This utilization of zero gravity data from the E-FIELD Flames experiment for beneficial purposes on earth and for a better understanding combustion control is compatible with the NASA objectives.

The E-FIELD Flames experiment was conceived and directed to create a data set that would permit the detailed assessment of the effects of ion-driven winds and charged species chemistry on non-premixed flames. A microgravity environment is required because the forces of ion driven winds and of gravity driven buoyancy winds are too close in magnitude to effectively separate their effects with any earth-gravity experiments where natural convection is a significant contributor to the combustion process. This rationale has been proven out by the E-FIELD flame experimental data now available as part of the PSI website. The Electric Field Effects on Laminar Diffusion Flames (E-FIELD Flames) experiment was completed in 2018, and with nearly 150 hours of operations over six months and the ignition of 250 electric field flames on the ISS, it was a very successful measurement campaign. The experiment aspect of this project is to extract and assemble a data portfolio that is ideally suited for validating numerical simulations. The simulation aspect of the project employs the PeleLM framework at the National Renewable Energy Laboratory (NREL), which was previously developed in the Lawrence Berkeley National Laboratory (LBNL). This open-source code is for evolving chemically reacting flows at low Mach number with block-structured adaptive mesh refinement (AMR) using the AMReX library. It is the extension of this code that we exploit for comparison with E-FIELD Flames data, with the ultimate goal of understanding and then controlling hydrocarbon flames under the influence of electric fields. PeleLMeX is the non-subcycling version of PeleLM, and this project will use the non-subcycling strategy at the same time understanding other aspects using the PeleLM.

Control of combustion using electric fields has for decades been proposed for extending flammability limits, reducing emissions, and preventing instability and blowoff, as well as modifying soot production. Implementation of the concept has been difficult, however, because of the complex coupling between buoyant driven flows, ion production in the flame, ion acceleration by the electric field, and the resulting forces on the neutral gas. With an efficient and accurate simulation capability, it will be possible to explore temporal and geometric variations of the electric field to manipulate flames in a wide range of beneficial ways.

As the GRI-Mech 3.0-based super reduced mechanism was implemented but without excited species and with a few ions, in the prior fiscal year, the main focus of this cycle was to implement a more complicated San Diego-based reduced mechanism with positive ions, negative ions, and excited species. The specific aims and objective of this project in this year include:

(1) Establishing flame images for fuel dilution cases, particularly for the 70% and 40% methane dilution in 1-g laboratory electric field flames (2) Understanding transient flame response, particularly the influence of the ion-driven wind in the 1-g in-laboratory electric field flames with the experimental flow conditions and mesh distance matching the microgravity E-FIELD Flames ISS experiments (3) Verify the strategy with understanding the detailed electric field with respect to the PeleLMeX function, and then implementing PeleLMeX CFD with a (newer) reduced chemical mechanism that includes more varieties of charged and excited species (4) Verify the fidelity of the PeleLMeX simulation of flames under the influence of electric fields by direct comparison using the previously validated reduced mechanism without electric fields in comparison to schlieren images of the flame under 1g conditions (5) Demonstrate validated PeleLM simulations of microgravity flames in steady conditions and during unsteady transitions in response to changes in the external electric field

Computational Modeling Verification (with and without electric fields) — This task continues from the previous computational modeling effort in 0g ISS condition of E-FIELD Flames, using a reduced mechanism to predict ion concentrations in flames and the effects of ion driven winds in coflow. The focus is to simulate in gravity condition and observe how the computational fluid dynamics (CFD) simulation captures the difference between microgravity, using the E-FIELD Flames PSI dataset for effective validation of the computational models and 1g experiments. [Ed. Note: PSI is the NASA Physical Sciences Informatics (PSI) data repository for Physical Sciences experiments performed on the ISS.]

Experiment of Fuel Dilution and Optical Imaging Diagnostics in 1g using ISS Conditions — This test is to establish the methane fuel dilution part in laboratory experiments similar to the conditions on the ISS, understand the flame dynamics using Schlieren imaging, and further validate the results of the computational models described above.

Computational Modeling Implementation Employed with a Set of Species — To predict ion concentrations in flames and the effects of ion driven winds in coflow and jet flames, a complete set of charged species and excited species to include in the computation is essential and the most challenging. This task is to find the best scientific method to include more charged species and excited species implementation into the CFD simulation with appropriate reduced chemistry to ensure robust simulation capability.

The previous PeleLM simulation computation capability validation using pressurized conditions also continues. Further details of the project progress through the third year are described in the following.

PeleLMeX Computational Results and Discussion

The progress in the PeleLM computational aspects of the project are described in a series of conference papers provided in the publication section of this report and in recent American Society for Gravitational and Space Research (ASGSR) presentations (Grabon et al., 2024, Deu Morel et al., 2023). Much of the work has so far concentrated on the cases without an electric field in order to assess the complications of the code implementation, its conversion to cylindrically symmetric geometry, to identify the choice of reduced models, and to ensure that the boundary conditions are appropriate to simulate an extruded fuel tube.

Using PeleLM with Various Pressures and Fuel Dilution with Different Species – no electric field

This part of the progress continues to investigate high pressure methane diffusion flames with water addition in comparison to CO2 addition using PeleLM simulation (Esquivias et al., 2024). The study used the same geometry coflow burner as for E-FIELD Flames and examined its behavior under the pressure conditions of previous tests of water addition simulation. By examining a similar flame geometry, but under different conditions, the study continues to evaluate the potential for PeleLM simulation to capture the chemico-physical details of small diffusion flames. This work is to study the behavior of a methane diffusion flame with various amounts of water vapor and CO2 added to the fuel and compares the flame behavior with various percentages of water vapor and CO2 content (0% to 65%). It provides the detailed temperature, CO, and CO2 and other species concentration distribution using the PeleLM adaptive mesh refinement (AMR) code. This work computes the combustion behavior of a coflow burner flame jet using PeleLM. The flame was simulated to pressures of 1.0, 1.4, 5.7, and 11.1 atm. The results extracted and analyzed include temperature profiles and various species mole fractions compared with their equilibrium state. The fundamentals between water and CO2 addition on the flame behaviors, peak temperature, as well as maximum mole concentration of CO2 and CO, was studied. The details of how the peak values of species vary for water and CO2 addition as compared to their equilibrium values are investigated.

The peak concentration for species CO, CO2, OH, and H2O were plotted and compared between water addition and CO2 at two pressure conditions, 11.1 and 1 atm. In 1 atm, CO2 addition has less CO and H2O production in comparison to water addition, while the peak value of OH is generally higher in all of the CO2 addition cases. At 11.1 atm, OH is very similar with both diluents added, this is due to the compact nature of the flame under the high pressure that does not allow enough room for the thermo-chemical reaction to reveal itself. Because of the addition of CO2, there is a distinct reduction in the production of CO, meaning that water addition participated in the reaction as both chemical and physical effect, but again, these are the peak value comparison only at atmospheric and high pressure conditions. Detailed reaction evaluations still are needed, as well as a sensitivity analysis. In general, the effect of pressure can be seen in the increase of species of CO, OH, and H2O.

The species of O, H, and OH of water and CO2 addition were compared to equilibrium conditions. The ratio between the OH, O, and H mole fraction species and the mole fraction at equilibrium are plotted against with the increase of the addition concentration. The mole fraction at equilibrium is computed using the Colorado State University online calculator, which is based on the original Stanjan Equilibrium Code (Dandy, 2019). The comparison conditions are 0%, 30% and 65% water addition and CO2 addition to a stoichiometric fuel/air mixture. The ratio of, for example, Hthis sstudy/ Heq., presents the ratio between the equilibrium prediction and the kinetics result for these atoms at all pressure conditions. The results show that H and O are consistently more than two times their equilibrium expectation for all water content at lower pressure, while OH stays near the equilibrium value throughout all the conditions in this study.

There is also relatively limited effect of diluent composition, as the thermal images of flames diluted with water and carbon dioxide are very similar under similar pressure and dilution conditions. The addition of water and carbon dioxide has similar results on the flame behavior. Future work will include an analysis of existing high-pressure experiment results and the superequilibrium prediction for the species (H and O) at different pressures.

These details show that the PeleLM framework is able to capture the proper boundary conditions, even though the absolute geometry is slightly different, because the experiment has a slightly extruded fuel tube and the simulation has a flat inlet. The flow inlet is varied to approximate the conditions expected for the extruded case. Further results appear in all of the conference presentations and publications listed in the publications section of this report.

Preliminary Results 1g Electric Fields Using ISS Conditions – PeleLMeX

The PeleLMeX computation of this research is implemented on the UCI High Performance Community Computing Cluster (HPC3) running on Linux 8.6. Its aim is to boost the effectiveness of simulations by utilizing a considerable number of cores. To handle user admission and resource allocation, the cluster employs entities like login nodes and the Slurm scheduler.

The goal of this project is to use the reduced mechanism Deu Morel et al. (2023, see reference) tested for microgravity condition with a newer version of the PeleLMeX as well as newer HPC3 hardware. This work initiated with reproducing the behavior of the coflow burner used in the ISS. The primary distinction between the numerical and experimental setups lies in the extruded fuel pipe. As the PeleLMeX code has not integrated the extruded geometry, it is unable to replicate the impact of the fuel pipe on the inflow velocities or temperature. The effect of the extruded tube is simulated using a customized inlet velocity profile. The primary assumption is that the Low Mach number regime is in effect. It is acknowledged that the velocities will remain within this regime. This is a widely accepted hypothesis in the study of combustion and was validated at the conclusion of the simulations.

The simulation was initiated with 19 (and 15) cm/s fuel flow velocity with 3.1 cm/s coflow air with 300K initial temperature with CFL 0.1. This case without the electric field can help explain the operation of PeleLMeX. Moreover, the cases without electric field are distinctive for shorter compilation time, due to the reduced number of equations/species to be solved resulting from the deactivation of the electric field. The reduced mechanism (by Di Renzo et al., 2022), comprising 26 species and 134 reactions with only two positive ions, H3O+ and HCO+, and one negative ion, O2- and electrons, was tested in microgravity and gravity conditions. This mechanism is currently used for simulating under 1g condition in comparison to the laboratory results to validate or invalidate the transient flame dynamic response to the ion wind in describing the E-FIELD Flames experiment.

As described in prior reports, for these simulations, the electric field implementation is in operation, requiring the use of HPC3, a supercomputer, due to its substantially longer computational time. Multiple convergence difficulties emerged during these simulations. After resolving these matters, we determined the necessary computation time to achieve a stable state. With the use of the electric field method, it takes 1 hour and 10 minutes to reach 1000 steps, in contrast to 3 minutes without it. Additionally, a stable state is achieved in the simulation within 60 milliseconds in simulation time. With the electric field, 1,000 steps correspond to a maximum of 10E-4 seconds. This requires over 30 days of consecutive simulations on the supercomputer to attain a steady state.

1g In-Laboratory Experiments

The 1g in-laboratory experiments include a coflow burner, rotameters for fuel and air, and a high voltage power supply (Trek 609A) with shunt current resistance for measuring ion current. The details of the experiment setups are described in Chien et al. (2018) Energies 11(5), 1235 and Chien et al. (2019) Combust. Flame 204(250-259). Both the High Voltage Power Supply (HVPS) output and ion current are monitored and controlled all along the running tests using data acquisition system NI 6341. The acquisition signals are sampled at a rate of 1000 samples per second and recorded with a LABVIEW code. The electrode potential is programmed before every recording allowing precise voltage sweeps. For safety, a High Voltage and Reverse Protection Circuit is used in addition to the HVPS. This specialized electronic arrangement is designed to protect electronic devices from potential damage caused by high voltage and reverse current. The High Voltage Protection circuit is composed of Zener diodes, which will open the circuit if the voltage reaches its threshold. This particular device safeguards the sensitive part of the circuit from being harmed by preventing electric arc formation, for instance. The Reverse Current Protection circuit is used to protect electronic components from being damaged by reverse polarity. If the circuit is incorrectly connected with the power supply by the user, diodes will block the current and avoid component destruction.

In this research, the mesh is set at 27 mm height above the burner due to the height limitation of the current experimental setup. The data collected includes ion current to applied voltages, natural flame images using a Nikon D90 camera, and also Shlieren imaging.

The measurements are made with two fuel flow speeds, 19 cm/s and 15 cm/s with airflow of 3.1 cm/s at 100% methane fuel concentration. The flame images were taken at ISO 200, aperture f/4.5, and a shutter speed of 1/10s, and the pictures were taken in the dark with a black fabric piece placed on the back to avoid reflections from the glass casing.

The ion current measurement is performed using a shunt resistor circuit. It consists of a 1 𝑀Ω resistor connected in parallel with two Zener diodes. This circuit is symmetrical, which means the polarity isn’t a concern here. The ion current results for positive electric fields are measured and recorded as a confirmation. The results indicate that a positive E-FIELD leads to the appearance of three regions. In the sub-saturation region, ion current increases with applied voltage until it reaches a plateau, known as the saturation region, which is reached at around 1.65 kV/cm for each flame speed. Finally, there is a final increase in ion current in the supersaturation region. The limits between the saturation regions are not the same depending on the flame velocity due to the total amount of carbon content input. Therefore, the saturation plateau appears at a lower electric field strength for lower velocity flames.

To understand how the ion-driven wind and the process of electric fields affect the flame dynamics, the Shlieren imaging optics diagnostic technique is performed. The Schlieren imaging is an optical setup composed of two parabolic convergent mirrors (f=48’’), a punctual light source, a razor blade and the Nikon D90 camera as the focal plane. The Z-type optical setup is used due to space limitations. The principle of Schlieren imaging is that the light rays reflected by the first mirror can be considered parallel because the light source is placed at its focal distance. If nothing is placed between the mirror, the light continues its path through the system and the image projected on the camera sensor is blank. However, if a change in density, which implies a change of refraction index, occurs between the two mirrors, then the ray light will not arrive parallel in the second mirror. These rays will then either be blocked by the razor blade and appear darker or pass above the razor blade and appear brighter.

Schlieren imaging is efficient when it comes to observing density changes and will be employed to validate the simulation. The changes in the flame shape are more remarkable on the 19 𝑐𝑚/𝑠 flame in comparison to the lower flow rate. The Schlieren imaging also shows that the thermal plume is wider for the higher flow rate. As the electric field strength increases for the positive field, the flame is concentrated and pulled upward to the mesh and the thermal plume is narrower. For the negative field configuration, the flame is pushed toward the burner and becomes wider; the thermal plume would be widened instead. When the saturation plateau is reached at around 1. 6 𝑘𝑉/𝑐𝑚, the flame shape doesn’t change significantly, nor does the thermal plume as the electric field strength increases. The electric field can affect the soot production (orange region of the flame), which can be a great potential for application as soot production can indicate a more efficient combustion.

[Ed. Note: Per the Principal Investigator (PI), the investigation gathered images of negative field and positive field with fuel dilution at 70% and 40%. These images resemble the fuel dilution conditions from the 0g experiments, and further analysis is underway.]

Comparison Between Simulation and Experiments

An initial comparison of the Schlieren imaging and the density profile from the PeleLMeX computation, specifically for flame widths for both flame velocities, are conducted. The result is computed from the reduced mechanism (Di Renzo 2022) under gravitational conditions, which was tested in microgravity. We understand that the comparison is qualitative as each image is normalized relative to its maximum intensity. A detailed Schlieren imaging is needed as well as a higher acquisition rate, for capturing the images between electric field charges can help further understand flame temperature and fully authenticate the chemical mechanism including more positive, negative ions and excited species from the PeleLMeX. The work is currently being carried out on a local high performance computing cluster, but there are limited hours available for that system so additional computational resources will be explored. Further results appear in the publications identified in that section of this report.

Preliminary Outcome for 0g Implementation with More Species – PeleLMeX

To bridge the gap between simulation outcomes and empirical data, this part of the work is to evaluate potential mechanism candidates: (A) 30 species and 86 reaction using Coffee mechanism as foundation (1993) by Pederson and C. Brown; (B) 26 species and 134 reactions but only 4 ionic species using GRI Mech 3.0 by Renzo et al.; and (C) the recently published 45 species and 216 reactions using San Diego reduced mechanism, by López-Cámara et al. This task is to implement the recently published reduced San Diego mechanism in microgravity. These mechanisms are chosen for their comprehensive coverage of ionic and exited species within the flame, compared to the mechanism currently used by Di Renzo et al. which is not able to capture the ion current behavior in microgravity condition as part of work from the previous fiscal year.

The main difference lies in the (C) reduced San Diego mechanism chemical model, which includes 11 charged species, including H3O+, HCO+, C2H3O+, CH5O+, O2−, OH−, e−, CO3−, CHO2−, O−, and CHO3− (in comparison to the previous used (B) H3O+, HCO+, O2− and e−) and included excited species CH and OH. This mechanism includes many species, and as the electric field is introduced, the calculation became extensive and the computation time increased substantially. With various tests and approaches for launching the mechanism at various stages, the simulation remains computing with extended calculation time, unresolved and/or unstable with inconsistent results. The heat release from the calculation can be massive at a single cell and that can lead to a relatively large local flow, and/or the chemical solver may fail to sort through the reaction and lead to confusion. Further tests on the mechanism, particularly on sorting through how the solver processes each step in the calculation, as well as how the mechanism works with a relatively newer version of the simulation, are needed.

Conclusions

The modeling work is progressing with the newer version of PeleLMeX using the reduced mechanism (with four ions, electron, and no excited species) in gravitational environment. The experiments with diluents are completed, and the ion current data and the Schlieren imaging are verified using ISS conditions. Locally with PeleLM non-subcycling version, the code has also continued to be exercised under different coflow flame conditions with different levels of fuel dilution with carbon dioxide and different pressures to demonstrate the validity from the previous work with water addition of the boundary condition implementation. Challenges remain for how to implement a relatively large San Diego based mechanism including more ions and charged species into the newer version of PeleLMeX. Continuing work on this project focuses more heavily on the computations with detailed assessment of the ion current prediction and the role the chemical mechanism plays in that prediction.

References

Esquivias, B. and Chien, Y.-C. (2024), “A comparison between water addition and CO2 addition to a diffusion jet flame,” Combustion Science and Technology, Special issue of the 29th ICDERS. doi: 10.1080/00102202.2024.2380083

Deu Morel, R., Day, M. and Chien, Y.-C. (2023) “Investigation of mechanisms for microgravity ion-driven wind using simulation,” Graduate Poster Presentation, American Society for Gravitational and Space Research Annual Meeting 2023, Nov. 14 –18, Washington DC.

Grabon, K., Wimer, N., Deak, N., Day, M, and Chien, Y.-C. (2024) “Ion-driven wind of E-FIELD Flames,” Combustion 2 Presentation, American Society for Gravitational and Space Research Annual Meeting 2024, Dec. 3 – 7, San Juan, Puerto Rico.

Invited Technical Lectures/Presentations

NASA Physical Sciences Informatics (PSI) User Group, “ACME — E-FIELD Flames,” April 11th, 2024 (Virtual).

Women in Technology at UCI (WITIN) inclusive networking, “NASA Keck, Diversity & Lovin’ UCI,” Irvine, California. May 17th, 2024 (Virtual).

PeleLM is a state-of-the-art reacting flow code tailored with a unique integration scheme that efficiently couples high-speed processes with the slower evolution of the flow structures using the adaptive mesh refinement (AMR) strategy. This project is to improve the ion chemistry (particularly employing reduced chemical mechanisms available in the published literature) to generalize the electric field configurations in PeleLM, and to then simulate zero-g coflow flames in the E-FIELD Flames dataset. Comparisons include V-I (voltage-ion current) curves and flame shapes (deduced from CH* images), varying with fuel types, compositions and flow rates. The data includes step function changes in voltage that can be used to evaluate the time response of the flame. It is important to also consider the impact of image and ion current collection timescales, which will be critical for model validation.

The work takes advantage of the unique E-FIELD Flames data, measured in microgravity, to validate the proposed PeleLM model augmentations, to permit a detailed investigation of the ion dynamics under the influence of electric fields, and to evaluate the subsequent interactions with non-charged species, as well as to assess the potential of electric field forcing to affect thermal transport and sooting in diffusion flames. Once the simulation is validated against the zero-g data, it will be possible to further extend the simulations to flames in earth gravity, where potential applications for improved heat release and emission reduction can be explored. This utilization of zero gravity data from the E-FIELD Flames experiment for beneficial purposes on earth and for a better understanding combustion control is compatible with the NASA objectives.

The E-FIELD Flames experiment was conceived and directed to create a data set that would permit the detailed assessment of the effects of ion-driven winds and charged species chemistry on non-premixed flames. A microgravity environment is required because the forces of ion driven winds and of gravity driven buoyancy winds are too close in magnitude to effectively separate their effects with any earth-gravity experiments where natural convection is a significant contributor to the combustion process. This rationale has been proven out by the E-FIELD flame experimental data now available as part of the PSI website. The Electric Field Effects on Laminar Diffusion Flames (E-FIELD Flames) experiment was completed in 2018, and with nearly 150 hours of operations over six months and the ignition of 250 electric field flames on the ISS, it was a very successful measurement campaign. The experiment aspect of this project is to extract and assemble a data portfolio that is ideally suited for validating numerical simulations. The simulation aspect of the project employs the PeleLM framework at the National Renewable Energy Laboratory (NREL), which was previously developed in the Lawrence Berkeley National Laboratory (LBNL). This open-source code is for evolving chemically reacting flows at low Mach number with block-structured adaptive mesh refinement (AMR) using the AMReX library. It is the extension of this code that we exploit for comparison with E-FIELD Flames data, with the ultimate goal of understanding and then controlling hydrocarbon flames under the influence of electric fields. PeleLMeX is the non-subcycling version of PeleLM, and this project will use the non-subcycling strategy at the same time understanding other aspects using the PeleLM.

Control of combustion using electric fields has for decades been proposed for extending flammability limits, reducing emissions, and preventing instability and blowoff, as well as modifying soot production. Implementation of the concept has been difficult, however, because of the complex coupling between buoyant driven flows, ion production in the flame, ion acceleration by the electric field, and the resulting forces on the neutral gas. With an efficient and accurate simulation capability, it will be possible to explore temporal and geometric variations of the electric field to manipulate flames in a wide range of beneficial ways.

As the prior fiscal year has tested a reduced mechanism with positive ions, the specific aims and objective of this project in this year include:

(1) Establishing 1-g in-laboratory electric field flames data with the experimental flow conditions and mesh distance matching the microgravity E-FIELD Flames ISS experiments (2) Verify the reduced mechanism with two positive ions, one negative ion, and electrons using the PeleLMeX CFD (3) Verify the strategy with understanding the detailed electric field with respect to the PeleLMeX function, and then implementing PeleLMeX CFD with a (newer) reduced chemical mechanism that includes more varieties of charged and excited species (4) Verify the fidelity of the PeleLMeX simulation of flames under the influence of electric fields by direct comparison with E-FIELD Flames data for both 0g and 1g conditions (5) Demonstrate validated PeleLM simulations of microgravity flames in steady conditions and during unsteady transitions in response to changes in the external electric field

Computational Modeling Verification (with and without electric fields) — This task continues from the previous computational modeling effort in 0g ISS condition of E-FIELD Flames using a reduced mechanism to predict ion concentrations in flames and the effects of ion driven winds in coflow. The focus is to simulate in gravity condition, and observe how the CFD simulation captures the difference between microgravity, using the E-FIELD Flames Physical Sciences Informatics (PSI) dataset for effective validation of the computational models and 1g experiments.

Experiment in 1g using ISS condition — This test is to establish the in-laboratory experiments similar to the conditions on the ISS, and validating the results of the computational models described above.

Computational Modeling Implementation Employed with a Set of Species — To predict ion concentrations in flames and the effects of ion driven winds in coflow and jet flames, a complete set of charged species and excited species to include in the computation is essential and the most challenging. This task is to find the best scientific method to include more charged species and excited species implementation into the CFD simulation with appropriate reduced chemistry to ensure robust simulation capability.

The previous PeleLM simulation computation capability validation using pressurized conditions also continues. Further details of the project progress through the second year in both the experimental and computational aspects are described in the following.

PeleLMeX Computations – Computational Results and Discussion

The progress in the PeleLM computational aspects of the project are described in a series of conference papers provided in the publication section of this report and in a recent publication (Chien, et al., 2023). [Ed. Note: See References and Cumulative Bibliography.] Much of the work has so far concentrated on the cases without an electric field in order to assess the complications of the code implementation, its conversion to cylindrically symmetric geometry, to identify the choice of reduced models, and to ensure that the boundary conditions are appropriate to simulate an extruded fuel tube.

Using PeleLM with various pressures and fuel dilution with different species – no electric field

This part of the progress investigates high pressure methane diffusion flames with water addition, in comparison to CO2 addition using PeleLM simulation (see reference). The study used the same geometry coflow burner as for E-FIELD Flames but examined its behavior under the pressure conditions of previous tests of water addition simulation. This work is to study the behavior of a methane diffusion flame with various amounts of water vapor and CO2 added to the fuel and to compare the flame behavior with various percentages of water vapor and CO2 content (0% to 65%). It provides the detailed temperature, CO, and CO2 and other species concentration distribution using the PeleLM adaptive mesh refinement (AMR) code. This work uses a coflow burner in PeleLM for simulating of the flame jet by computing the combustion behavior. The flame was simulated to pressures of 1.0, 1.4, 5.7, and 11.1 atm. The results extracted and analyzed include temperature profiles and various species mole fractions compared with their equilibrium state.

Regarding physical properties, both water and carbon dioxide present similar behaviors. At atmospheric pressure and high water content, the flame lifts off from the burner tip, while as the pressure rises the flame anchors back near the burner tip. In contrast, CO2 addition does not lift the flame from the burner tip. As pressure rises, the flame width reduces. Both water and CO2 addition decrease the flame peak temperature. This artificial horizontal expansion highlights the narrowing of the flame with changes in pressure while much smaller change is visible with changes in diluent fraction. There is also relatively limited effect of diluent composition as the thermal images of flames diluted with water and carbon dioxide are very similar under similar pressure and dilution conditions.

In order to provide more structure details, the mole fraction of oxygen, water, carbon dioxide, and carbon monoxide in the different conditions of the flames along the radius of the burner is observed. In both cases, the temperature profile follows the same path as water and carbon dioxide respectively. The maximum mole concentration of CO2 and CO were also investigated for all conditions. For the water addition, the CO2 and CO have an increase of maximum concentration when pressure rises. From 1 to 1.4 atm, there is a decrease in the maximum concentration of CO. As for the CO2 addition, there is an increase of CO in all conditions. H2O and CO2 concentration profile was also compared for all the conditions along the centerline of the diffusion flame. The rise of pressure produces an increase of water concentration approximately 0.4-0.6 mm from the burner, and it is observed that the water concentration decays slower than at lower pressure conditions. The addition of water and carbon dioxide has similar results on the flame behavior. Future work will include an analysis of existing high-pressure experiment results and the superequilibrium prediction for the species (H and O) at different pressures.

These details show that the PeleLM framework is able to capture the proper boundary conditions even though the absolute geometry is slightly different because the experiment has a slightly extruded fuel tube and the simulation has a flat inlet. The flow inlet is varied to approximate the conditions expected for the extruded case. Further results appear in all of the conference presentations and publications listed in the publications section of this report.

Preliminary Results 1g Electric Fields using ISS conditions - PeleLMeX

The PeleLMeX computation of this research is implemented on the UCI High Performance Community Computing Cluster (HPC3) running on Linux 8.6. Its aim is to boost the effectiveness of simulations by utilizing a considerable number of cores. To handle user admission and resource allocation, the cluster employs entities like login nodes and the Slurm scheduler.

The goal of this project is to use the reduced mechanism Donzeau et al. (2023, see reference) tested for microgravity condition with a newer version of the PeleLMeX as well as newer HPC3 hardware. This work initiated with reproducing the behavior of the coflow burner used in the ISS. The primary distinction between the numerical and experimental setups lies in the extruded fuel pipe. As the PeleLMeX code has not integrated the extruded geometry, it is unable to replicate the impact of the fuel pipe on the inflow velocities or temperature. The effect of the extruded tube is simulated using a customized inlet velocity profile. The primary assumption is that the Low Mach number regime is in effect. It is acknowledged that the velocities will remain within this regime. This is a widely accepted hypothesis in the study of combustion and was validated at the conclusion of the simulations. The chosen hypothesis pertains to the chemistry model. The mechanism was selected, comprising 26 species and 134 reactions. The project's goal is to validate or invalidate this chemical mechanism's efficacy in describing the E-FIELD experiment.

The simulation initiated with 15 and 19 cm/s fuel flow velocity with 3.1 cm/s coflow air with 300K initial temperature with CFL 0.1. This case without the electric field can help explain the operation of PeleLMeX. Moreover, the cases without electric field are distinctive for its shorter compilation time due to the reduced number of equations/species to be solved resulting from the deactivation of the electric field. The results depict the OH concentration in flames at different velocities. A comparison between these findings and the Abel inverse transform of the experimental flames is also conducted. A distinctive shape is observed for both flames with a reasonably uniform distribution of OH. The normalized colorbar facilitates comparison of concentrations. A greater concentration of OH is found in the 19 cm/s flame, which is expected due to the increased velocity leading to higher OH emission Into the flame. This effect is also apparent in the experiment in the following paragraphs. For an inlet velocity of 19 cm/s, velocity increases exponentially with respect to Y position. The velocity at y=0 is 19 cm/s and the maximum velocity of the mixture exceeds 86 cm/s, which is well above the flame. Velocity within the flame ranges between 19 and 35 cm/s.

For these simulations, the electric field implementation is in operation, requiring the use of HPC3, a supercomputer, due to its substantially longer computational time. Multiple convergence difficulties emerged during these simulations. The resolution is no longer straightforward, and there are residual values that struggle to converge during the solution of the Newton’s equations. To resolve these issues, the input includes a constant time step of 10E-6 seconds. Boundary condition issues on the cathode and anode are resolved by utilizing a Dirichlet condition in the X direction and a Neumann condition in the Y direction. To clarify, the Neumann condition is a flux condition, while the Dirichlet condition enforces a specific temperature at the boundary. After resolving these matters, we determined the necessary computation time to achieve a stable state. With the use of the electric field method, it takes 1 hour and 10 minutes to reach 1000 steps, in contrast to 3 minutes without it. Additionally, the preceding section illustrates that a stable state is achieved in the simulation within 60 milliseconds in simulation time. With the electric field, 1,000 steps correspond to a maximum of 10E-4 seconds. This requires over 30 days of consecutive simulations on the supercomputer to attain a steady state. Due to time constraints, the initiated simulations did not reach a steady state but are still being computed. Nonetheless, it is still possible to observe the variation in temperature profile between, with and without the electric field.

1g in-laboratory experiments

The 1g in-laboratory experiments include a coflow burner, rotameters for fuel and air, and a high voltage power supply (Trek 609A) with shunt current resistance for measuring ion current. The details of the experiment setups are described in Chien et al. (2018) Energies 11(5), 1235 and Chien et al. (2019) Combust. Flame 204(250-259). The mesh is set at 27 mm height above the burner. The data collected includes ion current to applied voltages, natural flame images using a Nikon D90 camera, and also OH* imaging from a PI-MAX 4 (Princeton Instruments).

The measurements are made with two fuel flow speeds, 15 cm/s and 19 cm/s, performed with d = 37 mm and airflow of 3.1 cm/s. Positive and negative electric fields are tested and analyzed. These results indicate that a positive E-FIELD leads to the appearance of three regions. In the sub-saturation region, ion current increases with applied voltage until it reaches a plateau, known as the saturation region, which is reached at around 1.65 kV/cm for each flame speed. Finally, there is a final increase in ion current at around 2.2 kV/cm in the supersaturation region. For the negative electric field, a plateau is reach around 0.6 kV/cm. However, it’s important to note that the presence and magnitude of the ion current in these flames may not provide a direct or straightforward indication of combustion efficiency.

Several factors can influence the ion current, making its interpretation more complex. Ionization processes are heavily influenced by temperature. The degree of ionization tends to rise with hotter temperatures. In a non-premixed flame, temperatures can fluctuate significantly within the combustion zone, resulting in regions with varying ionization levels. Regions with higher temperatures in the flame may display greater ionization and subsequently a higher ion current. This fluctuation in temperature poses a challenge to using ion current as a direct measure of combustion efficiency.

The flame visualization is captured by camera of each data point presented and related to the current–voltage curve (IV curve). There are four cases, two pertaining to speed and two to electric field.

When an electric field is positive, the flame tends to increase in width and appear bluer, indicating there is more air entrainment into the diffusion flame reacting with the fuel. Therefore, electric field can control the efficiency of the flame. Electric fields can indirectly accelerate the movement of ions and charged particles within the flame. This increases the ionization rate of the flame, promoting the formation of more ions. The increased ionization tends to enhance combustion by improving the mixing of fuel and oxidizer, which can result in a broader, more stable flame. The flame appears bluer, which is often associated with a higher combustion temperature and a more complete combustion process. Flames have an optimal or stoichiometric fuel-to-air ratio at which they burn most efficiently. This is the ideal balance of fuel and oxidizer that results in complete combustion.

Next, we observe how the introduction of electricity affects the formation of soot in the diffusion flame generated by the coflow burner. To increase the amount of soot particles in the flame, a high flow rate of methane up to 65 mL/min was utilized. Specifically, the fuel employed was CH4; using ethylene, for example, could result in even greater soot emission, which is why the flow rate was increased. Due to their yellow color, it is easier to track the greater volume of soot that results from the higher flow rate. A Nikon D90 camera was employed with the following settings: ISO 100, f/4.5, and a shutter speed of 1/10s. A positive electric field increases the velocity of the flame, causing it to go from 15,20 cm/s to over 100 cm/s. As a result, the soot residence time decreases and the amount of soot in the flame decreases as well. When a negative electric field is applied, the prevailing airflow towards the burner surface lengthens the time that soot particles spend in the flow, thus promoting soot accumulation. Results of this case were not taken with a higher flow rate since the previous negative field images have shown an increase of the soot volume in the flame.

The OH chemiluminescence images are captured and Abel inverted to reveal the local non-integrated chemiluminescence within the reaction sheet. While the quantitative value of the concentration of OH* cannot be determined, a qualitative comparison of the location of this excited species is possible after normalizing all the flame pictures. This information will be used for initial comparison with the computation results.

To analyze the variation in width with changes of electric field strength, a simple imaging method is used. This analysis is comparing the flame width variation with no a flame with no electric field applied. Wb is the width of the flame and Wb,0 is the width without electric field. The results are analyzed for the positive field for the two fuel flow speed. The result shows 3 sections for the width of the flame, similar to the ion current. It follows the IV curve, showing the first sub-saturation, the saturation zone, and then supersaturation. For a flame of 19 cm/s the width increases compared to Wb,0 but for the 15 cm/s flame the flame width decreases compared to Wb,0 but increases with the electric field. The use of an electric field is connected to stoichiometry. In slower flames, it decreases the initial width and then broadens it later. Conversely, higher velocity flames tend to expand beyond their initial width. A further analysis for the increase in flame width is in the sub-saturation zone for the two flowrates are also plotted with their linear regression. This may provide a means of comparing the development of flame width in relation to the electric field between the experiment and the simulation, particularly for the sub-saturation region.

Comparison between simulation and experiments

An initial comparison of the OH chemiluminescence and the OH species from the computation specifically for flame widths for both flame velocities are conducted. While we understand that the amount chemiluminescence is small and yet to be simulated, we would like to observe based on what we can get from the experiments and results from PeleLMeX. An example of images can present the comparison of (a) a flame from the experiment and (b) a simulated OH, both with a velocity of 19 cm/s and an electric field of 0 kV /cm.

The discrepancy in flame appearance is due to the scale, which may be misleading. At a similar scale, the experimental flame measures 3.15 mm, while the simulation displays a flame of 2.8 mm, resulting in a 11% error to the experiments. Although this may appear significant, uncertainties in measuring the experimental flame width must be considered. Even a slight variation in air flow during image capture can affect the width of the flame. Furthermore, the camera calibration has a relative uncertainty of 0.5 mm. A summary Table 1 of flame widths between the simulation and the experiment is presented.

OH* and OH are situated in the reaction sheets closed to each other (although the analysis is qualitative), while the flame width are closely within 10% difference. To fully authenticate the chemical mechanism, the next step is to compare the experiment and the simulation within the context of an electric field flame.

The work is currently being carried out on a local high performance computing cluster, but there are limited hours available for that system so additional computational resources will be explored. Further results appear in the conference presentations and publications identified in that section of this report.

Preliminary Outcome for 0g implementation with more species - PeleLMeX

To bridge the gap between simulation outcomes and empirical data, this part of the work is to evaluate potential mechanism candidates: (A) 30 species and 86 reaction using Coffee mechanism as foundation (1993) by Pederson and C. Brown, (B) 26 species and 134 reactions but only 4 ionic species using GRI Mech 3.0 by Renzo et al., and (C) the recent published 45 species and 216 reactions using San Diego reduced mechanism, by López-Cámara et al.. This task is to implement the recently published reduced San Diego mechanism in microgravity. These mechanisms are chosen for their comprehensive coverage of ionic and exited species within the flame compared to the mechanism currently used by Renzo et al. which is not able to capture the ion current behavior in microgravity condition as part of work from the previous fiscal year.

The main difference lies in the (C) reduced San Diego mechanism chemical model which includes 11 charged species, including H3O+, HCO+, C2H3O+, CH5O+, O2-, OH-, e-, CO3-, CHO2-, O-, and CHO3- (in comparison to the previous used (B) H3O+, HCO+, O2- and e-) and included excited species CH and OH. This mechanism includes many species and as the electric field is introduced the calculation became extensive and the computation time increased substantially. With various tests and approaches for launching the mechanism at various stages, the simulation remains computing with extended calculation time, unresolved and/or unstable with inconsistent results. The heat release from the calculation can be massive at a single cell and that can lead to a relatively large local flow, and/or the chemical solver may fail to sort through the reaction and lead to confusion. Further tests on the mechanism particularly on sorting through how the solver processes each step in the calculation, as well as how the mechanism works with a relatively newer version of the simulation is needed.

Conclusions

The modeling work are progressing with the newer version of PeleLMeX using the same reduced mechanism in gravitational environment. The experiments are conducted, and data are measured using ISS conditions. Locally with PeleLM non-subcycling version, the code has also continued to be exercised under different coflow flame conditions with different levels of fuel dilution with carbon dioxide and different pressures to demonstrate the validity from the previous work with water addition of the boundary condition implementation. Challenges remain for implementing a relatively large newer mechanism including more ions and charged species into the newer version of PeleLMeX. Continuing work on this project focuses more heavily on the computations with detailed assessment of the ion current prediction and the role the chemical mechanism plays in that prediction.

References:

Esquivias, B., Girodon, H. and Chien, Y.-C. (2023).“A Comparison Between Water Addition and CO2 Addition to a Diffusion Jet Flame,” abstract 263, International Colloquium on the Dynamics of Explosions and Reactive Systems, Siheung, Korea, July 23-28.

Donzeau, M., Esclapez, L., Day, M. S. and Chien, Y.-C. (2023). “Recent Progress on Numerical Modeling for Microgravity Electric Field Flames Results,” 13th U.S. National Combustion Meeting, March 19-22, 2023, College Station, Texas.

Invited Technical Lectures/Presentations for this PSI Project:

National Academies of Sciences Speaker, “Sustainability & ADEI – Research Scientist in Higher Education Academic Setting”, and Early-Career Panelist for the Committee on Biological and Physical Sciences in Space (CBPSS), Space Science Week 2023, Washington D.C., March 29th, 2023 (this project was presented).

PeleLM is a state-of-the-art reacting flow code tailored with a unique integration scheme that efficiently couples high-speed processes with the slower evolution of the flow structures using the adaptive mesh refinement (AMR) strategy. This project is to improve the ion chemistry (particularly employing reduced chemical mechanisms available in the published literature) to generalize the electric field configurations in PeleLM, and to then simulate zero-g coflow flames in the E-FIELD Flames dataset. Comparisons include V-I (voltage-ion current) curves and flame shapes (deduced from CH* images), varying with fuel types, compositions and flow rates. The data includes step function changes in voltage that can be used to evaluate the time response of the flame. It is important to also consider the impact of image and ion current collection timescales, which will be critical for model validation.

The work takes advantage of the unique E-FIELD Flames data, measured in microgravity, to validate the proposed PeleLM model augmentations, to permit a detailed investigation of the ion dynamics under the influence of electric fields, and to evaluate the subsequent interactions with non-charged species, as well as to assess the potential of electric field forcing to affect thermal transport and sooting in diffusion flames. Once the simulation is validated against the zero-g data, it will be possible to further extend the simulations to flames in earth gravity, where potential applications for improved heat release and emission reduction can be explored. This utilization of zero gravity data from the E-FIELD Flames experiment for beneficial purposes on earth and for a better understanding combustion control is compatible with the NASA objectives.

The E-FIELD Flames experiment was conceived and directed to create a data set that would permit the detailed assessment of the effects of ion-driven winds and charged species chemistry on non-premixed flames. A microgravity environment is required because the forces of ion driven winds and of gravity driven buoyancy winds are too close in magnitude to effectively separate their effects with any earth-gravity experiments where natural convection is a significant contributor to the combustion process. This rationale has been proven out by the E-FIELD flame experimental data now available as part of the PSI website. The Electric Field Effects on Laminar Diffusion Flames (E-FIELD Flames) experiment was completed in 2018, and with nearly 150 hours of operations over six months and the ignition of 250 electric field flames on the ISS, it was a very successful measurement campaign. The experiment aspect of this project is to extract and assemble a data portfolio that is ideally suited for validating numerical simulations. The simulation aspect of the project employs the PeleLM framework developed in the Lawrence Berkeley National Laboratory (LBNL), and recently transitioned to the National Renewable Energy Laboratory (NREL). This open-source code is for evolving chemically reacting flows at low Mach number with block-structured adaptive mesh refinement (AMR) using the AMReX library. It is the extension of this code that we exploit for comparison with E-FIELD Flames data, with the ultimate goal of understanding and then controlling hydrocarbon flames under the influence of electric fields.

Control of combustion using electric fields has for decades been proposed for extending flammability limits, reducing emissions, and preventing instability and blowoff, as well as modifying soot production. Implementation of the concept has been difficult, however, because of the complex coupling between buoyant driven flows, ion production in the flame, ion acceleration by the electric field, and the resulting forces on the neutral gas. With an efficient and accurate simulation capability, it will be possible to explore temporal and geometric variations of the electric field to manipulate flames in a wide range of beneficial ways. Hence, the specific aims and objective of this project include:

(1) Access, extract, and analyze ion current, voltage, and flame images from the PSI data base to serve as the benchmark microgravity data set for nonpremixed flames under the influence of an external electric field. (2) Modify PeleLM CFD software to approximate the physical boundary of an extruded tube coflow burner geometry with cylindrical symmetry (3) Extend PeleLM CFD software to include electrodynamics in the two-dimensional circular symmetric system (4) Implement PeleLM CFD with a reduced chemical mechanism that includes charged and excited species (5) Verify the fidelity of the PeleLM simulation of flames under the influence of electric fields by direct comparison with E-FIELD Flames data (6) Demonstrate validated PeleLM simulations of microgravity flames in steady conditions and during unsteady transitions in response to changes in the external electric field

Experimental Data (selection and analysis)--There are two major components of our using the E-FIELD Flames PSI dataset for effective validation of the computational models. First is the evaluation and regularizing the main features of the flame that will be used for the comparison. Second, is the selection of which data will serve as the baseline canonical dataset for comparison of computational findings.

Computational Modeling (with and without electric fields)--To predict ion concentrations in flames and the effects of ion driven winds in coflow and jet flames, a focus has been on a full CFD simulation with appropriate reduced chemistry, including first modeling the same burner geometry (but without electric fields) to ensure robust simulation capability.

Further details of the project progress through the first year in both the experimental and computational aspects are described in the following.

Experimental Results and Discussion

The experimental results and discussion part of the work (Chien, et al., 2022) described the preparatory and initial measurements of diffusion flames under the influence of an electric field aboard the International Space Station (ISS), as part of the Advanced Combustion via Microgravity Experiments (ACME) project. [Ed. Note: See References below and the Cumulative Bibliography in this Task Book record.] Intended as the foundation publication for the experimental effort, the work comprehensively included the space experiment methods, the capabilities of the ACME insert, experiment procedures, data, and limitations and constraints. The measurements presented included images and ion currents of small diffusion flames of methane (in air), subjected to a ramping electric field in microgravity. While there have been prior microgravity studies of Electric Field Effects on Laminar Diffusion Flames (E-FIELD Flames) using a drop tower, the published measurements represented the first mapping of the effects of an ion-driven wind on flame behavior under an electric field in the absence of the confounding influences of buoyancy. The paper described the challenges of remote measurement and manipulation of flames on the ISS and presented preliminary results from the first set of coflow flames. The results showed clearly that the flame is most compact at saturation while the measured voltage-current characteristic (VCC) curve demonstrates parabolic behavior after saturation which differs from observations in 1 g on Earth (shown in Chein, et al., 2022, figure 6). The flame images in microgravity of methane coflow conditions (in Chien, et al., 2022, figure 9) corresponding to the same comprehensive set of results. Identifying these images is nontrivial because the camera time stamp is not directly coincident with the scalar data information. Hence, one of the key accomplishments has been the linking of these two disparate data outcomes into a single unified comparison set. The flame images with different fuel flow rates 19 or 15 cm/s, and concentrations at 100%, 70% and 40% also show how different flames can look while still producing the same ion current as shown in Chein, et al., 2022 figure 11. The above scientific results provided the information that permitted the selection of 70% methane 30% nitrogen with a positive electric field applied as the most distinct condition to serve as the first comparison data set for the modeling.

The PSI data set is complete and the tools for extracting images and scalar data from that data set are achieved. Evaluating all of the different conditions showed that the 70% methane fuel case, with a positive electric field has the most consistent and distinct features while also providing a complete range of field strengths to use as the first comparison data set for the modeling efforts. The experiments also highlight the very clear first peak (called the acme point in Chien, et al., 2022) where the ion current reaches a saturation apex before decreasing, which is distinctly different from what happens in 1-g Earth gravity where a saturation plateau is seen. It is this feature that provides the clearest challenge to the computational model – to both demonstrate the validity of the model and to explain the physical processes that lead to this distinct acme point behavior in zero-g.

Computational Results and Discussion

The progress in the PeleLM computational aspects of the project are described in a series of conference papers provided in the publication section of this report and a recent publication (Chien, et al., 2023). Much of the work has so far concentrated on the cases without an electric field in order to assess the complications of the code implementation, its conversion to cylindrically symmetric geometry, to identify the choice of reduced models, and to ensure that the boundary conditions are appropriate to simulate an extruded fuel tube.

(a) Using PeleLM with various pressures and fuel dilution – no electric field

This part of the progress investigates high pressure methane diffusion flames with water addition using PeleLM simulation. The study used the same geometry coflow burner as for E-FIELD Flames but examined its behavior under the pressure conditions of previous ground-based experiments. This work provides a comprehensive preliminary understanding of how combustion changes when dilution (water vapor in this case) is introduced into the fuel flow from 0% - 60% mole fraction. It provides the detailed temperature, CO, and CO2 concentration distribution using the PeleLM adaptive mesh refinement (AMR) code. The results showed that the peak temperature in the flame reduces with water addition and increases with pressure increase. The general conclusion is that the pressure influence dominates over the dilution effect. This computational result helps to understand the scope of change when dilution is introduced into the high pressure flow and helps the next step of experiment planning (for possible future high pressure zero-g combustion studies). An example of these results appears in Chien, et al., 2023, figure 4.

This work concludes that the pressure has the most dominance over the water dilution conditions, while the flame is observed to lift at the 60% limit close to extinguishment. The work also plots the peak concentration of the species CO2, CO, O2 and OH at atmospheric pressure, to observe the relationships not only on OH + CO -> CO2, but also the O2 peak value change with the increase of dilution with water. The peak value of CO2 does not change throughout all the conditions but OH and CO are both decreased. O2 peak value remains similar with water dilution. Therefore, further analysis on how water vapor addition affects or participates in the reaction particularly on hydroxyl, carbon monoxide, water, and carbon dioxide is needed for understanding the influence from water addition, while the role of elevated pressures remains significant with the density change.

These details show that the PeleLM framework is able to capture the proper boundary conditions even though the absolute geometry is slightly different because the experiment has a slightly extruded fuel tube and the simulation has a flat inlet. The flow inlet is varied to approximate the conditions expected for the extruded case.

(b) Preliminary Results including Electric Fields - PeleLMeX

The first step was to select the conditions most appropriate for comparing simulations with experiments. The selected case from the 70% methane coflowing with air conditions starting from the acme point, already identified above as the 70% positive field situation. The first results of a zero gravity flame simulation with temperature plot, and the flame is subjected to an electric field influence is computed. The simulation shows clearly that the electric field pulls the flame upward and narrows it as would be expected. The detailed comparison of the ion current with the calculated ion current has the same order of magnitude but is not yet quantitatively matching. There are many avenues yet to explore, including the reduced chemistry mechanism. The current work is using the only published reduced mechanism for ions but that reduced mechanism was created for turbulent flames and may not provide the best performance for laminar flames. Continued work is needed on identifying the best mechanism to use. Another challenge with the simulations is the computational time needed. The work is currently being carried out on a local high performance computing cluster, but there are limited hours available for that system so additional computational resources will be explored. Further results appear in the conference presentations and publications identified in that section of this report.

Computational Conclusions

The computational progress has been substantial, with the conversion of PeleLM to a version that functions in cylindrical symmetry and one that includes the appropriate equations for the electric field and the ion mobility in that field (PeleLMeX). The code has also been exercised under different coflow flame conditions with different levels of fuel dilution and different pressures to demonstrate the validity of the boundary condition implementation, the ignition and steady state achievement in the code, and the capability of simulating an extruded tube geometry with a flat boundary by adjusting the inlet flow profile.

Continuing work on this project focuses more heavily on the computations with detailed assessment of the ion current prediction and the role the chemical mechanism plays in that prediction.

Reference:

Chien, Y.-C., Stocker, D., Hegde, U. and Dunn-Rankin, D. (2022) “Electric-Field Effects on Methane Coflow Flames Aboard the International Space Station (ISS): ACME E-FIELD Flames,” Combustion and Flame, Vol. 246. DOI: 10.1016/j.combustflame.2022.112443

Chien, Y.-C., Girodon, H., Esquivias Rodriguez, B. (2023), “Modeling of elevated pressure diffusion flames with water addition,” Combustion Science and Technology, Special issue of ICDERS 2022. doi: 10.1080/00102202.2023.2182205

Invited Technical Lectures/Presentations for this PSI Project:

1. National Academies of Sciences Speaker, “Sustainability & ADEI – Research Scientist in Higher Education Academic Setting,” and Early-Career Panelist for the Committee on Biological and Physical Sciences in Space (CBPSS), Space Science Week 2023, Washington D.C., March 29th, 2023 (this project was presented).

2. Technical talk and dinner at SWE Holiday Soirée (Society of Women Engineer – Orange County), Microgravity E-FIELD Flames Aboard the International Space Station (ISS): Research, Education and Embracing Diversity,” Irvine, California. December 8th, 2022.

3. Seminar in Institute Pprime - CNRS/ISAE-ENSMA (Centre national de la recherche scientifique - Institut P', Université de Poitiers – ISAE-ENSMA), Microgravity E-FIELD Flames on the ISS and Gas Hydrates for Combustion Research, November 24, 2022.

4. Seminar in ICARE - CNRS Orléans Campus (Centre national de la recherche scientifique - Institut de Combustion, Aérothermique, Réactivité et Environnement), Microgravity E-FIELD Flames on the ISS and Gas Hydrates for Combustion Research, Orléans, France, November 18, 2022.

5. Seminar for SIRiPods (Samueli Interdisciplinary Research in Pods) Junior Engineering Students at University of California, Irvine, hosted by Christine King — My Research Experience for E-FIELD Flames: Experiment Operation on the ISS from Earth, August 10, 2021.

6. Guest Speaker for the EXpanding Communities and Encouraging Leadership (EXCEL) students, 1st and 2nd year STEM undergrads in Bio/Chem at University of California, Irvine, hosted by Dr. Harris, the Director of CAMP, February 16, 2021.

PeleLM is a state-of-the-art reacting flow code tailored with a unique integration scheme that efficiently couples high-speed processes with the slower evolution of the flow structures using the adaptive mesh refinement (AMR) strategy. We have already completed simulations with PeleLM of a coflow burner, without an electric field, in Earth gravity, in comparison to experiments of OH fluorescence (Vicariotto, University of California Irvine dissertation, 2019). This proposal is to improve the ion chemistry and generalize the electric field configurations in PeleLM (with PeleLM’s primary developer, Dr. Day), and to then simulate zero-g coflow flames in the E-FIELD Flames dataset. Comparisons include V-I (voltage-ion current) curves and flame shapes (deduced from CH* images), varying with fuel types, compositions, and flow rates. The data includes step function changes in voltage that can be used to evaluate the time response of the flame. The Principal Investigator, Dr. Chien, was responsible for acquiring the data, and is an expert in the impact of image and ion current collection timescales, which will be critical for model validation.

The proposed work will take advantage of the unique E-FIELD Flames data, measured in microgravity, to validate the proposed PeleLM model augmentations, to permit a detailed investigation of the ion dynamics under the influence of electric fields, and to evaluate the subsequent interactions with non-charged species, as well as to assess the potential of electric field forcing to affect thermal transport and sooting in diffusion flames. Once the simulation is validated against the zero-g data, it will be possible to further extend the simulations to flames in Earth gravity, where potential applications for improved heat release and emission reduction can be explored. This utilization of zero gravity data from the E-FIELD Flames experiment for beneficial purposes on Earth and for a better understanding combustion control is compatible with the NASA objectives.

Reference: M. Vicariotto, Water vapor addition in high concentrations to the fuel side of a two-dimensional methane/air diffusion flame, Ph.D. Dissertation, UC Irvine, 2019.