CGAM Process¶

This tutorial demonstrates how to perform an exergy analysis of a CGAM process using ExerPy. The analysis is carried out with:

Ebsilon
Aspen Plus
tespy

The flowsheet of the system is shown below.

CGAM flowsheet — Figure: Flowsheet of the CGAM process model¶

Note

The tespy based flowsheet has slight deviations from the show flowsheet, due to not using the identical component definitions. For example, Ebsilon provides a component incorporating drum, recirculation pump and evaporator in a single part, while these are individual components in tespy.

The import of the exerpy dependency is the same for all simulators:

Import Necessary Modules

from exerpy import ExergyAnalysis

Ebsilon

Download the Ebsilon simulation model here: cgam.ebs

Initialize the Exergy Analysis

Create an instance of the ExergyAnalysis class using the from_ebsilon method. Since there is a combustion process in the power plant, it is necessary to specify the chemical exergy library to use. In this case, we will use the Ahrendts library. You can decided weather to split the physical exergy into its mechanical and thermal shares by setting the split_physical_exergy parameter to True or False.

model_path = 'cgam.ebs'

ean = ExergyAnalysis.from_ebsilon(model_path, chemExLib='Ahrendts', split_physical_exergy=False)

Define the exergy flows crossing the system boundaries

ExerPy requires you to specify the fuel (E_F), product (E_P), and loss (E_L) exergy streams in your system. These are defined using dictionaries with "inputs" and "outputs" keys, containing lists of connection IDs from your Ebsilon model.

fuel = {
    "inputs": ['1', '10'],
    "outputs": []
}

product = {
    "inputs": ['E1', '9'],
    "outputs": ['8']
}

loss = {
    "inputs": ['7'],
    "outputs": []
}

For this process, the exergetic fuel of the system (E_F) is the methane flow (10) and air flow (1). The exergetic product of the system (E_P) consists of the net power output (E_1) and the difference between the exergy of the vapour outlet (9) and of the water inlet (8). The exhaust stream (7) represents the exergy loss of the system (E_L).

Note

The dictionary labels must exactly match the connection labels defined in the Ebsilon simulation model. By default, Ebsilon assigns generic names to streams (e.g., “water”, “water_1”). It is strongly recommended to rename the stream labels in Ebsilon using consistent and meaningful labels.

For example: use “1”, “2”, “3” for material connections and “E1”, “E2”, “E3” for electrical connections, “H1”, “H2”, “H3” for heat connections, etc.

tespy

For the tespy model we have prepared the code to run the simulation in the dropdown below. To learn how to set up tespy models and what things to be aware of when working with tespy, we kindly refer to the online documentation of tespy.

Initialize the Exergy Analysis

After setting up the model, we set up the ExergyAnalysis instances using the from_tespy method. It takes the converged tespy.Network object along with ambient state and (optionally) the chemical exergy library as inputs.

Tip

TESPy can handle the splitting of physical exergy into its mechanical and thermal shares, therefore split_phyiscal_exergy can always be set to True when using tespy. In this instance it is set to False because ASPEN cannot handle this, and we wanted to cross validate the results of the examples for all three simulators.

Define the exergy flows crossing the system boundaries

For this power plant, the exergetic fuel of the system (E_F) is the methane flow (10) and air flow (1), which are the inputs to the combustion process. The exhaust stream (7) represents the exergy losses of the system (E_L). The product exergy (E_P) is the power output of the turbine minus the power requirement of the compressor e3 as well as the difference in total exergy between the flow of water entering the economizer 8`and the saturated steam leaving the drum :code:`9.

Aspen Plus

Download the Aspen simulation model here: cgam.bkp

Initialize the Exergy Analysis

Create an instance of the ExergyAnalysis class using the from_aspen method. Since there is a combustion process in the power plant, it is necessary to specify the chemical exergy library to use. In this case, we will use the Ahrendts library.

Note

At the moment, it is not possible to split the physical exergy into its mechanical and thermal shares when using Aspen Plus. Therefore, the split_physical_exergy parameter should be always set to False when using the from_aspen method.

model_path = 'cgam.bkp'

ean = ExergyAnalysis.from_aspen(model_path, chemExLib='Ahrendts', split_physical_exergy=False)

Define the exergy flows crossing the system boundaries

fuel = {
    "inputs": ['1', '10'],
    "outputs": []
}

product = {
    "inputs": ['E1', '9'],
    "outputs": ['8']
}

loss = {
    "inputs": ['7'],
    "outputs": []
}

Note

The dictionary labels must exactly match the connection labels defined in the Aspen simulation model. It is strongly recommended to rename the stream labels in Aspen using consistent and meaningful labels.

For example: use “1”, “2”, “3” for material connections and “E1”, “E2”, “E3” for electrical connections, “H1”, “H2”, “H3” for heat connections, etc.

Running the exergy analysis and working with the results is now independant for all simulators.

Perform the Exergy Analysis

Run the analysis by invoking the analyse method on the ExergyAnalysis instance, passing the defined fuel, product, and loss exergy streams.

ean.analyse(E_F=fuel, E_P=product, E_L=loss)

Retrieve and Display Results

After the analysis is complete, retrieve the results using the exergy_results method.

# Retrieve and display the results
df_components, df_material_connections, df_non_material_connections = ean.exergy_results()

# Print the components exergy results
print(df_components)

# Optionally, save the results to CSV files
df_components.to_csv('components_exergy_results.csv')
df_material_connections.to_csv('material_connections_exergy_results.csv')
df_non_material_connections.to_csv('non_material_connections_exergy_results.csv')