Loading tests/test_ThermoGISDoublet_Benchmark.py +24 −19 Original line number Diff line number Diff line Loading @@ -189,37 +189,42 @@ class ThermoGISDoubletBenchmark(TestCase): utc_properties = (self.setup_template_utc_properties_builder() .setOpexPerPower(100) .setOpexBase(0) .setHpDirectHeatInputTemp(70) .setUseHeatPump(True) .setHpApplicationMode(True) .setHpDirectHeatInputTemp(80) .setDhReturnTemp(35) .setViscosityMode(ViscosityMode.KESTIN) .build()) # Create an instance of a ThermoGISDoublet doublet = ThermoGISDoublet(Mockito.mock(Logger), RNG(0), utc_properties) # Setting input values thickness = 100 permeability = 175 transmissivity = thickness * permeability doublet.setReservoirTemp(50) doublet.setDepth(2000) doublet.setDhReturnTemp(35) doublet.setInjectionTemp(15) doublet.setNtg(1.0) doublet.setNoStimulation() temperature = 50 input = DoubletInput( -999.0, # unknowninput thickness, transmissivity, 0.0, # transmissivityWithNtg 1.0, # ntg 2000.0, # depth 0.0,# porosity temperature, None, # ates input ) doublet = ThermoGISDoublet(Mockito.mock(Logger), RNG(0), utc_properties) # Act doublet.calculateDoubletPerformance(-999, thickness, transmissivity, False) results = doublet.calculateDoubletPerformance(input) # Assert self.assertTrue(np.isclose(17499.99940, transmissivity, 0.001)) self.assertTrue(np.isclose(163.99771118164062, doublet.getFlowrate(), 1)) self.assertTrue(np.isclose(60, doublet.getDrawdownPressure() / 1e5, 0.001)) self.assertTrue(np.isclose(5.67, doublet.getHeatPowerPerDoublet(), 0.001)) self.assertTrue(np.isclose(955.27, doublet.getWellDistance(), 0.001)) self.assertTrue(np.isclose(3.51, doublet.getCop(), 0.001)) self.assertTrue(np.isclose(10.9, doublet.getUtcPeurctkWh(), 0.001)) self.assertTrue(np.isclose(17.56, doublet.getSumcapex(), 0.001)) self.assertTrue(np.isclose(163.99771118164062, results.flow(), 1)) self.assertTrue(np.isclose(60, results.pres(), 0.001)) self.assertTrue(np.isclose(4.97566556930542, results.power(), 0.001)) self.assertTrue(np.isclose(989.2578125, results.welld(), 0.001)) self.assertTrue(np.isclose(4.383212200392486, results.cop(), 0.001)) self.assertTrue(np.isclose(10.401010755009017, results.utc(), 0.001)) self.assertTrue(np.isclose(16.50359210062243, results.capex(), 0.001)) def test_calculateDoubletPerformance_ORC(self): Loading Loading
tests/test_ThermoGISDoublet_Benchmark.py +24 −19 Original line number Diff line number Diff line Loading @@ -189,37 +189,42 @@ class ThermoGISDoubletBenchmark(TestCase): utc_properties = (self.setup_template_utc_properties_builder() .setOpexPerPower(100) .setOpexBase(0) .setHpDirectHeatInputTemp(70) .setUseHeatPump(True) .setHpApplicationMode(True) .setHpDirectHeatInputTemp(80) .setDhReturnTemp(35) .setViscosityMode(ViscosityMode.KESTIN) .build()) # Create an instance of a ThermoGISDoublet doublet = ThermoGISDoublet(Mockito.mock(Logger), RNG(0), utc_properties) # Setting input values thickness = 100 permeability = 175 transmissivity = thickness * permeability doublet.setReservoirTemp(50) doublet.setDepth(2000) doublet.setDhReturnTemp(35) doublet.setInjectionTemp(15) doublet.setNtg(1.0) doublet.setNoStimulation() temperature = 50 input = DoubletInput( -999.0, # unknowninput thickness, transmissivity, 0.0, # transmissivityWithNtg 1.0, # ntg 2000.0, # depth 0.0,# porosity temperature, None, # ates input ) doublet = ThermoGISDoublet(Mockito.mock(Logger), RNG(0), utc_properties) # Act doublet.calculateDoubletPerformance(-999, thickness, transmissivity, False) results = doublet.calculateDoubletPerformance(input) # Assert self.assertTrue(np.isclose(17499.99940, transmissivity, 0.001)) self.assertTrue(np.isclose(163.99771118164062, doublet.getFlowrate(), 1)) self.assertTrue(np.isclose(60, doublet.getDrawdownPressure() / 1e5, 0.001)) self.assertTrue(np.isclose(5.67, doublet.getHeatPowerPerDoublet(), 0.001)) self.assertTrue(np.isclose(955.27, doublet.getWellDistance(), 0.001)) self.assertTrue(np.isclose(3.51, doublet.getCop(), 0.001)) self.assertTrue(np.isclose(10.9, doublet.getUtcPeurctkWh(), 0.001)) self.assertTrue(np.isclose(17.56, doublet.getSumcapex(), 0.001)) self.assertTrue(np.isclose(163.99771118164062, results.flow(), 1)) self.assertTrue(np.isclose(60, results.pres(), 0.001)) self.assertTrue(np.isclose(4.97566556930542, results.power(), 0.001)) self.assertTrue(np.isclose(989.2578125, results.welld(), 0.001)) self.assertTrue(np.isclose(4.383212200392486, results.cop(), 0.001)) self.assertTrue(np.isclose(10.401010755009017, results.utc(), 0.001)) self.assertTrue(np.isclose(16.50359210062243, results.capex(), 0.001)) def test_calculateDoubletPerformance_ORC(self): Loading