导体内部的稳定温度分布大物DTP.docVIP

导体内部的稳定温度分布 严鹏熙 （物理与电信工程学院 物理学三班 20082301105） 引言： 热学中一个物体的温度处于温恒态时各部分温度按一定规律分布，而导体内部的温度分布通常不可直接测量，有时我们需要掌握导体内部温度变化规律，则可引入勒让德函数来建立这个物理模型，只要测出导体表面的温度，导体内部的温度分布便可推知，同时我们可以画出导体内部温度分布变化图形，以此来形象表达该模型。 模型：一个半球形热良导体，在球坐标系（r,θ,ω）下，其表面维持温度为函数 ，底面维持为0，求这个导体内部的稳定温度分布。 解 球体表面温度分布函数 与ω无关，因此球体温度分布函数与ω无关，于是u满足下列Dirichlet边值问题 ， 由半球底面上的边界条件以及方程关于θ的奇偶不变性，可以对关于θ=做奇延拓为新未知函数，变成球面上的定解问题 . 设，作代换x=，并记为，分析可知R，P分别满足Euler方程和Legendre方程 它们有物理意义的特解分别为和P=。由叠加原理，可以设方程有物理意义的级数解为 在变换x=下，边界条件可以改写为。因此 又是n次多项式，而等式左边的次数为2，因此，当n时，。由legendre函数的正交关系可以求出 n=0，1，2. 是奇函数，所以。而 因此 从而 该半球形良导体内部的稳定温度分布解得图像如下： 结论：半球形导体的温度分布具有由球顶点向四周递增的特点。 参考文献：《数学物理方程与特殊函数》 华中科技大学 李元杰 《数学物理方程》 桂子鹏 康盛亮 编程程序： #include WizDemo.h #include math.h float step; Text title; float P_radius, P_omega,V; // 构造函数 demoApp::demoApp() : WizGLApp() { viewMode = VIEW_PERSPECTIVE2; // 观察视角模式 demoState = DEMO_STOP; // 初始动画状态 wndWidth = 1024; // 初始窗口宽度 wndHeight = 715; // 初始窗口高度 showParamWnd = true; // 初始是否显示参数窗口 showInfoWnd = false; // 初始是否显示信息窗口 showToolbar = true; // 初始是否显示工具栏 // wiz::Assign(bgColor, 0.32f, 0.72f, 0.54f, 0.0f); // 设置背景色 } // 析构函数 demoApp::~demoApp() { } double BesselFun(int n, float x) { double Result = 0.0f; double theta; int m, count; count = (abs(n)+1)*360; for (m=0; m=count; m++) { theta = m*PI/count; Result += cos(x*sin(theta) - abs(n)*theta); } Result /= count; if (n0) Result *= pow(-1, abs(n)); return Result; } double LegendreFun(int n, float x) { double Result = 0.0f; double den, acc; int i,m; if (n0) return Result; den=acc=1.0f; for (m=0; m=n/2; m++) { for (i=2*(n-m); i0; i--) den *= i; for (i=n-m; i0; i--) acc *= i; for (i=n-(2*m); i0; i--) acc *= i; for (i=m; i0; i--) acc *= i; Result += pow(-1, m) * den * pow(x, n-(2*m)) / (pow(2, n) * acc); } return Result; } /////////////////////////////////////////////////////////////////////// // 在此创建场景 void demoApp::InitializeScenes()