CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM (
-   -   helmholtzFoam / direct matrix solver (

mirko August 2, 2010 03:30

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 23:48.