Difference between revisions of "Extend-bazaar/template"

From OpenFOAMWiki
(Wray Agarwal 1-Equation Turbulence Model and Example Case)
 
(9 intermediate revisions by 2 users not shown)
Line 1: Line 1:
This template is only a suggestion aimed at saving your time. Please feel free to use any other format on the subpage for your contribution!
 
  
'''Click "Edit" to access and copy the page code.'''
 
  
 
== Description ==
 
== Description ==
This code implements  the Wray-Agarwal one-eqn turbulence closure model. [Developed on OpenFOAM V4.1]
+
This code solves the continuity and momentum equations along with the transport equations for cations and anions and the Poisson equation. The electrostatic equations are coupled to the fluid flow equations through the convective terms in ionic transport equations. The ionic transport equations are coupled to each other through the electrostatic Poisson equation. This solver is written to address the internal electrohydrodynamic flow. However, since this code is developed from the InterFOAM solver, it can easily be used for the external electrohydrodynamic flow and can also address the multiphase flow.
  
    Reference:
+
The fixed gradient boundary condition library in OpenFOAM is modified to create a customized gradient boundary condition for the ions that can be updated in each iteration. This boundary conditions accounts for the flux of counter-ions at the electrode surface due to the faradaic reactions.
    \verbatim
+
        Wray, T. J., Agarwal, R. K.,
+
        A New Low Reynolds Number One Equation Turbulence Model Based on a k-omega Closure,
+
        AIAA Journal, Vol. 53, No. 8, 2015, pp. 2216-2227
+
    \endverbatim
+
  
 
== Installation ==
 
== Installation ==
Download from https://openfoamwiki.net/images/4/4c/Contribution.zip
+
Download from https://www.dropbox.com/sh/p45ax8p4qvgbsyu/AABJt38xa5ghHRvOLorv5WU_a?dl=0
  
Then compile by following the instructions in the ReadMe file in the zip.
+
 
 +
 
 +
Then compile by running:
 +
<bash>cd EHDIonFOAM/
 +
./wmake
 +
</bash>
  
 
== Usage examples ==
 
== Usage examples ==
To use this code, select the model as the RAS model to be used in a simulation in the constant/turbulenceProperties dictionary.
+
To use this code, you have to create a channel that passes a fluid flow inside. A sample is given in page 31 of the following paper:
 +
[https://aip.scitation.org/doi/pdf/10.1063/1.5050793?class=pdf Kashir et al (2019)]
 +
 
 +
After creating the channel, you should specify the position of the electrodes in the channel. The discharge of counter-ions is already taken into account through the Boundary.H header file. The value of j here accounts for the flux of counter-ions and is found through dividing the measured electric current density at the electrode surface divided 
 +
 
  
The working variable was labeled RWA (to avoid possible conflict with the Reynolds stress tensor R).
 
  
An example test case, applying the model to the usual pitzDaly case is included in the zip.
+
Please have this mind that the channel flow is just a simple example. In reality electrostatic atomizers might have a hyperboloid structure. The only difference in complicated structures rather the sample channel flow presented in the paper is the complexity in grid creation. The usage of the EHDIonFOAM solver for understanding the electrification and charging mechanism inside electrostatic atomizers do not change.
  
 
== Documentation ==
 
== Documentation ==
This implementation is based on the Wray Agarwal model as presented in the NASA turbulence model research page (https://turbmodels.larc.nasa.gov/wray_agarwal.html).
+
This implementation is based on publication by [https://aip.scitation.org/doi/pdf/10.1063/1.5050793?class=pdf Kashir et al (2019)].
  
 
== History ==
 
== History ==
22 Mar 2018: First upload
+
23 Sep 2019: First upload

Latest revision as of 00:08, 23 September 2019


1 Description

This code solves the continuity and momentum equations along with the transport equations for cations and anions and the Poisson equation. The electrostatic equations are coupled to the fluid flow equations through the convective terms in ionic transport equations. The ionic transport equations are coupled to each other through the electrostatic Poisson equation. This solver is written to address the internal electrohydrodynamic flow. However, since this code is developed from the InterFOAM solver, it can easily be used for the external electrohydrodynamic flow and can also address the multiphase flow.

The fixed gradient boundary condition library in OpenFOAM is modified to create a customized gradient boundary condition for the ions that can be updated in each iteration. This boundary conditions accounts for the flux of counter-ions at the electrode surface due to the faradaic reactions.

2 Installation

Download from https://www.dropbox.com/sh/p45ax8p4qvgbsyu/AABJt38xa5ghHRvOLorv5WU_a?dl=0


Then compile by running:

cd EHDIonFOAM/
./wmake

3 Usage examples

To use this code, you have to create a channel that passes a fluid flow inside. A sample is given in page 31 of the following paper: Kashir et al (2019)

After creating the channel, you should specify the position of the electrodes in the channel. The discharge of counter-ions is already taken into account through the Boundary.H header file. The value of j here accounts for the flux of counter-ions and is found through dividing the measured electric current density at the electrode surface divided


Please have this mind that the channel flow is just a simple example. In reality electrostatic atomizers might have a hyperboloid structure. The only difference in complicated structures rather the sample channel flow presented in the paper is the complexity in grid creation. The usage of the EHDIonFOAM solver for understanding the electrification and charging mechanism inside electrostatic atomizers do not change.

4 Documentation

This implementation is based on publication by Kashir et al (2019).

5 History

23 Sep 2019: First upload