baram icon indicating copy to clipboard operation
baram copied to clipboard

force monitor: reference pressure issue

Open kjrstory opened this issue 1 year ago • 0 comments

문서페이지의 튜토리얼 중 2번째 문제(https://baramcfd.org/tutorials/2023/09/05/AhmedBody-post/) 인 Ahmed Body를 하던 중에 Force 모니터가 생각한대로 출력이 되지 않습니다. 그림을 보면 (아직 수렴이 다 되지 않았지만) Cd는 약 0.32로 실험값과 비슷한데 비해 Cl은 약 -2.5로 실험값과 부호와 크기가 큰 차이가 납니다.
Capture1 GUI에서 Reference Values 항목의 Operating Pressure값을 바꿔도 Cl이 동일합니다.
모니터 Plot과 postProcessing폴더의 텍스트 파일 둘 다 확인하였습니다. 이에 대한 원인을 분석 중 pRef가 잘못 들어가는 것이 아닌가 추정됩니다.
baram에 의해 만들어진 control dict의 force function을 살펴보면

 patches
       (
         ...
       );
     rho rho;
     Aref 0.056;
     lRef 1;
     magUInf 40;
     rhoInf 1.2;
     dragDir
 (1 0 0);
     liftDir
 (0 1 0);
     CofR
 (0 0 0);
 ...

이런식으로 되어있습니다.
여기에 pRef값이 없습니다.
텍스트 에디터로 강제로 아래 문구를 추가하면 (forces, forceCoeffs 둘 다)

patches
      (
        ...
      );
// Reference pressure [Pa]
    pRef            101325;
    rho rho;
...

다음 그림과 같이 Cl이 0.37정도로 다른 참고문헌들과 비슷해집니다.

Capture2

윈도우 v24.1.3에서 테스트 하였습니다. 리눅스에서는 값이 이상한 것만 확인하였습니다. libforce 자체가 이상한 것 같진 않고 controldict를 작성하는 부분만 잘못된 것 아닌가 생각합니다.

kjrstory avatar May 13 '24 11:05 kjrstory