This question as you might already known has been mostly addressed by means of physical models. In numerical models (CFD) vegetation is most usually represented as rigid elements, but this is used to understand the effect of vegetation on the flow field, which is I believe the inverse question to yours.
After a quick google search I came up with this papers, which seems to be going in the direction you want.
Article Eulerian-Lagrangian flow-vegetation interaction model using ...
They make use of OpenFoam as modelling platform, which would have been one of my first guesses to try to answer your question.