libmesh 实现流程

最新推荐文章于 2023-11-13 14:59:23 发布

北冥有鱼虚空充盈

最新推荐文章于 2023-11-13 14:59:23 发布

阅读量1.2k

点赞数

分类专栏： finite element method

finite element method 专栏收录该内容

34 篇文章 15 订阅

订阅专栏

libmesh 实现流程

标签： libmesh poisson 有限元

2014-09-19 22:56 371人阅读评论(0) 收藏举报

分类：

开源 libmesh（8）

介绍poisson方程的libmesh实现过程：

Poisson Equation $- \Delta u = f$

Weak Form $(F(u),v) = \int _{\Omega} \nabla u \cdot \nabla v dV + \int_{\partial \Omega} \nabla u \cdot \textbf{n} dS$

FEM Approximation $(F(u^h),v^h) = \int _{\Omega^h} \nabla u^h \cdot \nabla v^h dV = \sum_{e=1}^{Ne} \int_{\Omega^e} \nabla u^h \cdot \nabla v^h dV^e$

Basic Function/Shape Function $v^h|_{\Omega^e} = \sum_{i=1}^N v_i \phi_i$

Approximation in Element Level $\int_{\Omega^e}v^h f dV = v_i \sum_{i=1}^N f \phi_i \\ \int_{\Omega^e} \nabla u \nabla v dV = v_i \sum_{i=1}^N \nabla \phi_j \nabla \phi_i u_j$

Matrix Notation $K^e U^e = F^e$

Local To Global Element Map $F^e = \int_{\Omega_e } f \phi_i dV = \int_{\hat \Omega_e} f(x(\xi)) \phi_i |J| d \xi \\K^e = \int_{\Omega_e } \phi_j \cdot \phi_i dV = \int_{\hat \Omega_e} \hat \nabla\phi_j \hat \nabla\phi_i |J| d \xi \\where, \hat \nabla = J^{-1} \nabla$

Integrals Numerical Quadrature with Nq points at \xi_q with weight value w_q $F^e = \int_{\hat \Omega_e} f(x(\xi)) \phi_i |J| d \xi \approx \sum_{q=1}^N f(\xi_q) \phi(\xi_q) J(\xi_q) w_q\\K^e = \int_{\hat \Omega_e} \hat \nabla\phi_j \hat \nabla\phi_i |J| d \xi \approx \sum_{q=1}^N \phi _j(\xi_q) \phi _i(\xi_q) J(\xi_q) w_q$

Libmesh provides the following variables at each Quadrature Point JxW -> J(\xi_q) w_q phi[i][q] -> \phi_i{\xi_q) dphi[i][q] -> \nabla \phi_i(\xi_q) xyz[q] -> x(\xi_q) 至此，给出了poisson方程完整的libmesh分析流程。参考libmesh slides.

北冥有鱼虚空充盈

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
libmesh 实现流程

libmesh 实现流程标签： libmesh poisson 有限元 2014-09-19 22:56 371人阅读评论(0)收藏举报分类：开源 libmesh（8）版权声明：本文为博主原创文章，未经博主允许不得转载。介绍poisson方程的libmesh实现过程： Poisson Equation
复制链接

扫一扫