Zero solution by solving the Laplace problem on part of the mes

yes, you need to fix the boundary conditions to get an exact match.
See “natural boundary conditions” to get them for the road-only example.
I took some of this and your mesh and ended up with this code
Transport-mod.edp (3.9 KB)

https://community.freefem.org/t/normal-boundary-condition/341