A three-dimensional fluid simulation method is disclosed wherein: a first mesh of a flow domain non-uniformly split with respect to each of three degrees of freedom is defined; a second mesh uniformly sprit with respect to only one of the three degrees of freedom but non-uniformly with respect to other two degrees is defined; an object model is set in the first mesh and a motion equation is formed and calculated to obtain fluid velocity; based on the fluid velocity, flow imbalance is computed for each cell; based on the flow imbalance, fluid pressure correction equation is formed; the flow imbalance is mapped onto the second mesh and the fluid pressure correction is computed; the fluid pressure correction is mapped onto the first mesh; and until the flow imbalance and motion equations are converged the computation is repeated.