diff --git a/bc_80ns_T10Z6.csv b/bc_80ns_T10Z6.csv new file mode 100644 index 0000000..04bccba --- /dev/null +++ b/bc_80ns_T10Z6.csv @@ -0,0 +1,5 @@ +t (s),n (m^-3),u (m s^-1),T (eV),Z +0.000000E-000,1.111111E+025,0.000000E+000,1.000000E+001,6.000000E+000 +8.000000E-008,1.111111E+025,0.000000E+000,1.000000E+001,6.000000E+000 +8.500000E-008,1.000000E+020,0.000000E+000,5.000000E+000,6.000000E+000 +2.000000E-006,1.000000E+020,0.000000E+000,5.000000E+000,6.000000E+000 diff --git a/bc_80ns_T30Z11.csv b/bc_80ns_T30Z11.csv new file mode 100644 index 0000000..a23d865 --- /dev/null +++ b/bc_80ns_T30Z11.csv @@ -0,0 +1,5 @@ +t (s),n (m^-3),u (m s^-1),T (eV),Z +0.000000E-000,1.111111E+025,0.000000E+000,3.000000E+001,1.100000E+001 +8.000000E-008,1.111111E+025,0.000000E+000,3.000000E+001,1.100000E+001 +8.500000E-008,1.000000E+020,0.000000E+000,5.000000E+000,1.100000E+001 +2.000000E-006,1.000000E+020,0.000000E+000,5.000000E+000,1.100000E+001 diff --git a/bc_80ns_T60Z16.csv b/bc_80ns_T60Z16.csv new file mode 100644 index 0000000..63c934d --- /dev/null +++ b/bc_80ns_T60Z16.csv @@ -0,0 +1,5 @@ +t (s),n (m^-3),u (m s^-1),T (eV),Z +0.000000E-000,1.111111E+025,0.000000E+000,6.000000E+001,1.600000E+001 +8.000000E-008,1.111111E+025,0.000000E+000,6.000000E+001,1.600000E+001 +8.500000E-008,1.000000E+020,0.000000E+000,5.000000E+000,1.600000E+001 +2.000000E-006,1.000000E+020,0.000000E+000,5.000000E+000,1.600000E+001 diff --git a/bc_80ns_T6Z4.csv b/bc_80ns_T6Z4.csv new file mode 100644 index 0000000..f0bdb43 --- /dev/null +++ b/bc_80ns_T6Z4.csv @@ -0,0 +1,5 @@ +t (s),n (m^-3),u (m s^-1),T (eV),Z +0.000000E-000,1.111111E+025,0.000000E+000,6.000000E+000,4.000000E+000 +8.000000E-008,1.111111E+025,0.000000E+000,6.000000E+000,4.000000E+000 +8.500000E-008,1.000000E+020,0.000000E+000,5.000000E+000,4.000000E+000 +2.000000E-006,1.000000E+020,0.000000E+000,5.000000E+000,4.000000E+000 diff --git a/plasmaExpansion.f90 b/plasmaExpansion.f90 index 896629d..7aa0251 100644 --- a/plasmaExpansion.f90 +++ b/plasmaExpansion.f90 @@ -219,7 +219,8 @@ program plasmaExpansion call setTimeFormat(nt) t = 0 call writeOutputRef() - call writeOutputF(t, dt, nr, r, nv, v, f_i_old) + ! call writeOutputF(t, dt, nr, r, nv, v, f_i_old) + call writeOutputFCum(t, dt, r(rCum_index), nv, v, fCum_i) call writeOutputPhi(t, dt, nr, r, phi, E, n_e) call writeOutputMom(t, dt, nr, r, n_i, u_i, T_i, Zave) @@ -228,7 +229,7 @@ program plasmaExpansion time = t * dt + t0 call boundaryConditions%get(time, n_bc, u_bc, Temp_bc, Zave_bc) call writeOutputBoundary(t, dt, n_bc, u_bc, Temp_bc, Zave_bc) - u_bc = sqrt(Zave_bc * Temp_bc) + ! u_bc = sqrt(Zave_bc * Temp_bc) ! f0(j0:nv) = v(j0:nv)**2 / sqrt(PI*Temp_bc**3) * exp(-(v(j0:nv) - u_bc)**2 / Temp_bc) f0(j0:nv) = 1.0_dp / sqrt(PI*Temp_bc) * exp(-(v(j0:nv) - u_bc)**2 / Temp_bc) f0 = f0 * n_bc / (sum(f0)*dv) @@ -326,7 +327,7 @@ program plasmaExpansion ! ! Calculate new potential to ensure 0 current at the edge ! if (n_i(nr) > 1.0e-10_dp) then - ! phiF = phi0 + T_i(1) * log((2.0_dp*sqrt(pi)*Zave(nr)*n_i(nr)*u_i(nr)) / (Zave(1)*n_i(1)*sqrt(m_i*T_i(1)/m_e))) + ! 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))) ! ! else ! phiF = phi(nr-5) @@ -386,7 +387,7 @@ program plasmaExpansion ! Write output if (mod(t,everyOutput) == 0 .or. t == nt) then - call writeOutputF(t, dt, nr, r, nv, v, f_i_old) + ! call writeOutputF(t, dt, nr, r, nv, v, f_i_old) call writeOutputPhi(t, dt, nr, r, phi, E, n_e) call writeOutputMom(t, dt, nr, r, n_i, u_i, T_i, Zave) call writeOutputFCum(t, dt, r(rCum_index), nv, v, fCum_i) diff --git a/scripts_python/plotCumF.py b/scripts_python/plotCumF.py index e0896b5..8a72112 100644 --- a/scripts_python/plotCumF.py +++ b/scripts_python/plotCumF.py @@ -11,18 +11,20 @@ m_i = 1.9712e-25 # paths = ['../quasiNeutral_partialAblation/','../Poisson_partialAblation/','../PoissonTi_partialAblation/','../2024-11-05_13.56.23/'] # paths = ['../2024-11-11_13.58.54/'] # paths = ['../quasiNeutral_fullAblation/','../Poisson_fullAblation/', '../quasiNeutral_partialAblation/', '../Poisson_partialAblation/'] -paths = ['../2024-11-28_10.01.56/'] +# paths = ['../2024-12-02_21.07.52/', '../Poisson_50ns_T30Z11/'] +# paths = ['../2024-12-10_18.45.17/', '../Poisson_50ns_T30Z11/'] +paths = ['../2024-12-11_12.38.27/', '../Poisson_polytropic_fa_T30Z11/', '../Poisson_fa_T30Z11/'] # paths = ['../Poisson_partialAblation/','../Poisson_partialAblation_lowerT/','../Poisson_partialAblation_lowT/','../Poisson_partialAblation_highT/'] labels = [path[3:-1] for path in paths] for path, label in zip(paths, labels): filesCum_i = sorted(glob.glob(path+'time_*_fCum_i.csv')) - start = 0 - end = len(filesCum_i) - every = 20 - for fileCum_i in filesCum_i[start:end+1:every]: - time, x, v, f_i = readF.read(fileCum_i) - plt.plot(v**2*m_i*0.5/e, f_i[0]*e/m_i/v, label='t = {:.1f} ns'.format(time*1e9)) + # start = 0 + # end = len(filesCum_i) + # every = 20 + # for fileCum_i in filesCum_i[start:end+1:every]: + # time, x, v, f_i = readF.read(fileCum_i) + # plt.plot(v**2*m_i*0.5/e, f_i[0]*e/m_i/v, label='t = {:.1f} ns'.format(time*1e9)) time, x, v, f_i = readF.read(filesCum_i[-1]) plt.plot(v**2*m_i*0.5/e, f_i[0]*e/m_i/v, label=label) diff --git a/scripts_python/plotMom.py b/scripts_python/plotMom.py index 95ea48a..f980681 100644 --- a/scripts_python/plotMom.py +++ b/scripts_python/plotMom.py @@ -11,23 +11,23 @@ from scipy.constants import e, k # paths = ['../quasiNeutral_partialAblation/','../Poisson_partialAblation/'] # paths = ['../2024-10-02_14.30.44/'] # paths = ['../quasiNeutral_fullAblation/','../Poisson_fullAblation/'] -paths = ['../2024-11-28_10.01.56/'] +paths = ['../2024-12-10_18.45.17/'] labels = [path[3:-1] for path in paths] for path, label in zip(paths, labels): filesPhi = sorted(glob.glob(path+'time_*_phi.csv')) filesMom_i = sorted(glob.glob(path+'time_*_mom_i.csv')) - start = 0 - end = 50#len(filesMom_i) - every = 10 + start = 50 + end = len(filesMom_i) + every = 20 fig, ax = plt.subplots(4, sharex='all') for fileMom_i, filePhi in zip(filesMom_i[start:end+1:every], filesPhi[start:end+1:every]): time, r, phi, E, n_e = readPhi.read(filePhi) time, r, n_i, u_i, T_i, Zave = readMom.read(fileMom_i) - # ax[0].plot(r, phi, label='t = {:.1f} ns'.format(time*1e9)) + ax[0].plot(r, phi, label='t = {:.1f} ns'.format(time*1e9)) # ax[0].plot(r, E, label='t = {:.1f} ns'.format(time*1e9)) - ax[0].plot(r, (Zave*n_i - n_e), label='t = {:.1f} ns'.format(time*1e9)) + # ax[0].plot(r, (Zave*n_i - n_e), label='t = {:.1f} ns'.format(time*1e9)) ax[1].set_yscale('log') # ax[1].set_ylim(bottom=1e16) ax[1].plot(r, Zave*n_i)