Java 有人能帮我理解如何模拟流体吗？

我正在尝试制作一个程序，模拟加工过程中流体的物理特性。IDE中包含了一个示例：

    /**
     * Fluid 
     * by Glen Murphy. 
     * 
     * Click and drag the mouse to move the simulated fluid.
     * Adjust the "res" variable below to change resolution.
     * Code has not been optimised, and will run fairly slowly.
     */

    int res = 2;
    int penSize = 30;
    int lwidth;
    int lheight;
    int pnum = 30000;
    vsquare[][] v;
    vbuffer[][] vbuf;
    particle[] p = new particle[pnum];
    int pcount = 0;
    int mouseXvel = 0;
    int mouseYvel = 0;

    void setup() 
    {
      size(200, 200);
      noStroke();
      frameRate(30);
      lwidth = width/res;
      lheight = height/res;
      v = new vsquare[lwidth+1][lheight+1];
      vbuf = new vbuffer[lwidth+1][lheight+1];
      for (int i = 0; i < pnum; i++) {
        p[i] = new particle(random(res,width-res),random(res,height-res));
      }
      for (int i = 0; i <= lwidth; i++) {
        for (int u = 0; u <= lheight; u++) {
          v[i][u] = new vsquare(i*res,u*res);
          vbuf[i][u] = new vbuffer(i*res,u*res);
        }
      }
    }

    void draw() 
    {
      background(#666666);

      int axvel = mouseX-pmouseX;
      int ayvel = mouseY-pmouseY;

      mouseXvel = (axvel != mouseXvel) ? axvel : 0;
      mouseYvel = (ayvel != mouseYvel) ? ayvel : 0;

      for (int i = 0; i < lwidth; i++) {
        for (int u = 0; u < lheight; u++) {
          vbuf[i][u].updatebuf(i,u);
          v[i][u].col = 32;
        }
      }
      for (int i = 0; i < pnum-1; i++) {
        p[i].updatepos();
      }
      for (int i = 0; i < lwidth; i++) {
        for (int u = 0; u < lheight; u++) {
          v[i][u].addbuffer(i, u);
          v[i][u].updatevels(mouseXvel, mouseYvel);
          v[i][u].display(i, u);
        }
      }
    }

    class particle {
      float x;
      float y;
      float xvel;
      float yvel;
      int pos;
      particle(float xIn, float yIn) {
        x = xIn;
        y = yIn;
      }

      void updatepos() {
        float col1;
        if (x > 0 && x < width && y > 0 && y < height) {
          int vi = (int)(x/res);
          int vu = (int)(y/res);
          vsquare o = v[vi][vu];    

          float ax = (x%res)/res;
          float ay = (y%res)/res;

          xvel += (1-ax)*v[vi][vu].xvel*0.05;
          yvel += (1-ay)*v[vi][vu].yvel*0.05;

          xvel += ax*v[vi+1][vu].xvel*0.05;
          yvel += ax*v[vi+1][vu].yvel*0.05;

          xvel += ay*v[vi][vu+1].xvel*0.05;
          yvel += ay*v[vi][vu+1].yvel*0.05;

          o.col += 4;

          x += xvel;
          y += yvel;
        }
        else {
          x = random(0,width);
          y = random(0,height);
          xvel = 0;
          yvel = 0;
        }

        xvel *= 0.5;
        yvel *= 0.5;
      }
    }

    class vbuffer {
      int x;
      int y;
      float xvel;
      float yvel;
      float pressurex = 0;
      float pressurey = 0;
      float pressure = 0;

      vbuffer(int xIn,int yIn) {
        x = xIn;
        y = yIn;
        pressurex = 0;
        pressurey = 0;
        }

      void updatebuf(int i, int u) {
        if (i>0 && i<lwidth && u>0 && u<lheight) {
          pressurex = (v[i-1][u-1].xvel*0.5 + v[i-1][u].xvel + v[i-1][u+1].xvel*0.5 - v[i+1][u-1].xvel*0.5 - v[i+1][u].xvel - v[i+1][u+1].xvel*0.5);
          pressurey = (v[i-1][u-1].yvel*0.5 + v[i][u-1].yvel + v[i+1][u-1].yvel*0.5 - v[i-1][u+1].yvel*0.5 - v[i][u+1].yvel - v[i+1][u+1].yvel*0.5);
          pressure = (pressurex + pressurey)*0.25;
          }
        }
      }

    class vsquare {
      int x;
      int y;
      float xvel;
      float yvel;
      float col;

      vsquare(int xIn,int yIn) {
        x = xIn;
        y = yIn;
        }

      void addbuffer(int i, int u) {
        if (i>0 && i<lwidth && u>0 && u<lheight) {
          xvel += (vbuf[i-1][u-1].pressure*0.5
                  +vbuf[i-1][u].pressure
                  +vbuf[i-1][u+1].pressure*0.5
                  -vbuf[i+1][u-1].pressure*0.5
                  -vbuf[i+1][u].pressure
                  -vbuf[i+1][u+1].pressure*0.5
                  )*0.25;
          yvel += (vbuf[i-1][u-1].pressure*0.5
                  +vbuf[i][u-1].pressure
                  +vbuf[i+1][u-1].pressure*0.5
                  -vbuf[i-1][u+1].pressure*0.5
                  -vbuf[i][u+1].pressure
                  -vbuf[i+1][u+1].pressure*0.5
                  )*0.25;
          }
        }

      void updatevels(int mvelX, int mvelY) {
        if (mousePressed) {
          float adj = x - mouseX;
          float opp = y - mouseY;
          float dist = sqrt(opp*opp + adj*adj);
          if (dist < penSize) {
            if (dist < 4) dist = penSize;
            float mod = penSize/dist;
            xvel += mvelX*mod;
            yvel += mvelY*mod;
            }
          }

        xvel *= 0.99;
        yvel *= 0.99;
      }

      void display(int i, int u) {
        float tcol = 0;
        if (col > 255) col = 255;
        if (i>0 && i<lwidth-1 && u>0 && u<lheight-1) {
          tcol = (+ v[i][u+1].col 
                  + v[i+1][u].col 
                  + v[i+1][u+1].col*0.5
                  )*0.4;
          tcol = (int)(tcol+col*0.5);
          }
        else {
          tcol = (int)col;
          }
        fill(tcol, tcol, tcol);
        rect(x,y,res,res);
      }

}

/**
*流质
*格伦·墨菲。
* 
*单击并拖动鼠标以移动模拟流体。
*调整下面的“res”变量以更改分辨率。
*代码尚未优化，将运行得相当慢。
*/
int res=2；
int penSize=30；
国际宽度；
整数lheight；
int pnum=30000；
vsquare[]v；
vbuffer[]]vbuf；
粒子[]p=新粒子[pnum]；
int pcount=0；
int mouseXvel=0；
int mouseYvel=0；
无效设置（）
{
大小（200200）；
仰泳（）；
帧率（30）；
lwidth=宽度/分辨率；
lheight=高度/分辨率；
v=新的v方形[L宽度+1][L高度+1]；
vbuf=新的vbuffer[lwidth+1][lheight+1]；
对于（int i=0；i0&&i0&&u0&&i0&&u255）列=255；
如果（i>0&&i0&&u这篇论文是一个很好的起点，它将向您展示流体模拟背后的数学，并在第三章中描述流体解算器的实现。sourceforge中还有一个可用的工具（您需要使用cvs检查源）.
Hi@Miles.Java可能不是你想要的这种物理密集型应用程序的方式。此外，你可能在错误的地方问了这个问题。你可能想先学习流体物理。尝试在physics StackExchange网站上提问。（类似于StackOverflow，但适用于物理）。在开发之前，您应该首先了解您的域。希望这对您有所帮助！