helmholtzFoam / direct matrix solver
 

Hello,

I am writing a solver for the Helmholtz equation:

fvm::laplacian(T) + k2 * T

Right now, I am solving it in one dimension only.

I used the laplaceFoam solver code and tutorial as guidance.

I am having divergence problems using the PCG solver with the DIC conditioner. From iteration to iteration, the magnitude of T keeps increasing.

I can see that the k2 term influences the diagonal term. I would like to try solving the equation first with a direct matrix solver (since this is currently a one-dimensional problem), but I was unable to find any direct solvers.

Does OF have any direct solvers? I looked at the src directories, but nothing obvious jumped out to my untrained eyes.

Thanks,

Mirko