setPorousRegionSolidFields.H 613 Bytes
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
    fvMesh& mesh = porousSolidRegions[i];

    basicSolidThermo& thermo = porousSolidThermos[i];
    const volScalarField& betav = betavSolid[i];

    tmp<volScalarField> trho = thermo.rho();
    const volScalarField& rho = trho();

    tmp<volScalarField> tcp = thermo.Cp();
    const volScalarField& cp = tcp();

    tmp<volScalarField> tkappa = betav*thermo.K();
    //tmp<volSymmTensorField> tkappa = thermo.directionalK()*betav;

    const volScalarField& kappa = tkappa();
    //const volSymmTensorField& K = tK();

    volScalarField& T = thermo.T();

    IObasicSourceList& sources = solidHeatSources[i];