Modeling and simulating two-phase flows with moving contact lines pose significant challenges in fluid dynamics due to the problem's multi-scale and multi-physics nature. A continuum model must consider the nanoscale slipness near a contact line. While there are numerous microscopic models available, applying them to macroscopic two-phase flow problems is extremely difficult, especially when dealing with rough or chemically inhomogeneous solid boundaries. In this presentation, I will discuss recent progress in modeling and simulations for complex moving contact line problems. Specifically, I will demonstrate the effectiveness of utilizing the Onsager variational principle as a powerful approximation tool. We derive coarse-grained boundary conditions for moving contact line problems, both with and without contact angle hysteresis. These boundary conditions serve as the foundation for developing efficient numerical methods to solve macroscopic two-phase flow problems.