Advancements in Carbon Dioxide and Water Vapor Separations Using COMSOL
J. Knox1 K. Kittredge1 R.F. Coker1 R. Cummings1 C.F. Gomez1
1NASA - Marshall Space Flight Center, Huntsville, AL, USA
“NASA's Advanced Exploration Systems (AES) program is pioneering new approaches for rapidly developing prototype systems, demonstrating key capabilities, and validating operational concepts for future human missions beyond Earth orbit” . Under the new Atmosphere Revitalization Recovery and Environmental Monitoring (ARREM) project , efforts are focused on improving current state-of-the-art systems utilizing fixed beds of sorbent pellets by evaluating structured sorbents, seeking more robust pelletized sorbents, and examining alternate bed configurations to improve system efficiency and reliability. These development efforts combine testing of sub-scale systems and multi-physics computer simulations to evaluate candidate approaches, select the best performing options, and optimize the configuration of the selected approach, which is then implemented in a full-scale integrated atmosphere revitalization test. Reference 3 discusses the hardware design and sorbent screening and characterization effort in support of the ARREM project within the AES program. This paper describes the development of atmosphere revitalization models and simulations in support of the ARREM project. COMSOL Multiphysics has been used in the following technology development efforts: (1) Development of 1-D and 2-D Axisymmetric Packed Bed Models and Comparison with Breakthrough Curves. Adsorption in fixed beds of pelletized sorbents is the primary means of gas separation for atmosphere revitalization systems. For the bulk separation of CO2 and H2O, temperature changes due to the heat of adsorption are significant, requiring the modeling and simulation of the heat balance equations. For columns with small tube diameter to pellet diameter ratios, as encountered in internally heated columns, flow channeling along the column wall can have a strong influence on overall performance . In some cases, the influence is great enough to necessitate the use of 2-D simulations as shown in Figure 1. (2) Isothermal Bulk Desiccant Development. Initial models of the 4-column IBD with an aluminum (Al) foam lattice were conducted using the COMSOL Multiphysics code. For the porous flow part of the model, the code solves the low Mach number, compressible Brinkman equations modified by a Forchheimer drag term. Forchheimer drag adds a drag term that is proportional to the square of the fluid velocity, rather than just linear, as is the case in Darcy flow, with a coefficient given by βF = 1.75ρg*√(ε /150/κ). Heat transfer is modeled by solving the heat transfer equations separately for the porous medium (the beds) and the solid Al canister, with appropriate boundary conditions. Modeling results are showing in Figure 2. A second prototype test fixture, also shown in Figure 2, will be modeled following model validation. (3) Microlith Residual Desiccant Development. Figure 3 provides a cross-section of the Microlith concept and 3-D COMSOL velocity mapping of the jelly roll. This analysis provided important insight leading to successful corrective action for a flow channeling problem. Further uses of this model will be to map flow maldistribution following corrective action. With the addition of adsorption physics, this model will provide a means for optimization of cyclic parameters for this hardware, and allow for design optimization studies for new Microlith designs.