OpenFOAM

• F

  rhoCentralFoam版本更替的交流
  • Foamer24  

  3
  0
  Votes
  3
  Posts
  81
  Views

  F

  @sibo 厉害厉害，第一个问题理解啦
• F

  可压缩流求解器rhoEnergyFoam的下载以及代码适应 问题
  • Foamer24  

  2
  0
  Votes
  2
  Posts
  57
  Views

  F

  现在主要的问题是对于下面的式子编译时，会出现error: assignment of read-only location，这个百度上说不能通过该指针改变其指向的值。但是不知道下面的代码要怎么修改才可以。 const polyPatch& patch = mesh.boundaryMesh()[iPatch] ;
• L

  如何在每个时间步更新边界条件
  • liailei0627  

  5
  1
  Votes
  5
  Posts
  104
  Views

  L

  @Foamer24 感谢回答，您觉得Groovy可以解决我的问题， 一是添加扰动，二是基于流场状态反馈得到入口条件，并在每个时间步上添加。 对Swak4Foam没用过，还请赐教
• 红

  自定义函数求助
  • 红豆沙  

  2
  0
  Votes
  2
  Posts
  81
  Views

  红

  问题 稍微解决，，将symmetry边界换成wall边界，使用zeroGradient 可以算了，，先算 再看
• W

  openfoam中bounded和unbounded的区别
  • win  

  7
  0
  Votes
  7
  Posts
  267
  Views

  D

  代码中的描述的作用是为了使求解稳定，在原有格式的基础上减去了fvm::Sp(fvc::surfaceIntegrate(faceFlux), vf) ，而计算收敛的时候该项为0，所以理论上应该不会影响精度。
• J

  [菜鳥求救] K-omega SST 头几秒就崩溃了
  • jinjolee  

  18
  0
  Votes
  18
  Posts
  594
  Views

  D

  @jinjolee openFOam里的湍流模型有些时候确实不好用。（1）用你现在的条件计算出k、omega和nut的值，保证初始场initialValue也等于计算出的值；（2）检查壁面函数是否合适；（3）检查空间离散格式fvSchemes。
• 宝

  rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
  • 宝丁  

  6
  0
  Votes
  6
  Posts
  134
  Views

  宝

  @东岳 这个离SCI还很遥远，应该是两个求解器之间的差别吧，实验和计算对起来也很麻烦，我再做两个算例
• M

  twoPhaseEulerFoam 热传递
  • MX  

  21
  0
  Votes
  21
  Posts
  7111
  Views

  F

  @王孟飞 这个问题我也遇到了，不知道怎么解决
• 1

  使用rhoPimpleDyMFoam计算在开头就无法进行下去错误
  • 100yearsalone  

  7
  0
  Votes
  7
  Posts
  1424
  Views

  F

  @东岳 我非定常计算也出现pseudoTime：not converged within 20 iterations 每迭代一步，都要内迭代20次。不知道放宽收敛标准会不会容易发散
• D

  一个有毒的算例……
  • durm17  

  10
  0
  Votes
  10
  Posts
  1554
  Views

  D

  @东岳 后来就只用snappyHexMesh画的网格做了一些简单的模拟……
• Y

  求助：openfoam中的Immersed Boundary Methods
  • yuan_neu  

  6
  0
  Votes
  6
  Posts
  325
  Views

  红

  我也正在刚开始研究IBM，，不知道你现在做的如何了，可以交流交流么
• L

  overPimpleDyMFoam计算中出现报错
  • lv1995  

  1
  0
  Votes
  1
  Posts
  44
  Views

  L

  在使用overPimpleDyMFoam计算旋转机械 时，出现了如下报错： Create time Create mesh for time = 0 Selecting dynamicFvMesh dynamicOversetFvMesh Selecting motion solver: rigidBodyMotion Selecting rigidBodySolver Newmark PIMPLE: Operating solver in PISO mode Reading field p Reading field U Reading/calculating face flux field phi Creating cellMask field to block out hole cells inverseDistance : detected 2 mesh regions zone:0 nCells:29272 voxels:(193 193 1) bb:(-300.001077033 -200.001077033 1.99892296657) (700.001077033 200.001077033 3.00107703343) zone:1 nCells:8145 voxels:(193 193 1) bb:(-50.0001414007 -49.9939730307 1.99985859933) (50.0001414007 49.9720506607 3.00014140067) --> FOAM Warning : From function bool Foam::oversetPolyPatch::master() const in file oversetPolyPatch/oversetPolyPatch.C at line 149 The master overset patch is not the first patch. Generally the first patch should be an overset patch to guarantee consistent operation. Overset analysis : nCells : 37417 calculated : 37186 interpolated : 204 (interpolated from local:204 mixed local/remote:0 remote:0) hole : 27 Creating interpolatedCells field Selecting incompressible transport model Newtonian Selecting turbulence model type RAS Selecting RAS turbulence model kOmegaSST Selecting patchDistMethod meshWave RAS { RASModel kOmegaSST; turbulence on; printCoeffs on; alphaK1 0.85; alphaK2 1; alphaOmega1 0.5; alphaOmega2 0.856; gamma1 0.555555555556; gamma2 0.44; beta1 0.075; beta2 0.0828; betaStar 0.09; a1 0.31; b1 1; c1 10; F3 false; decayControl false; kInf 0; omegaInf 0; } Reading/calculating face velocity Uf No MRF models present No finite volume options present Courant Number mean: 0.0195783147649 max: 21.9831544632 Starting time loop Courant Number mean: 0.0195783147649 max: 21.9831544632 Time = 0.005 forces forces: Not including porosity effects --> FOAM FATAL ERROR: Could not find rho:rho From function void Foam::functionObjects::forces::initialise() in file forces/forces.C at line 204. FOAM exiting 使用的动网格文件如下： motionSolverLibs ("librigidBodyMeshMotion.so"); dynamicFvMesh dynamicOversetFvMesh; motionSolver rigidBodyMotion; // rigidBodyMotionCoeffs - not needed any more // { report on; solver { type Newmark; gamma 0.1; // Velocity integration coefficient beta 0.1; // Position integration coefficient } accelerationRelaxation 1.0; //- prescribedRotation requires some sub-iterations to converge faster // to desired value. nIter 3; bodies { rotor { type rigidBody; parent root; centreOfMass (-0.265 0 0); mass 1; rho 1; inertia (0.001 0 0 0.001 0 0.001); transform (1 0 0 0 1 0 0 0 1) (-0.265 0 0); patches (BLADE1); innerDistance 1000; outerDistance 2000; joint { type Rz; } } } restraints { rotorRotation { type prescribedRotation; body rudder; referenceOrientation (1 0 0 0 1 0 0 0 1); axis (0 0 1); // Axis of rotation omega sine; // Function1 entry omegaCoeffs { frequency 4; amplitude 0.2; scale (0 1 0); level (0 1 0); } } } 但是如果把动网格文件改为： dynamicFvMesh dynamicOversetFvMesh; dynamicOversetFvMeshCoeffs { // layerRelax 0.3; } solver multiSolidBodyMotionSolver; multiSolidBodyMotionSolverCoeffs { rotor { solidBodyMotionFunction rotatingMotion; rotatingMotionCoeffs { origin (0 0 0); axis (0 0 1); omega 1.08; } } } 就可以正常计算了。请问这是哪里出现错误了，跪求各位大神解答。
• C

  单个气泡传热理论解
  • CYW  

  6
  0
  Votes
  6
  Posts
  180
  Views

  桎

  @CYW 谢谢。在interFoam里加入过能量方程，但是在加物理量比如普康特数啥的出现了点问题。。。我对照着compressibleInterFoam研究研究，，有问题再向你请教。不好意思没帮到你什么忙
• 

  找到vof中interface的位置
  • linhan.ge  

  50
  0
  Votes
  50
  Posts
  9379
  Views

  W

  @linhan-ge 测试过，另外可以参看一下Johan Roenby的说法："For you surface tension people, I believe you will also need an improved estimate of the curvature, right? This information is somehow contained in the difference in surface normal direction between adjacent surface cells. Hmm... still thinking about how to cook this up in a consistent and robust way for general meshes..." 来自https://www.cfd-online.com/Forums/openfoam-news-announcements-other/180439-isoadvector-release.html
• 麦

  求助：液池热毛细对流算例，已编译好求解器和算例，但就是和文献自编程序不一样
  • 麦迪文  

  10
  0
  Votes
  10
  Posts
  353
  Views

  L

  @麦迪文 等着看云图
• 夏

  小白新人请教流固耦合是否适合使用OpenFOAM处理？
  • 夏雨天  

  2
  0
  Votes
  2
  Posts
  97
  Views

  T

  @夏雨天 在 小白新人请教流固耦合是否适合使用OpenFOAM处理？ 中说： immersed-boundary method openfoam-extend了解一下，immersed-boundary也有拓展包可以单独安装； https://sourceforge.net/p/openfoam-e...leaseNotes.txt https://openfoamwiki.net/index.php/Extend-bazaar
• F

  Openfoam如何高阶计算？
  • Foamer24  

  5
  0
  Votes
  5
  Posts
  226
  Views

  F

  上面的算法都有试过，差别不是很大，和实验数据还是有一些区别。现在在调试fvSolution里面的设置，希望有效果。
• L

  LamBremhorstKE模型
  • lambs  

  9
  0
  Votes
  9
  Posts
  1372
  Views

  

  @lambs 这个算例后来验证的如何了？
• 一

  LES介绍的文章的一个公式
  les • • 一二  

  10
  0
  Votes
  10
  Posts
  3570
  Views

  

  @东岳 东岳老师，有个小问题咨询一下。就是 \begin{equation} \nu_{\mathrm{SGS}}=\rho\left(C_{\mathrm{SGS}} \Delta \right)^2 \sqrt{\frac{1}{2}(\nabla\bar{\bfU}+\nabla\bar{\bfU}^{\mathrm{T}}):(\nabla\bar{\bfU}+\nabla\bar{\bfU}^{\mathrm{T}})} \label{lilly} \end{equation} 中的应变率张量数值格式在使用投影法计算时应该为显式吧？
• C

  yPlus在openfoam代码里面的实现
  • Calf.Z-DNS  

  17
  0
  Votes
  17
  Posts
  792
  Views

  C

  @东岳 谢谢老师！这里面的uTau为何直接取1.003呢？不需要植入它的公式吗？