扩散方程稳态扩散与非稳态扩散 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
数学物理学报2018,38A(4):750-769 h ttp://a cta m s.w ip m.a c.c n 一类具有变扩散系数的非局部反应-扩散方程解的 爆破分析* 赵元章马相如# (中国海洋大学数学科学学院山东青岛266100) 摘要：该文考虑了具有变扩散系数的反应-扩散方程D ir ic h le t初边值问题解的爆破现象.利 用辅助函数法和修正微分不等式技巧，对变扩散系数和非线性项给出适当的条件，以保证解整 体存在或有限时刻发生爆破，并在整体空间中（N>1)导出了爆破时间的界.同时，给出几 个应用举例. 关键词：反应-扩散方程；变扩散系数；爆破时间的界. M R(2010)主题分类：35K65; 35B30; 35B40 中图分类号：O175.29 文献标识码：A 文章编号：1003-3998(2018)04-750-20 1引言 我们考虑具有变扩散系数和非局部源项的反应--扩散方程 ut=d iv(a(x)V u(x, t))+f (u),(x, t)G Q x(0, t*),(U)给出齐次D ir ic h le t边界条件和初始条件 u(x, t)=0,(x, t)G dQ x(0, t*),(1.2) u(x, 0)=u〇 (x),:x G Q,(1.3)其中0c R n(N21)为具有光滑边界d n的有界区域，t*<+⑴表示可能发生爆破的时 间，反之t*=+⑴?变扩散系数a(x)为正的适当光滑函数，非线性项f(u)为非负连续函数并 满足非局部条件，比如，包含（u(x,t))p J^(u(x,t))qd x型非局部项，其中p + q>1?初值u〇(x) 为非负C1类函数且满足适当的相容性条件?因此，由经典拋物型理论知，问题（1.1)-(1.3) 存在唯一的非负局部解且充分光滑. 方程（1.1)出现在许多物理现象和生物种群理论.比如，热传导现象中温度，流体的流 动中浓度及某种生物种群密度的扩散等，见文献[1-3]及相关文献. 收稿日期：2017-07-11;修订日期：2017-12-11 E-mail: zhaoyz@https://www.doczj.com/doc/842125633.html,;xrmaouc@https://www.doczj.com/doc/842125633.html, *基金项目：山东省研究生创新计划项目（SDYY14127) Supported by Innovation Program for Graduates of Shandong Province (SDYY14127) **通讯作者
扩散方程的差分解法 在研究热传导过程、扩散过程、边界层现象时，我们常常遇到抛物型方程，这类方程中最典型、最简单的就是热传导方程。热传导方程中的自变量中包括时间t ，它是描述一种随时间变化的物理过程，即所谓不定常现象。这类问题的基本定解问题应是初值问题，即在初始时刻（t=0）时给定定解条件，求解t>0时的解。 本文主要运用有限差分法对一维扩散方程进行求解，并对差分解的适定性、相容性、收敛性及稳定性进行分析，同时与解析解进行对比。 1.扩散方程 一维扩散方程为： 22u u t x α??=?? （1） 式中，u 为因知量，α为扩散系数，x 为坐标，t 为时间。 其定解条件如下： 初始条件： (,0)() 0x u x f x L =≤≤ （2） 边界条件： 12(0,)() , (,)()u t f t u L t f t == （3） 一般假定函数()f x ，1()f t ，2()f t 满足连接条件，即1(0)(0) f f =，2()(0) f L f =。 2.有限差分法 有限差分法是数值计算解微分方程古老的方法之一，也是系统化地、数值地求解数学物理方法的方程。其控制方程中的导数用离散点上函数值的差商代替。 差分格式可以分为显格式和隐格式。所谓显格式是指在任一结点上因变量在新是时间层上的值可以通过之前的时间层上相邻结点变量的值显式解出来。由于这些层的变量值是已知的，当时间向前推进时，空间点上的新的变量值就只需逐点计算就行了，因此显格式计算起来比较省事。隐格式则是指任一结点上变量在新的时间层的值，不能通过之前的时间层上相邻结点的值显式解出来，它不仅与之前的时间层上的已知值有关，而且也与新时间层的相邻结点的变量值有关。因而一个差分方程常常包括几个相邻结点上的未知数，未知数的个数取决于格式的构成形式。为了解出这些未知数需要联立新的方程，而每引进一个新的方程往往又同时引进了新的未知数。因此，隐格式总是伴随着求解巨大的代数方程组。隐格式的主要缺点是计算工作量大，因而不如显格式计算得快，但这只是就时间步长一样的情况而言的。隐格式的主要优点是时间步长可以比显格式能够采用的最大步长大很多。显格式的时间步长受到稳定性条件的限制，而隐格式则几乎不受限制。 3.方程的离散 3.1 显格式 采用时间前差及第n 时间层的空间中心差，得一维扩散方程的显格式解： 111 2 2()n n n n n j j j j j u u u u u t x α ++---+=?? （4） 即 111(2) n n n n n j j j j j u u r u u u ++-=+-+ （5）