<div dir="ltr"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><span style="font-size:12.8px">There are more efficient split methods than SIMPLE, particularly for</span><br style="font-size:12.8px"><span style="font-size:12.8px">steady state or large time steps, but they are still low Mach solvers.</span></blockquote><div><br></div><div>I only mention SIMPLE as a canonical segregated/pressure-based method; segregated methods aren't really my area.</div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><span style="font-size:12.8px">Are you doing steady state or transient solves?  RANS? </span></blockquote><div><br></div><div>Varies -- I wrote a high order unstructured parallel DG code for coupled compressible NS that does fantastic with explicit time DNS, and I'm trying to get any kind of good implicit time or steady state solutions to attack the majority of problems where explicit is useless. Basically anything in the AIAA high order workshops is of interest to me. RANS, LES, and steady state are all of interest; I think steady state solutions in this area are probably doomed to use a pseudo-transient method to get there anyway so it's the same problem. I don't know if anyone in high order cfd has actually gotten a purely steady state formulation to converge on something nontrivial.</div><div><br></div><div><br style="font-size:12.8px"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><span style="font-size:12.8px"> There aren't </span><span style="font-size:12.8px">good AMG solvers available for this kind of system; not for lack of<br></span><span style="font-size:12.8px">trying but it remains a worthy area of research.</span></blockquote></div><div><br></div><div>Hence me using petsc as a solver library for research =), it's the best option for trying various kinds of algebraic solvers (seriously you devs have done fine work). Relatively generic AMG solvers have had massive success in the same problems for FV discretizations, so it seems hopeful they could extend to DG.</div><div><br></div><div>I recall you working in this area, so any suggestions welcome.</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Nov 8, 2017 at 11:12 PM, Jed Brown <span dir="ltr"><<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">Mark Lohry <<a href="mailto:mlohry@gmail.com">mlohry@gmail.com</a>> writes:<br>
<br>
>><br>
>><br>
>>>> You want to use the Schur complement PCs for NS. We have some support<br>
>>>> for PC where you give us a mass matrix.<br>
>>><br>
>>><br>
>>> I'll look into it, but I'm not aware of Schur complement being used for<br>
>>> compressible/coupled NS, only for incompressible/segregated.<br>
>>><br>
>><br>
>><br>
>> Really? OK, you've looked into it more than me.<br>
>><br>
><br>
> Correct me if I'm wrong, I haven't looked into it much at all. I just<br>
> haven't encountered Schur complement in CFD outside of segregated<br>
> SIMPLE-type algorithms (like in the ex70/PCFIELDSPLIT example.) Those<br>
> segregated pressure-correction methods tend to get progressively worse as<br>
> you become less elliptic / non-zero Mach.<br>
<br>
</span>There are more efficient split methods than SIMPLE, particularly for<br>
steady state or large time steps, but they are still low Mach solvers.<br>
Are you doing steady state or transient solves?  RANS?  There aren't<br>
good AMG solvers available for this kind of system; not for lack of<br>
trying but it remains a worthy area of research.<br>
</blockquote></div><br></div>