메쉬 의존적 안정성을 가진 요소의 유용성


15

3D Stokes 문제에서 요소의 안정성과 관련된 수학을 수행 한 후 이 임의의 4 면체 메쉬에 대해 안정적이지 않다는 것을 깨닫기 위해 약간 충격을 받았습니다 . 더 정확하게 말하면, 모든 노드와 4 개의 패싯 중 3 개의 패싯이 Dirichlet 조건이있는 도메인의 경계에있는 요소가있는 경우, 단일 행렬을 얻게됩니다. 스토크 스 시스템의 약한 형태로 결론을 내리는 것은 실제로 사소한 일입니다.피2−피1

(COMSOL)에 액세스 할 수있는 유일한 상용 Stokes 코드를 테스트했으며 그러한 메쉬를 만들 수있었습니다. 해결을 클릭하면 예상대로 '오류 : 특이 행렬'이 나타납니다. (COMSOL은 크리핑 흐름 모듈에 을 사용한다는 인상을 받았습니다 .)피2−피1

나쁜 메쉬.

문제가 다른 구성과 관련이 없다는 것을 추가로 테스트하기 위해 다음 메쉬를 시도하고 모든 것이 예상대로 작동합니다.

여기에 이미지 설명을 입력하십시오

질문 : (적응 형 또는 비 적응 형) 메시 생성기에서 이러한 종류의 제약이 고려됩니까? 다양한 연구 논문 에서이 요소가 꽤 인기있는 것으로 보입니다. 이러한 종류의 경계 불안정성은 일반적으로 사용할 방법을 선택할 때 중요하지 않은 것으로 간주됩니까? 더 중요한 것은 안정적인 유한 요소 를 갖는 것이 실제로 무엇을 의미 합니까? 즉, 어떤 종류의 메시 종속 불안정성이 너무 많아서 방법이 나쁘다는 결론을 내릴 수 있습니까?


5
재미있는 질문! 내가 아는 한, 이러한 요소는 일반적으로 큐브에서 구조화 된 사면체 메쉬 생성의 결과이며 구조화되지 않은 노드 화 알고리즘이있는 실제 응용 프로그램에서는 작은 역할 만합니다. 나는 얼마 전에 조금 시도했지만 완전히 구조화되지 않은 메쉬를 생성하는 메쉬 생성기로 이러한 메쉬를 생성 할 수 없었습니다. 나는 그들이 과도하게 구속 된 요소를 피하기위한 메커니즘을 사용하고 있다고 생각합니다. COMSOL에 액세스 할 수는 없지만 대부분의 솔버에서는이 요소가 큰 문제를 일으키지 않는다고 생각합니다.
— Christian Waluga 2016 년

이것이 MINI 요소에 문제가 있는지 궁금합니다.
— Daniel Shapero

1
MINI에는이 문제가 없다고 생각합니다. 단일 내부 DOF는 상황을 저장합니다. 스토크의 고유성 조건은 입니다. 보자 . 선택 여기서 일부 테트라에있는 기포. 이것은 가 국소 상수이며 연속성이 나머지를 처리 한다는 것을 의미 합니다. p ( x , y ) = a + b x + c y v = ( b ϕ , c ϕ ) ϕ p(∇⋅V,피)=0 ∀V∈Vh⇒피=글로벌 const.피(엑스,와이)=ㅏ+비엑스+씨와이V=(비ϕ,씨ϕ)ϕ피
— knl

답변:


1

이 문제는 삼각형의 모든 정점이 경계에있을 때 2D에서도 발생합니다. 스톡스 문제는 메쉬가 실패 할 수있는 유일한 문제는 아니다하는 의 일부 값 -Laplacian 문제 소자의 제로 기울기를 찾아 실패한다.p피피

발전기는 일반적으로이를 처리 할 수있는 옵션이 메쉬, 예를 들면 2 차원 메시 생성기 bamg의은 freefem++갖는 -splitpbedge경계에 양단을 갖는 임의의 에지의 중간 노드를 추가하는 옵션. bamg문서 에 따르면 비정형 메쉬 생성 은 이러한 삼각형을 반환 할 수 있습니다 .


2D Stokes의 Taylor-Hood와 같은 경우입니까? 내 직감은 가장자리와 관련된 DOF가 상황을 저장한다고 알려줍니다. 3D Taylor-Hood에서는 패싯과 관련된 DOF가 없으므로 불안정성이 발생합니다.
— knl

당신 말이 맞아요, 그럴 수도 있습니다 Taylor-Hood의 inf-sup 조건에 대한 Verfuhrt의 증거는 이것을 확인하기에 충분히 건설적이지만 지금은 시간이 없다고 생각합니다.
— Joce
당사 사이트를 사용함과 동시에 당사의 쿠키 정책과 개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.