data-raw/j0j0_examples.R

In mstejner/j0j0r: Calculate the j0j0 terms in the fluctuation model of Bindslev 1996

#script to generate j0j0_examples, distribution_examples
devtools::uninstall("C:/ctsr/j0j0r")
devtools::document("C:/ctsr/j0j0r")
devtools::install("C:/ctsr/j0j0r")

library(rlang)
library(magrittr)

phi <- 80
k <- 2 * pi / (j0j0r::const$c / 100e9)
B = 2.5
n = 4e19
A = 2
Z = 1

future::plan(strategy = "multisession", .skip = TRUE)

maxwellian = j0j0r::maxwellian_setup(
  n = n,
  T_eV = 2000,
  A = A,
  Z = Z,
  name = "maxwellian"
)

maxwellian_example <- j0j0r::j0j0(
  k = 2 * pi / (j0j0r::const$c / 100e9),
  phi = c(60, 80, 86),
  frequencies = seq(0, 400e6, by = 200e6),
  B = 2.5,
  particles = list(
    maxwellian = maxwellian
  ),
  integration_method = "hcubature"
)
usethis::use_data(maxwellian_example, overwrite = TRUE)


bimaxwellian <-  j0j0r::bimaxwellian_setup(
  n = n,
  T_eV_perp = 2000,
  T_eV_par = 1000,
  v_drift = 2e5,
  A = A,
  Z = Z,
  name =  "bimaxwellian"
)

lorentzian <- j0j0r::generalized_lorentzian_setup(
  n = n,
  T_eV = 2000,
  kp = 5,
  A = A,
  Z = Z,
  name = "lorentzian"
)

ring = j0j0r::maxwellian_ring_setup(
  n = n,
  v_width = 1.5e5,
  v_rad = 0.5e6,
  A = A,
  Z = Z,
  name = "ring"
)


center <- c(5e5, 3e5)
T_eV <- 2e3
v_term <- j0j0r::find_p_term(T_eV, A) / (A * j0j0r::const[["amu"]])
covariance <- rbind(c(v_term^2/4, v_term^2/3), c(v_term^2/5, (v_term)^2))

bvtnorm <- j0j0r::bvtnorm_setup(
  n = 4e19,
  center = center,
  covariance = covariance,
  A = A,
  Z = Z,
  name = "bivariate_normal"
)


slowdown_b5 = j0j0r::slowdown_setup(
  b = 5,
  n = n,
  A = A,
  Z = Z,
  birth_energy = 20e3,
  n_e = n,
  T_e_eV = 2e3,
  ions = data.frame(
    Z = c(1),
    A = c(2),
    n = n
  ),
  name = "deuterium_slowdown_b5"
)

slowdown_b0 = j0j0r::slowdown_setup(
  b = 0,
  n = n,
  A = A,
  Z = Z,
  birth_energy = 20e3,
  n_e = n,
  T_e_eV = 2e3,
  ions = data.frame(
    Z = c(1),
    A = c(2),
    n = n
  ),
  name = "deuterium_slowdown_b0"
)

distribution_examples <- j0j0r::calculate_distribution_data_frame(
  particles = list(
    maxwellian = maxwellian,
    bimaxwellian = bimaxwellian,
    lorentzian = lorentzian,
    ring = ring,
    bvtnorm = bvtnorm,
    slowdown_b5 = slowdown_b5,
    slowdown_b0 = slowdown_b0
  ),
  v_par = seq(-2.5e6, 2.5e6, length.out = 1000),
  v_perp = 0
)

dist_2D_examples <-
  dplyr::bind_rows(
    j0j0r::calculate_distribution_data_frame(
      particles = list(
        maxwellian = maxwellian,
        bimaxwellian = bimaxwellian,
        lorentzian = lorentzian,
        ring = ring,
        bvtnorm = bvtnorm
      ),
      v_par = seq(-1e6, 1e6, length.out = 50),
      v_perp = seq(0, 1e6, length.out = 50)
    ),
    j0j0r::calculate_distribution_data_frame(
      particles = list(
        slowdown_b5 = slowdown_b5,
        slowdown_b0 = slowdown_b0
      ),
      v_par = seq(-1.5e6, 1.5e6, length.out = 300),
      v_perp = seq(0, 1.5e6, length.out = 300)
    )
  )
usethis::use_data(dist_2D_examples, overwrite = TRUE)


j0j0_examples <- j0j0r::j0j0(
  k = k,
  phi = phi,
  frequencies = seq(0, 600e6, by = 2e6),
  B = B,
  particles = list(
    maxwellian = maxwellian,
    bimaxwellian = bimaxwellian,
    lorentzian = lorentzian,
    ring = ring,
    bvtnorm = bvtnorm,
    slowdown_b0 = slowdown_b0,
    slowdown_b5 = slowdown_b5
  ),
  integration_method = "stats"
)
usethis::use_data(j0j0_examples, distribution_examples, overwrite = TRUE)


phi_sweep <- j0j0r::j0j0(
  k = k,
  phi = c(60, 70, 80, 84, 86),
  frequencies = seq(0, 600e6, by = 2e6),
  directions = c("z"),
  B = B,
  particles = list(
    maxwellian = maxwellian,
    ring = ring,
    slowdown_b5 = slowdown_b5
  ),
  integration_method = "stats"
)
phi_sweep[["phi"]] <- as.factor(phi_sweep[["phi"]])
phi_sweep %>% plot_j0j0(wrap_by = "particle", color_by = "phi")

usethis::use_data(phi_sweep, overwrite = TRUE)