[Getdp] integral quantities

[André Nicolet]

Hello,

I would like to know how to compute integral quantities in order to be able to compute forces, energy, power, etc.

I have modified Magnetostatic.pro this way in the postprocessing :

{ Name wm ; Value { Local { [ SquNorm[{d a}] * nu[] / 2.] ; In Domain ; Jacobian JVol ; } } }

{ Name Wmag ; Value { Integral { [ SquNorm[{d a}] * nu[] / 2.] ; In Domain ; Jacobian JVol ; } } }


{ Name lapfor ; Value { Local { [ Cross [ js[], {d a} ] ] ; In Domain_S ; Jacobian JVol ; } } }

and then in test.pro

Print[ wm, OnElementsOf Domain, File "energy.pos"];

Print[ Wmag, OnRegion Domain, File "en.pos"];

Print[ lapfor, OnElementsOf Domain_S, File "force.pos"];

The first line computes the magnetic energy density and the second one is an attempt to compute the total energy. Unfortunately, it writes zeroes in the file en.pos.

The third line gives the Laplace force density. How to integrate it on the conductors to get the total force ?

How to integrate on the boundary of a region to get the forces by the Maxwell stress tensor or the energy with the Poynting vector ?

Thank you for your help !

********************************************
Prof. André Nicolet
Université d'Aix-Marseille III
Dom. Univ. de Saint-Jérôme
Service 162 
Institut Fresnel (UMR 6133)
Avenue Escadrille Normandie-Niemen
F-13397 Marseille cedex 20
France
tel +33 4 91 28 87 73
fax +33 4 91 28 88 80
secr +33 4 91 28 83 76
mailto:andre.nicolet at fresnel.fr
********************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/getdp/attachments/20011005/3b7d0015/attachment.html>