6 不可压缩牛顿流体流动 - 腾讯云

6 不可压缩牛顿流体流动. 2018-08-08 01:00:28阅读4880. 流体流动计算非常复杂，远远超出本教材。

最常见的SIMPLE算法编程实现也是比较复杂的，幸运的是MAC算法可以相对 ... 腾讯云备案控制台腾讯云开发者社区专栏问答沙龙团队主页TVP搜索搜索关闭创作写文章发视频提问登录注册展开腾讯云·社区登录首页专栏问答沙龙团队主页TVP返回腾讯云官网周星星9527156篇文章6不可压缩牛顿流体流动转到我的清单专栏首页javascript趣味编程6不可压缩牛顿流体流动00分享分享文章到朋友圈分享文章到QQ分享文章到微博复制文章链接到剪贴板海报分享海报分享6不可压缩牛顿流体流动发布于2018-08-0816:00:28阅读5000流体流动计算非常复杂，远远超出本教材。

最常见的SIMPLE算法编程实现也是比较复杂的，幸运的是MAC算法可以相对容易的实现一些简单算例。

Matlabfileexchange上一个顶驱方腔流动的例子，使用Matlab计算流体流动，代码如下：1.clc 2.clearvars 3.%parametersthatcanbemodified 4.nx=40;%numberofpoints 5.ny=40; 6.xmin=0;xmax=1;%domaindimentions 7.ymin=0;ymax=1; 8.dt=0.0005; 9.nstep=40000; 10.mu=0.00002; 11.un=1;us=0;ve=0;vw=0;%b.cfortangentialvelocity 12.maxit=850; 13.pit=150;%ploteach'pit'iterations 14.plotvorticity=0;%1toplotvorticity0tonotplot 15.%---------------------------------- 16.%--------initializing-------------- 17.%---------------------------------- 18.time=0.0; 19.dx=(xmax-xmin)/nx; 20.dy=(ymax-ymin)/ny; 21.[X,Y]=meshgrid(xmin:dx:xmax,ymin:dy:ymax); 22.u=zeros(nx+1,ny+2); 23.v=zeros(nx+2,ny+1); 24.p=zeros(nx+2,ny+2); 25.pp=zeros(nx+1,ny+1); 26.ut=zeros(nx+1,ny+2); 27.vt=zeros(nx+2,ny+1); 28.uu=zeros(nx+1,ny+1); 29.vv=zeros(nx+1,ny+1); 30.ifplotvorticity==1 31.w=zeros(nx+1,ny+1); 32.end 33. 34.foris=1:nstep 35.%enforcingb.cbyinterpolation 36.u(1:nx+1,1)=2*us-u(1:nx+1,2); 37.u(1:nx+1,ny+2)=2*un-u(1:nx+1,ny+1); 38.v(1,1:ny+1)=2*vw-v(2,1:ny+1); 39.v(nx+2,1:ny+1)=2*ve-v(nx+1,1:ny+1); 40. 41.i=2:nx;j=2:ny+1;%temporaryu-velocity 42.ut(i,j)=u(i,j)+dt*((-0.25)*... 43.(((u(i+1,j)+u(i,j)).^2-(u(i,j)+u(i-1,j)).^2)/dx...%DUUDX 44.+((u(i,j+1)+u(i,j)).*(v(i+1,j)+v(i,j))...%DUVDY1 45.-(u(i,j)+u(i,j-1)).*(v(i+1,j-1)+v(i,j-1)))/dy)...%DUVDY2 46.+(mu)*((u(i+1,j)-2*u(i,j)+u(i-1,j))/dx^2...%D2UDX 47.+(u(i,j+1)-2*u(i,j)+u(i,j-1))/dy^2));%D2UDY 48. 49.i=2:nx+1;j=2:ny;%temporaryv-velocity 50.vt(i,j)=v(i,j)+dt*((-0.25)*... 51.(((u(i,j+1)+u(i,j)).*(v(i+1,j)+v(i,j))...%DUVDX1 52.-(u(i-1,j+1)+u(i-1,j)).*(v(i,j)+v(i-1,j)))/dx...%DUVDX2 53.+((v(i,j+1)+v(i,j)).^2-(v(i,j)+v(i,j-1)).^2)/dy)...%DVVDY 54.+(mu)*((v(i+1,j)-2*v(i,j)+v(i-1,j))/dx^2...%D2VDX 55.+(v(i,j+1)+-2*v(i,j)+v(i,j-1))/dy^2));%D2VDY 56.pt=p; 57.forit=1:maxit%solveforpressure 58.%neumannb.cforpressure 59.pt(1,:)=pt(2,:); 60.pt(nx+2,:)=pt(nx+1,:); 61.pt(:,1)=pt(:,2); 62.pt(:,ny+2)=pt(:,ny+1); 63. 64.i=2:nx+1;j=2:ny+1; 65.pt(i,j)=(0.5/(dx^2+dy^2))*((dy^2)*(pt(i+1,j)+pt(i-1,j))... 66.+(dx^2)*(pt(i,j+1)+pt(i,j-1))... 67.-(dx*dy/dt)*(dy*(ut(i,j)-ut(i-1,j))... 68.+dx*(vt(i,j)-vt(i,j-1)))); 69. 70.Er=max(max(pt-p)); 71.p=pt; 72.ifEr<10^-5 73.break; 74.end 75.end 76.%correctthevelocity 77.u(2:nx,2:ny+1)=... 78.ut(2:nx,2:ny+1)-(dt/dx)*(p(3:nx+1,2:ny+1)-p(2:nx,2:ny+1)); 79.v(2:nx+1,2:ny)=... 80.vt(2:nx+1,2:ny)-(dt/dy)*(p(2:nx+1,3:ny+1)-p(2:nx+1,2:ny)); 81. 82.time=time+dt;%plottheresults 83.ifis==pit*ceil(is/pit) 84. 85.uu(1:nx+1,1:ny+1)=0.5*(u(1:nx+1,2:ny+2)+u(1:nx+1,1:ny+1)); 86.vv(1:nx+1,1:ny+1)=0.5*(v(2:nx+2,1:ny+1)+v(1:nx+1,1:ny+1)); 87.pp=0.5*(p(1:nx+1,1:ny+1)+p(2:nx+2,2:ny+2)); 88.holdoff,h=pcolor(X,Y,flipud(rot90(pp))); 89.set(h,'FaceAlpha',0.3,'EdgeAlpha',0); 90.holdon; 91.quiver(X,Y,flipud(rot90(uu)),flipud(rot90(vv)),2,'color',[000]); 92.title(time); 93.ifplotvorticity==1 94.w(1:nx+1,1:ny+1)=(u(1:nx+1,2:ny+2)-u(1:nx+1,1:ny+1))/(2*dx)... 95.+(v(1:nx+1,1:ny+1)-v(2:nx+2,1:ny+1))/(2*dy); 96.contour(X,Y,flipud(rot90(w)),10); 97.end 98.axisequal,pause(0.005),drawnow 99.end 100. 101.end复制计算结果如下：此部分内容为选读，仅供参考，有兴趣的读者可以将其移植到HTML5程序。

接下来的课程我们将进入能源与动力工程专业相关的设计类课程，例如换热器的设计、蒸发器、冷凝器、汽包锅炉水位控制等等。

文章分享自微信公众号：传输过程数值模拟学习笔记复制公众号名称本文参与腾讯云自媒体分享计划，欢迎热爱写作的你一起参与！原始发表时间：2018-07-28如有侵权，请联系[email protected]删除。

展开阅读全文其他举报点赞0分享登录后参与评论0条评论非牛顿流体不可压缩流动的混合高阶方法，具有类似动力的对流行为在这项工作中，我们设计并分析了一种混合高阶（HHO）离散化方法，用于非牛顿流体的不可压缩的流动，并具有类似动力的对流行为。

我们在粘度和对流规律的一般假设下工作，...用户8436237CFD初步（Matlab工具箱CFDTool试用）NavierStokes（纳维叶－斯托克斯）方程是流体力学中描述粘性牛顿流体的方程，是目前为止尚未被完全解决的方程，只有大约一百多个特解被解出来，是最复杂的方...云深无际使用javascript多快好省绘制简单CAD图纸[题引]：CAD（计算机辅助制图）是随着计算机技术发展而来的新技术，用于精确绘制。

图纸上承理论设计，下接生产制造，重要性不言而喻。

当前CAD软件种类繁多...周星星9527简单二维CAD图纸(使用javascript多快好省绘制)[题引]：君子性非异也，善假于物也。

javascript本身不能绘制CAD图纸，但借助第三方库就可以绘制一些简单的二维图纸了。

到被戏称为“全球最大的同性...周星星9527北大韦神出手，一天解决了六博士四个月没搞定的难题机器之心报道 编辑：泽南、蛋酱 韦东奕：太简单了，没必要要钱。

基础科学领域的研究总是让人感到晦涩难懂，但我们常常能从一些「小事」上看出其中意义。

5月6...机器之心Java学习历程之----提升篇（二）描述粘性不可压缩流体动量守恒的运动方程。

简称N-S方程。

因1821年由纳维和1845年由斯托克斯分别导出而得名。

是流体力学中描述粘性牛顿流体的方程，是尚未被完...用户5410712北大韦神出手，一天解决了六博士四个月没搞定的难题关注并星标 从此不迷路 计算机视觉研究院 公众号ID｜ComputerVisionGzq 学习群｜扫码在主页获取加入方式 计算机视觉研究院专栏 作者：E...计算机视觉研究院计算流体力学简介计算流体力学是利用高速计算机求解流体流动的偏微分方程组，主要研究内容是通过计算机和数值方法来求解流体力学的控制方程，对流体力学问题进行模拟和分析。

目的是为了更好...CAE之家综述：CFD的未来之路计算流体动力学（CFD）于近50年来兴起，是一门相对年轻的学科。

但远不如上个世纪80年代，过去的15年里，计算流体力学（CFD）的发展一直停滞不前。

CAE之家ProCAST有限元铸造模拟仿真软件CAEProCAST作为铸造模拟方案的先锋，是全球首屈一指的材料物理学数值模拟原型和制造流程供应商。

PROCAST成功的关键在于使用真实材料物理特性，能够进行更真实的...用户10047369为什么空调能制冷和滴水？儿子前两天问我，为什么空调可以制冷？空调出风口为什么会有水滴？对于这种求知欲，应该鼓励，当然不能乱讲，于是上网看了一些资料。

bisal多“维”优化——前端高并发策略的更深层思考（转载）说到优化，大家在收到“优化指标”任务的时候。

通常会做两件事情——分析“优化指标”对应的痛点、寻找解决痛点的技术方案并施行。

那这样是否就足够了呢？我的答案是否定的...ClearloveCAN-PINN:基于耦合自动数值微分法的快速物理信息神经网络在这项研究中，提出了一种新的物理信息神经网络（PINN）方法，通过泰勒级数展开将相邻支持点和自动微分（AD）耦合起来，以提高训练的效率和准确性。

在配置点处进行P...芽芽多“维”优化：前端高并发策略的更深层思考本文笔者从前端高并发优化这一具体点出发，逐步向大家阐述笔者在优化的“术”之上思维层面的一些思考。

希望能给各位带来共鸣和感悟。

WeTest质量开放平台团队如何应对访问量激增？前端高并发策略深层思考作 者Gary，腾讯web前端开发高级工程师商业转载请联系腾讯WeTest获得授权，非商业转载请注明出处。

WeTest导读一项指标的变好，总少不了...WeTest质量开放平台团队YOffleNet|YOLOV4基于嵌入式设备的轻量化改进设计最新的基于CNN的目标检测模型相当精确，但需要高性能GPU实时运行。

对于内存空间有限的嵌入式系统来说，它们在内存大小和速度方面依旧不是很好。

集智书童公众号深度学习准备「爆破」著名的欧拉方程250多年来，数学家一直试图「爆破」物理学中一些最重要的方程：那些描述流体如何流动的方程。

如果他们成功了，那么他们将会发现一种情况，在这种情况下，这些方程会被...脑机接口社区6/6阀门和调节器设计在液体火箭发动机系统中，将控制工质流动或调节工质参数的组件称为阀门（简称阀），其主要作用是保障发动机启动前的准备工作、过渡过程（启动、转级、关机）和调节过程等都...ClancyCC写给初学者的Tensorflow介绍Tensorflow是广泛使用的实现机器学习以及其它涉及大量数学运算的算法库之一。

Tensorflow由Google开发，是GitHub上最受欢迎的机器学习库之...云水木石更多文章周星星9527关注专栏文章156阅读量93.4K获赞307作者排名2412腾讯云原生专题云原生技术干货，业务实践落地。

云安全最佳实践-创作者计划火热征文中，发布文章赢千元好礼！立即查看腾讯云自媒体分享计划入驻腾讯云开发者社区，共享百万资源包。

立即入驻广告关闭社区专栏文章阅读清单互动问答技术沙龙技术视频团队主页腾讯云TI平台活动自媒体分享计划邀请作者入驻自荐上首页技术竞赛资源技术周刊社区标签开发者手册开发者实验室关于视频介绍社区规范免责声明联系我们友情链接归档问题归档专栏文章归档快讯文章归档关键词归档开发者手册归档开发者手册Section归档腾讯云开发者扫码关注腾讯云开发者领取腾讯云代金券热门产品域名注册云服务器区块链服务消息队列网络加速云数据库域名解析云存储视频直播热门推荐人脸识别腾讯会议企业云CDN加速视频通话图像分析MySQL数据库SSL证书语音识别更多推荐数据安全负载均衡短信文字识别云点播商标注册小程序开发网站监控数据迁移Copyright©2013-2022TencentCloud.AllRightsReserved.腾讯云版权所有京公网安备11010802017518粤B2-20090059-1扫描二维码扫码关注腾讯云开发者领取腾讯云代金券