유체 시뮬레이션 결과의 정확성을 어떻게 확인할 수 있습니까?


9

입자 기반 유체 시뮬레이션 프로그램을 작성했습니다. 내가 올바른 결과를 얻는 지 말하기가 어렵습니다. 시각화 된 결과는 합리적으로 보이지만 일부는 이상하게 보입니다. 그것이 유체의 특징인지는 모르겠습니다. 내 프로그램이 올바른지 확인하는 정확한 방법이 있습니까?

세부 사항 수정 :

내 프로그램은 2D 파티클 기반 시뮬레이션 프로그램입니다. 유체는 압축 가능합니다. 구현은 거의 고전적인 논문을 기반으로합니다.

ül 러, 마티아스, 데이비드 샤리 파, 마르쿠스 그로스. "대화 형 응용 분야를위한 입자 기반 유체 시뮬레이션." 2003 ACM SIGGRAPH의 절차

반복 방법으로 Navier-Stokes 방정식을 풀었습니다. 압력, 중력, 점도 및 표면 장력 만 고려했습니다.


숫자 미분법으로 NS 방정식의 항을 다시 계산하고 어떻게 취소하는지 확인할 수 있습니다.
Yves Daoust

답변:


2

다른 사람의 소프트웨어와 비교하십시오. 일부 표준화 된 테스트를 실행하고 다른 테스트와 거의 동일한 답변이 있는지 확인하십시오. 당신이 같은 대답을 얻는다면, 코드를 올바르게 가질 확률보다 상당히 높습니다.

일부 테스트 :

  1. 실린더를지나 흐름. 2D에서 직사각형 영역, 중앙의 실린더, 왼쪽의 유입, 싸움의 유출을 취하고 실린더의 힘을 계산하십시오. 다음은 소수의 코드를 비교하는 벤치 마크 입니다.

  2. 부력 흐름. 닫힌 상자, 하단의 핫 플레이트, 상단의 콜드 플레이트, 부력으로 인해 뜨거운 유체가 상승하기 시작합니다. 여기 벤치 마크가 있습니다.

  3. 기포 상승, 벤치 마크 .

그러나 불행히도 해당 벤치 마크에서 코드를 과학 코드와 비교하는 것은 매우 어려울 수 있습니다. 나는 당신이 SPH 또는 안정된 유체가 아닌 정확성을 위해 만들어진 것을 구현했다고 생각합니다.

실린더를 통과하는 흐름을 예로 들어 보겠습니다. 아주 작은 레이놀즈 수로 테스트를 시작한 다음 시뮬레이션 정확도를 높이면 (더 낮은 시간 단계, 세분화 또는 입자 수 증가) 실린더의 힘을 측정합니다. 힘이 어떤 수로 수렴합니까? 그렇지 않다면 문제가있는 것보다, 예인 경우 벤치 마크 용지를보고 다른 사람과 결과를 비교하십시오.


이 방법은 레이트 레이서 테스트에 사용하는 것과 매우 유사한 기술입니다. 다른 사람의 렌더러로 테스트 장면을 렌더링하고 내 결과와 비교합니다. 그들은 같은 결과에 수렴합니까? 맞다면 내가 옳은 것보다 아니라면, 틀린 것보다 그렇습니다.


알려진 실제 측정 및 유체 역학 벤치 마크에 대해 예스를 테스트하는 소프트웨어 대신. 그렇지 않으면 오류가 오염됩니다. 같은 질문을 stackexhange 네트워크 btw의 다른 곳에 게시했습니다.
joojaa

1
물리 측정이 올바른 경우 실제 측정에 대한 테스트가 테스트에 적합하다고 생각합니다. 프로그램 만 디버깅하려면 다른 코드를 테스트하는 것이 좋습니다. 또한 컴퓨터 시뮬레이션에서 실험에 영향을주지 않고 무엇이든 측정 할 수 있습니다. 예를 들어 어떤 시점에서 유체 속도를 측정하는 것은 실제 실험에서는 불가능하지만 컴퓨터 시뮬레이션에서는 쉽지 않습니다.
tom

그렇습니다. 그러나 당신도 그들의 솔버 문제를 물려받습니다. 나는 멀티 바디 시뮬레이터를 개발하고 MSC Adams의 결과를 확인하지만이 사실은
별다른

실제 실험을 확인하는 것이 더 좋습니까? 나는 그것을 의심하지만 잘못되었을 수 있습니다. 멀티 바디 물리의 상황은 유체 물리와는 상당히 다릅니다. 당구처럼 단순한 것조차 혼란스러운 행동을합니다. 더욱이 접촉이있는 강체 동역학은 수학 문제도 잘 제기되지 않습니다. Painlevé 역설을 알고 있습니까? 따라서 멀티 바디 물리학의 수치 시뮬레이션은 일반적으로 실패 할 것입니다. 참고 문헌 : plus.maths.org/content/chaos-billiard-table en.wikipedia.org/wiki/Painlev%C3%A9_paradox
tom

1
예, 나는 멀티 바디 다이나믹이 어떻게 작동하는지 알고 있습니다. 그러나 알려진 분석 솔루션을 확인하는 것이 쉽지 않았습니다. 그러나 실제 체액은 다체 역학과 유사하게 혼란입니다. 따라서 층류 상황 등을 점검 할 수 있어야합니다. 마찰은 암캐입니다.
joojaa
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.