SIF file in PIEZO example

General discussion about Elmer
annier
Posts: 1168
Joined: 27 Aug 2013, 13:51
Antispam: Yes

Re: SIF file in PIEZO example

Post by annier »

Hi Mancf,
You may change the source files of the stress solve and static current conduction in accordance to the mathematical requirement and compile them with elmerf90 wrapper. Then you can access them via user defined solver (Elmer Programmers Manual).

Yours Sincerely,
Anil KUnwar
Anil Kunwar
Faculty of Mechanical Engineering, Silesian University of Technology, Gliwice
mancf
Posts: 11
Joined: 15 Jul 2015, 18:58
Antispam: Yes

Re: SIF file in PIEZO example

Post by mancf »

Hi Annier,
Up to now I have not had to modify the code... Taking into account that I have never programmed in fortran, this can take quite a while for me. I have always used the compiled version for Windows.
In any case, I will try to do something with it during the next weeks.
Thanks in any case,
mancf
annier
Posts: 1168
Joined: 27 Aug 2013, 13:51
Antispam: Yes

Re: SIF file in PIEZO example

Post by annier »

Hi mancf,
In Elmer there is a template how you develop solvers in accordance to the structure of actual PDE.
1. Please see about chapter 7 about Partial Differential Equations in Elmer Programmers Manual (http://www.nic.funet.fi/pub/sci/physics ... torial.pdf). The example for this chapter ishttps://github.com/ElmerCSC/elmerfem/bl ... delPDE.F90.

2. The structure of a solver in Elmer module is described in Elmer Solver Manual (Chapter 18, section 18.4).http://www.nic.funet.fi/pub/sci/physics ... Manual.pdf
Now, you can inspect the contents of stress solver and static current conduction solver in accordance to the aspects descibed in the solver manual and can do modification inaccordance to the governing Partial Differential Equation (PDE).


Yours Sincerely,
Anil Kunwar
Anil Kunwar
Faculty of Mechanical Engineering, Silesian University of Technology, Gliwice
mancf
Posts: 11
Joined: 15 Jul 2015, 18:58
Antispam: Yes

Re: SIF file in PIEZO example

Post by mancf »

Hi again,
Sorry. Quite a long time ago, but I could not dedicate time to it.
What I realised is that the problem is the same in other coupled problems using the electric field. If you run the electrokinetics example, you get the same problem (it does not work). There is no effect of the electric field in the advection of the concentration. I guess that it worked sometime ago (the results are shown in the "ElmerTutorials_nonGUI.pdf" file, page 72), but there must have been a change in the code.
Maybe it could be possible to detect the problem easily by comparing the code with previous versions.
Regards,
mancf
Post Reply