Equilibrium structures, solvation forces, and conformational dynamics of thin confined films of n-hexadecane and squalane are investigated using a new grand canonical ensemble molecular dynamics method for simulations of confined liquids. The method combines constant pressure simulations with a computational cell containing solid surfaces and both bulk and confined liquid regions in equilibrium with each other. For both molecular liquids layered density oscillations in the confined films are found for various widths of the confining gap. The solvation force oscillations as a function of the gap width for the straight chain n-hexadecane liquid are more pronounced exhibiting attractive and repulsive regions, while for the branched alkane the solvation forces are mostly repulsive, with the development of shallow local attractive regions for small values of the gap width. Furthermore, the nature of the transitions between well-formed layered configurations is different in the two systems, with the n-hexadecane film exhibiting solid-like characteristics portrayed by step-like variations in the number of confined segments occurring in response to a small decrease in the gap width, starting from well-layered states of the film. On the other hand the behavior of the squalane film is liquid-like, exhibiting a monotonic continuous decrease in the number of confined segments as the gap width is decreased. These characteristics are correlated with structural properties of the confined films which, for n-hexadecane, exhibit enhanced layered ordering and in-plane ordered molecular arrangements, as well as with the relatively high tendency for interlayer molecular interdigitation in the squalane films. Reduced conformational (trans-guache) transition rates in the confined films, compared to their bulk values, are found, and their oscillatory dependence on the degree of confinement is analyzed, showing smaller transition rates for the well-formed layered states of the films.