文章编号:100926825(2007) 220362202二维渗流场的 MATLAB 仿真 收稿日期:20070423 作者简介:陶  承(1978) ,女 ,助理工程师 ,临安市水利水电局 ,浙江 临安  311300 凌前江(1965) ,男 ,助理工程师 ,临安市水利水电局 ,浙江 临安  311300 严  妹(19782) ,女 ,助理工程师 ,临安市水利水电局 ,浙江 临安  311300 陶 承  凌前江  严 妹 摘  要 :分析了二维渗流场与二维温度场之间的相似性 ,提出利用 MATLAB 软件的热分析功能进行非饱和土渗流的计算方法 ,阐述 MATLAB 软件中的 PDE工具箱在求解计算流体力学问题中的应用过程 ,算例结果表明应用 MATLAB 软件进行非饱和土渗流场的数值模拟是可行的。 关键词 :MATLAB ,渗流场 ,数值模拟 中图分类号 :TP391 文献标识码 :A   渗流问题的解法可分为四种 :解析法、数值解法、图解法、模型试验法 ,热比拟法是模型试验法的一种。鉴于非饱和土中水的渗流与温度场的热传导都遵循质量 (或能量) 守恒定律及微分方程满足连续性条件 ,可应用 MATLAB 软件的热分析功能进行非饱和土的渗流计算。 1  非饱和土渗流场与温度场的相似性分析 1. 1  二维非饱和土渗流的基本微分方程 依据质量守恒原理和非饱和流动的 Darcy 定律 ,当不考虑土骨架的变形并把土中水视为不可压缩流体时 ,以总水头为控制变量的基本微分方程[2]为 : x [ kwxhw x ] +y[ kwyhw y ] = chw 9t (1)其中 , hw 为总水头 ; kwx和 kwy分别为 x , y 方向的渗透系数 ; c 为比水容量。 1. 2  二维热传导的基本微分方程 假设物体中无内热源 ,根据能量守恒原理导出以温度 T 为控制变量的二维热传导基本微分方程[3] : λ x2 T x2 +λ y2 T y2 =ρ c′T 9t (2)其中 , λ x 和λ y 分别为 x , y 方向的热传导系数 ;λ为比热系 数 ;c′为单位质量物体每增加 (或降低) 1 ℃需吸收 (或释放) 的热量 ;ρ为介质密度。 1. 3  二维非饱和土渗流场与二维温度场相似性分析渗流场中存在水头梯度 ,水缓慢地从高水头的地方流向低水头的地方 ;温度场中的热能在温度梯度作用下从温度高的地方向温度低的地方传导。渗流与热传导的微分方程形式、控制变量以 及相关参数[4 ,5]如表 1 所示 ,具有很大的相似性。 表 1  渗流与热传导相关参数对照 渗流 热传导 项目 符号 单位 项目 符号 单位水头 hw m 温度 T ℃ 水量 V m3 热量 W J 水流量 Qw m3/ s 热流量 Q J/ s(或 W) 水流速度 v m/ s 热流密度 q W/ m2 渗透系数 kw m/ s 导热系数 λ W/ m·℃ 比水容量 c L/ m 比热系数 ×密度 c′ρ J/ m3·℃   渗流场中的水流量 Qw = V/ t ,温度场中热流量 Q = W/ t ;水流速度 v = Qw/ A ,其中 , A 为某断面的平均面积 ,m2 ,热流密度q = Q/ A ,由以上各式可见水流量与热量以及水流速度与热流密度具有很大的相似性 ,根据渗透系数与导热系数二者的定义 ,分别有 : kw = v hw/n (3)λ = q T/9n (4)其中 , n 为等温面或等水头面的法线方向 ,由式 (3) 和式 (4) 可见 ,渗透系数与导热系数具有相似的物理意

matlab油气田渗流,二维渗流场的MATLAB仿真相关推荐

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

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

  2. 利用matlab实现DMD动态模态分解(在一维信号或二维流场矢量中的应用)

    利用matlab实现DMD动态模态分解(在一维信号或二维流场矢量中的应用) 0 前言 0.1 特征根的计算与含义 1 DMD的基本思路 2 一维DMD算法 3 二维DMD算法 4 总结 (2020年9 ...

  3. 利用matlab实现POD分解(在一维信号或二维流场矢量中的应用)

    利用matlab实现POD分解(在一维信号或二维流场矢量中的应用) 0 前言 0.1 matlab中特征值计算 0.2 matlab中SVD分解计算 0.3 信号的正交性 1 一维信号POD分解 1. ...

  4. matlab 如何画二维图形,Matlab 学习 画图篇 一 二维图形

    matlab给绘制二维图形提供了很多的函数,把一些绘制二维图形的基本函数做成一张表,如下图所示: 我就按照表的顺序一一记录一些个函数的简略用法. 首先是 1.plot函数 plot函数有很多重载方法, ...

  5. 同轴全息matlab仿真,HoloSpec2D 二维全息谱的matlab程序,含有频谱校正 276万源代码下载- www.pudn.com...

    文件名称: HoloSpec2D下载  收藏√  [ 5  4  3  2  1 ] 开发工具: matlab 文件大小: 61 KB 上传时间: 2014-05-24 下载次数: 24 详细说明:二 ...

  6. matlab中饼图绘制程序,MATLAB中绘制二维饼图的函数是(? ? ) (1.5分)

    MATLAB中绘制二维饼图的函数是(? ? ) (1.5分) 答:pie 计算机的内存储器比外存储器读写速度 ______ 答:快 新民主主义革命时期,提出了反帝反封建这个任务就等于真正懂得了什么叫反 ...

  7. matlab中二维散点图,MATLAB实例:二维散点图

    MATLAB实例:二维散点图 用MATLAB实现简单的散点图的绘制,可以自行设置点的形状与颜色,可以去掉坐标轴边框. 数据来源:MATLAB高斯混合数据的生成 将二维数据保存为gauss.txt. 1 ...

  8. 高斯 二维 matlab,MATLAB实例:二维高斯分布图

    MATLAB实例:二维高斯分布图 问题描述:用MATLAB在一张图上同时画出二维高斯散点图与分布图. 数据来源:MATLAB高斯混合数据的生成 将二维数据保存为gauss_data.txt. 前期工作 ...

  9. MATLAB scatter 画二维/三维散点图时 用颜色表示数值大小 colorbar

    MATLAB scatter 画二维/三维散点图时 用颜色表示数值大小 自己在编程的时候,没有在网上找到很好的解释. 所以把我的代码贴出来,给遇到同样问题的盆友一个参考~ 带注释.看不懂的地方可以给我 ...

最新文章

  1. Linux-CentOS 7 增加root分区容量
  2. zigbee绑定 使用_遇见-果加智能锁F2——使用体验
  3. C# 6.0:Expression – Bodied Methods
  4. 清理linux 服务器的命令行,使用Linux上的Magic SysRq键修复冻结的X服务器,清理重新启动并运行其他低级命令 | MOS86...
  5. linux github代码仓库,centos7上使用git命令把代码放到github上
  6. android samba github,安卓手机访问树莓派samba文件共享出错解决
  7. 借助离散数学解决“哈弗大学智商测试”一题 --编程算法
  8. 小学用计算机画画的课件,用计算机画图画.ppt
  9. eplices如何导入外部代码_搜索引擎的蜘蛛是如何爬的
  10. equals 跟== 的区别
  11. 【Oracle】 expdp/impdp
  12. lopatkin俄大神精简中文系统Windows 8.1 Pro 18655 x86-x64 ZH-CN PIP
  13. 什么是UKey?Ukey在密评中的应用 双因素身份认证 安当加密
  14. Unity 查找重复图片资源以及引用 工具
  15. 中小企业生产信息化:私有系统还是云方案?
  16. 89c51 单片机 c语言 编写sszymmh 歌曲 文档,Proteus仿真51单片机生日快乐音乐播放器...
  17. 如何删除金山毒霸网址导航
  18. uniapp ios打包详细步骤
  19. TCP/UDP、TCP/IP区别
  20. 《计算机操作系统》(第四版)汤子瀛等编著——思维导图详细版本

热门文章

  1. icassp 2021 语音合成文章总结
  2. 使用 React Hooks实现一个返回顶部的按钮
  3. rib fib arp fdb
  4. oracle 添加输出数据集的函数
  5. 我的未来作文计算机行业,我的未来作文(三篇)
  6. Python学习手札01——二分法求值问题
  7. 【实验】:OSPF基本配置(OSPF路由协议、影响DR/BDR选举、配置接口开销值、被动接口、认证功能)
  8. OpenGL程序:实例练习
  9. python之扑克牌小游戏
  10. Sambamba: process your BAM data faster!