!Caculate Generation G In Gaas
Do K=Z_Gaas1,Z_Gaas2
Do J=Y_Gaas1,Y_Gaas2
Do I=X_Gaas1,X_Gaas2
If((((I-I0)*Dx/G_Dx*2)**2+((J-J0)*Dy/G_Dy*2)**2)<=1)Then !*** G Is Only Exist In The Pulse Circle
G(I,J,K)=G_I*G_A/G_H/G_V*Exp(-G_A*(K0-K)*Dz)*Exp(-((I-I0)*Dx)**2/G_Sigma_X**2-((J-J0)*Dy)**2/G_Sigma_Y**2)*Exp(-((T-T0)*Dt-(K0-K)*Dz*Sqrt(E*U0))**2/G_Sigma_T**2)
End If
End Do
End Do
End Do
Do J=Y_Gaas1,Y_Gaas2-1
Do I=X_Gaas1,X_Gaas2-1
If(Flag(I,J)/=0)Then
G(I,J,Z_Gaas2)=0.0
G(I+1,J,Z_Gaas2)=0.0
G(I,J+1,Z_Gaas2)=0.0
G(I+1,J+1,Z_Gaas2)=0.0
End If
End Do
End Do
!Caculate Recombination R In Gaas
Do K=Z_Gaas1,Z_Gaas2
Do J=Y_Gaas1,Y_Gaas2
Do I=X_Gaas1,X_Gaas2
If(N(I,J,K)*P(I,J,K)/=0) Then
R(I,J,K)=(N(I,J,K)*P(I,J,K))/(Tn0*N(I,J,K)+Tp0*P(I,J,K))
End If
End Do
End Do
End Do
!Carrier N In Gaas
Do K=Z_Gaas1,Z_Gaas2
Do J=Y_Gaas1,Y_Gaas2
Do I=X_Gaas1,X_Gaas2
N(I,J,K)=N(I,J,K)+Dt/Q*((Jnx(I,J,K)-Jnx(I-1,J,K))/Dx+(Jny(I,J,K)-Jny(I,J-1,K))/Dy+(Jnz(I,J,K)-Jnz(I,J,K-1))/Dz+Q*G(I,J,K)-Q*R(I,J,K))
End Do
End Do
End Do
!Carrier P In Gaas
Do K=Z_Gaas1,Z_Gaas2
Do J=Y_Gaas1,Y_Gaas2
Do I=X_Gaas1,X_Gaas2
P(I,J,K)=P(I,J,K)-Dt/Q*((Jpx(I,J,K)-Jpx(I-1,J,K))/Dx+(Jpy(I,J,K)-Jpy(I,J-1,K))/Dy+(Jpz(I,J,K)-Jpz(I,J,K-1))/Dz-Q*G(I,J,K)+Q*R(I,J,K))
End Do
End Do
End Do
!Caculate Jnx In Gaas
Do K=Z_Gaas1,Z_Gaas2
Do J=Y_Gaas1,Y_Gaas2
Do I=X_Gaas1,X_Gaas2-1
E_X=Edcx(I,J,K)+Ex(I,J,K)
Unx=Un0/Sqrt(1+(Un0*E_X/Vnsat)**2)
Dnx=Kb*Temp/Q*Unx
Jnx(I,J,K)=Q*Unx*(N(I+1,J,K)+N(I,J,K))/2*E_X +Q*Dnx*(N(I+1,J,K)-N(I,J,K))/Dx
End Do
End Do
End Do
!Caculate Jny In Gaas
Do K=Z_Gaas1,Z_Gaas2
Do J=Y_Gaas1,Y_Gaas2-1
Do I=X_Gaas1,X_Gaas2
E_Y=Edcy(I,J,K)+Ey(I,J,K)
Uny=Un0 /Sqrt(1+(Un0*E_Y/Vnsat)**2)
Dny=Kb*Temp/Q*Uny
Jny(I,J,K)=Q*Uny*(N(I,J+1,K)+N(I,J,K))/2.0*E_Y +Q*Dny*(N(I,J+1,K)-N(I,J,K))/Dy
End Do
End Do
End Do
!Caculate Jnz In Gaas
Do K=Z_Gaas1,Z_Gaas2-1
Do J=Y_Gaas1,Y_Gaas2
Do I=X_Gaas1,X_Gaas2
E_Z=Edcz(I,J,K)+Ez(I,J,K)
Unz=Un0/Sqrt(1+(Un0*E_Z/Vnsat)**2)
Dnz=Kb*Temp/Q*Unz
Jnz(I,J,K)=Q*Unz*(N(I,J,K+1)+N(I,J,K))/2*E_Z +Q*Dnz*(N(I,J,K+1)-N(I,J,K))/Dz
End Do
End Do
End Do
!Caculate Jpx In Gaas
Do K=Z_Gaas1,Z_Gaas2
Do J=Y_Gaas1,Y_Gaas2
Do I=X_Gaas1,X_Gaas2-1
E_X=Edcx(I,J,K)+Ex(I,J,K)
Upx=Up0/Sqrt(1+(Up0*E_X/Vpsat)**2)
Dpx=Kb*Temp/Q*Upx
Jpx(I,J,K)=Q*Upx*(P(I+1,J,K)+P(I,J,K))/2*E_X -Q*Dpx*(P(I+1,J,K)-P(I,J,K))/Dx
End Do
End Do
End Do
!Caculate Jpy In Gaas
Do K=Z_Gaas1,Z_Gaas2
Do J=Y_Gaas1,Y_Gaas2-1
Do I=X_Gaas1,X_Gaas2
E_Y=Edcy(I,J,K)+Ey(I,J,K)
Upy=Up0/Sqrt(1+(Up0*E_Y/Vpsat)**2)
Dpy=Kb*Temp/Q*Upy
Jpy(I,J,K)=Q*Upy*(P(I,J+1,K)+P(I,J,K))/2*E_Y -Q*Dpy*(P(I,J+1,K)-P(I,J,K))/Dy
End Do
End Do
End Do
!Caculate Jpz In Gaas
Do K=Z_Gaas1,Z_Gaas2-1
Do J=Y_Gaas1,Y_Gaas2
Do I=X_Gaas1,X_Gaas2
E_Z=Edcz(I,J,K)+Ez(I,J,K)
Upz=Up0/Sqrt(1+(Up0*E_Z/Vpsat)**2)
Dpz=Kb*Temp/Q*Upz
Jpz(I,J,K)=Q*Upz*(P(I,J,K+1)+P(I,J,K))/2*E_Z -Q*Dpz*(P(I,J,K+1)-P(I,J,K))/Dz
End Do
End Do
End Do
本文档为【Caculate Generation G In Gaas】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑,
图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
该文档来自用户分享,如有侵权行为请发邮件ishare@vip.sina.com联系网站客服,我们会及时删除。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。