[dune-pdelab] Accessing solution vector for MPFA

[Shubhangi Gupta]

Hi Christoph,

Thanks a lot for your answer.

I have another quick question:
Do you think I can use a non linear two-point flux approximation method in DUNE PDELab? Or will this also create a problem due to neighbouring pressure values?

Thanks again,
Best regards,
Shubhangi


----- Ursprüngliche Mail -----
Von: "Christoph Grüninger" <christoph.grueninger at iws.uni-stuttgart.de>
An: dune-pdelab at dune-project.org
Gesendet: Dienstag, 20. September 2016 15:42:43
Betreff: Re: [dune-pdelab] Accessing solution vector for MPFA

Hi Shubhangi,
does your MPFA method fulfil "FE and FV schemes requiring at most
face-neighbors" from [1]?
My experience is, that complicated FV schemes exceed this restriction.
You can try to program around this inside the local operator (like your
pointer to the global solution, that should be never needed), but the
resulting stencil generated by PDELab will remain too small.

I still try to answer your question:
The local solution is altered to calculate the finite difference.
Usually, the local solution is changed by 1e-6, 1e-8 or similar.

Bye
Christoph

[1] https://www.dune-project.org/modules/dune-pdelab/

-- 
The method has been developed for use on a high-speed electronic
computer and would be impractical for hand-solution purposes.
                                           [Harlow & Welch 1965]


_______________________________________________
dune-pdelab mailing list
dune-pdelab at dune-project.org
http://lists.dune-project.org/mailman/listinfo/dune-pdelab