Patent ID: 8296112

Claim:
A computer implemented method for simulating a flow of incompressible fluid, using a particle method, the method comprising: representing an incompressible fluid by a plurality of particles, the particles are grouped according to different particle sizes depending on a spatial resolution required at positions in a simulation domain; temporarily updating a velocity and a position coordinate of each of the particles to a first velocity and a first position coordinate by implicitly calculating a variation of the velocity of each of the particles due to a viscosity of the incompressible fluid in each of time steps having a predetermined time interval; and updating the first velocity and the first position coordinate to a second velocity and a second position coordinate of each of the particles at a next time step of the time steps by calculating a velocity correction of the first velocity due to a pressure gradient of the incompressible fluid using the first velocity, wherein the velocity correction is calculated with an incompressibility condition where a divergence of the first velocity of each of the particles is zero, wherein the variation due to the viscosity is calculated by a weight providing a magnitude of interaction between the particles using an interaction radius of each of the particles and a volume of each of the particles, the interaction radius defined as a radius of an interaction area in which the interaction between the particles extends, wherein the weight wij between particles i and j among the particles is zero when a distance A between the particles i and j is greater than an arithmetic average B of interaction radii of the particles i and j, and the weight wij is a product of a divided value of a difference between the distance A and the arithmetic average B by the distance A and a volume ratio of the particle j to the particle i when the distance A is equal to or less than the arithmetic average B.