Hi, I am back with my scattering problem. Setting the Dropping_Tolerance to some small values, I obtained convergence for lower (but still high) values of Nb_fill. The results are good using second order elements on a finer mesh (see attached pdf), but spending almost 1GB of RAM and 1 hour.

I don't manage to obtain convergence for other values of the parameter Freq of my .pro files (for example, Freq = 1e7 instead of 1e8). So I am still interested in any advice which could make the resolution easier.

<br><br>Thanks,<br><br>Laurent Bernard<br><br>PS: the .geo .par and .pro files are attached, here are the corresponding processing messages:<br><br>P r o c e s s i n g . . .<br>Operation : Generate[A]<br>Solver    : Loading parameter file '/home/lbernard/getdp/sphere_homo/solver.par'

<br>Info      : Setting System {A,b} to zero<br>Resources : cpu 14.5586 s<br>Operation : Solve[A]<br>Solver    : No scaling of system of equations<br>Solver    : RCMK algebraic renumbering<br>Solver    : N: 62588, NZ: 2629288, BW max/avg: 112/42, SW max: 62315

<br>Resources : cpu 16.8311 s<br>Solver    : ILUTP (Float, fill-in = 1500)<br>Solver    : N: 62588, NZ: 96125406, BW max/avg: 3000/1535, SW max: 10681<br>Resources : cpu 3432.68 s<br>Solver    : Generalized Minimum RESidual (GMRES)

<br>    1  2.4992055e+02  1.0000000e+00<br>    2  ...<br><br>   34  1.4337349e-08  5.7367627e-11<br>Solver    : 34 Iterations / Residual: 1.43364e-08<br>Resources : cpu 3544.23 s<br>Operation : SaveSolution[A]<br>Resources : cpu 

3544.41 s<br>E n d   P r o c e s s i n g<br><br><br><br><br><div><span class="gmail_quote">On 12/1/06, <b class="gmail_sendername">laurent bernard</b> <<a href="mailto:ldbernard@gmail.com">ldbernard@gmail.com</a>> wrote:

<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hello, I am trying to compute the field scattered by an homogeneous sphere exposed to an incident  electromagnetic plane wave.

<br>The unknown is the scattered field (electric field for example).<br>The source is the incident field in the volume of the sphere. (See for example 

P.D.Ledger's thesis §5.3 for the variational formulation of the problem <a href="http://www.sam.math.ethz.ch/%7Eledger/publication.html" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">http://www.sam.math.ethz.ch/~ledger/publication.html

</a>). PML are used to troncate the computational domain.

Only one-quarter of the sphere is represented (taking advantage of the symmetries). The mesh is generated using Gmsh (netgen algorithm with optimization) (18159 elements) Problems: Nb_Fill must be set to 3000 in order the obtain convergence. For lower Nb_Fill, the solver faces break-down. In consequence, the resolution is highly RAM and time consuming. So i don't manage to solve the problem on a finer mesh, and i can't use second order elements.

<br>The solution shows unexpected irregularities (see attached pdf file, which shows the components of the total electric field on a segment passing through the sphere).<br><br>Questions:<br>Is there any mistake in the .pro files? In the geommetry? (see attached files)

<br>If not, what should i do to improve the results? to solve the problem using less RAM?<br><br>thank you,<br>best regards<br><span class="sg"><br>Laurent Bernard<br><br><br>

</span><br clear="all"></blockquote></div><br>