Monotonicity and strict concavity

Hi all,

I’ve got a partial differential equation that I’m trying to approximate with a physics inspired neural network approach. The issue is that it doesn’t quite come out right. I know that the solution should be monotonic and strictly concave, and I think I have to impose those conditions, and I’m wondering if anyone has any ideas about how to go about it.

Many thanks,
Federico