[SeMar]

Dear foamers,

I'm trying to modify the icoIbFoam in the latest foam-extend release 3.2. The idea is, adding a temperature diffusion/convection-equation to the code in order to simulate a conjugate-heat-transfer problem. So far I've implemented a modified "DT" Field

DT = cellIbMask * lambdaF / (rhoF * CpF) + (1-cellIbMask) * lambdaS / (rhoS * CpS);

All the values xF, xS (fluid, solid) are lookup values from the transportPropertiesDict.

In the solverloop I added

fvScalarMatrix TEqn
        (
            fvm::ddt(T)
            + fvm::div(phi, T)
            - fvm::laplacian(DT, T)
        );
TEqn.solve();
T.correctBoundaryConditions();

While temperature is being calculated in the fluid cells, the value of the Ib-cells does not change. I tried several combinations of 'setDeadCellValue' and 'fixesValue' in the 0/T BC file and also erased the 'T.correctBoundaryConditions();' without success. Setting a temperature manually in the code however works. So I guess there must be something preventing the immersed boundary cells from being solved. Does anybody have a clue, where to find the corresponding lines of code?

Thanks in advance

Sebastian

[supersoldier]

Recently , I also want to do some work like yours, so have you had success?

[SeMar]

No, unfortunately not. But I'm still very interested in a solution.

[supersoldier]

Originally Posted by SeMar

Dear foamers,

I'm trying to modify the icoIbFoam in the latest foam-extend release 3.2. The idea is, adding a temperature diffusion/convection-equation to the code in order to simulate a conjugate-heat-transfer problem. So far I've implemented a modified "DT" Field

DT = cellIbMask * lambdaF / (rhoF * CpF) + (1-cellIbMask) * lambdaS / (rhoS * CpS);

All the values xF, xS (fluid, solid) are lookup values from the transportPropertiesDict.

In the solverloop I added

fvScalarMatrix TEqn
        (
            fvm::ddt(T)
            + fvm::div(phi, T)
            - fvm::laplacian(DT, T)
        );
TEqn.solve();
T.correctBoundaryConditions();

While temperature is being calculated in the fluid cells, the value of the Ib-cells does not change. I tried several combinations of 'setDeadCellValue' and 'fixesValue' in the 0/T BC file and also erased the 'T.correctBoundaryConditions();' without success. Setting a temperature manually in the code however works. So I guess there must be something preventing the immersed boundary cells from being solved. Does anybody have a clue, where to find the corresponding lines of code?

Thanks in advance

Sebastian

what do you mean by "Setting a temperature manually in the code however works" ? I posted a thread at http://www.cfd-online.com/Forums/ope...tml#post606287 . #67.

Can you give me some advance?

[SeMar]

By manually I mean assigning a value to the field.

T = 300.0;

This would result in every cell of field T including immersedboundary cells having a value of 300.0.