PDE with discontinous flux function

Question

Moritz 2013-1-28

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/60095-pde-with-discontinous-flux-function

Hi,

i have a strongly degenerate quasilinear PDE (sedimentation), describing the change of concentration depending on time and radius. The second order PDE degenerates to a hyperbolic PDE at certain concentrations.

dc/dt=d(w^2*r/g*f(c))/dr+d(dA(c)/dr)/dr

where c..concentration t...time w..omega (angular velocity) r..radius f..flux function (Richardson Zaki; discontinuous) A..primitive of the diffusion coefficient a (which itself is discontinuous) a.. power law function

considering the PDE pdepe solves, f and s are 0 at certain concentration intervals.

I do not have a strong mathematical background and will discuss this problem during a math.seminar but until then i am trying to figure out how deep i would have to go into discretization (there are quite a few numerical methods published) or if a matlab solver is applicable.

I considered the question from zhao qingyuan "How to set the coefficient of PDE equation as a user-defined matlab function?"

would this be the way to go ? To define an if else condition ?

Thank you for your help