luke
/
ngspice
mirror of https://git.code.sf.net/p/ngspice/ngspice
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

forgot some parameters that depend on "area" and T, but are needed as 

model variables in load
This commit is contained in:
Markus Mueller 2020-07-05 18:03:23 +02:00
parent 320a66c0a4
commit d7dd26880c
3 changed files with 28 additions and 15 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
src/spicelib/devices/hicum2/hicum2defs.h
View File

@ -140,11 +140,16 @@ typedef struct sHICUMinstance {
dual_double HICUMvds_t;
dual_double HICUMvpts_t;
//variables that depend on "area" and "m" but not on temperature
//model variables that depend on "area" and "m" and are needed in load
double HICUMcjci0_scaled;
double HICUMcjcx0_scaled;
double HICUMicbar_scaled;
double HICUMcbepar_scaled;
double HICUMcbcpar_scaled;
double HICUMcth_scaled;
double HICUMrth_scaled;
double HICUMrci0_scaled;
double HICUMibets_scaled;
double HICUMrbi;
double HICUMiavl;

22
src/spicelib/devices/hicum2/hicumL2.cpp
View File

@ -720,7 +720,7 @@ HICUMload(GENmodel *inModel, CKTcircuit *ckt)
FFitf_ick = itf/I_CK;
FFdTef = tef0_t*exp(model->HICUMgtfe*log(FFitf_ick));
FFdQef = FFdTef*itf/(1+model->HICUMgtfe);
if (here->HICUMicbar_scaled<0.05*(model->HICUMvlim/model->HICUMrci0)) {
if (here->HICUMicbar_scaled<0.05*(model->HICUMvlim/here->HICUMrci0_scaled)) {
FFdVc = 0;
FFdVc_ditf = 0;
} else {
@ -852,7 +852,7 @@ HICUMload(GENmodel *inModel, CKTcircuit *ckt)
std::function<duals::duald (duals::duald, duals::duald, duals::duald)> calc_ibet = [&](duals::duald Vbiei, duals::duald Vbpei, duals::duald T){
//Tunneling current
duals::duald ibet;
if (model->HICUMibets > 0 && (Vbpei <0.0 || Vbiei < 0.0)){ //begin : HICTUN
if (here->HICUMibets_scaled > 0 && (Vbpei <0.0 || Vbiei < 0.0)){ //begin : HICTUN
duals::duald pocce,czz, cje0_t, vde_t, ibets_t, abet_t;
double T_dpart = T.dpart();
ibets_t = here->HICUMibets_t.rpart;
@ -897,7 +897,7 @@ HICUMload(GENmodel *inModel, CKTcircuit *ckt)
// double v_bord_r;
// double volatile v_bord_r1;
// int use_aval;
// if ((model->HICUMfavl > 0.0) && (model->HICUMcjci0 > 0.0)) {
// if ((model->HICUMfavl > 0.0) && (model->HICUMcjci0_scaled > 0.0)) {
// use_aval = 1;
// } else {
// use_aval = 0;
@ -1196,24 +1196,24 @@ HICUMload(GENmodel *inModel, CKTcircuit *ckt)
for (; model != NULL; model = HICUMnextModel(model)) {
// Model_initialization
int selfheat = ((model->HICUMflsh > 0) && (model->HICUMrthGiven) && (model->HICUMrth > 0.0));
int selfheat = ((model->HICUMflsh > 0) && (model->HICUMrthGiven) && (here->HICUMrth_scaled > 0.0));
int nqs = ((model->HICUMflnqs != 0 || model->HICUMflcomp == 0.0 || model->HICUMflcomp == 2.1) && (model->HICUMalit > 0 || model->HICUMalqf > 0));
// Depletion capacitance splitting at b-c junction
// Capacitances at peripheral and external base node
C_1 = (1.0 - model->HICUMfbcpar) *
(model->HICUMcjcx0 + here->HICUMcbcpar_scaled);
(here->HICUMcjcx0_scaled + here->HICUMcbcpar_scaled);
if (C_1 >= here->HICUMcbcpar_scaled) {
cbcpar1 = here->HICUMcbcpar_scaled;
cbcpar2 = 0.0;
//cjcx01 = C_1 - here->HICUMcbcpar_scaled;
//cjcx02 = model->HICUMcjcx0 - cjcx01; //not needed herein
//cjcx02 = model->HICUMcjcx0_scaled - cjcx01; //not needed herein
}
else {
cbcpar1 = C_1;
cbcpar2 = here->HICUMcbcpar_scaled - cbcpar1;
//cjcx01 = 0.0;
//cjcx02 = model->HICUMcjcx0; //not needed herein
//cjcx02 = model->HICUMcjcx0_scaled; //not needed herein
}
// Parasitic b-e capacitance partitioning: No temperature dependence
@ -1239,7 +1239,7 @@ HICUMload(GENmodel *inModel, CKTcircuit *ckt)
}
// Turn avalanche calculation on depending of parameters
if ((model->HICUMfavl > 0.0) && (model->HICUMcjci0 > 0.0)) {
if ((model->HICUMfavl > 0.0) && (here->HICUMcjci0_scaled > 0.0)) {
use_aval = 1;
} else {
use_aval = 0;
@ -2491,10 +2491,10 @@ HICUMload(GENmodel *inModel, CKTcircuit *ckt)
hicum_diode(Temp,here->HICUMiscs_t,model->HICUMmsc, Vsici, &ijsc, &ijsc_Vsici, &ijsc_Vrth);
// Self-heating calculation
if (model->HICUMflsh == 1 && model->HICUMrth >= MIN_R) {
if (model->HICUMflsh == 1 && here->HICUMrth_scaled >= MIN_R) {
pterm = (Vbiei-Vbici)*it + (here->HICUMvdci_t.rpart-Vbici)*iavl;
pterm_dT = (Vbiei-Vbici)*it_dT + (here->HICUMvdci_t.rpart-Vbici)*iavl_dT + here->HICUMvdci_t.dpart*iavl;
} else if (model->HICUMflsh == 2 && model->HICUMrth >= MIN_R) {
} else if (model->HICUMflsh == 2 && here->HICUMrth_scaled >= MIN_R) {
pterm = (Vbiei-Vbici)*it + (here->HICUMvdci_t.rpart-Vbici)*iavl + ibei*Vbiei + ibci*Vbici + ibep*Vbpei + ijbcx*Vbpci + ijsc*Vsici;
pterm_dT = (Vbiei-Vbici)*it_dT + (here->HICUMvdci_t.rpart-Vbici)*iavl_dT + here->HICUMvdci_t.dpart*iavl +
ibei_dT*Vbiei + ibci_dT*Vbici + ibep_dT*Vbpei + ijbcx_dT*Vbpci + ijsc_Vrth*Vsici;
@ -2703,7 +2703,7 @@ HICUMload(GENmodel *inModel, CKTcircuit *ckt)
// Following code is an intermediate solution (if branch contribution is not supported):
// ******************************************
//if(model->HICUMflsh == 0 || model->HICUMrth < MIN_R) {
//if(model->HICUMflsh == 0 || model->HICUMrth_scaled < MIN_R) {
// I[br_sht] <+ Vrth/MIN_R;
//} else {
// I[br_sht] <+ Vrth/rth_t-pterm;

14
src/spicelib/devices/hicum2/hicumL2temp.cpp
View File

@ -164,6 +164,8 @@ int hicum_thermal_update(HICUMmodel *inModel, HICUMinstance *inInstance, double
double rcx_scaled ;
double rbi0_scaled ;
double rth_scaled ;
double cth_scaled ;
double ibets_scaled ;
// double cjci0_t, vdci_t, vptci_t, cjep0_t, vdep_t, ajep_t, vdcx_t, vptcx_t, cscp0_t, vdsp_t, vptsp_t, cjs0_t, vds_t, vpts_t;
@ -209,6 +211,7 @@ int hicum_thermal_update(HICUMmodel *inModel, HICUMinstance *inInstance, double
ireps_scaled = model->HICUMireps * here->HICUMm;
cjep0_scaled = model->HICUMcjep0 * here->HICUMm;
cbepar_scaled = model->HICUMcbepar * here->HICUMm;
ibets_scaled = model->HICUMibets * area_times_m;
//BC junction
ibcis_scaled = model->HICUMibcis * area_times_m;
cjci0_scaled = model->HICUMcjci0 * area_times_m;
@ -223,12 +226,17 @@ int hicum_thermal_update(HICUMmodel *inModel, HICUMinstance *inInstance, double
rcx_scaled = model->HICUMrcx / area_times_m;
rbi0_scaled = model->HICUMrbi0 / area_times_m;
//these variables depend only on scale, but not on temperature.
// These model variables depend on on scale.
// They are put into the here struct for usage in load routine.
here->HICUMicbar_scaled = icbar_scaled;
here->HICUMcbepar_scaled = cbepar_scaled;
here->HICUMcbcpar_scaled = cbcpar_scaled;
here->HICUMcth_scaled = cth_scaled;
here->HICUMibets_scaled = ibets_scaled;
here->HICUMrci0_scaled = rci0_scaled;
here->HICUMrth_scaled = rth_scaled;
here->HICUMcjcx0_scaled = cjcx0_scaled;
here->HICUMcjci0_scaled = cjci0_scaled;
Tnom = model->HICUMtnom;
k10 = model->HICUMf1vg*Tnom*log(Tnom);
@ -470,7 +478,7 @@ int hicum_thermal_update(HICUMmodel *inModel, HICUMinstance *inInstance, double
vdep_t.dpart(here->HICUMvdep_t.dpart);
//Tunneling current factors
if (model->HICUMibets > 0) { // HICTUN_T
if (ibets_scaled > 0) { // HICTUN_T
duals::duald a_eg,ab,aa;
ab = 1.0;
aa = 1.0;
@ -482,7 +490,7 @@ int hicum_thermal_update(HICUMmodel *inModel, HICUMinstance *inInstance, double
ab = (cjei0_t/cjei0_scaled)*sqrt(a_eg)*vdei_t*vdei_t/(model->HICUMvdei*model->HICUMvdei);
aa = (model->HICUMvdei/vdei_t)*(cjei0_scaled/cjei0_t)*pow(a_eg,-1.5);
}
a = model->HICUMibets*ab;
a = ibets_scaled*ab;
here->HICUMibets_t.rpart = a.rpart();
here->HICUMibets_t.dpart = a.dpart();
a = model->HICUMabet*aa;