  This work presents the application of the non-local multicontinuum method
(NLMC) for the Darcy-Forchheimer model in fractured media. The mathematical
model describes a nonlinear flow in fractured porous media with a high inertial
effect and flow speed. The space approximation is constructed on the
sufficiently fine grid using a finite volume method (FVM) with an embedded
fracture model (EFM) to approximate lower dimensional fractures. A non-local
model reduction approach is presented based on localization and constraint
energy minimization. The multiscale basis functions are constructed in
oversampled local domains to consider the flow effects from neighboring local
domains. Numerical results are presented for a two-dimensional formulation with
two test cases of heterogeneity. The influence of model nonlinearity on the
multiscale method accuracy is investigated. The numerical results show that the
non-local multicontinuum method provides highly accurate results for
Darcy-Forchheimer flow in fractured media.
