This paper presents a computational framework for modelling the biomechanics of human facial expressions. A detailed high-order (Cubic-Hermite) finite element model of the human head was constructed using anatomical data segmented from magnetic resonance images. The model includes a superficial soft-tissue continuum consisting of skin, the subcutaneous layer and the superficial Musculo-Aponeurotic system. Embedded within this continuum mesh, are 20 pairs of facial muscles which drive facial expressions. These muscles were treated as transversely-isotropic and their anatomical geometries and fibre orientations were accurately depicted. In order to capture the relative composition of muscles and fat, material heterogeneity was also introduced into the model. Complex contact interactions between the lips, eyelids, and between superficial soft tissue continuum and deep rigid skeletal bones were also computed. In addition, this paper investigates the impact of incorporating material heterogeneity and contact interactions, which are often neglected in similar studies. Four facial expressions were simulated using the developed model and the results were compared with surface data obtained from a 3D structured-light scanner. Predicted expressions showed good agreement with the experimental data.

Although many blood flow models have been constructed for cerebral arterial trees, few models have been reported for their venous counterparts. In this paper, we present a computational model for an anatomically accurate cerebral venous tree which was created from a computed tomography angiography (CTA) image. The topology of the tree containing 42 veins was constructed with 1-D cubic-Hermite finite element mesh. The model was formulated using the reduced Navier-Stokes equations together with an empirical constitutive equation for the vessel wall which takes both distended and compressed states of the wall into account. A robust bifurcation model was also incorporated into the model to evaluate flow across branches. Furthermore, a set of hierarchal inflow pressure boundary conditions were prescribed to close the system of equations. Some assumptions were made to simplify the numerical treatment, e.g., the external pressure was considered as uniform across the venous tree, and a vein was either distended or partially collapsed but not both. Using such a scheme we were able to evaluate the blood flow over several cardiac cycles for the large venous tree. The predicted results from the model were compared with ultrasonic measurements acquired at several sites of the venous tree and agreements have been reached either qualitatively (flow waveform shape) or quantitatively (flow velocity magnitude). We then discuss the significance of this venous model, its potential applications, and also present numerical experiments pertinent to limitations of the proposed model.