@李东岳 感谢老师回复！动量预测是关掉了的。

打开试了试，步长极小，残差也一直保持在一个值不下降，cumulative则很大，如下：

Courant Number mean: 7.953414508e-05 max: 0.4993036972 deltaT = 4.979675437e-35 Time = 6.00001e-05 PIMPLE: iteration 1 DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 0.04941709914, No Iterations 29 DICPCG: Solving for pcorr, Initial residual = 0.03122316034, Final residual = 0.001556737496, No Iterations 121 DICPCG: Solving for pcorr, Initial residual = 0.006804280077, Final residual = 0.0003314734696, No Iterations 359 DICPCG: Solving for pcorr, Initial residual = 0.002308369425, Final residual = 9.709852492e-11, No Iterations 911 time step continuity errors : sum local = 1.793718362e-20, global = 1.178518414e-23, cumulative = 1.037419916e+19 smoothSolver: Solving for Ux, Initial residual = 3.985498772e-05, Final residual = 2.227155701e-11, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 6.302253167e-05, Final residual = 2.385619185e-11, No Iterations 1 GAMG: Solving for p, Initial residual = 0.004038713667, Final residual = 0.0001722133227, No Iterations 18 GAMG: Solving for p, Initial residual = 0.0005583071565, Final residual = 2.673668049e-05, No Iterations 34 GAMG: Solving for p, Initial residual = 9.485975434e-05, Final residual = 4.53317788e-06, No Iterations 33 GAMG: Solving for p, Initial residual = 2.596708738e-05, Final residual = 1.231466526e-06, No Iterations 38 time step continuity errors : sum local = 7.929103562e-15, global = -8.302368146e-16, cumulative = 1.037419916e+19 GAMG: Solving for p, Initial residual = 0.003717218882, Final residual = 0.000164557501, No Iterations 17 GAMG: Solving for p, Initial residual = 0.0005378364039, Final residual = 2.53035408e-05, No Iterations 34 GAMG: Solving for p, Initial residual = 9.265726415e-05, Final residual = 4.535601234e-06, No Iterations 36 GAMG: Solving for p, Initial residual = 2.556961887e-05, Final residual = 9.918214454e-11, No Iterations 548 time step continuity errors : sum local = 6.482264028e-19, global = 6.603865839e-20, cumulative = 1.037419916e+19 PIMPLE: iteration 2 ... PIMPLE: iteration 50 smoothSolver: Solving for Ux, Initial residual = 1.226945138e-10, Final residual = 6.6691269e-16, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 1.651845516e-10, Final residual = 3.177048458e-16, No Iterations 1 GAMG: Solving for p, Initial residual = 0.000949721577, Final residual = 4.516299707e-05, No Iterations 41 GAMG: Solving for p, Initial residual = 0.0004457930337, Final residual = 2.206483819e-05, No Iterations 23 GAMG: Solving for p, Initial residual = 7.977007558e-05, Final residual = 3.747481321e-06, No Iterations 36 GAMG: Solving for p, Initial residual = 2.418915618e-05, Final residual = 1.162060008e-06, No Iterations 48 time step continuity errors : sum local = 7.482194647e-15, global = -7.221446192e-16, cumulative = 1.037419916e+19 GAMG: Solving for p, Initial residual = 0.003705959955, Final residual = 0.0001646823634, No Iterations 17 GAMG: Solving for p, Initial residual = 0.0005377574684, Final residual = 2.535476181e-05, No Iterations 34 GAMG: Solving for p, Initial residual = 9.241219827e-05, Final residual = 4.525895839e-06, No Iterations 36 GAMG: Solving for p, Initial residual = 2.542969083e-05, Final residual = 9.774400044e-11, No Iterations 548 time step continuity errors : sum local = 6.388632178e-19, global = 6.502735317e-20, cumulative = 1.037419916e+19 PIMPLE: not converged within 50 iterations ExecutionTime = 21923.2 s ClockTime = 21963 s Courant Number mean: 7.964991508e-05 max: 0.4989139005 deltaT = 4.99051583e-35 Time = 6.00001e-05

看到站内其它的帖子，试了试有量纲/无量纲计算（此前圆柱直径0.05，nu=50；改为圆柱直径0.0005，nu=0.5），也无法收敛。

走投无路试了一下StarCCM，多低雷诺数都能算，和OpenFOAM那些能收敛的算例比得还贼好

各位好，我现在设想一个工况是在水深非常浅，同时水深与船舶吃水非常小（大约1.2）的情况，在设置重叠部分也即船舶部分网格时，blockmesh需要画的比较小，底部与船舶底部很接近，这就造成了snappyhexmesh后船舶下部的网格质量很差，请问有没有什么办法或者思路来改善网格质量呢，如有相似问题的朋友欢迎讨论，谢谢!261e3ad06976a58e400a1fc15243dc1.png 78823dc1464db8ab58735e54dcd17d2.png efa88c4c50ed4254badb8dbc8db3089.png

@学流体的小明 好的，非常感谢解答

我也在算类似的问题，是不是给定进口水位和流量的设置方法不对

用interFoam计算带自由液面的倒虹吸，参考其他人边界设置，往往是入口给定流量（流速），出口给定水位（alpha.air)。这种设置方式如果通过手动设置（setFields)往往因初始流场不好极易发散，有没有办法能跟Fluent一样通过混合初始化得到一个初始场保证计算能够进行下去。（openFoam好像可以通过potentialFoam初始化场，但是这个对于interFoam也可以使用吗）79538e5a-b907-420a-8777-c7a72dd6011f-image.png

@cccrrryyy谢谢前辈

你好，请问你编写的NSRFG，在使用时计算域是多大，网格尺度是多大，能否在整个计算域上保持低耗散甚至无耗散么

参考之前大神发的帖子就可以解决了
https://www.cfd-china.com/topic/6750/mapfield报错

@jxzs_zy 在开发，好像还没正式发布。Added dynamicPolyRefinementFvMesh

问题找到了，对于这个算例，在滤波尺度的选择中，不能够使用剪切层调整亚格子长度尺度修正，但是别的文献中有人使用过，目前没有找到原因。

你好，请问这个问题解决了吗？

您好，您解决了吗

您好，您解决了吗

@BznW 您好。打扰一下，您解决了吗

前一段时间搭建了流化床求解器（燃烧、密相颗粒流动求解器，下称newSolver），目前需要对这个求解器进行冷态颗粒流动的验证（常温，关闭化学反应、燃烧和传热）。通过验证发现颗粒轴向速度和实验吻合较好，但是气相压降差距很大，而使用MPPICFoam计算相同的case可以得到准确的压降（稳定时newSolver比MPPICFoam压降高两倍多）。

为了找到问题所在，于是关闭颗粒计算纯气体流动的速度和压降，发现newSolver和MPPICFoam计算得到的气相速度相同但是压降有一定差距。因此猜测是气相压力方程计算的问题。

由于newSolver是基于coalChemistryFoam流体方程加入alpha的修改的，纯气体流动情况下alpha=1。为了排除我修改方程所导致的错误，因此用coalChemistryFoam（关闭化学反应、燃烧和传热）计算纯气体流动看能否得到准确的结果。结果发现压降仍然有差距。

难道coalChemisryFoam变密度的求解器无法计算常密度的情况吗？coalChemistryFoam是反应所导致的密度变化，关闭化学反应按理说密度应该是不变的，为何计算结果与实际差距大？常温纯气体流动压降都不准，加上化学反应燃烧后压降会准吗？

MPPICFoam压力方程

fvScalarMatrix pEqn ( fvm::laplacian(alphacf*rAUcf, p) == fvc::ddt(alphac) + fvc::div(alphacf*phiHbyA) );

coalChemistryFoam压力方程

fvScalarMatrix pEqn ( fvm::ddt(psi, p) + fvc::div(phiHbyA) - fvm::laplacian(rhorAUf, p) == coalParcels.Srho() + fvOptions(psi, p, rho.name()) );

纯气体流动验证的几何模型和网格如下图所示
图片1.png

请各位大佬指教！

这段用来统计网格单元内粒子个数的代码，编译的路径是什么，求东岳老师指导

基于interFoam求解器，添加了标量场S，在模拟了一段时间后，得到了两相流体分布，再使用swak4foam中的funkySetFields，把标量赋值到alpha1流体中，另一相不考虑，设置如下：
b8c93503-efa8-4a6d-9b18-7852dcb2d072-image.png
运行后报错:
15a01aae-36ee-42e4-8223-4fc9eba68dd9-image.png
请问各位大神该如何解决呢？

@liuyf 我在求解时候也遇到了这个问题，导致鼓泡塔的液面比实验高很多，因为上边部分存在液体速度导致求解一开始液体会一下子冲上去，最后液面高度就会变得非物理，您解决这个问题了吗

g++ -std=c++11 -m64 -Dlinux64 -DWM_ARCH_OPTION=64 -DWM_DP -DWM_LABEL_SIZE=32 -Wall -Wextra -Wold-style-cast -Wnon-virtual-dtor -Wno-unused-parameter -Wno-invalid-offsetof -Wno-attributes -O3 -DNoRepository -ftemplate-depth-100 -I/usr/local/include -I./windProfile/lnInclude -I./extrudedPatch/lnInclude -I/home/windlab_lzx/openfoam8/OpenFOAM-8/src/finiteVolume/lnInclude -I/home/windlab_lzx/openfoam8/OpenFOAM-8/src/meshTools/lnInclude -I/home/windlab_lzx/openfoam8/OpenFOAM-8/src/sampling/lnInclude -I/home/windlab_lzx/openfoam8/OpenFOAM-8/src/surfMesh/lnInclude -I/opt/intel/oneapi/mkl/2024.1/include -IlnInclude -I. -I/home/windlab_lzx/openfoam8/OpenFOAM-8/src/OpenFOAM/lnInclude -I/home/windlab_lzx/openfoam8/OpenFOAM-8/src/OSspecific/POSIX/lnInclude -fPIC -fuse-ld=bfd -Xlinker --add-needed -Xlinker --no-as-needed /home/windlab_lzx/openfoam8/OpenFOAM-8/platforms/linux64GccDPInt32Opt/applications/utilities/preProcessing/DFSRTurb/DFSRTurb.o -L/home/windlab_lzx/openfoam8/OpenFOAM-8/platforms/linux64GccDPInt32Opt/lib \ -L/usr/local/include -lwindProfile -lextrudedPatch -lmeshTools -lsampling -lsurfMesh -lfiniteVolume -lfftw3 -L/opt/intel/oneapi/mkl/2024.1/lib/intel64 -lmkl_intel_lp64 -L/opt/intel/oneapi/mkl/2024.1/lib/intel64 -lmkl_core -L/opt/intel/oneapi/mkl/2024.1/lib/intel64 -lmkl_sequential -lm -ldl -lOpenFOAM -ldl \ -lm -o /home/windlab_lzx/OpenFOAM/windlab_lzx-8/platforms/linux64GccDPInt32Opt/bin/DFSRTurb /usr/bin/ld.bfd: 找不到 -lwindProfile /usr/bin/ld.bfd: 找不到 -lextrudedPatch collect2: error: ld returned 1 exit status /home/windlab_lzx/openfoam8/OpenFOAM-8/wmake/makefiles/general:140: recipe for target '/home/windlab_lzx/OpenFOAM/windlab_lzx-8/platforms/linux64GccDPInt32Opt/bin/DFSRTurb' failed make: *** [/home/windlab_lzx/OpenFOAM/windlab_lzx-8/platforms/lDFSRinux64GccDPInt32Opt/bin/DFSRTurb] Error 1在编译DFSR方法最后一步wmake时候报错，为什么呀，找了很久没解决

@shirouchen 谢谢，我暂时是用matlab排列的

大家好！想请教各位在tecplot中对三维模型进行Mirror操作后，使用slices提取出一个面，可以观察到提取出来的面云图都还是对称的,但画出来的流线为什么不是对称的了？原计算案例对称轴所在的面本身使用的也是symmetryPlane型边界条件。
QQ图片20240407133422.png
$DR1%SP2~_Y1_A{0ZO8ZND8.png
TUG[A$3E@2}8%R85`77T]3K.png

那我不知道咋办了，要不你用apt删了重装一下，或者编译源代码安装。
或者你电脑提前改变了apt的安装路径，安装到了其它地方？
下面这个里面安装的openfoam9就是在/opt目录下面。
https://blog.csdn.net/sagjhdj/article/details/123435344

@学流体的小明
1、谢谢大佬！就是这个原因，mesh.Cf()改成mesh.C()就没问题了！
2、当时注释这个是因为我发现有问题后用Info输出坐标和mesh的size想看看问题出在哪，但是只输出了processer0的网格信息，所以有了这个疑问。后来我查到了是Info只会输出master processor的信息，我改成Pout就可以输出所有processor的信息。

最近在做一个模型，有一个面既按照一定的角速度旋转，同时又有流体射出，但是不知道该怎么设置这个面的边界条件，有大佬知道吗

@尚善若水 问题解决了，多谢多谢

已解决。

这么做尖角网格是不是看起来不太对呀，这类尖角网格有没有什么好一点的做法呢？
虽然看着网格质量还行，但总感觉不太对
Snipaste_2024-04-03_21-48-45.png

在轴对称模型中，可以只计算两个方向（轴向和径向）的量吗，在写方程时，如何不计算第三个方向的偏导数，可以像volTensorField一样定义一个Tensor2D的体积场吗

@coolhhh 太细致了！！！

@huangyuhui723 我明白你的意思了。

@huangyuhui723 在 关于在求解器中使用热物理库中的函数 中说：

通过其它计算更新了T之后想对he赋值

如果是求解TEqn，可以是得到温度T后更新he。OpenFOAM里的求解器一般都是求解EEqn，得到he后通过thermo.correct()更新T。

如果已知边界温度T后要对he赋值，应该在求解EEqn之前对he边界赋值，然后求解EEqn。大致代码是这样：

volScalarField& he = thermo.he(); //****************************************//加这一段 forAll(p.boundaryField(), patchi) { forAll(p.boundaryField()[patchi], facei) { he.boundaryFieldRef()[patchi][facei] = function(p.boundaryField()[patchi][facei], T.boundaryField()[patchi][facei]); } } //****************************************// fvScalarMatrix EEqn ( ...... );

这个function函数就是T与he的关系，类似上面，混合物的he通过各组分的显焓的质量分数加权得到，显焓用janaf里的关系确定。

供参考。

@学流体的小明 感谢老师，我找时间试试看

我在做单根不锈钢丝切割液滴的模拟，发现液滴触碰到固体后，会快速铺展覆盖固体，这明显与实际情况不符：
dc6202db-0a4b-4bee-8dfd-d8db7a8ba44e-image.png
9f370870-b0be-4f65-9075-57dd718ce5cd-image.png
求解器用的v2212的MPPICInterFOAM,固体壁面的边界条件为：U noSlip；P fixedFluxPressure；Alpha constantAlphaContactAngle，接触角为155疏水界面

文献和实际可视化拍摄中，液体不会有如此快速的铺展，还请各位指点一二
c833896f-62e9-4161-ae81-1b226e382477-image.png

@尚善若水 感谢，我去试试

1711709658497.png阿里云的VFC链接可以安装图形界面，按照这个教程来就可以，非常简单https://help.aliyun.com/zh/ecs/user-guide/install-a-gui-desktop-for-ubuntu-18

可以从边界上的face返回这个face属于的cell

const Foam::fvBoundaryMesh &procBoundary = mesh.boundary(); forAll(procBoundary, patch) { forAll(procBoundary[patch], facei) { const label bcell = boundaryMesh[patch].faceCells()[facei]; } }

你可以得到所有的边界单元编号，然后做你其它向做的事就应该比较方便了

@pengdi 没有收缩段，RANS结果就是很差的，LES即使入口只加平均风，也比RANS结果精确很多

已经搞清楚了，谢谢！

openfoam自带有interfaceHeight的function

报错的自定义边界如下，如图所示报错为：“‘incompressible’ does not name a type“，恳请论坛内各位老师解惑
/---------------------------------------------------------------------------\

\ / F ield OpenFOAM: The Open Source CFD Toolbox \ / O peration \ / A nd Copyright (C) 2011 OpenFOAM Foundation \/ M anipulation

License
This file is part of OpenFOAM.

OpenFOAM is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version. OpenFOAM is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.

*---------------------------------------------------------------------------*/

#include "adjointOutletPressureHeatFvPatchScalarField.H"
#include "addToRunTimeSelectionTable.H"
#include "fvPatchMapper.H"
#include "volFields.H"
#include "surfaceFields.H"
#include "RASModel.H"

// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //

Foam::adjointOutletPressureHeatFvPatchScalarField::
adjointOutletPressureHeatFvPatchScalarField
(
const fvPatch& p,
const DimensionedField<scalar, volMesh>& iF
)
:
fixedValueFvPatchScalarField(p, iF)
{}

Foam::adjointOutletPressureHeatFvPatchScalarField::
adjointOutletPressureHeatFvPatchScalarField
(
const adjointOutletPressureHeatFvPatchScalarField& ptf,
const fvPatch& p,
const DimensionedField<scalar, volMesh>& iF,
const fvPatchFieldMapper& mapper
)
:
fixedValueFvPatchScalarField(ptf, p, iF, mapper)
{}

Foam::adjointOutletPressureHeatFvPatchScalarField::
adjointOutletPressureHeatFvPatchScalarField
(
const fvPatch& p,
const DimensionedField<scalar, volMesh>& iF,
const dictionary& dict
)
:
fixedValueFvPatchScalarField(p, iF)
{
fvPatchField<scalar>::operator=
(
scalarField("value", dict, p.size())
);
}

Foam::adjointOutletPressureHeatFvPatchScalarField::
adjointOutletPressureHeatFvPatchScalarField
(
const adjointOutletPressureHeatFvPatchScalarField& tppsf,
const DimensionedField<scalar, volMesh>& iF
)
:
fixedValueFvPatchScalarField(tppsf, iF)
{}

// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //

void Foam::adjointOutletPressureHeatFvPatchScalarField::updateCoeffs()
{
if (updated())
{
return;
}

const fvsPatchField<scalar>& phip = patch().lookupPatchField<surfaceScalarField, scalar>("phi"); const fvsPatchField<scalar>& phiap = patch().lookupPatchField<surfaceScalarField, scalar>("phib"); const fvPatchField<vector>& Uap = patch().lookupPatchField<volVectorField, vector>("Ub"); const dictionary& transportProperties = db().lookupObject<IOdictionary>("transportProperties"); dimensionedScalar nu(transportProperties.lookup("nu")); scalarField Up_n = phip / patch().magSf();//Primal scalarField Uap_n = phiap / patch().magSf();//Adjoint const incompressible::RASModel& rasModel = db().lookupObject<Foam::incompressible::RASModel>("RASProperties"); scalarField nueff = rasModel.nuEff()().boundaryField()[patch().index()]; const scalarField& deltainv = patch().deltaCoeffs(); // distance^(-1) scalarField Uaneigh_n = (Uap.patchInternalField() & patch().nf()); //operator ==((Up_n * Uap_n) +2*nu.value()*deltainv*(Uap_n-Uaneigh_n));

operator ==((Up_n * Uap_n) +2*nueff.value()deltainv(Uap_n-Uaneigh_n));

fixedValueFvPatchScalarField::updateCoeffs();

}

void Foam::adjointOutletPressureHeatFvPatchScalarField::write(Ostream& os) const
{
fvPatchScalarField::write(os);
writeEntry("value", os);
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

namespace Foam
{
makePatchTypeField
(
fvPatchScalarField,
adjointOutletPressureHeatFvPatchScalarField
);
incompressible
}
f0c155d4-8455-4c68-be20-3db9aa49d6b6-image.png

@疏影横斜水清浅 在 topoSet, createPatch, createBaffle 在网格内部添加一个平面 求指教 中说：

@脉动通风 你好，能详细说一下么，你是用的哪种颗粒注入方式，内部画网格是指生成面还是cellzone，你使用的方式是连续注入么

内部画网格就是我自己在房间中间挖了一块正方体的区域，然后再使用ICEM画的网格。最后传入Openfoam中，再对那个正方体的一个面设置气体输入的浓度。
使用的方式是连续注入，我是注入气体。

各位老师，大佬晚上好！openfoam小白最近在阅读icoFoam源代码时遇到了一点问题，根据李老师对icoFoam求解器的理论推导，它的压力泊松方程如下：7.png
99.png
在源代码中首先定义了phiHbyA,应该是HbyA的面通量:
1.png
但是,它的压力泊松方程里面的用的是phiHbyA而不是HbyA,这里我不是很能理解：
3.png
另外按理说piso算法每算一次U、P，按理来说应该要更新HbyA一次：
8.png
但是源代码里面更新的好像是phiHbyA（速度变了，且是个通量）
在我的理解里，这段代码每次进行piso算法时，HbyA似乎没变，毕竟它是从之前的速度方程那里得到的？（对UEqn.A这个函数不太了解）
9.png
但是piso算完以后的U用的又是和理论推导的结果一致，用的是HbyA，而不是phiHbyA:
6.png
恳请各位老师，大佬解惑，谢谢各位！

@tidedrinker 前辈你好，我在使用olaFlow做一个二维的狭缝共振，但是把两个浮箱加上之后造出的波浪不向前传递是什么原因呀，详细的我放在了这个问题里面问题地址

参考openfoam自带的dambreach这个算例

@wangfei9088 好滴，感谢大佬:xinxin2:

@学流体的小明 嗯嗯，谢谢大佬

eb812a45-a24b-435b-a01c-246db6db8afc-250c0b40b51035fb0e43108730efee4.jpg
请教各位：
我想如图设置一个扇形的计算域，中间一段圆环设置为旋转动网格，其中涉及a-a' 与b-b'两个对应面之间设置的问题，因为整个过程中，两组对应面并不是时刻交互的，某些面区域会作为边界，这种情况应该如何设置呢？谢谢

各位前辈好，我请教一下一个问题，我在使用olaFoam来建立一个模型，是一个二维的双浮体狭缝共振，我单独造波的时候（波高0.024，周期2s）的时候是可以造波的
a1150344-d9bd-4c60-a03a-df521fec528d-image.png
但是当我把两个固定的浮体加上之后，波浪就几乎没有了，请问这是什么原因

@学流体的小明 我感觉也是时间太长，我目前打算检测一个点，观察温度的峰值和峰谷，取这个时间段作为周期

@wangfei9088 非常感谢大佬的指点！！ 完美解决了这个问题:140:

@Do1975 确实，可能还是看具体需求场景，我这是一方面DBS方法速度提升好像不多，然后在AB区域交界附近流场会出现间断，我没解决这个问题，就没用了