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

Answer 1

[Moritz]

Well, RTFM (Read the f... Manual).

So it is all described in the PDE Manual at 2.16 i guess.