boundary conditions effect on Matrix coefficients
Hi all,
I am trying to understand how openfoam generates the matrices for the differential operators. I have programmed a simple example for the laplacian operator: Code:
// Define a matrix that contains the coefficients of the Laplacian * Matrix Diagonal terms * Matrix out of diagonal terms * Indices of the out of diagonal terms * source terms (right hand side of the linear system) (zero in my example) But the effect of the boundary conditions is not yet included in the matrix and source term that I get. My question is: How can I get the contribution (coefficients and indices) of the Boundary Conditions to the matrix coefficients (diag and out-of-diag) and to the source term? Thanks in advance for your help. --CME |
same question here. any news on this?
|
Quote:
Not much on this unfortunaltely: I have use these two extra lines to get the information from the boundary conditions: Code:
Info <<"\n TEqnLaplaciano.boundaryCoeffs() \n " << TEqnLaplaciano.boundaryCoeffs() << endl; --Carlos |
Thanks Carlos,
this explains why I am getting incorrect results for the eigenvalues/functions of the laplacian when the mesh is non cartesian. does this mean that the additional non-orthogonal contribution are in TEqn.source()? If this is the case I could just try to put them back in the matrix. or how about using an uncorrected or fourth scheme? I was able to add the boundary condition contribution to the matrix with this code Code:
forAll(TEqn.internalCoeffs(), patchI) Quote:
|
All times are GMT -4. The time now is 12:23. |