目    录

目    录 1

1.  问题描述 1

2.  相关的有限元理论基础  1

2.1 二次泛函极值原理和里兹解法  1

2.2 伽辽金加权余数法  2

3.  二维圆柱绕流的有限元解法  3

3.1  圆柱绕流问题的数学模型  3

3.2 数学模型的有限元法求解  4

3.3 有限元法求解的 MATLAB 实现  8

3.4  结果讨论  21

参 考 文 献  23

1. 问题描述

图 1 二维圆柱绕流问题

本作业尝试探索二维圆柱绕流的流函数的有限元数值解法及并绘出其流线图。

2.相关的有限元理论基础

2.1 二次泛函极值原理和里兹解法

设微分方程

Au f ,(在区域中) ,                            (1)

式中 A 为线性微分算子,如拉普拉斯算子等,它作用于函数 u 后产生函数 f ,这里 f 是给定

的已知函数。如果在区域内任意两函数 u、v 满足如下内积等式,则称算子 A 为对称的,即

(u, Av) uAvd                              (2)

如果(u, Au) 0 ,那么称A 是正定的。

如果算子 A 是对称正定的,可以证明使二次泛函

J(u)=(Au, u)-2(f, u)                               (3)

为极小值的 u 是微分方程(1)的解。因此,解微分方程(1)可以转化为解泛函(3)的极小值问题。

而如果算子 A 是负定的,即(u, Au) 0 ,那么泛函(3)的极大值 u 即为微分方程(1)的解。

例如对于泊松方程第一边值问题。

2 2 

 

u u

-( 2 + 2 )=f(x,y)x, y, 

x y

                         (4)

u 0

s 

这里 S 为区域 Ω的边界。

1

可以证明,对于泊松方程第一边值问题,算子A 2 是对称的,正定的。

这样可以按照(3)式建立泛函并整理

J(u) (Au, u) 2(f, u)

2

matlab圆柱饶流,有限元法解二维圆柱绕流问题.pdf相关推荐

  1. 二维有限元方程matlab,有限元法求解二维Poisson方程的MATLAB实现

    有限元法求解二维 Poisson 方程的 MATLAB 实现 陈 莲a ,郭元辉b ,邹叶童a ( 西华师范大学 a. 数学与信息学院; b. 教育信息技术中心,四川南充 6437009) 摘 要: ...

  2. galerkin有限元法matlab实现,有限元法求解二维Poisson方程的MATLAB实现

    有限元法求解二维Poisson方程的MATLAB实现 陈莲a,郭元辉b,邹叶童a [摘要]文章讨论了圆形区域上的三角形单元剖分.有限元空间,通过变分形式离散得到有限元方程. 用MATLAB编程求得数值 ...

  3. fluent二维叶型仿真_ICEM划分嵌套网格之二维圆柱绕流

    首先,介绍一下嵌套网格.网络上关于嵌套网格的的内容大多数是关于直接利用软件进行计算的过程,而对于前处理过程中的网格生成过程并没有什么描述,其实这种技术已经在学术界流传已久,只是用的都是自己的程序算法, ...

  4. Algorithm之MC:基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动

    Algorithm之MC:基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动 目录 输出结果 设计代码 输出结果 基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动 设计代码 %基于Matlab实 ...

  5. LBM学习记录5 Python实现IB二维圆柱绕流 1.0

    浸润边界(IB)下二维圆柱绕流,bug1.0版 import numpy as np from math import sqrt from numpy import pi as PI from num ...

  6. Matlab绘制三维曲面(以二维高斯函数为例)

    原文地址为: Matlab绘制三维曲面(以二维高斯函数为例) 寒假学习了一下Python下的NumPy和pymatlab,感觉不是很容易上手.来学校之后,决定继续看完数字图像处理一书.还是想按照上学期 ...

  7. python下二进制流文件转二维码

    python下二进制流文件转二维码 之前做过的一个ctf中的misc,经过一系列操作得到一个flag.txt 里面是一串二进制流 突然束手无策 经简单计算 260^2=67600 所以尝试能否转化为一 ...

  8. 二维计算机动画,二维计算机动画制作新技术.pdf

    二维计算机动画制作新技术.pdf w年7月 浙江教育学院学报 2007 :uly 第4期 INSTIllFIE No.4 JoURNALOFZHF-JIANGEDUCATION 二维计算机动画制作新技 ...

  9. 生成二维码并将二维码附到pdf文件上,扫描二维码完成跳转

    依托e-iceblue工具完成 官网:https://www.e-iceblue.cn/tutorials.html 本文参考:https://www.e-iceblue.cn/pdf_java_ba ...

最新文章

  1. No module named ‘skimage.metrics‘在Anaconda3中的解决方法
  2. ThunderNet:国防科大、旷视提出首个在ARM上实时运行的通用目标检测算法
  3. iPhone第二季度在美销量降23%,但iPhone SE有亮点
  4. c语言实现1024点fft程序,C语言1024点快速傅里叶变换(FFT)程序,最好经过优化,执行速度快...
  5. psql json操作符合函数
  6. k8s核心技术-资源编排(yaml)的编写方式---K8S_Google工作笔记0019
  7. 动态向客户端注册脚本文件
  8. Windows下安装最新的Apache+PHP+MySQL方法--记录方便自己参考
  9. 洛谷——P1219 [USACO1.5]八皇后 Checker Challenge
  10. Axure的使用---原型图
  11. 大写日期转换器(大写日期转换器)
  12. ISO9001认证审核常见问题大汇总
  13. VS2010提示asp.net v4.0 尚未在web服务器上注册
  14. 微信小程序 图片上传预览删除
  15. 监听浏览器刷新事件,拦截浏览器返回,js监听移动端浏览器页面显示、隐藏
  16. iPad mini 6 评测:多功能升级,一手掌握的极佳手感,适合娱乐休闲
  17. 抖音服务器升级中无法修改名字,抖音名字改不了怎么回事
  18. 服装销售系统计算机毕业论文,大学本科计算机专业服装销售管理系统设计毕业论文.doc...
  19. 前端失业 3个月,尝试接私单的感触
  20. 30行代码实现蚂蚁森林自动偷能量

热门文章

  1. 交换两个变量的值的效率的测试
  2. android 豆瓣客户端 视频
  3. MYSQL函数应用----替换函数replace()用法
  4. context.Request.Files超过了最大请求长度
  5. JSON数据格式----- JavaScript与JSON、JavaScript的JSON对象、构建JSON格式数据
  6. JavaScript算法【特殊篇】数组去重的十种实用方案
  7. origin2016中怎么画多条曲线,并且分别给不同曲线设置标记
  8. linux 删除in链接,linux - 从TCP_IN删除端口不会从CSF上的外部流量关闭它 - 堆栈内存溢出...
  9. flutter不支持热更新_真当Flutter不能热更新?众能动态化Flutter
  10. KNN算法——分类部分