Computational Design of Intercooled Packing for CO 2 Absorbers Grigorios Panagakos, Research Engineer, Carnegie Mellon University, LRST-National Energy Technology Laboratory
Computational Design of Intercooled Packingfor CO2 Absorbers
Grigorios Panagakos, Research Engineer,Carnegie Mellon University, LRST-National Energy Technology Laboratory
Additively Manufactured Intensified Device for Enhanced Carbon Capture by Monoethanolamine
Background: Temperature Bulge of Liquid• Occurs as a result of the exothermic
reaction of CO2 with MEA
• The location of the temperature bulge depends on the solvent-to-gas ratio
• The T bulge as well as L/G can significantly affect the removal rate
• The location of where to provide the cooling is not obvious
(Kvaamsdal & Rochelle, 2008)
Rel
ativ
e lo
catio
n of
T_b
aldg
efro
m b
otto
m b
ed %
Tem
pera
ture
(K)
Rem
oval
rate
(%)
L/G (kg/kg)
L/G (kg/kg)
Motivation-New Subtask• Hypothesis
– Optimization of heat removal could promote more CO2 absorption but tradeoff between capitaland operating cost must be taken into consideration
• Background– Current capture equipment design: Decoupled unit operations with mass transfer and heat
transfer– Decoupled stages with external cooling
• Objective (3years)– Provide computational and theoretical tools to assist the development, prototyping and
validation of enhanced CO2 capture with intercooled/intensified devices– Develop new computational framework to
• Understand the effect of operating conditions and material properties on system performance• parameterize geometry and programmatically update CFD model
– Investigate the effect of geometry through detailed CFD modeling– Develop Framework to link CFD, process modeling and performance optimization
Absorber Model• Carbon Capture model developed
by the CCSI team1,2 formonoethanolamine (MEA).
• Model tested to be accurate at theNCCC and TCM
• Model includes mass transfer,heat transfer, vapor-liquidequilibrium, and chemistry of theMEA-H2O-CO2 system
Process Modeling
4
[1] A. Soares Chinen, J. C. Morgan, B. Omell, D. Bhattacharyya, C. Tong, and D. C. Miller, “Development of a Rigorous Modeling Framework for Solvent-Based CO2 Capture. 1. Hydraulic and Mass TransferModels and Their Uncertainty Quantification,” Ind. Eng. Chem. Res., vol. 57, no. 31, pp. 10448–10463, Aug. 2018.[2] J. C. Morgan et al., “Development of a Rigorous Modeling Framework for Solvent-Based CO2 Capture. Part 2: Steady-State Validation and Uncertainty Quantification with Pilot Plant Data,” Ind. Eng.Chem. Res., vol. 57, no. 31, pp. 10464–10481, Aug. 2018.[3] Debangsu Bhattacharyya1 and David C Miller , Post-combustion CO2 capture technologies — a review of processes for solvent-based and sorbent-based CO2 capture, Chemical Engineering 2017,17:78–92
3
Implementation of Embedded Cooling
5
20
30
40
50
60
70
80
0 0.2 0.4 0.6 0.8 1
Tem
pera
ture
(°C)
z/L
Solvent Gas Coolant
20
30
40
50
60
70
80
0 0.2 0.4 0.6 0.8 1
Tem
pera
ture
(°C)
z/L
Solvent Gas Coolant
20
30
40
50
60
70
80
0 0.2 0.4 0.6 0.8 1
Tem
pera
ture
(°C)
z/L
Solvent Gas Coolant
• Adaptable to multipleactive/inactive coolingsections
• Co-current or counter-current flow
3 Active Sections 5 Active Sections
2 Active Sections
J.C. Morgan et al., Development of a Rigorous Modeling Framework forSolvent-Based CO2 Capture. Part 2: Steady-State Validation and UncertaintyQuantification with Pilot Plant Data, Ind. Eng. Chem. Res. 2018, 57,10464−10481
• Investigating the trade-off between increase in absorption performance andthe increase in equipment size
• Geometry of shell and tube type heat exchanger used to relate heat transferarea to increase in absorber size
Preliminary Modeling Work II
6
0%
5%
10%
15%
20%
25%
0 2.5 5 7.5 10 15 20 30 500.0%
0.5%
1.0%
1.5%
2.0%
2.5%
Abso
rber
Size
Incr
ease
Specific Cooling Area (m2/m3)
Capt
ure
Incr
ease
Capture
Size
• Investigating the trade-off between increase in absorption performance andthe increase in equipment size
• Geometry of shell and tube type heat exchanger used to relate heat transferarea to increase in absorber size
Preliminary Modeling Work II
7
0%
5%
10%
15%
20%
25%
0 2.5 5 7.5 10 15 20 30 500.0%
0.5%
1.0%
1.5%
2.0%
2.5%
Abso
rber
Size
Incr
ease
Specific Cooling Area (m2/m3)
Capt
ure
Incr
ease
Capture
Size
How do you solve physics problems?
? ??
Fridolin Okkels and Henrik Bruus, Scaling behavior of optimally structured catalytic microfluidic reactors, Physical Review E 75, 016301 2007L.H. Olesen et al. A high-level programming-language implementation of topology optimization applied to steady-state Navier–Stokes flow, Int. J. Numer. Meth. Engng 2006; 65:975–1001
How do you solve physics problems?
? ??
Model Solver Solution
Fridolin Okkels and Henrik Bruus, Scaling behavior of optimally structured catalytic microfluidic reactors, Physical Review E 75, 016301 2007L.H. Olesen et al. A high-level programming-language implementation of topology optimization applied to steady-state Navier–Stokes flow, Int. J. Numer. Meth. Engng 2006; 65:975–1001
How do you solve physics problems?
? ??
Model Solver Solution
Governing Eq.+ Constitutive Laws+ Setup + BC
COMSOL u,v,p,c
Fridolin Okkels and Henrik Bruus, Scaling behavior of optimally structured catalytic microfluidic reactors, Physical Review E 75, 016301 2007L.H. Olesen et al. A high-level programming-language implementation of topology optimization applied to steady-state Navier–Stokes flow, Int. J. Numer. Meth. Engng 2006; 65:975–1001
How do you solve physics problems?
? ??
Model Solver Solution
Governing Eq.+ Constitutive Laws+ Setup + BC
COMSOL u,v,p,c
e.g. NS + CDR+ BC + Goal TOPOPT Best design
Fridolin Okkels and Henrik Bruus, Scaling behavior of optimally structured catalytic microfluidic reactors, Physical Review E 75, 016301 2007L.H. Olesen et al. A high-level programming-language implementation of topology optimization applied to steady-state Navier–Stokes flow, Int. J. Numer. Meth. Engng 2006; 65:975–1001
How do you solve physics problems?
? ??
Model Solver Solution
Governing Eq.+ Constitutive Laws+ Setup + BC
COMSOL u,v,p,c
e.g. NS + CDR+ BC + Goal TOPOPT Best design
Fridolin Okkels and Henrik Bruus, Scaling behavior of optimally structured catalytic microfluidic reactors, Physical Review E 75, 016301 2007L.H. Olesen et al. A high-level programming-language implementation of topology optimization applied to steady-state Navier–Stokes flow, Int. J. Numer. Meth. Engng 2006; 65:975–1001
Fundamentals of Topology Optimization
M.P.Bendsøe and O. Sigmund, Topology Optimization - Theory, Methods and Applications, Springer (2003)
γ=0, Solidγ=1, Void
Topology Optimization in Aircraft Design
D. Walker et al. Topology Optimization of an Aircraft Wing, 56th AIAA/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference
“…Airbus researches use of topology optimization on aircraft wing ribs. It is stated that usage of topology optimization results in around 1000 kg of weight savings per aircraft…”
https://topologyoptimization.wordpress.com/2011/03/11/airbus/
• Optimization Problem = Inverse Problem
Fundamentals of Optimization Problems
Tune the design-variables (γ) of a problem, such that an objective function (Φ) is minimized, under given constrains.
• Design variables: →
• Objective function:
• Constrains:
u[γ]γ →
Φ
Mathematical Formulation of the PDE problem
Dirichlet b.c.
Neumann b.c.
𝚪𝚪4 ≡−𝐽𝐽1−𝐽𝐽2
𝐹𝐹4 ≡ 𝒗𝒗 � 𝛻𝛻𝑐𝑐 + 𝑅𝑅
We can use the software package COMSOL to solve the reacting flow problem.Requires Partial Differential Equations to be expressed in divergence form
Then the reacting flow problem can be written as:
The Finite Element Method (FEM)
• Expand both u and γ on a finite basis set
• Use P1 basis for p(r), c(r) and γ(r)
• Use P2 basis for u(r)
Finite element mesh; basis function ϕj has small support
• Minimize objective function Φ(u) subject to constraints
Constrained optimization problem
Setup and procedure for micro-reactor
Setup and procedure for micro-reactor
COMSOL GUI
Modeling
Setup and procedure for micro-reactor
COMSOL GUI
Modeling
COMSOL Script
MMA.m(MATLAB)
Adj-Pr.m
TOPOPT.m
Method
L.H. Olesen, F. Okkels, og H. Bruus, Int. J. Num. Meth. Eng. 65, 975 (2006)
Setup and procedure for micro-reactor
COMSOL GUI
Modeling
COMSOL Script
MMA.m(MATLAB)
Adj-Pr.m
TOPOPT.m
Method
L.H. Olesen, F. Okkels, og H. Bruus, Int. J. Num. Meth. Eng. 65, 975 (2006)
COMSOL ScriptCOMSOL/MATLAB
Visualization/Analysis
General export
• Characteristic Examples
Topology Optimized Micro-Reactors
𝐃𝐃𝐃𝐃 = 𝟏𝟏𝟏𝟏−𝟓𝟓,𝐃𝐃 = 𝟏𝟏 � 𝟏𝟏𝟏𝟏−𝟖𝟖𝐦𝐦𝟐𝟐
𝐬𝐬 ,𝐝𝐝𝐝𝐝 = 𝟏𝟏.𝟓𝟓 𝐏𝐏𝐃𝐃,𝐤𝐤𝐃𝐃 = 𝟏𝟏𝐃𝐃𝐃𝐃 = 𝟏𝟏𝟏𝟏−𝟓𝟓,𝐃𝐃 = 𝟏𝟏 � 𝟏𝟏𝟏𝟏−𝟖𝟖𝐦𝐦𝟐𝟐
𝐬𝐬 ,𝐝𝐝𝐝𝐝 = 𝟏𝟏.𝟏𝟏 𝐏𝐏𝐃𝐃,𝐤𝐤𝐃𝐃 = 𝟏𝟏
Iter=227Iter=188
• Characteristic Examples
Topology Optimized Micro-Reactors
𝐃𝐃𝐃𝐃 = 𝟏𝟏𝟏𝟏−𝟓𝟓,𝐃𝐃 = 𝟏𝟏 � 𝟏𝟏𝟏𝟏−𝟖𝟖𝐦𝐦𝟐𝟐
𝐬𝐬 ,𝐝𝐝𝐝𝐝 = 𝟏𝟏.𝟓𝟓 𝐏𝐏𝐃𝐃,𝐤𝐤𝐃𝐃 = 𝟏𝟏𝐃𝐃𝐃𝐃 = 𝟏𝟏𝟏𝟏−𝟓𝟓,𝐃𝐃 = 𝟏𝟏 � 𝟏𝟏𝟏𝟏−𝟖𝟖𝐦𝐦𝟐𝟐
𝐬𝐬 ,𝐝𝐝𝐝𝐝 = 𝟏𝟏.𝟓𝟓 𝐏𝐏𝐃𝐃,𝐤𝐤𝐃𝐃 = 𝟏𝟏.𝟐𝟐𝟓𝟓
Iter=227Iter=121
NETLLocal Information- High fidelity
CFD
ORNLExperimentation, manufacturing and testing
WVUProcess modeling
Performance Optimization
Models and numericalcampaign to informexperimental work
Experimental campaign tocollect data to improvemodels
Partnership with ORNL and WVU
25
Moving Forward
26
Improve processmodel in terms ofCoolant model andimproved discreteor continuousdesign of coolinglocations
Link CFD withProcess Modeling
Study Processeconomics
MEA+Intensifiedreactor Goal: improve packing with in situ continuous cooling (absorber) and heating (stripper)
Develop CFD framework to simulate Multiphase flow Chemical reactions Thermal coupling
between chemistry and fluid flow on limiting mechanisms
Packing Fluid properties Operating conditions
MEA+Mellapakequivalent, to match literature performance
Calibrate modeland validatecomputational andtheoretical workagainst actual 3Dprinted structuredpacking
Acknowledgements
27
Carbon Capture Simulation for Industry Impact (CCSI2)
ORNL: Xin Sun, Costas Tsouris, Lonie Love, Canhai Lai
NETL: Michael Matuszewski, Benjamin Omell
WVU: Debangsu Bhattacharyya, Ryan Hughes, Goutham Kotamreddy
Disclaimer: This work is made available as an account of work sponsored by an agency of the United States Government. Neither the UnitedStates Government nor any agency thereof, nor any of their employees, makes any warranty, express or implied, or assumes any legal liability orresponsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its usewould not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark,manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United StatesGovernment or any agency thereof. The views and opinions of authors expressed herein do not necessarily state or reflect those of the UnitedStates Government or any agency thereof.
For more informationhttps://www.acceleratecarboncapture.org/
Grigorios Panagakos, Research EngineerCarnegie Mellon University-LRST/[email protected] [email protected]