@李东岳
老师好,用的interFoam和twoLiquidMixingFoam,是VOF?
merlinizlf
帖子
-
-
想设置一组边界条件,就是两相流压力驱动流,我设置的边界条件总是导致发散,想请教一下有没有稳定的边界条件?
这个边界条件算单向流的时候好像没有问题,怎么到两相流就不行了呢?Inlet: U: type zeroGradient; p: type fixedValue; value uniform 10; Outlet: U: type zeroGradient; p: type fixedValue; value uniform 0;
-
@OItoCFD 非常感谢,我们也准备加密网格,但是具体加密的程度还不知道怎么把握,您这边是使所有网格大小尽量均一吗?
-
-
@李东岳
老师好,是我没有表达清楚,是类似于这样
库朗数和deltaT都非常非常小,计算非常慢,到最后deltaT基本变成0了
我设置的最大库朗数是0.9,但是计算出来才0.0012,不应该自动缩小时间步长呀请问老师,这个应该怎么解决,如何提高计算速度?
Courant Number mean: 3.58164e-05 max: 0.901511 Interface Courant Number mean: 2.94865e-08 max: 0.00119605 deltaT = 1.7962e-05 Time = 2.01626 MULES: Solving for alpha.Sarin Phase-1 volume fraction = 6.16252e-05 Min(alpha.Sarin) = 0 Max(alpha.Sarin) = 0.109368 MULES: Solving for alpha.Sarin Phase-1 volume fraction = 6.16252e-05 Min(alpha.Sarin) = 0 Max(alpha.Sarin) = 0.109368 smoothSolver: Solving for alpha.Sarin, Initial residual = 5.46198e-10, Final residual = 5.46198e-10, No Iterations 0 smoothSolver: Solving for Ux, Initial residual = 0.00020754, Final residual = 9.85112e-08, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 0.000110383, Final residual = 4.3286e-08, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 0.000191943, Final residual = 1.79758e-06, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 0.00248148, Final residual = 1.41595e-05, No Iterations 6 time step continuity errors : sum local = 3.84475e-13, global = 1.89675e-14, cumulative = -3.04519e-09 GAMG: Solving for p_rgh, Initial residual = 0.00152746, Final residual = 5.47895e-06, No Iterations 7 time step continuity errors : sum local = 1.48826e-13, global = -6.48062e-15, cumulative = -3.04519e-09 ExecutionTime = 680.02 s ClockTime = 680 s Courant Number mean: 3.57663e-05 max: 0.90151 Interface Courant Number mean: 2.94446e-08 max: 0.00119381 deltaT = 1.79319e-05 Time = 2.01628
-
不好意思老师,是我说的不清楚,发现一直在计算同一个时间,比如一直计算0.1111秒
-
各位老师好:
写了一个非常简单的pisoFoamTemp,测试了一下发现一切正常,用elbow.msh测试的
然后在实际网格中,发现根本就没有计算出来,一直在入口处,都没有进入计算域,但是速度压力求解都正常,请问这是什么原因?
不好意思图片不是很清晰,但是可以看到温度就一直在壁面处没有变化,但是速度计算看上去比较正常,这个是计算150s的结果附上算例文件
链接: https://pan.baidu.com/s/1RaX15gCuychiMvbJBc9MMg 密码: 8gha -
各位老师好:
用自动调节时间步的时候经常发现一直在循环计算某一个时间,这是什么原因,是发散了吗?
因为一开始测试为了简单和利于收敛,设置的层流条件,实际上有部分地方应该是湍流,这个会影响吗?如果降低进口速度,deltaT,有助于改善吗?
-
各位老师好:
运行decomposePar的时候显示如下错误:
decomposePar: error while loading shared libraries: libscotch.so: cannot open shared object file: No such file or directory
此外,如果是单核运算的话,是没问题的。
请问各位老师,这是什么问题,怎样解决啊?
-
@bestucan 非常感谢老师的提醒,让网格变稀疏之后发现了两处可疑的质量较差的地方。尝试之后发现,有一些长方体网格划的太瘦了,ICEM的网格质量检查没有问题,修改之后OpenFOAM中就也没有问题了。
-
@bestucan 非常感谢老师的指点,非常明白,我会再检查一下模型转换的问题。
然后还是想了解一下,请问有没有什么能够快速 定位这些有问题的网格单元的位置 的方法呢?
-
@bestucan 老师您好,非常感谢老师您的回答。经过几天的检查发现,没封闭的原因在于把stl文件导入fluent的时候有的线条和面的格式出现了问题。
但是第一个网格质量较差的原因却没有查找出来,按照fluent对于skewness的定义,skewness应该不能超过1才对,而且我在icem中的网格质量最低是0.6,已经达到较高的标准了。
翻了翻资料说OpenFOAM检查的网格单元和 面元对应顶点的坐标、两个网格单元中心连线向量、面元中心与网格单元中心连线向量、面元的面法相向量 相关,但是没有找到具体的算法
请问老师这个问题怎么排查?或者说怎么才能快速的找到有质量问题的网格?
-
各位老师好,大家好:
使用ICEM画了网格之后checkMesh时显示了好几个问题,但是我不是很能理解是什么意思,也不太清楚往什么方向检查或者修改。但是,我之前在icem里面网格质量检查都显示所有网格都有较高的质量。。。第一个网格问题:
说有1296个高扭曲的面,但是同时我发现,这一点并不影响计算,请问这个问题是说有一些网格质量太差吗?是OpenFOAM对于网格质量的判定标准和icem不一样?Checking geometry... Overall domain bounding box (-72.0987 -6.78501 -2.15365) (59.2013 6.71664 3.24987) ...... Face pyramids OK. ***Max skewness = 4.54714, 1269 highly skew faces detected which may impair the quality of the results <<Writing 1269 skew faces to set skewFaces Coupled point location match (average 0) OK. Failed 1 mesh checks.
第二个问题:
这个就有些严重了,直接导致计算无法进行,请问这个是说边界没有封闭吗:Checking geometry... Overall domain bounding box (-89.7486 -38.8002 1.14635) (-27.4486 24.9498 21.3019) Mesh has 3 geometric (non-empty/wedge) directions (1 1 1) Mesh has 3 solution (non-empty) directions (1 1 1) ***Boundary openness (0.0173317 5.152e-10 -7.40693e-05) possible hole in boundary description. ***Open cells found, max cell openness: 1, number of open cells 2112 <<Writing 2112 non closed cells to set nonClosedCells Minimum face area = 0.00318294. Maximum face area = 0.230028. Face area magnitudes OK. Min volume = 0.000594917. Max volume = 0.0453085. Total volume = 7140.58. Cell volumes OK. ...... Failed 2 mesh checks.
还请各位老师不吝赐教,非常感谢!
-
@fangyuanaza 非常感谢!
-
@东岳 老师好,非常感谢您!
请问是要在输出网格的时候存成ASCII码格式吗?如果只有.msh文件可以用ICEM转换成ascii吗? -
东岳老师说过,fluent3DMeshToFoam对于转换3D网格更好一些,对于我自己使用ICEM画的网格,都很好,这还是我第一次遇到转换失败的问题
另外,这个网格文件打开都是乱码:(3010 (e 1 9211 2 3) ( @J窯? ?@ 鄒W? 棱G? @ef@ 纔
是因为比较大的原因吗?
我自己画的网格文件打开的话都是一行行数字啊 -
大家好,各位老师好:
最近网上下载了一个网格,在用ICEM观察时没有发现错误,
但是使用fluent3DMeshToFoam转换网格时报错:--> FOAM FATAL ERROR: Do not understand characters: on line 19 From function virtual int yyFlexLexer::yylex() in file fluent3DMeshToFoam.L at line 748. FOAM exiting
请问这会是什么原因?怎样才能解决呢?
请问会是用gambit生成和ICEM生成导致的吗?
另外,如果改成fluentMeshToFoam转换的话,还是会出现乱码:
Found end of section in unknown:) Embedded blocks in comment or unknown: ( @�@��@@��@@��������������@��Found end of section in unknown:) �@@��@@�����������@��@@�`Found end of section in unknown:# (����@@�Embedded blocks in comment or unknown:�> ��@@�Found end of section in unknown:> ��@@���@@����@@@����@@����@@��@@�}�@@����@@�Embedded blocks in comment or unknown:�> ���@@�Embedded blocks in comment or unknown:��
-
哎呀,名额满了吗?
没满加一个?
18994065595@163.com -
希望对大家有帮助。
实时更新sci-hub发帖混个积分
-
@东岳 非常感谢老师 !
就是这样迂回的话好像有些难,我还要再思考思考,再尝试尝试 -
你好,我也是课题组里唯一一个做流体的学生。
你这样的话应该非常容易发散,因为k和epsilon对于收敛影响很大,你可以根据估算方法自己算一组k-epsilon值,然后再慢慢调整。
另外,分享一个可能有用的k-epsilon经验计算器:
链接:https://pan.baidu.com/s/1AP6LuHwi9Gb4dfQ9fDgyyg
提取码:12zd希望对你有用
-
@东岳 非常感谢老师的指点,已经成功编译了UEqn
tmp<fvVectorMatrix> tUEqn//fvVectorMatrix UEqn ( fvm::ddt(rho, U) + fvm::div(rhoPhi, U) + turbulence->divDevRhoReff(rho, U) == fvOptions(rho,U)//add this ); fvVectorMatrix& UEqn = tUEqn.ref();//add UEqn.relax(); fvOptions.constrain(UEqn);//add this
但是,还想请问一下老师,怎样才能在rho方程上添加fvOptions?rho不是只由alpha1和rho1,2决定吗?
rho = alpha1*rho1 + alpha2*rho2;
-
CPU型号:i7-8750H 系统:虚拟机 OpenFOAM版本:OpenFOAM-6 4 910.12 2 883.78 1 1060.35
重新跑的结果,这次除了虚拟机什么都没运行。。。
希望能对统计有点用。。。 -
好的,我重新再算一遍
-
CPU型号:i7-8750H 系统:虚拟机 OpenFOAM版本:OpenFOAM-6 4 920.13 2 1098.04 1 572.75
为什么单核的时候反而最快?
-
对了,如果还有同伴也需要用到fvOptions的话,有一个网站推荐给大家:
https://caefn.com/openfoam/fvoptions里面写了在算例中添加fvOptions的一些方法,希望对大家有帮助
-
@马乔 老师您好,非常感谢老师您的指点,我仔细对比了一下,pisoFoam的createFields.H中本来就包含了createFvOptions.H这个文件,现在我加上createFvOptions.H后已经编译成功了。
但是,现在又出现了一个问题,alpha的fvOptions已经可以计算了,但是U在计算过程中却出现了单位不一致的问题, 但是我却检查不出来问题,还请老师您能再指点一下
报错:
--> FOAM FATAL ERROR: incompatible dimensions for operation [U[1 -2 -2 0 0 0 0] ] == [U[0 1 -2 0 0 0 0] ] From function void Foam::checkMethod(const Foam::fvMatrix<Type>&, const Foam::fvMatrix<Type>&, const char*) [with Type = Foam::Vector<double>] in file /home/dyfluid/OpenFOAM/OpenFOAM-6/src/finiteVolume/lnInclude/fvMatrix.C at line 1283.
UEqn.H代码:
fvVectorMatrix UEqn ( fvm::ddt(rho, U) + fvm::div(rhoPhi, U) + turbulence->divDevRhoReff(rho, U) == fvOptions(U)//add this ); UEqn.relax(); fvOptions.constrain(UEqn);//add this if (pimple.momentumPredictor()) { solve ( UEqn == fvc::reconstruct ( ( - ghf*fvc::snGrad(rho) - fvc::snGrad(p_rgh) ) * mesh.magSf() ) ); fvOptions.correct(U);//add this }
-
@东岳 老师您好,添加这个项就是想在不改变网格的条件下,在某个位置设置一个新的入口,设置fvOptions后就能通过只修改topoSet中的区域划分来不断修改新的入口的位置。不然的话,没改变一次新的入口的位置就需要重新划分一下网格,太麻烦了。(/ω\)
之前在pisoFoam中是成功实现了,但是在这里却失败了
pisoFoam中UEqn改的地方:fvm::ddt(U) + fvm::div(phi, U) + MRF.DDt(U) + turbulence->divDevReff(U) == fvOptions(U) ); UEqn.relax(); fvOptions.constrain(UEqn);//
-
各位老师好,大家好:
最近在twoLiquidMixingFoam中添加fvOptions项后wmake时出现了如下错误,显示我的fvOptions没有定义:In file included from my_twoLiquidMixingFoam.C:78:0: alphaDiffusionEqn.H: In function ‘int main(int, char**)’: alphaDiffusionEqn.H:11:10: error: ‘fvOptions’ was not declared in this scope == fvOptions(alpha1)//add ^~~~~~~~~ alphaDiffusionEqn.H:11:10: note: suggested alternative: ‘fvOption_H’ == fvOptions(alpha1)//add ^~~~~~~~~ fvOption_H
但是,我在C文件的开头中添加了fvOptions头文件,此外,Make文件夹中options中也有lfvoptions:
#include "pimpleControl.H" #include "fvOptions.H"//add this
EXE_LIBS = \ ...... -lfvOptions \ -lmeshTools
请问大家这是什么原因?有什么解决办法吗?
期待大家的回复和解答。此外,我还有一个问题,fvOptions和fvOption_H有什么区别呢?我没有找到fvOption_H的文件
另外,我的alphaDiffusionEqn.H是这样写的
{ fvScalarMatrix alpha1Eqn ( fvm::ddt(alpha1) - fvc::ddt(alpha1) - fvm::laplacian ( volScalarField("Dab", Dab + alphatab*turbulence->nut()), alpha1 ) == fvOptions(alpha1)//add ); fvOptions.constrain(alpha1Eqn);//add this alpha1Eqn.solve(); fvOptions.correct(alpha1);//add this alpha2 = 1.0 - alpha1; rhoPhi += alpha1Eqn.flux()*(rho1 - rho2); } rho = alpha1*rho1 + alpha2*rho2;
-
@Leon-Zhao 是的,东岳老师讲的课深入浅出,非常适合初学者。我的C语言和数学基础几乎为零,但是上过课之后也能初级的应用openfoam了,对CFD也有了一些初级的概念。现在我的目标是用of计算大量数据之后再融入神经网络,希望自己能成功吧。
-
新知识点
-
@Samuel-Tu 非常感谢您的回复。讲的非常明白,我之前没有想到这些。其实,我想实现的就是一个热源点慢慢加热最终提高出口温度,或许我应该试试改成热源强度再试一试看。非常感谢您。
-
@Samuel-Tu 非常感谢您的回复。T是瞬态,但是,TEqn跟U只是单向耦合。我的TEqn放在了pisoLoop的外面,这一点我没有考虑过,请问会造成什么影响?因为我看openfoamwiki上icoFoamTemp修改是放在了runTime.loop外面,所以我也放在了同一位置
非常期待您的回复
附上ofwiki的链接链接文本 -
大家好:
请教一个问题,添加温度源的时候,发现一段时间后就不再扩散,明显不符合常识啊,会是什么原因?在我的例子中t=6之后就不再变化,如图:
温度源是写在system/fvOptions里的,SRC是topoSet中分出来的一个区域,附上代码:
TSource { type scalarSemiImplicitSource; active true; scalarSemiImplicitSourceCoeffs { selectionMode cellSet; // all, cellSet, cellZone, points cellSet SRC; volumeMode absolute; // absolute or specific injectionRateSuSp { T (10 0); } } }
求解器是基于pisoFoam写的,其中添加的温度方程参照scalarTransportFoam写的,代码:
fvScalarMatrix TEqn ( fvm::ddt(T) + fvm::div(phi, T) - fvm::laplacian(DT, T) == fvOptions(T) ); TEqn.relax(); fvOptions.constrain(TEqn); TEqn.solve(); fvOptions.correct(T);
请大家不吝解惑,欢迎任何评论
-
东岳老师好,大家好:
我是正在学习OpenFOAM的新手
最近在学习模拟湍流时运行pisoFoam下自带算例porousBlockage时出现错误:Create time Create mesh for time = 0 Reading field p Reading field U Reading/calculating face flux field phi Selecting incompressible transport model Newtonian Selecting turbulence model type laminar Selecting laminar stress model Stokes No MRF models present Creating finite volume options from "constant/fvOptions" Selecting finite volume options model type explicitPorositySource Source: porosity1 - selecting cells using cellZone porousBlockage --> FOAM FATAL ERROR: Cannot find cellZone porousBlockage Valid cellZones are 0() From function void Foam::fv::cellSetOption::setCellSet() in file cellSetOption/cellSetOption.C at line 142. FOAM exiting
请问是什么原因呢?
此外,还想问一下,如果使用标准k-epsilon模型(RAS算例下)模拟湍流,进出口使用什么条件比较合适呢?
请教稳定的两相流压力驱动边界条件
请教稳定的两相流压力驱动边界条件
自动调节时间步长循环计算某一时间步
自动调节时间步长循环计算某一时间步
自动调节时间步长循环计算某一时间步
自动调节时间步长循环计算某一时间步
openfoam求解.msh网格问题
自动调节时间步长循环计算某一时间步
decomposePar出问题
checkMesh中发现的几个网格问题如何检查修改
checkMesh中发现的几个网格问题如何检查修改
checkMesh中发现的几个网格问题如何检查修改
checkMesh中发现的几个网格问题如何检查修改
.msh文件转换网格失败
.msh文件转换网格失败
.msh文件转换网格失败
.msh文件转换网格失败
Daniele Marchisio的CFD Talk登记
发现一个网站实时更新sci-hub
添加源项fvOptions编译出错
关于RAS Model里k-Epsilon中k和Epsilon的取值问题~
添加源项fvOptions编译出错
200万网格并行算力测试(OpenFOAM版本)
200万网格并行算力测试(OpenFOAM版本)
200万网格并行算力测试(OpenFOAM版本)
添加源项fvOptions编译出错
添加源项fvOptions编译出错
添加源项fvOptions编译出错
添加源项fvOptions编译出错
OKS课程买家秀
paraview显示问题
请教扩散现象不随时间变化的原因?
请教扩散现象不随时间变化的原因?
请教扩散现象不随时间变化的原因?
湍流算例porousBlockage出现错误