The development of a complex functional multicellular organism from a single cell involves tightly regulated and coordinated cell behaviours coupled through short- and long-range biochemical and mechanical signals. To truly comprehend this complexity, alongside experimental approaches we need mathematical and computational models, which can link observations to mechanisms in a quantitative, predictive, and experimentally verifiable way. In this talk I will describe our efforts to model aspects of embryonic development, focusing in particular on the planar polarised behaviours of cells in epithelial tissues, and discuss the mathematical and computational challenges associated with this work. I will also highlight some of our work to improve the reproducibility and re-use of such models through the ongoing development of Chaste (github.com/chaste), an open-source C++ library for multiscale modelling of biological tissues and cell populations.