search

mabarnes / moment_kinetics

Other
2 stars 4 forks source link

1D2V simulations with electric field switch #268

Open mrhardman opened 5 days ago

mrhardman commented 5 days ago

To support Mantas Abazorius's work, it would be useful to implement a switch to turn off Ez so that the DKE solved for is 

\frac{\partial F}{\partial t} + v_{\|}\frac{\partial F}{\partial z} = C[F,F]

@johnomotani @LucasMontoya4 Where should I implement the switch? Do you agree that this could be an interesting feature?

LucasMontoya4 commented 5 days ago

I'm not sure, but maybe it could belong in the composition section of the input file, considering this would be entirely different physics from kinetic or fluid ions? So I guess one of the new flags for ion_physics?

mrhardman commented 5 days ago

I think this is orthogonal to all other options if one just sets a flag that is used in em_fields.jl to set Ez and Er to zero.

LucasMontoya4 commented 5 days ago

Sounds much more sensible!

johnomotani commented 4 days ago

I'd be tempted to add this as an electron_physics flag - because electron_physics tells us how to calculate phi and Ez. Could we argue that this switch is equivalent to massless, zero-temperature electrons? So call it something like zero_temperature?

mrhardman commented 4 days ago

I view this as a physics hack for debugging and testing rather than a self-consistent model. The present version of the PR would let you calculate everything everywhere, just with Ez= Er = 0 https://github.com/mabarnes/moment_kinetics/pull/269. However, if you are sure that this is a physical limit, then I have no objection to your suggestion to make this as an electron model. Perhaps boltzmann_zero_temperature would make it clear that the electron response must be trivial in this option (if it is a real limit)? @MantasAbazorius What do you think?