PDE with discontinous flux function
1 次查看(过去 30 天)
显示 更早的评论
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
0 个评论
回答(1 个)
另请参阅
类别
在 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!