libmesh 实现流程

2024-04-06 17:48:59

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.

libmesh 实现流程相关推荐

淘宝获取单笔订单信息服务端调用API及流程
淘宝获取单笔交易接口(文档地址):https://open.taobao.com/api.htm?docId=54&docType=2 调用接口所需依赖(文档地址):https://devel ...
用伪代码模拟洗衣机的运转流程
今天的软导课又学到了不少"骚操作",其中就包括Pseudocode和Top-down design. 不如现在就借着介绍洗衣机的运转流程向大家介绍一下这两个简单的东西. 题目如下 ...
vue-devTools插件安装流程
vue-devTools插件安装流程本文主要介绍 vue的调试工具 vue-devtools 的安装和使用工欲善其事, 必先利其器, 快快一起来用vue-devtools来调试开发你的vue项目吧 ...
RPC 笔记（01）— RPC概念、调用流程、RPC 与 Restful API 区别
1. 基本概念 PRC 远程过程调用 Remote Procedure Call,其就是一个节点请求另外一个节点提供的服务.当两个物理分离的子系统需要建立逻辑上的关联时,RPC 是牵线搭桥的常见技术手 ...
etcd 笔记（05）— etcd 代码结构、各模块功能、整体架构、各模块之间的交互、请求和应答流程
1. etcd 项目结构和功能 etcd 项目代码的目录结构如下: $ tree ├── auth ├── build ├── client ├── clientv3 ├── contrib ├── ...
浅显易懂 Makefile 入门（01）— 什么是Makefile、为什么要用Makefile、Makefile规则、Makefile流程如何实现增量编译
1. 什么是 Makefile Makefile 文件描述了 Linux 系统下 C/C++ 工程的编译规则,它用来自动化编译 C/C++ 项目.一旦写编写好 Makefile 文件,只需要一个 ma ...
Go 学习笔记（57）— Go 第三方库之 amqp （RabbitMQ 生产者、消费者整个流程）
1. 安装 rabbitmq 的 golang 包 golang 可使用库 github.com/streadway/amqp 操作 rabbitmq .使用下面命令安装 RabbitMQ . go ...
简述Web程序开发流程
总体开发流程分析需求, 列出功能清单或写需求说明书. 设计程序功能, 写功能规格书和技术规格书. 进入开发与测试的迭代. 调试和性能等专项测试. 部署上线运维前端开发流程根据功能规格绘制页面草 ...
【Sql Server】Database-sql语言的流程控制语句
流程控制语句 If else语句 If (表.列) {语句|语句块begin,,,end} Else {语句|语句块begin,,,end } If else 语句嵌套 If Begin If ...
pytorch 调用forward 的具体流程
forward方法的具体流程: 以一个Module为例: 1. 调用module的call方法 2. module的call里面调用module的forward方法 3. forward里面如果碰到M ...

最新文章

热门文章