As discussed here:
You just need to set the material parameters, using a function like SetSiliconParameters.
The scripts in python_packages are not documented in the manual. You would have to look at the code. I will work on adding docstrings in the future.
https://github.com/devsim/devsim/issues/146
This function solves the Potential equation in a semiconductor region, assuming equilbrium.
This function sets the boundary condition at a contact for the potential equation in equilibrium.
Setups up the Potential, Electron, and Hole equations for drift-diffusion simulation.
Sets up an ideal ohmic contact for drift diffusion simulation.
Set the material parameters for your material. The existing equations should work for a single region device.
Please check the discussion here:
https://forum.devsim.org/t/simulation-of-simple-2d-sic-diode/328/2
One approach is to use the Step function. However, if you want to have an abrupt junction that is independent of meshing effects, you can have a P region and N region with an interface between them. Then the
CreateSiliconSiliconInterface function would be useful.