[petsc-users] issues with VecSetValues in petsc 3.19

[Matthew Knepley]

On Mon, Apr 17, 2023 at 6:37 AM Edoardo alinovi <edoardo.alinovi at gmail.com>
wrote:

> Sure thing, the solver I am working on is this one:
> https://gitlab.com/alie89/flubio-code-fvm.
>
> It is a 3D, collocated, unstructured, finite volume solver for
> incompressibility NS. I can run steady, unsteady and I can use SIMPLE, PISO
> and Factional step method (both explicit and fully implicit momentum). I
> can also solve for turbulence (k-omega, BSL, SST, Spalart-Allmars, LES). I
> have also implemented some kind of Immersed boundary (2D/3D) that I need to
> resume at some point.
>
> Hot topic of the moment, I am developing a fully coupled pressure based
> solver using field-split. What I have now is working ok, I have validated
> it on a lot of 2D problems and going on with 3D right now.  If all the
> tests are passed, I'll focus on tuning the field splitting which looks to
> be a quite interesting topic!
>

I think a very good discussion of the issues from the point of view of FEM
is

   https://arxiv.org/abs/1810.03315

There should be a similar analysis from the FVM side, although it might not
be possible to
find a pressure discretization compatible with the FVM velocity for this
purpose.

  Thanks,

     Matt


> Flubio is a project I have been carrying on since PhD days. The
> implementation is 99% on my shoulders, despite the fact I am
> collaborating with some people around. I am coding evenings and
> weekends/free time, it gives me a lot of satisfaction and also a lot of
> insights!
>

-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener

https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20230417/442e0447/attachment-0001.html>