造福来者，代码贴到下方，主要参考turbulentInlet这个边界，我的目的是每个时间步产生不同的脉动，所以种子用的是时间～

INLET
    {
        type            codedFixedValue;
        redirectType    rampedFixedValue;
        
        code
        #{
            const fvPatch& boundaryPatch = this->patch();
            vectorField U = boundaryPatch.Cf();
            vectorField u = boundaryPatch.Cf();
            
            scalar Uave = 9.65;
            scalar dU = 12.3;
            scalar deltaW = 0.0002;
            scalar C0 = 0.1;
            
            label index = this->db().time().timeIndex();
            Random perturbation(index);      

            forAll(boundaryPatch, faceI)
            {   
                vector deviation
                (
                    1.0 - 2.0*perturbation.scalar01(),
                    1.0 - 2.0*perturbation.scalar01(),
                    1.0 - 2.0*perturbation.scalar01()
                );
                //Info << "random: " << deviation << endl;
                
                scalar y = boundaryPatch.Cf()[faceI].y();
                U[faceI] = vector(Uave + dU*tanh(2*y/deltaW)/2, 0, 0);
                
                u[faceI] = C0*dU*deviation*exp(-sqr(y)/(2*deltaW));
            }

            operator==(U+u);
        #};
        
        codeInclude
        #{
            #include "Random.H"
        #};
        
        value           $internalField;

    }

@李东岳 前辈，解决了，原来是输出精度不够导致的，我的锅:lol:

@cfd-china 谢谢回复，我研究了源码已经搞懂了，顺便附上我总结的链接laplacian项离散解码

众所周知，用p_rgh计算会有较好的稳定性。如果计算结束，我们想要关注压力云图，必定是看p。
因为p_rgh = p - rho*gh
那么问题来了，rho*gh会影响p_rgh的值，所以它也叫虚假压力。
由于N-S方程求解的是静压，p_rgh和rho*gh将p分成了两个部分，两者在不同参考系下一般是不同的，它们是怎么相互作用组合成了我们要的压力p呢？

@random_ran “数学公式”选项，码命令即可

看样子是压缩过度了，如果界面太sharp，CSF模型计算曲率会出问题，代数重构方法的通病

植入后的cavity测试，与时间步无关

@xpqiu 好的，谢谢邱老师，我再试试😄

@xpqiu 好的，我的网格在观测点处太粗了，时间步长选择的和探针的采样频率一致，是不是需要比这个值还小呢？

@xpqiu 是的，定性的可以对上。计算中波动被磨平了很多，您认为是怎么回事呢？

邱老师好，我做了监测点的PSD。请问对比时域的绝对值是不是没意义呢？

水听器是没于水中测声压的探头。我有一组实验数据，我想把模拟的检测结果和实验对比。我不知道后处理的时候应该直接用监测点的瞬时值呢？还是应该用压力脉动值？请大家给些建议

如题，比如相体积分数，有人尝试过么？

如题，目前的方法是blockMesh先生成方形单层mesh，再利用extrudeMesh生成360的3D网格。
但是出现个问题，里面的面不会自动的成为interior，不太清楚怎么处理。

以下是extrudeMeshDict

constructFrom patch;
sourceCase "$FOAM_CASE";

sourcePatches (side1);
exposedPatchName side2;

extrudeModel    sector;

nLayers             10;//90;

sectorCoeffs
{
    point       (0 0 0);
    axis        (0 1 0);
    angle       20;//360;
}

flipNormals false;

mergeFaces  false;

@南风 他这个意思是群泡传质速率，不是单泡的。另外这个值与你其他参数取值也有关系

通过求解ode方程，得到了一系列采样点，想对其进行fft分析，计算出主频对应的幅值，不知道怎么用of里的fft

看到个帖子不错，记录一下。关键词：householder transform