CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Programming & Development

Integration issue

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

LinkBack Thread Tools Display Modes
Old   November 4, 2013, 11:45
Default Integration issue
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.

Secondly I'm having problems finding the differential length of my cells. I'm 99.9% sure that my mesh is uniform around the wire, so I'd probably only need the length of one cell.

I've been searching the source code info for any functions to solve for what I'm looking for but, so far I'm coming up dry. Any help that could be provided is much appreciated.
MechBoom is offline   Reply With Quote

Old   November 6, 2013, 22:56
New Member
Join Date: May 2013
Posts: 9
Rep Power: 6
MechBoom is on a distinguished road
I've been working on a modified electrostaticFoam solver. So far I've only added the following:

volVectorField E("E",-fvc::grad(phi)); //Solve For E-Field

Now I'm still trying to figure out how to solve for my integral. I've been looking and think that maybe I need to use nCells at the wire and then divide that by the circumference of my wire. If the mesh is uniform this will give me the dl I need. Having said that, I'm not sure how to actually go about doing it.

This is driving me a little crazy. It's so frustrating when you know what you need to do but not how to do it!

Any help you could give would be greatly appreciated!
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
How are the integration constants for the JANAFthermo class derived? tomf OpenFOAM Running, Solving & CFD 2 May 15, 2015 07:47
CFD-Post: Polyline Export at Integration Points minger ANSYS 4 February 22, 2014 10:23
Strange issue while launching ANSYS workbench in CentOS 6.4 Philip_C ANSYS 11 August 29, 2013 06:44
CyclicAMI Issue In OpenFOAM 2.2.0 prasant OpenFOAM Running, Solving & CFD 17 March 16, 2013 03:00
Fluent Integration Numerics RodriguezFatz FLUENT 11 February 21, 2013 07:35

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