了解纽曼 BC

计算科学 有限元 边界条件
2021-11-28 06:01:30

我理解 Neumann BC 的物理含义以及如何暗示它们。但是,我无法完全理解它在数学上表示为的方式u/n在哪里n是边界的外法向量,并且u是解决方案。

我知道法线向量对于计算表面积分是必不可少的。但是我们在这里做什么?牵引力是解向量wrt对法线向量的偏导的结果?这是什么意思?

任何人都可以尝试简单地简化这一点吗?

非常感谢!

2个回答

基本概念

法线向量对应一个物理参数,例如笛卡尔坐标x. 例如,沿x方向可以由达西定律定义如下:

ux(x,t)=kμPx

根据 Neumann 的 BC,通量在边界处应为常数或为零,即u(0,t)=u(L,t)=0

换一种说法,kμP(0,t)x=kμP(L,t)x=0

例子

考虑一个有 2 个网格的系统,如下所示:

在此处输入图像描述

如果我们将上述概念应用于该系统,我们可以得到以下方程:

PiPi1Δx=P1P0Δx=0

P1=P0

当我们写u/n这只是一个简写u˙n. 因此,我们采用梯度和单位法线的点积。以矩形区域为例:

在此处输入图像描述

左墙上的法线单位是 (1,0) (或 (-1,0) 取决于您的定义)。因此,诺伊曼条件将是u/x=0在左边墙上。类似地,对于底部边缘,单位法线为 (0,1),因此诺伊曼条件为u/y=0.

就表面积分而言,我认为维基百科的解释非常好。它部分说维基百科

或者,如果我们对向量场的法线分量进行积分,则结果是一个标量。想象一下,我们有一个流过 S 的流体,因此 v(x) 决定了流体在 x 处的速度。通量定义为每单位时间流过 S 的流体量。该图表明,如果矢量场在每个点都与 S 相切,则通量为零,因为流体仅与 S 平行流动,既不流入也不流出。