@小李的幸福生活 抱歉才看到消息。还需要吗？需要的话你留个邮箱吧，我发你。
H
Hope
@Hope
帖子
-
中性大气环境湍流动能的自保持 | 附有算例下载 -
cyclicAMI：旋转周期边界设置问题@李东岳 谢谢李老师，网格检查过了，没有问题，把cyclicAMI的边界设置为symmetry是可以计算的。刚才按照您的建议试过了，也还是报相同的错误，错误信息如下：
/*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2112 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : _6e1fca0e-20220610 OPENFOAM=2112 patch=220610 version=2112 Arch : "LSB;label=32;scalar=64" Exec : dbnsFoam Date : Jun 06 2024 Time : 18:13:40 Host : DESKTOP-H49IU2O PID : 137843 I/O : uncollated Case : /home/duan/JIA_Xiaotong/ICSFoam-master/tutorials/test_5 nProcs : 1 trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20) allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Steady-state analysis detected pseudoTime: tolerance = 1e-10, relTol = 1e-10 Local timestepping selected Reading thermophysical properties Selecting thermodynamics package { type hePsiThermo; mixture pureMixture; transport const; thermo hConst; equationOfState perfectGas; specie specie; energy sensibleInternalEnergy; } AMI: Creating addressing and weights between 18620 source faces and 18620 target faces AMI: Patch source sum(weights) min:0 max:1.84903 average:0.999732 AMI: Patch target sum(weights) min:0.030709 max:1.15763 average:0.994612 #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in /lib/x86_64-linux-gnu/libc.so.6 #3 ? at ??:? #4 ? at ??:? #5 Foam::psiThermo::addfvMeshConstructorToTable<Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > > >::New(Foam::fvMesh const&, Foam::word const&) at ??:? #6 ? at ??:? #7 Foam::psiThermo::New(Foam::fvMesh const&, Foam::word const&) at ??:? #8 ? in ~/OpenFOAM/duan-v2112/platforms/linux64GccDPInt32Opt/bin/dbnsFoam #9 ? in /lib/x86_64-linux-gnu/libc.so.6 #10 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6 #11 ? in ~/OpenFOAM/duan-v2112/platforms/linux64GccDPInt32Opt/bin/dbnsFoam Floating point exception老师，有没有可能是因为叶栅模型的旋转半径太大（120m)导致的呢？
-
cyclicAMI：旋转周期边界设置问题我在在设置旋转周期边界时遇到了以下问题：
我使用的是一个openfoam外开发的求解器，求解C3X叶栅的模拟，在使用相同的边界条件时，当设置叶片两个周期性边界为对称边界时，求解一切正常，但将其改为cyclicAMI时则计算无法进行，个人角度感觉是边界条件设置的问题，尝试过在·boundary文件中设置separationVector等，均会出下以下错误：// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Steady-state analysis detected pseudoTime: tolerance = 1e-10, relTol = 1e-10 Local timestepping selected Reading thermophysical properties Selecting thermodynamics package { type hePsiThermo; mixture pureMixture; transport const; thermo hConst; equationOfState perfectGas; specie specie; energy sensibleInternalEnergy; } AMI: Creating addressing and weights between 18620 source faces and 18620 target faces AMI: Patch source sum(weights) min:0 max:1.84903 average:0.999732 AMI: Patch target sum(weights) min:0.030709 max:1.15763 average:0.994612 #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in /lib/x86_64-linux-gnu/libc.so.6 #3 ? at ??:? #4 ? at ??:? #5 Foam::psiThermo::addfvMeshConstructorToTable<Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > > >::New(Foam::fvMesh const&, Foam::word const&) at ??:? #6 ? at ??:? #7 Foam::psiThermo::New(Foam::fvMesh const&, Foam::word const&) at ??:? #8 ? in ~/OpenFOAM/duan-v2112/platforms/linux64GccDPInt32Opt/bin/dbnsFoam #9 ? in /lib/x86_64-linux-gnu/libc.so.6 #10 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6 #11 ? in ~/OpenFOAM/duan-v2112/platforms/linux64GccDPInt32Opt/bin/dbnsFoam Floating point exception网格的质量没有问题，boundary处设置的条件如下：
( blade { type wall; inGroups 1(wall); nFaces 25088; startFace 6512669; } inlet { type patch; nFaces 8526; startFace 6537757; } outlet { type patch; nFaces 8526; startFace 6546283; } Lower_periodicity { type cyclicAMI; inGroups 1(cyclicAMI); nFaces 18620; startFace 6554809; matchTolerance 0.0001; transform rotational; neighbourPatch Upper_periodicity; rotationAxis (1 0 0); rotationCentre (0 0 0); } Upper_periodicity { type cyclicAMI; inGroups 1(cyclicAMI); nFaces 18620; startFace 6573429; matchTolerance 0.0001; transform rotational; neighbourPatch Lower_periodicity; rotationAxis (1 0 0); rotationCentre (0 0 0); } back { type wall; inGroups 1(wall); nFaces 22363; startFace 6592049; } front { type wall; inGroups 1(wall); nFaces 22363; startFace 6614412; } )同时也尝试过将cyclicAMI改为cyclic，但是两个边界的网格差异过大，有以下问题：
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 --> FOAM FATAL ERROR: (openfoam-2112 patch=220610) face 0 area does not match neighbour by 68.9656% -- possible face ordering problem. patch:Lower_periodicity my area:6.66444e-09 neighbour area:1.36797e-08 matching tolerance:0.1 Mesh face:6554809 fc:(-0.0309194 -0.0589349 119.957) Neighbour fc:(-0.0196249 0.058833 119.957) If you are certain your matching is correct you can increase the 'matchTolerance' setting in the patch dictionary in the boundary file. Rerun with cyclic debug flag set for more information. From void Foam::cyclicPolyPatch::calcTransforms(const primitivePatch&, const pointField&, const vectorField&, const pointField&, const vectorField&) in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 215. FOAM exiting望大佬帮忙解答下困惑。
-
中性大气环境湍流动能的自保持 | 附有算例下载@李东岳 李老师，yang2009这篇sci论文里面我没看到有壁面剪应力的相关设置。由于这部分内容是我师兄做的，我先依据硕士论文解释您提出的几个问题，具体的我还需要再查证之后再给您回复。
硕士论文中引用的入流条件也是出自这篇sci文献，有出入的是壁面函数的粗糙高度取值，硕士论文中的粗糙高度是依据文献（Blocken et al. 2007, CFD simulation of the atmospheric boundary layer: wall function problems）计算所得：
壁面剪应力的取值是依据硕士论文中的一些假定提出的，而非yang2009中相关参数。 -
中性大气环境湍流动能的自保持 | 附有算例下载@李东岳 李老师您好，我们这边之前在fluent里面有验证过该入流条件，具体结果可参见刘安彬的硕士论文《平衡大气边界层的CFD数值模拟》。但是在模拟时上部流场的湍动能并无明显衰减，只是近壁面有衰减，并在施加壁面剪应力之后可有效改善。您可以考虑下顶面边界条件、壁面剪应力这两方面的设置。
下图模拟结果截取自上述硕士论文，未加任何修正
下图是施加壁面剪应力后的模拟结果
-
算例比较upwind与linearUpwind对结果的影响你是否比较Upwind格式下的残差曲线和linearUpwind 格式的差别呢？这两方法在OpenFOAM中的收敛性是不一样的。或者设置Porbes探针观察下数据采样点最后的波动情况。我觉的是收敛性影响的。有复杂地形或者建筑群的边界层计算时，OpenFOAM高阶阶方案的收敛性并不好，会影响计算精度@李东岳
-
请教各位老师关于动量控制方程在数值求解过程中的问题请教各位老师两个问题（1）在采用有限体积法对方程进行离散之后，动量控制方程中左侧的惯性力项如何计算？
（2）在高湍流度的流动中，方程右侧的粘性力项很小，流体的受力状态如何理解？（左侧惯性力与右侧压力梯度的平衡吗？广义源项的加入如何影响流动的发展呢？） 谢谢~~~
-
边界条件tableFile读取错误问题@东岳 好的，我试一下。谢谢老师~
-
边界条件tableFile读取错误问题@东岳 嗯嗯，是在同一个文件夹0/下的
-
边界条件tableFile读取错误问题@东岳 李老师，我边界条件代码修改后如下：
INLET { type fixedProfile; profile tableFile; profileCoeffs { nHeaderLine 0; // Number of header lines refColumn 0; // Reference column index componentColumns (1 2 3); // Component column indices separator ","; // Optional (defaults to ",") mergeSeparators no; // Merge multiple separators file "inVelocity.txt";//"inletProfiles/inVelocity.csv"; outOfBounds clamp; // Optional out-of-bounds handling interpolationScheme linear; // Optional interpolation scheme } direction (0 1 0); origin 0; }同时inVelocity.txt数据按照下面格式给出，但还是有同样的报错信息
2.31474,0,0 3.07832,0,0 3.44144,0,0 3.64768,0,0 3.78464,0,0 3.88646,0,0 3.96812,0,0 4.03721,0,0 4.09802,0,0 4.15319,0,0 4.20449,0,0 4.25314,0,0 4.30008,0,0 4.34601,0,0 4.39149,0,0 4.43702,0,0 4.483,0,0 4.5298,0,0 -
边界条件tableFile读取错误问题模拟二维大气边界层流动，入口添加ESDU的对数律剖面时，采用如下方式定义，
INLET { type fixedProfile; profile tableFile; profileCoeffs { nHeaderLine 0; // Number of header lines refColumn 1; // Reference column index componentColumns (1); // Component column indices separator " "; // Optional (defaults to ",") mergeSeparators no; // Merge multiple separators file "inVelocity.txt";//"inletProfiles/inVelocity.csv"; outOfBounds clamp; // Optional out-of-bounds handling interpolationScheme linear; // Optional interpolation scheme } direction (0 1 0); origin 0; }运行
simpleFoam时，提示如下错误
inVelocity.txt文件为0.899251 3.78464 3.88646 3.96812 4.03721 4.09802 4.15319 4.20449 4.25314 4.30008 4.34601 4.39149 4.43702 4.483 4.5298 4.57773 4.62706 4.67807 4.73099 4.78608 4.84357 4.90368 4.96665 5.03268 5.10199 5.17476 5.25119 5.33143 5.41563 5.50394 5.59645 5.69327 5.79447 5.90012 6.01023不知道问题在哪里？还请各位指教