buoyantBoussinesqPimpleEHDFoam

An OpenFOAM solver to model enhancement of mixed convection heat transfer due to electrohydrodynamic (EHD) flow, developed by Navid Zehtabiyan-Rezaie.

Description

The developed solver handles governing equations in mixed convection heat transfer including continuity, momentum, energy, and electrical equations. The developed solver is able to handle the coupled equations. For a comprehensive understanding of this solver, see this publication.

Target platform

The code has been rigorously tested and verified to be fully compatible with OpenFOAM v-4.1, ensuring its smooth integration and reliable performance with this specific release.

How to compile and use the solver

1- Download the source code using the following command:

git clone https://github.com/nzhtbyn/buoyantBoussinesqPimpleEHDFoam.git

2- To compile the solver, go to your work directory via the following command:

cd $WM_PROJECT_USER_DIR

Copy the folder to your work directory, and compile the new solver with the following command

wmake

3- Copy the folder validationCase to your run directory. To initiate the simulation, execute the following command:

.//Allrun

Charge density boundary condition

One of the problems in solving the conservation of space charge density is the application of a boundary condition on the emitting electrode. Since all essential ionization processes next to the emitting electrode are neglected, strict boundary conditions for the charge density cannot be formulated. The problem can be solved by using the Kaptzov hypothesis. He suggested that if the corona discharge occurs at some point of the emitting electrode and charge is injected, the electric field at this point remains at the value it takes at the corona onset. According to the above-mentioned explanations, a common method is used in numerical studies to properly estimate the charge density boundary condition on the emitting electrode. First, the charge density at the wire is guessed and iterated until the electric field is sufficiently close to Peek’s law, which specifies the onset electric field on the emitting electrode. Read more, here.

Validation

The local heat transfer coefficient over a plate with a heat flux of 187 $\text{W}/\text{m}^2$ is compared with the experimental results of Owsenek and Seyed-Yagoobi. Three different voltages are applied to the wire placed at a distance of 2 cm above the heated plate. Results obtained from the developed solver show a very good agreement with the experimental data.

The performance of the solver has been further assessed through extensive testing against reference data in the study of Zehtabiyan-Rezaie et al. (2020).

How to cite

Please, cite this library as:

@article{ZEHTABIYANREZAIE_ehdHeat_JoE_2020,
author = {Navid Zehtabiyan-Rezaie and Majid Saffar-Avval and Kazimierz Adamiak},
title = {Forced convection heat transfer enhancement using a coaxial wire-tube corona system},
journal = {Journal of Electrostatics},
volume = {103},
pages = {103415},
year = {2020},
issn = {0304-3886},
doi = {https://doi.org/10.1016/j.elstat.2019.103415}
}