Matlab
求解偏微分的代码欧米茄方程
关于
Omega
方程是一个与时间无关的偏微分,它控制着水的垂直速度。
通常,控制流体的方程与时间无关,但一种特殊的近似,称为准地转近似,允许将复杂的
Navier-Stokes
方程简化为这种更简单的形式。
这种近似在历史上一直用于公海,但最近,Libe
Washburn
和
Chris
Gotschalk
发现浮游植物的存在远低于阳光允许它们生存的深度。
使用他们在圣巴巴拉海峡收集的数据,人们可以解出水柱垂直速度的欧米茄方程,作为浮游植物惊人深度的可能解释。
运行代码
只需克隆
repo
并运行
Matlab
scipt
omega_eqn_zero_neumann.m。
技术信息
*.mat
文件是运行代码所必需的数据。
omega方程是一个非常数系数的二阶椭圆pde,所以我们用GMRES迭代求解。
域为直角棱柱,每个轴方向的步长恒定。
两个matlab文件对应不同的边界条件:
omega_eqn_zero_dirichlet_fft_z.m
对应于盒子所有面上的零狄利克雷边界条件。
这在物理上是不合理的,包括用于测试目的以及与学术文献进
2022-08-24 22:08:00
6.55MB
系统开源
1