EHDFoam

Pavithra · June 28, 2019, 01:37

Hello Everyone,

I am new to OpenFoam.

Based on the solver posted in this thread, I have combined electrical equations with buoyantBoussinesqPimpleFoam.

My intention is to have a solver for single phase EHD flow. I am able to compile the solver successfully.

Could some one please suggest me a good benchmark case to validate the solver.

Thank You.

With Thanks.
Pavithra.

Nicolao · June 28, 2019, 10:39

Dear Pavithra,

I guess you probably already seen this paper (from Herrera),

http://gfs.sourceforge.net/papers/ehd.pdf

but in there they show several validation cases for each equation of the leaky dielectric model (for EHD flows). Of course they use a two-phase approach, but for the Poisson and the charge conservation equations you can use the same cases they used there.

"My intention is to have a solver for single phase EHD flow"

Now, specifically for your case I guess you already have something in mind. =)

Cheers

Nicolao · June 28, 2019, 10:47

Quote:

Originally Posted by Pavithra

Hello Everyone,

I am new to OpenFoam.

Based on the solver posted in this thread, I have combined electrical equations with buoyantBoussinesqPimpleFoam.

My intention is to have a solver for single phase EHD flow. I am able to compile the solver successfully.

Could some one please suggest me a good benchmark case to validate the solver.

Thank You.

With Thanks.
Pavithra.

In the work of Roghair you can also see some validation tests. I'm not quite sure, but I guess those are similar to the ones used by Herrera. He also included some MultiRegion validation that he used specifically for his work.

https://link.springer.com/article/10...404-015-1581-5

Pavithra · June 29, 2019, 01:12

Quote:

Originally Posted by Nicolao

In the work of Roghair you can also see some validation tests. I'm not quite sure, but I guess those are similar to the ones used by Herrera. He also included some MultiRegion validation that he used specifically for his work.

https://link.springer.com/article/10...404-015-1581-5

Dear Sir,

Thank you so much for your kind reply. I will read the research papers which you have shared.

https://www.sciencedirect.com/scienc...04388615000030

I have tried to validate the solver with respect to the above article.

I am trying to replicate the velocity profile presented in Fig 7 of the above paper.

I am able to see electroconvection of charge. I am able to obtain a velocity profile of the same shape. But, the magnitude is 10^-10 smaller.

One possible reason might be some mistake in non-dimensional parameters. I have cross-checked it several times.

I am stuck at this point.

Thank You.

With Thanks,
Pavithra.

CloudBird · September 22, 2019, 09:53

Quote:

Originally Posted by Ivooo

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.

Dear Ivo,

I am so excited to know that you once prepared to share your code and your cases. I am doing some work with two phase electrohydrodynamics with VOF model recently but I met some problems. I try to download the code but it seems that the URL is now unavailable. Would you mind to upload your code again in other ways?

I'll be very grateful for any reply from you.

Best regards,

Golden

Pavithra · September 23, 2019, 02:23

Quote:

Originally Posted by CloudBird

Dear Ivo,

I am so excited to know that you once prepared to share your code and your cases. I am doing some work with two phase electrohydrodynamics with VOF model recently but I met some problems. I try to download the code but it seems that the URL is now unavailable. Would you mind to upload your code again in other ways?

I'll be very grateful for any reply from you.

Best regards,

Golden

Hi Cloud Bird,

Please find the files of InterEHDFoam attached, herewith.

XinXin · October 22, 2019, 10:04

Hello，everyone. I am simulating a water droplet impact surface under electrical field based on the interFoamEHD solver which designed by ivook.My model is 2D-Axisymmetric. I find that the charge density distribution is weird, there is some charge appearing above the droplet along the axis. I am puzzled by this problem. The charge density distribution is related to the velocity Magnitude.

I think the charge conversation equation is calculated for all computation field,both water and air in the interFoam solver. However, for water, it will get polarized and ohmic migration of charged ions. For air,it can be regarded as a perfect dielectric, there is no charge convection. I think we should use different physical model for water and air. But can it be realized in interfoam. I am comfused by the interfoam solver. How does this sovler distinguish the phase1 and phase2? I am reading the twoPhaseMixture.H, incompressibleTwoPhaseMiture.H ,but still get stuck.Could you please give me some advices?
charge distribution 2.png

charge distribution.png

equation.jpg

mcfdma · January 8, 2020, 14:37

Hello,

Firstly, thank you for the solver. I am trying to update the solver on OF6. However, the droplet breakup is significantly different to the original solver.

I have uploaded my solver on Google Drive due to exceeding file size [1]. Could someone point out where exactly am I going wrong in terms of modifying my solver.

I have attached the screenshot for the last time step with OF231 and OF6.

OF231
last_OF231.png

OF6
last_OF6.jpg

[1] https://drive.google.com/open?id=1Jd...KiYCL0Wc3JHxyi

tom_flint2012 · February 23, 2020, 18:11

Hi All,

Thanks for sharing the code. I hope this is an appropriate thread to ask my questions as they are a bit more general than the electrohydrodynamics solver discussed here, and I think its more of a magnetohydrodynamics question really, but I was wondering if this thread was relevant to my problem.

Basically I currently have a VOF solver fro simulating multi-component metallic melting and vapourisation, with energy equations, source terms in the NS equations for solidification etc. I use this to simulate power beam welding applications and arc welding applications, but I would really like to incorporate the electric charge transport and magnetic effects to look at the effect of the Lorentz force on the molten pool behaviour.

I have a mhd + interfoam solver that I constructed based on this work:

https://www.researchgate.net/publica...etallic_liquid

As I understand it I could either use an mhd approach or an ehd approach, and I was wondering if people had any thoughts on which may be appropriate. Are the source terms in the NS equations in the interehdfoam solver equivalent to the JxB lorentz force term in the mhd (and my intermhdfoam solver)??

But I would like to be able to apply a alternating electric potential as a boundary condition. Would you guys have any pointers on the best way for me to progress? I could work out the time varying B field using Amperes law based on an alternating current density (I think), and put this as my B boundary condition in my mhd+interfoam solver. But I thought I'd try and get some advice from people more experienced than me before I go down this route to try and save myself some time in the long run.

Thanks in advance for any advice you guys can offer. I'm very grateful.

All the best,

Tom

Ivooo · November 18, 2020, 07:52

I shared my code a while back, but the link seems to stop working every now and then. So I did what I should have done a long time ago, create a github repository. You can find it here (but still based on OF2.1.1): https://github.com/iroghair/openFoamEHD

Although I am not focussing my research on EHD specifically at this time, I welcome pull/merge requests; a few people in this thread already took the source and updated it, I would be happy to update any updates or other changes to help other users.

Nicolao · November 18, 2020, 08:05

Quote:

Originally Posted by Ivooo

I shared my code a while back, but the link seems to stop working every now and then. So I did what I should have done a long time ago, create a github repository. You can find it here (but still based on OF2.1.1): https://github.com/iroghair/openFoamEHD

Although I am not focussing my research on EHD specifically at this time, I welcome pull/merge requests; a few people in this thread already took the source and updated it, I would be happy to update any updates or other changes to help other users.

Hello Ivo, how are?
Just making an observation here:
I've worked a lot with your solver and with some adaptations I've made to it. I can surely affirm here, that it can be easily compiled/adapted to any other version of OpenFOAM. I've compiled it on the most recent versions as a test, and it works just fine.

I'm glad you added a github repository.
Cheers to all

Marco Trofa · January 7, 2021, 13:20

Hi Nicolao,

could you please share your adapted most updated version of openFoamEHD?
It would be very usefull to compare the results with the original version.

Thanks

Nicolao · January 7, 2021, 14:51

Quote:

Originally Posted by Marco Trofa

Hi Nicolao,

could you please share your adapted most updated version of openFoamEHD?
It would be very usefull to compare the results with the original version.

Thanks

Hi Marco,

the last version of OpenFOAM I've tested was the OF1812

I attached here this adapted version. I haven't checked the numerical results, but qualitatively it seemed alright for me.

Hope it works well with you.

Cheers

res23 · February 27, 2021, 19:30

Hey Nicolao, many thanks for providing your updated version of the EHD solver.

When I try to compile the solver, I get a few errors that seem to be related to CorrectPhi.H. I'm very new with openFoam so I'm not sure what I'm doing wrong eactly. I would appreciate if you could provide any help.

Thanks!

Nicolao · March 1, 2021, 20:51

Quote:

Originally Posted by res23

Hey Nicolao, many thanks for providing your updated version of the EHD solver.

When I try to compile the solver, I get a few errors that seem to be related to CorrectPhi.H. I'm very new with openFoam so I'm not sure what I'm doing wrong eactly. I would appreciate if you could provide any help.

Thanks!

Hello Res,

can you share the error message and tell me what version of OpenFOAM are you trying to compile it?

Best

res23 · March 2, 2021, 08:27

Hey Nicolao,

I managed to solve the problem, thanks for offering help though!

I was wondering if you have by any chance any cases that I could use to test my solver. I tried using the original cases provided by Ivooo but they don't seem to be working.

I am using OpenFOAM 1812 and using the files you attached in a previous post here (of1812.zip).

Many thanks!

Nicolao · April 13, 2021, 16:10

Quote:

Originally Posted by res23

Hey Nicolao,

I managed to solve the problem, thanks for offering help though!

I was wondering if you have by any chance any cases that I could use to test my solver. I tried using the original cases provided by Ivooo but they don't seem to be working.

I am using OpenFOAM 1812 and using the files you attached in a previous post here (of1812.zip).

Many thanks!

Hi res23, how are you?
Sorry for the late response. Things are kinda crazy around here...

I am attaching here a case that I think it can work for you (or anyone else interested). Its an axisymmetric droplet subjected to an electric field - the classic case that you can find on my paper:

https://www.sciencedirect.com/scienc...77025714001426

I am not sure thou about what kind of conditions I used here, but you can use any other condition you like that the case will probably work just fine.

Hope it helps you out.

Cheers

YHP · May 21, 2021, 08:31

Quote:

Originally Posted by Nicolao

Hi Marco,

the last version of OpenFOAM I've tested was the OF1812

I attached here this adapted version. I haven't checked the numerical results, but qualitatively it seemed alright for me.

Hope it works well with you.

Cheers

Thank you very much for sharing. I successfully compiled your solver, but when I verified the droplet deformation in the electric field, I could only observe the long deformation of the droplet. I hope you can answer my question.

gsam · May 9, 2022, 14:25

If someone has the latest version of this solver could you upload it here? It will be much helpful for me... I'm kind of in an urgency

chi thuc · August 3, 2022, 00:50

Hi Ivoo and you guys, Thank for your upload file code on GitHub. Im totally newbie in CFD, I download interFoamEHD and interFoamEHDMR solver one month ago to try simulate some cases like paper of Ivoo. I complied both of solvers without prolems and can simulate the cases related to interFoamEHD solver. ' I use OpenFoam2.1.1, Ubuntu18.04.5, Window 10 '
However, i facing prolem like this,
1. When i try run solver inteFoamEHDMR for case case ' dielectric-dielectric_Multiregion ' . The error ' can not upload groovyBC ' and ' can not upload libcoxvoinov.so ' appeared on terminal. So how can I fix this prolem

2. I try to run blockMesh of case ' dielectric-dielectric_Multiregion ', but in the polymesh folder have two diffrent file blockMeshDict and blockMeshDict.bak with different geometry. So, how can I run blockMesh to get full geometry multi region solid-fluid-fluid like paper.

I hope will receive your help soon. I really appreciate for you help. Thank you guys so much.

June 28, 2019, 01:37		#41
Pavithra Member Join Date: Apr 2019 Location: India Posts: 81 Rep Power: 7	Hello Everyone, I am new to OpenFoam. Based on the solver posted in this thread, I have combined electrical equations with buoyantBoussinesqPimpleFoam. My intention is to have a solver for single phase EHD flow. I am able to compile the solver successfully. Could some one please suggest me a good benchmark case to validate the solver. Thank You. With Thanks. Pavithra. Last edited by Pavithra; June 28, 2019 at 01:38. Reason: adjust spacing

June 28, 2019, 10:39		#42
Nicolao New Member Nicolao Join Date: Jul 2012 Location: São Paulo, Brazil Posts: 13 Rep Power: 14	Dear Pavithra, I guess you probably already seen this paper (from Herrera), http://gfs.sourceforge.net/papers/ehd.pdf but in there they show several validation cases for each equation of the leaky dielectric model (for EHD flows). Of course they use a two-phase approach, but for the Poisson and the charge conservation equations you can use the same cases they used there. "My intention is to have a solver for single phase EHD flow" Now, specifically for your case I guess you already have something in mind. =) Cheers Pavithra likes this.

October 22, 2019, 10:04	water impact surface under electrical field	#47
XinXin New Member XinZhou Join Date: Apr 2018 Posts: 11 Rep Power: 8	Hello，everyone. I am simulating a water droplet impact surface under electrical field based on the interFoamEHD solver which designed by ivook.My model is 2D-Axisymmetric. I find that the charge density distribution is weird, there is some charge appearing above the droplet along the axis. I am puzzled by this problem. The charge density distribution is related to the velocity Magnitude. I think the charge conversation equation is calculated for all computation field,both water and air in the interFoam solver. However, for water, it will get polarized and ohmic migration of charged ions. For air,it can be regarded as a perfect dielectric, there is no charge convection. I think we should use different physical model for water and air. But can it be realized in interfoam. I am comfused by the interfoam solver. How does this sovler distinguish the phase1 and phase2? I am reading the twoPhaseMixture.H, incompressibleTwoPhaseMiture.H ,but still get stuck.Could you please give me some advices? charge distribution 2.png charge distribution.png equation.jpg

February 23, 2020, 18:11	intermhdfoam vs interehdfoam for lorentz force in arc welding scenarios	#49
tom_flint2012 Member Thomas Flint Join Date: Jan 2016 Posts: 60 Rep Power: 10	Hi All, Thanks for sharing the code. I hope this is an appropriate thread to ask my questions as they are a bit more general than the electrohydrodynamics solver discussed here, and I think its more of a magnetohydrodynamics question really, but I was wondering if this thread was relevant to my problem. Basically I currently have a VOF solver fro simulating multi-component metallic melting and vapourisation, with energy equations, source terms in the NS equations for solidification etc. I use this to simulate power beam welding applications and arc welding applications, but I would really like to incorporate the electric charge transport and magnetic effects to look at the effect of the Lorentz force on the molten pool behaviour. I have a mhd + interfoam solver that I constructed based on this work: https://www.researchgate.net/publica...etallic_liquid As I understand it I could either use an mhd approach or an ehd approach, and I was wondering if people had any thoughts on which may be appropriate. Are the source terms in the NS equations in the interehdfoam solver equivalent to the JxB lorentz force term in the mhd (and my intermhdfoam solver)?? But I would like to be able to apply a alternating electric potential as a boundary condition. Would you guys have any pointers on the best way for me to progress? I could work out the time varying B field using Amperes law based on an alternating current density (I think), and put this as my B boundary condition in my mhd+interfoam solver. But I thought I'd try and get some advice from people more experienced than me before I go down this route to try and save myself some time in the long run. Thanks in advance for any advice you guys can offer. I'm very grateful. All the best, Tom

November 18, 2020, 07:52	Updated link to source	#50
Ivooo New Member Ivo Join Date: Feb 2012 Posts: 26 Rep Power: 14	I shared my code a while back, but the link seems to stop working every now and then. So I did what I should have done a long time ago, create a github repository. You can find it here (but still based on OF2.1.1): https://github.com/iroghair/openFoamEHD Although I am not focussing my research on EHD specifically at this time, I welcome pull/merge requests; a few people in this thread already took the source and updated it, I would be happy to update any updates or other changes to help other users. Nicolao and LiedtkeJ like this.

January 8, 2020, 14:37		#48
mcfdma Member X Join Date: Jan 2019 Posts: 63 Rep Power: 7	Hello, Firstly, thank you for the solver. I am trying to update the solver on OF6. However, the droplet breakup is significantly different to the original solver. I have uploaded my solver on Google Drive due to exceeding file size [1]. Could someone point out where exactly am I going wrong in terms of modifying my solver. I have attached the screenshot for the last time step with OF231 and OF6. OF231 last_OF231.png OF6 last_OF6.jpg [1] https://drive.google.com/open?id=1Jd...KiYCL0Wc3JHxyi

January 7, 2021, 13:20		#52
Marco Trofa New Member Join Date: Jan 2021 Location: Italy Posts: 1 Rep Power: 0	Hi Nicolao, could you please share your adapted most updated version of openFoamEHD? It would be very usefull to compare the results with the original version. Thanks

February 27, 2021, 19:30		#54
res23 New Member Join Date: Feb 2021 Posts: 3 Rep Power: 5	Hey Nicolao, many thanks for providing your updated version of the EHD solver. When I try to compile the solver, I get a few errors that seem to be related to CorrectPhi.H. I'm very new with openFoam so I'm not sure what I'm doing wrong eactly. I would appreciate if you could provide any help. Thanks!

March 2, 2021, 08:27		#56
res23 New Member Join Date: Feb 2021 Posts: 3 Rep Power: 5	Hey Nicolao, I managed to solve the problem, thanks for offering help though! I was wondering if you have by any chance any cases that I could use to test my solver. I tried using the original cases provided by Ivooo but they don't seem to be working. I am using OpenFOAM 1812 and using the files you attached in a previous post here (of1812.zip). Many thanks!

May 9, 2022, 14:25	EHD foam latest version	#59
gsam New Member ghanashyam k c Join Date: May 2022 Posts: 8 Rep Power: 4	If someone has the latest version of this solver could you upload it here? It will be much helpful for me... I'm kind of in an urgency

August 3, 2022, 00:50	interFoamEHDMR prolem	#60
chi thuc New Member chi thuc Join Date: Jun 2022 Posts: 3 Rep Power: 4	Hi Ivoo and you guys, Thank for your upload file code on GitHub. Im totally newbie in CFD, I download interFoamEHD and interFoamEHDMR solver one month ago to try simulate some cases like paper of Ivoo. I complied both of solvers without prolems and can simulate the cases related to interFoamEHD solver. ' I use OpenFoam2.1.1, Ubuntu18.04.5, Window 10 ' However, i facing prolem like this, 1. When i try run solver inteFoamEHDMR for case case ' dielectric-dielectric_Multiregion ' . The error ' can not upload groovyBC ' and ' can not upload libcoxvoinov.so ' appeared on terminal. So how can I fix this prolem 2. I try to run blockMesh of case ' dielectric-dielectric_Multiregion ', but in the polymesh folder have two diffrent file blockMeshDict and blockMeshDict.bak with different geometry. So, how can I run blockMesh to get full geometry multi region solid-fluid-fluid like paper. I hope will receive your help soon. I really appreciate for you help. Thank you guys so much.