精华内容
下载资源
问答
  • 一维非稳态导热微分方程的数值求解matlab程序,《传热学》、《数值传热学》、《工程热力学》等课程上机作业。采用差分法、热平衡法建立的离散方程。
  • 非稳态导热问题由于有时间变量,其数值计算出现了一些新的特点。在非稳态导热微分方程中,与时间因素相关的非稳态项是温度对时间的一阶导数,这给差分离散带来...本文从个具体问题出发来研究非稳态导热问题的数值计算。
  • 用于一维平板的非稳态导热计算,其中包含三类边界条件,采用差分法将偏微分方程离散,计算得到平板内部温度分布随时间变化规律
  • 有关一维非稳态导热问题的数值解及相关的问题解,里面有很多关于这个问题的解法
  • 一维非稳态导热热传导Matlab程序,通过此可以解决偏微方程
  • 一维非稳态导热方程求解(附Matlab程序)》由会员分享,可在线阅读,更多相关《一维非稳态导热方程求解(附Matlab程序)(3页珍藏版)》请在人人文库网上搜索。1、使用差分方法求解下面的热传导方程 2 ( , )( , )(01,00....

    《一维非稳态导热方程求解(附Matlab程序)》由会员分享,可在线阅读,更多相关《一维非稳态导热方程求解(附Matlab程序)(3页珍藏版)》请在人人文库网上搜索。

    1、使用差分方法求解下面的热传导方程 2 ( , )( , )(01,00.2,1) txx T x ta Tx txta 初值条件: 2 ( ,0)44T xxx; 边值条件: (0, )0 (1, )0 Tt Tt ; 使用差分公式 1,1,2 22 (, )2 ( , )(, )2 ( , )() ijijijiji jij xxij T xh tT x tT xh tTTT Tx tO h hh ,1, ( ,)( , ) ( , )( ) ijiji ji j tij T x tkT x tTT T x tO k kk 上面两式带入原热传导方程 ,1,1,1, 2 2 i ji jiji。

    2、 jij TTTTT kh 令 2 2 4 k r h ,化简上式的 ,1,1,1, (1 2 )() i ji jijij Tr Tr TT 如下图: i x j t 1,ij rT , (1 2 ) i j r T 1,ij rT ,1i j T 编程 MATLAB 程序,运行结果如下 0 0.2 0.4 0.6 0.8 1 0 0.05 0.1 0.15 0.2 0 0.2 0.4 0.6 0.8 1 x t T function mypdesolution c=1; xspan=0 1; tspan=0 0.2; ngrid=100 10; f=(x)4*x-4*x.2; g1=(t)。

    3、0; g2=(t)0; T,x,t=rechuandao(c,f,g1,g2,xspan,tspan,ngrid); x,t=meshgrid(x,t); mesh(x,t,T); xlabel(x) ylabel(t) zlabel(T) function U,x,t=rechuandao(c,f,g1,g2,xspan,tspan,ngrid) % 热传导方程: % Ut(x,t)=c2*Uxx(x,t) axb tst0.5 error(为了保证算法的收敛,请增大步长 h 或减小步长 k!) end s=1-2*r; U=zeros(ngrid); % 边界条件 U(:,1)=g1(t); U(:,m)=g2(t); % 初值条件 U(1,:)=f(x); % 差分计算 for j=2:n for i=2:m-1 U(j,i)=s*U(j-1,i)+r*(U(j-1,i-1)+U(j-1,i+1); end end %本文来自互联网。

    展开全文
  • 使用差分方法求解下面的热传导方程2T (x, t) a T (x, t) (0 x 1, 0 t 0.2, a 1)t xx初值条件:T ( x,0) 4 x 4 x2 ;T (0, t) 0边值条件: ;T...

    使用差分方法求解下面的热传导方程

    2

    T (x, t) a T (x, t) (0 x 1, 0 t 0.2, a 1)

    t xx

    初值条件:T ( x,0) 4 x 4 x2 ;

    T (0, t) 0

    边值条件: ;

    T (1, t) 0

    使用差分公式

    T ( x h, t ) 2T ( x , t ) T ( x h, t ) T 2T T

    i j i j i j 2 i1, j i , j i1, j

    T ( x , t ) O(h ) 

    xx i j h2 h2

    T ( x , t k) T ( x , t ) T T

    T ( x , t ) i j i j O(k)  i, j1 i, j

    t i j k k

    上面两式带入原热传导方程

    T T T 2T T

    i, j1 i, j i1, j i, j i1, j

    k h2

    2

    4 k

    令r ,化简上式的

    h2

    T (12r)T r(T T )

    i, j1 i, j i1, j i1, j

    如下图:

    T

    i, j1

    rT (12r)T rT

    i1, j i, j i1, j

    t

    展开全文
  • 章熙民的第六版《传热学》里,较为简单的介绍了非稳态导热的数值计算,本文根据此书,以计算一个可视为无限大平壁的复合墙体传热过程为例,讨论一维非稳态导热问题数值求解的问题。 初始条件:一个由四层材料构成的...

    章熙民的第六版《传热学》里,较为简单的介绍了非稳态导热的数值计算,本文根据此书,以计算一个可视为无限大平壁的复合墙体传热过程为例,讨论一维非稳态导热问题数值求解的问题。

    初始条件:一个由四层材料构成的墙体(可视为无限大平壁),已知材料热工性质如下,该墙体的内表面换热系数为8.7W/(m^2*K),外表面换热系数为23W/(m^2*K),室内温度为20℃恒定不变,室外温度随昼夜变化。

    现给出一星期的室外温度数据(每小时),求墙体内部温度场。

    构造层材料干密度ρ-kg/m^3导热系数λ-W/(m*K)比热容c- J/(Kg*K)厚度mm

    热扩散率a=λ/(c*ρ)

    m^2/s

    1石膏板10500.361050123.2653E-07
    2钢筋混凝土25002.049201208.8696E-07
    3聚苯乙烯泡沫板200.04613801201.6667E-06
    4陶瓷红砖18000.78801204.4192E-07
    室外温度第1天第2天第3天第4天第5天第6天第7天
    0-25.4-24.5-20.5-19.8-21-23.5-15.5
    1-25.4-24.7-21.3-19.8-23.1-23.9-15.9
    2-25.4-25-22-19.8-25.1-24.4-16.4
    3-25.4-25.2-22.9-19.6-24.4-24.2-17.5
    4-25.4-25.5-23.8-19.4-23.6-24-18.6
    5-25.4-25.7-24.7-19.2-22.9-23.8-19.7
    6-25.4-26.4-25.1-19.2-23.4-23.2-18
    7-25-26.4-25-19-23.5-22.4-15.6
    8-24-25.4-24.2-18.5-22.9-21-13.2
    9-22-22.2-22.3-17.2-20.4-18.6-11.6
    10-19.8-18.9-20.2-15.7-17.6-15.8-10.3
    11-17.8-16.2-18.2-14.2-14.9-13.2-9.2
    12-16.5-15.6-16.9-13-13.3-11.2-8.1
    13-15.8-15.7-16.2-12.3-12.7-9.9-7.4
    14-15.9-16.3-16.2-12-12.8-9.5-7.2
    15-17-16.9-17.1-12.5-14-10.5-7.7
    16-18.4-17.7-18.3-13.3-15.2-12.1-8.6
    17-20-18.6-19.4-14.3-16.2-13.7-9.9
    18-21.3-19.4-19.8-15-16.2-13.9-11.5
    19-22-19.3-19.7-15.8-16.8-14.1-13.4
    20-23-19.6-19.8-16.5-17.1-14.4-15.1
    21-23.4-19.7-19.8-17.3-19.1-14.6-15.8
    22-23.8-19.7-19.8-18.2-21-14.8-16.4
    23-24.2-19.8-19.8-19-23-15-17.1
    • 前期准备

    在《传热学》第四章第三节中,对什么是非稳态导热的给出了如下解释:

    因此根据上述场景将该问题提炼为:一个两边是第三类边界条件的无限大平壁的一维非稳态导热问题。这里我们选用显示离散格式(求解上更为简单,但是对步长要求较为严苛)。

    内部节点方程式:

    稳定条件:

    边界节点方程式(第三类边界条件):

     稳定条件:

    •  问题分析

    在经过多次试验之后,将该墙体按间距3mm划分为132层,取时间步长为1s,此时满足稳定条件。由于这里是复合墙体,因此后续计算时也要注意节点的热物性参数会相应改变。

    还有一个重要的问题是墙体的初始温度场,显式离散格式的逻辑就是从初始温度出发逐个计算出各个时刻的温度。这里没有给,所以就用初始温度条件(室内20℃,室外-25.4℃)计算这一条件下的稳态温度。

     根据上式求得墙体内初始温度场。

    •  Matlab求解

    在完成上述工作之后,便可以在Matlab中求解了。

    clc;
    clf;
    format short e;
    global N delta_x delta_t TN;
    
    h1 = 8.7; %墙体内表面传热系数 W/(m^2*K)
    h4 = 23; %墙体外表面传热系数 W/(m^2*K)
    D = 0.332; %墙体厚度/m
    N = 124; %厚度划分数目
    delta_x = 0.003; %划分间距/m
    delta_t = 1; %时间步长/s
    day = 7; %计算天数
    TN = 3600*24*day; %计算时间节点数目
    Tn = 20; %室内温度
    T = zeros(TN,N+1); %预设数组大小
    Tf = zeros(1,TN);
    
    %设置室外空气温度 这一步是因为给的是每小时的温度 换成每秒均等增加的温度值
    for hour = 1:24*day %取值范围
        %将第hour小时和hour+1小时的温度等间距分成3600个
        %自定义数组outdoor_air为室外温度
        temp = linspace(outdoor_air(hour),outdoor_air(hour+1),3600); 
        for t = 1:3600
            Tf(3600*(hour-1)+t) = temp(t);%每秒均匀变化的温度
        end
    end
    
    %设置初始墙体内温度
    for I = 1:N+1
        %自定义数组wall为墙体初始温度
        T(1,I) = wall(I);
    end
    
    %第一层:石膏板
    a1 = 3.265e-7; %热扩散率 m^2/s
    lamb1 = 0.36; %热导率 W/(m*K)
    Fo1 = a1*delta_t/(delta_x)^2;
    Bi1 = h1*delta_x/lamb1;
    c1 = Fo1 - 1/(2*Bi1+2);
    
    %第二层:钢筋混凝土
    a2 = 8.870e-7; %热扩散率 m^2/s
    lamb2 = 2.04; %热导率 W/(m*K)
    Fo2 = a2*delta_t/(delta_x)^2;
    
    %第三层:聚苯乙烯泡沫板
    a3 = 1.667e-6; %热扩散率 m^2/s
    lamb3 = 0.046; %热导率 W/(m*K)
    Fo3 = a3*delta_t/(delta_x)^2;
    
    %第四层:陶瓷红砖
    a4 = 4.419e-7; %热扩散率 m^2/s
    lamb4 = 0.7; %热导率 W/(m*K)
    Fo4 = a4*delta_t/(delta_x)^2;
    Bi4 = h4*delta_x/lamb4;
    c4 = Fo4 - 1/(2*Bi4+2);
    
    while (Fo1>0.5 || Fo2>0.5 || Fo3>0.5 || Fo4>0.5 || c4>0  || c1>0)
        error('无法收敛');
    end

     第一部分主要是输入参数,设置初始条件,判断选择的步长是否满足稳定性条件。

    for t = 1:TN
        for i = 1:N+1      
            if i == 1
                T(t+1,i) = 2*Fo1*(T(t,i+1) + Bi1*Tn) + (1 - 2*Bi1*Fo1 - 2*Fo1)*T(t,i);
            end
            
            if (i>=2 && i<= 5)
                T(t+1,i) =  Fo1*(T(t,i-1) + T(t,i+1)) + (1 - 2*Fo1)*T(t,i);
            end
            
            if (i>=6 && i<= 45)
                T(t+1,i) =  Fo2*(T(t,i-1) + T(t,i+1)) + (1 - 2*Fo2)*T(t,i);
            end
            
            if (i>=46 && i<= 85)
                T(t+1,i) =  Fo3*(T(t,i-1) + T(t,i+1)) + (1 - 2*Fo3)*T(t,i);
            end
            
            if (i>=86 && i<= 124)
                T(t+1,i) =  Fo4*(T(t,i-1) + T(t,i+1)) + (1 - 2*Fo4)*T(t,i);
            end
            
            if i == N+1
                T(t+1,i) = 2*Fo4*(T(t,i-1) + Bi4*Tf(t)) + (1 - 2*Bi4*Fo4 - 2*Fo4)*T(t,i);
                %这里给了室外空气温度Tf一个时间函数关系
            end
        end
    end

     第二部分就是计算显式离散节点方程了,没有什么好说的,就用了最简单的迭代法,甚至不需要求解方程组,唯一需要注意的就是各个材料层的节点方程式使用的参数不一样。

    X = zeros(1,24*day*2);
    Y = zeros(1,N+1);
    TT = zeros(24*day*2,N+1);
    
    J = 1;
    for I = 1:24*day*2
        X(I) = 1800*I; %30分钟一个数
        for J = 1:N+1
            Y(J) = J;
            TT(I,J) = T(X(I),Y(J));
        end
    end
     
    figure(1)
    hold on;
     
    % Make a contour with the presence of isolines & texts
    [XX,YY] = meshgrid(Y,X);
    surfc(YY,XX,TT)
    
    xlabel('Time - AXIS')
    ylabel('Thickness - AXIS')
    zlabel('T (^{\circ}C)')
    title(colorbar,'^{\circ}C')

     第三部分画图,半个小时取一个数,因为一秒钟一个温度数太多了没有必要,画图也不好看。

    最后呈现出来的温度-时间-节点的图形就是这样:

    展开全文
  • 一维非稳态导热求解,绝热和热流边界条件,附Matlab程序,共同学习提高function oned% 根据网上一维算例改了下边界条件,边界条件为绝热和恒定热流 % by hxgclear all;clc;%%%%%%%%%%%%%%需要输入的物性参数Lambda=10;%...

    一维非稳态导热求解,绝热和热流边界条件,附Matlab程序,共同学习提高

    function oned

    % 根据网上一维算例改了下边界条件,边界条件为绝热和恒定热流 % by hxg

    clear all;clc;

    %%%%%%%%%%%%%

    %需要输入的物性参数

    Lambda=10;%导热系数

    cp=440;%热容

    rou=7800;%密度

    qw=500000;%壁面热流

    %%%%%%%%%%%%%

    a=Lambda/rou/cp;%定义中间系数

    c=qw/Lambda;%定义中间系数

    xspan=[0 0.012];%轴向坐标起止位置

    tspan=[0 10];%仿真时间起止

    ngrid=[1000 20];%空间网格数和时间网格数

    %%%%%%%%%%%%%

    %调用子函数

    [T,x,t]=rechuandao(a,c,xspan,tspan,ngrid);

    %画图

    [x,t]=meshgrid(x,t);

    figure(1)

    mesh(x,t,T);

    xlabel('x')

    ylabel('t')

    zlabel('T')

    function [U,x,t]=rechuandao(a,c,xspan,tspan,ngrid)

    % 热传导方程:

    % Ut(x,t)=c^2*Uxx(x,t) a

    %一维非稳态导热方程

    % 参数说明

    % c:方程中的系数

    % f:初值条件

    % g1,g2:边值条件

    % xspan=[a,b]:x的取值范围

    % tspan=[ts,tf]:t的取值范围

    % ngrid=[n,m]:网格数量,m为x网格点数量,n为t的网格点数量 % U:方程的数值解

    % x,t:x和t的网格点

    n=ngrid(1);

    m=ngrid(2);

    h=range(xspan)/(m-1);

    x=linspace(xspan(1),xspan(2),m);

    展开全文
  • 《传热学上机C程序源答案之一维非稳态导热的数值计算》由会员分享,可在线阅读,更多相关《传热学上机C程序源答案之一维非稳态导热的数值计算(4页珍藏版)》请在人人文库网上搜索。1、二维稳态导热的数值计算2.1物理...
  • 均匀杆的一维非稳态热传导数值求解
  • 一维稳态导热微分方程的数值求解,《传热学》、《数值传热学》、《工程热力学》等课程上机作业。采用差分法、热平衡法建立离散方程
  • 一维非稳态扩散算例一样的初始和边界条件 ,采用乘方格式 时间步长为0.001s,初始温度场为200,速度为2m/s,长度为2cm,t=0s时刻东侧温度突然降至0C。时间差分采用全隐式格式
  • 对平壁在导热第三类边界下的一维非稳态导热过程中的灯用传递进行计算,并将导温系数在定值与变值下的计算结果加以比较,对非稳态条件下灯用传递的时空分布规律进行分析。结果表明,灯用流密度随时间按指数规律递减,...
  • 一维非稳态导热通用程序(不变部分)

    千次阅读 2016-10-11 17:36:59
    一维非稳态导热通用程序(不变部分) ! This is a general purpose program to solve 1-D diffusion. ! problem in the form of: ! ρcdt/dz=1/a(x)d/dx(a(x)Γdt/dx)+s !***********
  • 为了研究非稳态导热问题中的常用有限体积方案(FVM1)和新提出的有限体积方案(FVM2)的精度差异问题,采用理论推导和实例验证的方法,建立基于有限体积法原理的热平衡积分方程,利用矩形网格进行离散化,建立了平面非稳态热...
  •  二维稳态导热的数值计算主要采用了热平衡法。用差分法建立节点的热平衡方程, 将节 点所在的单元体的四个方向传递的热流密度,内热源在单元体产生的热流密度, 根据能量守 恒的原则建立方程,可以得到每个节点...
  • 种二维非稳态导热问题的数值解法,可作为计算机数值分析的参考
  • 一维非稳态导热求解,绝热和热流边界条件,附Matlab程序,共同学习提高function oned% 根据网上一维算例改了下边界条件,边界条件为绝热和恒定热流 % by hxgclear all;clc;%%%%%%%%%%%%%%需要输入的物性参数Lambda=10;%...
  • 数值传热学非稳态导热,采用显式、隐式、C-N格式求解
  • 种二维非稳态导热问题的数值解法维普资讯第 9卷 第 2f/I 石 油 化 工 高 等 学 校 学 报 Vo1. 9No.21996~ 6月 JOURNAL OF PETROCHEM ICAL UNIVERSITIESOF SINOPEC J衄. 1996、a0/ 种二维非稳态导热问题的...
  • 传热学基本概念三维非稳态导热微分方程导热微分方程计算机在材料科学与工程中的应用 叶卫平 本 章 要 点 4.1材料学主要物理场 4.1材料学主要物理场 4.2 传热学基本概念 4.2 传热学基本概念 4.2 传热学基本概念 4.2 ...
  • matlab一维导热模型

    2018-12-18 18:48:44
    matlab一维导热模型,适合初学者,非常有借鉴意义,是个好东西
  • 维非稳态导热

    千次阅读 2020-03-17 18:27:41
    [外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-Uv8DxZfk-1584440806238)(geo.png)]{#img width=".4\textwidth"} 数学物理模型 物理模型 二有内热源的稳态导热模型。 数学模型 如图1所...
  • 复合圆筒壁导热问题有广泛的工程应用背景。例如,注汽和采油生产中的井筒由隔热油管、环空、套管、水泥环和地层等构成;地面油气输送的多层保温管道等也属复合圆筒壁。这些结构边界条件复杂,分析求解难度大甚至无法...
  • §4-3非稳态导热问题的数值解法由前可知:非稳态导热和稳态导热...一、一维非稳态导热时间——空间区域的离散化1 、基本概念如图 4-8 所示, x 为空间坐标, τ 为时间坐标。1 )时间步长:指从一个时间层到下一个时...
  • 一维非稳态非线性导热数值解。 有限差分法离散求解区域,TDMA法求解离散方程组。
  • matlab解稳态导热问题

    2021-04-20 01:34:12
    有限差分和Matlabpde求解...计算传热学程序报告 姓名: 题目:一维非稳态导热问题的数值解学号: 学院:能源与动力工程学院 专业:工程热物理 日期:2014 年 5 月 25 日 一维非稳态导热问题数值......算例表明, MATLAB 对...
  • 一维传热方程

    2021-04-23 19:47:49
    精品文档 第五次作业(前三题写在作业纸上) 一、用有限差分方法求解一维非定常热传导方程,初始条件和边界条件见说明.pdf 文件, 热扩散系数α =const, ?T ? ...传热基本方程及传热计算_物理_自然科学_专业资料。. 第三...

空空如也

空空如也

1 2 3 4 5 ... 11
收藏数 205
精华内容 82
关键字:

一维非稳态导热

友情链接: 55593383Excel.RAR