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中文网

R

random_ran

@random_ran
大神
关于
帖子
206
主题
12
群组
1
粉丝
35
关注
0

帖子

最新

  • 有关CFDEM安装
    R random_ran

    官方的安装指南通常来说是最详细的。

    我自己试着安装了一下这种需要联合编译的软件,这不是一个成功的安装日志只是记录一下自己的遇到的问题。

    • linux 服务器:

    LSB Version: :base-4.0-amd64:base-4.0-noarch:core-4.0-amd64:core-4.0-noarch:graphics-4.0-amd64:graphics-4.0-noarch:printing-4.0-amd64:printing-4.0-noarch
    Distributor ID: CentOS
    Description: CentOS release 6.4 (Final)
    Release: 6.4
    Codename: Final

    • OpenFOAM
      预安装,module load 加载,版本: openfoam/3.0.1
    • 无root权限

    这个安装网页,没有提到其它的依赖库,似乎只要有O.F. 相应的版本就可以了。

    下面我一步步来尝试安装。

    1. 首先加载OpenFOAM
    2. 通过git 下载 LIGGGHTS
    $ cd
    $ mkdir LIGGGHTS
    $ cd LIGGGHTS 
    $ git clone git://github.com/CFDEMproject/LIGGGHTS-PUBLIC.git LIGGGHTS-PUBLIC
    

    同样安装 CFDEM(R)

    $ cd
    $ mkdir LIGGGHTS
    $ cd LIGGGHTS 
    $ git clone git://github.com/CFDEMproject/LIGGGHTS-PUBLIC.git LIGGGHTS-PUBLIC
    $ git clone git://github.com/CFDEMproject/CFDEMcoupling-PUBLIC.git CFDEMcoupling-PUBLIC-$WM_PROJECT_VERSION
    

    获取库中最新文件

    cd $HOME/CFDEM/CFDEMcoupling-PUBLIC-$WM_PROJECT_VERSION
    git stash  //小心这个操作会擦出本地修改的代码
    git pull
    

    设置环境变量:

    $ cd
    $ emacs .bashrc
    

    在.bashrc中写入

    #================================================#                                                                                                                                               
    #- source cfdem env vars                                                                                                                                                                         
    export CFDEM_VERSION=PUBLIC
    export CFDEM_PROJECT_DIR=$HOME/CFDEM/CFDEMcoupling-$CFDEM_VERSION-$WM_PROJECT_VERSION
    export CFDEM_SRC_DIR=$CFDEM_PROJECT_DIR/src
    export CFDEM_SOLVER_DIR=$CFDEM_PROJECT_DIR/applications/solvers
    export CFDEM_DOC_DIR=$CFDEM_PROJECT_DIR/doc
    export CFDEM_UT_DIR=$CFDEM_PROJECT_DIR/applications/utilities
    export CFDEM_TUT_DIR=$CFDEM_PROJECT_DIR/tutorials
    export CFDEM_PROJECT_USER_DIR=$HOME/CFDEM/$LOGNAME-$CFDEM_VERSION-$WM_PROJECT_VERSION
    export CFDEM_bashrc=$CFDEM_SRC_DIR/lagrangian/cfdemParticle/etc/bashrc
    export CFDEM_LIGGGHTS_SRC_DIR=$HOME/LIGGGHTS/LIGGGHTS-PUBLIC/src
    export CFDEM_LIGGGHTS_MAKEFILE_NAME=fedora_fpic
    export CFDEM_LPP_DIR=$HOME/LIGGGHTS/mylpp/src
    export CFDEM_PIZZA_DIR=$HOME/LIGGGHTS/PIZZA/gran_pizza_17Aug10/src
    . $CFDEM_bashrc
    #================================================#         
    

    重新登录,信息提示:

    using default CFDEM_ADD_LIBS_DIR=/home/userID/CFDEM/CFDEMcoupling-PUBLIC-3.0.1/src/lagrangian/cfdemParticle/etc/addLibs_universal
    using default CFDEM_ADD_LIBS_NAME=additionalLibs_3.0.1
    !!! ERROR !!!: CFDEM_ADD_LIBS_DIR/CFDEM_ADD_LIBS_NAME=/home/userID/CFDEM/CFDEMcoupling-PUBLIC-3.0.1/src/lagrangian/cfdemParticle/etc/addLibs_universal/additionalLibs_3.0.1 does not exist.
    make new dirs /home/userID/CFDEM/userID-PUBLIC-3.0.1 ? (y/n)
    

    检查安装:

    $CFDEM_PROJECT_DIR
    $CFDEM_SRC_DIR
    $CFDEM_LIGGGHTS_SRC_DIR 
    

    出现三次:

    /home/userID/CFDEM/CFDEMcoupling-PUBLIC-3.0.1: is a directory
    

    检查pizza的地址?

    $ cd $CFDEM_SRC_DIR/lagrangian/cfdemParticle/etc/
    $ sh cfdemSystemTest.sh
    

    检查结果

    
    *********************************************
    * C F D E M (R) c o u p l i n g             *
    *                                           *
    * by DCS Computing GmbH                     *
    * www.dcs-computing.com                     *
    *********************************************
    
    *********************************
    CFDEM(R)coupling system settings:
    *********************************
    CFDEM_VERSION=PUBLIC
    couple to OF_VERSION=3.0.1
    compile option=Opt
    
    check if paths are set correctly
    valid:yes critical:yes - $CFDEM_PROJECT_DIR = /home/userID/CFDEM/CFDEMcoupling-PUBLIC-3.0.1
    valid:NO  critical:no - $CFDEM_PROJECT_USER_DIR = /home/userID/CFDEM/userID-PUBLIC-3.0.1 does not exist
    valid:yes critical:yes - $CFDEM_SRC_DIR = /home/userID/CFDEM/CFDEMcoupling-PUBLIC-3.0.1/src
    valid:yes critical:yes - $CFDEM_SOLVER_DIR = /home/userID/CFDEM/CFDEMcoupling-PUBLIC-3.0.1/applications/solvers
    valid:yes critical:yes - $CFDEM_TUT_DIR = /home/userID/CFDEM/CFDEMcoupling-PUBLIC-3.0.1/tutorials
    valid:yes critical:yes - $CFDEM_LIGGGHTS_SRC_DIR = /home/userID/LIGGGHTS/LIGGGHTS-PUBLIC/src
    valid:NO  critical:yes - $CFDEM_LPP_DIR = /home/userID/LIGGGHTS/mylpp/src does not exist
    valid:yes critical:yes - $CFDEM_ADD_LIBS_DIR = /home/userID/CFDEM/CFDEMcoupling-PUBLIC-3.0.1/src/lagrangian/cfdemParticle/etc/addLibs_universal
    valid:NO  critical:no - $CFDEM_PIZZA_DIR = /home/userID/LIGGGHTS/PIZZA/gran_pizza_17Aug10/src does not exist
    valid:NO  critical:no - $CFDEM_TEST_HARNESS_PATH = /home/userID/CFDEM/userID-PUBLIC-3.0.1/log/logFilesCFDEM-PUBLIC-3.0.1 does not exist
    valid:NO  critical:no - $C3PO_SRC_DIR =  does not exist
    
    library names
    $CFDEM_LIGGGHTS_LIB_NAME = lmp_fedora_fpic
    $CFDEM_LIB_NAME = lagrangianCFDEM-PUBLIC-3.0.1
    $LD_LIBRARY_PATH  = /global/scratch/userID/OpenFOAM/userID-3.0.1/platforms/linux64IccDPInt32Opt/lib:/global/software/OpenFOAM-3.0.1/ThirdParty-3.0.1/platforms/linux64IccDPInt32Opt/lib/openmpi-system:/global/software/OpenFOAM-3.0.1/platforms/linux64IccDPInt32Opt/lib/openmpi-system:/global/software/OpenFOAM-3.0.1/ThirdParty-3.0.1/platforms/linux64IccDPInt32Opt/lib:/global/software/OpenFOAM-3.0.1/platforms/linux64IccDPInt32Opt/lib:/global/software/OpenFOAM-3.0.1/platforms/linux64IccDPInt32Opt/lib/dummy:/global/software/intel/composerxe/mkl/lib/intel64:/global/software/intel/composerxe/lib/intel64:/global/system/globus-5.2.5/lib64
    $WM_NCOMPPROCS  = 
    $WM_LABEL_SIZE = 32
    
    Additional lib settings
    .Makefile_vtk_tmp:1: /home/userID/CFDEM/CFDEMcoupling-PUBLIC-3.0.1/src/lagrangian/cfdemParticle/etc/addLibs_universal/additionalLibs_3.0.1: No such file or directory
    make: *** No rule to make target `/home/userID/CFDEM/CFDEMcoupling-PUBLIC-3.0.1/src/lagrangian/cfdemParticle/etc/addLibs_universal/additionalLibs_3.0.1'.  Stop.
    
    *******************
    g++:
    /usr/bin/g++
    g++ (GCC) 4.4.7 20120313 (Red Hat 4.4.7-3)
    Copyright (C) 2010 Free Software Foundation, Inc.
    This is free software; see the source for copying conditions.  There is NO
    warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
    
    gcc:
    /usr/bin/gcc
    gcc (GCC) 4.4.7 20120313 (Red Hat 4.4.7-3)
    Copyright (C) 2010 Free Software Foundation, Inc.
    This is free software; see the source for copying conditions.  There is NO
    warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
    
    mpic++:
    /global/software/openmpi-1.6.5/intel/bin/mpic++
    icpc (ICC) 14.0.2 20140120
    Copyright (C) 1985-2014 Intel Corporation.  All rights reserved.
    
    mpirun:
    /global/software/openmpi-1.6.5/intel/bin/mpirun
    mpirun (Open MPI) 1.6.5
    
    Report bugs to http://www.open-mpi.org/community/help/
    **********************
    additional packages...
    
    

    结果显示

    $CFDEM_PROJECT_USER_DIR
    $CFDEM_PIZZA_DIR
    $CFDEM_TEST_HARNESS_PATH
    $C3PO_SRC_DIR
    

    不存在,不知道有什么影响?

    编译LIGGGHTS(R) and CFDEM(R)coupling:
    回到home:

    git clone git://github.com/CFDEMproject/CFDEMcoupling-PUBLIC.git CFDEMcoupling-PUBLIC-$WM_PROJECT_VERSION
    

    注意,这个git clone从手册中似乎不全,我是到他的rep查到的准确地址。

    编译LIGGGHTS(R)

    安装过程很快,但是我注意到有一个处有错误

    using CFDEM_LAMMPS_LIB_DIR=/home/userID/LIGGGHTS/LIGGGHTS-PUBLIC/src/../lib defined by user.
    using CFDEM_ADD_LIBS_DIR=/home/userID/CFDEM/CFDEMcoupling-PUBLIC-3.0.1/src/lagrangian/cfdemParticle/etc/addLibs_universal defined by user.
    using CFDEM_ADD_LIBS_NAME=additionalLibs_3.0.1 defined by user.
    !!! ERROR !!!: CFDEM_ADD_LIBS_DIR/CFDEM_ADD_LIBS_NAME=/home/userID/CFDEM/CFDEMcoupling-PUBLIC-3.0.1/src/lagrangian/cfdemParticle/etc/addLibs_universal/additionalLibs_3.0.1 does not exist. 
    

    有一个lib似乎不存在。

    接着编译CFDEM(R)coupling

    $ cd
    $ cfdemCompCFDEM
    

    然后就出现了很多error,由于太多,我就不贴了,注意到几处细节:

    Please provide the utilities to be compiled in the /global/home/userID/CFDEM/CFDEMcoupling-PUBLIC-3.0.1/src/lagrangian/cfdemParticle/etc/utilities-list.txt file.
    structure:
    path  to provide the path relative to CFDEM_UT_DIR
    
    example:
    cfdemPostproc/dir
    

    似乎要编译的话,要在

    /global/home/userID/CFDEM/CFDEMcoupling-PUBLIC-3.0.1/src/lagrangian/cfdemParticle/etc/utilities-list.txt
    

    指定一些具体的求解器, 安装说明书到这里就没有多写了。这次不成功的安装也就到这吧。


  • 圆柱扰流计算与物理现象问题思考
    R random_ran

    @yhdthu

    3900 已经是sub-critical了,这个区间的雷诺数已经有3D效应了,建议extrude你的2D 网格,在轴向上拉长 pi 的长度并布置48或者更高层网格。

    对于这个St从Cd和Cl的计算产生的差异,我觉得这样理解比较好:

    想象钟表平面是那个圆柱截面,3点钟方向是流体流动方向,也就是Cd的方向,12点钟是升力系数为正的方向。这样规定好座标轴以后,我们在讨论作用在圆柱体表面的力的周期性变化。

    你可以想象钟表上的时针代表合力的方向,合力的周期性变化表现在时针的来回转动。

    • 当1点钟的时候,我们假定一个周期开始,这个时候Cl出现了最大正值。
    • 当3点钟的时候,Cd出现最大正值。
    • 5点钟的时候Cl,最大负值。

    从1点到5点,经历了半个周期,你会发现:
    Cl从最大正值逐渐变为最小负值,这就是正好是vortex shedding 的半个周期。

    再来看看Cd:
    Cd在1点钟的时候,我们要把Cd的真值是1点钟投影到3点钟上的大小,所以Cd此时是最小值,从1点钟到3点钟,Cd逐渐增大到 最大值(此时此刻Cd不用投影)。

    接下来从3点钟到5点钟,Cd又逐渐减小到最小值,同样需要在3点钟方向上做投影。

    总结一下:

    从1点到5点,这个半个vortex shedding 周期内:

    • Cl 也是半个周期
    • Cd 经历了一个全周期

    所以依据Cd计算的St大小应该是依据Cl的两倍。而基于Cl得出的St,和vortex shedding的频率值正好吻合。


  • sampleDict 是否可以切取圆柱面采集数据?
    R random_ran

    或许可以等计算完成之后,用ParaVIEW后处理:

    Filter -> PlotOnIntersectionCurve

    我做圆柱界面的压力分布历史就用的这个功能。


  • 并行效率疑问
    R random_ran

    @xpqiu

    非常感谢你的建议, renumberMesh 在一个800万网格 96核的模拟中,时间缩短了40%。

    以后多核计算前,我都会renumberMesh 的,也建议大家用上这个好用的功能,加速你的算例。:D


  • Cd计算不准的问题-LES-Re3900-三维圆柱绕流-pisoFoam
    R random_ran

    如果单纯从均值Cd的大小上看,这两个边界条件,对于我所用的那个网格并没有太大影响。

    我最开始的时候也是用过symmetry的边界条件,但是我觉得slip-wall 的边界条件可能更适合描述我当初对这个问题的理解。

    后来我发现这个研究方向上的文献中,绝大多数人在span-wise都是用cyclic的边界条件。用cyclic边界的条件的优势之一就是能把模拟的有限长度的圆柱,当成无限长来出来。换句话说就是end-effect上有很强的优势。这个即便是做风洞实验都无法克服的问题,因为风洞中都要用end-plate来矫正端部对流场的影响。

    Yeo, DongHun, and Nicholas P. Jones. “Investigation on 3-D Characteristics of Flow around a Yawed and Inclined Circular Cylinder.” Journal of Wind Engineering and Industrial Aerodynamics 96.10–11 (2008): 1947–1960. Web.

    这篇文章中3.1节讨论了 slip-wall 和 cyclic wall 对涡结构的影响。作者发现在cyclic边界条件受边界的影响要小,因此cyclic边界条件能更好模拟三维圆柱绕流这一经典问题。


  • 并行效率疑问
    R random_ran

    @cfd-china

    你说得很对。从hierarchical到scotch是我要考虑的一个方向。不过由于decompose的时间花费,现在我在处理一些solver上的选择,这样可以不decompse。这个表格是我现在做过的一些case。

    我的目标计算时间是1个小时算出0.5second。

    #case p-solver #核心数 decompose 方法 数值精度 求解器 1小时后,模拟时间 湍流模型 #PISO 矫正次数
    1 GAMG+symGaussSeidel 576 hiera (32 9 2) -6 pisoFoam - LES-kqE 3
    2 - - - - - - - -
    3 PCG+DIC 576 hiera (32 9 2) -6 pisoFoam 0.133 LES-WALE 3
    4 PCG+DIC 576 hiera (32 9 2) -6 pisoFoam 0.129 LES-kqE 3
    5 PCG+DIC 576 hiera (32 9 2) -6 pisoFoam 0.125 LES-TKE 3
    6 PCG+DIC 576 hiera (32 9 2) -7 pisoFoam 0.125 LES-TKE 3
    7
    8
    9 GAMG+smootherGaussSeidel 576 hiera (4 4 36) -7 pisoFoam 0.015(2914s) LES-TKE 3
    10 GAMG+smootherGaussSeidel 288 hiera (16 9 2) -6 pisoFoam 0.171 LES-TKE 2
    11 GAMG+smootherGaussSeidel 288 hiera (16 9 2) -6 pisoFoam 0.132 LES-TKE 3

  • Cd计算不准的问题-LES-Re3900-三维圆柱绕流-pisoFoam
    R random_ran

    最开始确实是DES。 dynamicKEqution 的一大特点就是波动性,我第一次看到这个湍流模型运行出来的Cd历史的时候也被惊讶到了。可以参考:

    Lysenko, Dmitry A., Ivar S. Ertesvåg, and Kjell Erik Rian. “Large-Eddy Simulation of the Flow over a Circular Cylinder at Reynolds Number 3900 Using the Openfoam Toolbox.” Flow, Turbulence and Combustion 89.4 (2012): 491–518. Web.

    如果湍流结果是一个前置条件的话,特别是你的雷诺数上千以后,2D的模拟总会计算出一个较高的Cd结果。 我不知道你是怎么计算Cd,O.F. 自己有专门计算Cd的object:

    FoamFile
    {
        version     2.0;
        format      ascii;
        class       dictionary;
        object      forceCoeffsIncompressible;
    }
    

    你可以按自己的需求调整输出频率。另外,文献中对采样点的多少也是各有各的看法。有的人(Franke 2002)采样150个周期(shedding period), 有的人(Kravchenko and Moin 2000)说采样10几个周期的结果和自己的实验结果很吻合,而更多的周期则会出现偏差。个人建议30~50个周期。

    我这里有最开始那个DES的Cp图,传统的plot方法。后面的因为考察的方向不同,没有用传统的方式画,对你也没有参考价值。Lysenko 2012 年那篇里有Cp的plot,你可以参考他的。

    Cp的做法,我是用ParaView 后处理得到的。先在OF设置好:在一个周期内,输出,比如20个点 (包含p,U场)。计算完成后在圆柱的中截面截取出压力值用来计算Cp (这个Cp是截面Cp值),这和Cd (整个圆柱)的算法是有区别的。得到了一个个的输出的时间点的Cp值之后,再对所有的采样的周期,比如30个周期,求出平均值。

    ps 论坛的回复没有办法插入图片,如果你有需要我那个Cp,可以发我的邮箱.

    个人习惯:不用微信以及一切即时通信软件。这是并行计算的哲学告诉我的:尽量减少cpu之间的通信以提升效率。 :joking:


  • 大网格-奇怪报错
    R random_ran

    @cfd-china

    非常感谢你的建议!非常管用,起码现在能跑起来了。

    
    Time = 0.005
    
    Courant Number mean: 0.000119441 max: 0.017568
    DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 2.34151e-11, No Iterations 2
    DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 1.61665e-10, No Iterations 2
    DILUPBiCG:  Solving for Uz, Initial residual = 1, Final residual = 1.63882e-10, No Iterations 2
    GAMG:  Solving for p, Initial residual = 1, Final residual = 0.035832, No Iterations 36
    time step continuity errors : sum local = 8.55925e-06, global = -9.25811e-08, cumulative = -9.25811e-08
    GAMG:  Solving for p, Initial residual = 0.000128631, Final residual = 6.05672e-06, No Iterations 39
    time step continuity errors : sum local = 4.27446e-07, global = -5.85326e-09, cumulative = -9.84343e-08
    GAMG:  Solving for p, Initial residual = 8.87864e-06, Final residual = 9.28164e-08, No Iterations 74
    time step continuity errors : sum local = 6.57138e-09, global = 9.95505e-13, cumulative = -9.84333e-08
    DILUPBiCG:  Solving for k, Initial residual = 0.297752, Final residual = 7.56559e-07, No Iterations 4
    bounding k, min: -2.24189e-08 max: 2e-05 average: 2.54114e-09
    ExecutionTime = 42.97 s  ClockTime = 43 s
    

    我改动的fvSolution 部分是:

        "(U|k|B|nuTilda|s)"
        {
            solver          PBiCG;
            preconditioner  DILU
            tolerance       1e-07;
            relTol          0;
        }
    

    目前运行了半个小时了。

    我知道后面还会有的新的问题。


  • 新手求指导
    R random_ran

    最基本的问题在论坛上问不是特别有效率的。熟悉官方教程,你就可以找到答案。

    另外,熟悉一些最基本的 linux 语言会对你学习OpenFOAM大有帮助。

    cp -r $FOAM_TUTORIALS/incompressible/icoFoam/cavity/cavity .  
    

    然后在 /cavity 目录下

    blockMesh
    

    祝你好运。


  • 稳态求解器simpleFoam的迭代问题
    R random_ran

    大神不敢当,我也是和你一样刚接触CFD的小白。

    不好意思对于 continuity error 的问题没能讲得很明白,这是因为我对这个问题没有理解到“代码”级别,不过很多时候,这个问题没有给我带来太多麻烦,所以就停留在这个地方。如果有熟悉这方面的朋友,还希望分享一下更好的见解。

    对于一个物理用瞬态求解还是稳态求解,也是一个很大的问题。可以说,我至今也没完全搞明白。不过最基本的区分,或许能帮助到你:

    1. 瞬态: 关注的对象是非常依赖时间的。pisoFoam 和 pimpleFoam 在此类。最好把时间上的离散用二阶以上精度。

    2. 稳态问题: 关注的物理量和时间无关。

    我在用大涡模拟的时候就用的是pisoFoam求解器。因为旋涡脱落过程中,有一个参量叫Strouhal数是我重点考察的一个对象之一。

    对于稳态计算,通常可以选择RNAS的湍流模型。一般来讲,时间步长要高于顺态求解,所以计算速度更快,计算消耗资源更小。

    个人理解,对于真实的物理世界,或多或少都和时间有关系。只是我们在研究问题的时候,关注点不一样,用不同的方法,高效地处理不同的问题。

    至于算法级别,归结到代码上的问题,实在无能为力。我也还在学习中。
    :sad:


  • ParaView-streamline 问题
    R random_ran

    #问题1 自我解决,感觉像是bug.
    0_1488141087886_streamline.png


  • Cd计算不准的问题-LES-Re3900-三维圆柱绕流-pisoFoam
    R random_ran

    @sorrby 改了很多地方,比如fvm 的scheme, 时间步长,收敛标准。

    对于之前那套网格,实在是没有精力了。一个二维的网格就16W的cell了,已经做得很密了。 反倒是圆形的这个计算域,只用到4W的cell就把很多参数算得很准了。要是有时间,再试试那个矩形的计算网格。

    另外,我的意思是这样的,对于相同的雷诺数,不同的流体性质,会不会对一些测量参数产生一定范围内的影响?比如10%以内。


  • LES初始速度场设置和壁面函数设置
    R random_ran

    @hoversoar 台兄过奖了。我只是好奇罢了。自己玩的很杂,可是遇到能一起玩的,实在不易,所以慢慢来。我在找 of 的脚本,自动生成。最好简单可控,小 mesh 大家可以互相交流,大 mesh 大家自己留着做研究就好。


  • LES初始速度场设置和壁面函数设置
    R random_ran

    看到有人讨论圆柱绕流就很激动。

    可是我已经挖了太多坑还没填。关注中,加油。


  • Cd计算不准的问题-LES-Re3900-三维圆柱绕流-pisoFoam
    R random_ran

    @Regina

    谢谢你对我的论文感兴趣。

    我自己来说,一开始是卡着 CFL =1 来的。不过也看到有人用 2。

    0.3 比较小的了,可以迈开步子!

    比如我目前在运行的一个例子
    Courant Number mean: 0.0522079 max: 0.977915

    够可以的吧?哈哈..


  • OpenFOAM后处理与Fluent结果比较
    R random_ran

    也想比对Fluent和O.F.的计算效率与准度的问题,但是发现实际操作起来还挺麻烦的。边界条件的处理,离散,数值积分方面的处理要做到非常接近才能让对比结果有说明意义。这篇文献供大家参考。

    Robertson, E., Choudhury, V., Bhushan, S. & Walters, D. K. Validation of OpenFOAM numerical methods and turbulence models for incompressible bluff body flows. Comput. Fluids 123, 122–145 (2015).


  • 最快速的上google的方式是什么?
    R random_ran

    不是秀网速, 怎么样科学上网才能达到这样的效果?
    0_1516073318925_test.png


  • [分享] 一个 LaTeX 论文模板
    R random_ran

    @l-j刘侃

    如果必须要用邮箱联系的时候,我会电邮你的。QQ已经3年没用过了。有问题讨论的话,我会在论坛看到的。有时候没有及时回复,我也会收藏起来,等稍微有空的时候回复的。

    尽量降低产生这种涡流的成因 或者是降低涡流的强度和范围

    一个设备中 当两股流体相碰的时候 比如 一股上升流和一个反弹回来的反弹流 能够在近壁面处形成一个非常强烈的垂直向的涡流

    我还是不太清楚这个涡流的形态,有没有什么关键词?或者一些重要文献?或者视频演示?

    如果没有搞清问题是什么的前提下,很难给你出有效的解决办法的。

    于研究不同形状表面对于圆柱后的涡流区的影响

    Zhang Kai 一个曾在日本的博士生,他对这个问题有过深入的探讨 在 Journal of Wind Engineering and Industrial Aerodynamics 也发过文章。他的博士值得看看。

    表面呈锯齿状 就好比高尔夫球

    我觉得,你指得是 drag crisis, 对于这种分离点很敏感的几何体, 层流边界层还是湍流边界层,对拽力升力有巨大的影响。

    改变几何体,很多情况下是把分离点固定,从而达到稳定流体的目的。

    要实现这个目的还有很多方式...

    我自己的论文中也对这个问题稍有涉足:

    Chapter 6 中,通过添加一个 分割板
    在 FIGURE 6.6 and 6.7 可以看到 升力和阻力随着板的特征的变化而变化的趋势。

  • 登录

  • 没有帐号? 注册

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