From 73cb0a4afe889d49653621bd4e16f7d9c5997d2b Mon Sep 17 00:00:00 2001 From: JGonzalez Date: Mon, 18 Nov 2024 15:16:44 +0100 Subject: [PATCH] Trying to fix polytropic electrons. Still not working... --- plasmaExpansion.f90 | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/plasmaExpansion.f90 b/plasmaExpansion.f90 index 47ac293..cd103d1 100644 --- a/plasmaExpansion.f90 +++ b/plasmaExpansion.f90 @@ -286,6 +286,8 @@ program plasmaExpansion ! Diagonal matrix for Newton integration scheme db_dphi = n_e / T_e ! Isotropic + ! db_dphi = Zave(1) * n_i(1) / (gamma_e * T_e) * & + ! (1.0_dp + ((gamma_e - 1.0_dp)/gamma_e*(phi_old-phi0)/T_e)**((2.0_dp - gamma_e)/(gamma_e - 1.0_dp))) ! Polytropic diag = -2.0_dp / dr**2 - db_dphi diag_low = 1.0_dp / dr**2 - 1.0_dp / (r(2:nr) * dr) diag_high = 1.0_dp / dr**2 + 1.0_dp / (r(1:nr-1) * dr) @@ -313,7 +315,7 @@ program plasmaExpansion ! Calculate distribution of electrons n_e = Zave(1) * n_i(1) * exp((phi- phi0) / T_e) ! Isothermal (Boltzmann) - ! n_e = Zave(1) * n_i(1) * (1.0_dp + ((gamma_e - 1.0_dp)/gamma_e*(phi-phi0)/T_e)**(1.0_dp/(gamma_e - 1.0_dp))) ! Not working + ! n_e = Zave(1) * n_i(1) * (1.0_dp + ((gamma_e - 1.0_dp)/gamma_e*(phi-phi0)/T_e)**(1.0_dp/(gamma_e - 1.0_dp))) ! Polytropic ! Check if the solution has converged phiConv = maxval(abs(Res),1)