看看我的command文件对不对？

xiaxiaoxia

我要算的是铝电解槽中的气液两相流，液体装在槽子中。气泡体积分数为10%，在固体壁面上产生（把此壁面当作气体的入口），再进入液体中运动，搅动液体循环，最后气泡由出口（设为Press型patch）逸出。
看看我的command文件对不对？还有那些要注意的问题吗？我用的是CFX4.3
>>CFX4
  >>SET LIMITS
    TOTAL INTEGER WORK SPACE 50000000
    TOTAL CHARACTER WORK SPACE 1000000
    TOTAL REAL WORK SPACE 100000000
    MAXIMUM NUMBER OF BLOCKS 200
    MAXIMUM NUMBER OF PATCHES 100
    MAXIMUM NUMBER OF INTER BLOCK BOUNDARIES 100
  >>OPTIONS
    TURBULENT FLOW
    BUOYANT FLOW
    NUMBER OF PHASES 2
>>MODEL DATA
  >>AMBIENT VARIABLES
    PHASE NAME 'PHASE1'
    VOLUME FRACTION 9.0000E-01
  >>AMBIENT VARIABLES
    PHASE NAME 'PHASE2'
    VOLUME FRACTION 1.0000E-01
  >>DIFFERENCING SCHEME
    U VELOCITY 'HYBRID'
    V VELOCITY 'HYBRID'
    W VELOCITY 'HYBRID'
    PRESSURE 'CENTRAL'
    VOLUME FRACTION 'UPWIND'
    K 'HYBRID'
    EPSILON 'HYBRID'
  >>HYSICAL PROPERTIES
    >>BUOYANCY PARAMETERS
      GRAVITY VECTOR 0.000000E+00 -9.810000E+00 0.000000E+00
      BUOYANCY REFERENCE DENSITY 2.1300E+03
    >>FLUID PARAMETERS
      PHASE NAME 'PHASE1'
      VISCOSITY 3.0000E-03
      DENSITY 2.1300E+03
    >>FLUID PARAMETERS
      PHASE NAME 'PHASE2'
      VISCOSITY 1.8773E-04
      DENSITY 4.0000E-01
    >>MULTIPHASE PARAMETERS
      >>HASE DESCRIPTION
        PHASE NAME 'PHASE1'
        LIQUID
        CONTINUOUS
      >>HASE DESCRIPTION
        PHASE NAME 'PHASE2'
        GAS
        DISPERSE
        MEAN DIAMETER 2.0000E-02
      >>MULTIPHASE MODELS
        >>MOMENTUM
          INTER PHASE TRANSFER
          IPSAC
      >>INTER PHASE TRANSFER MODELS
        >>MOMENTUM
          FIRST PHASE NAME 'PHASE1'
          SECOND PHASE NAME 'PHASE2'
          SURFACE TENSION COEFFICIENT 2.5000E-02
          >>PARTICLE DRAG MODEL
            FLOW REGIME 'DISTORTED'
          >>NON DRAG FORCES
            LIFT FORCE
            LIFT COEFFICIENT 3.0000E-02
            WALL LUBRICATION FORCE
            LUBRICATION COEFFICIENTS -1.000000E-02 5.000000E-02
            TURBULENT DISPERSION FORCE
            DISPERSION COEFFICIENT 1.0000E-01
    >>TURBULENCE PARAMETERS
      >>TURBULENCE MODEL
        PHASE NAME 'PHASE1'
        TURBULENCE MODEL 'K-EPSILON'
        PARTICLE INDUCED TURBULENCE 'SATO'
      >>TURBULENCE MODEL
        PHASE NAME 'PHASE2'
        TURBULENCE MODEL 'K-EPSILON'
        PARTICLE INDUCED TURBULENCE 'SATO'
      >>TURBULENCE CONSTANTS
        CMUPI 6.0000E-01
>>SOLVER DATA
  >>PROGRAM CONTROL
    MAXIMUM NUMBER OF ITERATIONS 2000
    MASS SOURCE TOLERANCE 1.0000E-04
  >>EQUATION SOLVERS
    U VELOCITY 'STONE'
    V VELOCITY 'STONE'
    W VELOCITY 'STONE'
    PRESSURE 'ICCG'
    VOLUME FRACTION 'STONE'
    K 'STONE'
    EPSILON 'STONE'
  >>FALSE TIMESTEPS
    U VELOCITY 1.0000E-03
    V VELOCITY 1.0000E-03
    W VELOCITY 1.0000E-03
    VOLUME FRACTION 1.0000E-03
    K 1.0000E-03
    EPSILON 1.0000E-03
  >>REDUCTION FACTORS
    U VELOCITY 1.0000E-01
    V VELOCITY 1.0000E-01
    W VELOCITY 1.0000E-01
    PRESSURE 1.0000E-02
    VOLUME FRACTION 1.0000E-01
    K 1.0000E-01
    EPSILON 1.0000E-01
  >>SWEEPS INFORMATION
    >>MAXIMUM NUMBER
      U VELOCITY 100
      V VELOCITY 100
      W VELOCITY 100
      PRESSURE 500
      VOLUME FRACTION 100
      K 100
      EPSILON 100
  >>UNDER RELAXATION FACTORS
    U VELOCITY 1.0000E-01
    V VELOCITY 1.0000E-01
    W VELOCITY 1.0000E-01
    PRESSURE 1.0000E-01
    VOLUME FRACTION 1.0000E-01
    K 1.0000E-03
    EPSILON 1.0000E-03
>>MODEL BOUNDARY CONDITIONS
  >>INLET BOUNDARIES
    PHASE NAME 'PHASE2'
    PATCH NAME 'INLET1'
    U VELOCITY 0.0000E+00
    V VELOCITY -2.0000E-03
    W VELOCITY 0.0000E+00
    K 1.0000E-04
    EPSILON 1.0000E-04
>>INLET BOUNDARIES
    PHASE NAME 'PHASE2'
    PATCH NAME 'INLET2'
    U VELOCITY 0.0000E+00
    V VELOCITY -2.0000E-03
    W VELOCITY 0.0000E+00
    K 1.0000E-04
    EPSILON 1.0000E-04
  >>PRESSURE BOUNDARIES
    PHASE NAME 'PHASE2'
    PATCH NAME 'PRESS1'
    PRESSURE 0.0000E+00
  >>PRESSURE BOUNDARIES
    PHASE NAME 'PHASE2'
    PATCH NAME 'PRESS2'
    PRESSURE 0.0000E+00
>>STOP

yanhj

没有给出出口，气泡如何出去？
好象有个脱气边界条件，你看看。

xiaxiaoxia

    我设置的时候是把出口设Press类的Patch，再在这个Patch上设压力边界，用这个文件计算的结果倒是气体可以从Press类的Patch中出去。
    我试过把出口定义为OUTLET类型的Patch，可是在边界条件里，只有“入口、压力、质量流”三种边界条件，没办法设出口压力。（压力边界条件只能为Press型的Patch设置压力）。
     “脱气边界条件”？出口设为OUTLET，再在“MODEL BOUNDARY CONDITIONS”里设置吗？我还没见过， 试一下先。
    谢谢 yanhj！请多指教！！：）