Anisotropy of thermal conductivity

[Anne Schulz]

Hello everyone,
I would like to incorporate the anisotropy of thermal conductivity into my simulations.

For this purpose, I have structured the specification of the material properties as follows:

MATERIAL_PROPERTY rotliegend
    ID 2
    POROSITY 0.0794
    ROCK_DENSITY 2440
    PERMEABILITY
        PERM_HORIZONTAL 5.5e-15                           
        VERTICAL_ANISOTROPY_RATIO 5.500000000000001e-16
    /
    HEAT_CAPACITY 755
    CHARACTERISTIC_CURVES default
    THERMAL_CHARACTERISTIC_CURVES rotliegend
    /
    THERMAL_CHARACTERISTIC_CURVES rotliegend                                    
        THERMAL_CONDUCTIVITY_FUNCTION LINEAR_RESISTIVITY
            THERMAL_CONDUCTIVITY_DRY 2.95
            THERMAL_CONDUCTIVITY_WET 2.49
            LINEAR_RESISTIVITY_COEFFICIENTS 0.7167d0 2.949D4                   
            ANISOTROPY_RATIO_X  1.0000D+0
            ANISOTROPY_RATIO_Y  1.0000D+0
            ANISOTROPY_RATIO_Z  0.85                    
        END
        TEST
    END

I used the specifications for the THERMAL_CHARACTERISTIC_CURVES according to an example from the user documentation, but unfortunately, I am receiving the following error message:

ERROR: Cannot combine material-based thermal conductivity input format with thermal characteristic curves. Use TCC with "DEFAULT" specification instead.

I then changed the THERMAL_CONDUCTIVITY_FUNCTION to DEFAULT and commented out the LINEAR_RESISTIVITY_COEFFICIENTS, but I still get the same error message.

MATERIAL_PROPERTY rotliegend
    ID 2
    POROSITY 0.0794
    ROCK_DENSITY 2440
    PERMEABILITY
        PERM_HORIZONTAL 5.5e-15                           
        VERTICAL_ANISOTROPY_RATIO 5.500000000000001e-16
    /
    HEAT_CAPACITY 755
    CHARACTERISTIC_CURVES default
    THERMAL_CHARACTERISTIC_CURVES rotliegend
    /
    THERMAL_CHARACTERISTIC_CURVES rotliegend                                    
        THERMAL_CONDUCTIVITY_FUNCTION DEFAULT
            THERMAL_CONDUCTIVITY_DRY 2.95
            THERMAL_CONDUCTIVITY_WET 2.49
            # LINEAR_RESISTIVITY_COEFFICIENTS 0.7167d0 2.949D4                   
            ANISOTROPY_RATIO_X  1.0000D+0
            ANISOTROPY_RATIO_Y  1.0000D+0
            ANISOTROPY_RATIO_Z  0.85                    
        END
        TEST
    END

Can anyone help me resolve this issue?
Best regards,
Anne

[Eduardo Basurto]

Hi Anne,

I ran into this same issue recently. The solution is to also set THERMAL_CHARACTERISTIC_CURVES for ALL of your material properties even though you won't necessarily be setting the anisotropy ratio x,y, z value to anything other than 1.

-Eddie

[Anne Schulz]

Thanks, Eddie. 

I'll give it a try and hope there aren't any more hidden bugs in my code. =)

Best regards

[Eduardo Basurto]

No problem, let me know if this resolves it, otherwise I can look over your input deck.

-Eddie