  We report a hybrid numerical method for the solution of the model H
fluctuating hydrodynamic equations for binary mixtures. The momentum
conservation equations with Landau-Lifshitz stresses are solved using the
fluctuating lattice Boltzmann equation while the order parameter conservation
equation with Langevin fluxes are solved using the stochastic method of lines.
Two methods, based on finite difference and finite volume, are proposed for
spatial discretisation of the order parameter equation. Special care is taken
to ensure that the fluctuation-dissipation theorem is maintained at the lattice
level in both cases. The methods are benchmarked by comparing static and
dynamic correlations and excellent agreement is found between analytical and
numerical results. The Galilean invariance of the model is tested and found to
be satisfactory. Thermally induced capillary fluctuations of the interface are
captured accurately, indicating that the model can be used to study nonlinear
fluctuations.
