用FreeFem++求解二维Boussinesq系统:平底情形
数值分析
2012-07-10 v2 偏微分方程分析
摘要
我们考虑二维空间中的不同Boussinesq系统族。这些系统近似三维欧拉方程,由三个耦合的非线性色散波动方程组成,描述理想流体在水平底床上小振幅长表面波的传播。我们给出一个FreeFem++代码,用于数值求解这些系统,其中采用P1有限元空间离散,时间上采用二阶Runge-Kutta格式。我们详细给出代码,使用网格自适应技术。对所用算法进行了优化,并对不同Boussinesq族的解进行了比较。所得结果与文献一致。
引用
@article{arxiv.1205.3147,
title = {Solution of 2D Boussinesq systems with FreeFem++: The flat bottom case},
author = {Georges Sadaka},
journal= {arXiv preprint arXiv:1205.3147},
year = {2012}
}