Charge relaxation problem

Numerical methods and mathematical models of Elmer
Post Reply
PauloFonte
Posts: 3
Joined: 24 Aug 2018, 23:47
Antispam: Yes

Charge relaxation problem

Post by PauloFonte » 25 Aug 2018, 23:04

Dear Elmer experts

I would like to solve a charge relaxation problem. That is
div(sigma*grad(V)) = div(d/dt (epsilon*grad(V)))

The related Elmer solver seems to be the StatCurrentSolver
div(sigma*grad(V)) = d/dt (rho)

Searched the forum and the most relevant discussion I could found was
http://elmerfem.org/forum/viewtopic.php ... =10#p16162
It seems that the time derivative is not included in StatCurrentSolver and there is a suggestion of a modification, but I have no idea how to use it (and don't want to become an expert myself...).

I also dreamed that the rhs term could be calculated in a parallel solver, but couldn't find a way to calculate the time derivative. Also, have no clue about the chances of convergence of such contraption.

Any suggestions welcome. Thanks in advance

raback
Site Admin
Posts: 3321
Joined: 22 Aug 2009, 11:57
Antispam: Yes
Location: Espoo, Finland
Contact:

Re: Charge relaxation problem

Post by raback » 02 Sep 2018, 22:12

Hi

As written before this would be just a few lines of code. We may be soon writing an optimized version of the StatCurrentSolver and could include that term in the process.

-Peter

raback
Site Admin
Posts: 3321
Joined: 22 Aug 2009, 11:57
Antispam: Yes
Location: Espoo, Finland
Contact:

Re: Charge relaxation problem

Post by raback » 18 Sep 2018, 19:46

Hi All,

I'm in the process of writing a new optimized version of StatCurrentSolver. At the same time I could augment it to be more complete. Hence I thought about adding this relaxation term. Here \sigma is electric conductivity and \eps is permittivity and the source term is given.
NewStatCurrent.PNG
NewStatCurrent.PNG (5.93 KiB) Viewed 290 times
Would this solve the charge relaxation problem. Any other ideas while we are at it?

-Peter

PauloFonte
Posts: 3
Joined: 24 Aug 2018, 23:47
Antispam: Yes

Re: Charge relaxation problem

Post by PauloFonte » 18 Sep 2018, 21:43

For me this will be fine.
Don't know how this fits the general philosophy of Elmer, but in some cases the frequency domain version (d/dt -> Iw) may be more useful (and, I guess, simpler to calculate).
Thanks in advance.
Paulo

raback
Site Admin
Posts: 3321
Joined: 22 Aug 2009, 11:57
Antispam: Yes
Location: Espoo, Finland
Contact:

Re: Charge relaxation problem

Post by raback » 09 Nov 2018, 13:14

Hi Paulo

This solver has been there for some time. Would you have some test case to try it out?

-Peter

PauloFonte
Posts: 3
Joined: 24 Aug 2018, 23:47
Antispam: Yes

Re: Charge relaxation problem

Post by PauloFonte » 09 Nov 2018, 21:20

Dear Peter

Sorry, but unfortunately right now I'm too busy with classes to tackle this. Must first learn Elmer...

Hope to be able to work on this in January. Will keep you informed.

Kjnd regards
Paulo

Post Reply