This paper proposes a quantitative model of the reaction-diffusion type to examine thedistribution of interferon-α (IFNα) in a lymph node(LN). The numerical treatment of the model is based on using an original unstructured meshgeneration software Ani3D and nonlinear finite volume method for diffusion equations. Thestudy results in suggestion that due to the variations in hydraulic conductivity ofvarious zones of the secondary lymphoid organs the spatial stationary distribution ofIFNα is essentially heterogeneous across the organs. Highly protecteddomains such as sinuses, conduits, co-exist with the regions in which where the stationaryconcentration of IFNα is lower by about 100-fold. This is the first studywhere the spatial distribution of soluble immune factors in secondary lymphoid organs ismodelled for a realistic three-dimensional geometry.