Exchange Problem with PDEPE solver

3 次查看（过去 30 天）
Thomas
Thomas 2015-1-14
评论： Thomas 2015-1-15
I want to use the PDEPE (or whatever else) solver to describe an exchange model between two different layers. Within the layers I want diffusion to occurr.
I want to start with a "hot" sphere surrounded by a cold spherical shell. Between those two reservoirs I want an exchange of the form d/dt U(leftborder) = -k*(U(leftborder)-U(rightborder). Within the layers I want a propagation (of heat for example) via diffusion.
My problems are:
(1) how to i solve a problem of two partial differential equations that are only coupled at their border in matlab. As far as i understand i can only solve coupled equations that have the same domain, i.e. the same xmesh.
(2) how can I implement a time dependent boundary condition like the one stated above.
Thanks for any help!
Best regards.

请先登录，再回答此问题。

采纳的回答

Torsten
Torsten 2015-1-14
There is no ready-to-use solver in MATLAB for this kind of problem.
You will have to discretize the PDEs in space and solve the resulting system of ODEs for the temperatures in the grid points by a solver for ordinary differential equations (e.g.ODE15s).
What second transmission condition should apply at the interface between the two layers ? Continuity of heat flux ?
Best wishes
Torsten.
  2 个评论
Thomas
Thomas 2015-1-14
Hello Torsten and thank you for your answer! I do not fully understand what you mean by second transmission condition. The model should contain heat transfer towards the border by diffusion and then (kind of a rude discontinuity in the diffusion) an exchange at the border. In the second layer this should propagate by diffusion again. So in principle is a diffuse, hop over barrier (and also backwards hop) then keep diffusing type of problem.
Thank you very much, your answer already helped me a lot!
Torsten
Torsten 2015-1-14
You have two unknowns at the interface: U(leftborder) and U(rightborder). So you need two equations to solve for them.
Usually, if e.g. material properties change at the interface, the transmission conditions are continuity of temperature (Tleft=Tright) and continuity of heat flux (Lambda_left*dT_left/dx=Lambda_right*dT_right/dx).
In your case, only one condition is stated.
Best wishes
Torsten.

请先登录，再进行评论。

更多回答（1 个）

Bill Greene
Bill Greene 2015-1-14
Torsten said:
"Usually, if e.g. material properties change at the interface, the transmission conditions are continuity of temperature (Tleft=Tright) and continuity of heat flux (Lambda_left*dT_left/dx=Lambda_right*dT_right/dx)."
If, in fact, these are the conditions at the material interface, solution using the pdepe function is straightforward. You have a single xmesh that spans both layers. In your function that calculates the pde coefficients, you use the value of the x-variable that is passed in to calculate the flux, i.e, if x < x_interface, f= k1*dudx, otherwise f = k2*dudx.
Bill
  4 个评论
显示 2更早的评论
Torsten
Torsten 2015-1-15
Well, I repeat that PDEPE is not suited for your problem, but maybe its worth for gaining experience with your problem.
Best wishes
Torsten.
Thomas
Thomas 2015-1-15
Hm, okay it looks like i will have to do a numeric solution then! Thanks for the help!

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Structural Mechanics 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by