Questions tagged «fenics»

我才刚刚开始使用FEniCS。我正在用三阶元素求解泊松，并想将结果可视化。但是，当我使用plot（u）时，可视化只是结果的线性插值。当输出到VTK时，我会得到相同的结果。在我正在使用的另一个代码中，我编写了一个VTK输出器，该输出器将对高阶元素进行升采样，以使它们实际上在Paraview中看起来较高阶。FEniCS中是否有类似（或更好）的东西？

14 finite-element  fenics 

我具有由泊松方程在两个维度支配的物理问题 我有两个梯度分量的测量 ∂ ü / ∂ X和 ∂ ú / ∂ ý沿边界的某些部分， Γ 米，所以想强加 ∂ û- ∇2u = f（x ，y），我ñΩ−∇2u=f(x,y),inΩ -\nabla^2 u = f(x,y), \; in \; \Omega ∂你/ ∂X∂u/∂x\partial{u}/\partial{x}∂你/ ∂ÿ∂u/∂y\partial{u}/\partial{y}Γ米Γm\Gamma_m 和传播到远场。∂ü∂X一世0= 克米，Ø ñΓ米∂u∂xi0=gm,onΓm \frac{\partial{u}}{\partial{x_i}}_0 = g_m, \; on \; \Gamma_m 切向梯度分量，，我可以整合，然后通过狄利克雷条件执行，使得 ＆Integral;Γ米∂ù∂ü∂X（t ，0 ）∂u∂x(t,0)\frac{\partial{u}}{\partial{x}}_{(t,0)} 为了同时施加法向分量， ∂ ù∫Γ米∂ü∂X（t ，0 ）ds …

11 finite-element  python  fenics 

我正在尝试使用Discontinuous Galerkin方法（DG）和以下离散化方法求解二维Poisson方程（抱歉，我有一个png文件，但不允许上传）： 方程： ∇⋅(κ∇T)+f=0∇⋅(κ∇T)+f=0\nabla \cdot( \kappa \nabla T) + f = 0 新的方程： q=κ∇T∇⋅q=−fq=κ∇T∇⋅q=−fq = \kappa \nabla T\\\nabla \cdot q = -f 弱形式与数值通量Ť和q：T^T^\hat{T}q^q^\hat{q} ∫q⋅wdV=−∫T∇⋅(κw)dV+∫κT^n⋅wdS∫q⋅∇vdV=∫vfdV+∫q^⋅nvdS∫q⋅wdV=−∫T∇⋅(κw)dV+∫κT^n⋅wdS∫q⋅∇vdV=∫vfdV+∫q^⋅nvdS\int q \cdot w dV = - \int T \nabla \cdot (\kappa w) dV + \int \kappa \hat{T} n \cdot w dS\\ \int q \cdot \nabla v …

10 pde  finite-element  fenics 

基本上，FEM似乎是一个“已解决”的问题。现有许多强大的框架，例如Trilinos，PETSc，FEniCS，Libmesh或MOOSE。 它们的共同点是：它们非常“重”。首先，安装通常非常麻烦。其次，它们的接口/ API很繁琐-您必须将整个想法转化为相应库的思想。这也意味着很难满足特殊要求或现有代码的互操作性和可扩展性。 其他项目（例如随机示例）Boost，LibIGL，Aztec（线性求解器），Eigen或CGAL证明，绝对有可能编写功能强大的库，这些接口可以通过非常精简的界面无缝集成到C ++或Python代码中，而无需安装超重的框架。 FEM是否有真正轻巧的包装？我不是在寻找简单，自动的求解器-我是在寻找一个提供强大功能，同时保持精益接口，与常见数据结构（例如C ++ STL）互操作性以及轻量级安装（例如仅标题）的库。

9 finite-element  python  c++  petsc  fenics