This paper describes the extension of a recently developed numerical solver for the Landau-Lifshitz Navier-Stokes (LLNS) equations to binary mixtures in three dimensions. The LLNS equations incorporate thermal fluctuations into macroscopic hydrodynamics by using white-noise fluxes. These stochastic PDEs are more complicated in three dimensions due to the tensorial form of the correlations for the stochastic fluxes and in mixtures due to couplings of energy and concentration fluxes (e.g., Soret...

Under the key assumption of finite $\rho$-variation, $\rho \in [1,2)$, of the covariance of the underlying Gaussian process, sharp a.s. convergence rates for approximations of Gaussian rough paths are established. When applied to Brownian resp. fractional Brownian motion (fBM), $\rho =1$ resp. $\rho =1/\left(2H\right)$, we recover and extend the respective results of (Trans. Amer. Math. Soc.361 (2009) 2689–2718) and (Ann. Inst. Henri Poincasé Probab. Stat.48(2012) 518–550). In particular, we establish an a.s. rate $k^{-(1/\rho -1/2-\epsilon )}$, any $\epsilon gt;0$, for Wong–Zakai and Milstein-type...