We present numerical simulations of a directional coupler based on three-dimensional waveguides made of a nematic liquid crystal, acting as the waveguide core, infiltrated in polydimethysiloxane channels. Modeling is based on the combination of minimization of Oseen-Frank energy of the liquid crystal molecules with a beam propagation algorithm. Design of the coupler waveguides is optimized to minimize coupling lengths and maximise efficiencies. Such components can be made at low cost on flexible plastic substrates and can be also integrated with optofluidic devices for biomedical applications.