불연속 Galerkin 방식의 CFL 조건


9

나는 다음과 같은 유형의 보존 법칙의 선형 시스템을 해결하기 위해 ADER-Discontinuous Galerkin 체계를 구현했습니다. tU+AxU+ByU=0CFL 조건이 매우 제한적이라는 것을 관찰했습니다. 참고 문헌에서 시간 단계의 상한Δthd(2N+1)λmax 어디에서 찾을 수 있습니까 h 셀 크기입니다. d 차원의 수이며 N 다항식의 최대 차수입니다.

이 문제를 피할 수있는 방법이 있습니까? WENO-ADER Finite 볼륨 구성표로 작업하고 있었고 CFL 제한이 훨씬 완화되었습니다. 예를 들어, 5 차 체계의 경우 DG를 사용할 때 0.04보다 낮은 CFL을 부과해야하지만 WEFL-ADER FV 체계에서 CFL = 0.4를 계속 사용할 수 있습니다.

예를 들어, 계산 에어로 아쿠 스틱 (선형 오일러 방정식) 또는 유사한 응용 (가스 역학, 얕은 물, 자기 유체 역학)에서 ADER-FV 대신 DG 방식을 사용하는 이유는 무엇입니까? 훨씬 낮은 시간 단계에도 불구하고 체계의 전체 계산 비용이 ADER-FV의 계산 비용과 비슷합니까?

이에 대한 생각과 제안을 환영합니다.

답변:


6

DG 방식의 제한적 CFL은 일반적으로 높은 정확도와 소형 스텐실의 조합에서 비롯됩니다 ( 예 : 이 참조 참조 ). CFL은 다양한 형태의 경계에 따라L2도함수와 미량의 다항식에 의존하는 솔루션의 표준. Bernstein 또는 Markov 형제 불평등과 불연속 미량 불평등을 사용하여 이러한 각 수량의 경계는 상수에 반비례합니다.h 그리고 이차적으로 N전체 CFL이 O(h/N2).

참고로, 이전에 언급 한 CFL을 보았지만 입증 된 위치를 기억할 수 없습니다. 나는 그들이 이차 의존성을 피하는 방법을 알고 싶습니다.N 그들의 경계에.

유한 한 차이와 WENO 스킴 (주기적 메쉬의 B- 스플라인 기반 유한 요소 방법)은 CFL 조건이 느슨합니다. 유사 경계의 상수는 N. 이는 스텐실 크기가 순서에 따라 증가하는 경향이 있기 때문입니다.N이 문제 중 일부가 줄어 듭니다.

DG 방법은 더 비싸지 만 구조화되지 않은 메시를 쉽게 처리 할 수 ​​있으며 효율적으로 구현할 수 있습니다. 구조화되지 않은 그리드에는 WENO (또는 유사한 재구성)의 상위 버전이 있지만 추가 수학 또는 구현상의 복잡성을 유발할 수 있습니다.


Jesse에게 자세한 답변을 보내 주셔서 감사합니다.이 문제에 대한 폭 넓은 견해를 제시했습니다. DG-ADER를 사용한 수치 시험에서 구조적 사변형 메쉬 (예 : 사각형, 사다리꼴 또는 평행 사변형 등 임의의 사변형)를 사용할 때 수치 솔루션이 진동하지 않고 정확한 솔루션에 수렴한다는 것을 알았습니다. 그러나 구조화되지 않은 메시로 이동할 때 구조화 된 메시의 노드를 조금씩 변위시켜 준 구조화 된 메시에서도 진동이 나타납니다. 이것이 예상되는 동작입니까?
Adr

1
@Adrian-균일 한 메쉬에서 멀어지면 진동이 나타나는 것이 일반적입니다. 일반 메쉬를 사용하면 메쉬 크기로 정확히 무엇을 의미하는지 더 이상 명확하지 않습니다.h. 셀 지름, 가장 짧은 가장자리의 길이, 면적의 제곱근 (2d) 또는 "메쉬 크기"를 정의하는 다른 방법 일 수 있습니다.
Wolfgang Bangerth 님이

설명을 해주셔서 감사합니다. Wolfgang. 지금까지 나는 설정했다h가장 짧은 가장자리의 길이로. 그러나 어쨌든 공식으로 주어진 처방 된 CFL에서 CFL 수를 1 배 이상 줄이면 여전히 진동합니다.
Adr
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.