CFD中文网

coolhhh

@coolhhh 如果还有问题，就用最笨的方法，用probes，人工生成监测点位置，这样结果就没问题。但probes没法监测壁面边界的数据，比如wallShearStress，见帖子probes和sample提取wallShearStress问题

coolhhh

@lizhisongsjtu
我之前是用of2206， axis用distance试下。另外网格划分有1000个吗，1000个点是否会很密了？

    x_0_1 
    { 
        type        uniform; 
        axis        distance; 
        start       (0.00010000 0.00000100 -1.72500000); 
        end         (0.00010000 1.99999900 -1.72500000); 
        nPoints     100; 
    }

学

取一个或多个平行于壁面的面（Kim的图也是不同高度的，y^+=5.39和149.23），使用surfaces工具输出这个面上的速度，然后每个面计算速度的展向、流向空间相关函数就行了。
可以直接计算二维的相关。
也可以固定z坐标，一条线一条线地计算流向的相关函数，再做个平均。

学

@bcp 我搞错了，之前没看对你计算rField的公式。我以为是你rField的计算公式不对。

你把报错代码发一下。
debug，分别输出rho、U、rField的量纲看一看

学

UEqn当中各项的量纲是 $[\mu][U]/L^2=M/(L^2T^2)$
rField的量纲应该是$T^{-1}$才对，但看起来你写出来的并不是这样
你可以在定义场的时候显式地指定量纲

volScalarField B                                            
(
    IOobject
    (
        "B",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    mesh,
    dimensionedScalar(dimless, 0.0)
);

学

平均流速就是用瞬时的流场速度U进行三维空间上的积分。可以看看代码
https://develop.openfoam.com/Development/openfoam/blob/OpenFOAM-v2012/src/fvOptions/sources/derived/meanVelocityForce/meanVelocityForce.C

和重力不挂勾的，重力产生的加速度是它自己的效果，记为A；这个体积源项本身会产生一个加速度，记为B。
A不变，B会一直变化以使得平均速度保持为你设定的那个值。

coolhhh

@lb0331 参考这个帖子如何控制使用pisoFoam计算时的库朗数？，计算时候严格让CFL完全小于1是很难的，可以适当大一点，只要最后结果计算是合理的。

coolhhh

@郑学习好的，感谢回复

coolhhh

@郑学习现在我也遇到了buff/cache过大导致计算变慢问题。想问下现在把swap拉高到256，后面还会出现buff/cache过大问题吗？

CFD中文网

讲师

帖子

成员列表