Commit d16af64e authored by Caneses Marin, Juan Francisco's avatar Caneses Marin, Juan Francisco
Browse files

Using the mean of Tpar and Tper for collision operator. using vx insread of vperp in P22

parent e8f03111

src/CoulombCollisions.f90

+15 −15
Original line number Diff line number Diff line
@@ -119,7 +119,7 @@ np0 = plasma%np(i) 
!np0  = params%ne0

Up0  = plasma%Up(i)

!Tp0  = params%Te0

Tp0  = abs(plasma%Tparp(i))

Tp0  = 0.5*abs(plasma%Tparp(i) + plasma%Tperp(i))





! Test particle properties:

src/MoveParticlePack.f90

+8 −5
Original line number Diff line number Diff line
@@ -363,7 +363,7 @@ TYPE(paramsTYP), INTENT(IN) :: params 
! Define local variables:

INTEGER(i4) :: i, ix, frame, NZ

INTEGER(i4), DIMENSION(3) :: ixLeft, ixRight

REAL(r8) :: vpar, Ma, Ep, alpha, xip, vper, v, a

REAL(r8) :: vpar, Ma, Ep, alpha, xip, vper, v, a, phi, vx

REAL(r8), DIMENSION(mesh%NZmesh + 4) :: n, nU, unU, P11, P22, nUE, cf



! Initialize local mesh quantities:
@@ -395,6 +395,9 @@ DO i = 1,params%NC 
		xip  = plasma%xip(i)

		vpar = v*xip

		vper = v*sqrt(1 - xip**2.)

		call random_number(phi)

		phi = phi*2*pi

		vx = vper*cos(phi)



		! Mesh point with ghost cells:

		ix = plasma%m(i) + 2
@@ -420,9 +423,9 @@ DO i = 1,params%NC 
		P11(ix+1)  = P11(ix+1)  + plasma%wR(i)*a*Ma*vpar**2.



		! Stress tensor 22:

		P22(ix-1)  = P22(ix-1)  + plasma%wL(i)*a*Ma*vper**2.

		P22(ix)    = P22(ix)    + plasma%wC(i)*a*Ma*vper**2.

		P22(ix+1)  = P22(ix+1)  + plasma%wR(i)*a*Ma*vper**2.

		P22(ix-1)  = P22(ix-1)  + plasma%wL(i)*a*Ma*vx**2.

		P22(ix)    = P22(ix)    + plasma%wC(i)*a*Ma*vx**2.

		P22(ix+1)  = P22(ix+1)  + plasma%wR(i)*a*Ma*vx**2.



		! Energy flux density:

		nUE(ix-1)  = nUE(ix-1)  + plasma%wL(i)*a*vpar*Ep

src/run_linFP.sh

+1 −1
Original line number Diff line number Diff line
@@ -9,7 +9,7 @@ INPUT_FILE="xp_TestDrag.in" 


# Set number of threads:

# ===================================================

export OMP_NUM_THREADS=48

export OMP_NUM_THREADS=96



# Set processor binding for openMP:

# ===================================================