Difference between revisions of "Contrib/EHDFoam"

From OpenFOAMWiki
m (model Eqution)
m (Wyldckat moved page Contrib electromagnetics to Contrib/EHDFoam: Renamed since the original title was misleading)
 
(12 intermediate revisions by 3 users not shown)
Line 1: Line 1:
 
=EHDFoam=
 
=EHDFoam=
== introduction ==
+
== Introduction ==
  
Electrohydrodynamics deal with fluid motion induced by electric fields. In the mid 1960s G.I. Taylor introduced the leaky dieletric model to explain the behaviour of droplets deformed by steady field {cite}, and J.R. Melcher used it extensively to develop electrohydrodynamic{cite}. Here, we developed a Numerical model of electrohydrodynamics, the model concides the Taylor-Melcher Leaky Dielectric Model with VOF (Volume of Fraction) method. The Numerical solver was created in OpenFOAM-1.5.
+
Electrohydrodynamics deal with fluid motion induced by electric fields. In the mid 1960s G.I. Taylor introduced the leaky dieletric model to explain the behaviour of droplets deformed by steady field, and J.R. Melcher used it extensively to develop electrohydrodynamic. Here, we developed a Numerical model of electrohydrodynamics, the model concides the Taylor-Melcher Leaky Dielectric Model with VOF (Volume of Fraction) method. The Numerical solver was created in OpenFOAM-1.5.
  
 
== Model Equations ==
 
== Model Equations ==
 +
:<math>\vec{E}=-\nabla U_E</math>
 +
:<math>\varepsilon\nabla\cdot\vec{E}= - \varepsilon\Delta U_E=\rho_E</math>
 +
:<math>\frac{\partial \rho_E}{\partial t} + \nabla\cdot\left(\rho_E\vec{V}\right) + \frac{\sigma}{\epsilon}\rho_E =0</math>
 +
:<math>\vec{F}_b=\rho_E\vec{E} - \frac{1}{2}\left|\vec{E}\right|^2\nabla\varepsilon</math>
 +
:<math>\frac{\partial\vec{V}}{\partial t} + \nabla\cdot\left(\vec{V}\vec{V}\right)= -\frac{1}{\rho}\nabla p + \frac{1}{\rho}\nabla\cdot\vec{\vec{\tau}} + \frac{1}{\rho}\vec{F}_b + \vec{g}</math>
 +
:<math>\nabla\cdot\vec{V}=0</math>
 +
:<math>\vec{\vec{\tau}}=\mu\left(\nabla\vec{V} + \nabla\vec{V}^T\right)</math>
 +
:<math>\frac{\partial \alpha}{\partial t} + \nabla\cdot\left(\alpha\vec{V}\right) = 0</math>
 +
:<math>\rho=\rho_1\alpha + \rho_2\left(1-\alpha\right) </math>
 +
:<math>\mu=\mu_1\alpha + \mu_2\left(1-\alpha\right) </math>
 +
:<math>\varepsilon=\varepsilon_1\alpha + \varepsilon_2\left(1-\alpha\right) </math>
 +
:<math>\sigma=\sigma_1\alpha + \sigma_2\left(1-\alpha\right) </math>
  
[[Image:Equation.jpg]]
+
<!--  [[Image:Equation.jpg]]       the formula shown in the jpg were horrible and inconsistent (different notation was used). i left the link here though -->
  
== source code ==
+
== Source Code ==
  
sources can be downloaded below:
+
Sources can be downloaded below:
[[Media:EHDFoam.rar]]
+
[[Media:EHDFoam.zip]]
  
 +
==How to Install==
  
==how to install==
+
Compile the EHDTwoPhaseMixture first and then compile EHDFoam.
  
compile the  EHDTwoPhaseMixture first and then EHDFoam
+
== Test Case ==
  
 +
[[Media:EHDdroplet.zip]]
  
 
+
[[Category:Electromagnetics solvers]]
== test case ==
+
 
+
[[Media:EHDdroplet.zip]]
+

Latest revision as of 19:53, 21 March 2015

1 EHDFoam

1.1 Introduction

Electrohydrodynamics deal with fluid motion induced by electric fields. In the mid 1960s G.I. Taylor introduced the leaky dieletric model to explain the behaviour of droplets deformed by steady field, and J.R. Melcher used it extensively to develop electrohydrodynamic. Here, we developed a Numerical model of electrohydrodynamics, the model concides the Taylor-Melcher Leaky Dielectric Model with VOF (Volume of Fraction) method. The Numerical solver was created in OpenFOAM-1.5.

1.2 Model Equations

\vec{E}=-\nabla U_E
\varepsilon\nabla\cdot\vec{E}= - \varepsilon\Delta U_E=\rho_E
\frac{\partial \rho_E}{\partial t} + \nabla\cdot\left(\rho_E\vec{V}\right) + \frac{\sigma}{\epsilon}\rho_E =0
\vec{F}_b=\rho_E\vec{E} - \frac{1}{2}\left|\vec{E}\right|^2\nabla\varepsilon
\frac{\partial\vec{V}}{\partial t} + \nabla\cdot\left(\vec{V}\vec{V}\right)= -\frac{1}{\rho}\nabla p + \frac{1}{\rho}\nabla\cdot\vec{\vec{\tau}} + \frac{1}{\rho}\vec{F}_b + \vec{g}
\nabla\cdot\vec{V}=0
\vec{\vec{\tau}}=\mu\left(\nabla\vec{V} + \nabla\vec{V}^T\right)
\frac{\partial \alpha}{\partial t} + \nabla\cdot\left(\alpha\vec{V}\right) = 0
\rho=\rho_1\alpha + \rho_2\left(1-\alpha\right)
\mu=\mu_1\alpha + \mu_2\left(1-\alpha\right)
\varepsilon=\varepsilon_1\alpha + \varepsilon_2\left(1-\alpha\right)
\sigma=\sigma_1\alpha + \sigma_2\left(1-\alpha\right)


1.3 Source Code

Sources can be downloaded below: Media:EHDFoam.zip

1.4 How to Install

Compile the EHDTwoPhaseMixture first and then compile EHDFoam.

1.5 Test Case

Media:EHDdroplet.zip