boutproject / hermes Goto Github PK

Sorry not sure if this is an issue, but Div_par_diffusion seems to only calculate fluxes at the upper boundary:

https://github.com/boutproject/hermes/blob/master/div_ops.cxx#L529

What about the lower boundary?

I'm curious because I'm trying to use this operator elsewhere and it's causing some issues.

Cell centred differencing is susceptible to zig-zag / chequerboard modes unless some form of dissipation is used. Currently this is somewhat ad-hoc, and does not always work well.

Methods developed for SD1D should be applied to Hermes: minmod limiter and flux splitting for parallel advection. The skew-symmetric central differencing formulation may also be helpful for parallel current divergence, to suppress potential-jpar zig-zags.

Hi,

I am trying to run the Hermes and did 'run'. However, I've got a error as below. I would appreciate if you provide me a solution.

Thanks!

$ make
Compiling hermes-1.cxx
In file included from hermes-1.cxx:22:0:
hermes-1.hxx:30:10: fatal error: bout/invert/laplace3d.hxx: No such file or directory
#include <bout/invert/laplace3d.hxx>
^~~~~~~~~~~~~~~~~~~~~~~~~~~
compilation terminated.
../../make.config:318: recipe for target 'hermes-1.o' failed
make: *** [hermes-1.o] Error 1

Many compiler error:

hermes-1.cxx:298:21: error: ‘class Mesh’ has no member named ‘g11’
g11norm = mesh->g11 / averageY(mesh->g11);
^~~
hermes-1.cxx:298:42: error: ‘class Mesh’ has no member named ‘g11’
g11norm = mesh->g11 / averageY(mesh->g11);
^~~
hermes-1.cxx:383:33: error: ‘class Mesh’ has no member named ‘ngx’
for (int i = 0; i < mesh->ngx; i++) {
^~~
hermes-1.cxx:384:35: error: ‘class Mesh’ has no member named ‘ngy’
for (int j = 0; j < mesh->ngy; j++) {
^~~
hermes-1.cxx:398:15: error: ‘class Solver’ has no member named ‘addToRestart’; did you mean ‘addMonitor’?
solver->addToRestart(density_error_integral, "density_error_integral");
^~~~~~~~~~~~
addMonitor
hermes-1.cxx:401:15: error: ‘class Solver’ has no member named ‘addToRestart’; did you mean ‘addMonitor’?
solver->addToRestart(pe_error_integral, "pe_error_integral");
^~~~~~~~~~~~
addMonitor
hermes-1.cxx:436:11: error: ‘class Mesh’ has no member named ‘dx’
mesh->dx /= rho_s0 * rho_s0 * Bnorm;
^~
hermes-1.cxx:437:11: error: ‘class Mesh’ has no member named ‘Bxy’

The adaptive source currently tries to maintain the core profiles, but can result in source profiles which are peaked or otherwise contribute to numerical problems.

An adaptive source is needed which has a fixed profile, but whose amplitude is varied to maintain some measure of core density, temperature. Core boundary value? Or perhaps an average over a core region.