plasma.h

来自「等离子体覆盖目标的电磁特性的FDTD分析」· C头文件 代码 · 共 34 行

#define CollisionalFreq 150.0e9*2.0*pii//100e6  

double PlasmaFreq(int z)


{

return   sqrt(5.0e17*1.6e-19*1.6e-19/9.11e-31/8.854e-12*(z-MZ/2+thick)/(2.0*thick));
}

double chi(int m,int z)
{
	double wp,vc;
	wp=PlasmaFreq(z);
	vc=CollisionalFreq;

	return (epsn*wp*wp/vc)*(1.0-exp(-1.0*vc*dt))*exp(-1.0*m*vc*dt);

}

double xi(int m,int z)
{
	double wp,vc;
	wp=PlasmaFreq(z);
	vc=CollisionalFreq;
	return (epsn*wp*wp/(vc*vc*dt))*(1.0-(1.0+vc*dt)*exp(-1.0*vc*dt))*exp(-m*vc*dt);

	
}

inline double dChi0(int z){	return chi(0,z)-chi(1,z);}
inline double dXi0(int z){	return xi(0,z)-xi(1,z);}