Właśnie zacząłem bawić się w FEniCS. Rozwiązuję Poissona z elementami trzeciego rzędu i chciałbym wizualizować wyniki. Kiedy jednak używam wykresu (u), wizualizacja jest po prostu liniową interpolacją wyników. To samo pojawia się, gdy przesyłam dane do VTK. W innym kodzie, z którym pracuję,...
Mam problem fizyczny rządzony równaniem Poissona w dwóch wymiarach Mam pomiary dwóch składników gradientu ∂ u / ∂ x i ∂ u / ∂ y wzdłuż pewnej części granicy, Γ m , więc chciałbym narzucić ∂ u−∇2u=f(x,y),inΩ−∇2u=fa(x,y),janΩ -\nabla^2 u = f(x,y), \; in \; \Omega ∂u /...
Próbuję rozwiązać równanie Poissona 2D za pomocą nieciągłej metody Galerkina (DG) i następującej dyskretyzacji (mam plik png, ale nie mogę go przesłać, przepraszam): Równanie: ∇ ⋅ ( κ ∇ T) + f= 0∇⋅(κ∇T)+f=0\nabla \cdot( \kappa \nabla T) + f = 0 Nowe równania: q= κ ∇ T∇ ⋅ q= - fq=κ∇T∇⋅q=−fq =...
Zasadniczo FEM wydaje się problemem, który jest w zasadzie „rozwiązany”. Istnieje wiele potężnych frameworków, takich jak Trilinos, PETSc, FEniCS, Libmesh lub MOOSE. Łączy je jedno: są wyjątkowo „ciężkie”. Po pierwsze, instalacja zwykle jest bardzo bolesna. Po drugie, ich interfejs / API jest...