【摘 要】
:
Poisson-Nernst-Planck(PNP)方程是由Poisson方程和Nernst-Planck方程组合而成的强耦合非线性偏微分方程组.此类方程广泛用于描述生物化学的静电扩散反应过程、半导体的离子输运以及生物细胞膜间的离子转换等应用领域.有限元方法是求解PNP方程的一种流行离散化方法,因此,研究PNP方程的有限元误差估计及其快速算法具有重要的理论意义与实际应用价值.本文主要开展了以下三个
论文部分内容阅读
Poisson-Nernst-Planck(PNP)方程是由Poisson方程和Nernst-Planck方程组合而成的强耦合非线性偏微分方程组.此类方程广泛用于描述生物化学的静电扩散反应过程、半导体的离子输运以及生物细胞膜间的离子转换等应用领域.有限元方法是求解PNP方程的一种流行离散化方法,因此,研究PNP方程的有限元误差估计及其快速算法具有重要的理论意义与实际应用价值.本文主要开展了以下三个方面的研究工作.针对描述离子质量守恒的经典含时PNP方程,通过引入一种合适的有限元投影算子,首先证明了其半离散线性有限元离散格式的最优2误差估计;接着对一种应用广泛且更能保持PNP方程物理特性的向后Euler全离散线性有限元格式建立了最优2误差估计理论.数值实验验证了理论结果的正确性.针对含时PNP方程的全离散有限元方法,构造了半解耦和全解耦格式的两网格算法.该算法基于解耦思想,通过使用粗网格空间的有限元解作为细网格解的一个可靠逼近解,在每个时间层上将原始的全离散耦合系统解耦成更小的独立系统,与标准有限元方法相比提高了求解PNP方程的计算效率.基于本文推导的有限元解的最优2误差估计,获得了静点势的两网格有限元解在1范数下的最优误差估计;对浓度的两网格有限元解分别给出了2和1误差估计.数值实验表明,当网格尺寸和?满足=(?1/2),两网格法具有与标准有限元法相同的误差收敛阶.包括实际应用的离子通道问题在内的一系列数值实验验证了两网格算法的有效性,更少的CPU时间表明两层网格法大大地提高了有限元方法求解PNP方程的效率.针对一类改进的非线性稳态PNP模型的线性有限元格式,讨论了其梯度重构型后验误差估计.首次对该模型问题进行了严格的后验误差分析,并获得了静电势和浓度的后验误差估计子的上界和下界估计,设计了基于非线性PNP模型后验误差估计子的自适应有限元算法.特别地,本文推导的后验误差估计子依然适用于经典的稳态PNP方程.若干数值实验验证了后验误差估计子的可靠性和有效性,自适应计算提高了奇性PNP模型问题的计算效率.
其他文献
本文我们考虑几类常见的流体方程,研究它们的强解及相关极限问题,也就是,局部解的粘性消失极限和整体解的衰减这两类问题。更确切地说,粘性消失极限问题是指,当粘性系数或扩散系数趋于零时,粘性流体方程的解收敛到无粘性或理想流体方程的解。在有界区域,边界条件将是一个关键,我们主要考虑的是Slip边界下粘性消失极限问题。而整体解的衰减问题是一个大时间行为,是指当时间趋于无穷大时,能量趋于零,本文也包括衰减率和
E.T.A.霍夫曼是德国18世纪末、19世纪初的著名作家、音乐家、指挥、舞台设计。他被誉为歌德和海涅之间最具世界影响的德国作家。他的很多作品都表现了神秘力量对人的控制,使得人产生分裂和异化,因此他很早就被称为“幽灵霍夫曼”。但他的作品中也确实存在体现和谐、统一的因素,而这些因素都可以归到他最有名的童话小说《金罐》里的“阿特兰提斯”这一概念之中。籍由“阿特兰提斯”表达出来的思想内涵贯穿于霍夫曼的所有
作为计算流体力学研究的一个重要内容,双曲守恒律方程的数值解法在流体力学发展过程中占据着非常重要的地位。在层出不穷的数值计算方法之中,高精度、高分辨率的数值计算方法因为其具有良好的特性,在计算流体力学的发展中占据着重要的地位。本文的主要目的是研究几类具有高分辨率、高精度的数值格式。具体内容如下:首先,我们基于有限体积法思想,通过增加光滑因子中非光滑部分的权重,提出了能有效提高CWENO-Z格式分辨率
随着社会经济的飞速发展和科学技术水平的日益提高,人类面临着越来越复杂的实际决策环境,而决策者认知的模糊性和决策因素的不确定性,导致了决策者往往难以获得确定的决策信息;并且,在实际决策过程中,影响决策者的决策信息越来越多,由此产生的决策数据的维数和量级也越来越大;致使决策者难以进行有效且理性的决策。作为现代决策理论与实践的重要组成部分,直觉模糊多属性决策能够有效的模拟复杂决策环境,而原有基于矩阵理论
弱有限元方法(weak Galerkin finite element methods,简称WG方法)是最近发展起来的求解偏微分方程的有效数值方法.它的主要思想是利用弱微分算子代替传统意义下的微分算子,然后把其应用到通常的变分形式中以数值求解偏微分方程.弱有限元方法的逼近函数为分片间断多项式,逼近函数在单元与单元之间的联系则通过单元边界上的特定多项式实现.自从弱有限元方法在2011年被王军平和叶秀
两相渗流驱动模型多用于石油资源的运移聚集数值模拟,描述在盆地发育中油水运移聚集演化的历史,它对于油田的勘探和合理开发有着极其重要的价值.近年来流体动力学在油藏模拟及地下水污染等重要工程领域应用的研究取得了重大的进展,在模拟油藏生成发展过程,尤其是在不同热量和应变压力作用下的进化过程有重要意义,其数学模型是一组多层对流扩散非线性耦合系统的动边值问题,由于这些方程具有强非线性且相互耦合,标准的有限差分
代数Riccati方程是一类特殊的矩阵方程,在科学计算和工程应用中发挥了重要作用。通过对各种不同的现实问题进行建模,包括最优控制、队列模型、输运理论、与粒子束的传输有关的应用和Markov过程,都可以发现与代数Riccati方程有关。加倍算法是用来求解代数Riccati矩阵方程其及相关方程的一种先进有效的方法,主要利用矩阵和矩阵束的特定结构,以及由不变子空间或压缩子空间来得到矩阵方程的解。特征选择
有限元方法因其具有完善的数学理论及对不规则几何区域较强的适应性等特点,被广泛应用于科学与工程计算领域。虽然对有限元方法已经有大量的研究工作,但仍有一些问题值得进一步探讨。标准有限元先验误差估计只给出了网格尺寸与有限元误差之间的渐近关系,但没有体现网格质量(如单元形状和大小、网格对称性)对有限元解逼近精度的影响。本文借助单元分析,构造了两个可计算量Ge和Gv来刻画网格质量和有限元误差之间的关系,从而
大量研究表明具有非局部特性的分数阶微分算子非常适用于描述具有记忆特性和遗传性质的材料.因此,近年来分数阶微分方程得到了广泛的关注和应用.然而很多分数阶微分方程的解析解是很难得到的,于是在实际应用中数值模拟成为研究分数阶微分方程的一个重要手段.本文致力于二维Riesz空间分数阶扩散方程及分数阶Sine-Gordon方程的有效数值格式及快速算法的研究.第二章中,采用ADI-CN格式将二维Riesz空间
Allen-Cahn方程作为描述相场模型最基本的方程之一,是用于模拟在一定温度下二元合金相位分离的模型。Allen-Cahn方程在实际问题中的应用非常广泛,如晶体生长、相变、图像分析、晶粒生长、材料科学中的界面动力学等。而且在实际问题中复杂的Allen-Cahn方程不易求得精确解,故而,如何数值求解该方程就变得尤为重要。在第三章中,为了数值求解Allen-Cahn方程,我们采用了线性元对其进行空间