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

started derivatives of hicum2 minority charge

This commit is contained in:
mariok 2020-04-20 09:59:49 +02:00
parent c7e5df27db
commit 9da946bf2d
2 changed files with 58 additions and 35 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
.gitignore vendored
View File

@ -73,3 +73,5 @@ Makefile.in
# Visual Studio user options files # Visual Studio user options files
**/*.vcxproj.user **/*.vcxproj.user
# Visual Studio Code user options files
.vscode/

91
src/spicelib/devices/hicum2/hicum2load.c
View File

@ -502,6 +502,9 @@ HICUMload(GENmodel *inModel, CKTcircuit *ckt)
double a_h,Q_pT,d_Q; double a_h,Q_pT,d_Q;
double Qf,Cdei,Qr,Cdci; double Qf,Cdei,Qr,Cdci;
double ick,vc,cjcx01,cjcx02; double ick,vc,cjcx01,cjcx02;
double ick_Vciei;
double Qf_Vbiei, Qf_Vbici, Q_bf_Vbiei, Q_bf_Vbici, Qr_Vbiei, Qr_Vbici, Q_pT_Vbiei, Q_pT_Vbici;
int l_it; int l_it;
//NQS //NQS
@ -1226,6 +1229,7 @@ HICUMload(GENmodel *inModel, CKTcircuit *ckt)
//Critical current for onset of high-current effects //Critical current for onset of high-current effects
{ // HICICK { // HICICK
double Ovpt,a,d1,vceff,a1,a11,Odelck,ick1,ick2,ICKa; double Ovpt,a,d1,vceff,a1,a11,Odelck,ick1,ick2,ICKa;
double d1_Vciei, vceff_Vciei, a1_Vciei, ick1_Vciei, ick2_Vciei, ICKa_Vciei;
Ovpt = 1.0/model->HICUMvpt; Ovpt = 1.0/model->HICUMvpt;
a = vc/here->HICUMvt; a = vc/here->HICUMvt;
d1 = a-1; d1 = a-1;
@ -1239,6 +1243,18 @@ HICUMload(GENmodel *inModel, CKTcircuit *ckt)
ICKa = (vceff-here->HICUMvlim_t)*Ovpt; ICKa = (vceff-here->HICUMvlim_t)*Ovpt;
ick = ick2*(1.0+0.5*(ICKa+sqrt(ICKa*ICKa+model->HICUMaick))); ick = ick2*(1.0+0.5*(ICKa+sqrt(ICKa*ICKa+model->HICUMaick)));
// derivative: maxima
d1_Vciei = 1/here->HICUMvt;
// vceff_Vciei = (((d_1*d1_Vciei/sqrt(d_1^2+1.921812)+d1_Vciei)*v_t)/2
vceff_Vciei = (d1*d1_Vciei/sqrt(d1*d1+1.921812)+d1_Vciei)*here->HICUMvt/2;
a1_Vciei = vceff_Vciei/here->HICUMvlim_t;
// ick1_Vciei = (exp(log(exp(d_elck*log(a1))+1)/d_elck+d_elck*log(a1))*a1_Vciei)/(a1*(exp(d_elck*log(a1))+1))
ick1_Vciei = (exp(Odelck*log(1+exp(model->HICUMdelck*log(a1))) + model->HICUMdelck*log(a1))*a1_Vciei)/(a1*(exp(model->HICUMdelck*log(a1))+1));
ick2_Vciei = -1.0*a11*ick1_Vciei/ick1/ick1;
ICKa_Vciei = vceff_Vciei*Ovpt;
// ick_Vciei = (0.5*(sqrt(I_CKa(x)^2+a_ick)+I_CKa(x))+1.0)*('diff(i_ck2(x),x,1))+0.5*i_ck2(x)*((I_CKa(x)*('diff(I_CKa(x),x,1)))/sqrt(I_CKa(x)^2+a_ick)+'diff(I_CKa(x),x,1))
ick_Vciei = (0.5*(sqrt(ICKa*ICKa+model->HICUMaick)+ICKa)+1.0)*ick2_Vciei+0.5*ick2*((ICKa*(ICKa_Vciei))/sqrt(ICKa*ICKa+model->HICUMaick)+ICKa_Vciei);
} }
//Initialization //Initialization
@ -1283,7 +1299,6 @@ HICUMload(GENmodel *inModel, CKTcircuit *ckt)
Tf = T_f0; Tf = T_f0;
Qf = T_f0*itf; Qf = T_f0*itf;
HICQFF(here, model, itf,ick,&Tf,&Qf,&T_fT,&Q_fT,&Q_bf); HICQFF(here, model, itf,ick,&Tf,&Qf,&T_fT,&Q_fT,&Q_bf);
//TODO: itf=f(Vbiei,Vbici) -> Qf, Q_bf Ableitungen nach Vbiei, Vbici
//Initial formulation of reverse diffusion charge //Initial formulation of reverse diffusion charge
Qr = Tr*itr; Qr = Tr*itr;
@ -1293,7 +1308,6 @@ HICUMload(GENmodel *inModel, CKTcircuit *ckt)
//Iteration for Q_pT is required for improved initial solution //Iteration for Q_pT is required for improved initial solution
Qf = sqrt(T_f0*itf*Q_fT); Qf = sqrt(T_f0*itf*Q_fT);
Q_pT = Q_0+Qf+Qr; Q_pT = Q_0+Qf+Qr;
//TODO: Q_pT_Vbiei, Vbici
d_Q = Q_pT; d_Q = Q_pT;
while (fabs(d_Q) >= RTOLC*fabs(Q_pT) && l_it <= l_itmax) { while (fabs(d_Q) >= RTOLC*fabs(Q_pT) && l_it <= l_itmax) {
double a; double a;
@ -1337,6 +1351,21 @@ HICUMload(GENmodel *inModel, CKTcircuit *ckt)
} //if } //if
// TODO calculate here the derivatives once! They are not needed inside the newton... // TODO calculate here the derivatives once! They are not needed inside the newton...
// Qf = Tf*itf;
Qf_Vbiei = 0;
Qf_Vbici = 0;
// Qr = Tr*itr;
Qr_Vbiei = 0;
Qr_Vbici = 0;
// Q_pT = Q_0+Qf+Qr;
Q_pT_Vbiei = Q_0_Vbiei + Qf_Vbiei + Qr_Vbiei;
Q_pT_Vbici = Q_0_Vbici + Qf_Vbici + Qr_Vbici;
// Q_bf = Q_0+Qf+Qr;
Q_bf_Vbiei = 0;
Q_bf_Vbici = 0;
// itf, itr
//TODO: itf=f(Vbiei,Vbici) -> Qf, Q_bf Ableitungen nach Vbiei, Vbici
itf_Vbiei = itf/VT_f; // TODO: missing the derivatives of Qf itf_Vbiei = itf/VT_f; // TODO: missing the derivatives of Qf
itr_Vbici = itr/here->HICUMvt; // TODO: missing the derivatives of Qpt itr_Vbici = itr/here->HICUMvt; // TODO: missing the derivatives of Qpt
@ -1395,10 +1424,10 @@ HICUMload(GENmodel *inModel, CKTcircuit *ckt)
sq_smooth = sqrt(denom*denom+0.01); sq_smooth = sqrt(denom*denom+0.01);
hl = 0.5*(denom+sq_smooth); hl = 0.5*(denom+sq_smooth);
iavl = itf*avl/hl; iavl = itf*avl/hl;
iavl_Vbici = itf*avl_Vbici/hl; iavl_Vbici = itf*avl_Vbici/hl; // TODO itf_Vbici
} else { } else {
iavl = itf*avl; iavl = itf*avl;
iavl_Vbici = itf*avl_Vbici; iavl_Vbici = itf*avl_Vbici; // TODO itf_Vbici
} }
} }
} else { } else {
@ -1447,6 +1476,7 @@ HICUMload(GENmodel *inModel, CKTcircuit *ckt)
// Consideration of peripheral charge // Consideration of peripheral charge
if(Qf > 0.0) { if(Qf > 0.0) {
rbi = rbi*(Qjei+Qf*model->HICUMfqi)/(Qjei+Qf); rbi = rbi*(Qjei+Qf*model->HICUMfqi)/(Qjei+Qf);
// Qf_Viei and Qf_Vbici
rbi_Vbiei = (Qjei+Qf*model->HICUMfqi)*rbi_Vbiei/(Qjei+Qf) + rbi*Cjei/(Qjei+Qf) - (Qjei+Qf*model->HICUMfqi)*rbi*Cjei/(Qjei+Qf)/(Qjei+Qf); rbi_Vbiei = (Qjei+Qf*model->HICUMfqi)*rbi_Vbiei/(Qjei+Qf) + rbi*Cjei/(Qjei+Qf) - (Qjei+Qf*model->HICUMfqi)*rbi*Cjei/(Qjei+Qf)/(Qjei+Qf);
rbi_Vbici = rbi_Vbici*(Qjei+Qf*model->HICUMfqi)/(Qjei+Qf); rbi_Vbici = rbi_Vbici*(Qjei+Qf*model->HICUMfqi)/(Qjei+Qf);
} }
@ -1574,22 +1604,16 @@ HICUMload(GENmodel *inModel, CKTcircuit *ckt)
// Excess Phase calculation // Excess Phase calculation
if ((model->HICUMflnqs != 0 || model->HICUMflcomp == 0.0 || model->HICUMflcomp == 2.1) && Tf != 0 && (model->HICUMalit > 0 || model->HICUMalqf > 0)) { if ((model->HICUMflnqs != 0 || model->HICUMflcomp == 0.0 || model->HICUMflcomp == 2.1) && Tf != 0 && (model->HICUMalit > 0 || model->HICUMalqf > 0)) {
// Vxf1 = Vbxf1; Vxf1 = Vbxf1;
// Vxf2 = Vbxf2; Vxf2 = Vbxf2;
// Ixf1 = (Vxf2-itf)/Tf*model->HICUMt0; Ixf1 = (Vxf2-itf)/Tf*model->HICUMt0;
// Ixf2 = (Vxf2-Vxf1)/Tf*model->HICUMt0; Ixf2 = (Vxf2-Vxf1)/Tf*model->HICUMt0;
// Qxf1 = model->HICUMalit*model->HICUMt0*Vxf1; Qxf1 = model->HICUMalit*model->HICUMt0*Vxf1;
// Qxf1_Vxf1 = model->HICUMalit*model->HICUMt0; Qxf1_Vxf1 = model->HICUMalit*model->HICUMt0;
// Qxf2 = model->HICUMalit*model->HICUMt0*Vxf2/3; Qxf2 = model->HICUMalit*model->HICUMt0*Vxf2/3;
// Qxf2_Vxf2 = model->HICUMalit*model->HICUMt0/3; Qxf2_Vxf2 = model->HICUMalit*model->HICUMt0/3;
// Itxf = Vxf2; Itxf = Vxf2;
Ixf1 = Vbxf1;
Ixf2 = Vbxf2;
Qxf1 = 0;
Qxf2 = 0;
Qxf1_Vxf1 = 0;
Qxf2_Vxf2 = 0;
// TODO derivatives of Ixf1 and Ixf2 // TODO derivatives of Ixf1 and Ixf2
@ -2316,7 +2340,7 @@ c Branch: sis, Stamp element: Rsu
/* /*
c Branch: xf1-ground, Stamp element: Ixf1 c Branch: xf1-ground, Stamp element: Ixf1
*/ */
// rhs_current = (Ixf1 - Ixf1_Vrth*Vrth - Ixf1_Vbiei*Vbiei - Ixf1_Vbici*Vbici - Ixf2_Vxf2*Vxf2); // rhs_current = (Ixf1 - Ixf1_Vxf1*Vxf1 - Ixf1_Vrth*Vrth - Ixf1_Vbiei*Vbiei - Ixf1_Vbici*Vbici - Ixf1_Vxf2*Vxf2); // TODO
rhs_current = Ixf1; rhs_current = Ixf1;
*(ckt->CKTrhs + here->HICUMxf1Node) += rhs_current; // into xf1 node *(ckt->CKTrhs + here->HICUMxf1Node) += rhs_current; // into xf1 node
// *(here->HICUMxf1TempPtr) += -Ixf1_Vrth; // *(here->HICUMxf1TempPtr) += -Ixf1_Vrth;
@ -2324,33 +2348,30 @@ c Branch: xf1-ground, Stamp element: Ixf1
// *(here->HICUMxf1EmitEIPtr) += +Ixf1_Vbiei; // *(here->HICUMxf1EmitEIPtr) += +Ixf1_Vbiei;
// *(here->HICUMxf1BaseBIPtr) += -Ixf1_Vbici; // *(here->HICUMxf1BaseBIPtr) += -Ixf1_Vbici;
// *(here->HICUMxf1CollCIPtr) += +Ixf1_Vbici; // *(here->HICUMxf1CollCIPtr) += +Ixf1_Vbici;
// *(here->HICUMxf1Xf2Ptr) += +Ixf1_Vxf2; // TODO // *(here->HICUMxf1Xf2Ptr) += +Ixf1_Vxf2;
/* /*
c Branch: xf1-ground, Stamp element: Qxf1 c Branch: xf1-ground, Stamp element: Qxf1 // TODO Test in AC simulation!
*/ */
rhs_current = Iqxf1 - Iqxf1_Vxf1*Vxf1; // rhs_current = Iqxf1 - Iqxf1_Vxf1*Vxf1;
*(ckt->CKTrhs + here->HICUMxf1Node) += rhs_current; // into ground // *(ckt->CKTrhs + here->HICUMxf1Node) += rhs_current; // into ground
*(here->HICUMxf1Xf1Ptr) += Iqxf1_Vxf1; // *(here->HICUMxf1Xf1Ptr) += Iqxf1_Vxf1;
/*
c Branch: xf1-ground, Stamp element: Rxf1
*/
*(here->HICUMxf1Xf1Ptr) += 1; // current Ixf1 is normalized to Tf
/* /*
c Branch: xf2-ground, Stamp element: Ixf2 c Branch: xf2-ground, Stamp element: Ixf2
*/ */
// rhs_current = (Ixf2 - Ixf2_Vxf1*Vxf1 - Ixf2_Vrth*Vrth - Ixf2_Vbiei*Vbiei - Ixf2_Vbici*Vbici - Ixf2_Vxf2*Vxf2); // TODO
rhs_current = Ixf2; rhs_current = Ixf2;
*(ckt->CKTrhs + here->HICUMxf2Node) += rhs_current; // into xf node *(ckt->CKTrhs + here->HICUMxf2Node) += rhs_current; // into xf2 node
// *(here->HICUMxf2TempPtr) += -Ixf2_Vrth; // *(here->HICUMxf2TempPtr) += -Ixf2_Vrth;
// *(here->HICUMxf2BaseBIPtr) += -Ixf2_Vbiei; // *(here->HICUMxf2BaseBIPtr) += -Ixf2_Vbiei;
// *(here->HICUMxf2EmitEIPtr) += +Ixf2_Vbiei; // *(here->HICUMxf2EmitEIPtr) += +Ixf2_Vbiei;
// *(here->HICUMxf2BaseBIPtr) += -Ixf2_Vbici; // *(here->HICUMxf2BaseBIPtr) += -Ixf2_Vbici;
// *(here->HICUMxf2CollCIPtr) += +Ixf2_Vbici; // *(here->HICUMxf2CollCIPtr) += +Ixf2_Vbici;
/* /*
c Branch: xf2-ground, Stamp element: Qxf2 c Branch: xf2-ground, Stamp element: Qxf2 // TODO Test in AC simulation!
*/ */
rhs_current = Iqxf2 - Iqxf2_Vxf2*Vxf2; // rhs_current = Iqxf2 - Iqxf2_Vxf2*Vxf2;
*(ckt->CKTrhs + here->HICUMxf2Node) += rhs_current; // into ground // *(ckt->CKTrhs + here->HICUMxf2Node) += rhs_current; // into ground
*(here->HICUMxf2Xf2Ptr) += Iqxf2_Vxf2; // *(here->HICUMxf2Xf2Ptr) += Iqxf2_Vxf2;
/* /*
c Branch: xf2-ground, Stamp element: Rxf2 c Branch: xf2-ground, Stamp element: Rxf2
*/ */
@ -2367,7 +2388,7 @@ c Branch: xf-ground, Stamp element: Ixf
// *(here->HICUMxfBaseBIPtr) += -Ixf_Vbici; // *(here->HICUMxfBaseBIPtr) += -Ixf_Vbici;
// *(here->HICUMxfCollCIPtr) += +Ixf_Vbici; // *(here->HICUMxfCollCIPtr) += +Ixf_Vbici;
/* /*
c Branch: xf-ground, Stamp element: Qxf c Branch: xf-ground, Stamp element: Qxf // TODO Test in AC simulation!
*/ */
// rhs_current = model->HICUMtype * (Iqxf - Iqxf_Vxf*Vxf); // rhs_current = model->HICUMtype * (Iqxf - Iqxf_Vxf*Vxf);
// *(ckt->CKTrhs + here->HICUMxfNode) += rhs_current; // into ground // *(ckt->CKTrhs + here->HICUMxfNode) += rhs_current; // into ground