抛物型方程的时空有限元方法

来源 :数值计算与计算机应用 | 被引量 : 0次 | 上传用户:jmzsren1
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
一、 引言 抛物型方程的数值求解,流行的是有限差分法和半离散Galerkin方法,本文在最近发展起来的活动边界问题的时空有限元方法的(见文[1]和[2])基础上,将这种方法移用来求解固定边界的抛物型方程,数值实验和结果分析证明,本方法精度高,稳定性好,相当经济有效。
其他文献
近二十年来发展起来的分裂法在求解高维数学物理问题方面得到了广泛的应用,在实际应用中通常按坐标方向进行分裂或按物理问题进行分裂,我们对于球头驻点线高超音速粘性绕流的计算,按问题的物理性质构造了分裂隐式差分格式SISOMS,计算结果与
期刊
设f:R~n→R为R~n上的二次可微函数,且已知其二阶导数Hesse矩阵▽~2f具有某种稀疏、对称型结构,[1]已叙述了求解这类问题的各种稀疏拟牛顿方法,它们仅利用目标函数及其梯度的变化来估计二阶导数的近似海色矩阵,并在每次校正的迭代中保持与海色
期刊
一、前言 变分问题的近似解法,由于实际上的需要,曾为许多学者所研究,例如,欧拉、里兹、米赫林、康托罗维奇等人都对变分问题的近似解法作过研究,并给出了一些方法。但是,在过去已有的方法中,有的精度不高(如欧拉的折线方法);有的在具体执行中甚为困难
期刊
一.Lewis算法存在的问题 1982年J.Lewis在[1]中给出了不对称三对角阵求逆的算法,设三对角阵为
期刊
一、问题的提出 在托卡马克装置上,要估计快波加热的效果,需要计算天线的耦合阻抗,在计算椭圆截面半环天线等效负载阻抗时,要用到马丢函数,即要计算马丢方程
期刊
一、前言 对于超越方程 f(x)=0 (1)到目前为止已有许多求解的方法,对于各种方法,当所求的根为单根时,收敛性质一般都较好;当所求的根为p(p≥2)重根时,各种方法的收敛速度都要减慢,且p越大收敛越
期刊
§1.引言 我们讨论整数Bottleneck问题其中c_i>0,a_i>0(i=1,2,…,n),b>0,且设(否则目标函数值只能取另值),记,整数},并假设Τ≠Φ,因而(1.1)一定有解。 引理1.1 (1.1)的解x~*满足下列估计:
期刊
一、前言 气藏模拟,意指从气藏模型的特性推断出实际气藏特点的过程,一旦涉及具体气藏,往往由于底水锥进和气层静、动态参数分布及其变化不清楚,难以编制出科学合理的气藏开发调整方案,为了改善气藏开发现状,提出保藉气藏稳产和提高最终采收率的措施,需研究气藏静、动态特征,掌握气藏(#)水浸规律,为合理开发气田提供依据,时至今日,数
期刊
一、问题的提出 关于无约束最优化问题 min f(x), x∈R~n在[1]中曾提到曲线线性搜索是一个求解这类问题的令人感兴趣的研究课题,采用的一般曲线为 x(α)=x+Φ_1(α)8+Φ_2(α)p(a≥0参数),其中Φ_1(0)=Φ_2(0)=0,且Φ_1,Φ_2应满足条件 Φ_1~′(0)=0,Φ_2~′(0)>0,Φ_1~″(0)>0.最简单的是取
期刊
程序实现的是Stiff稳定的1—6阶的方法族M_p(ε)(下标p表示方法的阶,ε是方法族依赖的参变量),详见[1]。现将各阶方法M_p(ε)相应的生成多项式及表征其Stiff稳定性特征的参数α,D(含义见图1)列于下面: 二阶方法M_2(ε):
期刊