Poisson-Nernst-Planck方程的两网格法与自适应有限元方法

论文摘要

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模型问题的计算效率.

论文目录

• 摘要
• Abstract
• 第一章 绪论
•   1.1 研究背景与现状
•   1.2 本文主要工作
• 第二章 预备知识
•   2.1 Sobolev空间
•   2.2 有限元估计中的基本不等式
•     2.2.1 有限剖分与有限元空间
•     2.2.2 基本不等式
•   2.3 Cl′ement插值与梯度重构算子
• 第三章 含时Poisson-Nernst-Planck方程的有限元误差估计
•   3.1 PNP模型与变分形式
•   3.2 有限元离散与投影算子
• 2 误差估计'>  3.3 最优 l²误差估计
•     3.3.1 半离散误差估计
•     3.3.2 全离散误差估计
•   3.4 数值实验
•   3.5 本章小结
• 第四章 含时Poisson-Nernst-Planck方程的两网格有限元方法
•   4.1 半解耦算法
•   4.2 全解耦算法
•   4.3 数值实验
•   4.4 离子通道问题中的应用
•   4.5 本章小结
• 第五章 非线性稳态Poisson-Nernst-Planck方程的自适应有限元方法
•   5.1 模型方程与变分问题
•   5.2 MNSPNPE的有限元逼近
•   5.3 MNSPNPE的后验误差估计
•     5.3.1 上界估计
•     5.3.2 下界估计
•     5.3.3 自适应算法
•   5.4 数值实验
•   5.5 本章小结
• 总结与展望
• 参考文献
• 致谢
• 个人简历、在读博士期间发表的学术论文和研究成果

文章来源

类型: 博士论文

作者: 沈瑞刚

导师: 舒适,阳莺

关键词: 方程,有限元误差估计,两网格方法,梯度重构,后验误差估计,自适应有限元方法

来源: 湘潭大学

年度: 2019

分类: 基础科学

专业: 数学,数学

单位: 湘潭大学

分类号: O241.82

DOI: 10.27426/d.cnki.gxtdu.2019.001727

总页数: 136

文件大小: 6848k

下载量: 10