#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;氘質(zhì)量m_D=3.348e-27;vD = sqrt(2.0 * Ti * qi / m_D); //氘離子熱速度,離子溫度Ti=20,離子電荷qi=1.6021865314e-19;//printf("%lf \n", vD); stime();//電子質(zhì)量me=9.1093897e-31; ve=sqrt(-2.0*qe*Te/me);sqrt(2.0*qe*Te/me); 電子電荷qe=-1.6021865314e-19; 磁場(chǎng)B0=5.3;lamd_e = me * ve / fabs(qe * B0);//電子回旋半徑lamd_D = m_D * vD / fabs(qi * B0);//氘離子回旋半徑//Ne=1e19;電子數(shù)密度 //真空介電常數(shù)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;時(shí)間步長(zhǎng) //dz_plas=lz_plas/double(Nz_plas-1);z方向空間步長(zhǎng)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("德拜長(zhǎng)度=%10.6e,空間步長(zhǎng)=%10.6e\n", deby_leng, dz_plas);printf("電子振蕩周期=%10.6e,時(shí)間步長(zhǎng)=%10.6e\n", peri_e, dt);printf("wpe*dt=%10.6f,v*dt/dz=%10.6f\n", wpe_dt, vdt_dz);}

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

?

?

?

總結(jié)

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

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。