|February 1, 2015, 01:55||
How can H() return values if itself contains unknowns?
Join Date: Jul 2012
Location: Baton Rouge Louisiana
Posts: 12Rep Power: 6
Hi there, I suddenly get quite confused by the H() function . The meaning of H() is explained here:
and in this excellent post:
From where I quote:
Therefore, the linear system becomes (A + H')U = R, which is the same as AU = R - H'U. The right hand side is simply called H, so H = R - H'U. Therefore, we know have AU = H.
In the above explanation we are trying to solve for U, if we define H = R-H'U, then H itself contains the unknown, this doesn't make sense. On the other hand, if the above operation is understood to be Jacobi iteration (http://en.wikipedia.org/wiki/Jacobi_method). Then this does not make sense either: we don't use Jacobi iteration for anything, and solvers have no error control of it.
In icoFoam, the H() is called after "solve(UEqn == -fvc::grad(p));". There is no return value or assignment for this solve(), then what's the effect of solve() anyway? Does it modify the UEqn itself, then how? If it does modify the matrix, then which U (U at time n-1, U* newly obtained, or a mix of the two) contribute to results of H()?
I would appreciate any explanation from gurus, and people who have the same puzzle are welcomed to post so we can discuss.
|Thread||Thread Starter||Forum||Replies||Last Post|
|Loading Cell Values in Tecplot||P1361||OpenFOAM Post-Processing||3||September 16, 2014 04:02|
|return values from udf||garvitgoel||FLUENT||4||June 16, 2014 13:14|
|cell values vs node values||reversemermaid||FLUENT||0||March 13, 2014 19:06|
|strange node values @ solid/fluid interface - help||JB||FLUENT||2||November 1, 2008 13:04|
|Plotting raw data values||Wilesco||CD-adapco||0||January 5, 2006 06:34|