Fluid flows in porous media are ubiquitous in the shallow Earth environments. The resistivity of the porous medium to the flow of either a single or multiphase fluid does not only effect the fluid discharge trough the porous medium, but also the transport of additional dissolved components or heat. In many instances, chemicals or heat drive reactions at the fluid-solid interface. These reactions can affect the porous media geometry at the pore-scale leading to a local change in permeability and phases saturation. Reactions can therefore have a drastic effect on discharge of mass and advected quantities. In this thesis we aim to better understand and quantify the complex coupling between the different physical processes that take place at the pore-scale for both single and multiphase fluids in isothermal and thermal reactive environments. We do that by performing pore-scale numerical simulations for both single and multiphase fluids in synthetic porous media samples. The numerical algorithm we choose is the lattice Boltzmann method.