TimeDependentResults
时间PDE解决方案和派生的数量
描述
一个TimeDependentResults对象包含PDE的解决方案及其梯度的形式方便策划和后处理。
一个
TimeDependentResults对象包含解决方案及其梯度计算节点的三角形或四面体网格,生成的generateMesh。出现在解决方案值的节点
NodalSolution财产。解决方案*中出现
SolutionTimes财产。解决方案的梯度值的三个组件的节点出现在
XGradients,YGradients,ZGradients属性。数组的尺寸
NodalSolution,XGradients,YGradients,ZGradients使您能够提取解决方案和梯度值指定的时间指标,和方程指数在PDE系统。
插入解决方案或其梯度定义网格(例如,指定的meshgrid),使用interpolateSolution或evaluateGradient。
创建
有几种方法可以创建一个TimeDependentResults对象:
解决时间问题使用
solvepde函数。这个函数返回一个PDE作为一个解决方案TimeDependentResults对象。这是推荐的方法。解决时间问题使用
抛物线或双曲函数。然后使用createPDEResults函数来获得一个TimeDependentResults返回的对象从PDE的解决方案抛物线或双曲。请注意,抛物线和双曲遗留功能。他们不建议解决PDE的问题。
属性
对象的功能
evaluateCGradient |
评估通量PDE的解决方案 |
evaluateGradient |
评估在任意点梯度PDE的解决方案万博 尤文图斯 |
interpolateSolution |
插入PDE解决任意点 |
例子
一个抛物问题的解决方案
解决一个抛物线的二维几何问题。
创建和查看几何图形:正方形和一个圆域。
%广场集中在(1)通过rect1 = [3、4 0; 2; 2。0, 0, 0; 2; 2);%的圆心在(1.5,0.5)circ1 =(1; 1.5;综合成绩;0.25);%添加额外的零循环circ1 = [circ1; 0(长度(通过rect1)长度(circ1), 1)];gd =(通过rect1, circ1);ns = char (“通过rect1”,“circ1”);ns = ns ';科幻小说=“通过rect1 + circ1”;(dl, bt) = decsg (gd、科幻、ns);pdegplot (dl,“EdgeLabels”,“上”,“FaceLabels”,“上”)轴平等的ylim ([-0.1, 2.1])
在PDE模型包括几何。
模型= createpde ();geometryFromEdges(模型、dl);
设定边界条件,较高的温度和左边缘10。
applyBoundaryCondition(模型,“边界条件”,“边缘”(2、3),“u”10);
设置初始条件,正方形在温度0,和圆在温度100。
setInitialConditions(模型中,0);100年setInitialConditions(模型,“面子”2);
定义模型的系数。
specifyCoefficients(模型,“m”0,“d”,1“c”,1“一个”0,“f”,0);
解决问题的步骤0.01中通过1/2 * 0。
generateMesh(模型,“Hmax”,0.05);tlist = 0:0.01:0.5;结果= solvepde(模型、tlist);
情节的解决方案* 0.02,0.04,0.1和0.5。
索尔= results.NodalSolution;次要情节(2 2 1)pdeplot(模型,“XYData”索尔(:3)标题(“0.02”次要情节(2,2,2)pdeplot(模型,“XYData”索尔(:5)标题(“0.04”次要情节(2,2,3)pdeplot(模型,“XYData”索尔(:11)标题(“0.1”次要情节(2,2,4)pdeplot(模型,“XYData”索尔(:51)标题(“0.5”)
版本历史
介绍了R2016a
