inclusionsMultiscale.m modified

[Junqiang Wang]

Dear MRST Community,

I tried to use the MSFV method to run inclusionsMultiscale.m instead of MsRSB. However I met the error called " assert(isfield(CG, 'dual'), 'The MsFV method requires a dual grid!');". Then I add the code and met another error as follows,

The completed script file is  attached. How to modify the example script if I want to use the msfv method. 

Looking forward to your reply.

Best regards.

J.WANG

[Olav Møyner]

Hi,

This example is specifically made to demonstrate one of the known issues with the MsFV method: Setting up a valid dual grid is difficult even for a Cartesian grid when there are many inactive cells. So, I'm afraid there is no quick fix to setting up a valid dual grid to use the MsFV method here unless you want to draw it manually. You can read the paper that this example is taken from for more context:

https://doi.org/10.1016/j.jcp.2015.10.010

My co-author Knut-Andreas has a preprint archived here: https://folk.ntnu.no/andreas/papers/msrsb-incomp.pdf

Best regards,

Olav

---

From: sinte...@googlegroups.com <sinte...@googlegroups.com> on behalf of Junqiang Wang <junqiang...@gmail.com>
Sent: Thursday, April 11, 2024 4:59
To: MRST-users: The Matlab Reservoir Simulation Toolbox User Group <sinte...@googlegroups.com>
Subject: [MRST Users] inclusionsMultiscale.m modified

 

--
You received this message because you are subscribed to the Google Groups "MRST-users: The Matlab Reservoir Simulation Toolbox User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email to sintef-mrst...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/sintef-mrst/c75fc560-4323-4d5e-819b-e6e19e3eaf37n%40googlegroups.com.

[Junqiang Wang]

Hello，

I read https://doi.org/10.1016/j.jcp.2015.10.010 and have some questions for the mrst community. Attached is a screenshot of the formula

Is there a specific derivation and reference proof of how the coarsening correction in that formula ensures that the update does not remove the coarser scale conservation.

Looking forward to your reply.

Best regards.

J.WANG

[Olav Møyner]

Hi,

If you move the term y^k into the parenthesis by multiplying with the inverse of R A_c^{-1} P you will see that the right hand side is simply P A_c^{-1} R d^k, i.e. the defect itself, solved with the conservative coarse scale multiscale pressure operator.

The notation may be different, but this trick exists in the earlier iterative MsFV papers to preserve conservation (divergence free velocity fields after reconstruction) for incompressible flow.

Best regards,

Olav

---

From: sinte...@googlegroups.com <sinte...@googlegroups.com> on behalf of Junqiang Wang <junqiang...@gmail.com>

Sent: Thursday, June 13, 2024 11:18

To: MRST-users: The Matlab Reservoir Simulation Toolbox User Group <sinte...@googlegroups.com>

Subject: Re: [MRST Users] inclusionsMultiscale.m modified

 

To view this discussion on the web visit https://groups.google.com/d/msgid/sintef-mrst/4fc76c8b-39d0-4d88-8dc0-7f04dbf73335n%40googlegroups.com.