In[]:=
Clear["Global`*"]
In[]:=
​
In[]:=
P=0.18Pa;ρ=88.7
-3
cm
;T=190K;kb=1.38*
-23
10
kg*
2
m
*
-2
s
*
-1
K
;​​m=kb*rho/P*T
$RecursionLimit
:Recursion depth of 1024 exceeded.
Out[]=
1.45667×
-20
10
kgrho
2
TerminatedEvaluation[RecursionLimit]
Pa
2
s
In[]:=
​​
In[]:=
La=3*10
-2
cm
*
-1
s
;​​σ=2.02*
-18
10
2
cm
;​​nNO=1*
8
10
-3
cm
;​​deltaz=900km;​​H=10km;​​k=3.5*
-7
10
3
cm
;​​ne=5.23*
5
10
-3
cm
;​​ni=5.09*
5
10
-3
cm
;​​​​​​(*Ionproduction*)​​N[q=La*σ*nNO*Exp[-(deltaz/H)-1]]​​​​(*Ionloss*)​​N[l=k*ne*ni]​​
Out[]=
1.82673×
-48
10
3
cm
s
Out[]=
93172.5
3
cm
In[]:=
​
In[]:=
Plot
93172.5
3
cm
,{cm,-8,8}
Out[]=
-5
5
-15000
-10000
-5000
5000
10000
15000
In[]:=
N[q=La*σ*nNO]
Out[]=
6.06×
-9
10
3
cm
s
In[]:=
Plot3D
6.06×
-9
10
3
cm
s
,{cm,-1.01099,1.01099},{s,-1.,1.}
Out[]=
In[]:=
​​​​​​​​Ti=1129.8K;Te=3500K;mi=1.66*
-27
10
kg*16;g=9.8m*
-2
s
;me=9.109382*
-31
10
kg;
k
;​​k=1.38*
-23
10
kg*
2
m
*
-2
s
*
-1
K
;​​(*k=Quantity["BoltzmannConstant"]*)​​​​N[Hi=k*Ti/mi*g]​​N[Hi=k*Te/me*g]
$RecursionLimit
:Recursion depth of 1024 exceeded.
$RecursionLimit
:Recursion depth of 1024 exceeded.
Out[]=
5.75279×
6
10
3
TerminatedEvaluation[RecursionLimit]
4
s
Out[]=
5.19618×
11
10
3
TerminatedEvaluation[RecursionLimit]
4
s