Coupled ODE to PDE with Method of lines
1 次查看(过去 30 天)
显示 更早的评论
Dear all,
i do have a PDE which describes sedimentation. I do discretize it in space and end up with a system of ode´s. Fairly simple.
Now i use the same PDE (settling of particles) inside a syringe (no particles leave the syringe) and move the piston (pressure filtration). Therfore i do have a moving boundary at the piston. So i have to solve for the piston height too.
The system looks like
dh/dt=(height(n+1)-height(n))/dt
dy/dt= dh/dt*(something discretized) + h(n)/h(n+1)*(something discretized)
My question is if i can put dh/dt and dy/dt in a vector and treat it as a coupled system of ode´s ? Could someone give me a hint how to code that ?
E.g
i=1:N
dydt=dydt(N+1)*(something discretized) + ???*(something discretized)
i=N+1
dydt=myOdePistonheight
So if the height variable is part of my input array y0 how do i get h(n+1)/h(n)
Sincerly
Moritz
回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Boundary Conditions 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!