Discretization of FVM 1D particle model, heat transport

6 次查看(过去 30 天)
I'm trying to solve a boundary value problem using an ode solver.
It solves heat transport by conduction in 1 dimension of a spherical particle. The boundary condition in the center a symmetry condition. The boundary condition on the particle surface is given by a heat flux by both convection and conduction (algebraic eqs).
I am having some trouble implementing the boundary conditions. Using FVM all information is stored in the center of a cell, thus giving problems when I need the information on the node. Introducing ghost cells would give similar problems?
Any ideas?

回答(1 个)

Torsten
Torsten 2015-1-6
I don't understand your problem.
At the center, u_(-1/2)=u_(1/2).
At the outer boundary, (u_(N+1/2)-u_(N-1/2))/deltar is approximately du/dr at r=R.
Best wishes
Torsten.

类别

Help CenterFile Exchange 中查找有关 Numerical Integration and Differential Equations 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by