[lzhou]

Hi,

I am using Finite area method to simulate the scour on river bed.

When calculating the divergence of an "areaVectorField" Qsa (fac::div(Qsa)), I found out that the program suppose the Qsa value at the boundary of the finite area mesh is (0 0 0), even though I set the boundary condition of Qsa as zeroGradient in file 0/Qsa. So the divergence calculated is not correct at the boundary of finite area mesh.

I would like to ask how to solve this problem ? Or is there a way to get the label of the face cell at the boundary of the finite area mesh and correct it ? Thanks !

Lu ZHOU

[tiat]

HI LU,

I am having the same problem as you described and I wonder have you figured out a solution yet? Thanks!

[lzhou]

Hello,

I asked someone else and the answer I got was that the boundary of finite area method did obey the zeroGradient boundary condition. So I think the problem is when we use vsm.mapToSurface(), only the internal field is mapped but not the boundary values. So I set the boundary value by myself because in my case I just need to set the boundary value to zero.

I would like to find someway to modify the boundary values as a function in finite volume mesh like U.correctBoundaryCondition() but I couldn't find one. If you get some better ways to solve this please let me know.

Best regards !

Lu