<div dir="ltr">This has been a great discussion to follow. Regarding<div><br></div><div>> when time stepping, you have enough mass matrix that cheaper preconditioners are good enough</div><div><br></div><div>I'm curious what some algebraic recommendations might be for high Re in transients. I've found one-level DD to be ineffective when applied monolithically or to the momentum block of a split, as it scales with the mesh size. For high Re boomeramg is ineffective perhaps until <a href="https://gitlab.com/petsc/petsc/-/issues/1362">https://gitlab.com/petsc/petsc/-/issues/1362</a> is resolved. I should try fiddling around again with Pierre's work in HPDDM, but curious if there are other PETSc PC recs, or if I need to overcome my inertia/laziness and move beyond command line options.</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, May 8, 2023 at 6:46 AM Jed Brown <<a href="mailto:jed@jedbrown.org">jed@jedbrown.org</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">Sebastian Blauth <<a href="mailto:sebastian.blauth@itwm.fraunhofer.de" target="_blank">sebastian.blauth@itwm.fraunhofer.de</a>> writes:<br>
<br>
> Hello everyone,<br>
><br>
> I wanted to briefly follow up on my question (see my last reply).<br>
> Does anyone know / have an idea why the LSC preconditioner in PETSc does <br>
> not seem to scale well with the problem size (the outer fgmres solver I <br>
> am using nearly scale nearly linearly with the problem size in my example).<br>
<br>
The implementation was tested on heterogeneous Stokes problems from geodynamics, and perhaps not on NS (or not with the discretization you're using).<br>
<br>
<a href="https://doi.org/10.1016/j.pepi.2008.07.036" rel="noreferrer" target="_blank">https://doi.org/10.1016/j.pepi.2008.07.036</a><br>
<br>
There is a comment about not having plumbing to provide a mass matrix. A few lines earlier there is code using PetscObjectQuery, and that same pattern could be applied for the mass matrix. If you're on a roughly uniform mesh, including the mass scaling will probably have little effect, but it could have a big impact in the presence of highly anistropic elements or a broad range of scales.<br>
<br>
I don't think LSC has gotten a lot of use since everyone I know who tried it has been sort of disappointed relative to other methods (e.g., inverse viscosity scaled mass matrix for heterogeneous Stokes, PCD for moderate Re Navier-Stokes). Of course there are no steady solutions to high Re so you either have a turbulence model or are time stepping. I'm not aware of work with LSC with turbulence models, and when time stepping, you have enough mass matrix that cheaper preconditioners are good enough. That said, it would be a great contribution to support this scaling.<br>
<br>
> I have also already tried using -ksp_diagonal_scale but the results are <br>
> identical.<br>
<br>
That's expected, and likely to mess up some MG implementations so I wouldn't recommend it.<br>
</blockquote></div>