[ChasingNeutrons]

Hi all,

I have made a solver which apparently works quite well. I calculate the scalar flux of neutrons by tracing rays across the geometry and summing up each ray's contribution to each cell. However, in doing this I've played it pretty loose with boundary conditions; a solver such as this typically has boundary conditions on the rays whereas OF (as far as I'm aware) applies them to the fields themselves, the scalar flux in this case. I have implemented boundary conditions on my rays although in a more roundabout manner which, frankly, I'm fairly happy with for the moment.

However, I gather I still must apply boundary conditions to my scalar flux field. Initially I presumed this wouldn't cause a problem as I thought BCs has to be explicitly applied either via a command or automatically during matrix computations (which I certainly don't do in my solver). Sadly it seems I haven't avoided this problem - rather these boundary conditions are applied and do affect the solution somewhat.

Hence, my question: is there any way to avoid applying BCs in the manner of zeroGradient, uniformValue etc.? Alternately, is there a non-intrusive BC I can easily apply instead?

Cheers for any help.

[kmooney]

Have you tried the 'calculated' type?