1. State precisely the mathematical problem describing axially symmetric heat conduction 
    in a hollow cylinder Ω:  Rin ≤ r ≤ Rout,   0 ≤ z ≤ Z, 
    with imposed temperature boundary conditions on all boundaries, 
    starting at a uniform temperature uinit.

2. Set up a 2-dimensional (r,z) mesh for the axially symmetric heat transfer problem above, 
    by subdividing [Rin, Rout] into Mr subintervals and [0,Z] into Mz subintervals.
    Determine (formulas for) nodes:  ri , zj ,  areas of radial faces:  Ari-½,j ,  
    areas of axial faces:  Azi,j-½ ,  and control volumes:  ΔVij.
    Make sure you specify the range of indices for each array.

3. Derive the discrete heat balance for axially symmetric heat conduction in a (hollow) cylinder using the above mesh.  
   [This is the conservation law applied on each control volume, in terms of U and fluxes, 
    where  volumes=..., areas=...., fluxes=....,  as it will be coded up, do NOT plug these into the expression.  
    Do NOT write code, only the formulas needed for update, use half-index notation.]

4. Outline the explicit Finite Volume algorithm (formulas, NOT code)
    for the problem with convective BC at r=Rin, imposed flux at r=Rout,
    and insulated top and bottom (at z=0, z=Z), starting at uniform initial temperature uinit. 
    For simplicity, assume constant properties and mesh uniform in each direction.

5. Find the CFL condition and the stable Δtexpl for the explicit scheme (for internal nodes only).
   [For this you DO need to plug in formulas for fluxes, areas, vol, to see the coefficients]