TNO Intern

Commit 1a4447dc authored by Florian Knappers's avatar Florian Knappers
Browse files

setup doublet test (fails for now)

parent 8acb04e7
Loading
Loading
Loading
Loading
+20 −10
Original line number Diff line number Diff line
import numpy as np

from pythermogis.workflow.utc.doublet import calculate_doublet_performance, DoubletInput, DoubletOutput
from pythermogis.workflow.utc.utc_properties import UTCConfiguration

@@ -5,23 +7,31 @@ from pythermogis.workflow.utc.utc_properties import UTCConfiguration

def test_calculate_doublet_performance_runs_with_default_props():
    # Arrange: instantiate default UTCConfiguration
    props = UTCConfiguration()
    props = UTCConfiguration(
        viscosity_mode="kestin",
        dh_return_temp=40,
    )

    # Create a minimal valid DoubletInput
    input_data = DoubletInput(
        unknown_input_value=0.0,
        unknown_input_value=-999.0,
        thickness=100.0,
        transmissivity=5.0,
        transmissivity_with_ntg=3.0,
        ntg=0.1,
        depth=1500,
        porosity=0.2,
        temperature=80.0,
        transmissivity=17500.0,
        transmissivity_with_ntg=0.0,
        ntg=1.0,
        depth=2000,
        porosity=0.0,
        temperature=76.0,
    )

    # Act
    result = calculate_doublet_performance(props, input_data)
    print(result)

    # Assert
    assert isinstance(result, DoubletOutput)
    assert np.isclose(result.flow, 227.2757568359375, rtol=1)
    assert np.isclose(result.pres, 60, rtol=0.001)
    assert np.isclose(result.utc, 6.616096470753937, rtol=0.001)
    assert np.isclose(result.welld, 1159.17968, rtol=0.001)
    assert np.isclose(result.power, 8.636903762817383, rtol=0.001)
    assert np.isclose(result.cop, 13.627557754516602, rtol=0.001)