We present an efficient scheme within the phase field framework for imposing dynamic contact angle boundary conditions for wall-bounded flows of two immiscible incompressible fluids with large density ratios. First, we develop an algorithm for imposing the dynamic contact angle boundary conditions to the Cahn–Hilliard equation. Our algorithm consists of two components: (i) we ignore the boundar...