On Sep 19, 2008, at 5:11 PM, Zhiwen Liang wrote:
The higher-order discretization sounds interesting.
I knew I was going to regret even mentioning that as a possibility. A new discretization is not simple and will not be soon. I am also not at all convinced that it is necessary.
I am solving the elastic force balance equations. I compared the results that I got from my fipy code with those from ABAQUS (a finite element software). The results are very similar except they are offset somehow (depending on the problem). If I use the "getFaceGrad()" instead of "getGrad()" when calculating the strain from the displacement, I get closer results.
Is this offset affected by changes to the mesh? Has your solution converged?
Maybe the higher-order interpolation in finite element method makes the differences.
It's possible, but first we would need to know that your solutions are well converged and that they are posed correctly. What equations are you solving? How are you representing them in FiPy?
