Skip to content

Gk sheath bc with mu dependent vcut - #95

Open
Antoinehoff wants to merge 43 commits into
mainfrom
gk_sheath_vcut_mu
Open

Gk sheath bc with mu dependent vcut#95
Antoinehoff wants to merge 43 commits into
mainfrom
gk_sheath_vcut_mu

Conversation

@Antoinehoff

Copy link
Copy Markdown
Collaborator

This PR is related to these Gkeyll PR and DR

…ernels, add_apar and add_apardot kernels. We conserve the electrostatic part of the scripts.
…ity of generating separated EM/ES kernels. This will be cleaned once the design in Gkeyll is fully decided.
…lock instance. This did not change a iota of the output kernels but better be safe than sorry.

There is an error in the addapardot volume kernel I think because my first tests explode when I include it, will work on that tomorrow 🤓 have a good night folks 😘
…s much as possible the precomputed geometric quantities like normcurlbhat_quad for surface flux and bioverJB_list for the volume computation. A comment in each script is added to point out the terms responsible to oscillations observed in the 3x2v Alfven testcase. The volume term builds an oscillation in n_e along y at longer time whereas the surface flux term creates a discontinuity in n_e along y at the first frame of the regression test.

It is still unknown if these terms are the source of the 3x2v alfven issue or if they are triggered unphysically from other part of the code like the Helmholtz solver.
Comment thread maxima/g0/gk-sheath/ms-gk-sheath.mac Outdated
Comment thread maxima/g0/gk_collisionless/verify_no_leak.mac Outdated
Comment thread maxima/g0/gk-sheath/sheath-1x2v.mac
Comment thread maxima/g0/gk-sheath/sheath-2x2v.mac
Comment thread maxima/g0/gk-sheath/sheath-3x2v.mac
Comment thread maxima/g0/gk_collisionless/dg_gk-vol.mac Outdated
Comment thread maxima/g0/gk_collisionless/ms-dg_gyrokinetic-header.mac Outdated
Comment thread maxima/g0/gk_collisionless/ms-dg_gyrokinetic-vol.mac Outdated
Comment thread maxima/g0/gk_collisionless/ms-gk_collisionless_flux-header.mac Outdated
Comment thread maxima/g0/gk_collisionless/ms-gk_collisionless_flux.mac Outdated
@Antoinehoff

Copy link
Copy Markdown
Collaborator Author

thanks @manauref for your review! I think I've addressed everything, only a small question if float(expan( instead of float(fullratsimp( is really better. I checked the unit test and it still passes with the new version so it's an aesthetic decision I guess.

@Antoinehoff
Antoinehoff requested a review from manauref May 15, 2026 13:28
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