[Getdp] Newton Raphson method with a time dependent problem

Olivier Castany castany at quatramaran.ens.fr
Tue Aug 3 20:30:04 CEST 2010


Hello,
 
>         JacNL[ dlambdadT[{T}] * Dof{T} * {d T} , {d T} ]
>         [ lambda[{T}] * Dof{d T} , {d T} ]

This looks good.

> But I must confess I was surprised to see that you have also a "Dof" in the second line. I looked at the "MagnetoStaticsNonLinear<https://geuz.org/trac/getdp/wiki/MagnetoStaticsNonLinear>" example from the Wiki:
> 
>       Galerkin { JacNL[ - dbdh[ {d phi} ] * Dof{d phi} , {d phi} ] ;
>                  In Steel ; Jacobian MyJac ; Integration MyInt ; }
> 
>       Galerkin { [ b[ - {d phi} ] , {d phi} ]  ;
>                 In Steel ; Jacobian MyJac ; Integration MyInt ; }
> 
> There is no Dof used for the non-JacNL part.
> The interesting thing is that in your example it is necessary to have it whereas in the Wiki example it does not work if I use Dof (by replacing "b[-{d phi}]" by "mu_Steel[SquNorm[{d phi }]] * Dof{d phi}").

Strange. Can you check with the opposite sign ?  
Maybe you made a mistake with the sign.
We have : h = - Grad phi
so : b = mu * h = mu * (- Grad phi)

A+