Skip to content

Conversation

@Maxwell-Rosen
Copy link
Collaborator

@Maxwell-Rosen Maxwell-Rosen commented Nov 19, 2025

$\mu_{\rm left} = \max_{i<k}\left(\frac{\frac{1}{2}m_s v_\parallel^2 + q_s(\phi(z) - \phi_{mi})}{B_{mi} - B(z)}\right)$
$\mu_{\rm right} = \max_{i>k}\left(\frac{\frac{1}{2}m_s v_\parallel^2 + q_s(\phi(z) - \phi_{mi})}{B_{mi} - B(z)}\right)$
$\mu = \min\left(\mu_{\rm left} ,\mu_{\rm right} \right)$

image

pgkyl gk_mirror_boltz_elc_poa_1x2v_p1-ion_fdot_multiplier_0.gkyl sel --z2 0.5 pl
image

… of B instead of the min/max operation. This lets us find the value of B at the wall as well
…he grid the position map operates on for the optimize function is always 3x but during the mask it's only 1x because of deflation
… functionality so that it gets the wall too. Pass the wall coordinate and Bmag into some updaters. It's valgrind clean on my laptop, which is why I'm making the commit here. Things should be passed into the right places and it's a good checkpoint. Changes are made to the cuda scripts, but they are not tested with gpu-sanitizeer
…rror throat. It's general for species of all charges. Valgrind clean. Update GPU kernels to compute this as well.
…lation. For some reason, the time step for the simulation in the OAP is still tiny. Must debug, but the cfl rate diagnostics look very reasonable and as we would expect
…at we actually should multiply it. Testing will show if this is a good ida. I think it makes sense because of the v_th in the omegaH formula.
…ding a regression test for looking at tandem mirrors. When the code is run in single mirror cell configurations, the tandem components are there to duplicate the outer cell. This implementation assumes symmetry in the magnetic field around z=0, which is a good assumption for the input files we are working with.
… symmetry check and modify mu_bound calculation in loss_cone_mask_gyrokinetic_advance for correct boundary handling. Tandem mirror testing is working
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants