with pdepe, how to solve coupled ode-pde with ode citing pde's BC

3 views (last 30 days)
Hi everyone,
I have coupled pde-ode equations, where pde's one boundary condition is ode's input parameter. I want to write both equations as PDE so that I can use pdepe, but cannot figure out how to get boundary condition of pde and pass that to ode in pdefun. Could anyone tell me how to code it in order to use pdepe? Thanks!
for example, on x from 0 to 1, i have
dC1/dt = d^2 C1 / d x^2
dC2/dt = C2 - C1(x = 1)
pde's BC : dC1/dx = C1(x=1) - C2
best, Rui
  2 Comments
Bill Greene
Bill Greene on 5 Oct 2013
Hi,
Let me ask a couple of questions first. In the line: pde's BC : dC1/dx = C1(x=1) - C2 is this the BC at x=1? That is, are dC1/dx and C2 both values at x=1? Since, the equation for C1 is second-order, you also need a BC at x=0. What is that?
Bill
Rui
Rui on 2 Nov 2013
Bill, sorry for the slow feedback. The system did not send me the notification.
at X = 0, dC1/dt = 0, no flux BC.
C1 diffuses to the surface at x=1, where some of it goes across interface between C2/C1.

Sign in to comment.

Answers (0)

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!