Implementation of Robin BC by groovyBC : difference wrt. theoretic solution
1 Attachment(s)
Hello everyone, we consider the following stationary heat equation:
Laplacian(T)=0 Code:
T(x) a bar with length L=0.4 m, mesh Delta x=1 mm At x=0 : T(0)=T0=273 K At x=L : kT'(x)=h*(T(L)Tinf) with k=237 W/mK, h=100 W/m^2K, Tinf=873 K We have implemented this BC using swak4Foam, in 0/T for the patch outlet: Code:
type groovyBC; However, the OF solution differs a little from the theoretic solution : Attachment 12409 
Quote:

After we changed the gradient scheme to leastSquares, everything is fine !! However we still don't know why...

All times are GMT 4. The time now is 23:10. 