1. A coronary circulation simulation method executed by a first processor and a plurality of second processors, each of the first and the second processors being coupled to a storage device, the coronary circulation simulation method comprising: storing a geometric model representing a shape of a heart as a collection of elements formed from a plurality of nodes and a plurality of connections among the nodes in the storage device; receiving a first vessel network model representing a network of first vessels whose diameters are larger than or equal to a specified threshold, wherein the first vessel network model is a macro model; storing the first vessel network model in the storage device; receiving a plurality of second vessel network models each representing a network of second vessels whose diameters are smaller than the specified threshold, the plurality of second vessel network models being independent of each other, wherein one or more of the second vessel network models are connected to each of the nodes, wherein the second vessel network models are micro models; storing the plurality of second vessel network models in the storage device; first analyzing, by the first processor, first hemodynamics in the first vessels belonging to the first vessel network model, based on the geometric model of the heart; second analyzing, by the plurality of second processors in parallel without interactions, second hemodynamics in the second vessels belonging to the second vessel network models connected to the plurality of nodes, by using blood pressure data in the first vessels at the nodes, the blood pressure data being included in result data of the first analyzing which indicates the first hemodynamics in the first vessels at each of the nodes; and outputting the first hemodynamics in the first vessels and the second hemodynamics in the second vessels, wherein: the first analyzing and the second analyzing include analyzing the first and the second hemodynamics in the first and second vessels reflecting motion of the heart by using a Newton-Raphson method for nonlinear analysis, and are performed alternately; the second analyzing includes calculating tentative blood pressure increments in the second vessel network models by using an equation from which first blood pressure increments in the first vessel network model have been eliminated; the first analyzing includes calculating the first blood pressure increments in the first vessel network model by using the tentative blood pressure increments in the second vessel network models; and the second analyzing includes calculating second blood pressure increments in the second vessel network models by using the first blood pressure increments in the first vessel network model.