CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Post-Processing

Finding my current.

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   November 6, 2013, 19:31
Default Finding my current.
New Member
Join Date: May 2013
Posts: 9
Rep Power: 6
MechBoom is on a distinguished road
My problem is that I know what I need to do, I'm just no sure how to code it.

Here's the integral I'm working with:

I=-ρ b∫ ∂V/∂n dl

rho and b are known constants. ∂V/∂n is the partial of Voltage normal to the plane and dl is obviously my differential length. With this integral I'm finding the corona current around a wire with fairly high current.

I'm basically having issues with two parts. The first is figuring out how to find the partial of voltage normal to my plane (wire). I have solved for the voltage field throughout my mesh but don't know how to single out the values around the wire. I know that the values at the wire itself all are normal to the plane so all I need to find are the magnitudes of each mesh cells voltage at the wire.

Any help is much appreciated!!
MechBoom is offline   Reply With Quote

Old   November 7, 2013, 04:44
Senior Member
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 15
Bernhard is on a distinguished road
The integral that you wrote here, is this a line integral?

For the surface normal gradient, you have to snGrad() functions available.
Bernhard is offline   Reply With Quote

Old   November 7, 2013, 18:01
New Member
Join Date: May 2013
Posts: 9
Rep Power: 6
MechBoom is on a distinguished road
Since I'm only analyzing a 2D surface, it is essentially a line integral. I'm just integrating around the circumference of a circle (the wire).

And thanks for the snGrad() tip, I'm looking at that function now.
MechBoom is offline   Reply With Quote

Old   November 8, 2013, 15:21
New Member
Join Date: May 2013
Posts: 9
Rep Power: 6
MechBoom is on a distinguished road
Basically all I've done so far is modify the electrostatic solver to solve for the E-field throughout the duct.

    while (runTime.loop())
        Info<< "Iteration = " << runTime.timeName() << nl << endl;

        solve(fvm::laplacian(phi) + rho/epsilon0);

        rhoFlux = -k*mesh.magSf()*fvc::snGrad(phi);

        solve(fvm::ddt(rho) + fvm::div(rhoFlux, rho));

    volVectorField E("E",-fvc::grad(phi)); //Solve For E-Field
now that I've got the e-field, I need to find the corona current around the wire. The integral above does that. So far I'm diggin into the snGrad() function and can't figure out how to get it work. My programming knowledge is somewhat limited. I've been teaching my self C++ and OpenFOAM at the same time. The reason is for my undergrad research. I'm not looking for someone to write my code! I just need some assistance in figuring out:

a) how to find the magnitude of my e-field normal to the wire. As stated above I just need the values of the e-field at the wire since all E-vectors at the wire are normal.

b) how to find the length of each mesh cell at the wire. I could always hard wire a value in for this, but that's not exactly what I would like to do.

Any help that anyone could provide would be GREATLY appreciated. I've been stuck at this point for a month now, and don't know where else to turn.
MechBoom is offline   Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Current density visualisation (PEM fuel cell add-on module) pchoopanya FLUENT 6 July 14, 2016 15:26
how is parasitic current now? houkensjtu OpenFOAM Programming & Development 18 June 8, 2016 05:13
Turbidity current julien.n FLUENT 0 June 25, 2012 19:12
What is the difference between current time step and current time djing FLUENT 4 May 1, 2012 16:18
dilute gravity current at very flow base doronzo FLUENT 0 October 5, 2011 15:38

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