科学网

 找回密码
  注册

tag 标签: 粒子模拟

相关帖子

版块 作者 回复/查看 最后发表

没有相关内容

相关日志

等离子体粒子模拟及其应用 第11讲-粒子模拟的并行化
YUNJU 2020-5-26 19:59
个人分类: 科研笔记|1680 次阅读|0 个评论
[转载]等离子体粒子模拟及其应用 第56讲-电磁模型
YUNJU 2020-5-26 19:05
第一节 基本方程与无量纲化 第二节 粒子速度的推进 第三节 一维电磁模式(1-1/2) 第四节 哨声波的一维粒子模拟 二维电磁模型
个人分类: 科研笔记|1718 次阅读|0 个评论
等离子体粒子模拟及其应用课程大作业暨2020年第4次组会报告(附教学视频)
YUNJU 2020-5-26 00:07
目录 程序变量说明 nss:两种成分的电子 np:粒子数目 idimv:每个电子有两个坐标,三个速度 nx,ny:分别是x,y方向的网格点 nloop:计算的总步数 dt:时间步长 part(idimv,np,nss):所有电子坐标和速度的数组 dns(nss):两种电子成分的权重 part(idimv,np):所有离子坐标和速度的数组 vbe(idimp),vte(idimp):电子初始的平均速度和热速度 vbi(idimp),vti(idimp):离子初始的平均速度和热速度 fx(nxv,ny),fy(nxv,ny):网格上电场的x,y分量(实空间) fxc(nxv,ny),fyc(nxv,ny):网格上电场的x,y分量(复空间) qs(nxv,ny,nss):两种电子成分在网格上的电荷 q(nxv,ny):网格上的电荷密度(实空间) qc(nxv,ny):网格上的电荷密度(复空间) we:通过电势算的电场能量 eke:电子动能 eki:粒子动能 wt:系统中的总能量 fx:电场在x方向的分量 fy:电场在y方向的分量 wx:x方向上的电场能量 wy:y方向上的电场能量 验证数值不稳定性 网格自加热效应验证 束流不稳定性结果(附MATLAB代码) clc;clear all;load D:\\Fortran\\shuliu\\ee.dat; loop=15;Cxy=256;qq=zeros(Cxy,Cxy,loop); ex=zeros(Cxy,Cxy,loop);ey=zeros(Cxy,Cxy,loop); for i=1:loop for j=1:Cxy for k=1:Cxy temp=Cxy*(j-1)+k+Cxy*Cxy*(i-1);qq(j,k,i)=ee(temp,3); ex(j,k,i)=ee(temp,4);ey(j,k,i)=ee(temp,5); end end end x=1:Cxy;y=1:Cxy; =meshgrid(x,y) figure(1) for i=1:loop subplot(5,3,i);contour(X,Y,qq(:,:,i)) end figure(2) for i=1:loop subplot(5,3,i);contour(X,Y,ex(:,:,i)) end figure(3) for i=1:loop subplot(5,3,i);contour(X,Y,ey(:,:,i)) end 附教学视频(有两个大于4G的视频未上传): https://pan.baidu.com/s/1caKhAKCmdKRibdgF2aJJ0A
个人分类: 科研笔记|121 次阅读|1 个评论
从粒子到网格:粒子云分室法(附代码)
YUNJU 2020-5-23 20:20
例子:网格50*50,网格距为1,每个网格里有均匀分布的100个粒子,粒子电量为1,使用粒子云分室法计算格点的电荷量。 结果图 MATLAB程序 function =electron_density() q=zeros(51,51); for i=1:50 for j=1:50 for m=1:10 for n=1:10 q(i,j)=q(i,j)+(11-m)*(11-n)/(11*11); q(i+1,j)=q(i+1,j)+m*(11-n)/(11*11); q(i,j+1)=q(i,j+1)+(11-m)*n/(11*11); q(i+1,j+1)=q(i+1,j+1)+m*n/(11*11); end end end end raw=zeros(1,51);column=zeros(52,1) q_raw= ;q_column= figure(1);pcolor(q_column);colorbar; end
个人分类: 科研笔记|1882 次阅读|0 个评论
[转载]等离子体粒子模拟及应用 第234讲-静电模型
YUNJU 2020-5-12 13:11
第一节基本方程 计算步骤 第二节 粒子的处理方法 泊松方程的求解 第三节 粒子移动格式的精度 第四 节 网格-粒子云-蛙跳动力论 例子Electrostatic Solitary Structures in 2D plasma 来自一维PIC仿真的电场的堆叠轮廓。 一维PIC模拟得出的电子的电场和相空间图。 多维的结构 2D PIC simulations电子双流不稳定性:最初两个电子成分具有相同的密度和温度
个人分类: 科研笔记|1689 次阅读|0 个评论
等离子体粒子模拟及应用 第1讲-绪论
YUNJU 2020-5-10 18:59
参考书目 主要内容 第一章绪论 等离子体维持的条件(易忽略验证) 等离子体模拟的方法(上学期学了基础的MHD模拟) 粒子模拟方法概述 和MHD模拟的不同之处 粒子模拟历史 粒子模拟的一些概念 粒子模拟的基本方程 有限粒子大小对碰撞的影响 一些定量的结果 原文链接: https://asipp.yunjunet.cn/58232.html
个人分类: 科研笔记|2178 次阅读|0 个评论

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-6-16 17:37

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部