Acabo de empezar a jugar con FEniCS. Estoy resolviendo Poisson con elementos de tercer orden y me gustaría visualizar los resultados. Sin embargo, cuando uso plot (u), la visualización es solo una interpolación lineal de los resultados. Me sale lo mismo cuando salgo a VTK. En otro código con el que...
Tengo un problema físico regido por la ecuación de Poisson en dos dimensiones Tengo mediciones de los dos componentes de gradiente ∂ u / ∂ x y ∂ u / ∂ y a lo largo de alguna parte del límite, Γ m , por lo que me gustaría imponer ∂ u- ∇2u = f( x , y) ,yo nΩ−∇2u=f(x,y),inΩ -\nabla^2 u = f(x,y),...
Estoy tratando de resolver la ecuación de Poisson en 2D usando el método Discontinuous Galerkin (DG) y la siguiente discretización (tengo un archivo png pero no puedo cargarlo, lo siento): Ecuación: ∇ ⋅ ( κ ∇ T) + f= 0∇⋅(κ∇T)+F=0 0\nabla \cdot( \kappa \nabla T) + f = 0 Nuevas ecuaciones: q= κ ∇...
Básicamente, FEM parece ser un problema que está prácticamente "resuelto". Existen numerosos marcos potentes, como Trilinos, PETSc, FEniCS, Libmesh o MOOSE. Una cosa que tienen en común: son extremadamente "pesados". Primero, la instalación normalmente es muy dolorosa. En segundo lugar, su...