GMRES方法求解二维不可压Stokes方程(附MATLAB代码)

最新推荐文章于 2024-09-29 08:52:28 发布

waitingwinter

最新推荐文章于 2024-09-29 08:52:28 发布

阅读量7.2k

点赞数 5

分类专栏：计算数学文章标签： MATLAB GMRES 预优GMRES方法

本文链接：https://blog.csdn.net/waitingwinter/article/details/103863358

版权

一、问题描述

考虑如下不可压的Stokes方程
$\left\{\begin{array}{cc} -\lambda \Delta \textbf{u} +\nabla p=\textbf{f}& in \;\Omega,\\ \nabla \textbf{u} = 0& in \; \Omega . \end{array}\right.$

通常设 $\textbf{f}=0,$ 令 $$\textbf{u}=(u,v)^T,\; \Gamma=\partial \Omega.$\\$
定义流函数 $$\psi$,$ 满足

$u=\frac{\partial \psi}{\partial y}, \; v=-\frac{\partial \psi}{\partial x}.$$

显然有 $\nabla \textbf{u}=0.$
定义窝度 $\omega$: $\omega = -\frac{\partial u}{\partial y}+\frac{\partial v}{\partial x},$$ 经计算, stokes 方程可化为

$\left\{\begin{array}{cc} -\Delta u -\frac{\partial \omega}{\partial y}=0& in \;\Omega,\\ -\Delta v +\frac{\partial \omega }{\partial x}=0& in\;\Omega,\\ \lambda \Delta \omega =0& in\;\Omega. \end{array}\right.$

x相应的边界条件为

$\left\{\begin{array}{cc} u=1,\;v=0,\;\omega=-\frac{\partial u}{\partial y}& y=1,\; 0\leq x\leq 1,\\ u=0,\;v=0,\;\omega=-\frac{\partial u}{\partial y}& y=0,\; 0\leq x\leq 1,\\ u=0,\;v=0,\;\omega=\frac{\partial v}{\partial x}& x=1\;or\;0,\; 0\leq y\leq 1.\\ \end{array}\right.$