Craniosynostosis is a condition defined by premature closure of cranial vault sutures, which is associated with abnormalities of the brain and skull. Many causal relationships between discovered mutations and premature suture closure have been proposed but an understanding of the precise mechanisms remains elusive. This article describes a computational framework of biological processes underlying cranial growth that will enable a hypothesis driven investigation of craniosynostosis phenotypes using reaction-diffusion-advection methods and the finite element method. Primary centers of ossification in cranial vault are found using activator-substrate model that represents the behavior of key molecules for bone formation. Biomechanical effects due to the interaction between growing bone and soft tissue is investigated to elucidate the mechanism of growth of cranial vault.