Discretization of FVM 1D particle model, heat transport

1 view (last 30 days)
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?

Answers (1)

Torsten
Torsten on 6 Jan 2015
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.

Categories

Find more on Numerical Integration and Differential Equations in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!