Convective heat transfer in boundary

Hi, I am trying to calculate the heat transfer from a solid to a liquid. I am simplifying the boundary by assuming a Neumann bc and getting the convective heat transfer value from an experimental reference.
So my problem looks like this:

−k dT/dx n = h(Ts − Tl)

Where k is the heat conduction of the solid, h is the conductive heat transfer coefficient.

I am not sure how to implement this condition as my heat transfer will depend on the values of T on both sides of the interface (Ts = Tsolid)_at the interface and (Tl = Tliquid)_at the interface.

Any suggestions?

Look up for “Robin boundary conditions” in https://fenicsproject.org/pub/tutorial/sphinx1/._ftut1005.html.