清空筛选条件
清空筛选条件

Neumann boundary condition in a first order PDE

4 次查看（过去 30 天）
Antonio
Antonio 2012-9-25
I'm trying to solve the following equation using PDEPE:
dC/dt + v * dC/dx = constant
With the boundary conditions:
C(t,0)=Cin dC(t,L)/dx=0
My question is how can I incorporate the second BC in the PDEPE syntax, if I should define f = [-v]. Is there any posibility to call the penultime value and make it equal to ur, so dC/dx=0? u(x_n) = u(x_(n-1))
Thanks for your cooperation!
Antonio

请先登录，再回答此问题。

回答（1 个）

Tom
Tom 2012-9-25
编辑：Tom 2012-9-25
In this case you want to set
pr = 0;
qr = -1/v; %to cancel out f
  11 个评论
显示 9更早的评论
Tom
Tom 2012-9-25
I'm not sure how PDEPE can deal with two time terms, maybe it would be better to solve numerically in a for loop.
Antonio
Antonio 2012-9-25
About the last term (dq/dt), ill solve it differently. My problem now is this bc!

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Partial Differential Equation Toolbox 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by