solar-copper-extraction/main_final.py at main · Clare2D/solar-copper-extraction · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
# -*- coding: utf-8 -*-
"""
Created on Tue Sep 02 17:51:39 2014

@author: cmurray
"""

import start
import TrnsysCaller
import HeatCalcs
import results
import HeapSim
from matplotlib.pyplot import close

close('all')

dt = 0.25
time = 8760
h_heap = 6.
cu_concentration = 0.005     # 0.5%
pyr_concentration = 0.02     # 2%
t_initial=288.
t_pond1=t_pond2=t_heap_base=t_heap_top=t_initial
a_heap=200000.
case=0

# VARIABLES
g1 = 10.                      # kg/m2hr   Flow rate of fluid through the heap
collector_area = 0;
boiler_on = False
pond_cover = True

if g1==5.:
    t_req=50.+273
if g1==10.:
    t_req=55.+273

results_store={}
initial = start.InitialCalcs()
initial.calc_weather_array()
initial.initial_properties(h_heap,g1,a_heap)
initial.input_definition(results_store,case,boiler_on,pond_cover,cu_concentration,pyr_concentration,g1,h_heap)

initial_run = 1
# Initial TRNSYS Run
trn = TrnsysCaller.TrnsysCaller(time,dt)
if collector_area==0:
    solar_on=False
    trn.call_trnsys(0)
else:
    solar_on=True
    trn.edit_deck_file(collector_area)                                          # Sets the collector area
    trn.write_trnsys_input_file(initial_run,results_store,temp_i=t_pond2)     # Creates the file to export to Trnsys
    trn.call_trnsys(1)                                                          # Creates the file from Trnsys

# Initial HeapSim Run
heapsim = HeapSim.HeapSimProcessing(cu_concentration,time,initial=initial)
if boiler_on == True:
    t_hs=t_req
else:
    t_hs=t_initial
heapsim.heapsim_inputs(initial_run,boiler_on,g1,t_req,results_store,temp_i=t_hs)
heapsim.heapsim_run()
heapsim.heapsim_results(results_store)
temps_heap_base = heapsim.heap_temp_base(results_store)
temps_heap_top = heapsim.heap_temp_top(results_store)

# Initial Pond instances
pond1 = HeatCalcs.Pond(time, dt, initial.mfr_heap, 0, pond_cover, results_store, initial = initial, trn = trn)
if solar_on == False:
    pond2 = HeatCalcs.Pond(time, dt, initial.mfr_heap, 0, pond_cover, results_store, initial = initial, trn = trn)
else:
    pond2 = HeatCalcs.Pond(time, dt, initial.mfr_heap, 1, pond_cover, results_store, initial = initial, trn = trn)

# Energy Calculations
trn_convergence_checker = True
heapsim_convergence_checker = True
while ((trn_convergence_checker == True or heapsim_convergence_checker == True)):
    for i in range(int(time/dt)):
        t_pond1_new = pond1.p_calculations(i,t_pond1,t_heap_base,t_pond2)
        t_pond2_new = pond2.p_calculations(i,t_pond2,t_pond1,t_heap_top)

        t_heap_base = temps_heap_base[int(i*dt)]
        t_heap_top = temps_heap_top[int(i*dt)]

        t_pond1 = t_pond1_new
        t_pond2 = t_pond2_new

    pond1.condense_results('pond1')
    pond2.condense_results('pond2')

    # Check for Convergence: TRNSYS
    if solar_on == False:
        trn_convergence_checker = False
    else:
        trn_convergence_checker = trn.convergence_check(0.5,pond2.results_pond)
        if trn_convergence_checker == True:
            t_pond1 = t_pond2 = t_heap_base = t_heap_top = t_initial
            trn.write_trnsys_input_file(0,results_store)      # Creates the file to export to Trnsys
            trn.call_trnsys(1)
            pond1 = HeatCalcs.Pond(time, dt,initial.mfr_heap, 0,pond_cover,results_store, initial = initial, trn = trn)
            pond2 = HeatCalcs.Pond(time, dt,initial.mfr_heap, solar_on,pond_cover,results_store, initial = initial, trn = trn)

    # Check for Convergence: HeapSim
    if boiler_on == True:
        heapsim_convergence_checker = False
    else:
        heapsim_convergence_checker = heapsim.convergence_check(0.5,results_store)
        if heapsim_convergence_checker == True:
            heapsim.heapsim_inputs(0,boiler_on,g1,t_req,results_store)
            heapsim.heapsim_run()
            heapsim.heapsim_results(results_store)
            temps_heap_base = heapsim.heap_temp_base(results_store)
            temps_heap_top = heapsim.heap_temp_top(results_store)

    # Calculate heap generation and boiler requirements
    heapsim.heap_generation(results_store,initial.mfr_heap,pond2.cp_fluid,boiler_on,t_req)
    heapsim.boiler(t_req,results_store,boiler_on)
    t_pond1 = t_pond2 = t_heap_base = t_heap_top = t_initial

trn.condense_results(results_store)

# Results
res= results.ShowResults(time, dt,a_heap, initial = initial)
res.print_results(results_store)
res.plot_results(results_store,collector_area)
res.daily_heat_profile(100,results_store)