#include <stdio.h> #include "define.h" /* 坐標變換 Michael Komm,Interakce plazmatu se stenou tokamaku,2007,pp49 */ void init(particle* ptr_e, particle* ptr_D, particle* ptr_C3) {double dz_e, dz_D, dz_C3;double vx, vy, vz;double v_para, v_ver1, v_ver2;int i;vx = 0.0;vy = 0.0;vz = 0.0;dz_e = lz_plas / (numb_e + 1); //均勻放置粒子,不在邊界放置粒子dz_D = lz_plas / (numb_D + 1);dz_C3 = lz_plas / (numb_C3 + 1);/* 初始化電子 */for (i = 0; i < numb_e; i++){maxw(&vx, &vy, &vz, Te, me);ptr_e[i].z = dz_e * (i + 1);ptr_e[i].vx = vx;ptr_e[i].vy = vy;ptr_e[i].vz = vz;ptr_e[i].ener = 0.5 * (vx * vx + vy * vy + vz * vz) * me / qi;ptr_e[i].weig = weig_e;if (i == numb_e - 1)printf("z_e=%10.6e,", ptr_e[i].z);}/* 初始化氘離子 */for (i = 0; i < numb_D; i++){/* 存在定向速度的情況maxw(&v_para,&v_ver1,&v_ver2,Ti,m_D);vx=(v_para-vb_D)*bx-(v_ver1*by+v_ver2*bx*bz)/sqrt(bx*bx+by*by); //坐標系變換vy=(v_para-vb_D)*by+(v_ver1*bx-v_ver2*by*bz)/sqrt(bx*bx+by*by);vz=(v_para-vb_D)*bz+v_ver2*sqrt(bx*bx+by*by); */maxw(&vx, &vy, &vz, Ti, m_D);ptr_D[i].z = dz_D * (i + 1);ptr_D[i].vx = vx;ptr_D[i].vy = vy;ptr_D[i].vz = vz;ptr_D[i].ener = 0.5 * (vx * vx + vy * vy + vz * vz) * m_D / qi;ptr_D[i].weig = weig_D;if (i == numb_D - 1)printf("z_D=%10.6e,", ptr_D[i].z);}/* 初始化C3+離子 */for (i = 0; i < numb_C3; i++){/* 存在磁場及定向速度maxw(&v_para,&v_ver1,&v_ver2,Tc,m_C);vx=(v_para-vb_C)*bx-(v_ver1*by+v_ver2*bx*bz)/sqrt(bx*bx+by*by); //坐標系變換vy=(v_para-vb_C)*by+(v_ver1*bx-v_ver2*by*bz)/sqrt(bx*bx+by*by);vz=(v_para-vb_C)*bz+v_ver2*sqrt(bx*bx+by*by); */maxw(&vx, &vy, &vz, Tc, m_C);ptr_C3[i].z = dz_C3 * (i + 1);ptr_C3[i].vx = vx;ptr_C3[i].vy = vy;ptr_C3[i].vz = vz;ptr_C3[i].ener = 0.5 * (vx * vx + vy * vy + vz * vz) * m_C / qi;ptr_C3[i].weig = weig_C3;if (i == numb_C3 - 1)printf("z_C3=%10.6e\n", ptr_C3[i].z);}}

?

maxw.cpp

#include <stdlib.h> #include <math.h> #include "define.h" #include<stdio.h> /* James Joseph Szabo.Fully Kinetic Numerical Modeling of a Plasma Thruster.2011,p149 */void maxw(double* vx, double* vy, double* vz, double kt, double mass) {double ran1, vt, v;double theta, psi;/* 粒子熱速度 */vt = sqrt(2.0 * kt * qi / mass);//kt粒子溫度，qi離子電荷，mass粒子質(zhì)量//printf("電子熱速度%f \n", vt); ran1 = rand() / (double)RAND_MAX;//產(chǎn)生0-1之間的隨機數(shù)v = vt * sqrt(-log(ran1));//?ran1 = rand() / (double)RAND_MAX;theta = 2.0 * pi * ran1;ran1 = rand() / (double)RAND_MAX;psi = 2.0 * pi * ran1;*vx = v * sin(theta);*vy = v * cos(theta);*vz = v * sin(psi); //? }

總結

以上是生活随笔為你收集整理的[C++调试笔记]初始化3种粒子数据initmaxw的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

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