AbstractOxygen field evaluation is important in modeling and simulation of many important physiological processes of animals, such as angiogenesis. However, numerical simulation of the oxygen field in animal tissue is usually limited by the unusual coupling of different mechanisms, the nonlinearity of the model, and the complex geometry of refined blood vessel networks. In this work, a fast numerical method is designed for the simulation of oxygen supply in tissue with a large-scale complex vessel network. This method employs an implicit finite-difference scheme to compute the oxygen field. By virtue of an oxygen source distribution technique from vessel center lines to mesh points and a corresponding post-processing technique that eliminate the local numerical error induced by source distribution, square mesh with relatively large mesh sizes can be applied while sufficient numerical accuracy is maintained. The new method has computational complexity which is slightly higher than linear with respect to the number of mesh points and has a convergence order which is slightly lower than second order with respect to the mesh size. As an example, the oxygen field of a tissue irrigated by a blood vessel network with more than four thousand blood vessels can be accurately computed within one minute with our new method. The new method will definitely promote further researches based on evaluation of oxygen field, such as modeling of angiogenesis and pathogenesis of many cardiovascular diseases.