解决平流扩散pde

43岁的观点(30天)

显示旧的评论

Deepa Maheshvare 2018年12月25日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/437177-solving-advection-diffusion-pde

评论道: Deepa Maheshvare 2019年1月2日

我想解决上述pde给定的边界和初始条件。我遇到pdepe在MATLAB函数。

我有机会看看这个例子在这里。我不能理解pdex1pde函数定义为我的情况。

有人能帮忙吗?

1评论
显示没有隐藏没有

Deepa Maheshvare 2018年12月27日

你好,

以下是实现正确吗?

                               函数DiffusionConvection
                              
                               m = 0;
                              
                               x = linspace (0、62、10);
                              
                               t = linspace (0, 10100);
                              
                               索尔= pdepe (m, @pdefun, @icfun @bcfun x, t)
                              
                               函数[g、f、s] = pdefun (x、t、c、DcDx)
                              
                               D = 900;
                              
                               v = 10;
                              
                               g = 1;
                              
                               f = D * DcDx;
                              
                               s = - v * DcDx;
                              
                               结束
                              
                               函数c0 = icfun (x)
                              
                               c0 = 80;
                              
                               结束
                              
                               函数[pl, ql,公关,qr] = bcfun (xl, cl, xr、铬、t)
                              
                               pl = cl -10;
                              
                               ql = 1;
                              
                               公关= cr;
                              
                               qr = 1;
                              
                               结束
                              
                               结束

登录置评。

在回答这个问题。

答案(1)

比尔•格林 2018年12月27日

1
链接

直接链接到这个答案

//www.tianjin-qmedu.com/matlabcentral/answers/437177-solving-advection-diffusion-pde answer_354140

唯一的问题,我看到你的代码在边界条件。我纠正你 bcfun 功能和附加我的下面的代码的版本。

                             函数DiffusionConvection
                            
                             m = 0;
                            
                             x = linspace (0、62、10);
                            
                             t = linspace (0, 10100);
                            
                             D = 900;
                            
                             索尔= pdepe (m, @pdefun, @icfun @bcfun x, t);
                            
                             函数[g、f、s] = pdefun (x、t、c、DcDx)
                            
                             v = 10;
                            
                             g = 1;
                            
                             f = D * DcDx;
                            
                             s = - v * DcDx;
                            
                             结束
                            
                             函数c0 = icfun (x)
                            
                             c0 = 80;
                            
                             结束
                            
                             函数[pl, ql,公关,qr] = bcfun (xl, cl, xr、铬、t)
                            
                             pl = -10 * D;
                            
                             ql = 1;
                            
                             公关= 0;
                            
                             qr = 1;
                            
                             结束
                            
                             结束