## Display Body Force

Post processing utility for Elmer
Naorgk
Posts: 2
Joined: 19 Aug 2015, 13:40
Antispam: Yes

### Display Body Force

Hello everyone,

I've recently started using Elmergui to solve magneto dynamic problems, for this I have defined a current density as a body force.
I was wondering if there is a way to display the body force in ElmerPost to see if it is defined correctly.

mzenker
Posts: 1945
Joined: 07 Dec 2009, 11:49
Location: Germany

### Re: Display Body Force

Hi,

as the current density is a vector, you either have to calculate its magnitude and display that one (you need to enter some math commands for that), or display it as vector - but I don't remember how that worked. There is a ElmerPost help which explains many things.
If all else fails, you can change the output file extension from .ep to .vtu in the sif file, rerun the simulation and use ParaView as postprocessor.

HTH,

Matthias

jay
Posts: 14
Joined: 05 Jun 2015, 12:46
Antispam: Yes

### Re: Display Body Force

hi Naorgk

If your problem is still not solved then try new release of Elmer. The current density that you define as body force can be seen in the post processing. Thanks to CSC team.

For eg. you can see the attached screen shots of my post processed results in paraview.
Attachments
current density_2.jpg (64.21 KiB) Viewed 4469 times
current density1.jpg (35.92 KiB) Viewed 4469 times

jay
Posts: 14
Joined: 05 Jun 2015, 12:46
Antispam: Yes

### Re: Display Body Force

So with refence to my previous comment, I have some queries about the results obtained for current density.
For the model attached in my prvious comment I have used body force as:

Current Density 1 = Variable coordinate,
Real MATC "-(6.25e6)*tx(1)/sqrt((tx(1)*tx(1))+(tx(0)*tx(0)))"
Current Density 2 = Variable coordinate,
Real MATC "(6.25e6)*tx(0)/sqrt((tx(1)*tx(1))+(tx(0)*tx(0)))"

So here the absolute value is 6.25e6 A/m^2

Now according to post processed results, in one of clipped plan I can see in some patches of coil cross section, current densities are above 9E6 A/m^2.
So my question is paraview showing me absolute value or value of conponent of particular tangential vector of current in that place or current in the mesh element their.?

As per my understanding the for appied current, its direction in each element will be as shown in the attached image. So if i clip the model from some place, what quantity of current will paraview show me ?
Personally I wish to see absolute value of current which should be uniformly distributed over entire cross section of coil. How can it be acheived?
The simulation is carried out for static model.

I will be glad if some one can help me in understanding the nformation obtained in Paraview.

With Kindest Regards
jay panchal
Attachments
current direction in each first order element.jpg (48.31 KiB) Viewed 4460 times

pavel
Posts: 34
Joined: 17 Apr 2014, 17:08
Antispam: Yes

### Re: Display Body Force

Hi jay,

Have you tried to set keyword Discontinuous Galerkin = True as it is recommended in the Elmer Models Manual?

BR, Pavel

jay
Posts: 14
Joined: 05 Jun 2015, 12:46
Antispam: Yes

### Re: Display Body Force

pavel wrote:Hi jay,

Have you tried to set keyword Discontinuous Galerkin = True as it is recommended in the Elmer Models Manual?

BR, Pavel
Hi Pavel