different PDEs to solve through events in pdepe
3 次查看(过去 30 天)
显示 更早的评论
There is a model defined by a pointwise PDE solving for u(x,t), i.e. it is PDE1 when u(x,t)>0 and PDE2 when u(x,t)<=0. Can one use pdepe events to solve this pointwise PDE? I don't think the following code works because once umesh=0 anywhere the whole PDE will be switched to the other (e.g. PDE2 on the global domain rather than locally will be triggered if PDE1 is terminated by events).
function [value, isterminal, direction] = pdevents(m,t,xmesh,umesh)
value = umesh;
isterminal = ones(size(umesh));
direction = zeros(size(umesh));
4 个评论
Torsten
2024-2-8
Then you should try to manage these discontinuities within the function that defines the PDE equation, not by using the "events" facility.
"Events" might be useful in this context to stop an integration (e.g. if the PDE has reached a steady state or if some value has reached a certain limit level).
回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 PDE Solvers 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!