<!DOCTYPE HTML PUBLIC "-//W3C//DTD W3 HTML//EN">
<HTML>
<HEAD>

<META content=text/html;charset=iso-8859-1 http-equiv=Content-Type>
<META content='"MSHTML 4.72.3110.7"' name=GENERATOR>
</HEAD>
<BODY bgColor=#ffffff>
<DIV><FONT face="Courier New" size=2>
<P>Hello,</P>
<P>I would like to know how to compute integral quantities in order to be able 
to compute forces, energy, power, etc.</P>
<P>I have modified Magnetostatic.pro this way in the postprocessing :</P>
<P>{ Name wm ; Value { Local { [ SquNorm[{d a}] * nu[] / 2.] ; In Domain ; 
Jacobian JVol ; } } }</P>
<P>{ Name Wmag ; Value { Integral { [ SquNorm[{d a}] * nu[] / 2.] ; In Domain ; 
Jacobian JVol ; } } }</P>
<P></P>
<P>{ Name lapfor ; Value { Local { [ Cross [ js[], {d a} ] ] ; In Domain_S ; 
Jacobian JVol ; } } }</P>
<P>and then in test.pro</P><FONT face="Courier New" size=2>
<P>Print[ wm, OnElementsOf Domain, File "energy.pos"];</P>
<P>Print[ Wmag, OnRegion Domain, File "en.pos"];</P>
<P>Print[ lapfor, OnElementsOf Domain_S, File "force.pos"];</FONT></P>
<P></FONT><FONT color=#000000 face="Courier New" size=2>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.</FONT></P>
<P><FONT color=#000000 face="Courier New" size=2></FONT><FONT 
face="Times New Roman"></FONT><FONT face="Courier New" size=2>The third line 
gives the Laplace force density. How to integrate it on the conductors to get 
the total force ?</FONT></P>
<P><FONT face="Courier New" size=2>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 ?</FONT></P>
<P><FONT face="Courier New" size=2>Thank you for your help !</FONT></P></DIV>
<DIV><FONT color=#000000 
size=2>********************************************<BR>Prof. André 
Nicolet<BR>Université d'Aix-Marseille III<BR>Dom. Univ. de 
Saint-Jérôme<BR>Service 162 <BR>Institut Fresnel (UMR 
6133)<BR>Avenue Escadrille Normandie-Niemen<BR>F-13397 Marseille cedex 
20<BR>France<BR>tel +33 4 91 28 87 73<BR>fax +33 4 91 28 88 80<BR>secr +33 4 91 
28 83 76<BR><A 
href="mailto:andre.nicolet@fresnel.fr">mailto:andre.nicolet@fresnel.fr</A><BR>********************************************</FONT></DIV></BODY></HTML>