The foundation of numerical codes used in engineering analyses of two-phase flows is the two-fluid model (TFM). However, the TFM codes use artificial regularization to remove the high frequency ill-posed instability in the numerical solution. This work demonstrates that incorporating variational inertial coupling terms to the numerical two-fluid model code, makes it more complete, objective and well-posed without the need for any regularization. The variational TFM is implemented in an industrial CFD code and the two-fluid Burgers problem is used to verify the numerical TFM against analytic solutions.