«finite-element» 태그된 질문

나는 방금 FEniCS를 엉망으로 만들기 시작했습니다. 3 차 요소로 Poisson을 해결하고 결과를 시각화하고 싶습니다. 그러나 plot (u)를 사용할 때 시각화는 결과의 선형 보간입니다. VTK로 출력 할 때도 같은 결과가 나타납니다. 내가 사용하고있는 다른 코드에서는 VTK 출력기를 작성하여 고차원 요소를 업 샘플링하여 실제로 Paraview에서 고차원으로 보였습니다. FEniCS에 이와 같은 (또는 더 …

14 finite-element  fenics 

나는 Hesthaven / Warburton 책을 사용하여 DG-FEM 방법에 대한 이론을 배우고 있으며 '숫자 플럭스'의 역할에 대해 약간 혼란 스럽습니다. 이것이 기본적인 질문이라면 사과하지만, 그것에 대한 만족스러운 답변을 찾지 못했습니다. 선형 스칼라 웨이브 방정식을 고려하십시오. 선형 플럭스로 주어진다F(U)=U.∂유∂티+ ∂에프( u )∂엑스= 0∂u∂t+∂f(u)∂x=0\frac{\partial u}{\partial t} + \frac{\partial f(u)}{\partial x} = 0에프( U …

13 finite-element  discontinuous-galerkin 

FEM 시뮬레이션을 개발 중입니다. 초기 테스트를 위해 간단한 자체 작성 메쉬 및 메쉬 그래프 시각화를 사용합니다. 그러나 기존 메셔에서 생성 한 데이터를 사용하여 기존 시각화 도구로 출력하도록 프로그램을 준비하고 싶습니다. 파일 형식에 권장 (준) 표준과 (FEM) 메쉬에 대한 내부 데이터 형식이 있습니까?

13 finite-element  mesh  graph-theory 

양식의 문제로 시작하자 (L+k2)u=0(L+k2)u=0(\mathcal{L} + k^2) u=0 주어진 경계 조건 세트 ( Dirichlet , Neumann , Robin , Periodic , Bloch-Periodic ). 이는 일부 기하학 및 경계 조건 에서 일부 연산자 대한 고유 값과 고유 벡터를 찾는 것과 일치 합니다. 예를 들어, 음향, 전자기, 탄성 역학, 양자 역학에서 이와 같은 …

13 pde  finite-element  eigensystem  verification 

파동 방정식에서 : 씨2∇ ⋅ ∇ u ( x , t ) - ∂2u ( x , t )∂티2= f( x , t )c2∇⋅∇u(x,t)−∂2u(x,t)∂t2=f(x,t)c^2 \nabla \cdot \nabla u(x,t) - \frac{\partial^2 u(x,t)}{\partial t^2} = f(x,t) 통합하기 전에 먼저 테스트 함수 v ( x , t )v(x,t)v(x,t) 를 곱하는 이유는 무엇 입니까?

13 finite-element 

질문 : 유한 요소 행렬의 희소 구조를 정확하고 효율적으로 계산하기 위해 어떤 방법을 사용할 수 있습니까? 정보 : 저는 2 중 라그랑주 단위로 Galerkin의 방법을 사용하고 C로 작성하고 희소 행렬 저장 및 KSP 루틴에 PETSc를 사용하여 포아송 압력 방정식 솔버를 작업하고 있습니다. PETSc를 효율적으로 사용하려면 전역 강성 매트릭스를위한 메모리를 미리 …

13 matrix  finite-element  petsc  performance 

\newcommand{\v}[1]{\boldsymbol{#1}} 다음과 같은 스톡 흐름 모델 방정식이 있다고 가정합니다. { − d i v ( ν∇ u ) + ∇ pD 나 v에 U를= f= 0{−div(ν∇u)+∇p=fdivu=0 \tag{1} \left\{ \begin{aligned} -\mathrm{div}(\nu \nabla \v{u}) + \nabla p &= \v{f} \\ \mathrm{div} \v{u} &= 0 \end{aligned} \right. 여기서 점도 ν( x )ν(x)\nu(x) 는 …

13 finite-element 

압축 가능한 오일러 방정식에 대한 자체 솔버를 작성하고 싶습니다. 가장 중요한 것은 모든 상황에서 강력하게 작동하기를 원합니다. 나는 그것을 FE 기반으로하고 싶습니다 (DG는 괜찮습니다). 가능한 방법은 무엇입니까? 나는 0 차 DG (유한 체적)를하고 있다는 것을 알고 있으며 매우 강력하게 작동해야합니다. 기본 FVM 솔버를 구현했으며 훌륭하게 작동하지만 수렴이 매우 느립니다. 그러나 …

13 pde  hyperbolic-pde  finite-element 

나는 사람들이 종종 일관된 질량 행렬을 덩어리 대각선 행렬로 대체한다는 것을 알고 있습니다. 과거에는로드 벡터가 FEM 일관된 방식이 아닌 일괄 방식으로 조립되는 코드도 구현했습니다. 그러나 나는 왜 우리가 처음에 이것을 할 수 있는지 조사하지 않았습니다. 질량 및 하중 벡터에 적용 할 수있는 덩어리의 직관은 무엇입니까? 그것에 대한 수학적 정당성은 무엇입니까? …

13 finite-element 

이것은 SE의 소프트웨어 권장 사항 측면에 더 적합한 질문 일 수 있지만 SE의이 부분을 자주 이용하는 사람들이이 질문에 대답 할 가능성이 더 높다고 생각합니다. Comsol Multiphysics 의 무료 (자유뿐만 아니라) 대안을 찾고 있습니다. 까다로운 부분은 다음과 같습니다. 나는 단순히로드가있는 모델링 및 시뮬레이션 패키지를 찾는 것이 아니라 가능한 한 Comsol과 비슷한 …

12 finite-element  comsol  octave  openmodelica  modeling 

내 시스템은 lagrange multipliers와 같은 대칭 FE 문제입니다 (예 : 비 압축 Stokes 흐름). ( A비비티씨)(ABTBC)\begin{pmatrix}A & B^T \\ B & C\end{pmatrix} 여기서 이 일반적인 경우입니다 (라그랑 지 승수가 마지막에 나타나도록 방정식의 번호를 매겼습니다). 시스템은 상당히 큽니다 (+ 100k 라인).씨= 0C=0C = 0 이 질문에 대한 답변을 읽은 후 혼합 …

12 finite-element  petsc  krylov-method  preconditioning 

FEM은-이산화하고, 반응 확산 문제, 예를 들어, 풀 때 0 < ε ≪ 1 (단수 변동)을 갖는 ∂ Ω 에서 u = 0 인 경우, 이산 문제의 해는 일반적으로 경계에 가까운 진동 층을 나타냅니다. 함께 Ω = ( 0 , 1 ) , ε = 10 - 5 및 선형 유한 …

12 finite-element  discretization  numerical-analysis  singular-perturbation 

3 노드의 4 노드 요소에 다항식을 통합하고 싶습니다. FEA에 관한 몇 권의 책은 임의의 평평한 4-noned 요소에 대해 통합이 수행되는 경우를 다룹니다. 이 경우 일반적인 절차는 Jacobi 행렬을 찾아 적분 기준을 더 간단한 적분 한계 [-1; 1]가 있고 가우스 레전드 르 직교 기법이 쉽게 사용되는 정규화 된 것으로 변경하는 결정 …

12 finite-element  quadrature 

레이놀즈 수가 많으면 매우 얇은 경계층이 생성됩니다. Large Eddy Simulation에서 벽 해상도를 사용하는 경우 종횡비는 정도일 수 있습니다 10610610^6. 이 정권에서는 종횡비의 제곱근이 나빠질수록 주입 상수가 저하되기 때문에 많은 방법이 불안정 해집니다. inf-sup 상수는 선형 시스템의 조건 수와 이산 솔루션의 근사 특성에 영향을주기 때문에 중요합니다. 특히, 다음과 같은 불연속 오류 …

12 pde  finite-element  fluid-dynamics  stability  incompressible 

유한 요소 방법에 관한 문헌을 읽을 때 종종 "매달려있는 노드"라는 용어가 나타날 수 있습니다. 누구든지 실제로 교수형 노드가 무엇인지 말해 줄 수 있습니까?

11 finite-element  mesh-generation