PDEPE solver with x-dependent diffusion constant

2 次查看（过去 30 天）
Thomas
Thomas 2015-3-17
编辑： Ronald 2016-5-4
Hello,
according to the documentation of the PDEPE solver one can use diffusion constant values that are dependent on x. However, if I calculate a simple 1 D model and set the c value with an if statement in the solver, I find that the system converges to entirely different values than it did before. This should not be the case, the diffusion should go slower in the more rigid parts, but the endvalues should be exactly the same.
This also happens if I do not use an if statement buf for example c=sin(x). Is this a major bug? How can i cirumvent this?
Greetings, Thomas

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

采纳的回答

Torsten
Torsten 2015-3-18
The diffusion coefficient is under the d/dx-operator - so how can you put it in the c-array ?
You must define f=Diff*DuDx where Diff is the diffusion coefficient (which may depend on x).
If this does not answer your question, we will have to see the MATLAB code to give further advice.
Best wishes
Torsten.
  2 个评论
Thomas
Thomas 2015-3-18
Thank you, that was exactly what was going wrong. Splendid!
Ronald
Ronald 2016-5-4
编辑：Ronald 2016-5-4
Torsten: If I understand it corrctly, the constant Diff under the d/dx-operator refers to diffusivity (physically)? Since in d^2/dx^2 T - c d/dt T=0, the constant 1/c refers to diffusivity.
Nevertheless it should be possible to make the parameter c in the pdepe framework time-dependent, right? Do you maybe know how to realize this? Because the solver supposes c has the diagonal elements for each c_i of the pde-system, it is quite confusing.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Mathematics and Optimization 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by