非定常线性双曲方程的差分_间断Galerkin方法的精度估计.docVIP

———间断 Galer kin 方法的精度估计赵建丛 康 彤( 河北农业大学) ( 北京广播学院 ,中科院计算数学与科学工程计算研究所)【摘要】 本文对非定常线性双曲方程的差分间断 ———Galer kin 格式进行了误差分析 ,作 了收敛阶估计 ,得到了 L ∞(0 , T ; Hr + 1 (Ω) ) 模的次最优估计 。【关键词】 差分 —间断 Galer kin 方法 精度估计1 引言在5 中 , 曾建立了非定常线性双曲方程的差分 ———间断 Galer kin 格式 ( Finite Difference Disco ntinuo us Galer kin scheme , 以下简称 FD - D G 格式) ,这种方法简化了传统的时 —空剖分 方法的计算难度 ,本文进一步对 FD - D G 格式进行了误差分析 ,得到 L ∞( 0 , T ; Hr + 1 (Ω) ) 模的 次最优估计 。假设Ω R2 是有界区域 ,具有逐片光滑边界5Ω的有界区域 。So bolev 空间 W1 , ∞(Ω) = { v| v ,DV ∈L ∞ (Ω) } , 其范数 ‖v ‖1 , ∞,Ω = sup|α| ≤1‖Dαv ‖ , 半范数为 | v | = sup0 , ∞,Ω 10 , ∞,Ω 1 , ∞,Ω‖Dαv‖0 , ∞,Ω ,其中 , ‖v ‖0 , ∞,Ω = ess sup| v| , Hm (Ω) 表示 P = 2 的 m 阶 So blev 空间 。ΩL (0 , T ; x)假设 X 为 Banach 空间 ,其范数用 ‖1 ‖X 表示 ,映射φ( x ,t ) : 0 , T →L (0 , T ; x)∫T 21/ 2 = ‖ ‖2 ∞0 ‖φ‖x ( t ) dtφ L ( x) ; ‖φ‖L (0 , T ; x) = sup { ‖φ‖x ( t ) } = ‖φ‖L ∞ ) 。( x0 ≤t ≤T2 问题的陈述考虑非定常一阶双曲方程的第一边值问题 (A) :σ Ω5 σ Ω5 t + a ( x , t ) A u + ( x , t ) u = f ( x , t ) ( x , t ) ∈ ×[ 0 , T ] ( 211)u = g ( x , t ) x ∈ г- ( t ) t ∈[ 0 , T ] ( 212)u ( x , 0) = u0 ( x ) x ∈Ω ( 213)其中 β,( x ,t ) = (β1 ( x ,t ) , (β2 ( x ,t ) ,г- ( t ) 称为内流边界 ,定义为 г- ( t ) = { x ∈5Ωβ,( x ,t ) n ( x) 0} ,n ( x) 为Ω的外法向 。类似地 ,定义 г+ ( t ) = { x ∈5Ωβ,( x ,t ) n ( x) ≥0}为外流边界 。显然 ,г-关于问题 (A) 的解 、系数及右端项 ,假定满足条件 Ⅰ:a) 对任意 t ∈0 , T , βi ( x ,t ) ∈W1 , ∞(Ω) ,i = 1 ,2 ;b) 对任意 t ∈0 , T , σ( x ,t ) ∈W1 , ∞(Ω) ∩C (Ω) ;c) f ( x ,t ) ∈L ∞(0 , T ;L 2 (Ω) ) ,g ( x ,t ) ∈L ∞(0 , T ;L 2 (г- ( t ) ) ,初始条件 u0 ( x) ∈Hr + 1 (Ω) ;d) 问题 (A) 的解 u 存在且唯一 ,并满足 u ∈L ∞(0 , T ; Hr + 1 (Ω) ∩C (Ω ×0 , T ) ,5u ∈L ∞(0 ,5tT ; Hr + 1 (Ω) ) 。e) 对任意 t ∈0 , T , (βi ( x ,t ) i = 1 ,2) ,关于 x 满足 Lip schitz 条件 ,即存在常数 L 0 ,使| βi ( x , t ) - βi ( x′, t ) | ≤ L | x - x′|由于 г- ( t ) 随着 t 的不同而变化 ,因而随着 t 的改变 ,需要对Ω采取不同的剖分 。首先在 t = t 0h= 0 时 ,对Ω采取拟一致三角剖分 ,用 Tt0 表示 ,其网格参数为 h0h,假定对任一 t ∈0 , T 及 г-h(t ) ,对Ω均采用拟一致三角剖分 Tt ,并要求每次剖分的网格参数 hth均满足 P0≤ht / h0≤P1 ,其h中 ,P0 , P1 为两个确定的正常数 ,为了记号简单 ,对于每次剖分 Tt ,用 K 表示其剖分单元 ,5 Kh表示 K 的边界 γ,( x) = (γ1 ( x)