Solve PDE system - interface problem between two domains

5 次查看（过去 30 天）
MFV
MFV 2017-5-2
评论： Torsten 2023-12-12
I am trying to solve a PDE system in one space dimension using pdepe function. I have two equations, each of them applied to different domains and I need to calculate the time evolution in the interface, I am attaching the equations, boundary conditions and initial conditions.
I do not see how could I define a boundary condition in the interface and also different domains for each equations. Is it possible to do it in this way? If not, could anybody indicate me an alternative to solve the problem?
Thanks
  1 个评论
Bill Greene
Bill Greene 2017-5-3
No boundary or interface condition is needed where the material properties are discontinuous. The first argument to pdefun is x; you simply evaluate the coefficients based on the value of x. The pdepe documentation discusses this issue specifically: "Discontinuities in c and/or s due to material interfaces are permitted provided that a mesh point is placed at each interface."

请先登录，再进行评论。

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

回答（1 个）

Torsten
Torsten 2017-5-3
编辑：Torsten 2017-5-3
If the condition you state at x=delta is correct, you can solve the two PDEs one after the other since u1 and u2 are not coupled. I assume that u_i(t) is a function that is known a priori.
Best wishes
Torsten.
  5 个评论
显示 3更早的评论
Torsten
Torsten 2017-5-4
编辑：Torsten 2017-5-4
If those two conditions are your transmission conditions, you can proceed as Bill Greene suggests: You can solve both equations in the union of Domain 1 and Domain 2 and set the PDE coefficients depending on the actual x-position:
if x<=Delta
c=...;
s=...;
f=...;
else
...
end if
One grid point should be placed at x=Delta.
The transmission conditions will be automatically satisfied.
Best wishes
Torsten.
Torsten
Torsten 2023-12-12
In your code, I guess f is your function, but what about c and s?
No. c, f and s is the usual notation for parts of the PDE equation that is to be solved by "pdepe".
Take a look at the documentation under
https://uk.mathworks.com/help/matlab/ref/pdepe.html
Example
"Solve Heat Equation in Cylindrical Coordinates"
Section
Code Equation

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Boundary Conditions 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by