Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新
    1. Home
    2. random_ran
    R
    • Profile
    • Following
    • Followers
    • Topics
    • Posts
    • Groups

    random_ran

    @random_ran

    副教授

    205
    Posts
    3452
    Profile views
    25
    Followers
    0
    Following
    Joined Last Online

    random_ran Follow
    副教授

    Best posts made by random_ran

    • RE: 有关CFDEM安装

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

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

      • 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
      

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

      posted in OpenFOAM
      R
      random_ran
    • RE: 圆柱扰流计算与物理现象问题思考

      @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的频率值正好吻合。

      posted in OpenFOAM
      R
      random_ran
    • RE: sampleDict 是否可以切取圆柱面采集数据?

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

      Filter -> PlotOnIntersectionCurve

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

      posted in OpenFOAM
      R
      random_ran
    • RE: 并行效率疑问

      @xpqiu

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

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

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

      如果单纯从均值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边界条件能更好模拟三维圆柱绕流这一经典问题。

      posted in OpenFOAM
      R
      random_ran
    • RE: 并行效率疑问

      @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
      posted in OpenFOAM
      R
      random_ran
    • RE: Cd计算不准的问题-LES-Re3900-三维圆柱绕流-pisoFoam

      最开始确实是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:

      posted in OpenFOAM
      R
      random_ran
    • RE: 大网格-奇怪报错

      @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;
          }
      

      目前运行了半个小时了。

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

      posted in OpenFOAM
      R
      random_ran
    • RE: 新手求指导

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

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

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

      然后在 /cavity 目录下

      blockMesh
      

      祝你好运。

      posted in OpenFOAM
      R
      random_ran
    • RE: 稳态求解器simpleFoam的迭代问题

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

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

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

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

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

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

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

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

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

      posted in OpenFOAM
      R
      random_ran

    Latest posts made by random_ran

    • RE: Fluent虹吸VOF模拟

      跑了一个虹吸的模拟,挺有意思的,能模拟一些生活中的流体现象。

      https://github.com/randomwangran/capillary-will-rise

      以把重力翻过来的话,就和楼主的模拟比较接近了。

      posted in Fluent
      R
      random_ran
    • RE: k-omega 模型的网格无关性验证

      Sorry I had hard time to type Chinese.

      Thanks for sharing this. Your finding is consistent with what I've
      read this in Wilcox's book (P98).

      In free shear flow, k-epsilon model is unaffected by the numerical
      value of $k$ and $\epsilon$, but k-omega is very sensitive to the
      choice of $\omega$.

      But the book does not provide any suggestions on how to achieve mesh
      independence on k-omega model.

      Speaking to mesh independence study. The idea is simple: refining the
      mesh at the region with high gradient. In reality, it depends on
      researchers to define where are these regions. Even in the simple
      context like flow over a circular cylinder. Different people choose
      different types of the computational domain, different strategies to
      generate the mesh, different parameters to build the block.

      There are something I would investigate:

      1. continue refine the mesh in the "slot exit"
      2. expand the computational domain
      3. expand the "slot exit" region, especially where the shear layers
        are.
      4. use unsteady solver
      5. try LES if you have computational resource

      Wilcox, D. C., & others, (1998). Turbulence modeling for cfd. : DCW
      industries La Canada, CA.

      posted in Fluent
      R
      random_ran
    • RE: 分享:OpenFOAM不同离散格式的简单研究

      @波流力 非常感谢,我放在了这个仓库中,会做进一步步的研究,有值得分享的我会回来盖楼的.

      posted in OpenFOAM
      R
      random_ran
    • RE: OpenFOAM如何统计涡脱落频率?

      求 St 的方法有很多,比如对阻力或者升力系数 fft,或者在wake区域对某个参考点速度检测,然后在垂直于来流方向做fft。需要注意的是对Cd做fft所得的St应该是对Cl的两倍才对.

      ref: https://cfd-china.com/topic/774/圆柱扰流计算与物理现象问题思考/2

      posted in OpenFOAM
      R
      random_ran
    • RE: 圆柱扰流计算与物理现象问题思考

      @hoversoar

      还在,:).

      时间离散对结果影响大吗,还有Crank-Nicolson后面的系数应该怎么取?

      Please see: https://cfd-china.com/topic/3834/分享-openfoam不同离散格式的简单研究

      posted in OpenFOAM
      R
      random_ran
    • RE: 分享:OpenFOAM不同离散格式的简单研究

      感谢楼主分享,这些数据会让大家少走很多弯路.

      我有一个想法,不知楼主能否分享这个工做中 OpenFOAM 的一个算例文件?
      后来的人可以基于楼主研究,进一步研究相关问题.比如,在网格质量比较差情的况下,会有怎么样的影响.或者计算效率之类的.

      posted in OpenFOAM
      R
      random_ran
    • RE: CFD中文网很多功能取消了?

      可以先用谷歌代替:

      site:https://cfd-china.com/ pisoFoam

      posted in C斯达克
      R
      random_ran
    • RE: 未来的路该怎么走呢?

      深感这三年的付出如同泡沫一般付诸东流,自己无论怎么折腾都只是在瞎蹦跶。

      可是,写过的每一行代码,读过的每一本书,推过的每个一个公式,都是你的啊。

      国外的世界诚然宽阔,可是不及你心中的世界。

      胡思乱想的时候,跑跑步吧,刷 500 圈 400 米跑道;游 500圈 25 米的泳池;骑 500 公里的山路,之后,

      我发现,哪怕平平淡淡的生活,都是无比幸福的。

      posted in C斯达克
      R
      random_ran
    • RE: LES初始速度场设置和壁面函数设置

      @hoversoar

      上次算到step2的中心差分的50.5s处,说可能会阻力又下降。

      建议楼主定一个准则.比如计算了 N 个 Vortex Shedding 周.之后又计算了 2N 个周期,然后比较一下所感兴趣的物理量.

      那之后我把它继续算到了59.5s,然后速度对流项的linear换到了
      limitedLinear。把结果贴上来一下。

      不清楚楼主研究的重点是?层主个人来说,是很好奇, linear -> limitedLinear 的区别.

      整体看起来的话,三个阶段的定量CdCl差距不是很大。和文献也对的比较上。
      初 步认定是速度对流项选择了迎风格式造成的。跟之前的猜想有关,在55s处残
      差的Uy震荡明显,造成了Cd上升,Cl振幅也增大. 换成limitedLinear 0.5 之后,
      可以发现残差有明显的下降。

      不知楼主有没有对这些格式的比较?

      UDS 对网格要求更高.

      "Peaks or rapid variations in the variables will be smeared out and,
      because the rate of error reduction is only first order,
      very fine grids are required to obtain accurate solutions."

      楼主还可以探索一下这些:

      • Linear Interpolation (CDS)
      • Quadratic Upwind Interpolation (QUICK)
      • Higher-Order Schemes

      @random_ran 是否可以回答前面的问题,虽然P的量级并没有发生变化。但是Uz
      的曲线震荡明显(相比中心差分),猜想可能是由于TVD 里头那个0.5造成的,
      混入了一部分迎风格式,造成耗散?

      我没法准确回答你,因为我自己没有太多探索过.

      我这样中途切换差分格式是否ok?

      为什么不可以?如果是研究格式的问题,我可能也会这么做,省很多时间.但是如果是分析流体机理,最好从头开始用一种固定格式.

      posted in OpenFOAM
      R
      random_ran
    • RE: LES初始速度场设置和壁面函数设置

      今天才把case文件整理了一下,下一步找合适的网格.

      @东岳 岳兄过奖了.
      @hoversoar 今天才整理了一下你的 case.

      乃有新发现的话留言在这,大家都能看到,或者帮助你,我只能尽力,还是要靠氖自己的.

      posted in OpenFOAM
      R
      random_ran