JGonzalez/vlaplex
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Minor changes: Rename T_e to T_e0

Browse source
This commit is contained in:
Jorge Gonzalez 2025-07-03 20:48:55 +02:00
commit 4d44d8adf0
1 changed files with 12 additions and 12 deletions
Download patch file Download diff file Expand all files Collapse all files

24
src/vlaplex.f90
View file

@ -53,7 +53,7 @@ program VlaPlEx
real(dp), allocatable, dimension(:):: phi, phi_old, E, db_dphi real(dp), allocatable, dimension(:):: phi, phi_old, E, db_dphi
real(dp):: phiConv real(dp):: phiConv
real(dp):: phi0 real(dp):: phi0
real(dp):: T_e real(dp):: T_e0
! real(dp):: phiF ! real(dp):: phiF
integer:: k integer:: k
@ -145,7 +145,7 @@ program VlaPlEx
E_i = 0.0_dp E_i = 0.0_dp
T_i = 0.0_dp T_i = 0.0_dp
n_e = 0.0_dp n_e = 0.0_dp
T_e = 0.0_dp T_e0 = 0.0_dp
Zave = 0.0_dp Zave = 0.0_dp
Zave_bc_old = 0.0_dp Zave_bc_old = 0.0_dp
phi = 0.0_dp phi = 0.0_dp
@ -211,7 +211,7 @@ program VlaPlEx
! Get boundary conditions for specific time ! Get boundary conditions for specific time
call boundaryConditions%get(time, n_bc, u_bc, Temp_bc) call boundaryConditions%get(time, n_bc, u_bc, Temp_bc)
! Find new \bar{Z}_i based on T_e = Temp_bc and n_e = n_bc ! Find new \bar{Z}_i based on T_e0 = Temp_bc and n_e = n_bc
call Tene_to_Z%get(Temp_bc, n_bc, Zave_bc) call Tene_to_Z%get(Temp_bc, n_bc, Zave_bc)
! Assign Z(T,n) to bin ! Assign Z(T,n) to bin
z_inj = minloc(abs(Zlist - Zave_bc),1) z_inj = minloc(abs(Zlist - Zave_bc),1)
@ -229,7 +229,7 @@ program VlaPlEx
f_i_old(z_inj,1,j0:nv) = f0 f_i_old(z_inj,1,j0:nv) = f0
f_i(:,1,j0:nv) = f_i_old(:,1,j0:nv) f_i(:,1,j0:nv) = f_i_old(:,1,j0:nv)
T_e = Temp_bc T_e0 = Temp_bc
! r = rf, v<0 ! r = rf, v<0
f_i_old(:,nr,1:j0-1) = 0.0_dp f_i_old(:,nr,1:j0-1) = 0.0_dp
@ -248,12 +248,12 @@ program VlaPlEx
! Advect negative velocity ! Advect negative velocity
if (i < nr) then if (i < nr) then
f_i(iz,i,1:j0-1) = f_i_old(iz,i,1:j0-1) - v(1:j0-1)*dt/dr/r(i)**2*(r(i+1)**2*f_i_old(iz,i+1,1:j0-1) - & f_i(iz,i,1:j0-1) = f_i_old(iz,i,1:j0-1) - v(1:j0-1)*dt/dr/r(i)**2*(r(i+1)**2*f_i_old(iz,i+1,1:j0-1) - &
r(i )**2*f_i_old(iz,i ,1:j0-1)) r(i )**2*f_i_old(iz,i ,1:j0-1))
end if end if
! Advect positive velocity ! Advect positive velocity
if (i > 1) then if (i > 1) then
f_i(iz,i,j0:nv) = f_i_old(iz,i, j0:nv) - v( j0:nv)*dt/dr/r(i)**2*(r(i )**2*f_i_old(iz,i , j0:nv) - & f_i(iz,i,j0:nv) = f_i_old(iz,i, j0:nv) - v( j0:nv)*dt/dr/r(i)**2*(r(i )**2*f_i_old(iz,i , j0:nv) - &
r(i-1)**2*f_i_old(iz,i-1, j0:nv)) r(i-1)**2*f_i_old(iz,i-1, j0:nv))
end if end if
n_i(iz,i) = sum(f_i(iz,i,:))*dv n_i(iz,i) = sum(f_i(iz,i,:))*dv
@ -308,12 +308,12 @@ program VlaPlEx
! phi0=phi(1) ! Neumann ! phi0=phi(1) ! Neumann
! Calculate distribution of electrons ! Calculate distribution of electrons
! n_e = sum_ni(1) * exp((phi- phi0) / T_e) ! Isothermal (Boltzmann) ! n_e = sum_ni(1) * exp((phi- phi0) / T_e0) ! Isothermal (Boltzmann)
n_e = sum_ni(1) * (1.0_dp + (gamma_e - 1.0_dp)/gamma_e*(phi-phi0)/T_e)**gamma_e_exp !Polytropic n_e = sum_ni(1) * (1.0_dp + (gamma_e - 1.0_dp)/gamma_e*(phi-phi0)/T_e0)**gamma_e_exp !Polytropic
! Diagonal matrix for Newton integration scheme ! Diagonal matrix for Newton integration scheme
! db_dphi = n_e / T_e ! Isothermal (Boltzmann) ! db_dphi = n_e / T_e0 ! Isothermal (Boltzmann)
db_dphi = sum_ni(1) / (gamma_e * T_e) * & db_dphi = sum_ni(1) / (gamma_e * T_e0) * &
(1.0_dp + (gamma_e - 1.0_dp)/gamma_e*(phi-phi0)/T_e)**gamma_e_dexp !Polytropic (1.0_dp + (gamma_e - 1.0_dp)/gamma_e*(phi-phi0)/T_e0)**gamma_e_dexp !Polytropic
! Check if the solution has converged ! Check if the solution has converged
phiConv = maxval(abs(Res),1) phiConv = maxval(abs(Res),1)
@ -324,7 +324,7 @@ program VlaPlEx
! ! Calculate new potential to ensure 0 current at the edge ! ! Calculate new potential to ensure 0 current at the edge
! if (n_i(nr) > n_epsilon) then ! if (n_i(nr) > n_epsilon) then
! phiF = phi0 + T_e * log((2.0_dp*sqrt(pi)*Zave(nr)*n_i(nr)*u_i(nr)) / (Zave(1)*n_i(1)*sqrt(m_i*T_e/m_e))) ! phiF = phi0 + T_e0 * log((2.0_dp*sqrt(pi)*Zave(nr)*n_i(nr)*u_i(nr)) / (Zave(1)*n_i(1)*sqrt(m_i*T_e0/m_e)))
! !
! else ! else
! phiF = phi(nr-5) ! phiF = phi(nr-5)