I’m very sorry that some of my questions may have been too basic, but I do feel that I’m making progress and your replies are really very important to me, so thank you again for your patience in replying!
Sorry for the delay, but I have had a very busy week.
Convergence problems can be tricky to diagnose.
If it is diverging as you go to a higher bias, you may need to reduce the voltage step. If it is diverging near dc, there could be a problem with the mesh.
Is it not converging with the mesh that is included in the projects, or is this one that you had created yourself? Did you change many of the parameters or simulation scripts?
Hello, author:
The problem occurs when I run the bjt_dd.py program to generate the bjt_dd_0.msh file, which at this point in time should not have added the applied bias yet, but just the equations, and there is too much error in the continuity equations.
The mesh file I used was exactly the bjt.msh given in your example, and I changed all the parameters that I had previously modified (copied the code you shared on github).
Here is the bjt.msh, is it possible that the one given in your example is not the final msh file? Or something else?
I got someone else to look at your example for me and found that using your bjt.msh file, I really couldn’t run the bjt_dd.py file, and adding up the number of iterations to 1000 didn’t help convergence.
I didn’t change any parameters in your code.
Can you run the script bjt_dd.py in your spare time? I’ve been stuck on this for quite some time.
Hello ,author:
I saw it as soon as you replied to me, but I’m very sorry I’m in class in the morning and don’t have my computer with me, I’ll try to verify your replies first in the future.
The error did get smaller after I changed it, but with only two iterations,Is this normal?Is this too soon?
In TCAD device simulation, a full drift-diffusion simulation will not converge without an initial guess.
A potential-only simulation, with flat fermi levels, is equivalent to the solution of the full drift-diffusion simulation. This is since the electrons and holes are determined by the equilibrium potential solution.
Potential was solved in the potential-only simulation. When creating it again, it was reset to zero. This was a poor initial guess, even though the electrons and holes were correctly initialized here:
The file format has changed over the years. If you compare the files using a diff program, you can see that some unimportant data was removed from the file format.
You explained it in great detail and thank you very much for your patience.
I asked you earlier about grid refinement and I am very sorry that my description may have been too sketchy at the time.
I experimented many times and changed the sentence you mentioned and found that:
“mincl, maxcl”, these two parameters control should be the grid density between the emitter-base-collector region two by two.
“pdiff” has no effect on the generated grid.
After many attempts I still fail, I think the grid generated at each step looks very reasonable, but when running bjt_refine.py for the second time still “devsim_py3.error: DEVSIM FATAL” occurs.
If you have time could you take a look at what the problem is?
Unfortunately, devsim aborted when it determined there was something structually wrong with the mesh. I would suggest modifying the mesh and trying again.
