[marc_wattenberg]

Hi together,

I am currently working on understanding the basics of the solver code. I am a litte bit familar with programming, but do not understand e.g.

UEqn().A()
UEqn().H()
mesh.Sf()
etc

For me its clear that some fuction A,H,AU... will act on the momentum equation, but where is the function defined. I already searched in the source guide but can not find A or H. I am sure someone can tell me the right procedure for understanding this

Thanks a lot
Marc

[grtabor]

This is part of the PISO algorithm. The momentum equation can be written as

a_p U_p = H(U) - grad(p)

for a point p. We need to be able to rearange this to solve it. .A() and .H() give access to the relevant components; roughly the diagonal of the matrix and the off-diagonal non-pressure components, respectively.

.Sf() gives the face normal vectors IIRC

Gavin

[marc_wattenberg]

Thanks Gavin,

but this was clear to me. I want to know where to find the definitions of A()... inside the code

Marc

[grtabor]

Dear Marc,

Try rummaging around in OpenFOAM-1.5/src/finiteVolume/fvMatrices

Gavin