Leakage pipeline with lattice Boltzmann

I want to do a simulation of a pipe with one or more leaks, and i not sure about how to implement the boundary conditions in the outlets of the flows, specifically in relation to implementing the leak itself (for example in the middle of the domain of the simulation) and the output of the flow at the end of the pipe.
Has anyone done something similar or has any idea how to implement it?