<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
<div>Hi Frederic,</div>
<div><br>
</div>
<div>You can inspire yourself from the magnetostatic demo you get when downloading GetDP.</div>
<div><br>
</div>
<div>There is a complete Newton-Raphson example, i.e. a term for completing the Jacobian matrix is included (JacNL).</div>
<div>You can just remove or comment it for getting a fixed point iteration. The resolution does not change.</div>
<div><br>
</div>
<div>HTH,</div>
<div>Ruth</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<br>
<div>
<div>On 14 Jul 2014, at 18:29, Frederic Trillaud Pighi <<a href="mailto:ftrillaudp@pumas.iingen.unam.mx">ftrillaudp@pumas.iingen.unam.mx</a>> wrote:</div>
<br class="Apple-interchange-newline">
<blockquote type="cite">
<div>
<div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;">
<div>Dear all,<br>
</div>
<div><br>
</div>
<div>I have a non-linear magnetodynamics problem that involves a fixed point iteration. Would someone have an example to share or some recommendation how to implement the algorithm?<br>
</div>
<div><br>
</div>
<div>The problem is the following (A formulation, sig: electrical conductivity, Je: applied current):<br>
</div>
<div><br>
</div>
<div>Curl[nu*Curl(A)]+sig*dA/dt = Je<br>
</div>
<div><br>
</div>
<div>with,<br>
</div>
<div><br>
</div>
<div>sig = f(-dA/dt+Je/sig])<br>
</div>
<div><br>
</div>
<div><br>
Best,</div>
<div><br>
</div>
<div><br>
</div>
<div>Frederic<br>
</div>
</div>
</div>
_______________________________________________<br>
getdp mailing list<br>
<a href="mailto:getdp@geuz.org">getdp@geuz.org</a><br>
http://www.geuz.org/mailman/listinfo/getdp<br>
</blockquote>
</div>
<br>
<div apple-content-edited="true">—<br>
Prof. Ruth V. Sabariego<br>
KU Leuven <br>
Dept. Electrical Engineering ESAT/Electa<br>
<a href="http://www.esat.kuleuven.be/electa">http://www.esat.kuleuven.be/electa</a><br>
<br>
<br>
<br>
<br>
<br>
<br>
</div>
<br>
</body>
</html>