使用matlab求解二维浅水方程的数值解（一）

最近在读《ocean modelling for beginners》这本书，对于做海洋数值模拟工作的小白来说，这绝对是一本好书。强烈推荐给理论基础较弱的学习者，这本书循序渐进，由简入繁的讲解了物理海洋方面的物理现象、理论基础、算法实现等等。

第五章讲的是二维浅水方程，书中进行了一些案例的模拟，并给出了相关代码。其中，模拟部分的代码使用Fortran，后处理部分使用Scilab。我试着用matlab，按照书中的有限差分方法，写了该部分的模拟及后处理代码。首先，原始方程组如下：

离散方法如下：

模拟案例设置如下：
与书中练习8一样，模拟区域为500*500m，空间步长10m，时间步长0.1s，重点内容初始条件为网格中心处初始水位1m。
模拟结果如图所示：

这个模拟案例比较简单，模拟的是水位起伏引起的波动。下篇博客将模拟波浪的折射，看看水深地形将会对波浪造成什么样的影响，敬请期待。

使用matlab求解二维浅水方程的数值解（一）—浅水波相关推荐

使用matlab求解二维浅水方程的数值解（二）—波浪的折射
如果大家去过海边,会有这样的感受:如果你面向大海,不管海岸是平直还是蜿蜒曲折的,你感觉到海浪总是迎着你传过来.这是波浪传播中的一种物理现象--波浪的折射.波浪由远海传入近岸的过程中,随着水深变浅,波浪 ...
二维有限元方程matlab,有限元法求解二维Poisson方程的MATLAB实现
有限元法求解二维 Poisson 方程的 MATLAB 实现陈莲a ,郭元辉b ,邹叶童a ( 西华师范大学 a．数学与信息学院; b．教育信息技术中心,四川南充 6437009) 摘要: ...
galerkin有限元法matlab实现,有限元法求解二维Poisson方程的MATLAB实现
有限元法求解二维Poisson方程的MATLAB实现陈莲a,郭元辉b,邹叶童a [摘要]文章讨论了圆形区域上的三角形单元剖分.有限元空间,通过变分形式离散得到有限元方程. 用MATLAB编程求得数值 ...
二维对流方程matlab求解,二维对流扩散方程的有限元计算方法
冯立伟+张成+屈福志 " " " 摘要:针对二维对流扩散方程边值问题,采用三角形剖分,使用二维线性有限元进行计算分析.采用matlab编写了计算程序,使用算例进行了数值实 ...
有限元方法求解二维拉普拉斯方程C++实现
文章目录前言问题区域方程程序设计几何区域网格单元刚度矩阵组装数值结果问题区域网格 u 值图像(结果导出借助Matlab画图) 总结前言本文利用C++语言实现在二维任意区域(内部 ...
水深平均的二维浅水方程推导
浅水方程推导将三维的基本方程沿水深积分平均,即可得到沿水深平均的平面二维流动基本方程. 定义水深为,.为基准面下液面水位和河床高程: 定义沿水深平均流速为: 引用莱布尼兹公式自由表面及底部运动学条 ...
matlab求解二维矩阵并画图,Matlab教程2_ 绘图 _ 二维（2）
(作者:lcc) 二维曲线绘图的基本操作 n plot指令的基本调用格式 (1)plot(x) n x为向量时,以该元素的下标为横坐标.元素值为纵坐标绘出曲线 n x为实数二维数组时,则按列绘制 ...
[计算流体力学][Matlab] 使用 A,B,C 格式与蛙跳格式求解二维对流问题
1. 题目 2. 转述原题目要求可以简化为: 对于二维对流方程: ∂u/∂t+∂u/∂x+∂u/∂y=0 u(x,y,0)={█(1,when-4≤x≤4,-4≤y≤4@0,other)┤ 取计算范 ...
有限元方法基础-以二维拉普拉斯方程为例(附程序)
文章目录前言问题描述问题区域偏微分方程变分问题(弱形式) 有限元离散二维线性有限元 : 参考基函数 2D linear finite element : affine mapping 有限 ...

使用matlab求解二维浅水方程的数值解（一）—浅水波

使用matlab求解二维浅水方程的数值解（一）—浅水波相关推荐

最新文章

热门文章