https://linen.dev logo
Join Discord
Powered by
# incompressible
  • z

    Zino

    08/27/2025, 8:00 AM
    no clue what the important one was, but start with that, and then one at a time you can put your schemes back to second order
  • z

    Zino

    08/27/2025, 8:00 AM
    until it breaks
  • z

    Zino

    08/27/2025, 8:00 AM
    then you know the issue
  • l

    littleB123

    08/27/2025, 8:01 AM
    THANKS! Ill try these!! Again - Really appreciate your time and effort!!
  • z

    Zino

    08/27/2025, 8:01 AM
    Oh, I also set maxCO to 0.4 instead of 0.9
  • l

    littleB123

    08/27/2025, 8:02 AM
    IVe tried settigns that to 0.1, 0.2 and 0.5 but that didnt solve it alone. If i may ask - why did you remove the underrelaxation of p? ncorrectors -> makes it runs pimpe instead of piso right?
  • z

    Zino

    08/27/2025, 8:03 AM
    There's the file https://cdn.discordapp.com/attachments/726377664213942323/1410172845295927326/Discord.tgz?ex=68b00d4e&is=68aebbce&hm=7fa4247196b22dfca28bfcd44840ca3ce31f617dcacc5a981dfaca906d97dd9b&
  • l

    littleB123

    08/27/2025, 8:04 AM
    Thanks!
  • z

    Zino

    08/27/2025, 8:05 AM
    Oh, uh, force of habit, I normally use Pimple, probably right to leave underrelaxation on p 😆
  • z

    Zino

    08/27/2025, 8:05 AM
    I was speedrunning
  • z

    Zino

    08/27/2025, 8:06 AM
    But no, nCorrectors > 1 is still piso
  • z

    Zino

    08/27/2025, 8:06 AM
    nOuterCorrectors is what makes it pimple
  • z

    Zino

    08/27/2025, 8:08 AM
    Actually wait, if you're using piso you shouldn't have any relaxation factors
  • z

    Zino

    08/27/2025, 8:09 AM
    Yeah you should either get rid of them, or put ".*Final" 1; in there
  • z

    Zino

    08/27/2025, 8:10 AM
    though I guess it doesn't hurt to have while troubleshooting, will keep things a little more stable
  • z

    Zino

    08/27/2025, 8:10 AM
    Anyway, let me know if you're able to figure out your issue 🙂
  • l

    littleB123

    08/27/2025, 8:10 AM
    id look into it!! And ill give you an update!! Also if i manage to get it working with cyclic! 🙂
  • l

    littleB123

    08/27/2025, 8:11 AM
    Im not entirely sure i understand why - is it because they are applied to the outer loop (this case nOuterCorrectors)?
  • z

    Zino

    08/27/2025, 8:12 AM
    Relaxation factors are a thing with SIMPLE
  • z

    Zino

    08/27/2025, 8:12 AM
    For piso, if your courant number is <1 you don't need them
  • l

    littleB123

    08/27/2025, 8:13 AM
    Ahh okay - i didnt know that!
  • l

    littleB123

    08/27/2025, 8:14 AM
    But in piso it could be wrong then - then a timestep will be manipulated my relaxtion right? so even thoguh it makes it more stable then its not a true timestep?
  • z

    Zino

    08/27/2025, 8:15 AM
    Yep
  • l

    littleB123

    08/27/2025, 8:15 AM
    Makes sense thou
  • z

    Zino

    08/27/2025, 8:15 AM
    some of your variables will be simulated in slow motion
  • z

    Zino

    08/27/2025, 8:16 AM
    I make sure I always have:
    Copy code
    c++
    ".*Final"    1;
  • z

    Zino

    08/27/2025, 8:16 AM
    in my relaxation factors, avoids making that mistake
  • l

    littleB123

    08/27/2025, 8:17 AM
    I will - thats a pretty important detail
  • l

    littleB123

    08/27/2025, 10:34 AM
    I think i managed to get it working too. So the changing DDT from backward to euler did it... I used some of the other settings as well afterwards because i believe that is more correct. I still think it is strange that ddt scheme was causing this?
  • z

    Zino

    08/27/2025, 12:33 PM
    Could try letting it run with Euler, then swap to backward when settled