A phase field model for cubic to orthorhombic martensitic transformation (MT) at the nanoscale in a β titanium (Ti) alloy Ti–24Nb–4Zr–8Sn (in wt.%) is investigated by finite element simulation. The approach is based on phase field theory, time-dependent Ginzburg-Landau theory, and mechanical equilibrium equations. Partial differential equations (PDEs) were solved using the commercial software COMSOL Multiphysics. The morphology of the product phase exhibits plate-like or needle-like shapes that reduce the elastic strain energy of the system. The simulation result for random initial order parameters is in agreement with previous experimental observations. The final volume fractions of two different orthorhombic martensitic variants are not dependent on the initial conditions.