A new Matrix Equivalent Discrete Fracture Network model for variable density flow in fractured porous media: Modeling variable density flow (VDF) in fractured reservoirs is a challenging task because of (i) the high contrast of permeability between the fractures and the matrix, (ii) the complexity of the fracture network and (iii) the nonlinearities induced by density variation. The dual-porosity model (DPM) and the discrete fracture-matrix (DFM) model are among the most used models for the simulation of flow in fractured reservoirs. The DPM is efficient as it uses an implicit representation of the fra