|March 26, 2010, 09:52||
Does FOAM converge to exact solution in Laplace equation?
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 396Rep Power: 12
Recently Nisi and me have been working in understanding the nonOrthogonal correction method implemented in FOAM. As have been explained in another threads Over-Relaxed approach is used to correct the solution iteratively for non-orthogonal meshes. We carried out several runnings with the same problem but incresing non-orthogonality. The mesh was simple, a bar with only two triangular (wedge) elements, and different large-height ratio. Since the equation was the Laplace one and temperatures fixed at both extremes, analytical solution is a line and we could compare it with the numerical data. We've found that numerical solution converges but not to the analytical solution, and this phenomenon is so big as the non-orthogonality is. Nisi have found the following errors for different mesh ratios (for 9 nonOrthogonalCorrections):
Mesh ratio %e
2:1 0.1135 %
5:1 0.53375 %
20:1 1.2209 %
Is that interpretation correct? If it is correct is the approach based in any kind of compromise solution or what is the idea behind it (remember that even solution is fairly good, this approach would violate one the cornerstones of numerical analysis --> convergence to analytic solution) ?
Santiago MÁRQUEZ DAMIÁN, Ph.D.
Research Center for Computational Mechanics (CIMEC) - CONICET/FICH-UNL
T.E.: 54-342-4511594 Ext. 1005
Güemes 3450 - (3000) Santa Fe
Santa Fe - Argentina
|Thread||Thread Starter||Forum||Replies||Last Post|
|StarToFoam error||Kart||OpenFOAM Meshing & Mesh Conversion||1||February 4, 2010 05:38|
|Calculation of the Governing Equations||Mihail||CFX||5||July 25, 2008 18:29|
|Exact 2D NS solution for benchmarking?||Tony||Main CFD Forum||4||July 31, 2001 15:50|
|CFL Condition||Matt Umbel||Main CFD Forum||14||January 12, 2001 15:34|
|Analytical solution to the 3d Laplace equation||Chrys Correa||Main CFD Forum||4||September 11, 2000 19:50|