扩散方程的三种导出方法

扩散方程稳态扩散与非稳态扩散

扩散方程稳态扩散与非稳态扩散 1．稳态扩散下的菲克第一定律（一定时间内，浓度不随时间变化dc/dt=0） 单位时间内通过垂直于扩散方向的单位截面积的扩散物质流量（扩散通量）与该面积处的浓度梯度成正比 即J=－D（dc/dx） 其中D：扩散系数，cm2/s，J：扩散通量，g/cm2〃s ，式中负号表明扩散通量的方向与浓度梯度方向相反。 可见，只要存在浓度梯度，就会引起原子的扩散。 x轴上两单位面积1和2，间距dx，面上原子浓度为C1、C2 则平面1到平面2上原子数n1=C1dx ，平面2到平面1上原子数n2=C2dx 若原子平均跳动频率f, dt时间内跳离平面1的原子数为 n1f〃dt 跳离平面2的原子数为n2fdt，但沿一个方向只有1/2的几率，则单位时间内两者的差值即扩散原子净流量。 令，则上式 2．扩散系数的测定：

其中一种方法可通过碳在γ-Fe中的扩散来测定纯Fe的空心园筒，心部通渗碳气氛，外部为脱碳气氛，在一定温度 下经过一定时间后，碳原子从内壁渗入，外壁渗出达到平衡，则为稳态扩散单位时单位面积中碳流量： A：圆筒总面积，r及L：园筒半径及长度，q：通过圆筒的碳量 则： 即： 则： q可通过炉内脱碳气体的增碳求得，再通过剥层法测出不同r处的碳含量，作出C-lnr曲线可求得D。 第一定律可用来处理扩散中浓度不因时间变化的问 3．菲克第二定律：解决溶质浓度随时间变化的情况，即dc/dt≠0

两个相距dx垂直x轴的平面组成的微体积，J1、J2为进入、流出两平面间的扩散通量，扩散中浓度变化为，则单元体积中溶质积累速率为 （Fick第一定律） （Fick第一定律） （即第二个面的扩散通量为第一个面注入的溶质与在这一段距离内溶质浓度变化引起的扩散通量之和） 若D不随浓度变化，则 故： 4．Fick第二定律的解：很复杂，只给出两个较简单但常见问题的解 a. 无限大物体中的扩散

一类反应扩散方程解的长时间行为

I 一类反应扩散方程解的长时间行为 摘 要 本文主要在一个有界光滑区域中讨论了一类带有齐次Dirichlet 边值条件的反应扩散方程解的长时间行为，其方程的形式如下： 其中 偏微分算子是一致抛物的， ，满足一定条件。 对于以上方程，我们首先定义了该方程的弱解，之后我们在有限维空间中构造了一系列该方程的近似解，并证明了在维数趋于无穷时，存在子列收敛于该方程的弱解。最后，我们利用先验估计得到了该方程弱解的存在唯一性。 在获得方程弱解的存在唯一性后，我们便能定义伴随方程的解半群，并由此研究伴随方程解半群的全局吸引子。 为了证明解半群在 中存在全局吸引子，我们证明 了伴随方程的解半群在 与中有界吸收集的存在性，并利用Sobolev 紧嵌入定理得到了全局吸引子的存在性。 关 键 词：反应扩散方程；Galerkin 方法；全局吸引子；弱解

II ABSTRACT In this thesis, we mainly consider the long-time behavior of solutions for the following reaction-diffusion equation with homogeneous Dirichlet boundary condition in a bounded smooth domain : where The partial differential operator is uniformly parabolic, and satisfies some additional assumptions. First of all, we give the definition of weak solutions, and then, we construct a sequence of approximate solution in a n dimension subspace and show that there exists a subsequence will convergent to a weak solution of this problem when n goes to infinite. Finally, we establish the existence and uniqueness of weak solution by some aprior estimates. With the help of the existence and the uniqueness of weak solutions, we define the solution semigroup associate with the problem and investigate the existence of a global attractor for the semigroup. To prove the existence of a global attractor, we show that there exist bounded absorbing sets in and and obtain existence of a global attractor in by using the Sobolev compactness embedding theorem. KEY WORDS: Reaction-diffusion equation; Galerkin’s method ; Global attractor; Weak solution

热传导方程及其定解问题的导出

第一章 热传导方程 本章介绍最典型的抛物型方程—热传导方程,在研究热传导,扩散等物理现象时都会遇 到这类方程. §1 热传导方程及其定解问题的导出 1.1热传导方程的导出 物理模型 在三维空间中,考虑一均匀,各向同性的物体Ω,假定它内部有热源,并且与周围介质有热交换,需要来研究物体内部温度的分布和变化. 以函数),,,(t z y x u 表示物体Ω在位置),,(z y x 及时刻t 的温度.物体内部由于各部分温度不同,产生热量的传递,它们遵循能量守恒定律. 能量守恒定律 物体内部的热量的增加等于通过物体的边界流入的热量与由物体内部的热源所生成的热量的总和 . 在物体Ω内任意截取一块D .现在时段],[21t t 上对D 使用能量守恒定律. 设),,,(t z y x u u =是温度(度),c 是比热(焦耳∕度·千克),ρ是密度(千克/米3), q 是热流密度(焦耳/秒·米2),0f 是热源强度(焦耳/千克·秒). 注意到在dt 时段内通过D 的边界D ?上小块dS 进入区域D 的热量为dSdt n q ?-(n 是 D ?的外法向),从而由能量守恒律,我们有 ,)||(21 21 120??????????+?-=-?==t t D t t D D t t t t dxdydz f dt ds n q dt dxdydz u u c ρρ （1.1） 大家知道,热量流动的原因是因为在物体内部存在温差.依据传热学中的傅立叶实验定律,在一定条件下,热流向量与温度梯度成正比 ,u k q ?-= （梯度? ?? ? ????????==?z u y u x u gradu u ,,） （1.2） 这里负号表明热量是由高温向低温流动,k 是物体的导热系数.

热传导方程的求解

应用物理软件训练 前言 MATLAB 是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。 MATLAB是矩阵实验室（Matrix Laboratory）的简称，和Mathematica、Maple 并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其

他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。本部分主要介绍如何根据所学热传导方程的理论知识进行MATLAB数值实现可视化。本部分主要介绍如何根据所学热传导方程的理论知识进行MATLAB数值实现可视化。本部分主要介绍如何根据所学热传导方程的理论知识进行MATLAB数值实现可视化。 本部分主要介绍如何根据所学热传导方程的理论知识进行MATLAB数值实现可视化。

题目：热传导方程的求解 目录 一、参数说明 (1) 二、基本原理 (1) 三、MATLAB程序流程图 (3) 四、源程序 (3) 五、程序调试情况 (6) 六、仿真中遇到的问题 (9) 七、结束语 (9) 八、参考文献 (10)

一、参数说明 U=zeros(21,101) 返回一个21*101的零矩阵 x=linspace(0,1,100);将变量设成列向量 meshz(u)绘制矩阵打的三维图 axis([0 21 0 1]);横坐标从0到21，纵坐标从0到1 eps是MATLAB默认的最小浮点数精度 [X,Y]=pol2cart(R,TH);效果和上一句相同 waterfall(RR,TT,wn)瀑布图 二、基本原理 1、一维热传导问题 （1）无限长细杆的热传导定解问题 利用傅里叶变换求得问题的解是： 取得初始温度分布如下 这是在区间0到1之间的高度为1的一个矩形脉冲，于是得 （2）有限长细杆的热传导定解问题

热传导方程

前言 本文只是针对小白而写，可以使新手对热传导理论由很浅到不浅的认识，如想更深学习热传导知识，请转其它文档。 一、概念与常量 1、温度场： 指某一时刻下，物体内各点的温度分布状态。 在直角坐标系中：； 在柱坐标系中：； 在球坐标系中：。 补充：根据温度场表达式，可分析出导热过程是几维、稳态或非稳态的现象，温度场是几维的、稳态的或非稳态的。 2、等温面与等温线： 三维物体内同一时刻所有温度相同的点的集合称为等温面； 一个平面与三维物体等温面相交所得的的曲线线条即为平面温度场中的等温线。 3、温度梯度： 在具有连续温度场的物体内，过任意一点P温度变化率最大的方向位于等温线的法线方向上。称过点P的最大温度变化率为温度梯度(temperature gradient)。用grad t表示。 定义为： 补充：温度梯度表明了温度在空间上的最大变化率及其方向，是向量，其正向与热流方向恰好相反。对于连续可导的温度场同样存在连续的温度梯度场。

在直角坐标系中： 3、导热系数 定义式：单位 导热系数在数值上等于单位温度降度(即1)下，在垂直于热流密度的单位面积上所传导的热流量。导热系数是表征物质导热能力强弱的一个物性参数。 补充：由物质的种类、性质、温度、压力、密度以及湿度影响。 二、热量传递的三种基本方式 热量传递共有三种基本方式：热传导；热对流；热辐射 三、导热微分方程式（统一形式：） 直角坐标系： 圆柱坐标系： 球坐标系： 其中，称为热扩散系数，单位，为物质密度，为物体比热容，为物体导热系数，为热源的发热率密度，为物体与外界的对流交换系数。 补充： 1处研究的对象为各向同性的、连续的、有内热源、物性参数已知的导热物体。 2稳态温度场，即则有：，此式称为泊松方程。 3无内热源的稳态温度场，则有：，此式称为拉普拉斯方程。 四、单值条件 导热问题的单值条件通常包括以下四项： 1几何条件：表示导热物体的几何形状与大小（一维、二维或三维）

传热基本方程及传热计算

第三节传热基本方程及传热计算 可知，要强化传热过程主要应着眼于增加推动力和减少热阻， 也就是设法增大 t m 或者 增大传热面积A 和传热系数K 。 在生产上，无论是选用或设计一个新的换热器还是对已有的换热器进行查定，都是建 立在上述基本方程的基础上的， 传热计算则主要解决基本方程中的 Q ，A,K, tm 及有关量的 计算。传热基本方程是传热章中最主要的方程式。 、传热速率Q 的计算 冷、热流体进行热交换时，当热损失忽略，则根据能量守恒原理，热流体放出热 量Qh ，必等于冷流体所吸收的热量 Qc ，即Qn Qc ，称之热量衡算式。 i.i. 无相变化时热负荷的计算 (1) ( 1)比热法 Q m h c ph T 1 T 2 m c C pc t 2 11 式中 Q ――热负荷或传热速率， J .S 1或W ； mh ， mc ――热、冷流体的质量流量， kg.s -1; Cpc,Cph ――冷、热流体的定压比热，取进出口流体温度的算术平均值下的比热， k J . (kg.k ) -1; T 1 ，T 2——热流体进、出口温度，K(° C ); t 1 ，t 2 —冷流体的进出口温度，K(° C )。 (2) 热焓法 Q m(l 1 I 2) (4 — 13) 式中 丨 1 ――物料始态的焓，k J .kg -1； I 2 ――物料终态的焓，k J .kg -1。 2 ?有相变化时热负荷计算 Q Gr (4—14) 式中 G ――发生相变化流体的质量流量， kg.s -1; r ---- 液体汽化(或蒸汽冷凝)潜热， k J .kg -1。 注意：在热负荷计算时，必须分清有相变化还是无相变化， 然后根据不同算式进行计算。 对蒸汽的冷凝、冷却过程的热负荷，要予以分别计算而后相加。 当要考虑热损失时，则有： 从传热基本方程 或 Q kA t m t Q m 1 kA 传热推动力 传热热阻 (4-11) (4-lla) (4-12)

热传导方程抛物型偏微分方程和基本知识

1. 热传导的基本概念 1.1温度场 一物体或系统内部，只要各点存在温度差，热就可以从高温点向低温点传导， 即产生热流。因此物体或系统内的温度分布情况决定着由热传导方式引起的传热速率（导热速率）。 温度场：在任一瞬间，物体或系统内各点的温度分布总和。 因此，温度场内任一点的温度为该点位置和时间的函数。 〖说明〗 若温度场内各点的温度随时间变化，此温度场为非稳态温度场，对应于非稳 态的导热状态。 若温度场内各点的温度不随时间变化，此温度场为稳态温度场，对应于稳态 的导热状态。 若物体内的温度仅沿一个坐标方向发生变化，且不随时间变化，此温度场为 一维稳态温度场。 1.2 等温面 在同一时刻，具有相同温度的各点组成的面称为等温面。因为在空间同一点不可能同时有两个不同的温度，所以温度不同的等温面不会相交。 1.3 温度梯度 从任一点起沿等温面移动，温度无变化，故无热量传递；而沿和等温面相交 的任一方向移动，温度发生变化，即有热量传递。温度随距离的变化程度沿法向最大。 温度梯度：相邻两等温面间温差△t与其距离△n之比的极限。 〖说明〗 温度梯度为向量，其正方向为温度增加的方向，与传热方向相反。 稳定的一维温度场，温度梯度可表示为：grad t = dt/dx

2. 热传导的基本定律——傅立叶定律 物体或系统内导热速率的产生，是由于存在温度梯度的结果，且热流方向和 温度降低的方向一致，即与负的温度梯度方向一致，后者称为温度降度。 傅立叶定律是用以确定在物体各点存在温度差时，因热传导而产生的导热速率大小的定律。 定义：通过等温面导热速率，与其等温面的面积及温度梯度成正比： q = dQ/ds = -λ·dT/dX 式中：q 是热通量（热流密度），W/m2 dQ是导热速率，W dS是等温表面的面积，m2 λ是比例系数，称为导热系数，W/m·℃ dT / dX 为垂直与等温面方向的温度梯度 “－”表示热流方向与温度梯度方向相反 3. 导热系数 将傅立叶定律整理，得导热系数定义式： λ= q/(dT/dX) 物理意义：导热系数在数值上等于单位温度梯度下的热通量。因此，导热系 数表征物体导热能力的大小，是物质的物性常数之一。其大小取决于物质的组成结构、状态、温度和压强等。 导热系数大小由实验测定，其数值随状态变化很大。 3.1 固体的导热系数 金属：35～420W/(m·℃)，非金属：0.2～3.0W/ (m·℃) 〖说明〗

第八章 热传导和扩散问题的傅里叶解

第八章 热传导方程的傅里叶解 第一节 热传导方程和扩散方程的建立 8.1.1 热传导方程的建立 推导热传导方程和前面弦振动所用的数学方法完全相用，不同之处在于具体的物理规律不同。这里用到的是热学方面的两个基本规律，即能量守恒和热传导的傅里叶实验定律。 热传导的傅里叶实验定律：设有一块连续的介质，选定一定的坐标系，并用(,,,)u x y z t 表示介质内空间坐标为的一点在t 时刻的温度。若沿x 方向有一定的温度差，在x 方向也就一定有热量的传递。从宏观上看，单位时间内通过垂直x 方向的单位面积的热量q 与温度的沿x 方向的空间变化率成正比，即 x u q k x ?=-? （8-1.1） q 称为热流密度，k 称为导热系数。公式中的负号表示热流的方向和温度变化的方向正好相 反，即热量由高温流向低温。 研究三维各向同性介质中的热传导，在介质中三个方向上存在温度差，则有 x u q k x ?=-?，y u q k y ?=-?，z u q k z ?=-? 或 q k u =-?r 即热流密度矢量q r 与温度梯度u ?成正比。 下面以一维均匀细杆为例，根据傅里叶实验定律和能量守恒定律推导介质中的热传导方程。 第一步，定变量。研究介质x 位置处在t 时刻的温度(,)u x t 。 第二步，取局部。在介质内部隔离出从x 到x x +?一段微元长度，在t 到t t +?时间内温度的变化(,)(,)u u x t t u x t ?=+?-。 第三步，立假设。假设均匀介质的横截面积为A ，质量密度为ρ，比热为c ，热传导系数为k 。

第四步，找规律。隔离出来的微元长度在t 到t t +?时间内吸收的热量为： Q c m u c A x u ρ=????=???? （8-1.2） 在t 到t t +?时间内，同过x 位置处的横截面的热量为： 1x x x Q q A t k u A t =???=-?? （8-1.3） 在t 到t t +?时间内，同过x x +?位置处的横截面的热量为： 2x x x x x Q q A t k u A t +?+?=???=-?? （8-1.4） 如果在微元段内有其他的热源，假设在单位时间单位体积内产生的热量为(,)F x t ，则该热源在微元内产生的热量为： (,)3Q F x t t A x =???? （8-1.5） 第五步，列方程。根据能量守恒定律，净流入的热量应该等于介质在此时间内温度升高所需要的热量。 123Q Q Q Q =-+ 即 (,)x x x x x c A x u k u A t k u A t F x t t A x ρ+?????=-??+??+???? (,)x x x x x u u u c k F x t t x ρ+?-??? =+?? 得到： (,)t xx k F x t u u c c ρρ = + 令 a = (,)(,)F x t f x t c ρ= 则得到热传导方程为 (,)2t xx u a u f x t =+ （8-1.6） 当介质内部无其他热源时，热传导方程是齐次的，为 2t xx u a u = （8-1.7） 8.1.2 扩散方程的建立

第四章导热问题的数值解法

第四章导热问题的数值解法 1 、重点内容：①掌握导热问题数值解法的基本思路； ②利用热平衡法和泰勒级数展开法建立节点的离散方程。 2 、掌握内容：数值解法的实质。 3 、了解内容：了解非稳态导热问题的两种差分格式及其稳定性。 §4—1导热问题数值求解的基本思想及内节点方程的建立由前述 3 可知，求解导热问题实际上就是对导热微分方程在定解条件下的积分求解，从而获得分析解。但是，对于工程中几何形状及定解条件比较复杂的导热问题，从数学上目前无法得出其分析解。随着计算机技术的迅速发展，对物理问题进行离散求解的数值方法发展得十分迅速，并得到广泛应用，并形成为传热学的一个分支——计算传热学（数值传热学），这些数值解法主要有以下几种： （１）有限差分法（ 2 ）有限元方法（ 3 ）边界元方法 数值解法能解决的问题原则上是一切导热问题，特别是分析解方法无法解决的问题。如：几何形状、边界条件复杂、物性不均、多维导热问题。 一．分析解法与数值解法的异同点： ?相同点：根本目的是相同的，即确定① t=f(x ， y ， z) ；② 。 ?不同点：数值解法求解的是区域或时间空间坐标系中离散点的温度分布代替连续的温度场；分析解法求解的是连续的温度场的分布特征，而不是分散点的数值。 数值求解的基本思路及稳态导热内节点离散方程的建立 二．解法的基本概念 ?实质 对物理问题进行数值解法的基本思路可以概括为：把原来在时间、空间坐标系中连续的物理量的场，如导热物体的温度场等，用有限个离散点上的值的集合来代替，通过求解按一定方法建立起来的关于这些值的代数方程，来获得离散点上被求物理量的值。该方法称为数值解法。 这些离散点上被求物理量值的集合称为该物理量的数值解。 2 、基本思路：数值解法的求解过程可用框图 4-1 表示。 由此可见： 1 ）物理模型简化成数学模型是基础； 2 ）建立节点离散方程是关键； 3 ）一般情况微分方程中，某一变量在某一坐标方向所需边界条件的个数等于该变量在该坐标方向最高阶导数的阶数。 ?数值求解的步骤 如图 4-2 （ a ），二维矩形域内无内热源、稳态、常物性的导热问题采用数值解法的步骤如下：（1）建立控制方程及定解条件 针对图示的导热问题，它的控制方程（即导热微分方程）为：（ a ） 边界条件： x=0 时， x=H 时， 当 y=0 时，

偏微分方程数值解法

《偏微分方程数值解法》 课程设计 题目：六点对称差分格式解热传导方程的初边值 问题 姓名：王晓霜 学院：理学院 专业：信息与计算科学 班级： 0911012 学号： 091101218 指导老师：翟方曼 2012年12月14日

一、题目 用六点对称差分格式计算如下热传导方程的初边值问题 222122,01,01(,0),01 (0,),(1,),01x t t u u x t t x u x e x u t e u t e t +???=<<<≤?????=≤≤??==≤≤??? 已知其精确解为 2(,)x t u x t e += 二、理论 1．考虑的问题 考虑一维模型热传导方程 (1.1) )(22x f x u a t u +??=??，T t ≤<0 其中a 为常数。)(x f 是给定的连续函数。（1.1）的定解问题分两类： 第一，初值问题（Cauchy 问题）：求足够光滑的函数()t x u ,，满足方程（1.1）和初始条件： （1.2） ()()x x u ?=0,, ∞<<∞-x 第二，初边值问题（也称混合问题）：求足够光滑的函数()t x u ,，满足方程（1.1）和初始条件： ()13.1 ()()x x u ?=0,, l x l <<- 及边值条件 ()23.1 ()()0,,0==t l u t u ， T t ≤≤0 假定()x f 和()x ?在相应的区域光滑，并且于()0,0，()0,l 两点满足相容条件,则上述问题有唯一的充分光滑的解。 现在考虑边值问题（1.1），（1.3）的差分逼近 取 N l h = 为空间步长，M T =τ为时间步长，其中N ，M 是自然数， jh x x j ==, ()N j ,,1,0Λ=； τk y y k ==, ()M k ,,1,0Λ=

传热基本方程及传热计算

第三节 传热基本方程及传热计算 从传热基本方程 m t kA Q ?= （４－１１） 或 传热热阻传热推动力= ?=kA t Q m 1 （４－１１ａ） 可知，要强化传热过程主要应着眼于增加推动力和减少热阻，也就是设法增大m t ?或者 增大传热面积Ａ和传热系数Ｋ。 在生产上，无论是选用或设计一个新的换热器还是对已有的换热器进行查定，都是建立在上述基本方程的基础上的，传热计算则主要解决基本方程中的m t K A Q ?,,,及有关量的 计算。传热基本方程是传热章中最主要的方程式。 一、传热速率Ｑ的计算 冷、热流体进行热交换时，当热损失忽略，则根据能量守恒原理，热流体放出热量 h Q ，必等于冷流体所吸收的热量c Q ，即c n Q Q =，称之热量衡算式。 1． 1． 无相变化时热负荷的计算 （1） （1） 比热法 () ()1221t t c m T T c m Q pc c ph h -=-= （4-12） 式中 Q ——热负荷或传热速率，J.s -1或W ； c h m m ,——热、冷流体的质量流量，kg.s -1； ph pc c c ,——冷、热流体的定压比热，取进出口流体温度的算术平均值下的比热， k Ｊ.（kg.k ）-1； 21,T T ——热流体进、出口温度，Ｋ（°C ）； 21,t t －冷流体的进出口温度，Ｋ（°C ）。 （２）热焓法 )(21I I m Q -= （４－１３） 式中 1I ——物料始态的焓，k Ｊ.kg -1； 2I ——物料终态的焓，k Ｊ.kg -1。 ２．有相变化时热负荷计算 Gr Q = （４－１４） 式中 G ——发生相变化流体的质量流量，kg.s -1； r ——液体汽化（或蒸汽冷凝）潜热，k Ｊ.kg -1。 注意：在热负荷计算时，必须分清有相变化还是无相变化，然后根据不同算式进行计算。对蒸汽的冷凝、冷却过程的热负荷，要予以分别计算而后相加。 当要考虑热损失时，则有：

热传导方程的导出及其定解问题的导出

热传导方程的导出及其定解问题的导出 1. 热传导方程的导出 考察空间某物体G 的热传导问题。以函数(,,,)u x y z t 表示物体G 在位置(,,)x y z 及时刻t 的温度。 依据传热学中的Fourier 实验定律，物体在无穷小时段dt 内沿法线方向n 流过一个无穷小面积dS 的热量dQ 与物体温度沿曲面dS 法线方向的方向导数 u n ??成正比，即 (,,) u dQ k x y z dSdt n ?=-? (1-1) 其中(,,)k x y z 称为物体在点(,,)x y z 处的热传导系数，它应取正值。(1-1)式中负号的出现是由于热量总是从温度高的一侧流向低的一侧，因此dQ 应和 u n ??异号。 在物体G 内任取一闭曲面Γ，它所包围的区域记为Ω，由(1-1)式，从时刻1t 到2t 流进此闭曲面的全部热量为 21(,,)t t u Q k x y z dS dt n Γ?? ?=??????? ? (1-2) 这里 u n ??表示u 沿Γ上单位外法线方向n 的方向导数。 流入的热量使物体内部的温度发生变化，在实践间隔12(,)t t 中物体温度从1(,,,)u x y z t 变化到2(,,,)u x y z t ，它所应该吸收的热量是 21(,,)(,,)[(,,,)(,,,)]c x y z x y z u x y z t u x y z t dxdydz ρΩ -??? 其中c 为比热，ρ为密度。因此就成立 21 21(,,) (,,)(,,)[(,,,)(,,,)]t t u k x y z dS dt c x y z x y z u x y z t u x y z t dxdydz n ρΓΩ ??? =-??????????? (1-3) 假设函数u 关于变量,,x y z 具有二阶连续偏导数，关于t 具有一阶连续偏导数，利用格林公式，可以把(1-3)化为 2 21 1t t t t u u u u k k k dxdydzdt c dt dxdydz x x y y z z t ρΩΩ????????????????? ++=?????????? ? ? ? ?????????????? ???? 交换积分次序，就得到 2 1 0t t u u u u c k k k dxdydzdt t x x y y z z ρΩ?? ?????????????---=?????? ? ? ??????????????? ? (1-4) 由于12,,t t Ω都是任意的，我们得到

偏微分方程数值解法的MATLAB源码

[原创]偏微分方程数值解法的MATLAB源码【更新完毕】 说明：由于偏微分的程序都比较长，比其他的算法稍复杂一些，所以另开一贴，专门上传偏微分的程序谢谢大家的支持！其他的数值算法见：..//Announce/Announce.asp?BoardID=209&id=82450041、古典显式格式求解抛物型偏微分方程（一维热传导方程） function [U x t]=PDEParabolicClassicalExplicit(uX,uT,phi,psi1,psi2,M,N,C) %古典显式格式求解抛物型偏微分方程 %[U x t]=PDEParabolicClassicalExplicit(uX,uT,phi,psi1,psi2,M,N,C) % %方程：u_t=C*u_xx 0 <= x <= uX,0 <= t <= uT %初值条件：u(x,0)=phi(x) %边值条件：u(0,t)=psi1(t), u(uX,t)=psi2(t) % %输出参数：U -解矩阵，第一行表示初值，第一列和最后一列表示边值，第二行表示第2层…… % x -空间变量 % t -时间变量 %输入参数：uX -空间变量x的取值上限 % uT -时间变量t的取值上限 % phi -初值条件，定义为内联函数 % psi1 -边值条件，定义为内联函数 % psi2 -边值条件，定义为内联函数 % M -沿x轴的等分区间数 % N -沿t轴的等分区间数 % C -系数，默认情况下C=1 % %应用举例： %uX=1;uT=0.2;M=15;N=100;C=1; %phi=inline('sin(pi*x)');psi1=inline('0');psi2=inline('0'); %[U x t]=PDEParabolicClassicalExplicit(uX,uT,phi,psi1,psi2,M,N,C); %设置参数C的默认值 if nargin==7 C=1; end %计算步长 dx=uX/M;%x的步长 dt=uT/N;%t的步长 x=(0:M)*dx; t=(0:N)*dt;

第一章导热理论和导热微分方程

第一章 导热理论和导热微分方程 相互接触的物体各部分之间依靠分子、原子和自由电子等微观粒子的热运动而传递热量的过程称为导热。在纯导热过程中物体各部分之间没有宏观运动。 与固体物理的理论研究方法不同，传热学研究导热问题时不是对导热过程的微观机理作深入的分析，而是从宏观的、现象的角度出发，以实验中总结出来的基本定律为基础进行数学的推导，以得到如温度分布、温度-时间响应和热流密度等有用的结果。这种处理方法的物理概念简单明了，但所要求的数学知识和技能仍是复杂和困难的。本书在材料的选取上，注意在介绍有重要应用价值的结果的同时，也给予求解导热问题的典型数学方法以足够的重视，以培养和发展读者独立解决问题的能力。 1-1 导热基本定律 1-1-1 温度场 由于传热学以宏观的、现象的方式来研究导热问题，团此必须引入连续介质假定，以便用连续函数来描述温度分布。温度场就是在一定的时间和空间域上的温度分布。它可以表示为空间坐标和时间的函数。由于温度是标量，温度场是标量场。常用的空间坐标系有三种：直角坐标系、柱坐标系和球坐标系。在直角坐标系中，温度场可以表示为 (,,,)t f x y z τ= （1-1-1） 式中：t 表示温度；x 、y 、z 为三个空间坐标；τ表示时间。 若温度场各点的温度均不随时间变化，即0t τ??=，则称该温度场为稳态温度场，否则为非稳态温度场。若温度场只是一个空间坐标的函数，则称为一维温度场；若温度场是两个或三个空间坐标的函数，则称为二维或三维温度场。 1-1-2 等温面与温度梯度 物体内温度相同的点的集合所构成的面叫做等温面。对应不同温度值的等温面构成等温面族。等温面与任一截面的交线形成等温线。由于等温线具有形象直观的优点，二维温度场常用等温线来表示温度分布。 由于在同一时刻物体的一个点上只能有一个温度值，所以不同的等温面不可能相交。它们或者在域内形成封闭曲线，或者终止于物体的边界。 如图1-l 所示，在物体内某一点P 处，沿空间某一方向l 的温度的变化率 图1-l 等温线和温度梯度