混合边界条件的Poisson方程实例.docVIP

带混合边界条件的Poisson方程 现在我们考虑带混合边界条件的Poisson方程问题: 对于给定的函数 和,求得函数使得 (1) 其中,该问题的解在存在奇点. 预备知识: 如当求解区域为时, 的奇点就为,并且可以表示成如下形式 . 其中为常数, 为极坐标系参数,我们有极坐标变换 实例演算: 设其中为精确解.运用Freefem++解上述问题. 程序分析如下: border N(t=0,1) { x=-1+t; y=0; label=1; }; border D1(t=0,1){ x=t; y=0; label=2;}; border D2(t=0,1){ x=1; y=t; label=2; }; border D3(t=0,2){ x=1-t; y=1; label=2;}; border D4(t=0,1) { x=-1; y=1-t; label=2; }; //形成五条边构成边界(如下图所示) mesh T0h = buildmesh(N(10)+D1(10)+D2(10)+D3(20)+D4(10)); //给定一种划分,为了和后面的调整型划分区分开记此划分为T0h. plot(T0h,wait=true); //画出该三角划分,见图一. fespace V0h(T0h,P1); //形成P1元有限元空间. V0h u0, v0; //定义有限元空间中的两个函数. func f=-2*30*(x^2+y^2); //定义已知的函数. func us = sin(atan2(y,x)/2)*sqrt( sqrt(x^2+y^2) ); //根据表达式 定义. real K=10.; //有已知函数的表达式可知K=10. func ue = K*us + 30*(x^2*y^2); // 根据表达式定义 V0h fh=ue; //定义V0h空间的一个函数为真解表达式 plot(fh,wait=true); //画出真解ue的图,见图四(为了与数值解图三作比较). solve Poisson0(u0,v0) =int2d(T0h)( dx(u0)*dx(v0) + dy(u0)*dy(v0) ) //双线性部分 - int2d(T0h)( f*v0 ) //线性部分 + on(2,u0=ue) ; //边界条件 mesh Th = adaptmesh(T0h,us); //对原来的划分进行了调整,调整后形成的新划分能更好的处理解出现局部震荡(陡峭)的问题,在电子书5.5节有详细介绍. for (int i=0;i 5;i++) //通过for循环生成调整5次后的新的划分,使数值 { 解能更好的逼近真解(后面通过误差估计比较 mesh Th=adaptmesh(Th,us); 能体现这一点). } ; fespace Vh(Th, P1); //在新的划分下形成新的有限元空间Vh. Vh u, v; //在空间Vh中定义两个函数. solve Poisson(u,v) =int2d(Th)( dx(u)*dx(v) + dy(u)*dy(v) ) //双线性部分 - int2d(Th)( f*v ) //线性部分 + on(2,u=ue) ; //边界条件 plot(Th,wait=true,ps=adaptDNmix.ps); //画出新的三角划分,见图二. plot(u,wait=true); //画出所得数值解,见图三. Vh uue = ue; //定义Vh空间一个函数(真解