Difference between revisions of "Contrib/EHDFoam"

From OpenFOAMWiki
m
m (Wyldckat moved page Contrib electromagnetics to Contrib/EHDFoam: Renamed since the original title was misleading)
 
(5 intermediate revisions by 2 users not shown)
Line 5: Line 5:
  
 
== Model Equations ==
 
== Model Equations ==
:<math>\nabla\cdot\vec{V}=0</math>
+
:<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>\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>\vec{\vec{\tau}}=\mu\left(\nabla\vec{V} + \nabla\vec{V}^T\right)</math>
:<math>\vec{F}_b=\rho_E\nabla\vec{E} - \frac{1}{2}\left|\vec{E}\right|^2\nabla\varepsilon</math>
 
 
:<math>\frac{\partial \alpha}{\partial t} + \nabla\cdot\left(\alpha\vec{V}\right) = 0</math>
 
:<math>\frac{\partial \alpha}{\partial t} + \nabla\cdot\left(\alpha\vec{V}\right) = 0</math>
:<math>\varepsilon\nabla\vec{E}=\rho_E</math>
 
:<math>\vec{E}=-\nabla U_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>\rho=\rho_1\alpha + \rho_2\left(1-\alpha\right) </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>\mu=\mu_1\alpha + \mu_2\left(1-\alpha\right) </math>
Line 18: Line 18:
 
:<math>\sigma=\sigma_1\alpha + \sigma_2\left(1-\alpha\right) </math>
 
:<math>\sigma=\sigma_1\alpha + \sigma_2\left(1-\alpha\right) </math>
  
<!--   [[Image:Equation.jpg]]     the formula shown in the jpg were horrible and inconsistent (different notation was used). i left the link here though -->
+
<!-- [[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 ==
Line 32: Line 32:
  
 
[[Media:EHDdroplet.zip]]
 
[[Media:EHDdroplet.zip]]
 +
 +
[[Category:Electromagnetics solvers]]

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