温度场有限容积法程序入门之三: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一下。可以预见,如果将SxFix、SyFix,VolFix设置为Node类的成员变量,计算会更快。这里给出初步的计算结果(迭代100s的结果)。目前笔者没有贴出所有代码,这时按照笔者提供的程序是无法运行的,读者想想,还缺点什么?

      将其对称得到整个界面:

    有点样子了,这还不是最终的计算结果,凝固潜热还没有考虑进去,可以使用物理意义明确的温度回升法处理。另外,我们没有离散偏微分方程,但是我们的方法和离散偏微分方程殊途同归的。也许读者可以理解有限差分和有限容积的连续与区别了。

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



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

相关文章

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#代码计算数学表达式该方法很长,因此我将分段描述下面的代码片段显示了下一步以下代码显示该方法如

如何用Java结合经纬度位置计算目标点的日出日落时间详解

《如何用Java结合经纬度位置计算目标点的日出日落时间详解》这篇文章主详细讲解了如何基于目标点的经纬度计算日出日落时间,提供了在线API和Java库两种计算方法,并通过实际案例展示了其应用,需要的朋友... 目录前言一、应用示例1、天安门升旗时间2、湖南省日出日落信息二、Java日出日落计算1、在线API2