Elmer Discussion Forum

Hello everyone,
I am having some trouble solving a magnetodynamics problem.
The thing is that it appears the message of solution trivially zero, even though previous simulations worked well because I had done the modifications said in another post about this topic. I post my case file in case anyone can see the error and help me.
Moreover, another difference with previous cases that I have simulated, the mesh is bigger and in the solver log appears the message of GCR: failed to allocate memory of size: 1868765 x 500.
I really don't know what to do to solve this problems.
Thank you
Guillem

Hi

I don't see any boundary conditions for the vector potential as in
https://github.com/ElmerCSC/elmerfem/bl ... e/case.sif

GCR allocates vectors for all 500 iterations. Rather try BiCGStabl or Idrs and then you can use much more iterations too. The BCs may help in convergence.

-Peter

Hi Peter,
I really appreciate your response.
I have added BCs to my case and the objective is to introduce current into the wire (the case is a cube of air and thisconductor in the middle), but I still get solution trivially zero, with BiCGStabl.
If yo know what my error is it woul be really helpful.
Thank you and sorry for the inconveniences.
Guillem

I believe he is saying that the Electric Infinity BC = True is defined in the StatElecSolver
but not in the WhitneyAVSolver, the outer BC is the WhitneyAVSolver is like:

Boundary Condition 6
Name = "AirEnd"
Target Boundaries(1) = 6

AV re {e} 1 = Real 0.0
AV re {e} 2 = Real 0.0

AV im {e} 1 = Real 0.0
AV im {e} 2 = Real 0.0
End

And how do you specify to the Outer BC that it is an Infinite BC, like in the Electrostatic Solver?

My understanding is that is what AV = 0 is.

Perfect, thank you!
Now it works fine , but only introducing current density in real numbers.
Does anyone know how to introduce imaginary components of current density in the WhitneyAVSolver?
Guillem

Hi

The "infinite bc" is a Robin type of boundary condition that may be derived when you assume that your field goes down as in spherical symmetry. It is a useful approximation for diffusion type of equations far from the source.

Now for AV solver you cannot use this approximation. This is because 1) the \curl\curl operator behaves differently that the \div\grad operator (laplace) and 2) for the edge element formulation you don't even have any dof normal to the surface that we could set.

Hence, what Kevin says is the best you can do.

-Peter

Perfect!
It works perfectly now, thank you so much.
Now I just need to discover if it is possible to introduce current with real and imaginary components.
I really appreciate your answers!
Guillem