节点离散温度场有限差分(有限容积)程序入门之三:2D温度场显式迭代计算(暂不考虑潜热)...

本文主要是介绍节点离散温度场有限差分(有限容积)程序入门之三:2D温度场显式迭代计算(暂不考虑潜热)...,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

首先声明,我是一个菜鸟。一下文章中出现技术误导情况盖不负责

          我们首先分析度温场的求解吧,假设边界条件和初始条件经已设定。在贴码代之前,我们先谈谈这个类要需什么属性和行为:节点组数用于存储盘算变量、网格小大、度维定义、盘算函数,也就这么多了。如何盘算某节点的度温?盘算其东南西北方位相接节点对该节点的穿导热之和便可,读者这里可以虑考一下如何添加源相和对流换热进去。

package Soong.Solver
{public class TSolver{public  var Tlist:Vector.<Node>;private var xGridNum:uint = 1;//Number of Grid Allocated in X Directionprivate var yGridNum:uint = 1;//Number of Grid Allocated in X Directionpublic var dx:Number = 1;//Grid Size in X Directionpublic var dy:Number = 1;//Grid Size in Y Directionpublic var Sx:Number = 0;//Area of Heat Interface in X Directionpublic var Sy:Number = 0;//Area of Heat Interface in X Directionpublic var cellVol:Number = 0;//Volume of Control Volumepublic var Freezing:Boolean=false;//If Time to Freezepublic function TSolver(xGridNum:uint,yGridNum:uint,dx:Number,dy:Number){this.xGridNum = xGridNum;this.yGridNum = yGridNum;this.dx = dx;this.dy = dy;this.Sx = dy * 1;this.Sy = dx * 1;this.cellVol = dx * dy * 1;}public function Step(timeStep:Number):void{var col:uint = 0;var row:uint = 0;var node:Node = null;for (col = 2; col < xGridNum - 2; col++ ){for (row = 2; row < yGridNum-2; row++ ){node = Tlist[Index(col, row)] as Node;CalTnext(timeStep,node,col,row);node.T0=node.T;}}}public function CalTnext(timeStep:Number,node:Node,col:uint,row:uint):void{var conner:Boolean=false;var node_Adj:Node = null;var conductionHeat:Number = 0;//For Node on/in Connor or Edgevar SxFix:Number=1;//Area Fix Factor For Non-Interior Region in X Directionvar SyFix:Number=1;//Area Fix Factor For Non-Interior Region in Y Directionvar VolFix:Number=1;//Volume Fix Factor For Non-Interior Region in Y Directionif(((col==2)&&(row==2))||((col==2)&&(row==yGridNum-3))||((col==xGridNum-3)&&(row==2))||((col==xGridNum-3)&&(row==yGridNum-3))){SxFix=1/2.0;SyFix=1/2.0;conner=true;}if((col==2)||(col==xGridNum-3)){VolFix/=2;if(!conner){SyFix=1/2.0;}}if((row==2)||(row==yGridNum-3)){			VolFix/=2;if(!conner){SxFix=1/2.0;}}node_Adj = Tlist[Index(col+1, row)] as Node;conductionHeat+=node.eHeatExchangeFactor*(node_Adj.T0-node.T0)*Sx*SxFix;node_Adj = Tlist[Index(col-1, row)] as Node;conductionHeat+=node.wHeatExchangeFactor*(node_Adj.T0-node.T0)*Sx*SxFix;node_Adj = Tlist[Index(col, row+1)] as Node;conductionHeat+=node.nHeatExchangeFactor*(node_Adj.T0-node.T0)*Sy*SyFix;node_Adj = Tlist[Index(col, row - 1)] as Node;conductionHeat+=node.sHeatExchangeFactor*(node_Adj.T0-node.T0)*Sy*SyFix;var dT:Number = conductionHeat * timeStep;dT /= cellVol * VolFix * node.Rho * node.Cp;node.T = node.T0 + dT;}public function LatentHeatRelease(node:Node):void{}//Apply the Boundary Conditionpublic function ApplyBC():void{}private function Index(col:uint=0,row:uint=0):uint{return row * xGridNum + col;}}
}
每日一道理
因为自信,在呀呀学语时,我靠着纤嫩的双腿,迈出人生的第一步;因为自信,我一次次将第一名的奖状高高举起;因为自信,我毫不吝惜地剪掉飘逸的长发,在运动场上展现风采……感谢自信,它给了我一双翅膀,让我在电闪雷鸣中去飞翔,在风雨中去搏击人生!

    
        简略吧,要需意注的是不同置位的节点其传热面积以及控制体体积不尽相同,要需Fix一下。这里给出开端的盘算结果(代迭100s的结果)。前目笔者没有贴出全部码代,这时按照笔者供提的序程是法无运行的,读者想一想,还点缺什么?

    

      将其称对到得整个界面:

    

    有点子样了,这还不是终究的盘算结果,凝结潜热还没有虑考进去,后续会补上。另外,我们没有离散偏微分方程,但是我们的方法和离散偏微分方程归同途殊的。或许读者可以解理无限差分和无限容积的连续与区别了。

文章结束给大家分享下程序员的一些笑话语录: 一个合格的程序员是不会写出 诸如 “摧毁地球” 这样的程序的,他们会写一个函数叫 “摧毁行星”而把地球当一个参数传进去。

转载于:https://www.cnblogs.com/jiangu66/archive/2013/04/26/3045541.html

这篇关于节点离散温度场有限差分(有限容积)程序入门之三:2D温度场显式迭代计算(暂不考虑潜热)...的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/301947

相关文章

Spring Boot + MyBatis Plus 高效开发实战从入门到进阶优化(推荐)

《SpringBoot+MyBatisPlus高效开发实战从入门到进阶优化(推荐)》本文将详细介绍SpringBoot+MyBatisPlus的完整开发流程,并深入剖析分页查询、批量操作、动... 目录Spring Boot + MyBATis Plus 高效开发实战:从入门到进阶优化1. MyBatis

Mybatis从3.4.0版本到3.5.7版本的迭代方法实现

《Mybatis从3.4.0版本到3.5.7版本的迭代方法实现》本文主要介绍了Mybatis从3.4.0版本到3.5.7版本的迭代方法实现,包括主要的功能增强、不兼容的更改和修复的错误,具有一定的参考... 目录一、3.4.01、主要的功能增强2、selectCursor example3、不兼容的更改二、

如何用java对接微信小程序下单后的发货接口

《如何用java对接微信小程序下单后的发货接口》:本文主要介绍在微信小程序后台实现发货通知的步骤,包括获取Access_token、使用RestTemplate调用发货接口、处理AccessTok... 目录配置参数 调用代码获取Access_token调用发货的接口类注意点总结配置参数 首先需要获取Ac

Python FastAPI入门安装使用

《PythonFastAPI入门安装使用》FastAPI是一个现代、快速的PythonWeb框架,用于构建API,它基于Python3.6+的类型提示特性,使得代码更加简洁且易于绶护,这篇文章主要介... 目录第一节:FastAPI入门一、FastAPI框架介绍什么是ASGI服务(WSGI)二、FastAP

基于Python开发PDF转Doc格式小程序

《基于Python开发PDF转Doc格式小程序》这篇文章主要为大家详细介绍了如何基于Python开发PDF转Doc格式小程序,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 用python实现PDF转Doc格式小程序以下是一个使用Python实现PDF转DOC格式的GUI程序,采用T

将java程序打包成可执行文件的实现方式

《将java程序打包成可执行文件的实现方式》本文介绍了将Java程序打包成可执行文件的三种方法:手动打包(将编译后的代码及JRE运行环境一起打包),使用第三方打包工具(如Launch4j)和JDK自带... 目录1.问题提出2.如何将Java程序打包成可执行文件2.1将编译后的代码及jre运行环境一起打包2

Python如何计算两个不同类型列表的相似度

《Python如何计算两个不同类型列表的相似度》在编程中,经常需要比较两个列表的相似度,尤其是当这两个列表包含不同类型的元素时,下面小编就来讲讲如何使用Python计算两个不同类型列表的相似度吧... 目录摘要引言数字类型相似度欧几里得距离曼哈顿距离字符串类型相似度Levenshtein距离Jaccard相

在不同系统间迁移Python程序的方法与教程

《在不同系统间迁移Python程序的方法与教程》本文介绍了几种将Windows上编写的Python程序迁移到Linux服务器上的方法,包括使用虚拟环境和依赖冻结、容器化技术(如Docker)、使用An... 目录使用虚拟环境和依赖冻结1. 创建虚拟环境2. 冻结依赖使用容器化技术(如 docker)1. 创

C#使用yield关键字实现提升迭代性能与效率

《C#使用yield关键字实现提升迭代性能与效率》yield关键字在C#中简化了数据迭代的方式,实现了按需生成数据,自动维护迭代状态,本文主要来聊聊如何使用yield关键字实现提升迭代性能与效率,感兴... 目录前言传统迭代和yield迭代方式对比yield延迟加载按需获取数据yield break显式示迭

使用C#代码计算数学表达式实例

《使用C#代码计算数学表达式实例》这段文字主要讲述了如何使用C#语言来计算数学表达式,该程序通过使用Dictionary保存变量,定义了运算符优先级,并实现了EvaluateExpression方法来... 目录C#代码计算数学表达式该方法很长,因此我将分段描述下面的代码片段显示了下一步以下代码显示该方法如