CFD中文网

之前每个场的Dc值一样，这里我想设置每一个field的Dc值不一样，求问一下这里代码怎么修改一下：

{    
    
    surfaceScalarField phichems( phi*1.0 );
    forAll(fluidPhase.field, i)
    {
        fvScalarMatrix chemTransportEqn
        (
            fvm::ddt(fluidPhase.field[i])
            + fvm::div(phichems,fluidPhase.field[i],"div(phi,chems)")
            - fvm::laplacian(Dc,fluidPhase.field[i])  //想修改这一行的Dc
        );

        chemTransportEqn.relax();
        chemTransportEqn.solve();
        fluidPhase.field[i].correctBoundaryConditions();
    }
}

!
A物质在设置的反应区域有初始值，A的源项与B有关，计算B的分布。

那么A与B的反应应该发生在红色框，但是计算的B的值在设置区域后面

请问是对流项太强吗，反应后生成的B都被冲到后面的区域了,我将流速降低还是这样的。

这样啊谢谢，我看论文中测试的时间都没有这么久，就看着能加速反应速率，时间加快些不

求解对流扩散方程，P^k是变量， S^k是源项， S^k形式如下， k^on, k^off, k^cat, 是系数，z7, z2, [z2:z10]等为反应变量。
现在反应太慢了，我增大 2倍k^on, k^off, k^cat, 还可以算，再增大计算就发散了。

请问这里可以怎么处理呢？谢谢大家！！
求解器设置

p
{
    solver          PCG;
    preconditioner  DIC;
    tolerance       1e-06;
    relTol          0;
}

pFinal
{
    $p;
    relTol          0;
}

U 
{
    solver          GAMG;
    preconditioner  DILU;
    smoother        GaussSeidel;
    tolerance 1e-05;
    relTol 0.1;
};