CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (http://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   EHDFoam (http://www.cfd-online.com/Forums/openfoam-programming-development/105182-ehdfoam.html)

Nicolao July 24, 2012 13:57

EHDFoam
 
Hey everyone,

I'm trying to create a solver for electrohydrodynamics in OpenFoam. I've found one solver at this link

http://openfoamwiki.net/index.php/Co...ectromagnetics

but I'm having some problems about it.

The charge density in the case (EHDdroplet) is not calculated. And I don't understand why the electric force has been set up at the surface (I guess that it is at the surface) at the right side of the N-S equation.

I'm wondering if there is anyone working on a solver for EHD cases that could help me with it.

Just in case, I'm working with atomization process.

Thanks for any help!!

salehda November 16, 2012 16:33

Hi Nicola ,

any updates about these ?? I am trying to use it too , but I am stuck with compiling it

Thanks
David

Ivooo November 20, 2012 05:12

Working on it & almost there...
 
Hi,

I am currently working on an extended version of the EHDFoam model, using the version posted on the wiki as a base, but also using the JCP article of Lopez Herrera et al. [1] as a guide. This paper outlines how a charge density is generated on the interface between two fluids due to an electric field. As such, I am now able to simulate both dielectric and conducting fluids (any combination of the two), and I was able to verify most cases outlined in the paper. I'll put it on the wiki as soon as we have everything worked out in detail, but I am willing to collaborate with the last steps of development. Let me know if you're interested.

[1] Journal of Computational Physics 230:1939-1955, 2011

Ivooo November 20, 2012 05:15

Quote:

Originally Posted by Nicolao (Post 373295)
Hey everyone,

I'm trying to create a solver for electrohydrodynamics in OpenFoam. I've found one solver at this link

http://openfoamwiki.net/index.php/Co...ectromagnetics

but I'm having some problems about it.

The charge density in the case (EHDdroplet) is not calculated. And I don't understand why the electric force has been set up at the surface (I guess that it is at the surface) at the right side of the N-S equation.

I'm wondering if there is anyone working on a solver for EHD cases that could help me with it.

Just in case, I'm working with atomization process.

Thanks for any help!!

Just to answer your question, the charge density in the original EHDFoam model is not calculated, it is only transported so you'll have to initialise the charge densities e.g. using setFields. The electric force is solved explicitly (RHS), but I dont see a problem in that; it is the normal way of incorporating an additional force in the NS eqs.

Ivooo November 20, 2012 05:20

Quote:

Originally Posted by salehda (Post 392617)
Hi Nicola ,

any updates about these ?? I am trying to use it too , but I am stuck with compiling it

Thanks
David

I'm not sure what goes wrong in your case, but the EHDFoam has been written using OF1.5. I can't tell what happens if you compile it with a more recent version of OF, I started off by updating EHDFoam (both the library and the executable) to match the 2.1 version of OF.

Nicolao December 4, 2012 12:55

Quote:

Originally Posted by Ivooo (Post 393144)
Hi,

I am currently working on an extended version of the EHDFoam model, using the version posted on the wiki as a base, but also using the JCP article of Lopez Herrera et al. [1] as a guide. This paper outlines how a charge density is generated on the interface between two fluids due to an electric field. As such, I am now able to simulate both dielectric and conducting fluids (any combination of the two), and I was able to verify most cases outlined in the paper. I'll put it on the wiki as soon as we have everything worked out in detail, but I am willing to collaborate with the last steps of development. Let me know if you're interested.

[1] Journal of Computational Physics 230:1939-1955, 2011


Hey Ivooo, I also used the paper from Herrera as a base. And my solver is actually working quite fine. The fact is that I didn't implemented the equations from 13 to 18 (from the paper) on the solver. I want to know if it is necessary to put them on the solver, or if they are inherent to the solution of the Maxwell's equations.
(I don't know if that was clear enough...)

If you could help me on implementing them (equations 13 to 18) it would be of great help.

About the charge density, I initialized it on the 0 folder, and set as 0 as initial condition. And it worked. It appears only at the interface and the electric force too.

Thanks!

Rob Bacon January 23, 2013 13:24

Any success?
 
Hi, I'm looking to use a working EHD solver for analysis of Taylor cone ejection and I was wondering if either of your solvers is at a working stage yet? If so I'd be very interested in 'beta testing' it for you!

Ivooo February 5, 2013 09:37

Quote:

Originally Posted by Rob Bacon (Post 403599)
Hi, I'm looking to use a working EHD solver for analysis of Taylor cone ejection and I was wondering if either of your solvers is at a working stage yet? If so I'd be very interested in 'beta testing' it for you!

Sorry for the late response. I'll publish the code as soon as I have permission from my supervisors/sponsors (I will get that permission, shouldnt worry about it). Our intention is to publish it including several validation and example cases so people can start using it right away.

Rob Bacon February 16, 2013 20:00

Hi Ivooo, that would be great, looking forward to trying it!

Mojtaba.a July 20, 2013 07:46

Quote:

Originally Posted by Ivooo (Post 406103)
Sorry for the late response. I'll publish the code as soon as I have permission from my supervisors/sponsors (I will get that permission, shouldnt worry about it). Our intention is to publish it including several validation and example cases so people can start using it right away.

Dear Ivo,
have you got permission from you supervisor to publish the code?
I really look forward to it.

Thank you

Ivooo August 2, 2013 07:49

Quote:

Originally Posted by Mojtaba.a (Post 440947)
Dear Ivo,
have you got permission from you supervisor to publish the code?
I really look forward to it.

Thank you

Hi,

We are preparing an article about it now, which also involves a multi-region solver for EHD. If we're confident that the paper gets accepted I can publish it.


Cheers

chirs February 3, 2014 16:03

Hi David,

Currently, I am trying to set up EHD solver in openFOAM 2.2.2. but I can't get compile it, please could you give some advise; how to do it? how did you do it ?

Thanks in advance,
Christian

salehda May 29, 2014 12:44

Hi Chris,

I did not compile it , I left openFOAM and then I am getting back to it. I tried compiling a new library o include the two phase mixture, every thing goes well until it gives me an error

In file included from EHDTwoPhaseMixture/EHDtwoPhaseMixture.C:27:0:
EHDTwoPhaseMixture/EHDtwoPhaseMixture.H:190:39: error: conflicting return type specified for ‘virtual const volScalarField& Foam::EHDtwoPhaseMixture::nu() const’
../EHDincompressible/transportModel/transportModel.H:84:37: error: overriding ‘virtual Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::transportModel::nu() const’
make: *** [Make/linux64GccDPOpt/EHDtwoPhaseMixture.o] Error 1

Any idea, how to solve this ??

chirs June 4, 2014 05:58

Hi David,
Sorry for my late answer, I fix that problem as follows:
changing the 190 line for:
//- Return the kinematic laminar viscosity
virtual tmp<volScalarField& nu() const
{
return nu_;
}

salehda June 6, 2014 11:58

Thanks Chris,

I fixed this but now I am getting this ,
gammaEqnSubCycle.H: In function ‘int main(int, char**)’:
gammaEqnSubCycle.H:3:15: error: ‘piso’ was not declared in this scope
In file included from gammaEqnSubCycle.H:22:0,
from EHDFoam.C:75:
gammaEqn.H:6:26: error: ‘class Foam::interfaceProperties’ has no member named ‘cGamma’
In file included from gammaEqnSubCycle.H:30:0,
from EHDFoam.C:75:
gammaEqn.H:6:26: error: ‘class Foam::interfaceProperties’ has no member named ‘cGamma’
/opt/openfoam222/src/finiteVolume/lnInclude/readPISOControls.H:3:15: warning: unused variable ‘nOuterCorr’ [-Wunused-variable]
/opt/openfoam222/src/finiteVolume/lnInclude/readPISOControls.H:15:16: warning: unused variable ‘transonic’ [-Wunused-variable]
/opt/openfoam222/src/finiteVolume/lnInclude/readPISOControls.H:3:15: warning: unused variable ‘nOuterCorr’ [-Wunused-variable]
/opt/openfoam222/src/finiteVolume/lnInclude/readPISOControls.H:6:15: warning: unused variable ‘nCorr’ [-Wunused-variable]
/opt/openfoam222/src/finiteVolume/lnInclude/readPISOControls.H:12:16: warning: unused variable ‘momentumPredictor’ [-Wunused-variable]
/opt/openfoam222/src/finiteVolume/lnInclude/readPISOControls.H:15:16: warning: unused variable ‘transonic’ [-Wunused-variable]
/opt/openfoam222/src/finiteVolume/lnInclude/readTimeControls.H:38:8: warning: unused variable ‘maxDeltaT’ [-Wunused-variable]
make: *** [Make/linux64GccDPOpt/EHDFoam.o] Error 1

It seems that all the variable names has changed from the different versions , Can you please help me with that or can you upload the updated code if you were able to compile it.

Thank you

salehda August 9, 2014 13:29

Quote:

Originally Posted by Ivooo (Post 443480)
Hi,

We are preparing an article about it now, which also involves a multi-region solver for EHD. If we're confident that the paper gets accepted I can publish it.


Cheers

Hi Ivoo,

Any updates on that ?

Ivooo August 11, 2014 06:16

Quote:

Originally Posted by salehda (Post 505182)
Hi Ivoo,

Any updates on that ?

Yes, the paper has been submitted and is currently under review. It has been about a month now, so I guess I will hear something soon.

chirs April 20, 2015 13:16

hi Ivo,
I was reading your conference paper "An OpenFOAM-based electro-hydrodynamical model", it is interesting, I am working in something similar, but I could not set up the validation case for charge relaxation, please could you give me an idea how to set up it?

Best regards,

Chris

Ivooo April 21, 2015 09:28

Quote:

Originally Posted by chirs (Post 542805)
hi Ivo,
I was reading your conference paper "An OpenFOAM-based electro-hydrodynamical model", it is interesting, I am working in something similar, but I could not set up the validation case for charge relaxation, please could you give me an idea how to set up it?

Best regards,

Chris

If I recall correctly, I have used the funkySetFields in order to init a Gaussian charge bump. Perhaps that's where you are stuck?

Anyway, perhaps it is coincidence but today I received notification that our article on interFoamEHD including multi-region is available online. That means I can also upload the code and validation cases I have used.

I am currently digging through the files (I want the package with solvers, utilities, scripts and cases to be clean and tested once more). I will get back to this soon.

chirs April 21, 2015 10:01

Quote:

Originally Posted by Ivooo (Post 543006)
If I recall correctly, I have used the funkySetFields in order to init a Gaussian charge bump. Perhaps that's where you are stuck?

Anyway, perhaps it is coincidence but today I received notification that our article on interFoamEHD including multi-region is available online. That means I can also upload the code and validation cases I have used.

I am currently digging through the files (I want the package with solvers, utilities, scripts and cases to be clean and tested once more). I will get back to this soon.

Ivo, thanks to answer. Yes, the Gaussian charge bump is the problem. I will try your suggestion.
Congratulations for your paper, it will be great if you can upload the code and the validations.

Best Regards,

Christian Narvaez


All times are GMT -4. The time now is 18:18.