CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   helmholtzFoam / direct matrix solver (http://www.cfd-online.com/Forums/openfoam/78768-helmholtzfoam-direct-matrix-solver.html)

mirko August 2, 2010 03:30

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


All times are GMT -4. The time now is 13:38.