rhoSimpleFoam：k-epsion 湍流Pr数

[Exthan]

@cccrrryyy 你好，我将EddyDiffusivity里将dimensionedScalar Prt改成了volScalarField Prt,修改完成之后，编译也完成，但是想查看新的Prt是否参与计算，所以想输出结果查看，可在reconstruct算例时出现错误，错误信息是T的边界条件segmentation（如下图）

请问修改完EddyDiffusivity文件之后，是否还需要修改其他地方呢？

[Exthan]

@cccrrryyy 老师，不可压缩里通常是在solver里面直接用新的alpha+alphat来代替原来的alphaEff。可压缩哪种方便点,我的prt与mut，mu，nu,alpha,有关.所以感觉在eddyDiffusivity文件中修改更为方便。但是编译后，如何去查看是否调用新的prt，我在solver/creatFields中IOonbject了Prt,但并未生成文件，而且externalWallHeatFlux在映射生成温度T文件时候，也出错了。请问除了修改EddyDiffusivity文件是否还需要修改其他地方？2.如果想要判断新Prt是否成功读入应如何查看？

[cccrrryyy]

@exthan 你的目的是把Prt改成所说的和mut，mu，nu，alpha等等有关的那个表达式对吧？你改的eddyDiffusivity文件路径是啥，我瞅一眼看看，好久没搞可压缩啦。

你判断修改Prt是否成功的方法应该没问题，但是没有输出有点奇怪，按道理说即便是报错了应该也会输出的。如果这种间接的方式还不行，我一般会直接在程序里面让它输出，就是 Info<< 这样子，缺点是这样会输出每一步每一个cell上你要的数值，优点是你一下就能看出来有没有问题。可以把网格弄很小，反正只是测试程序嘛。

[Exthan]

@cccrrryyy
源文件在src/TurbulenceModels/compressible/EddyDiffusivity
我不知道如何传输文件。我可以把我的文件发送给你。这是我改写的Prt

[Exthan]

@cccrrryyy 你好，我刚在湍流模型里面Info,但是结果并未输出，您说的在.C文件中直接输出要如何写？我看了之前的OF文件读取，但是并不理解。

如果直接在.c中如上图输入，编译显示Prt无定义，请问该如何输出

[Exthan]

@cccrrryyy 我在.C文件中直接使用了
Prt.write(),结果生成的文件

打开之后却都是边界条件，这是什么原因？

[cccrrryyy]

@exthan 因为你的Prt是定值呀，这个文件不止包含边界条件，内部场也是有的，只不过是uniform 0.85。

你用info输出那里是不是应该直接改为Prt_，而不是Prt()。

[cccrrryyy]

@exthan 改EddyDiffusivity这一部分，可否直接自己建一个变量，只不过它的计算基于mut，mu这些，然后求解器那块全部使用这个自建的变量。就是说，原有的Prt全部保留，只是增加一个东西，然后求解器用的是新加的这个。

[Exthan]

@cccrrryyy 请问如何在.c文件中调用mu.nu这些量呢？还有网格，代码是如何表现，我试了调用，但总是报错，显示未定义，但是我已经在creatFields定义场了。

[cccrrryyy]

@exthan 这些量本身就是有的吧，mu_，nu_这样子？

你试过调用但显示未定义这个过程是什么样的？

[Exthan]

@cccrrryyy 编译solver的时候直接报错，显示没定义