Spatial source for diffusion equation
1 次查看(过去 30 天)
显示 更早的评论
Hi all, I would like to solve a diffusion equation -d^2n/dx^2 = S(x) in 1D between -L < x < +L, with a spatial source term, for instance "1" between -0.1L < x < 0.1L and 0 elsewhere. I'm a bit confused by the structure of pdepe. And I can't find an example with a source term.
Can anybody help me ?
Thanks a lot
Stéphane
0 个评论
回答(1 个)
Torsten
2016-6-14
Just take pencil and paper and integrate n''(x)=S(x) twice with respect to x.
If you don't succeed, use bvp4c instead of pdepe since your equation is stationary.
Best wishes
Torsten.
2 个评论
Torsten
2016-6-14
I don't know if the equation you have in mind is more complex than the one you wrote above. But integrating a piecewise constant function (S(x)) twice is a task I usually accompish without computer support.
Best wishes
Torsten.
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!