タグ付けされた質問 「fenics」

2
FeniCS:高次要素の視覚化
FEniCSをいじり始めたばかりです。3次要素でポアソンを解いているので、結果を視覚化したいと思います。ただし、plot(u)を使用すると、視覚化は結果の単なる線形補間になります。VTKに出力するときにも同じことが起こります。私が作業している別のコードで、高次の要素をアップサンプリングして、実際にParaviewで高次に見えるようにするVTK出力機能を作成しました。FEniCSにはこのような(またはそれ以上の)ものがありますか?

3
ポアソン方程式:ラグランジュ乗数による境界条件として完全勾配を課す
私は、2次元のポアソン方程式によって支配物理的な問題を有していて 私は2つの勾配成分の測定値を有する ∂ U / ∂ X及び ∂ U / ∂ Y境界の一部に沿って、 ΓのMよう課すしたい、 ∂ uと−∇2u=f(x,y),inΩ−∇2u=f(x,y),inΩ -\nabla^2 u = f(x,y), \; in \; \Omega ∂u/∂x∂u/∂x\partial{u}/\partial{x}∂u/∂y∂u/∂y\partial{u}/\partial{y}ΓmΓm\Gamma_m あり、遠方界に伝播します。∂u∂xi0=gm,onΓm∂u∂xi0=gm,onΓm \frac{\partial{u}}{\partial{x_i}}_0 = g_m, \; on \; \Gamma_m 接線方向の傾き成分、、私は統合することができ、次いで、ディリクレ条件を介して強制よう ∫ΓM∂U∂あなたは∂バツ(t 、0 )∂u∂x(t,0)\frac{\partial{u}}{\partial{x}}_{(t,0)} 同時に法線成分を課すために、 ∂ U∫Γm∂あなたは∂バツ(t 、0 )ds = u0∫Γm∂u∂x(t,0)ds=u0 \int_{\Gamma_m}\frac{\partial{u}}{\partial{x}}_{(t,0)} \, ds = u_0 、私はラグランジュ乗数を経由する必要があります集まりました。∂あなたは∂バツ(n 、0 …

2
不連続ガラーキン/ポアソン/フェニックス
不連続ガラーキン法(DG)と次の離散化を使用して2Dポアソン方程式を解こうとしています(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 数値フラックスおよび弱い形式: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 dV …

3
「軽量」のFEMパッケージはありますか?
基本的に、FEMはかなり「解決」された問題のようです。Trilinos、PETSc、FEniCS、Libmesh、MOOSEなど、数多くの強力なフレームワークが存在します。 彼らが共通していることの1つは、それらが非常に「重い」ことです。まず、インストールは通常非常に骨の折れる作業です。第二に、それらのインターフェース/ APIは分厚くて重いです-あなたはあなたの考え全体をそれぞれのライブラリの考えに翻訳しなければなりません。これは、特別な要件や既存のコードの相互運用性と拡張性が困難であることも意味します。 (ランダムな例)Boost、LibIGL、Aztec(線形ソルバー)、Eigen、またはCGALなどの他のプロジェクトは、C ++またはPythonコードにシームレスに統合する強力なライブラリーを、非常に無駄のないクリーンなインターフェースで、インストールの必要なしに作成することが絶対に可能であることを示しています超重いフレームワークの。 FEM用の本当に軽量なパッケージはありますか?簡単なオートマジックソルバーを探しているのではありません-無駄のないインターフェース、一般的なデータ構造(たとえばC ++ STLなど)との相互運用性、および軽量のインストール(たとえばヘッダーのみ)を維持しながら強力な機能を提供するライブラリを探しています。
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.