Skip to content
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]
皮肤
  • Light
  • Cerulean
  • Cosmo
  • Flatly
  • Journal
  • Litera
  • Lumen
  • Lux
  • Materia
  • Minty
  • Morph
  • Pulse
  • Sandstone
  • Simplex
  • Sketchy
  • Spacelab
  • United
  • Yeti
  • Zephyr
  • Dark
  • Cyborg
  • Darkly
  • Quartz
  • Slate
  • Solar
  • Superhero
  • Vapor

  • 默认(不使用皮肤)
  • 不使用皮肤
折叠
CFD中文网

CFD中文网

  1. CFD中文网
  2. OpenFOAM
  3. 请教关于进口条件的设置的问题

请教关于进口条件的设置的问题

已定时 已固定 已锁定 已移动 OpenFOAM
32 帖子 6 发布者 37.0k 浏览
  • 从旧到新
  • 从新到旧
  • 最多赞同
回复
  • 在新帖中回复
登录后回复
此主题已被删除。只有拥有主题管理权限的用户可以查看。
  • R 离线
    R 离线
    random_ran 大神
    写于 最后由 编辑
    #17

    我注意到你的log 文件里面从第二个时间步长开始, specific dissipation rate (omega)就已经发散了(10^9)。

    omega 文件里

    FoamFile
    {
        version     2.0;
        format      ascii;
        class       volScalarField;
        object      omega;
    }
    
    你的 internalField uniform $omegaInlet = 630 这个神奇的数字是怎么算的?
    从你的log file 看出,你是从2000 s 开始算的,贴出2000s的初始条件或许能帮助解决问题?
    
    对于你的的进,出口边界条件不是特别理解。用进口速度边界条件,出口用压力边界条件不能满足你的需求吗?气动问题大家通常都是用这个方法。

    Yours in CFD,

    Ran

    小 1 条回复 最后回复
  • 小 离线
    小 离线
    小葱小虫
    在 中回复了 random_ran 最后由 编辑
    #18

    @random_ran 我开始没有开湍流模型,想算个初场,2000部之前都是没有湍流模型的,omega用这个公式算的
    omega = k^0.5/(C^0.25 * l) c=0.09
    我的是进口条件是进口已知总压,和气流方向出口知道静压,算出来速度大约是0.67Ma

    1 条回复 最后回复
  • 小 离线
    小 离线
    小葱小虫
    写于 最后由 CFD中文网 编辑
    #19

    这是最新算的几部,麻烦老师们给看一下

    | =========                 |                                                 |
    | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
    |  \\    /   O peration     | Version:  3.0.1                                 |
    |   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
    |    \\/     M anipulation  |                                                 |
    \*---------------------------------------------------------------------------*/
    Build  : 3.0.1-119cac7e8750
    Exec   : simpleFoam
    Date   : Mar 25 2017
    Time   : 02:33:53
    Host   : "ubuntu"
    PID    : 44090
    Case   : /home/cc/OpenFOAM/cc-3.0.1/run/tutorials/incompressible/simpleFoam/nasa65
    nProcs : 1
    sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
    fileModificationChecking : Monitoring run-time modified files using timeStampMaster
    allowSystemOperations : Allowing user-supplied system call operations
    
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    Create time
    
    Create mesh for time = 0
    
    
    SIMPLE: convergence criteria
        field U	 tolerance 1e-05
        field p	 tolerance 1e-05
        field "(k|epsilon|omega)"	 tolerance 1e-05
    
    Reading field p
    
    AMI: Creating addressing and weights between 15340 source faces and 15340 target faces
    AMI: Patch source sum(weights) min/max/average = 0.999911, 1.00059, 1.00001
    AMI: Patch target sum(weights) min/max/average = 0.999839, 1.00646, 1.00006
    Reading field U
    
    Reading/calculating face flux field phi
    
    Selecting incompressible transport model Newtonian
    Selecting turbulence model type RAS
    Selecting RAS turbulence model kOmegaSST
    Selecting patchDistMethod meshWave
    kOmegaSSTCoeffs
    {
        alphaK1         0.85;
        alphaK2         1;
        alphaOmega1     0.5;
        alphaOmega2     0.856;
        gamma1          0.555556;
        gamma2          0.44;
        beta1           0.075;
        beta2           0.0828;
        betaStar        0.09;
        a1              0.31;
        b1              1;
        c1              10;
        F3              false;
    }
    
    No MRF models present
    
    No finite volume options present
    
    
    Starting time loop
    
    Time = 1
    
    smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 0.0813711, No Iterations 6
    smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.0714855, No Iterations 5
    smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 0.0771402, No Iterations 4
    GAMG:  Solving for p, Initial residual = 1, Final residual = 0.0876464, No Iterations 5
    time step continuity errors : sum local = 10.8537, global = 4.01707, cumulative = 4.01707
    smoothSolver:  Solving for omega, Initial residual = 0.259035, Final residual = 0.0221475, No Iterations 3
    smoothSolver:  Solving for k, Initial residual = 1, Final residual = 0.0890898, No Iterations 5
    ExecutionTime = 34.31 s  ClockTime = 36 s
    
    Time = 2
    
    smoothSolver:  Solving for Ux, Initial residual = 0.886786, Final residual = 0.058124, No Iterations 6
    smoothSolver:  Solving for Uy, Initial residual = 0.979345, Final residual = 0.0631154, No Iterations 6
    smoothSolver:  Solving for Uz, Initial residual = 0.810242, Final residual = 0.0670691, No Iterations 6
    GAMG:  Solving for p, Initial residual = 0.986904, Final residual = 0.0565546, No Iterations 4
    time step continuity errors : sum local = 1345.89, global = -10.0278, cumulative = -6.0107
    smoothSolver:  Solving for omega, Initial residual = 0.0850778, Final residual = 0.00220061, No Iterations 1
    smoothSolver:  Solving for k, Initial residual = 0.814193, Final residual = 0.0238861, No Iterations 2
    ExecutionTime = 45.36 s  ClockTime = 47 s
    
    Time = 3
    
    smoothSolver:  Solving for Ux, Initial residual = 0.0193303, Final residual = 0.00152468, No Iterations 6
    smoothSolver:  Solving for Uy, Initial residual = 0.011814, Final residual = 0.00102405, No Iterations 6
    smoothSolver:  Solving for Uz, Initial residual = 0.422504, Final residual = 0.0419746, No Iterations 6
    GAMG:  Solving for p, Initial residual = 0.791662, Final residual = 0.0642592, No Iterations 4
    time step continuity errors : sum local = 957.279, global = -143.421, cumulative = -149.432
    smoothSolver:  Solving for omega, Initial residual = 0.00274129, Final residual = 4.02132e-05, No Iterations 2
    smoothSolver:  Solving for k, Initial residual = 0.169936, Final residual = 0.0157838, No Iterations 1
    ExecutionTime = 56.05 s  ClockTime = 58 s
    
    Time = 4
    
    smoothSolver:  Solving for Ux, Initial residual = 0.0983965, Final residual = 0.00679746, No Iterations 6
    smoothSolver:  Solving for Uy, Initial residual = 0.0657851, Final residual = 0.00457164, No Iterations 6
    smoothSolver:  Solving for Uz, Initial residual = 0.206983, Final residual = 0.0189299, No Iterations 6
    GAMG:  Solving for p, Initial residual = 0.980449, Final residual = 0.0634521, No Iterations 4
    time step continuity errors : sum local = 1261.26, global = -92.2404, cumulative = -241.672
    smoothSolver:  Solving for omega, Initial residual = 0.000609973, Final residual = 3.54731e-05, No Iterations 1
    smoothSolver:  Solving for k, Initial residual = 0.10451, Final residual = 0.0098634, No Iterations 1
    ExecutionTime = 66.79 s  ClockTime = 69 s
    
    Time = 5
    
    smoothSolver:  Solving for Ux, Initial residual = 0.0261376, Final residual = 0.00214716, No Iterations 6
    smoothSolver:  Solving for Uy, Initial residual = 0.0187723, Final residual = 0.00145964, No Iterations 6
    smoothSolver:  Solving for Uz, Initial residual = 0.0626224, Final residual = 0.00595206, No Iterations 6
    GAMG:  Solving for p, Initial residual = 0.921146, Final residual = 0.0735434, No Iterations 4
    time step continuity errors : sum local = 1305.03, global = -174.17, cumulative = -415.842
    smoothSolver:  Solving for omega, Initial residual = 0.000273909, Final residual = 9.03935e-06, No Iterations 1
    smoothSolver:  Solving for k, Initial residual = 0.0349887, Final residual = 0.00218001, No Iterations 3
    ExecutionTime = 78.4 s  ClockTime = 81 s
    
    Time = 6
    
    smoothSolver:  Solving for Ux, Initial residual = 0.0566235, Final residual = 0.00458913, No Iterations 6
    smoothSolver:  Solving for Uy, Initial residual = 0.0400099, Final residual = 0.00302616, No Iterations 6
    smoothSolver:  Solving for Uz, Initial residual = 0.0384566, Final residual = 0.00324406, No Iterations 6
    GAMG:  Solving for p, Initial residual = 0.916166, Final residual = 0.0603851, No Iterations 4
    time step continuity errors : sum local = 668.506, global = -94.6009, cumulative = -510.443
    smoothSolver:  Solving for omega, Initial residual = 0.000292328, Final residual = 1.62378e-05, No Iterations 1
    smoothSolver:  Solving for k, Initial residual = 0.046623, Final residual = 0.00262596, No Iterations 2
    ExecutionTime = 89.53 s  ClockTime = 92 s
    
    Time = 7
    
    smoothSolver:  Solving for Ux, Initial residual = 0.233958, Final residual = 0.0176323, No Iterations 6
    smoothSolver:  Solving for Uy, Initial residual = 0.128911, Final residual = 0.00936957, No Iterations 6
    smoothSolver:  Solving for Uz, Initial residual = 0.0442786, Final residual = 0.00376373, No Iterations 6
    GAMG:  Solving for p, Initial residual = 0.952315, Final residual = 0.0728294, No Iterations 4
    time step continuity errors : sum local = 1263.9, global = -195.294, cumulative = -705.737
    smoothSolver:  Solving for omega, Initial residual = 0.000206228, Final residual = 1.22213e-05, No Iterations 1
    smoothSolver:  Solving for k, Initial residual = 0.0167436, Final residual = 0.00109875, No Iterations 4
    ExecutionTime = 100.52 s  ClockTime = 103 s
    
    Time = 8
    
    smoothSolver:  Solving for Ux, Initial residual = 0.064942, Final residual = 0.00538914, No Iterations 6
    smoothSolver:  Solving for Uy, Initial residual = 0.0316677, Final residual = 0.00284742, No Iterations 6
    smoothSolver:  Solving for Uz, Initial residual = 0.0769102, Final residual = 0.00692843, No Iterations 6
    GAMG:  Solving for p, Initial residual = 0.791177, Final residual = 0.0625584, No Iterations 6
    time step continuity errors : sum local = 645.649, global = -319.458, cumulative = -1025.19
    smoothSolver:  Solving for omega, Initial residual = 0.000258477, Final residual = 1.59141e-05, No Iterations 1
    smoothSolver:  Solving for k, Initial residual = 0.0181236, Final residual = 0.00171421, No Iterations 3
    ExecutionTime = 113.31 s  ClockTime = 116 s
    
    Time = 9
    
    smoothSolver:  Solving for Ux, Initial residual = 0.219725, Final residual = 0.0183762, No Iterations 6
    smoothSolver:  Solving for Uy, Initial residual = 0.158121, Final residual = 0.0111089, No Iterations 7
    smoothSolver:  Solving for Uz, Initial residual = 0.303528, Final residual = 0.0257045, No Iterations 6
    GAMG:  Solving for p, Initial residual = 0.94797, Final residual = 0.064793, No Iterations 4
    time step continuity errors : sum local = 1113.34, global = -309.45, cumulative = -1334.64
    smoothSolver:  Solving for omega, Initial residual = 0.000122954, Final residual = 1.00576e-05, No Iterations 1
    smoothSolver:  Solving for k, Initial residual = 0.0197785, Final residual = 0.00176744, No Iterations 3
    ExecutionTime = 126.36 s  ClockTime = 129 s
    
    Time = 10
    
    smoothSolver:  Solving for Ux, Initial residual = 0.0821976, Final residual = 0.00721542, No Iterations 6
    smoothSolver:  Solving for Uy, Initial residual = 0.03459, Final residual = 0.00247223, No Iterations 7
    smoothSolver:  Solving for Uz, Initial residual = 0.1429, Final residual = 0.0115741, No Iterations 6
    GAMG:  Solving for p, Initial residual = 0.861997, Final residual = 0.084008, No Iterations 5
    time step continuity errors : sum local = 1083.47, global = 389.776, cumulative = -944.868
    smoothSolver:  Solving for omega, Initial residual = 0.000136441, Final residual = 2.43071e-06, No Iterations 2
    smoothSolver:  Solving for k, Initial residual = 0.0235868, Final residual = 0.00172782, No Iterations 3
    ExecutionTime = 138.69 s  ClockTime = 141 s
    
    Time = 11
    
    smoothSolver:  Solving for Ux, Initial residual = 0.163168, Final residual = 0.0136899, No Iterations 6
    smoothSolver:  Solving for Uy, Initial residual = 0.0793402, Final residual = 0.00752288, No Iterations 5
    smoothSolver:  Solving for Uz, Initial residual = 0.386931, Final residual = 0.0294783, No Iterations 6
    GAMG:  Solving for p, Initial residual = 0.870129, Final residual = 0.0711586, No Iterations 6
    time step continuity errors : sum local = 919.618, global = 516.016, cumulative = -428.851
    smoothSolver:  Solving for omega, Initial residual = 0.000106145, Final residual = 1.01946e-05, No Iterations 1
    smoothSolver:  Solving for k, Initial residual = 0.0196961, Final residual = 0.0018465, No Iterations 3
    ExecutionTime = 151.28 s  ClockTime = 154 s
    
    Time = 12
    
    smoothSolver:  Solving for Ux, Initial residual = 0.274151, Final residual = 0.0245354, No Iterations 6
    smoothSolver:  Solving for Uy, Initial residual = 0.186447, Final residual = 0.0165338, No Iterations 5
    smoothSolver:  Solving for Uz, Initial residual = 0.156144, Final residual = 0.0136401, No Iterations 6
    GAMG:  Solving for p, Initial residual = 0.887478, Final residual = 0.070432, No Iterations 6
    time step continuity errors : sum local = 1119.32, global = 661.637, cumulative = 232.786
    smoothSolver:  Solving for omega, Initial residual = 0.000100293, Final residual = 2.24457e-06, No Iterations 2
    smoothSolver:  Solving for k, Initial residual = 0.0248731, Final residual = 0.00212404, No Iterations 3
    ExecutionTime = 164.01 s  ClockTime = 167 s
    
    Time = 13
    
    smoothSolver:  Solving for Ux, Initial residual = 0.205525, Final residual = 0.0187076, No Iterations 6
    smoothSolver:  Solving for Uy, Initial residual = 0.0923904, Final residual = 0.00894338, No Iterations 5
    smoothSolver:  Solving for Uz, Initial residual = 0.196898, Final residual = 0.0175523, No Iterations 6
    GAMG:  Solving for p, Initial residual = 0.817651, Final residual = 0.0790382, No Iterations 8
    time step continuity errors : sum local = 1128.57, global = 728.044, cumulative = 960.83
    smoothSolver:  Solving for omega, Initial residual = 0.00010952, Final residual = 2.80076e-06, No Iterations 2
    smoothSolver:  Solving for k, Initial residual = 0.0218195, Final residual = 0.00158256, No Iterations 4
    ExecutionTime = 178.57 s  ClockTime = 181 s
    
    Time = 14
    
    smoothSolver:  Solving for Ux, Initial residual = 0.209248, Final residual = 0.0174402, No Iterations 6
    smoothSolver:  Solving for Uy, Initial residual = 0.234617, Final residual = 0.0220479, No Iterations 6
    smoothSolver:  Solving for Uz, Initial residual = 0.330718, Final residual = 0.0281098, No Iterations 6
    GAMG:  Solving for p, Initial residual = 0.871195, Final residual = 0.0759581, No Iterations 10
    time step continuity errors : sum local = 1172.79, global = 858.035, cumulative = 1818.86
    smoothSolver:  Solving for omega, Initial residual = 0.000111802, Final residual = 3.26755e-06, No Iterations 2
    smoothSolver:  Solving for k, Initial residual = 0.0256925, Final residual = 0.00171392, No Iterations 4
    ExecutionTime = 193.33 s  ClockTime = 196 s
    
    Time = 15
    
    smoothSolver:  Solving for Ux, Initial residual = 0.173376, Final residual = 0.0142151, No Iterations 6
    smoothSolver:  Solving for Uy, Initial residual = 0.141625, Final residual = 0.0112524, No Iterations 7
    smoothSolver:  Solving for Uz, Initial residual = 0.203415, Final residual = 0.0168199, No Iterations 6
    GAMG:  Solving for p, Initial residual = 0.850352, Final residual = 0.0820699, No Iterations 8
    time step continuity errors : sum local = 1261.82, global = 852.091, cumulative = 2670.96
    smoothSolver:  Solving for omega, Initial residual = 0.000119946, Final residual = 3.04546e-06, No Iterations 2
    smoothSolver:  Solving for k, Initial residual = 0.0523867, Final residual = 0.00451596, No Iterations 2
    ExecutionTime = 206.27 s  ClockTime = 209 s
    
    Time = 16
    
    smoothSolver:  Solving for Ux, Initial residual = 0.189402, Final residual = 0.0147132, No Iterations 6
    smoothSolver:  Solving for Uy, Initial residual = 0.172538, Final residual = 0.0148313, No Iterations 6
    smoothSolver:  Solving for Uz, Initial residual = 0.216101, Final residual = 0.0175322, No Iterations 6
    GAMG:  Solving for p, Initial residual = 0.831111, Final residual = 0.0800892, No Iterations 13
    time step continuity errors : sum local = 1124.03, global = -888.775, cumulative = 1782.18
    smoothSolver:  Solving for omega, Initial residual = 0.000132821, Final residual = 2.71341e-06, No Iterations 2
    smoothSolver:  Solving for k, Initial residual = 0.0429148, Final residual = 0.00279213, No Iterations 3
    ExecutionTime = 220.4 s  ClockTime = 223 s
    
    Time = 17
    
    smoothSolver:  Solving for Ux, Initial residual = 0.240906, Final residual = 0.0185245, No Iterations 6
    smoothSolver:  Solving for Uy, Initial residual = 0.481738, Final residual = 0.0406032, No Iterations 5
    smoothSolver:  Solving for Uz, Initial residual = 0.306464, Final residual = 0.0227411, No Iterations 6
    GAMG:  Solving for p, Initial residual = 0.829137, Final residual = 0.0797512, No Iterations 5
    time step continuity errors : sum local = 1176.95, global = 495.085, cumulative = 2277.27
    smoothSolver:  Solving for omega, Initial residual = 0.000148683, Final residual = 1.27859e-05, No Iterations 1
    smoothSolver:  Solving for k, Initial residual = 0.0542324, Final residual = 0.0033722, No Iterations 3
    ExecutionTime = 231.68 s  ClockTime = 235 s
    
    Time = 18
    
    smoothSolver:  Solving for Ux, Initial residual = 0.206106, Final residual = 0.0161844, No Iterations 6
    smoothSolver:  Solving for Uy, Initial residual = 0.255402, Final residual = 0.0223374, No Iterations 5
    smoothSolver:  Solving for Uz, Initial residual = 0.341787, Final residual = 0.0275786, No Iterations 6
    GAMG:  Solving for p, Initial residual = 0.827099, Final residual = 0.082525, No Iterations 7
    time step continuity errors : sum local = 1134.39, global = 819.278, cumulative = 3096.54
    smoothSolver:  Solving for omega, Initial residual = 9.46476e-05, Final residual = 2.25188e-06, No Iterations 2
    smoothSolver:  Solving for k, Initial residual = 0.0505738, Final residual = 0.00348548, No Iterations 3
    ExecutionTime = 242.31 s  ClockTime = 246 s
    
    Time = 19
    
    smoothSolver:  Solving for Ux, Initial residual = 0.199931, Final residual = 0.0167285, No Iterations 6
    smoothSolver:  Solving for Uy, Initial residual = 0.361021, Final residual = 0.0275666, No Iterations 6
    smoothSolver:  Solving for Uz, Initial residual = 0.234612, Final residual = 0.0186392, No Iterations 6
    GAMG:  Solving for p, Initial residual = 0.829252, Final residual = 0.0781392, No Iterations 18
    time step continuity errors : sum local = 1182.87, global = 1048.32, cumulative = 4144.86
    smoothSolver:  Solving for omega, Initial residual = 0.000103443, Final residual = 9.28303e-06, No Iterations 1
    smoothSolver:  Solving for k, Initial residual = 0.160832, Final residual = 0.0122526, No Iterations 1
    ExecutionTime = 255.21 s  ClockTime = 259 s
    
    Time = 20
    
    smoothSolver:  Solving for Ux, Initial residual = 0.203207, Final residual = 0.0177404, No Iterations 6
    smoothSolver:  Solving for Uy, Initial residual = 0.328249, Final residual = 0.0269088, No Iterations 6
    smoothSolver:  Solving for Uz, Initial residual = 0.224519, Final residual = 0.0164899, No Iterations 6
    GAMG:  Solving for p, Initial residual = 0.838448, Final residual = 0.0774968, No Iterations 16
    time step continuity errors : sum local = 1282.19, global = 1084.48, cumulative = 5229.34
    smoothSolver:  Solving for omega, Initial residual = 0.000137159, Final residual = 1.18779e-05, No Iterations 1
    smoothSolver:  Solving for k, Initial residual = 0.105212, Final residual = 0.00733032, No Iterations 2
    ExecutionTime = 268.59 s  ClockTime = 273 s
    
    Time = 21
    
    smoothSolver:  Solving for Ux, Initial residual = 0.195891, Final residual = 0.0166656, No Iterations 6
    smoothSolver:  Solving for Uy, Initial residual = 0.295498, Final residual = 0.0264212, No Iterations 6
    smoothSolver:  Solving for Uz, Initial residual = 0.256556, Final residual = 0.019235, No Iterations 6
    GAMG:  Solving for p, Initial residual = 0.820159, Final residual = 0.0809549, No Iterations 27
    time step continuity errors : sum local = 1210.13, global = -1110.62, cumulative = 4118.72
    smoothSolver:  Solving for omega, Initial residual = 0.000209472, Final residual = 1.33081e-05, No Iterations 1
    smoothSolver:  Solving for k, Initial residual = 0.0963909, Final residual = 0.00779254, No Iterations 2
    ExecutionTime = 283.95 s  ClockTime = 288 s
    
    Time = 22
    
    smoothSolver:  Solving for Ux, Initial residual = 0.210533, Final residual = 0.0180219, No Iterations 5
    smoothSolver:  Solving for Uy, Initial residual = 0.451871, Final residual = 0.0368704, No Iterations 6
    smoothSolver:  Solving for Uz, Initial residual = 0.267514, Final residual = 0.0242343, No Iterations 5
    GAMG:  Solving for p, Initial residual = 0.828406, Final residual = 0.0758993, No Iterations 30
    time step continuity errors : sum local = 1183.8, global = 1080.35, cumulative = 5199.06
    smoothSolver:  Solving for omega, Initial residual = 0.00040392, Final residual = 1.58342e-05, No Iterations 1
    smoothSolver:  Solving for k, Initial residual = 0.116199, Final residual = 0.00839545, No Iterations 2
    ExecutionTime = 299.9 s  ClockTime = 304 s
    
    Time = 23
    
    smoothSolver:  Solving for Ux, Initial residual = 0.223249, Final residual = 0.0207695, No Iterations 5
    smoothSolver:  Solving for Uy, Initial residual = 0.479522, Final residual = 0.0438371, No Iterations 5
    smoothSolver:  Solving for Uz, Initial residual = 0.287509, Final residual = 0.0233311, No Iterations 5
    GAMG:  Solving for p, Initial residual = 0.819938, Final residual = 0.0816971, No Iterations 21
    time step continuity errors : sum local = 1364, global = -1206.8, cumulative = 3992.26
    smoothSolver:  Solving for omega, Initial residual = 0.000308337, Final residual = 2.31198e-05, No Iterations 1
    smoothSolver:  Solving for k, Initial residual = 0.183739, Final residual = 0.0103415, No Iterations 2
    ExecutionTime = 314.31 s  ClockTime = 319 s
    1 条回复 最后回复
  • C 离线
    C 离线
    CFD中文网
    写于 最后由 编辑
    #20

    质量严重不守恒,为何不把网格做到10万以下,轻量级跑跑试试?

    CFD中国标准用户测试帐号
    目前由徐笑笑登录

    小 1 条回复 最后回复
  • 小 离线
    小 离线
    小葱小虫
    在 中回复了 CFD中文网 最后由 编辑
    #21

    @cfd-china 好的,我这么设置周期边界没有问题吧,两个面的距离是22mm

    | =========                 |                                                 |
    | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
    |  \\    /   O peration     | Version:  3.0.1                                 |
    |   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
    |    \\/     M anipulation  |                                                 |
    \*---------------------------------------------------------------------------*/
    FoamFile
    {
        version     2.0;
        format      ascii;
        class       polyBoundaryMesh;
        location    "constant/polyMesh";
        object      boundary;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    
    8
    (
        hub
        {
            type wall;
            nFaces          30748;
            startFace       5380555;
        }
        inlet
        {
            type patch;
            nFaces          5546;
            startFace       5411303;
        }
        outlet
        {
            type patch;
            nFaces          4543;
            startFace       5416849;
        }
        per1
        {
              type            cyclicAMI;
            nFaces          15340;
            startFace       5421392;
            neighbourPatch   per2;
            matchTolerance    0.01;
    	transform         translational;
    	separationVector (0 0.022 0);
        }
        per2
        {
             type            cyclicAMI;
           nFaces          15340;
           startFace       5436732;
           neighbourPatch   per1;
           matchTolerance    0.01;
           transform         translational;
           separationVector (0 -0.022 0);
        }
        ps
        {
            type            wall;
            nFaces          10266;
            startFace       5452072;
        }
        ss
        {
            type            wall;
            nFaces          11151;
            startFace       5462338;
        }
        sym
        {
            type symmetryPlane;
            nFaces          30748;
            startFace       5473489;
        }
    )
    
    // ************************************************************************* //
    1 条回复 最后回复
  • C 离线
    C 离线
    CFD中文网
    写于 最后由 编辑
    #22

    为何不用cyclic?

    per1
        {
              type            cyclic;
            nFaces          15340;
            startFace       5421392;
            neighbourPatch   per2;
        }
        per2
        {
             type            cyclic;
           nFaces          15340;
           startFace       5436732;
           neighbourPatch   per1;
        }
    
    
    
    

    CFD中国标准用户测试帐号
    目前由徐笑笑登录

    小 2 条回复 最后回复
  • 小 离线
    小 离线
    小葱小虫
    在 中回复了 CFD中文网 最后由 编辑
    #23

    @cfd-china 我感觉这个我理解的可能有问题,是应该给平移距离的大小(0.022m),还是(0,1,0)这么给。当时从这里看见的0_1490441899214_cyclic.jpg

    separationVector (0 -0.022 0);
    1 条回复 最后回复
  • 小 离线
    小 离线
    小葱小虫
    在 中回复了 CFD中文网 最后由 编辑
    #24

    @cfd-china 使用cyclic会有这个问题,因为两侧网格并不完全匹配

    /*---------------------------------------------------------------------------*\
    | =========                 |                                                 |
    | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
    |  \\    /   O peration     | Version:  3.0.1                                 |
    |   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
    |    \\/     M anipulation  |                                                 |
    \*---------------------------------------------------------------------------*/
    Build  : 3.0.1-119cac7e8750
    Exec   : simpleFoam
    Date   : Mar 25 2017
    Time   : 05:33:37
    Host   : "ubuntu"
    PID    : 48984
    Case   : /home/cc/OpenFOAM/cc-3.0.1/run/tutorials/incompressible/simpleFoam/cynasa
    nProcs : 1
    sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
    fileModificationChecking : Monitoring run-time modified files using timeStampMaster
    allowSystemOperations : Allowing user-supplied system call operations
    
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    Create time
    
    Create mesh for time = 0
    
    
    
    --> FOAM FATAL ERROR: 
    face 47 area does not match neighbour by 10.4554% -- possible face ordering problem.
    patch:per1 my area:1.63789e-10 neighbour area:1.47515e-10 matching tolerance:0.1
    Mesh face:5421439 fc:(0.0411928 0.00446869 5e-07)
    Neighbour fc:(0.040347 0.0264296 5e-07)
    If you are certain your matching is correct you can increase the 'matchTolerance' setting in the patch dictionary in the boundary file.
    Rerun with cyclic debug flag set for more information.
    
        From function cyclicPolyPatch::calcTransforms()
        in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 221.
    
    FOAM exiting
    1 条回复 最后回复
  • 小 离线
    小 离线
    小葱小虫
    写于 最后由 编辑
    #25

    @cfd-china 我的算例可以算了,但是不知道为什么,壁面的分离算的特别小

    小 1 条回复 最后回复
  • 小 离线
    小 离线
    小葱小虫
    在 中回复了 小葱小虫 最后由 编辑
    #26

    @小葱小虫 我只是换了一套网格,其他没动

    1 条回复 最后回复
  • 小 离线
    小 离线
    小葱小虫
    写于 最后由 编辑
    #27

    @cfd-china 老师您好,由于我的算例进口马赫数为0.7,所以要考虑可压缩性,但是算着算着就发散了,用的rhosimplefoam,请问,能否帮我调试一下,由于我的上面是对称面,下面是壁面,所以网格10万很难满足壁面要求,能否帮忙调试,谢谢

    I 1 条回复 最后回复
  • I 离线
    I 离线
    Izumi
    在 中回复了 小葱小虫 最后由 编辑
    #28

    @小葱小虫 你好!我在计算压缩机时,也出现了time step continuity errors过大导致发散的情况,参考苏军伟的博客,我调整了非正交修正次数,目前没发散,得到了结果。你可以参考一下。
    主要就是修改fvsolution文件。

    SIMPLE
    {
        //nNonOrthogonalCorrectors 0;
        nNonOrthogonalCorrectors 1;
        residualControl
        {
            p               1e-2;
            U               1e-3;
            "(k|epsilon)"   1e-3;
        }
    }
    小 1 条回复 最后回复
  • 小 离线
    小 离线
    小葱小虫
    在 中回复了 Izumi 最后由 编辑
    #29

    @Izumi 是增加nNonOrthogonalCorrector这个数吗,谢谢我试试

    H 1 条回复 最后回复
  • H 离线
    H 离线
    Haining LUO
    在 中回复了 小葱小虫 最后由 编辑
    #30

    @小葱小虫
    @Izumi 应该说的就是nNonOrthogonalCorrectors

    另外你两个网格之间差别是什么,我也是用pointwise画网格,出现aspect ratio那么大可能是手误也有可能是中间有breakpoint,你用examine aspect ratio看看那些畸形的网格,用connector->distribution->delete breakpoint去掉。因为如果一个方向的网格意外的小,虽然你的求解器是simpleFoam没有CFL数,但也对稳定性有影响吧?

    1 条回复 最后回复
  • 宝 离线
    宝 离线
    宝丁
    写于 最后由 编辑
    #31

    我也跟你遇到了相同的问题,不知道是不是计算设置的问题,总压入口是不是还有别的需要注意的地方

    1 条回复 最后回复
  • 宝 离线
    宝 离线
    宝丁
    写于 最后由 编辑
    #32

    你好,你现在解决这个问题了吗?我现在也遇到了这个问题,不过是因为pressureDirectedInletVelocity 这个进口条件的问题。
    我也在做类似的研究,可以交流一下吗

    1 条回复 最后回复

  • 登录

  • 登录或注册以进行搜索。
  • 第一个帖子
    最后一个帖子
0
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]