创建一个流体动力学效果着色器，包含流动、湍流和涡旋效果，支持流速、粘度调节，提供三种流体类型可选。

uniform float uTime;
uniform float uFlowSpeed;
uniform float uViscosity;
uniform float uWaveFrequency;

void main() {
    vUv = uv;
    vPosition = position;
    vNormal = normal;

    // 流体流动效果
    float flow = sin(position.x * 3.0 + uTime * uFlowSpeed) *
                cos(position.y * 4.0 + uTime * uFlowSpeed * 0.8) *
                sin(position.z * 3.0 + uTime * uFlowSpeed * 0.6);

    // 粘度效果
    float viscosity = sin(position.x * 5.0 + uTime * uViscosity * 2.0) *
                     cos(position.z * 5.0 + uTime * uViscosity * 2.0);

    // 波动频率
    float wave = sin(position.x * uWaveFrequency * 6.0) *
                cos(position.y * uWaveFrequency * 8.0);

    vec3 newPosition = position + normal * (flow + viscosity * 0.5 + wave * 0.3) * 0.15;

    gl_Position = projectionMatrix * modelViewMatrix * vec4(newPosition, 1.0);
}

float fluidFlow(vec2 uv, float speed) {
    float flowField = 0.0;

    // 多层流体流动
    flowField += sin(uv.x * 4.0 + uTime * 0.3) *
                cos(uv.y * 6.0 + uTime * 0.5);
    flowField += sin(uv.x * 6.0 + uTime * 0.7) *
                cos(uv.y * 8.0 + uTime * 0.9) * 0.7;
    flowField += sin(uv.x * 8.0 + uTime * 1.1) *
                cos(uv.y * 10.0 + uTime * 1.3) * 0.4;

    return flowField * speed;
}

float turbulence(vec2 uv, float intensity) {
    float turbField = 0.0;
    for(int i = 0; i < 20; i++) {
        vec2 turbPos = vec2(
            noise(vec2(float(i) * 0.123, 0.456)),
            noise(vec2(float(i) * 0.456, 0.789))
        );
        float dist = distance(uv, turbPos);
        float turbSize = noise(vec2(float(i) * 0.789, 0.123)) * 0.04 + 0.02;
        float turbLife = mod(uTime * 1.5 + float(i) * 8.0, 6.28);
        float turbIntensity = sin(turbLife) * 0.7 + 0.3;
        turbField += exp(-dist * 45.0) * turbIntensity * intensity;
    }
    return turbField;
}

float vortices(vec2 uv, float intensity) {
    if(!uShowVortices) return 0.0;

    float vortexField = 0.0;

    // 涡旋流动
    float angle = atan(uv.y - 0.5, uv.x - 0.5);
    float radius = distance(uv, vec2(0.5, 0.5));
    vortexField += sin(angle * 8.0 + uTime * 2.0) *
                  cos(radius * 12.0 - uTime * 1.5) * intensity;

    return vortexField;
}

uniforms: {
    uTime: { value: 0.0 },
    uFlowSpeed: { value: 0.8 },
    uViscosity: { value: 0.6 },
    uTurbulence: { value: 0.7 },
    uGlowIntensity: { value: 1.2 },
    uWaveFrequency: { value: 0.8 },
    uFluidType: { value: 0 },
    uShowVortices: { value: true }
}