Computational fluctuating hydrodynamics aims at understanding the impact of thermal fluctuations on fluid motions small scales through numerical exploration. These are modeled as stochastic flux terms and incorporated into classical Navier-Stokes equations, which need to be solved numerically. In this paper, we present a novel projection-based method for solving incompressible (FHD) equations. ...