<!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>