CFD中文网

@李东岳 多谢李老师回复，我在cfd online（https://www.cfd-online.com/Forums/openfoam-post-processing/107491-reconstructpar-error.html
）上看到了跟我类似的问题，说是RAM内存不够的问题，我不是很能理解这个解释，李老师您清楚吗？

各位老师好，我最近在并行计算后用reconstructPar重组流场出现了报错：

我原本以为是存储空间不够导致写入不了，后面发现不是。机器的存储空间剩余了很多。。（但是删掉一些算例文件夹的大体积文件却莫名其妙地可以继续reconstructPar）我不是很清楚这个是为什么。。
请问各位老师知道怎么解决这个问题吗？

各位老师好，我最近在做自然对流的模拟时发现温度场出现了比较大的振荡，感觉是发散了（因为流场最高温度已经大于了底部边界给定的温度）：

但是Co数是没问题的。。
请问有老师知道这个问题是什么原因造成的，该怎么解决吗？

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

@尚善若水 你好，controlDict里面的代码是这样的：

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application     buoyantPimpleFoam;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         500;

deltaT          0.001;

writeControl    runTime;

writeInterval   10;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

adjustTimeStep  true;

maxCo           0.5;

各位老师，请问有谁知道在自适应步长设置（runTimeModifiable true;）下统一输出的精度吗？

不统一的话可能会导致程序续跑时读不了latestTime的问题，因为timePrecision为6时，计算时间步可能非常小，但是输出有效位就会变得特别长，而程序续跑又只能读到前6个有效数字。

烦请各位老师指教。

@李东岳 李老师，我发现应该是dev2(T(……)）这一小段出的问题，删掉这个转置操作就能编译成功。因为我是在传热求解器的基础上进行修改的，T在createFields.H里声明了是温度场，所以是不是程序将T识别成了温度而进行不了转置操作了呢？

@李东岳 李老师好，我的nuEff定义在了构建动量方程之前：

	volScalarField nu = turbulence->nu();
	volScalarField nut = turbulence->nut();
	volScalarField nuEff = nu + nut;

而且，报错信息中没有指出nuEff未被声明，所以我感觉不是这个的原因

@SuperBoy >divDevReff(U)的代码参考自linearViscousStress.C：

template<class BasicTurbulenceModel>
Foam::tmp<Foam::fvVectorMatrix>
Foam::linearViscousStress<BasicTurbulenceModel>::divDevRhoReff
(
    volVectorField& U
) const
{
    return
    (
      - fvc::div((this->alpha_*this->rho_*this->nuEff())*dev2(T(fvc::grad(U))))
      - fvm::laplacian(this->alpha_*this->rho_*this->nuEff(), U)
    );
}

各位老师好， 我最近将动量方程中turbulence->divDevReff(U)这一项还原成-fvc::div((nuEff*dev2(T(fvc::grad(U)))))-fvm::laplacian(nuEff, U)，在编译时出现了报错的情况，如下图所示：

我不是很清楚报错的原因，请问有老师知道该怎么解决吗？

@coolhhh 感谢您的回复，应该就是求解器中那些非线性模型导致的。感觉openfoam做无量纲化有不少trick。。

各位老师，我最近在计算自然对流的时候，发现在压力参考点（pRefCell = 0，pRefValue = 0.0）处，随着计算推进，压力跟速度的数值会变得异常的奇怪：压力固定为0，速度也非常小，但是其附近网格的压力速度的绝对值却很大，也就是参考点处的压力速度梯度很大，导致计算也特别容易发散。我不是很清楚这是为什么？
如下图所示，一开始几个迭代步没有问题，但是慢慢地，参考点处的压力速度就变得奇怪了起来。
30多个时间步之后的速度（上下边界为no slip）：

压力（fixedFluxPressure）：

温度（上下边界分别为固定190K和300K）：

有老师遇到过类似的情况知道怎么解决么？

@李东岳 就是长度，长度缩小10倍，体积缩小1000倍

@李东岳 东岳老师好，我做的是一个三维案例，所以线尺度缩小10倍，计算域整个体积就缩小了1000倍。我是这么想的，比方说下面的离散方程（忽略对流项），如果取alpha是线尺度缩放因子，那么，为了求解得到与缩放前相同的的T，只要将时间相应的缩小alpha^2倍即可。但是实际求解起来得到的结果却不一样，不知道是哪一点出现了问题

各位老师好，
我最近在做网格尺寸的缩放，但是碰到了求解结果不一致的情况。
比方我将计算域（单位网格）的线尺寸缩小10倍，那么网格体积和表面积将分别缩小10^3和10^2倍。对比离散后的控制方程，如果要得到缩放前的计算结果（速度或温度），那么对应的时间长度应该取缩小10^2倍，但是，这样求解出来结果与缩放前并不一致，请问有老师知道这是为什么吗？