helmholtzFoam / direct matrix solver
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.
|All times are GMT -4. The time now is 17:37.|