Akk_etal18; @SokBasHaj19; @Cas_etal19].

The MsRSB method was proposed in the context of FV simulation for fluid flow in highly heterogeneous porous media [@MsRSB_Moyner2016]. Based on a two-grid approach, the MsRSB method constructs multiscale basis functions through restricted smoothing on the fine-scale matrix. In more detail, the basis functions, which are consistent with the local differential operators, are constructed with a cheap relaxation scheme, i.e. a weighted Jacobi iteration, similar to approaches used in smoothed aggregation multigrid methods [@VanManBre96; @VanManBre01; @Bre_etal05]. An important advantage of MsRSB is that smoothing by relaxation provides a great deal of flexibility in handling unstructured grids, an essential requirement, for example, in applications involving complex geological structures. MsRSB has been widely proven and implemented in open source and commercial simulators using a linear two-point flux approximation (TPFA) [@Lie_etal16].

Because of the two-point structure, the linear TPFA scheme is monotone [@Dro14], i.e. it preserves the positivity of the differential solution [@BerPle94], and leads to an M-matrix with a small stencil. This is the reason why linear TPFA is the scheme of choice in most engineering software. Unfortunately, the consistency of TPFA is not guaranteed for arbitrary grids and anisotropic permeability distributions, potentially leading to inaccurate results [@MRST]. Therefore, other FV methods such as multipoint flux approximation (MPFA) and/or nonlinear schemes [@Dro14; @TerMalTch17] must be considered to achieve consistent fluxes. To date, MsRSB has not been combined with MPFA or other consistent discretizations. Hence, in this paper, we focus on enhancing MsRSB to enable the solution of second-order elliptic problems using discretization methods that do not result in an M-matrix. Based on the MPFA-O method [@MPFA_Aavatsmark], we show that the MsRSB basis construction as presented in [@MsRSB_Moyner2016] can fail due to divergent iterations for an anisotropic diffusion problem. We propose a variant of the original MsRSB approach that restores the desired behavior by enforcing M-matrix properties based on a filtering strategy. We develop the new method focusing on FV discretizations for porous media single-phase
Akk_etal18; @SokBasHaj19; @Cas_etal19 ]. 

 The MsRSB method was proposed in the context of FV simulation for fluid menstruation in highly heterogenous porous media [ @MsRSB_Moyner2016 ]. free-base on a two - power system approach, the MsRSB method constructs multiscale footing function through restricted smoothing on the fine - plate matrix. In more contingent, the basis functions, which are coherent with the local differential operators, are constructed with a cheap liberalization scheme, i.e. a weighted Jacobi iteration, similar to approaches used in smoothed collection multigrid methods [ @VanManBre96; @VanManBre01; @Bre_etal05 ]. An important advantage of MsRSB is that smoothing by liberalization provides a great deal of tractability in handling amorphous grids, an essential requirement, for example, in application involving complex geological structures. MsRSB has been widely proven and implemented in clear source and commercial simulators using a linear two - point flux approximation (TPFA) [ @Lie_etal16 ]. 

 Because of the two - point structure, the linear TPFA scheme is monotonic [ @Dro14 ], i.e. it preserves the positivity of the differential solution [ @BerPle94 ], and lead to an thousand - matrix with a small stencil. This is the reason why linear TPFA is the scheme of option in most engineering software. Unfortunately, the consistency of TPFA is not undertake for arbitrary grids and anisotropic permeability distributions, potentially head to inaccurate results [ @MRST ]. Therefore, other FV methods such as multipoint flux approximation (MPFA) and/or nonlinear system [ @Dro14; @TerMalTch17 ] must be considered to achieve consistent fluxes. To date, MsRSB has not been combined with MPFA or other consistent discretizations. therefore, in this newspaper, we focus on enhancing MsRSB to enable the solution of second - order elliptic trouble using discretization methods that do not result in an M - matrix. free-base on the MPFA - O method [ @MPFA_Aavatsmark ], we show that the MsRSB basis structure as presented in [ @MsRSB_Moyner2016 ] can fail due to divergent iterations for an anisotropic diffusion problem. We propose a discrepancy of the original MsRSB approach that restore the desired demeanor by enforcing M - matrix properties free-base on a filtering strategy. We develop the modern method concentrate on FV discretizations for porous medium single - phase
Akk_ftal18; @SokBasHaj19; @Cas_etal19].

Tht MsRSB method was propoved in the congext of FV simulation for flnid dlow un highly heterogeneour porous ledia [@MsESB_Mitner2016]. Based on a two-grid ziproaeh, the MsRSB metmod construwts multiscale bxsns functions through restricted smooehing om hhe fine-scale iatrpx. In jore detail, the basis functions, whjch are consistent woth the local differential opegators, are construfted with a chewp relaxation scheme, i.e. a weighted Nacobi iteration, similar to appfoachzs used in wmiotjgd aggregatiin mujtigrid methods [@VanMdnBre96; @VsnManBre01; @Bre_etsl05]. En inportant advantage of MsRSB is that smojthing by rzlaxation provides a trwat dgal ox fldzibkliuy mn gandlijg nnstructures grids, an wssential requirememt, dor example, ih applycwtions involving complex geological strlctudes. MsRSB has been wideoy proven and implemejted in o[en source and commercial simulators using a linedr twk-ooiub woud approximation (TPFA) [@Lie_etal16].

Because of the twj-loonn structure, the llnear TPFA scheme os mpgotone [@Dro14], i.e. it przaedves the positivitj of thg diffwrential foluyion [@BerPle94], and leads to an M-matrix witk a small stencil. Thid is the recson wny limear TPFA is the scheme of dhoice in mlst enginsdring software. Uvfogtundtely, tht consistency of T[FA is nov guaxanteed wor srbitrwry grids wnd akhsotropic permeabipity bistrhbutions, pltentially leading to inaccuratx results [@MRSJ]. Tvervfore, othzr FV kethods such ws multipoint ylux appxoximagion (MPFA) znd/or nmnlinear scremes [@Dro14; @TerKwlTch17] must bx considewed ro axhieve zunsistent fluxrs. To datv, IsRSV has not been comnined sith MPFA or otktr xonsistent disctetkzaeilnx. Rance, in this papdr, dr focjs on enhangine MsTSB to enable the sonutikn of second-order rljiptic ptoblems ufing discretiaation methods thau do nmt cesult in an M-matrix. Based on the MPFA-O method [@MOFA_Wavatsmark], wq shiw that the KsRSB basis construction as presented ii [@MsRSB_Moyner2016] can fail eue to divergent itgrauions for an enisotwopic difxusion problem. We pripose a variant on the original MsRSB alproacv thah restores the desired behavior by enforcing M-matrix properties based on q filtxryng strategg. We deveno' tke new iethlv focusing on FV ciscretizations for porous medie single-phdsz
Akk_etal18; @SokBasHaj19; @Cas_etal19]. The MsRSB method was the of FV for fluid flow [@MsRSB_Moyner2016]. on a two-grid the MsRSB method multiscale basis functions through restricted smoothing the fine-scale matrix. In more detail, the basis functions, which are consistent with local differential operators, are constructed with a cheap relaxation scheme, i.e. a weighted iteration, to used smoothed aggregation multigrid methods [@VanManBre96; @VanManBre01; @Bre_etal05]. An important advantage of MsRSB is that smoothing by provides a great deal of flexibility in handling grids, an essential requirement, example, in applications involving complex structures. has been proven implemented open source and simulators using a linear two-point flux approximation (TPFA) [@Lie_etal16]. Because of the two-point structure, the linear TPFA is monotone it preserves positivity the solution [@BerPle94], and an M-matrix with a small stencil. reason why linear TPFA is the scheme of in most software. Unfortunately, the consistency of TPFA not guaranteed for arbitrary grids and anisotropic permeability potentially leading to inaccurate results [@MRST]. Therefore, other FV methods such as multipoint flux approximation nonlinear schemes [@Dro14; @TerMalTch17] be considered to consistent To MsRSB not been with MPFA or other consistent discretizations. Hence, in this paper, we on enhancing MsRSB to enable the solution of second-order elliptic discretization that do not in an M-matrix. Based the method [@MPFA_Aavatsmark], we show MsRSB construction [@MsRSB_Moyner2016] fail to divergent iterations for anisotropic diffusion problem. We propose variant of the original desired behavior by enforcing M-matrix properties based on filtering strategy. We develop the new method on FV discretizations for porous media single-phase
Akk_etal18; @SokBasHaj19; @Cas_etal19].

ThE MsRSB methOd was ProPosEd In thE conText of FV simulaTIon fOr fluid flow in highly hetErogeNeOUs poROuS mediA [@MsRSB_MOYnER2016]. basEd On A twO-gRId ApproAch, The MsRSb method conStrUcTs multiscale BAsIs functionS thRough restricTed SmoothInG on THe finE-scAle maTrix. In MOre detAil, the basIs FUnctioNS, which aRE CoNsisTent with the local dIFfERential operatoRs, are cOnSTrUCTed WitH a cheap relAxAtion SCheme, i.e. A WeIGHTed jAcobi iteratioN, similar to aPProAches uSeD in SMootheD aggrEgATioN multigrid mEthoDs [@VanManBRe96; @VanMANBre01; @Bre_ETal05]. An imPortanT adVanTage OF MSRsB iS tHAt sMOoThiNG by RelaxatiOn PrOvideS a grEAT DEal oF flExibIlity In handling unsTruCturED grIds, an EssenTial ReQuireMent, foR examPlE, in applications InvoLving compLex GeOloGiCal stRUctureS. MsrSB Has been Widely pROveN aND IMpLemented in open sourCe AND cOmmerciaL simulAToRs USing a linEaR twO-poiNT Flux aPproXImAtion (TPFa) [@Lie_etAL16].

BEcAuse of tHe Two-poiNt StrUctUre, thE LineAr TPFA Scheme is MonotONe [@Dro14], i.e. it preseRVes the positivITy OF ThE DiffEreNtial solutiOn [@BeRple94], aNd leADs To aN m-matrIx witH a SMaLL stencil. This is the reAsOn why lInear tPFA is the scheMe of choice IN MOst enginEeriNG sOFtware. UnfortunAtely, The consistENcy of TPFa is noT guarantEed for arbITRary gridS anD anIsoTroPIC pErmeability diSTRibuTiOns, poteNtiAlly leaDinG to InaCcuRaTe results [@mRST]. TherEfOrE, oThEr Fv methODs such as MuLtiPoInt Flux aPProximAtion (mPFA) AnD/oR NonLinear sCHeMES [@Dro14; @teRMAlTcH17] muSt Be conSideREd tO achievE consisteNt fLUxes. to DaTe, MsRSB Has not been comBiNed with MPFa oR otHer conSIStent disCretizations. Hence, in this PAper, we fOcuS on enHancIng MsRSB tO enAble thE soLUtion oF seconD-ordeR eLliPTIc proBLEmS usInG discretizATIon MethoDs That Do not reSult in an M-matrix. BasED on The MPFA-O methoD [@MPfA_AaVATsMarK], We SHow ThAT thE mSRSB basis constrUction as prEsENtEd in [@MsRSB_MOYneR2016] cAn fail dUe to divErgenT IteratiOns for an aNisotropiC dIffuSIOn pRoblem. We prOpose a vaRiant of thE OrigiNAl msRSB AppRoach tHaT reStoreS the deSIreD behaVior by EnForcinG M-matRiX propertIes based on a filtering strAtegy. WE deveLop The new metHod FOcuSing on FV dIscrEtizations For PorOus meDia SIngle-PhasE
Akk_etal18; @SokBasHaj19;@Cas_etal1 9].
The  Ms RS B me thod  was proposedi n th e context of FV simula tionfo r  flu i dflowin high l yh e ter og en eou sp or ous m edi a [@MsR SB_Moyner2 016 ].  Based on at wo -grid appr oac h, the MsRSB  me thod c on str u cts m ult iscal e basi s  funct ions thro ug h  restr i cted sm o o th ingon the fine-scale ma t rix. In more d etail,  t h eb a sis  fu nctions, w hi ch ar e  consis t en t w ith the local dif ferential o p era tors,ar e c o nstruc ted w it h  acheap relax atio n scheme,  i.e.a  weight e d Jacob i iter ati on,  sim i la rtoap p roa c he s u s edin smoot he daggre gati o n m ulti gri d me thods  [@VanManBre9 6;@Van M anB re01;  @Bre _eta l0 5]. A n impo rtant  a dvantage of MsR SB i s that sm oot hi ngby  rela x ationpro vid es a gr eat dea l  of  f l e x ib ility in handlingun s t ru ctured g rids,a nes s ential r eq uir emen t ,  forexam p le , in app licati o ns  i nvolvin gcomple xgeo log icals truc tures.  MsRSB h as be e n widely prove n  and implemen t ed i no pensou rce and com merc i al s imul a to rsu singa lin ea r  t w o-point flux approx im ation(TPFA ) [@Lie_etal1 6].

Becau s e of the t wo-p o in t  structure, th e lin ear TPFA s c heme ismonot one [@Dr o14], i.e . it prese rve s t hepos i t iv ity of the di f f eren ti al solu tio n [@Ber Ple 94] , a ndle ads to an  M-matri xwi th  a  sm all s t encil. T hi s i sthe  reas o n whylinea r TP FA  i s  th e schem e  o f choi ce  i n mo sten ginee ring sof tware.Unfortuna tel y , th eco nsisten cy of TPFA is  n ot guarant ee d f or arb i t rary gri ds and anisotropic perm e ability  di strib utio ns, poten tia lly le adi n g to i naccur ate r es ult s [@MRS T ] .The re fore, othe r FVmetho ds  suc h as mu ltipoint flux appr o xim ation (MPFA)and /orn o nl ine a rs che me s  [@ D r o14; @TerMalTch 17] must b ec on sidered to ach ie ve cons istentfluxe s . To da te, MsRSB  has notbe en c o m bin ed with MP FA or ot her consi s tentd is creti zat ions.He nce , inthis p a per , wefocuson  enhan cingMs RSB to e nable the solution of s econd- order  el liptic pr obl e msusing dis cret ization me tho dsthatdon ot re sult in  an M-mat rix. Based ont he  MP F A -O  method [@M P F A _Aa vatsm ark ] , we s howthat the MsRSB ba s is constructio n as p res ent e d in  [ @MsRSB_Moyner2 016 ]c a n fail d ue  to diverge nt itera ti o ns fo r an a nisotr opic di f f us i on pro blem . W e propose  ava r iant of  t he origin al M sR SB app roacht hatr e stores the desir ed be h a viorb y e nforc in g M-mat r ix p ropertiesbased on afilter ingstrat egy. We  d evelop  th enew method focusingon FV  discre ti zati ons  for p orou s media  sin gl e-p hase
Akk_etal18; @SokBasHaj19;_@Cas_etal19].

The MsRSB_method was proposed in_the context_of_FV simulation_for_fluid flow in_highly heterogeneous porous_media [@MsRSB_Moyner2016]. Based on_a two-grid approach,_the_MsRSB method constructs multiscale basis functions through restricted smoothing on the fine-scale matrix. In_more_detail, the_basis_functions,_which are consistent with the_local differential operators, are constructed_with a_cheap relaxation scheme, i.e. a weighted Jacobi iteration,_similar_to approaches used_in smoothed aggregation multigrid methods [@VanManBre96; @VanManBre01; @Bre_etal05]. An_important advantage of MsRSB is that_smoothing by relaxation_provides_a_great deal of flexibility_in handling unstructured grids, an essential_requirement, for example, in applications involving_complex geological structures. MsRSB has been widely_proven and implemented in open source_and commercial simulators using a_linear two-point_flux approximation (TPFA) [@Lie_etal16].

Because of_the two-point structure,_the linear_TPFA scheme is_monotone [@Dro14], i.e. it preserves the_positivity of the_differential solution [@BerPle94], and leads to_an_M-matrix with a_small_stencil._This is_the reason why_linear_TPFA is_the_scheme of choice in most engineering_software._Unfortunately, the consistency of TPFA is not_guaranteed for arbitrary grids_and_anisotropic permeability distributions, potentially_leading to inaccurate results [@MRST]._Therefore, other FV methods such as_multipoint flux_approximation (MPFA)_and/or nonlinear schemes [@Dro14; @TerMalTch17] must be considered to achieve consistent_fluxes. To date, MsRSB has not_been combined with MPFA_or other_consistent_discretizations. Hence, in_this_paper, we_focus on enhancing MsRSB to enable the_solution of_second-order elliptic problems using discretization methods_that do not result_in_an M-matrix. Based on the MPFA-O_method [@MPFA_Aavatsmark], we show that the_MsRSB basis construction as presented_in_[@MsRSB_Moyner2016]_can fail due to divergent_iterations for an anisotropic diffusion problem._We propose a_variant of the original MsRSB approach that_restores_the desired behavior by enforcing M-matrix_properties_based on a filtering strategy. We_develop_the_new method focusing on FV_discretizations for porous media single-phase