解四阶抛物型方程高精度两层显格式
四阶抛物方程——一类新的交替分组显格式
四阶抛物方程——一类新的交替分组显格式
郭阁阳;陈超
【期刊名称】《天津职业技术师范大学学报》
【年(卷),期】2007(017)004
【摘要】给出了逼近四阶抛物方程一组新的Saul'yev非对称差分格式,利用这组非对称格式构造了一类新的交替分组显格式,并证明了该算法的绝对稳定性.数值实验表明,该格式具有良好的收敛性、较高的误差精度和绝对稳定性.
【总页数】4页(P32-35)
【作者】郭阁阳;陈超
【作者单位】天津工程师范学院数理与信息科学系,天津,300222;天津工程师范学院数理与信息科学系,天津,300222
【正文语种】中文
【中图分类】O246
【相关文献】
1.四阶抛物型方程的一个新的高精度显格式 [J], 陈世平
2.四阶抛物方程一类新的并行交替分段隐格式 [J], 郭阁阳;刘播
3.一类变系数抛物方程的交替分段显-隐差分格式 [J], 郝涛;王怡
4.色散方程的一类新的高精度交替分组显隐算法 [J], 张青洁;王文洽
5.四阶抛物方程的一类交替分组方法 [J], 冯青华
因版权原因,仅展示原文概要,查看原文内容请购买。
解四阶抛物型方程高精度紧致差分格式
7 2
大 学 数 学
第2 6卷
其 中 D , q 依次 为关 于 的一 阶偏 微分算 子 , 移算 子与一 阶 中心差分 算子 , 面建 立 中心差 分 T - 位 下 算子 和微分 算子 D 的关 系式. T yo 展开 , 得 由 a lr 可
一
+ 1'? : ,)
了一个 三层 隐式 差分 格式 , 是格 式 的精度 比较 低 ; [ ] 造 了一个 三层 显式 差分 格 式 , 稳定 性条 件 但 文 4构 其
和局 部截 断误 差 阶分别 为 f f 1 8和 o(Z f。 Z )) 文 [ ] < / r ( ) +(X ; X 5 构造 了一个 两层 恒稳 隐 式格 式 和 一
因此 , 文针对 四阶抛物 型方程 ( ) 本 1 的周期 初值 问题 , 造 出了一 个两 层 高精度 紧致 差分 格 式和一 个 构 三层高精 度紧致 隐格 式 , 其截 断误差 阶分别为 O(△£ + ( z 和 o(a£ +( )zz △ . ( ) z )) 5 ( ) S △£( ) +( )) X
一c<z × 0 ≤T × <C,≤f , 3 。
…
1 (2 “ 3+L,) t 一“( £ , 一 ∞ < < ∞ , ≤ £ T, T,) O ≤ 一o %x o. o % o
一
对 于这 类 四阶抛 物 型方程 的数值解 求 解 , a l e S u ’v在 文 [ ] 出 了一 类 含 权 因子 a的两 层 差 分 格 式 , 1提 当 a 一0时 为显 式格式 , 其稳定 性 条件 为 f f 1 2一 文 [ ] 造 了一族 三层 ( 殊 情况 下 为两 层 ) 含双 参 < / 。 ; 2构 r 特 、 数、 绝对 稳定 、 精度 、 对角 线型 的 隐式差 分 格 式 , 局 部 截 断误 差 为 O( z +( ) ) At△z分 高 五 其 ( ) 5 z z , , 5 别 为时 间及 空间 步长 ; 后 , 随 曾文平 针对 四阶抛 物型 方程 提 出 了一系列 的差分 格式 ]其 中文 [ ] 造 , 3构
解高维抛物型方程的一个高精度显式差分格式
V0 1 . 3 1 N o . 2
Ma r . 2 O1 3
文章编号 : 1 0 0 8—1 4 0 2( 2 0 1 3 ) D 2— 0 2 7 5一 o 4
解 高维 抛 物 型 方 程 的 一个 高精 度 显 式 差分 格式①
沈高峰 , 谷淑敏
( 1 . 郑州轻工业学院计算 机与通 信工程学院 . 河南 郑州 4 5 0 0 0 2; 2 . 中原工学院信 息商务学院基础部 。 河南 郑州 4 5 0 0 0 7 )
摘 要: 构造和研 究了五雏抛物型方程的高精度显式差分格式. 首先给 出了含参 变量的差分方 程, 并用待定 系数法适 当地选取 了这些参数的表 示式, 以使差分方程的截 断误差阶尽可能高地达 到了o ( z  ̄ t +△ ); 其 次 用稳定 性分 析 的 F o u r i e r 方 法给 出了所得 格 式 的稳 定 性 条件 ; 接 着确 定 了高精度 显 式差分 格 式的稳 定性 条件 为 r<2 / 5; 最后 给 出了数 值例 子 , 数值 结果 表 明 了本文格
件 r<1 / 6 或 r<1 / 4 . 本 文对 五 维热 导方 程构造 了
一
△ [ ‰ +( 叩 +叩 : ( 1 ◇)+7 7 3口
+7 7 田
+ o) ]=( + + + + )
( 7 7 6 +7 7 7 n 胁 - 1 )+[ ( 1 2 ◇ +3 ( 口
第3 1 卷 第 2期
2 0 1 3 年 O 3月
佳 木 斯 大 学 学 报 ( 自 然 科 学 版 ) J o u na r l o f J i a mu s i U n i v e r s i t y( N a t u r a l S c i e n c e E d i t i o n )
抛物型方程的一族双参数高精度恒稳格式
第 2卷 3
第 4期
华 侨 大 学 学 报 (自 然 科 学 版 )
J u na o a i o Un v r iy ( t r lS i n e o r l fHu q a i e st Na u a ce c )
V o1 .23 N o. 4
考 虑 抛 物型 方 程周 期初 值 问题
收 稿 日期 2 0 —60 0 20 —9 作者简 介 曾 文 平 ( 9 0 ) 男 , 授 14一, 教
基 金 项 目 国务 院 侨 务 办 公 室 自然 科 学 基 金 资 助 项 目
维普资讯
中图分类号
众所 周 知 , 人们 对 抛物 型 方程
a
一
“
的 差 分 解 法 的 研 究 由 来 已久 .常 见 的差 分 格 式 “ , 如 L ao e ,C a k Ni os n和 Du 诸 as n n rn — c lo h
fr— rn e 等格 式 , 都是 绝对 稳定 的 , otF a k l 虽 但它 们 截断误 差 阶均 较低 . 两 种分 别 是 O( ) 前 ( +
』
8  ̄q xP t
一
ap q z +z
、 ( ’ 一 1, … ) ” 2, .
( 2)
议 时 1 步 长 为 At 至 1J 长 为 △ 稻 域 由 点 集 ( t 八 一 0 1 2, , n一 0 1 … ) 组 日 J , 日步 , z , , , … J; , , 所
一
寺, 一0时, 得到一个两层高精度 的差分格式 . 同时, 证明了该族差分格式 , 对任意选取 的
非负 实参 数 n 0 ≥ O都是绝 对 稳定 的 .数 值例 子 表明 , ≥ , 本文 所提 格 式是 有效 的 , 论分 析 与 理 实际计 算 相 吻合 .
《抛物型方程的高精度时空有限体积元方法》范文
《抛物型方程的高精度时空有限体积元方法》篇一一、引言抛物型方程是一类重要的偏微分方程,广泛应用于物理学、工程学和金融学等多个领域。
为了解决这类方程的数值解问题,本文提出了一种高精度的时空有限体积元方法。
该方法在时间和空间上均采用离散化处理,能够有效地捕捉到抛物型方程的动态变化过程,并提高解的精度。
二、抛物型方程的描述抛物型方程通常描述了热量传导、扩散等现象。
其基本形式为:u_t = a u_xx + f(x,t)其中,u表示因变量,t表示时间,x表示空间坐标,a为扩散系数,f(x,t)为源项。
三、时空有限体积元方法本节将详细介绍抛物型方程的时空有限体积元方法。
1. 空间离散化处理空间离散化是将连续的空间划分为有限个离散的空间单元。
在每个空间单元上,抛物型方程的解可以近似表示为该空间单元的平均值。
通过对空间单元进行适当的剖分,可以得到一个离散的空间网格。
2. 时间离散化处理时间离散化是将连续的时间划分为有限个离散的时间步长。
在每个时间步长内,可以采用合适的数值方法来近似求解抛物型方程。
为了获得较高的解精度,本方法采用高阶的时间离散化技术。
3. 有限体积元的构建在空间和时间离散化处理的基础上,可以构建有限体积元。
每个体积元都包含一定的空间和时间范围,可以用于近似求解抛物型方程。
通过合理选择体积元的形状和大小,可以有效地提高解的精度。
四、高精度求解策略为了提高解的精度,本文采用以下策略:1. 采用高阶的空间离散化技术,以减小空间误差;2. 采用高阶的时间离散化技术,以减小时间误差;3. 优化有限体积元的构建过程,以提高近似解的精度;4. 采用迭代法或自适应网格法等数值优化技术,进一步提高解的精度。
五、算法实现与结果分析本节将通过具体的数值实验来验证所提方法的有效性。
首先,给出具体的算法实现步骤;然后,通过与其它数值方法进行比较,分析所提方法的优越性;最后,给出具体的数值结果,并进行分析和讨论。
六、结论本文提出了一种高精度的时空有限体积元方法来求解抛物型方程。
四阶抛物型方程基于子域精细积分方法的五次非多项式样条解
n r cib e a d pat a l. c
Ke r s ou r e a a oi q a o y wo d :f r o rp r b l e u t n;s b d ma n p e ie i tg ai n o - o y o a u ni p n ;sa i t n lss r r d c i u ・ o i rc s ne r t ;n n p l n mi q i t s h e tb l y a ay i ;e r o l c i o
解四阶抛物型方程新的高精度显式差分格式
0 引言
考虑下列四阶抛物型方程初边值问题
维普资讯
第 3期
单 双荣 等 :解 四阶抛 物型方程新 的高精 度显式差分格 式
・2 9・ 7
微 分方 程 ( ) 具有 尽 可能 高 阶的离 散误 差 ,而且 有 较好 的稳 定性 . 1
表 示 方 向 的 四阶 中心 差 分算
子,即:
0 u 3
一 一
+++) 如 弛窘+ ]2( +[4 h 0 1 +) +( ] + 一 ) 12 一 一 h’+ 一 2 _l2一 一 一2豢 一r (一 一 ’一 ’ ( 一h2+ 2 h+ + f( 一 0 r ) 2 2 + [ ” 2O ) 2 +l o ] i ̄ l2 。14 0 +( 雾] , ++2 + + 02) 5 ) 0 u
得:
( + + + + + ) 一 + 2 u+( + ) 0 ]一 (o + 一 一 ) u + h 2 02 ( +
一c o r l( +
一
+++一 軎 + 0 o 孑 1 3 )
+ [ ( + +
其中 0 , O , ,2 ,0, 0 为待定参数.适 当选取这些参数 ,可以使差分格式 ( )逼近 五 0 o0 0 , , 0 , 6 2
[ 收稿 日期 ]2 0 —1 0 1 0—1 0 [ 作者简介 ]单双荣 (9 6一) 15 ,男 ,讲 师 ,从 事计 算数学研究 .
维普资讯
四阶抛物型方程的一个高精度差分格式
() 2
边 值条件 : ( ,) M L,) 0 0≤ t T。 () “ 0 t一 ( t一 , ≤ 3 利用加 耗散项 的思想来构 造偏微 分方程 的差分
本 文对 高阶抛物 型方 程 ( ) 出如 下 的三 层 多 1提 参数差 分格式 :
收 稿 臼期 :0 8 9~0 2 0 ~0 5
科
学 m
版
A i h Ac u a y Di f r n e S h me f r S l i g H g c r c fe e c c e o ov n
t t d r Pa a o i r i lDif r nta he 4 h Or e r b lc Pa ta f e e i IEqu to a in
W A NG a - n LI Pi g, ANG Xi o Fe g, U n W Bo
( . p r me to t e t s Xi x a g Un v r iy Xi ̄a g 4 3 0 Ch n ; 1 De a t n fMa h ma i , n in ie st , n n 5 0 3, i a c
0问题 的提 出
格 式是 一个 重要 而有 效 的方法 , 文献 [ ,]对 四 阶 12
抛 物型方 程 ( ) 造 了若 干显式 、 1构 隐式和半显 式差分
在渗流 、 扩散 、 热传 导等领域 中经 常会 遇 到求 解 四阶抛物型方程 的 问题 , 一维情 形 , 在 其模 型为 如下 初 边值 问题 :
合 的。
关键 词 : 物 型 方 ; 稳 粮一类 号 : 24程8待 定 系数 法 ; 志 码 : . 抛 O 1.2 文 献 标 定性 中 图分 A
自
文 章 编 号 :6 4 3 6 2 0 ) 3 0 0 - 0 1 7 —3 2 ( 0 8 0 - 0 6 3
高阶抛物型偏微分方程的一个高精度差分格式
关键词 : 抛物型方程 ;截断误差 ; 稳定性 中图分类号 :O241. 82 文献标志码 :A 文 章编 号 :167423326 (2008) 042 00012 02
A High Accuracy D iffer encin g Scheme f or Sol vin g t he Para bol ic Par t ia l D iff er en t ial Equat ion
新乡学 院学 报 (自然科学版) Jour nal of Xinxia ng Univer sity( Na tural Science Edi 2008
高阶抛物型偏微分方程的一个高精度差分格式
蒋菊霞1 , 王 波2
3
(1. 新乡学院 数学系 ,河南 新乡 453003 ;2. 三门峡职业技术学院 机电工程系 ,河南 三门 峡 472000 )
J IANG J u2xia 1 , WANG Bo2
( 1. Depa rtment of Mathema tics , Xinxiang Univer sit y , Xinxia ng ,453003 , China ; 2. Mec hanical and Elect rical e nginee ring Depar tme nt of Sa nmenxia Polytechinc , Sanmenxia 472000 , China) Abstract : This paper pr esents a high accur ac y three2level diff erence scheme fo r solving parabolic equa tio ns. The sta2 bility condition a nd local tr uncation e rror for the sc heme are η < 0 , r ≤[ - 10 ( m - 6 η ) 2 + 360 (η- 1) 2 + 17 m 360 ]/ [ 180( 1 - 2η ) 4 m ] , re spectively. Key wor ds : para bolic diff erential equa tion ;truncation e rror ;stability
解四阶抛物型方程含参数绝对稳定差分格式
;8 , ) 詈一 或 1 及 ≤< 时所出 且< 圭一 亏 > 以 。 ,给 二 主 ( ) , 0 时 1 1( 1 1
.
的格 式绝 对稳 定 , 截断误差 阶数为 O( : . 当参数 适 当选取 时 ,其截断误 差阶数 最高 r +h ) 而 2
可达 of+ . 到 i
维普资讯
第 1 5卷
第 3期
漳 州 师 范学院 学报 ( 自然 科 学版 )
J un l f h n zo ec es l g ( t c. o ra o Z a g h uT a h r Col e Na. i e S )
voII N O. . 5 3
将 格 式 ‘)中各 节点 上 的 2 经整理后得 : 在 网点
,,
() 3
) 或 ( , ) 处展 开的 y。 f l 级数代 入
’ 岍 詈 十 ) + ) 2
”
一
2- +I  ̄ -
4 )
一
㈣
1
+
M :
+
2
2 舞 蛳 + 等 M 1 哪2 ㈤
满 足 一定 的条件 时 , 差 分格 式绝对 稳 定 .蜀部截 断误 差 阶 数 最高 可达 子 说 明对 稳定性 所 作 的分 析 是正 确 的 . 关键 词 四阶抛 物 型 方 程 、隐式 差 分 格式 ,绝对 稳 定 .
中 图 分 类 号 : 4 .2 0 2 18 文 献 标 识 玛 :A 文 章 编 号 : 0 8 7 2 【0 2 0 .“ 30 1 0 .8 6 2 0 ) 3O 4 .6
l— l0 8 = 82 一 i 0 0 0+ 一 一 2
lo , 0 =0】 = 0 1 一= {2 0: 12 / 一 0 = 一= / 一 2 13 — / 4 0 =1 2 ,0 = a
4段抛物线方程
4段抛物线方程
我们要找到一个4段抛物线方程。
首先,我们需要理解什么是抛物线方程。
一个抛物线方程通常可以表示为 y = ax^2 + bx + c,其中a, b和c是常数。
但是,题目要求的是4段抛物线,这意味着我们需要4个不同的抛物线段来组成整个图形。
我们可以考虑以下几种情况:
1. 四个不同的抛物线,每个都有不同的a, b和c值。
2. 两个抛物线,每个有两个部分(例如,y = ax^2 + bx + c 和 y = -ax^2 - bx - c)。
3. 一个抛物线,它被切割成四部分。
为了简化问题,我们假设每一段都是一个完整的抛物线,并且每一段的顶点都是 (x, y)。
这样,我们可以为每一段抛物线分别写出方程。
现在我们要来解这个问题,找出这四段抛物线的方程。
四段抛物线的方程分别为:
y = Eq(a1*x**2 + b1*x + c1, 0)
y = Eq(a2*x**2 + b2*x + c2, 0)
y = Eq(a3*x**2 + b3*x + c3, 0)
y = Eq(a4*x**2 + b4*x + c4, 0)
注意:您需要提供具体的顶点坐标 (x, y) 来解出 a, b 和 c 的具体值。
抛物型方程的O(h 4)精度新交替分段显隐格式
( .Sho o  ̄ l ts n yt i cs hnogU i rt, i n200 ,Sadn ,C i ; 1 col f ln i dSs m S e e,Sadn n e i J a 10 hnog hn Ma e ac a e c n v sy n 5 a
w ih Wa rv d t e hg l e u a y n me ia e p r n s h c s p o e b i y a c _ ̄ b u r l x ei o h re c me t.
Ke yw
:8 1 at , i s g e tsh me o v cin df s n p be 1 ̄ l e  ̄ , e ̄ n c e ;c n e t iu i rl m;p rl l loi m 1 r  ̄ n o o o aal g r h ea t
c m n e a ll m u t n w i a e eucnio l al. h f u c t i O h ) s e e a e sdi prl o pt o , h hws rvdt b nod i aysbe T esaM acr y a ( , h c bu n ae c a i c p o o tnl t pi c a r e s
维普资讯
第4卷 2
V0 . 2 14
第 l 期 2
N .2 o 1
山
东
大
学
学
报
( 理
学
版)
20 年 l O7 2月
De c.2 r 007
Junl f hn ogU i r t( a r cec) ora o adn nv sy N t a Si e S ei ul n
解四阶抛物型方程的两层显式差分格式
解四阶抛物型方程的两层显式差分格式近年来,解决抛物型方程(Parabolic equation)的方程及其解的研究内容越来越受到重视。
它是由于几何的变形、材料的变化、流体的流动等一系列连续空间和时间的概念而定义的方程。
另外,在自然界中,抛物型方程也被广泛应用于描述极限状态,例如平衡、流动、摩擦、热传导等问题。
针对抛物型方程,研究者已经利用显式差分格式建立了一系列解决方案,以更高的准确率来解决这些方程。
显式差分法(Explicit Difference Method)是一种有效的数值求解方法,它是根据抛物型方程的求解公式,用一些简单的数学公式来得出最终要求的解。
在过去的几年中,研究者开发了许多不同的抛物型方程的显式差分格式,其中最被广泛应用的是两层显式差分格式(2-Layer Explicit Difference Scheme)。
两层显式差分格式使用外层循环来设定抛物型方程的未知量,内层循环求解抛物型方程的未知量。
两层显式差分格式有助于减少抛物型方程的计算量,提高解的准确度。
因此,这种显式差分格式是解决抛物型方程的首选方案。
首先,两层显式差分格式可以将抛物型方程转变为一组常微分方程,以计算出正确的未知量。
然后,使用该显式差分格式为抛物型方程设定了一组相应的数值解,例如使用抛物型方程为一维、二维及更多维度时的未知量。
最后,两层显式差分格式可以根据特定的方程设定相应的解算公式,获得更准确的数值解,而不会因为细节而出现误差。
两层显式差分格式可以归纳为下列步骤:第一步,根据抛物型方程确定未知量,以及相应的微分方程,即根据抛物型方程的解算公式计算未知量;第二步,使用两层显式差分格式求解出抛物型方程的未知量;第三步,根据解算出的未知量,计算抛物型方程的解。
总而言之,两层显式差分格式是一种有效的数值计算技术,它可以有效地求解抛物型方程的未知量,从而获得更准确的数值解。
通过部分解和完全解来求解抛物型方程,两层显式差分格式可以有效地提高解的准确率,为研究人员提供更强大的工具。
解抛物型偏微分方程的一个高精度格式
解 抛 物 型 偏 微 分 方 程 的一 个 高 精 度 格 式
王 晓 峰 王 波 ,
( . 乡 学 院 数 学 系 , 南 新 乡 430 ; . 门 峡 职 业 技 术 学 院 机 电工 程 系 , 南 三 门 峡 420 ) 1新 河 50 32 三 河 700
摘
要: 用待定 系数法 给出 了解 一维抛物 型偏微分 方程初边 值 问题 的两层显 格 式 , 此格 式 的截 断误 差为 0( r
格 式 当 属 H dii o 格式 【 , 截 断误 差 为 0( + aj m s d 3其 ] r h ) 此 后 周顺 兴 又 构 造 了一个 精度 更 高 的格 式 , ,
即
Ox Ot 粤 , —x — ’ t O a O
: :
、‘ (2 2) ’ .
( .) 2_ L 3 2・ j
( .) 2 1
其 中 —— 坐 标 变 量 , — 的 已知 函数 ,—— 时 — t 间变 量 , — — 、 的 已知 函数 。 u t
在解 上 述 问题 的 差 分 格 式 中 , 典 差 分 格 式 [ 古 1 】
其 中 , 示 u在 节 点 ( , )= ( k )上 的值 , u表 , t j r h,
0引言
两 层 显格 式 , 精 度与 周 顺 兴 的三 层 隐格 式 相 同 , 其 也 为 O( +h ) 而 稳 定性 条 件 为 29< r≤2 3 其 中 r , / 1,
r= rh / 。
在 渗 流 、 散 、 传导 筝 领 域 中经 常会 遇到 求解 扩 热
抛物 型 方 程 的 问题 , 一 维情 形 , 模 型 为 如下 初边 在 其
— —2 m 、n O m+ ̄,m、 刀 非 负 整数 。 n为非 贝 拦 , 。 X 烈
解二阶抛物型方程含参数高精度两层差分格式
解二阶抛物型方程含参数高精度两层差分格式
解二阶抛物型偏微分方程是许多领域中非常重要的问题,例如热传导、扩散、化学反应等。
然而,在实际计算中,由于参数的存在,往往需要使用高精度的差分格式来保证数值计算的准确性。
以下是一种含参数的高精度两层差分格式,可以用于解决二阶抛物型偏微分方程:
1. 先验估计
为了使用高精度两层差分格式,我们首先需要进行先验估计,以确定合适的时间步长和空间步长。
具体来说,我们可以使用稳定性分析来确定时间步长和空间步长的上限值。
2. 差分格式
在确定了时间步长和空间步长之后,我们可以开始使用高精度两层差分格式来求解二阶抛物型偏微分方程。
该差分格式通常包括以下几个步骤:
(1)先用向前差分公式求解第一层,得到一个中间解。
(2)再采用Crank-Nicolson格式对第二层进行差分,同时使用前一步得到的中间解进行修正。
(3)最后,将得到的数值解反推回到未知函数的值域中,得到方程的
数值解。
需要注意的是,在使用这种高精度差分格式进行计算时,我们需要使
用高精度的算法来保证计算的准确性。
3. 参数调节
由于实际问题中经常存在参数不确定性的情况,因此,在进行数值计
算时,我们需要对参数进行调节和优化。
具体来说,我们可以通过多
次求解不同的二阶抛物型偏微分方程,来不断调节参数并逐步优化计
算结果。
以上是一个含参数的高精度两层差分格式,可以用于解决二阶抛物型
偏微分方程的计算问题。
该方法能够保证数值计算的高精度和准确性,同时也能够应对实际问题中的参数不确定性。
解抛物型方程的一族高精度差分格式
A [ l + 2; 3“+ + 1 / +1 “+ + ;1 2 = [5? +(一 5 7) 7“一] , “ “ 一 + ( l “一 )2 7(n} u2) ] 口 r“ ) 4 - / / 1 —1 “ +1 ; 6 6
() 2
其 中 A, 关 于 t的 一阶 向前 差商 , 是 关 于 的 二 阶中心 差商 , ~ 是 待定 参数 . 是 6 l 6 将 () 2 中各节 点 上的 “以其 在节 点处 展开 的 T y r 数代入 , 利用方 程 ( ) al 级 o 并 1可得
均达 0( £ + △ ) 在 不增 加节 点所 在 区间宽 度 的情况 下 , . 是否 存在精 度更 高 的- , 线型 隐格 式 ?文 [ ] -x 角  ̄ 3 曾给 出 了一个 截 断 误差 为 0( + △£ ) 隐 格 式, 由于 作 者 算 错 了 一个 误 差 项【 格 式 的 精 度 仍 为 t hsp pr T el a t n ainerri O ( t ee t l u t n r r ne i ti ae . h cl r ct r A + ae o s dn o u o o s
a s l e y sa l nd c n a iy s l e o bl we p n t o b ut l t b e a a bee sl v d by d u e s e i g me h d. o o Ke r : ne d me in lp r b i q a in; i h a c a y;mpl i ifr nc c me y wo ds o — i nso a a a ol e u to h g c ur c i c i td fe e es he c
维普资讯
第 1 第 5期 8卷
数学中的抛物型方程
数学中的抛物型方程抛物型方程(parabolic equation)是数学中一类重要的偏微分方程,它在物理学、工程学和社会科学等领域中具有广泛的应用。
本文将从抛物型方程的定义、特征和解法等方面进行论述,以帮助读者更好地理解和应用抛物型方程。
一、抛物型方程的定义在数学中,抛物型方程是一类二维或三维偏微分方程,其形式可以表示为:∂u/∂t = a∇²u + bu + c其中,∂u/∂t 表示函数 u 对时间 t 的偏导数,∇²u 表示函数 u 对空间坐标的拉普拉斯算子,a、b、c 是常数。
抛物型方程通常描述了某一物理现象随时间变化的规律,比如热传导、扩散等。
通过解抛物型方程,我们可以预测和分析这些物理现象。
二、抛物型方程的特征1. 热传导方程抛物型方程在热传导方程中的应用是最常见的。
热传导方程描述了物体内部温度随时间和空间的变化情况。
在一维情况下,热传导方程具有以下形式:∂u/∂t = α∂²u/∂x²其中,u(x, t) 表示在时刻 t 位置为 x 的温度,α 是热扩散系数。
2. 扩散方程抛物型方程在扩散方程中的应用也是非常重要的。
扩散方程描述了物质在浓度梯度驱动下的扩散过程。
在一维情况下,扩散方程具有以下形式:∂u/∂t = D∂²u/∂x²其中,u(x, t) 表示在时刻 t 位置为 x 的物质浓度,D 是扩散系数。
三、抛物型方程的解法对于抛物型方程,我们通常采用偏微分方程的求解方法,如分离变量法、格林函数法等。
1. 分离变量法分离变量法是一种常用的求解抛物型方程的方法。
它的基本思想是将多元函数分解为几个一元函数的乘积,并利用分离后的一元函数满足各自的方程来求解。
以热传导方程为例,我们可以将其分离变量为时间部分和空间部分:u(x, t) = X(x)T(t)代入原方程,得到两个方程:X''(x)T(t)/X(x) = T'(t)/T(t) = -λ²其中,λ² 是常数。
四阶抛物方程的crank-nicolson全离散格式
2019年12月第19卷第4期廊坊师范学院学报(自然科学版)Journal of Langfang Normal University(Natural Science Edition)Dec.2019Vol.19No.4四阶抛物方程的Cran k-Nicolson全离散格式杨晓侠(平顶山学院,河南平顶山467000)【摘要】对一类四阶抛物方程利用非协调EQJ01元和零阶Raviart-Thomas元提出了一个低阶Crank-Nicolson全离散逼近格式。
首先,证明该格式逼近解的稳定性,其次,基于上述两个单元的高精度分析,并借助插值后处理技术,导出了原始变量”的刃-模意义下,中间变量v=-△"的时空能量模意义下以及流量p=-Vu的-模意义下O(h2+T2)阶的超逼近性质和超收敛结果。
【关键词】四阶抛物方程;非协调混合元方法;Crank-Nicolson全离散格式;分裂技巧;超收敛Crank-Nicolson Full-discrete Scheme for A Class ofFourth-order Parabolic EquationsYANG Xiao-xia(Pingdingshan University,Pingdingshan467000,China )【Abstract]In this paper,with the help of nonconforming EQ^f element and zero order Raviart-Thomas element,a low order Crank-Nicolson full-discrete scheme is proposed for a class of fourth-order parabolic equations.Firstly,the stability of the approximation solution is proved.Secondly,based on the high accuracy analysis of the about two elements,with interpolated postprocessing approach,the superclose properties and superconvergence results with order O世+/)for the primitive solution u in H1-norm,the intermediate variable v=-A«in space-time energy norm,and flux p=-Vw in Z2-norm are derived,respectively. [Key words]fourth—order parabolic equations;nonconforming mixed finite element method;Crank-Nicol—son full-discrete schemes;splitting technique;superconvergence〔中图分类号〕0242.21〔文献标识码〕A〔文章编号J1674-3229(2019)04-0005-060引言考虑如下四阶抛物方程的初边值问题u t+A:A2m=/(m),(X,t)e Qx(0,T]■u(X,t)=Au(X,t)=0,(XOe30x(0,71(1)u(X,O)=u o(X),XeQ.其中a是一个矩形区域,oa为a的边界,T g(0,+oo)为一定值,X=(x,y),m0(x),f(X,t)为已知光滑函数,丘为正常数,/'(•)有界。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
中图分类号 :2 18 O 4 .2
收 稿 日期 :0 8 o —o 20 — 5 4
文献标 志码 : A
文章编号 :0 9 5 2 (0 8 o —0 l—0 10 - 18 20 )5 0 4 2
基金项 目: 陕西省教育厅专项研究基金项 目(5K 3 ) 西安建筑科技大学基 础研究基金 (3 R 2 0J 2 9 ; 0B 0 )
第2 3卷 第 5期
解 四阶抛物型方程高精度两层显格式
孙 雪莉 , 曲小钢
( 安建筑科技大学 理学院 , 西 西安 70 5 ) 105
摘
=
要: 给出了一个求解 四阶抛物型方程高精度两层显式差分格式 , 明了其 截断误差 为 O( + 。 , 证 h ) 稳定性 条件为 r
寺 . ≤
差 分 格 式 , 断 误 差 分 别 为 0( +h ) D( +h ), 需 求 解 联 立 方 程 组 , 算 量 较 大 . 对 这 一 问 题 , 截 r 及 . r 但 计 针 曾 文 平 提 出 了一 个 三 层 含 参 数 的 显 式 差 分 格 式 , 截 断 误 差 为 D( +h ), 度 得 到 了 一 定 的 提 高 , 却 其 r 精 但 使 用 了较 多 的 网格 结 点 . 文 在 上 述 研 究 的 基 础 上 , 造 了 一 种 新 的 两 层 显 式 差 分 格 式 , 截 断 误 差 为 0 本 构 其 ( 2+h ), 稳 定 性 条 件 为 r=三 r t 其
作者简介 : 孙雪莉 (9 3 )女 , 18一 , 陕西扶风人 , 西安建筑科技大学硕士研究生 ; 曲小钢 (9 9 ) 男 , 15 一 , 山西定襄人 , 西安建
筑科技大学理学院教授 , 硕士生导师.
考 虑 F列 四 阶 抛 物 型 方 程 的 初 边 僵 问 题 : I 等 :0, 等
=
- + 【 号
,
㈣
同理 , 用 T yo 公 式 可 以 得 到 4阶 和 8阶 微 商 的 下 列 差 分 近 似 : 利 al r
【 = 【 4(, 喀) ( 吐) 吐一 +- r 一1吐一 + 吐一 + ( ) 1丽 6 2 + 2 一 2 n 3
+ ( 咄 一 +吐 : 16 ( 2 )一面 9 6
h4c 6 - 2  ̄ = 4
,
进 一 步 扩 充 了这 一 问 题 的 已有 研 究 成 果 .
1 差 分 格 式 的 构 适
设 h和 . 别 表 示 空 间 步 长 和 时 间 步 长 , 域 由求 解 区 域 上 的 网格 点 集 { , } 成 , 中 = h, = r分 网 t 组 其 j t
这 是 一 个 两 层 显 格 式 , 截 断 误 差 为 0( +r 注 意 到 0( h ) mz 0( ), h )} 其 . r h ), r x{ 0( 。 .
故 格 式 ( 的 截 断 误 差 实 为 0( +r 。 . 8) h)
2 稳 足 性 分 析
用 分 离 变 量 法 研 究 差 分 格 式 的稳 定 性 , n e 1 1 令 = ( ≤仃 , = = ) 入 差 分 格 式 ( ) 可 得 传 0 i 『 代 8 ,
播 子 =+ls 0 9i詈 2 s 导 soJ 因 1【6 — s + 6 。 + i 02 3i 2 6 5i n n n n ,1 r
[丽 一 +2 + _1 竽2吐 一 + )【1 22u 一 + ) 丁2 咛)【 一 r 2 U + 4 一 r(,2 , 4r 1 1 r 一 n 瓦 r 】 r 3]] 0 5 ( n 5 +
一62u 一 n + 一1r(+ 2 + 。 8] ̄ 2 [ r . ) (: 22u 一 吐 ) 6  ̄ ] 。 + ( 8 )
。2+ ) o 。 一 ; 】 () u + h
( 6 )
() 7
【 【( 2+ + ( 2 + 一 (, 吐) ] 一吐 一 )1u 一 吐)6u 一 + , 3T + 9T 2 +
+2 4( o n . 2—2 + 2 )一3 8 l + j ]+D( 7 ( 一2 ) h)
n (, r j n为 整 数 ) 又 设 r= 为 网 格 比 , , 令 =u(h,r) 由 T yo 公 式 可 得 ㈨ : j n , a lr
=
【】 r ̄J D ) 詈 "+( r1
( 4 )
注 到 程-又 詈 一 】 【 人4,: 意 方 c ,有 = 【 , 】代 ( 得 】【 = , )
+
-
0≤ ≤ 1, l o≤ ≤
( )
v
--
{ = “ 0≤ ( ) ≤1 【 , 宣 “, t 2 =1 = ()_ () 0= 三 t
= o≤ o ≤ , t
( 2 ) ( 3 )
16 9 0年 , a], 构 造 了一个 显式差 分格式 … , 截 断误差 仅 为 O r+h ) 他 所提 出的 另外 两个 隐式 C y ̄ B ie 但 ( ,
维普资讯
20 0 8年第 5期
孙 雪莉 , : 等 解四阶抛物型 方程 高精 度两层显格式
・1 5・
将 ( ), 7 代 入 ( ), 去 高 阶 误 差 项 , 可 得 方 程 ( ) 下 列 差 分 格 式 : 6 () 5 略 即 1的
“ =
维普资讯
20 0 8年 9 月
渭南师 范学院学报
J u a f en n T a h r Unv r i o r lo i a e c e s n W ie t s y
S p . 2 08 et 0
Vo _ 3 NO. l2 5