scalarTransportFoam中如何实现随空间变化的rhoc密度场？

[大喵]

问题：在scalarTransportFoam中，目前是把rhoc当做常数，并除到K下面，得到DT

请问如果rho随空间变化（比如说现有 volScalarfield rho），该如何实现该方程的求解呢？

[李东岳]

在代码里加入

volScalarField rho
    (
        IOobject
        (
            "rho",
            runTime.timeName(),
            mesh,
            IOobject::MUST_READ,
            IOobject::AUTO_WRITE
        ),
        mesh
    );

然后直接用就行