一个2D的算例更好。可以调节进口速度，来测试是否是正常的。
1D算例由于通量守恒，不好测试

估计你下载的不是v11的版本。你要卸载互相适配的版本。

哥们有结果没，我最近也要使用这个求解器

@cccrrryyy 老师，想请教一下，如果对于mixtureKEpsilon这种两相都需要给个湍流模型初始值的该怎么弄呢？有k.air, k.water, km。。。之前估单相的用这个还没啥问题，现在上两相不知道咋整了:135:

屌爆了老铁，聊天记录已收藏，明天就发JFM

constrainPressure主要是更新压力第二类边界条件，公式如下：

\begin{equation}
\left( \nabla p_{rgh} \right)_f \cdot\bfn_f=
\frac{\left(\mathbf{HbyA}_f^{*}- \frac{1}{{{A^n_{\mathrm{P},f}}}}(\bfg\cdot\bfh\nabla\rho)_f - \mathbf{U}_f \right)\cdot\bfS_f}
{
|\bfS_f|
\frac{1}{{{A^n_{\mathrm{P},f}}}}
}
\end{equation}
上述公式与代码并不一致。在OpenFOAM中，constrainPressure为

forAll(pBf, patchi) { if (isA<fixedFluxPressureFvPatchScalarField>(pBf[patchi])) { refCast<fixedFluxPressureFvPatchScalarField> ( pBf[patchi] ).updateCoeffs ( ( phiHbyABf[patchi] - rho.boundaryField()[patchi] *MRF.relative(SfBf[patchi] & UBf[patchi], patchi) ) /(magSfBf[patchi]*rhorAUBf[patchi]) ); } }

多乘了一个密度。应该改为：

forAll(pBf, patchi) { if (isA<fixedFluxPressureFvPatchScalarField>(pBf[patchi])) { refCast<fixedFluxPressureFvPatchScalarField> ( pBf[patchi] ).updateCoeffs ( ( phiHbyABf[patchi] - MRF.relative(SfBf[patchi] & UBf[patchi], patchi) ) /(magSfBf[patchi]*rhorAUBf[patchi]/rho.boundaryField()[patchi]) ); } }

在非常老的OpenFOAM版本中，看起来是正确的，与公式一致

setSnGrad<fixedFluxPressureFvPatchScalarField> ( p_rgh.boundaryField(), ( phiHbyA.boundaryField() - fvOptions.relative(mesh.Sf().boundaryField() & U.boundaryField()) )/(mesh.magSf().boundaryField()*rAUf.boundaryField()) );

@yzwang 直接toposet就行，然后进行积分

@xuqiming 噢噢，这个问题解决了，需要厚度刚好覆盖这一层网格。

可以试试看symmetry边界

@李东岳 老师，我试过了，还是发散的，但是换更高质量的网格也可以收敛，（不收敛的网格icem中的quality为0.6，收敛的为0.75）这种情况是不是应该去尝试找一些新的收敛性更好的算法来实现呢

安装 geometricVofExt，遇到CGAL的链接和指向问题。

问题语句为 include $(GENERAL_RULES)/CGAL

请大佬指点一下，如何解决？！

安装过程如下：

CGAL 安装
安装后，检查是否安装成功： ****************************** Check CGAL ****************************** cgal=true flavour=library root="/home/ABCD/OpenFOAM/ThirdParty-v2506/platforms/linux64Gcc/CGAL-4.14.3" include="/home/ABCD/OpenFOAM/ThirdParty-v2506/platforms/linux64Gcc/CGAL-4.14.3/include" library="/home/ABCD/OpenFOAM/ThirdParty-v2506/platforms/linux64Gcc/CGAL-4.14.3/lib64" ********************************* Done *********************************

其中，ABCD表示 home 目录。

src/CGALVof 目录编译: 使用原始的 options 文件

无法编译 libCGALVof.so

>> wmake libso CGALVof make: *** /home/ABCD/OpenFOAM/OpenFOAM-v2506/wmake/rules/General/CGAL: Is a directory. Stop. make: *** /home/ABCD/OpenFOAM/OpenFOAM-v2506/wmake/rules/General/CGAL: Is a directory. Stop. wmake error: file 'Make/linux64GccDPInt32Opt/sourceFiles' could not be created in /home/ABCD/OpenFOAM/ABCD-v2506/modules/geometricVofExt/src/CGALVof

其中，原始 options 文件为

CGAL_EXACT = /*-DCGAL_DONT_USE_LAZY_KERNEL*/ CGAL_INEXACT = -DCGAL_INEXACT include $(GENERAL_RULES)/CGAL c++CGALWARN = -Wno-old-style-cast EXE_INC = \ -I$(LIB_SRC)/finiteVolume/lnInclude \ -I$(LIB_SRC)/fileFormats/lnInclude \ -I$(LIB_SRC)/surfMesh/lnInclude \ -I$(LIB_SRC)/dynamicFvMesh/lnInclude \ -I$(LIB_SRC)/meshTools/lnInclude \ -I$(LIB_SRC)/surfMesh/lnInclude \ -ICGALVof/lnInclude \ ${CGAL_EXACT} \ ${CGAL_INEXACT} \ ${CGAL_INC} \ ${c++CGALWARN} \ -fopenmp \ -frounding-math EXE_LIBS = \ ${CGAL_LIBS} \ -lfiniteVolume \ -lfileFormats \ -lsurfMesh \ -ldynamicFvMesh \ -lmeshTools

**根据错误提示，定位问题出在语句 include $(GENERAL_RULES)/CGAL **.

src/CGALVof 目录编译: 编译前导入环境变量

注释 options 文件中的语句 /*include $(GENERAL_RULES)/CGAL*/ **.

编译前导入的环境变量为

#!/bin/bash export CGAL_INC=-I${WM_THIRD_PARTY_DIR}/platforms/linux64Gcc/CGAL-4.14.3/include export CGAL_LIBS=-L${WM_THIRD_PARTY_DIR}/platforms/linux64Gcc/CGAL-4.14.3/lib64

编译命令为

source bashrc # bashrc 文件即为导入的环境变量 wclean wmake

成功编译出 libCGALVof.so

applications/utilities/setVofField 编译

根据第 3 步的经验，注释 options 文件中的语句 /*include $(GENERAL_RULES)/CGAL*/ **.

编译结果为

Making dependencies: setVofField.C g++ -std=c++17 -m64 -pthread -DOPENFOAM=2506 -DWM_DP -DWM_LABEL_SIZE=32 -Wall -Wextra -Wold-style-cast -Wnon-virtual-dtor -Wno-unused-parameter -Wno-invalid-offsetof -Wno-attributes -Wno-unknown-pragmas -O3 -DNoRepository -ftemplate-depth-100 -I/home/ABCD/OpenFOAM/OpenFOAM-v2506/src/finiteVolume/lnInclude -I/home/ABCD/OpenFOAM/OpenFOAM-v2506/src/fileFormats/lnInclude -I/home/ABCD/OpenFOAM/OpenFOAM-v2506/src/surfMesh/lnInclude -I/home/ABCD/OpenFOAM/OpenFOAM-v2506/src/dynamicFvMesh/lnInclude -I/home/ABCD/OpenFOAM/OpenFOAM-v2506/src/dynamicMesh/lnInclude -I/home/ABCD/OpenFOAM/OpenFOAM-v2506/src/overset/lnInclude -I/home/ABCD/OpenFOAM/OpenFOAM-v2506/src/meshTools/lnInclude -I/home/ABCD/OpenFOAM/ABCD-v2506/modules/geometricVofExt/src/CGALVof/lnInclude -I/home/ABCD/OpenFOAM/ABCD-v2506/modules/geometricVofExt/src/SimPLIC/lnInclude -I/home/ABCD/OpenFOAM/OpenFOAM-v2506/src/phaseSystemModels/twoPhaseInter/incompressibleInterPhaseTransportModel/lnInclude -I/home/ABCD/OpenFOAM/OpenFOAM-v2506/src/sampling/lnInclude -DCGAL_INEXACT -I/home/ABCD/OpenFOAM/ThirdParty-v2506/platforms/linux64Gcc/CGAL-4.14.3/include -Wno-old-style-cast -fopenmp -frounding-math -iquote. -IlnInclude -I/home/ABCD/OpenFOAM/OpenFOAM-v2506/src/OpenFOAM/lnInclude -I/home/ABCD/OpenFOAM/OpenFOAM-v2506/src/OSspecific/POSIX/lnInclude -fPIC -c setVofField.C -o Make/linux64GccDPInt32Opt/setVofField.o In file included from /usr/include/boost/type_index.hpp:25, from /usr/include/boost/variant/variant.hpp:21, from /usr/include/boost/variant.hpp:17, from /home/ABCD/OpenFOAM/ThirdParty-v2506/platforms/linux64Gcc/CGAL-4.14.3/include/CGAL/Object.h:38, from /home/ABCD/OpenFOAM/ThirdParty-v2506/platforms/linux64Gcc/CGAL-4.14.3/include/CGAL/kernel_basic.h:34, from /home/ABCD/OpenFOAM/ThirdParty-v2506/platforms/linux64Gcc/CGAL-4.14.3/include/CGAL/basic.h:47, from /home/ABCD/OpenFOAM/ThirdParty-v2506/platforms/linux64Gcc/CGAL-4.14.3/include/CGAL/Cartesian/Cartesian_base.h:29, from /home/ABCD/OpenFOAM/ThirdParty-v2506/platforms/linux64Gcc/CGAL-4.14.3/include/CGAL/Simple_cartesian.h:29, from /home/ABCD/OpenFOAM/ThirdParty-v2506/platforms/linux64Gcc/CGAL-4.14.3/include/CGAL/Exact_predicates_inexact_constructions_kernel.h:29, from /home/ABCD/OpenFOAM/ABCD-v2506/modules/geometricVofExt/src/CGALVof/lnInclude/CGALKernel.H:55, from /home/ABCD/OpenFOAM/ABCD-v2506/modules/geometricVofExt/src/CGALVof/lnInclude/CGALSurfaceMesh.H:41, from setVofField.C:36: /usr/include/boost/function_output_iterator.hpp:14:1: note: ‘#pragma message: This header is deprecated. Use <boost/iterator/function_output_iterator.hpp> instead.’ 14 | BOOST_HEADER_DEPRECATED("<boost/iterator/function_output_iterator.hpp>") | ^~~~~~~~~~~~~~~~~~~~~~~ /usr/include/boost/bind.hpp:36:1: note: ‘#pragma message: The practice of declaring the Bind placeholders (_1, _2, ...) in the global namespace is deprecated. Please use <boost/bind/bind.hpp> + using namespace boost::placeholders, or define BOOST_BIND_GLOBAL_PLACEHOLDERS to retain the current behavior.’ 36 | BOOST_PRAGMA_MESSAGE( | ^~~~~~~~~~~~~~~~~~~~ In file included from setVofField.C:43: /home/ABCD/OpenFOAM/ABCD-v2506/modules/geometricVofExt/src/SimPLIC/lnInclude/reconstruction.H:318:13: warning: type qualifiers ignored on function return type [-Wignored-qualifiers] 318 | const bool willMapAlphaField() const | ^~~~~ In file included from setVofField.C:98: setAlphaField.H: In function ‘int main(int, char**)’: setAlphaField.H:32:35: warning: unused variable ‘aridCells’ [-Wunused-variable] 32 | const DynamicList<label>& aridCells(cellLocs.aridCells()); | ^~~~~~~~~ In file included from setVofField.C:136: setAlphaField.H:32:35: warning: unused variable ‘aridCells’ [-Wunused-variable] 32 | const DynamicList<label>& aridCells(cellLocs.aridCells()); | ^~~~~~~~~ g++ -std=c++17 -m64 -pthread -DOPENFOAM=2506 -DWM_DP -DWM_LABEL_SIZE=32 -Wall -Wextra -Wold-style-cast -Wnon-virtual-dtor -Wno-unused-parameter -Wno-invalid-offsetof -Wno-attributes -Wno-unknown-pragmas -O3 -DNoRepository -ftemplate-depth-100 -I/home/ABCD/OpenFOAM/OpenFOAM-v2506/src/finiteVolume/lnInclude -I/home/ABCD/OpenFOAM/OpenFOAM-v2506/src/fileFormats/lnInclude -I/home/ABCD/OpenFOAM/OpenFOAM-v2506/src/surfMesh/lnInclude -I/home/ABCD/OpenFOAM/OpenFOAM-v2506/src/dynamicFvMesh/lnInclude -I/home/ABCD/OpenFOAM/OpenFOAM-v2506/src/dynamicMesh/lnInclude -I/home/ABCD/OpenFOAM/OpenFOAM-v2506/src/overset/lnInclude -I/home/ABCD/OpenFOAM/OpenFOAM-v2506/src/meshTools/lnInclude -I/home/ABCD/OpenFOAM/ABCD-v2506/modules/geometricVofExt/src/CGALVof/lnInclude -I/home/ABCD/OpenFOAM/ABCD-v2506/modules/geometricVofExt/src/SimPLIC/lnInclude -I/home/ABCD/OpenFOAM/OpenFOAM-v2506/src/phaseSystemModels/twoPhaseInter/incompressibleInterPhaseTransportModel/lnInclude -I/home/ABCD/OpenFOAM/OpenFOAM-v2506/src/sampling/lnInclude -DCGAL_INEXACT -I/home/ABCD/OpenFOAM/ThirdParty-v2506/platforms/linux64Gcc/CGAL-4.14.3/include -Wno-old-style-cast -fopenmp -frounding-math -iquote. -IlnInclude -I/home/ABCD/OpenFOAM/OpenFOAM-v2506/src/OpenFOAM/lnInclude -I/home/ABCD/OpenFOAM/OpenFOAM-v2506/src/OSspecific/POSIX/lnInclude -fPIC -Xlinker --add-needed -Xlinker --no-as-needed Make/linux64GccDPInt32Opt/setVofField.o -L/home/ABCD/OpenFOAM/OpenFOAM-v2506/platforms/linux64GccDPInt32Opt/lib \ -L/home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib -L/home/ABCD/OpenFOAM/ThirdParty-v2506/platforms/linux64Gcc/CGAL-4.14.3/lib64 -lfiniteVolume -lfileFormats -lsurfMesh -ldynamicFvMesh -loverset -lmeshTools -lwaveModels -lCGALVof -lSimPLIC -limmiscibleIncompressibleTwoPhaseMixture -lsampling -lOpenFOAM -ldl \ -lm -o /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/bin/setVofField /usr/bin/ld: Make/linux64GccDPInt32Opt/setVofField.o: in function `CGAL::Interval_nt<false>::Test_runtime_rounding_modes::Test_runtime_rounding_modes()': setVofField.C:(.text._ZN4CGAL11Interval_ntILb0EE27Test_runtime_rounding_modesC2Ev[_ZN4CGAL11Interval_ntILb0EE27Test_runtime_rounding_modesC5Ev]+0xc9): undefined reference to `CGAL::assertion_fail(char const*, char const*, int, char const*)' /usr/bin/ld: setVofField.C:(.text._ZN4CGAL11Interval_ntILb0EE27Test_runtime_rounding_modesC2Ev[_ZN4CGAL11Interval_ntILb0EE27Test_runtime_rounding_modesC5Ev]+0xe8): undefined reference to `CGAL::assertion_fail(char const*, char const*, int, char const*)' /usr/bin/ld: Make/linux64GccDPInt32Opt/setVofField.o: in function `CGAL::Interval_nt<true>::Test_runtime_rounding_modes::Test_runtime_rounding_modes()': setVofField.C:(.text._ZN4CGAL11Interval_ntILb1EE27Test_runtime_rounding_modesC2Ev[_ZN4CGAL11Interval_ntILb1EE27Test_runtime_rounding_modesC5Ev]+0xc9): undefined reference to `CGAL::assertion_fail(char const*, char const*, int, char const*)' /usr/bin/ld: setVofField.C:(.text._ZN4CGAL11Interval_ntILb1EE27Test_runtime_rounding_modesC2Ev[_ZN4CGAL11Interval_ntILb1EE27Test_runtime_rounding_modesC5Ev]+0xe8): undefined reference to `CGAL::assertion_fail(char const*, char const*, int, char const*)' /usr/bin/ld: Make/linux64GccDPInt32Opt/setVofField.o: in function `_GLOBAL__sub_I_setVofField.C': setVofField.C:(.text.startup+0x2e6): undefined reference to `CGAL::Random::Random()' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `mpfr_get_emin' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `__gmpn_com' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `__gmpq_swap' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `__gmpq_init' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `CGAL::precondition_fail(char const*, char const*, int, char const*)' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `__gmpn_copyi' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `__gmpq_cmp' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `mpfr_set_q' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `CGAL::postcondition_fail(char const*, char const*, int, char const*)' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `__gmpn_add_n' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `__gmpq_set_d' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `mpfr_get_d' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `__gmpq_sub' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `__gmpz_init_set_ui' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `__gmpq_clear' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `__gmpn_sqr' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `mpfr_set_emin' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `__gmpn_sub_n' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `__gmpq_mul' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `__gmpq_div' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `CGAL::Random::Random(unsigned int)' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `__gmpz_init_set_si' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `mpfr_subnormalize' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `__gmpq_set' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `__gmpn_mul' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `CGAL::warning_fail(char const*, char const*, int, char const*)' /usr/bin/ld: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/lib/libCGALVof.so: undefined reference to `__gmpq_add' collect2: error: ld returned 1 exit status make: *** [/home/ABCD/OpenFOAM/OpenFOAM-v2506/wmake/makefiles/general:182: /home/ABCD/OpenFOAM/ABCD-v2506/platforms/linux64GccDPInt32Opt/bin/setVofField] Error 1

setShapedFields 的编译错误与上面的类似。

@studentTan nCellsBetweenLevels 过渡层设置多一点，试下10层或20层

nCellsBetweenLevels 1;

画成一个平的圆形就行

这问题老大难了，基本要一事一议

有什么好的方法展示粒子的运动轨迹吗？

大家好
我想在SOWFA模拟过程中让风力机偏航角度随时间变化（动态偏航），请问有什么方法可以实现吗？我看SOWFA里面的偏航角设置是一个定值。

当时可真能写...

哎呀我勒个去 当时自己可真能写 现在看了都想删了 不嫌害臊