<div dir="ltr"><div>Hey Jonas,</div><div><br></div><div>welcome to PDELab! I wished I had a better welcoming gift to you...</div><div><br></div><div>I am assuming that the scheme you are trying to implement is something like Figure 1b from</div><div>"Nonlinear two-point flux approximation for modelingfull-tensor effects in subsurface flow simulations"</div><div>by Chen et al. - correct me if I am wrong.</div><div><br></div><div>We get this question and variations of it a lot, and unfortunately the most adequate answer is:</div><div>The assumption that a facet integral only depends on the two neighboring cells</div><div>is built quite deeply into the PDELab assembler. That being said, there is no clean solution to</div><div>your problem in PDELab.</div><div><br></div>I know that DuMux (<a href="http://www.dumux.org">www.dumux.org</a>, a downstream project) has treated that problem before.<div>I am not sure they are relying on the PDELab assembler in their implementation, but it is worth</div><div>asking them.</div><div><br></div><div>Of course, things can always be hacked. Most probably one can do something with precomputation</div><div>of all flux values at all intersection mid points, but I am not sure about the details (e.g. I have no</div><div>clues about how jacobians for such a system could be computed).<br></div><div><br></div><div>Best,</div><div>Dominic<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Jul 6, 2020 at 3:30 PM Jonas Beddrich <<a href="mailto:jonas.beddrich@gmail.com">jonas.beddrich@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi everyone, </div><div><br></div><div>I'm looking for a method to access the values of cells in the alpha_skeleton function of my local operator in addition to the values of the inside and outside cell of the intersection. <br></div><div><br></div><div>The nonlinear Two Point Flux Approximation I'm currently working on requires information about 4 cells to compute the flux through each intersection. So in addition to the inside and outside cell two additional cells m and k are required. The cells m and k are dependent on the intersection, the inside and outside cell as well as problem parameters. Right now I can determine the index of those cells and then get the corresponding elements by iterating over the grid. <br></div><div><br></div><div>Sadly, I haven't found an efficient way to access the values of these cells, since alpha_skeleton only gets the local vectors and local function spaces of the inside and outside cell. But I guess I'm not the only person encountering this problem. Can anyone give me a hint on how to solve this? <br></div><div><br></div><div>Best regards, <br></div><div>Jonas <br></div></div>
_______________________________________________<br>
dune-pdelab mailing list<br>
<a href="mailto:dune-pdelab@lists.dune-project.org" target="_blank">dune-pdelab@lists.dune-project.org</a><br>
<a href="https://lists.dune-project.org/mailman/listinfo/dune-pdelab" rel="noreferrer" target="_blank">https://lists.dune-project.org/mailman/listinfo/dune-pdelab</a></blockquote></div>