Process Development Heat Integration

Author

Student1, Student2, Student3

Published

February 1, 2024

source("../rosmose-setup.R", local = knitr::knit_global())

Introduction

This report presents the results of your heat integration simulations. Each run will generate an output json file stored in the “Rosmose/result” folder and a “frontend.html” file stored in “Rosmose” folder.

The json files will be numbered as 1, 2, 3,,,,,, for each run. However, the “frontend.html” file is overwritten each time hence to keep previous versions you must rename the file before running your simulation again.

once your models and units are ready you can run the rsomose simulation by pressing the “Preview” button in the top right corner or using the command line “quarto preview” in the terminal at the correct directory location.

A ready example is given to demonstrate extraction of data from an Aspen flowsheet, changing the process parameters in Aspen through your Rmd file and displaying your integration results.

all Aspen related files are found in this folder “model/AspenModel/.”

The “Aspen_methanol_process.rmd” file is coupled with the “methanol_process_PD1.bkp” file. The same example Aspen process flowsheet used in all tutorials throughout the course.

Methanol Process Aspen model

MODEL Methanol
Software Location Comment
ASPEN model_process_PD1.bkp
MODEL INPUTS Methanol
Name Path Value Units Comments
Electrolyzer_size /Data/Streams/POWER/Input/POWER -20000 kW
MODEL OUTPUTS Methanol
Name Path Units Comments
T_AIR_1 /Data/Streams/AIR1/Output/TEMP_OUT/MIXED C
T_AIR_2 /Data/Streams/AIR2/Output/TEMP_OUT/MIXED C
heater_1 /Data/Blocks/HEATER1/Output/QCALC kW
T_FUEL_1 /Data/Streams/FUEL1/Output/TEMP_OUT/MIXED C
T_FUEL_2 /Data/Streams/FUEL2/Output/TEMP_OUT/MIXED C
heater_3 /Data/Blocks/HEATER3/Output/QCALC kW
T_FUEL_3 /Data/Streams/FUEL3/Output/TEMP_OUT/MIXED C
heater_4 /Data/Blocks/HEATER4/Output/QCALC kW
T_FUEL_4 /Data/Streams/FUEL4/Output/TEMP_OUT/MIXED C
heater_5 /Data/Blocks/HEATER5/Output/QCALC kW
T_AIR_OUT /Data/Streams/AIROUT/Output/TEMP_OUT/MIXED C
T_AIR_4 /Data/Streams/AIR4/Output/TEMP_OUT/MIXED C
cooler_1 /Data/Blocks/COOLER1/Output/QCALC kW
T_FUEL_OUT /Data/Streams/FUELOUT/Output/TEMP_OUT/MIXED C
T_LTFUEL /Data/Streams/LTFUEL/Output/TEMP_OUT/MIXED C
cooler_2 /Data/Blocks/COOLER2/Output/QCALC kW
MSC_cooler_1_Tin /Data/Blocks/C1/Output/B_TEMP/1 C
MSC_cooler_1_Tout /Data/Blocks/C1/Output/COOL_TEMP/1 C
MSC_cooler_1_Duty /Data/Blocks/C1/Output/QCALC/1 kW
MSC_cooler_2_Tin /Data/Blocks/C1/Output/B_TEMP/2 C
MSC_cooler_2_Tout /Data/Blocks/C1/Output/COOL_TEMP/2 C
MSC_cooler_2_Duty /Data/Blocks/C1/Output/QCALC/2 kW
MSC_power_tot /Data/Blocks/C1/Output/WNET kW
T_S1 /Data/Streams/S1/Output/TEMP_OUT/MIXED C
T_S2 /Data/Streams/S2/Output/TEMP_OUT/MIXED C
HX1_Duty /Data/Blocks/HX1/Output/QCALC kW
T_S4 /Data/Streams/S4/Output/TEMP_OUT/MIXED C
T_S5 /Data/Streams/S5/Output/TEMP_OUT/MIXED C
HX2_Duty /Data/Blocks/HX2/Output/QCALC kW
T_S7 /Data/Streams/S7/Output/TEMP_OUT/MIXED C
T_S8 /Data/Streams/S8/Output/TEMP_OUT/MIXED C
HX3_Duty /Data/Blocks/HX3/Output/QCALC kW
T_S3 /Data/Streams/S3/Output/TEMP_OUT/MIXED C
R1 /Data/Blocks/R1/Output/QCALC kW
cond_ti /Data/Blocks/T1/Output/B_TEMP/2 C
cond_to /Data/Blocks/T1/Output/B_TEMP/1 C
cond_Q /Data/Blocks/T1/Output/COND_DUTY kW
reb_ti /Data/Blocks/T1/Output/B_TEMP/23 C
reb_to /Data/Blocks/T1/Output/B_TEMP/24 C
reb_Q /Data/Blocks/T1/Output/REB_DUTY kW
name value unit description
heater_1 6370.16251 kW
heater_5 2221.69628 kW
MSC_power_tot 1350.56296 kW

Problem Definition

Here we define the Energy Technology (ET) we want to solve.

This project will use the Methanol ET.

Biomass ET

This ET will use the following Layers

OSMOSE LAYERS Methanol
Layer Display name shortname Unit Color
ELEC Electricity elec kW yellow

The methanol ET contains the following units

OSMOSE UNIT Methanol
unit name type
Methanol Process

Methanol Unit

cost1 cost2 cinv1 cinv2 imp1 imp2 fmin fmax
0 0 0 0 0 0 1 1

Methanol Unit Streams

Note to Meire you can add the electricity resources streams here directly as inputs to the process

After importing the powers of your compressors and pumps in your Aspen model. you can use this ET to sum up everything and report your net electricity consumption.

layer direction value
ELEC in 21350.56296

Heat Streams

name Tin Tout Hin Hout DT min/2 alpha
heater_1 25.0 700.0 0 6370.16251 2.5 1
heater_3 25.0 33.7975644 0 76.0102859 2.5 1
heater_4 33.7975644 94.9629521 0 3443.99453 2.5 1
heater_5 94.9629521 700.0 0 2221.69628 2.5 1
cooler_1 700.02157 25.0 0 -7224.99235 2.5 1
cooler_2 700.02157 25.0 0 -2296.65435 2.5 1
MSC_cooler_1 207.329046 25.0 0 -448.438476 2.5 1
MSC_cooler_2 207.401247 25.0 0 -450.076111 2.5 1
HX_1 207.656056 250.0 0 106.222888 2.5 1
HX_2 237.559676 50.0 0 -2546.423 2.5 1
HX_3 52.5240658 250.0 0 1445.43097 2.5 1
R_1 237.559676 237.559676 0 -1968.45722 2.5 1
reb 76.8962779 90.7934626 0 1469.00502 2.5 1
cond 63.9201383 58.7161669 0 -1409.85994 2.5 1

Solver configuration chunck for student VMs

OSMOSE OPTIONS mathProg
Property Value Comments
language ampl
solver gurobi
solver_options {}

Here we solve the optimization problem.

Your objective function can either be MER (Minimum energy requirement) or TotalCost. To run in MER mode you need to specify only the ET you want to solve however to run in TotalCost mode, sufficient utilities must be supplied for the optimization problem to be solve-able.

When running in MER mode you can display the cc and gcc curves with the minimum energy requirements chunk however when solving for TotalCost you should report the CAPEX/OPEX chunk and display the icc and carnot plots to evaluate the utility integration.

Optimization Results

Energy Requirements

Minimum Heating Requirement: 65.26 kW 
Minimum Cooling Requirement: 1278 kW

Plots

Grand Composite Curve

Composite Curves