[mcfdma]

Hi.

I'm working on a two-phase flow simulation and have modified the interFoam solver to include the electrostatics component to get an electrohydrodynamic solver (similar to http://openfoamwiki.net/index.php/Co...ectromagnetics)

I'm calculating the electric field using the gradient of electric potential given by the equation , which in my header file looks like:

E = -fvc::grad(Ue);

For my simulation, I'm interested in finding the normal component of E inside and outside the interface.


I've tried Paraview, where I applied contour to get the isosurface followed by normal glyphs. This gave me the normal component of E outside of the interface, and I was unable to get the normal component of E inside. Also, this isn't a quick process, and I want to know how I can createFields for both the normal components in my solver.

What I plan to understand is the following:
1. How to acquire(&code) the surface normal vector. This is what I've done (can be wrong, and if it is, do let me know).

	surfaceVectorField phiNormal
	(   mesh.Sf()/mesh.magSf()//(fvc::reconstruct(phi))/mag(fvc::reconstruct(phi))
	);

2. How to use this surface normal vector to find the 2 normal components.

I appreciate your time and would be grateful if I can get some assistance on this area.

[mcfdma]

I have attached an image to explain my question better. Normal inside and outside of the surface.
1200px-Normal_vectors2.svg.jpg