[aviatorr]

Dear Sweet People,

I have developed a EHDsimpleFoam code fused from simpleFoam and electrostaticFoam.

The development encompassed 4 things:
1) adding coulomb force to UEqn, as seen in attachments
2) writing EEqn based on electrostatic foam, as seen in attachments
3) adding #include EEqn in the main file, as seen in attachments
4) fused the CreateFields.C from both simpleFoam and electrostaticFoam.

The Pressure equations I left unchanged.
I am running the sim with RAS kepsilon, turbulence off (so laminar pretty much).

The validation case that I am doing is at zero background velocity. All velocity is induced by ionic wind. The setup is half-domain, with symmetryPlane. It has emitter and collector, both being semi-circles (cause half-domain).

THE PROBLEM:
The velocity field seems nice, as attached. But the pressure field..... God Please.... what is this ??? Even the forceCoeffs outputs near 0 values for Cd at the big semi-circle.... I have been trying to figure it out why that happended but cant find a solution. If I run the same sim without charge, but set the inlet velocity to the value that was induced by ionic wind, then the pressure looks alright.

QUESTION:
How to fix the pressure ??