search

ammarhakim / gkyl

This is the main source repo for the Gkeyll 2.0 code. Please see gkeyll.rtfd.io for details.
https://gkeyll.readthedocs.io/en/latest/
54 stars 14 forks source link

For Poisson solve, output is showing electric potential for all 4 components. Last three components should be magnetic potential (0 for no B field case) #149

Closed crskolar closed 8 months ago

crskolar commented 8 months ago

I'm running sims with this field solver in the input file: 

   -- Field solver.
   field = Plasma.Field {
      epsilon0 = epsilon_0, mu0 = mu_0,
      evolve   = true, -- Evolve field?
      hasMagneticField = false,
      bcLowerPhi = {{T="D", V=V_left}},
      bcUpperPhi = {{T="D", V=0.0}},
   },

This should output a 4 component output with phi and 3 components of 0s (since no B field being included). However, all the outputs are the electric potential.

manauref commented 8 months ago

@crskolar, it's an electrostatic simulation without external B. Why should it output a (zero) vector potential?

JunoRavin commented 8 months ago

The output in field is (phi_tot, A_ext) always no? It's just the last 3 are zero in a pure Poisson calculation

manauref commented 8 months ago

@crskolar I think you are in an older version of the code please pull the latest in gkylzero and gkyl/g0-merge and build again, preferably from scratch. If this fixes it please close the issue.

crskolar commented 8 months ago

@manauref You were correct. g0 and g0-merge branch fixes it. I'll work with Kolter to get the emissionBC branches all merged in so that I don't keep having these issues.

I'm gonna close this issue now.