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

finally: avalanche at very high T fix

This commit is contained in:
Markus Mueller 2020-07-23 16:22:46 +02:00
parent 64704fd53e
commit 03c3efd762
4 changed files with 53 additions and 48 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -75,6 +75,7 @@ typedef struct sHICUMinstance {
double HICUMicVCS; /* initial condition inner C-S branch */
double HICUMtemp; /* instance temperature */
double HICUMtemp_Vrth;/* derivative device temperature to temperature at thermal node */
double HICUMdtemp; /* instance delta temperature */
double HICUMdtemp_sh; /* instance delta temperature because of self-heating */
double HICUMm; /* multiply factor for the hicum */

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

@ -484,6 +484,7 @@ HICUMload(GENmodel *inModel, CKTcircuit *ckt)
double Q_0_Vbiei, Q_0_Vbici, Q_0_dT;
double Temp;
double Tdev_Vrth; //derivative device temperature to Vrth
double Cjei_Vbiei,Cjci_Vbici,Cjep_Vbpei,Cjep_dT,Cjs_Vsici;
double Cjei_dT, Cjci_dT;
@ -1877,7 +1878,7 @@ HICUMload(GENmodel *inModel, CKTcircuit *ckt)
if (selfheat) { // Thermal_update_with_self_heating
Temp = here->HICUMtemp+Vrth;
_iret = hicum_thermal_update(model, here, Temp);
_iret = hicum_thermal_update(model, here, &Temp, &Tdev_Vrth);
here->HICUMdtemp_sh = Temp - here->HICUMtemp;
} else {
@ -2502,7 +2503,7 @@ HICUMload(GENmodel *inModel, CKTcircuit *ckt)
if(!selfheat) {
Ith = 0;
} else if (model->HICUMrth_de==1) {
Ith = -dT/here->HICUMrth_t.rpart+pterm; //Current from gnd to T
Ith = -Vrth/here->HICUMrth_t.rpart+pterm; //Current from gnd to T
} else {
Ith = -Vrth/here->HICUMrth_t.rpart+pterm; //Current from gnd to T
@ -3324,125 +3325,125 @@ load:
// #############################################################
// Stamp element: Ibiei f_Bi = + f_Ei = -
rhs_current = -Ibiei_Vrth*Vrth;
rhs_current = -Ibiei_Vrth*Tdev_Vrth*Vrth;
*(ckt->CKTrhs + here->HICUMbaseBINode) += -rhs_current;
*(ckt->CKTrhs + here->HICUMemitEINode) += rhs_current;
// with respect to Potential Vrth
*(here->HICUMbaseBItempPtr) += Ibiei_Vrth;
*(here->HICUMemitEItempPtr) += -Ibiei_Vrth;
*(here->HICUMbaseBItempPtr) += Ibiei_Vrth*Tdev_Vrth;
*(here->HICUMemitEItempPtr) += -Ibiei_Vrth*Tdev_Vrth;
// finish
// Stamp element: Ibpei f_Bp = + f_Ei = -
rhs_current = -Ibpei_Vrth*Vrth;
rhs_current = -Ibpei_Vrth*Tdev_Vrth*Vrth;
*(ckt->CKTrhs + here->HICUMbaseBPNode) += -rhs_current;
*(ckt->CKTrhs + here->HICUMemitEINode) += rhs_current;
// with respect to Potential Vrth
*(here->HICUMbaseBPtempPtr) += Ibpei_Vrth;
*(here->HICUMemitEItempPtr) += -Ibpei_Vrth;
*(here->HICUMbaseBPtempPtr) += Ibpei_Vrth*Tdev_Vrth;
*(here->HICUMemitEItempPtr) += -Ibpei_Vrth*Tdev_Vrth;
// finish
// Stamp element: Ibici f_Bi = + f_Ci = -
rhs_current = -Ibici_Vrth*Vrth;
rhs_current = -Ibici_Vrth*Vrth*Tdev_Vrth;
*(ckt->CKTrhs + here->HICUMbaseBINode) += -rhs_current;
*(ckt->CKTrhs + here->HICUMcollCINode) += rhs_current;
// with respect to Potential Vrth
*(here->HICUMbaseBItempPtr) += Ibici_Vrth;
*(here->HICUMcollCItempPtr) += -Ibici_Vrth;
*(here->HICUMbaseBItempPtr) += Ibici_Vrth*Tdev_Vrth;
*(here->HICUMcollCItempPtr) += -Ibici_Vrth*Tdev_Vrth;
// finish
// Stamp element: Iciei f_Ci = + f_Ei = -
rhs_current = -Iciei_Vrth*Vrth;
rhs_current = -Iciei_Vrth*Vrth*Tdev_Vrth;
*(ckt->CKTrhs + here->HICUMcollCINode) += -rhs_current;
*(ckt->CKTrhs + here->HICUMemitEINode) += rhs_current;
// with respect to Potential Vrth
*(here->HICUMcollCItempPtr) += Iciei_Vrth;
*(here->HICUMemitEItempPtr) += -Iciei_Vrth;
*(here->HICUMcollCItempPtr) += Iciei_Vrth*Tdev_Vrth;
*(here->HICUMemitEItempPtr) += -Iciei_Vrth*Tdev_Vrth;
// finish
// Stamp element: Ibpci f_Bp = + f_Ci = -
rhs_current = -Ibpci_Vrth*Vrth;
rhs_current = -Ibpci_Vrth*Vrth*Tdev_Vrth;
*(ckt->CKTrhs + here->HICUMbaseBPNode) += -rhs_current;
*(ckt->CKTrhs + here->HICUMcollCINode) += rhs_current;
// with respect to Potential Vrth
*(here->HICUMbaseBPtempPtr) += Ibpci_Vrth;
*(here->HICUMcollCItempPtr) += -Ibpci_Vrth;
*(here->HICUMbaseBPtempPtr) += Ibpci_Vrth*Tdev_Vrth;
*(here->HICUMcollCItempPtr) += -Ibpci_Vrth*Tdev_Vrth;
// finish
// Stamp element: Rcx f_Ci = + f_C = -
rhs_current = -Icic_Vrth*Vrth;
rhs_current = -Icic_Vrth*Vrth*Tdev_Vrth;
*(ckt->CKTrhs + here->HICUMcollCINode) += -rhs_current;
*(ckt->CKTrhs + here->HICUMcollNode) += rhs_current;
// with respect to Potential Vrth
*(here->HICUMcollCItempPtr) += Icic_Vrth;
*(here->HICUMcollTempPtr) += -Icic_Vrth;
*(here->HICUMcollCItempPtr) += Icic_Vrth*Tdev_Vrth;
*(here->HICUMcollTempPtr) += -Icic_Vrth*Tdev_Vrth;
// finish
// Stamp element: Rbx f_B = + f_Bp = -
rhs_current = -Ibbp_Vrth*Vrth;
rhs_current = -Ibbp_Vrth*Vrth*Tdev_Vrth;
*(ckt->CKTrhs + here->HICUMbaseNode) += -rhs_current;
*(ckt->CKTrhs + here->HICUMbaseBPNode) += rhs_current;
// with respect to Potential Vrth
*(here->HICUMbaseTempPtr) += Ibbp_Vrth;
*(here->HICUMbaseBPtempPtr) += -Ibbp_Vrth;
*(here->HICUMbaseTempPtr) += Ibbp_Vrth*Tdev_Vrth;
*(here->HICUMbaseBPtempPtr) += -Ibbp_Vrth*Tdev_Vrth;
// finish
// Stamp element: Re f_Ei = + f_E = -
rhs_current = - Ieie_Vrth*Vrth;
rhs_current = - Ieie_Vrth*Vrth*Tdev_Vrth;
*(ckt->CKTrhs + here->HICUMemitEINode) += -rhs_current;
*(ckt->CKTrhs + here->HICUMemitNode) += +rhs_current;
// with respect to Potential Vrth
*(here->HICUMemitEItempPtr) += +Ieie_Vrth;
*(here->HICUMemitTempPtr) += -Ieie_Vrth;
*(here->HICUMemitEItempPtr) += +Ieie_Vrth*Tdev_Vrth;
*(here->HICUMemitTempPtr) += -Ieie_Vrth*Tdev_Vrth;
// finish
// Stamp element: Rbi f_Bp = + f_Bi = -
rhs_current = -Ibpbi_Vrth*Vrth;
rhs_current = -Ibpbi_Vrth*Vrth*Tdev_Vrth;
*(ckt->CKTrhs + here->HICUMbaseBPNode) += -rhs_current;
*(ckt->CKTrhs + here->HICUMbaseBINode) += rhs_current;
// with respect to Potential Vrth
*(here->HICUMbaseBPtempPtr) += Ibpbi_Vrth;
*(here->HICUMbaseBItempPtr) += -Ibpbi_Vrth;
*(here->HICUMbaseBPtempPtr) += Ibpbi_Vrth*Tdev_Vrth;
*(here->HICUMbaseBItempPtr) += -Ibpbi_Vrth*Tdev_Vrth;
// finish
// Stamp element: Isici f_Si = + f_Ci = -
rhs_current = -Isici_Vrth*Vrth;
rhs_current = -Isici_Vrth*Vrth*Tdev_Vrth;
*(ckt->CKTrhs + here->HICUMsubsSINode) += -rhs_current;
*(ckt->CKTrhs + here->HICUMcollCINode) += rhs_current;
// with respect to Potential Vrth
*(here->HICUMsubsSItempPtr) += Isici_Vrth;
*(here->HICUMcollCItempPtr) += -Isici_Vrth;
*(here->HICUMsubsSItempPtr) += Isici_Vrth*Tdev_Vrth;
*(here->HICUMcollCItempPtr) += -Isici_Vrth*Tdev_Vrth;
// finish
// Branch: bpsi, Stamp element: Its
rhs_current = - Ibpsi_Vrth*Vrth;
rhs_current = - Ibpsi_Vrth*Vrth*Tdev_Vrth;
*(ckt->CKTrhs + here->HICUMbaseBPNode) += -rhs_current;
*(ckt->CKTrhs + here->HICUMsubsSINode) += rhs_current;
// f_Bp = + f_Si = -
// with respect to Vrth
*(here->HICUMbaseBPtempPtr) += Ibpsi_Vrth;
*(here->HICUMsubsSItempPtr) += -Ibpsi_Vrth;
*(here->HICUMbaseBPtempPtr) += Ibpsi_Vrth*Tdev_Vrth;
*(here->HICUMsubsSItempPtr) += -Ibpsi_Vrth*Tdev_Vrth;
// finish
if (nqs) {
// Stamp element: Ixf f_xf = +
rhs_current = -Ixf_dT*Vrth;
rhs_current = -Ixf_dT*Vrth*Tdev_Vrth;
*(ckt->CKTrhs + here->HICUMxfNode) += -rhs_current;
// with respect to Potential Vxf
*(here->HICUMxfTempPtr) += Ixf_dT;
*(here->HICUMxfTempPtr) += Ixf_dT*Tdev_Vrth;
// finish
// Stamp element: Ixf1 f_xf1 = +
rhs_current = -Ixf1_dT*Vrth;
*(ckt->CKTrhs + here->HICUMxf1Node) += -rhs_current;
// with respect to Potential Vxf1
*(here->HICUMxf1TempPtr) += Ixf1_dT;
*(here->HICUMxf1TempPtr) += Ixf1_dT*Tdev_Vrth;
// finish
// Stamp element: Ixf2 f_xf2 = +
rhs_current = -Ixf2_dT*Vrth;
rhs_current = -Ixf2_dT*Vrth*Tdev_Vrth;
*(ckt->CKTrhs + here->HICUMxf2Node) += -rhs_current;
// with respect to Potential Vxf2
*(here->HICUMxf2TempPtr) += Ixf2_dT;
*(here->HICUMxf2TempPtr) += Ixf2_dT*Tdev_Vrth;
// finish
}

2
src/spicelib/devices/hicum2/hicumL2.hpp
View File

@ -8,7 +8,7 @@ extern "C" {
void hicum_diode(double T, double IS, double UM1, double U, double *Iz, double *Gz, double *Tz);
void hicum_qjmodf(double T, double c_0, double u_d, double z, double a_j, double U_cap, double *C, double *C_dU, double *C_dvt, double *Qz, double *Qz_dU, double *Qz_dvt);
static double HICUMlimitlog( double deltemp, double deltemp_old, double LIM_TOL, int *check);
int hicum_thermal_update(HICUMmodel *, HICUMinstance *, double Temp);
int hicum_thermal_update(HICUMmodel *, HICUMinstance *, double * Temp, double * Tdev_Vrth);
int HICUMload(GENmodel *inModel, CKTcircuit *ckt);
int HICUMtemp(GENmodel *inModel, CKTcircuit *ckt);
#ifdef __cplusplus

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

@ -120,14 +120,14 @@ HICUMtemp(GENmodel *inModel, CKTcircuit *ckt)
if(here->HICUMdtempGiven) here->HICUMtemp = here->HICUMtemp + here->HICUMdtemp;
iret = hicum_thermal_update(model, here, here -> HICUMtemp);
iret = hicum_thermal_update(model, here, &here -> HICUMtemp, &here->HICUMtemp_Vrth);
}
}
return(OK);
}
int hicum_thermal_update(HICUMmodel *inModel, HICUMinstance *inInstance, double HICUMTemp)
int hicum_thermal_update(HICUMmodel *inModel, HICUMinstance *inInstance, double * HICUMTemp, double * Tdev_Vrth)
{
HICUMmodel *model = (HICUMmodel *)inModel;
HICUMinstance *here = (HICUMinstance *)inInstance;
@ -158,14 +158,17 @@ int hicum_thermal_update(HICUMmodel *inModel, HICUMinstance *inInstance, double
zetasct = mg-1.5;
// Limit temperature to avoid FPEs in equations
if(HICUMTemp < TMIN + CONSTCtoK) {
HICUMTemp = TMIN + CONSTCtoK;
*(Tdev_Vrth) = 1;
if(*(HICUMTemp) < TMIN + CONSTCtoK) {
*(HICUMTemp) = TMIN + CONSTCtoK;
*(Tdev_Vrth) = 0;
} else {
if (HICUMTemp > TMAX + CONSTCtoK) {
HICUMTemp = TMAX + CONSTCtoK;
if (*(HICUMTemp) > TMAX + CONSTCtoK) {
*(HICUMTemp) = TMAX + CONSTCtoK;
*(Tdev_Vrth) = 0;
}
}
temp = HICUMTemp+1_e; // dual number valued temperature
temp = *(HICUMTemp)+1_e* *(Tdev_Vrth); // dual number device temperature
vt = temp*CONSTKoverQ; // dual valued temperature voltage
here->HICUMvt0 = Tnom * CONSTKoverQ;