流形上的散度公式证明.docx
- 文档编号:24214174
- 上传时间:2023-05-25
- 格式:DOCX
- 页数:22
- 大小:482.92KB
流形上的散度公式证明.docx
《流形上的散度公式证明.docx》由会员分享,可在线阅读,更多相关《流形上的散度公式证明.docx(22页珍藏版)》请在冰豆网上搜索。
流形上的散度公式证明
流形上的散度公式证明
杨科
中国成都610017
E-mail:
more2010e@
摘要:
散度公式(又称Остроградский-Gauss公式)是现代数学、物理体系的核心公式之一.传统的散度公式证明逻辑体系,建立了基于空间直角坐标系投影法(简称投影法)的曲面积分与三重积分的公式关联,确立了投影法为曲面积分的根本方法.但是投影法存在诸多明显的缺陷(例如计算过程繁琐,不适用于不对称、不规则曲面等),以致于物理、工程领域的许多重要问题(例如电磁学领域的Maxwell方程组实例化和流体力学领域的任意不规则控制面积分)的解决途径,均建立在直角坐标系或其它坐标系的偏微分方程组求解基础上.一个多世纪以来的数学、物理和工程实践已经证明,通过投影法、直角坐标系或其它坐标系的偏微分方程组,难于甚至不能获得关于复杂几何对象[流形]的解析解、数值解;
传统的流形微积分学,用外微分形式推导出Green公式,Остроградский-Gauss公式,Stokes公式,乃至关于n维空间积分的广义Stokes公式[20],即
但是这类用外微分形式推导出的公式只具有抽象的理论意义,并没有揭示积分的具体实现过程,更无具体数值模型可言;
本稿件通过建立与具体几何对象(流形)匹配的个性化坐标系(即有什么样的几何形体,就建立什么样几何形体的坐标系,使用什么样几何形体的微元系数;而不再依赖于已有的少数几个直角坐标系、球面坐标系、柱面坐标系、广义球面坐标系及其相关微元系数等),用积分以及和式极限的方法,证明散度公式在无穷多个任意参数曲面(流形)坐标系[包括单连通可定向闭合曲面坐标系(基于Poincare猜想)和复连通可定向闭合曲面坐标系(环面坐标系)]的存在,使散度公式超越传统的直角坐标系框架,建立基于参数化空间点积法的曲面积分与基于个性化微元系数的三重积分之间的新公式关联,并且在无限丰富、绚丽的公式数值模型运算中实现两种类型积分相互验证,确立两种新型的积分方法的理论逻辑依据和数值模型.
"证明流形上的散度公式"本身不是唯一目的,"建立基于参数化空间点积法的曲面积分与基于个性化微元系数的三重积分之间的新公式关联,确立两种新型的积分方法的理论逻辑依据和数值模型"是根本目的.
本系列稿件相关的数值模型表明,使用基于参数化空间点积法的曲面积分以及基于个性化微元系数的三重积分,能够获得关于复杂几何形体[流形](尤其是不对称、不规则曲面及其包含空间区域)的解析积分值或任意精度浮点积分值;实现任意曲面积分以及任意空间区域三重积分,实现向量场(电场、磁场、流体场、引力场等)和数量场(电位场、温度场、密度场等)在任意自由曲面及其包含空间区域的精确积分计算,确立两种类型积分的逻辑关联关系,实现流形上的散度公式和工程意义上的流形积分.
关键词:
微积分学拓扑学物理学Poincare猜想向量场数量场自由参数曲面坐标系
单连通可定向闭合参数曲面坐标系复连通可定向闭合参数曲面坐标系
基于参数化空间点积法的曲面积分基于个性化微元系数的三重积分
流形上的散度公式证明数值模型和式极限两种类型积分的新公式关联
工程意义上的流形积分解析积分值任意精度浮点数积分值
中图分类号:
O17/O412.3
目录
引言1...........................................................2
引言2证明的前提条件—-单连通可定向闭合曲面坐标系的建立...................4
流形上的散度公式证明...............................................10
总结...........................................................13
参考书籍........................................................14
引言1
散度公式(又称Остроградский-Gauss公式)是现代数学、物理体系的核心公式之一[1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16][17].传统的散度公式证明逻辑体系,建立了基于空间直角坐标系投影法(简称投影法)的曲面积分与三重积分的公式关联,确立了投影法为曲面积分的根本方法.投影法的基本思路是将三维欧氏空间区域中的曲面积分,转化为某一空间直角坐标平面上的二重积分,以间接的方式达到目的.
投影法的缺陷是明显的:
第一,积分曲面在某一坐标平面的投影区域不能有重叠,这就决定了积分曲面只能是非常简单的函数曲面;在现实世界和物理、工程领域更为普遍存在复杂参数曲面,投影法则无能为力;
第二,投影法通常要求积分曲面具有某种对称性(点对称、轴对称和面对称等[2]),计算诸如"以三维坐标原点为中心的球面上侧、下侧、左侧、右侧曲面"类型的简单曲面积分,再乘以某一常数,得到整个球面的积分值;在现实世界和物理、工程领域更为普遍存在的不对称、不规则曲面,用投影法计算非常繁琐,在绝大多数情况下不能计算[1][2][3][4][5][6][7]
[8][9];以致于物理、工程领域的许多重要问题(例如电磁学领域的Maxwell方程组实例化和流体力学领域的任意不规则控制面积分)的解决途径,均建立在直角坐标系或其它坐标系(例如贴体坐标系等)的偏微分方程组求解(例如有限元法、边界元法、有限差分法等[18])基础上.一个多世纪以来的数学、物理、工程实践已经证明,通过投影法、直角坐标系或其它坐标系的偏微分方程组,难于甚至不能获得关于复杂几何对象(流形)的解析解、数值解;
第三,因不同积分曲面的几何差异,投影的方向、投影的次数千差万别[尤其是分面投影法].有100计算实例,就可能有100种投影方案.计算过程不可能标准化、模块化,不利于电子计算机编程;
第四,不论积分曲面复杂程度,投影法实际计算过程普遍繁琐;
第五,在物理、数学分析领域至关重要的Остроградский–Gauss公式,Stokes公式(在某种意义上也包括Green公式),投影法几乎没有直接计算实例(即使有,也是极个别的特例,没有代表性.如正方体、长方体表面外观的闭合“曲面”[12][13][15][16][17],数条正方体不同表面截线段围成的闭合“曲线”[12][13][15][16],平面“x+y+z=1”与三个直角坐标平面的相交三角形构成的闭合“曲线”[12][13][14][15][16]等).通常的思路是,先用符号逻辑的方式证明在直角坐标系中三大公式的存在,然后是如何应用这三大公式简化计算.非常遗憾、困惑的是没有这三大公式的丰富而绚丽的直接计算实例.
在球面坐标系,有向量场参数曲面积分[13]、球体空间区域三重积分[即通过三阶Jaccobi行列式变量变换];在极坐标系,有平面区域二重积分[即通过二阶Jaccobi行列式变量变换]等计算方法[12][13][14][16][17].但是存在下列问题:
第一,在球面坐标系内,向量场闭合曲面积分与数量场闭合空间区域三重积分彼此孤立,没有通过散度公式关联,并且两者的计算结果不能相互验证;
第二,向量场闭合曲面积分与数量场闭合空间区域三重积分局限于正交曲线坐标系[即球面坐标系和柱面坐标系][12][13][14][16][17]或广义球面坐标系[14],没有扩展到无穷多个任意参数曲面坐标系;
第三,在极坐标系内,向量场环路积分与数量场平面区域二重积分彼此孤立,没有通过Green公式关联,并且两者的计算结果也不能相互验证;
第四,数量场[二元函数]闭合平面区域二重积分局限于正交曲线坐标系[即极坐标系][11],没有扩展到无穷多个任意单连通闭合曲线坐标系.
传统的流形微积分学,用外微分形式推导出Green公式,Остроградский-Gauss公式,Stokes公式,乃至关于n维空间积分的广义Stokes公式[20],即
但是这类用外微分形式推导出的公式只具有抽象的理论意义,并没有揭示积分的具体实现过程,更无具体数值模型可言.
认识,无止境.
引言2证明的前提条件
——单连通可定向闭合曲面坐标系的建立
(一)
考察证明的对象---散度公式:
“散度公式设空间闭区域Ω是由光滑或分片光滑的闭曲面S围成,函数
P(x,y,z),Q(x,y,z),R(x,y,z)[构成向量场A]及其偏导数在空间闭区域Ω上连续,则
(1)
其中曲面S为空间闭区域Ω的整个边界曲面外侧,n为曲面S的单位外法向量,divA为向量场A的散度”.
在公式的定义中,强调空间闭区域Ω的边界闭合曲面S必须是能够区分其”内侧”、”外侧”的可定向曲面.
在传统的直角坐标系Остроградский-Gauss公式证明中,”抽象可定向闭合曲面∑”是这样定义的:
抽象可定向闭合曲面∑由三个子曲面
分片包围而成,其中曲面
皆为抽象二元函数.
(参见《高等数学(第六版)》(下册)同济大学数学系高等教育版2007P168-170)
也就是说,散度公式客观上要求,不论在空间直角坐标系,或者在其它坐标系,被证明的相关曲面S必须具有两种属性:
(1)闭合性;
(2)可定向性.
离开传统的空间直角坐标系,怎样刻画抽象的、具有普遍意义的”可定向闭合曲面”并且进一步建立”可定向闭合曲面坐标系”?
并没有现成的答案.
Poincare猜想[19]断定"任何与n维球面同伦的n维闭合流形必定同胚于n维球面",在散度或旋度公式涉及的三维欧氏空间,其对应的判断为"任何单连通、可定向2维闭合流形必定同胚于2维球面".
也就是说,根据Poincare猜想,在散度或旋度公式涉及的三维欧氏空间,任何单连通、可定向的闭合曲面(虽然仅仅是单连通),不论其几何外观如何千变万化,必定有同胚于”球面”这一普遍属性.
进一步的问题自然是”在三维欧氏空间,能否根据Poincare猜想这一普遍属性,定义单连通、可定向的闭合曲面的抽象的、普遍意义的表达式?
”这也正是本”引言”讨论的中心内容.
在空间解析几何学中,上述"2维球面"的参数表达式为[sin(u)cos(v),sin(u)sin(v),
cos(u)],其中参数u的变化范围[0,Pi],参数v的变化范围[0,2*Pi](在严格意义上,该参数表达式是”2维球面”在”空间直角坐标系”和”球面坐标系”之间的转换式;”二维球面”在球面坐标系的表达式是常数1).
在拓扑学领域,"同胚"的定义为"两个流形,如果可以通过弯曲、延展、剪切等操作把其中一个变为另一个,则认为两者是同胚的".
从解析几何学和拓扑学的角度再理解Poincare猜想,既然"2维球面"的参数方程为[sin(u)cos(v),sin(u)sin(v),cos(u)],其中参数变化范围u[0,Pi],v[0,2*Pi],则其变形[a*sin(u)cos(v),b*sin(u)sin(v),c*cos(u)],u[0,Pi],v[0,2*Pi]
(其中待定系数a,b,c为任意非零常数)即为任意椭球面的参数方程.在三维欧氏空间,任意椭球面皆同胚于球面,这是拓扑学的常识,无需讨论;
如果a,b,c为任意"一阶可导连续函数",可能出现怎样的情况?
参见下列图例:
图例1:
假设任意待定系数a=sin(u)+cos(v),b=cos(u),c=cos(v/2),则目标参数曲面[a*sin(u)*cos(v),b*sin(u)*sin(v),c*cos(u)](其中u∈[0,π],v∈[0,2π])
为[(sin(u)+cos(v))*sin(u)*cos(v),cos(u)*sin(u)*sin(v),
cos(v/2)*cos(u)](其中u∈[0,π],v∈[0,2π])
其实际参数图形为:
图例1由待定系数a,b,c输入任意"一阶可导连续函数",
输出参数曲面呈非单连通、非闭合状态,
与”Poincare猜想”及”流形上的散度或旋度公式”讨论的内容无关
图例2:
假设任意待定系数a=sin(u+v)+cos(v),b=cos(v),c=cos(v/2),
则目标参数曲面[a*sin(u)*cos(v),b*sin(u)*sin(v),c*cos(u)](其中u∈[0,π],v∈[0,2π])
为[(sin(u+v)+cos(v))*sin(u)*cos(v),cos(v)*sin(u)*sin(v),
cos(v/2)*cos(u)](其中u∈[0,π],v∈[0,2π])
其实际参数图形为:
图例2由待定系数a,b,c输入任意"一阶可导连续函数",
输出参数曲面呈非单连通、不可定向状态,
与”Poincare猜想”及”流形上的散度或旋度公式”讨论的内容无关
图例3:
假设任意待定系数a=sin(u),b=cos(u)+cos(u+3*v)/3,c=cos(u),则目标参数曲面[a*sin(u)*cos(v),b*sin(u)*sin(v),c*cos(u)]
(其中u∈[0,π],v∈[0,2π])
为[sin(u)*sin(u)*cos(v),(cos(u)+cos(u+3*v)/3)*sin(u)*sin(v),
cos(u)*cos(u)](其中u∈[0,π],v∈[0,2π])
其实际参数图形为:
图例3由待定系数a,b,c输入任意"一阶可导连续函数",
输出曲面呈单连通可定向闭合状态
可以作为”Poincare猜想”及”流形上的散度或旋度公式”讨论的对象
实验数据从原始现象表明,同样属于参数曲面[asin(u)cos(v),bsin(u)sin(v),
ccos(u)],u[0,Pi],v[0,2*Pi],因待定系数a,b,c的不同取值,一部份曲面属于单连通、可定向闭合曲面,一部分曲面则例外.
也就是说,参数曲面[asin(u)cos(v),bsin(u)sin(v),ccos(u)],u[0,Pi],v[0,2*Pi]存在两种情况:
(1)在待定系数a,b,c为任意非零常数的情况下,参数曲面为椭球面(自然同胚于球面);
(2)在待定系数a,b,c为任意一阶可导连续函数的情况下,参数曲面可以为单连通可定向闭合曲面(同胚于球面),也可以为非单连通可定向闭合曲面(不同胚于球面).
进一步的问题自然是“在参数曲面[asin(u)cos(v),bsin(u)sin(v),ccos(u)],u[0,Pi],v[0,2*Pi]模式中,能否通过某种定义将非单连通可定向闭合曲面(不同胚于球面)的情况排除?
”
(二)
设定“任意曲面”为一集合,则“任意单连通、可定向闭合曲面”是前者的子集合.Poincare猜想是这一子集合的属性,本论文“流形上的散度或旋度公式证明”及其”和式极限证明”则讨论散度或旋度公式是否适用于这一子集合.
Poincare猜想为用参数方程方法描述“任意单连通、可定向闭合曲面”的某种属性(即
“同胚于2维球面”这一属性)提供了实现途径.参考丘成桐院士2006年观点:
庞加莱猜想和三维空间几何化的问题是几何领域的主流,它的证明将会对流形性质的认识,甚至用数学语言描述宇宙空间产生重要影响.
基于上述情况,将无数具体的单连通、可定向闭合曲面抽象化为一个统一的表达式:
[asin(u)cos(v),bsin(u)sin(v),ccos(u)],u[0,Pi],v[0,2*Pi]
(其中待定系数a,b,c不能任意指定,而必须服从于曲面的”单连通、可定向闭合”的拓扑学属性)
也就是说,如果待定系数a,b,c能够任意指定,则目标曲面[a*sin(u)cos(v),
b*sin(u)sin(v),c*cos(u)],u[0,Pi],v[0,2*Pi]可能是”单连通、可定向闭合曲面”,也可能不是;
如果预先设定目标曲面[asin(u)cos(v),bsin(u)sin(v),ccos(u)],
u[0,Pi],v[0,2*Pi]本身就是”单连通、可定向闭合曲面”,则待定系数a,b,c就不能任意指定了.
从几何意义解释上述现象---在空间直角坐标系,球面(即[sin(u)cos(v),
sin(u)sin(v),cos(u)],u[0,Pi],v[0,2*Pi])沿x,y,z轴三个方向任意连续变化(即[asin(u)cos(v),bsin(u)sin(v),ccos(u)],u[0,Pi],v[0,2*Pi],其中待定系数a,b,c为任意一阶可导连续函数),不一定产生单连通、可定向闭合曲面;
反过来,在空间直角坐标系,任一单连通、可定向闭合曲面--必定由球面(即[sin(u)cos(v),sin(u)sin(v),cos(u)],u[0,Pi],v[0,2*Pi])沿x,y,z轴三个方向连续变化而成(也必定能够沿x,y,z轴三个方向连续变回球面)--Poincare猜想为依据.
例如,正六面体(即正方体)表面也可以被视为单连通、可定向闭合曲面---但是正六面体表面难于甚至不能用参数方程描述---但是不能否认,根据Poincare猜想,正六面体表面必定同胚于球面,必定由球面(即[sin(u)cos(v),sin(u)sin(v),cos(u)],u[0,Pi],
v[0,2*Pi])沿x,y,z轴三个方向连续变化而成(也必定能够沿x,y,z轴三个方向连续变回球面);根据Poincare猜想,正六面体表面同样可以用[asin(u)cos(v),bsin(u)sin(v),
ccos(u)],u[0,Pi],v[0,2*Pi]参数模式描述.
用[asin(u)cos(v),bsin(u)sin(v),ccos(u)],u[0,Pi],v[0,2*Pi]模式描述抽象的、具有普遍意义的单连通、可定向闭合曲面,实际上是用Poincare猜想来描述单连通、可定向闭合曲面的某种内在结构和属性(即同胚于球面这一属性),为进一步的公式推导设定一个恰当的前提条件.
需要特别指出,在具体曲面为“复连通、可定向闭合曲面”(例如环面及其同胚曲面)情况下还不能实现抽象化,还没有相关的理论依据为支持.
在实际操作层面,用Plot3D[属于WaterlooMaple计算机代数系统指令]指令绘画出某一参数曲面,必须在直观视觉上判定该曲面是否为单连通、可定向闭合曲面以后,才能决定是否适用于流形上的散度或旋度公式数值模型;从参数表达式本身无法判断曲面是否为单连通、可定向闭合曲面;
“参数曲面是否为单连通、可定向闭合曲面”的决定因素在拓扑学领域而不在解析几何领域;单凭解析几何的参数方程方法并不能够推导、演绎出某一曲面的单连通、可定向闭合属性.
(三)
[asin(u)cos(v),bsin(u)sin(v),ccos(u)],u[0,Pi],v[0,2*Pi]只是基于Poincare猜想定义的抽象的、普遍意义的单连通可定向闭合曲面表达式,还不属于坐标系;抽象单连通可定向闭合曲面坐标系为[rasin(u)cos(v),rbsin(u)sin(v),
rccos(u)],r[0,∞],u[0,Pi],v[0,2*Pi],其中r为向径,a,b,c为待定系数(因为a,b,c既可以为非零常数,也可以为一阶可导连续函数),具有不确定性.
实际上,抽象单连通可定向闭合曲面表达式[asin(u)cos(v),bsin(u)sin(v),
ccos(u)],u[0,Pi],v[0,2*Pi]与椭球面表达式[asin(u)cos(v),bsin(u)sin(v),
ccos(u)],u[0,Pi],v[0,2*Pi]在形式上是完全一致的,只是两者对待定系数a,b,c的解释不同:
前者将a,b,c解释为"任意非零常数或一阶连续可导函数(非任意,受曲面的单连通可定向闭合属性限制)”,而后者将a,b,c解释为只是"任意非零常数";故抽象单连通可定向闭合曲面坐标系与直角坐标系的对应关系是x=r*asin(u)cos(v),y=r*bsin(u)sin(v),z=r*ccos(u)(与椭球面坐标系-直角坐标系转换式是相同的).
流形上的散度公式证明:
散度公式设空间闭区域Ω是由光滑或分片光滑的闭曲面S围成,函数
P(x,y,z),Q(x,y,z),R(x,y,z)[构成向量场A]及其偏导数在空间闭区域Ω上连续,则
(1)
其中曲面S为空间闭区域Ω的整个边界曲面外侧,n为曲面S的单位外法向量,divA为向量场A的散度
证明:
定义任意单连通、可定向闭合曲面S的参数表达式:
[asin(u)cos(v),bsin(u)sin(v),ccos(u)]
(2)
其中a,b,c为非零常数或一阶可导连续函数表达式,单连通、可定向闭合曲面S决定a,b,c的取值;设定参数u,v的变化范围[0,
],[0,2
],使曲面S闭合.(参见Poincare猜想:
"任何与n维球面同伦的n维闭合流形必定同胚于n维球面")[19]
根据曲面参数表达式
(2),定义并计算第一偏导数矩阵,获取曲面S的切平面法向量:
=
(3)
从(3)式分别提取i,j,k项系数,获得曲面S的切平面法向量:
[
](4)
向量场A与曲面S的切平面法向量(4)的空间点积对曲面参数u,v的积分:
(5)
将曲面S的参数表达式
(2)各项通乘以向径r(设定r>0),将x,y,z轴方向上的曲面坐标参数转化为空间区域坐标参数:
[rasin(u)cos(v),rbsin(u)sin(v),rccos(u)](6)
根据空间区域坐标参数(6),定义并计算第二偏导数矩阵,获取空间闭区域Ω微元系数的一般表达式:
=
(7)
计算向量场A的散度,并将其从直角坐标形式(8)转变为空间闭区域Ω坐标形式(9):
(8)
(9)
散度(9)与空间闭区域Ω微元的乘积对变量r,u,v的三重积分:
(10)
其中,
即设定空间闭区域Ω微元本身对参数r,u,v的三重积分不能为零,也可以理解为设定空间闭区域Ω不能为零体积
即(5)式=(10)式:
=
亦可表述为
(1),证毕
环面散度公式证明和流形上的散度公式数值模型,参见”附件1流形上的散度公式证明和数值模型(分析与说明)”;
流形上的散度公式和式极限证明及其数值模型,参见”流形上的散度公式和式极限证明和数值模型[附件3分析与说明]”
总结
传统的散度公式证明逻辑体系,建立了基于空间直角坐标系投影法(简称投影法)的曲面积分与三重积分的公式关联,确立了投影法为曲面积分的根本方法.但是投影法存在诸多明显的缺陷(例如计算过程繁琐,不适用于不对称、不规则曲面等),以致于物理、工程领域的许多重要问题(例如电磁学领域的Maxwell方程组实例化和流体力学领域的任意不规则控制面积分)的解决途径,均建立在直角坐标系或其它坐标系的偏微分方程组求解基础上.一个多世纪以来的数学、物理和工程实践已经证明,通过投影法、直角坐标系或其它坐标系的偏微分方程组,难于甚至不能获得关于复杂几何对象(流形)解析解、数值解;
传统的流形微积分学,用外微分形式推导出Green公式,Остроградск
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 流形 公式 证明
![提示](https://static.bdocx.com/images/bang_tan.gif)