On 8/27/25 02:55, Dmitriy Samosvat wrote:
I am trying to solve the 30-band kp model for a silicon quantum well. I
managed to get the ground state for the electron wave function. However,
it has oscillations. I think a more physical basis for this problem
would be the sines (standing waves) basis (instead of Lagrangian finite
elements). However, when I try to use the sines basis, it gives
incorrect eigenvalues. I am sending you the file with the wave function
and the sines basis. I would be grateful if you could tell me how to fix
my code.
I think this is asking too much of the people on this forum -- you'll
have to debug things yourself first.
But a good starting point would be to take your element, a mesh with one
or perhaps four cells, create a DoFHandler with it, create a solution
vector that corresponds to a function you know, and then output the
solution with DataOut::build_patches() using a high value for the
subdivisions. This way you get to see the solution on a subdivided mesh.
Compare that to what you expect and make sure the shape functions are
correct. Then repeat this by using a DataPostprocessor object that
outputs the gradient of your shape functions. Chances are that the
mistake is in one of these two steps.
Best
W.
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To view this discussion visit
https://groups.google.com/d/msgid/dealii/7eb8984e-bdcf-4789-a120-3c0a3302b869%40colostate.edu.
