PISO Algorithm 관련 질문

Main Forums Solver PISO Algorithm 관련 질문

This topic contains 2 replies, has 2 voices, and was last updated by  geonhong 3 years, 10 months ago.

Viewing 3 posts - 1 through 3 (of 3 total)
  • Author
    Posts
  • #655

    geonhong
    Participant

    안녕하세요,

    MRFInterFoam을 이용하여 자유수면을 포함한 MRF 유동 해석을 하는 중 OpenFOAM의 PISO 알고리즘에 대하여 궁금한 점이 생겼네요.

    PISO Algorithm을 보면 Outer Correction Step이 있는데, Outer Correction이 진행될수록 pressure equation의 initial residual을 줄어들게 됩니다.

    그러다가 Final Iteration에 들어가면 갑자기 initial residual이 증가하네요.

    예를들어,

    fvSolution 파일에서 다음과 같이 설정을 해주었습니다.
    nOuterCorrectors 5;

    그리고 로그를 들여다보면,

    PIMPLE: iteration 1
    GAMG: Solving for p_rgh, Initial residual = 0.2246, Final residual = 0.0008629, No Iterations 3

    PIMPLE: iteration 2
    GAMG: Solving for p_rgh, Initial residual = 0.1631, Final residual = 0.0005971, No Iterations 3

    PIMPLE: iteration 3
    GAMG: Solving for p_rgh, Initial residual = 0.08544, Final residual = 0.0003218, No Iterations 3

    PIMPLE: iteration 4
    GAMG: Solving for p_rgh, Initial residual = 0.04923, Final residual = 0.0002085, No Iterations 3

    PIMPLE: iteration 5
    GAMG: Solving for p_rgh, Initial residual = 0.1326, Final residual = 0.0005350, No Iterations 3

    Outer Correction이 반복될수록 initial residual이 떨어지다가 마지막 corretion step에서는 되려 증가를 하는군요.

    재미있는 것은 nOuterCorrectors 옵션을 10으로 늘리고 수행해도 iteration 1~9는 initial residual이 줄어들다가 iteration 10에서는 다시 늘어난다는 것입니다.

    Solver를 뜯어봐도 final iteration이라고 특별히 residual이 증가해야 할 이유를 잘 모르겠는데,
    혹시 같은 문제를 경험하신 유저 계신가요?

    이거 왜 이럴까요?

    #770

    pig
    Participant

    우선, 답변이 아니라 죄송합니다..
    흥미로운 질문이라 저도 결과가 궁금하네요..

    #783

    geonhong
    Participant

    Momentum equation의 relaxation의 문제였네요

    pimpleControl.H에서 finalIteration의 경우 fvMesh data에 “finalIteration = true” 를 등록하는데,

    UEqn.relax() 과정에서 mesh의 data에 finalIteration이 등록되어 있으면 relax 시키지 않는 것으로 보입니다.

    • This reply was modified 3 years, 10 months ago by  geonhong.
Viewing 3 posts - 1 through 3 (of 3 total)

You must be logged in to reply to this topic.