24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

xingxiao

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 101.5
帖子: 110
在线: 62.9小时
虫号: 712761
注册: 2009-03-02
专业: 理论和计算化学

[求助] 已算出整个空间的流，下一步想得到在空间中某一位置的流

我已经实现在整个放电空间的粒子流密度，接下来想要得到在电极极板处了流密度，请问怎样实现？？？
下面是我编译的整个放电空间的粒子流密度：
implicit none
integer,intent(in)::Ng
real,intent(in)::Sbc,gone,q
   integer,dimension(1:Ng),intent(in)::ink
real,dimension(1:Ng,1:3000,1:7),intent(in)::AN
real,dimension(0:Ng),intent(out)::Jf
integer::i,j,m
real::Vq,qV,qv1,qv2
real,parameter::br=0.064,pi=3.14159
do i=0,Ng  整个空间包含的网格数
   Jf(i)=0.0 整个空间的粒子流
end do
do i=1,Ng
   do j=1,ink(i) 在某一个网格里德粒子数目
   if(q<0.0) then    如果是负电荷
         qV=abs(AN(i,j,7))*AN(i,j,4)       其中AN(i,j,7)是粒子的电荷量，
                                                                           AN(i,j,4)是离子在轴向的速度
      else if(q>0.0) then                      如果是正电荷
         qV=-abs(AN(i,j,7))*AN(i,j,4)
      end if
   call particleweighting(AN(i,j,1),Sbc,qV,m,qv1,qv2) 粒子加权到网格点上
      if(m==Ng) then    如果在最后一个网格，即右极板处
         Jf(m)=Jf(m)+qv1
      else
         Jf(m)=Jf(m)+qv1
         Jf(m+1)=Jf(m+1)+qv2
      end if
      end if
   end do
end do
Vq=abs(q)*pi*br**2*Sbc  每个网格的体积
   Jf(0)=Jf(0)/(0.5*Vq);Jf(Ng)=Jf(Ng)/(0.5*Vq)
do i=1,Ng-1
   Jf(i)=Jf(i)/Vq
end do

现在我想实现粒子在左极板处的流密度，它应该是一个数，而不是数组，但是上述算出来的Jf是一个数组，请问我该怎么办？？？

回复此楼

» 猜你喜欢

导师想让我从独立一作变成了共一第一已经有9人回复
博士读完未来一定会好吗已经有23人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有11人回复
读博已经有4人回复
JMPT 期刊投稿流程已经有4人回复
心脉受损已经有5人回复
Springer期刊投稿求助已经有4人回复
小论文投稿已经有3人回复
申请2026年博士已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼2012-02-06 11:14:04

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 xingxiao 的主题更新

返回列表