#include "define.h" #include <stdio.h> #include <math.h> void diag(void) {double vD; //氘離子熱速度double lamd_e, lamd_D;double w_pe, peri_e;double wpe_dt, vdt_dz;//離子溫度Ti=20,離子電荷qi=1.6021865314e-19;氘質量m_D=3.348e-27;vD = sqrt(2.0 * Ti * qi / m_D); //氘離子熱速度,離子溫度Ti=20,離子電荷qi=1.6021865314e-19;//printf("%lf \n", vD); stime();//電子質量me=9.1093897e-31; ve=sqrt(-2.0*qe*Te/me);sqrt(2.0*qe*Te/me); 電子電荷qe=-1.6021865314e-19; 磁場B0=5.3;lamd_e = me * ve / fabs(qe * B0);//電子回旋半徑lamd_D = m_D * vD / fabs(qi * B0);//氘離子回旋半徑//Ne=1e19;電子數密度 //真空介電常數eps0=8.854187817e-12;w_pe = sqrt(Ne * qe * qe / (eps0 * me));//電子震蕩頻率//printf("電子回旋半徑%lf \n", lamd_e); //printf("氘離子回旋半徑%lf \n", lamd_D); //printf("電子震蕩頻率%lf \n", w_pe); stime();//dt=1.0e-12;時間步長 //dz_plas=lz_plas/double(Nz_plas-1);z方向空間步長peri_e = 2.0 * pi / w_pe; //電子震蕩周期周期T=2π/ωwpe_dt = w_pe * dt; vdt_dz = ve * dt / dz_plas;printf("ne=%10.6e,Te=%10.6f\n", Ne, Te);printf("ve=%10.6e,v_D=%10.6e\n", ve, vD);printf("德拜長度=%10.6e,空間步長=%10.6e\n", deby_leng, dz_plas);printf("電子振蕩周期=%10.6e,時間步長=%10.6e\n", peri_e, dt);printf("wpe*dt=%10.6f,v*dt/dz=%10.6f\n", wpe_dt, vdt_dz);}

氘離子熱速度、電子回旋半徑、氘離子回旋半徑、電子震蕩頻率

?

?

?

總結

以上是生活随笔為你收集整理的[C++调试笔记]diag.cpp的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。