Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic Modelling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Electron heating in kinetic-Alfvén-wave turbulence

We report analytical and numerical investigations of subion-scale turbulence in low-beta plasmas using a rigorous reduced kinetic model. We show that efficient electron heating occurs and is primarily due to Landau damping of kinetic Alfvén waves, as opposed to Ohmic dissipation. This collisionless damping is facilitated by the local weakening of advective nonlinearities and the ensuing unimpeded phase mixing near intermittent current sheets, where free energy concentrates. The linearly damped energy of electromagnetic fluctuations at each scale explains the steepening of their energy spectrum with respect to a fluid model where such damping is excluded (i.e., a model that imposes an isothermal electron closure). In conclusion, the use of a Hermite polynomial representation to express the velocity-space dependence of the electron distribution function enables us to obtain an analytical, lowest-order solution for the Hermite moments of the distribution, which is borne out by numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Ignition delay time and H 2 O measurements during methanol oxidation behind reflected shock waves

To improve detailed chemical kinetics models, the oxidation of methanol was investigated behind reflected shock waves in shock tubes. Ignition delay times of methanol–air mixtures, with Ar as diluent, were studied between 940 and 1540 K in a heated shock tube, for pressures up to 14.9 atm and for equivalence ratios of 0.5, 1.0, and 2.0. Water profiles were measured by utilizing a laser absorption technique in the 1350-to-1600-K temperature range, at an average pressure of 1.3 atm and for similar equivalence ratios. The present study shows here the ignition delay times of methanol to be in very good agreement with results from the literature (Fieweger et al., 1997), whereas the other conditions have never been investigated before. The ignition delay time data are also in good agreement with modern detailed kinetics mechanisms such as the AramcoMech 3.0 model. The water time-history profiles were modeled using well-known literature mechanisms. Discrepancies were observed between these kinetics mechanisms, and poor predictions were observed for the lower temperatures investigated. Sensitivity and rate-of-production analyses were performed using 3 literature mechanisms (namely, AramcoMech 3.0, Princeton, and JetSurfII). Discrepancies were found among the models when predicting important reactions dominating the oxidation of methanol as well as the rate-of-production of H 2 O.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Experimental and computational investigation of the influence of iso -butanol on autoignition of n -decane and n -heptane in non-premixed flows

Experimental and computational investigation is carried out to elucidate the influence of iso-butanol on critical conditions of autoignition of n-decane and n-heptane. The counterflow configuration is employed. In this configuration an axisymmetric stream of air is directed over the surface of an evaporating pool of a liquid fuel. A stagnation plane is established. The temperature of the air is increased until autoignition is achieved in the mixing layer in the vicinity of the stagnation plane. The temperature of the air stream at autoignition, T ig , is measured at various values of the strain rate, which is defined as the axial gradient of the axial component of the flow velocity at the stagnation plane. Fuels tested are n-decane, n-heptane, iso-butanol and various mixtures of n-decane/iso-butanol and n-heptane/iso-butanol. Kinetic modeling is carried out using the recently updated version of the comprehensive CRECK chemical-kinetic mechanism. Critical conditions of autoignition are predicted for all fuels and fuel mixtures and the results are compared with the measurements. Computations show that low-temperature chemistry plays a significant role in promoting autoignition of n-decane and n-heptane, and the influence of low-temperature chemistry decreases with increasing strain rates because there is insufficient residence for the low-temperature reactions to take place. Experimental data and numerical simulations show that addition of even small amounts of iso-butanol to n-decane or n-heptane increases the value of T ig at low strain rates indicating that iso-butanol strongly inhibits the low-temperature chemistry of n-decane and n-heptane. Furthermore, the simulations show that when iso-butanol is present in the liquid fuel, only a small amount of n-decane is available in the gas phase and its low concentration cannot effectively sustain the low temperature degenerate branching oxidation process. Here, predicted flame structures show that the peak values of mole fraction of ketohydroperoxide are significantly reduced when iso-butanol is added to n-decane indicating that the kinetic pathway to low temperature ignition is blocked. This observation is confirmed from sensitivity analysis.

42 ENGINEERING↗
DOE OSTI2023

RAM-SCB Runs on Request at CCMC

The Ring current–Atmosphere interactions Model (RAM) with Self-Consistent magnetic (B) field (SCB) combines a large-scale kinetic model of ring current plasma with a three-dimensional (3-D) force-balanced model of the terrestrial magnetic field. RAM-SCB simulates the evolution of major ion species (H + , O + , and He + by default) and electrons as a function of azimuth, radial distance, energy, and pitch angle. Simulation outputs include the Dst index, and pressure and differential flux of the modeled species, thus providing benefit as a science code and to inform surface charging hazard within the model domain. Version 2.2 of this open-source simulation code has now—as of 8 August 2023—been made available for Runs-On-Request via the Community Coordinated Modeling Center.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Enhanced reactivity at the oil–water interface accelerates the synthesis of zymonic acid in microemulsions

Chemical reactions in microscale compartments, such as aerosols and emulsions, can exhibit significantly faster reaction rates relative to macroscale containers. This enhancement in chemistry is often due to the elevated importance of surfaces as reaction vessels are reduced to picoliter volumes. While most studies have focused on the air-water interface of droplets, there are comparably fewer studies of reactions in micron-scale aqueous solutions encapsulated by oil. Here we investigate the condensation reaction of pyruvic acid (PA) to form zymonic acid (ZA) and water. Using microfluidics and optical trapping, chemical kinetics are measured in monodisperse micron-sized emulsion droplets in situ via Raman spectroscopy. Relative to a macroscopic bulk solution, which exhibits little to no reaction over many days, we find efficient production of ZA over the same time period. A kinetic model is developed to elucidate the role of the interface in accelerating the microdroplet reaction kinetics. After quantifying the surface partitioning of PA from interfacial tension measurements, the rate coefficient for the condensation reaction at the oil-water interface is determined to be 1.8 × 10-2 M-1 s-1. This rate coefficient is estimated to be 105 larger than the reaction rate in bulk aqueous solutions. Compared to previous studies of accelerated ZA formation at the air-water interface on nanodroplets, we find that the reaction at the oil-water interface is 20 times more efficient. Despite this difference, the overall ZA formation rate in emulsions is significantly slower than in the same-sized aerosols, which arises from the weaker partitioning of PA to the oil-aqueous relative to air-water interface. These results highlight the interplay between interfacial partitioning and reactivity in accelerating chemistry in microcompartments and provides new insights into how interfacial composition influences condensation reactions.

Wallace, Brandon J↗
DOE OSTI2021

Modeling of Glycolate Destruction in the Recycle Collection Tank

The Savannah River Site’s DWPF is being upgraded with the introduction of the NG flowsheet. Glycolic acid has been shown to be a superior alternative to formic acid for sludge processing. The new flowsheet improves or maintains necessary parameters such as 1) reduction of mercury, 2) adjustment of feed rheology, 3) pH stability, and 4) adjustment of melter oxidation/reduction potential. Further, the use of glycolic acid virtually eliminates the potential for catalytic hydrogen generation in DWPF processing DWPF process condensates are collected and returned to the Savannah River Site (SRS) CSTF. The RCT collects off-gas condensate during chemical processing, vitrification, and other unit operations performed in DWPF and is the singular return vessel delivering recycle effluent back to CSTF. Each batch of recycle will have a small amount of glycolate from chemical processing and melter off-gas condensates. To avoid potential flammability issues due to thermolysis of glycolate in the CSTF, Savannah River National Laboratory (SRNL) provided to Savannah River Remediation (SRR) at their request a Task Technical and Quality Assurance Plan (TTQAP) to quantify and mitigate glycolate returns via DWPF’s recycle stream. The request included testing of a process to oxidize glycolate and other organic species that are responsible for hydrogen generation from thermolysis. Following that work SRR provided a Task Technical Request (TTR) that requested process modeling. In 2021 a TTQAP was issued to cover the modeling work. Modeling draws data from laboratory scale studies using chemical simulants and radioactive waste samples. Chemical kinetic modeling was performed to evaluate the feasibility of using sodium permanganate to destroy glycolate in the RCT. The results from the laboratory studies were summarized in a series of reports. Reference 9 is a report of lab scale processing of actual DWPF Slurry Mix Evaporate Condensate Tank (SMECT) and Offgas Condensate Tank (OGCT) samples in the SRNL Shielded Cells. Tests at caustic conditions demonstrated sodium permanganate was effective in converting glycolate to oxalate, and permanganate (Mn 7+ ) is reduced to manganate (Mn 6+ ) with no significant formation of carbon dioxide or carbonate. Equation (1) was found to best describe the observed reaction of glycolate with permanganate under nominal (60 to 145 mg/L in RCT) glycolate entrainment conditions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2019

Impact of Ethanol additions on Autoignition Characteristics of a Full Boiling Range Gasoline and its Surrogates at Advanced Engine Conditions

Ethanol is widely blended into commercial gasolines, as an oxygenate to mitigate pollutant formation, as well as an octane enhancer to improve knock resistance. Ethanol can also facilitate the use of lower-grade, or less-refined gasoline in modern internal combustion (IC) engines. However, the chemical kinetic interactions of ethanol with petroleum-based gasoline are not well understood, especially towards autoignition propensity over a range of conditions. These interactions are critical for improving knock resistance, as well as towards performance in low temperature combustion schemes, such as gasoline compression ignition. This study investigates the chemical kinetic interactions of ethanol with a full boiling range gasoline, where blending ratios of 0 to 30 vol./vol. are covered. A mid-octane range (Anti Knock Index 91.5) Fuels for Advanced Combustion Engines gasoline is employed, while a variety of multi-component surrogate blends are also studied. Autoignition data are acquired within a rapid compression machine at temperature, pressure, and fuel loading conditions representative of advanced IC engines. Changes in low-temperature behavior, including first-stage ignition times and extents of low temperature heat release, as well as main ignition times are quantified for the neat fuels and the ethanol-blended fuels. A recently-updated, detailed chemical kinetic model for the gasoline surrogates is used to interpret the measurements, and better understand the kinetic pathways.

autoignition↗
DOE OSTI2025

Integrated Molten Salt Reactor Modeling Capabilities in NEAMS Thermal Hydraulics Tools

The DOE neams program supports a full range of computational thermal fluids analysis capabilities and code developments for a broad range of advanced reactor concepts. The research and development approach under the thermal fluids technical area synergistically combines three length and time scales in a hierarchical multi-scale approach. To enable multi-scale thermal fluids capability using these codes, a key joint effort has been underway to develop an integrated system- and engineering-scale thermal fluids analysis capability, through integration of SAM and Pronghorn codes, both based on the MOOSE framework. This report summarizes recent advances in developing an integrated system- and engineering-scale modeling capability for the msr concept, which has gained significant interest in recent years. A consistent framework was established by coupling Pronghorn and SAM through the Saline interface, with thermophysical properties provided by the Molten Salt Thermal Property Database (MSTDB-TP). Further improvements were made to the coupling schemes and domain-overlapping strategies, enhancing the stability and robustness of multi-code simulations. Verification and validation efforts demonstrate the accuracy of this integration across a range of benchmark problems, including one-dimensional heated pipe flows, three-dimensional natural convection loops with evolving isotopic compositions, and \gls{msre} demonstration cases. Within Pronghorn, new capabilities were introduced to model corrosion and noble-metal plating phenomena, supported by an extended thermal-hydraulics framework and refined turbulence treatments. To capture two-phase flow behavior, a multiphase Euler–Euler model was implemented in Pronghorn, including advanced closure relations, high-resolution advection techniques, and capillary force reconstruction. Preliminary verification cases confirm the fidelity of the approach, while planned validation efforts target canonical multiphase benchmarks and application to msr components such as the msre pump bowl. Finally, updates to SAM’s msr mass transfer modeling were extended to consider noble gas migration into porous structures like graphite. The point kinetics model was updated to include reactivity feedback contributions from any defined species, such as xenon. The gas transport model was expanded for applicability to gas mixtures, bubble efflux phenomena, and species transport between liquid and gas phases. A selection of multi-scale Sherwood number correlations from MOSCATO/NekRS and multi-phase correlations from literature have been added for improved accuracy in calculating mass transfer coefficients. A companion effort on developing system-level redox corrosion has also been incorporated into SAM. Collectively, these enhancements strengthen the predictive capability of SAM and Pronghorn for simulating MSR thermal-hydraulics, corrosion, multiphase behavior, and fission-product transport, providing a more complete toolset for design, safety analysis, and licensing support of next-generation \gls{msr}s.

42 - ENGINEERING↗
DOE OSTI2018

Fractional analog of a chemical system inspired by Braess’ paradox

In this paper, we introduce the fractional analog of a chemical model arouse from a mathematical paradox attributed to Dietrich Braess. Two basic examples which serve fractional kinetic models as better suited models to the real data sets than the integer-order counterparts are given. Existence and uniqueness of the rebuilt model’s solutions are proved. It is shown that asymptotic stability conditions of the solutions are provided. A comparison is made between two different solution methods and numerical simulations are also presented to exemplify the mathematical outcomes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

An evaluation of the efficacy of various coal combustion models for predicting char burnout

Coal combustion is comprised of several subprocesses including devolatilization and heterogeneous reactions of the coal char with O 2 , CO 2 , H 2 O and potentially several other species. Much effort has been put forth to develop models for these processes which vary widely in both complexity and computational cost. This work investigates the efficacy of models for devolatilization and char reactions at either end of the complexity and cost spectrums for a range of particle sizes and furnace temperatures and across coal types. The overlap of simulated devolatilization and char consumption is also examined. In the gas phase, a detailed kinetics model based on a reduced version of the GRI 3.0 mechanism is used. The Char Conversion Kinetics and an n th -order Langmuir-Hinshelwood models are considered for char oxidation. The Chemical Percolation and Devolatilization and a two-step model are considered for devolatilization. Results indicate that high-fidelity models perform better at representing particle temperature and mass data across a wide range of O 2 concentrations as well as coal types. A significant overlap in devolatilization and char consumption is observed for both char chemistry and devolatilization models.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2015

Interstellar hydrogen fluxes measured by IBEX-LO in 2009: Numerical modeling and comparison with the data

In this paper, we perform numerical modeling of the interstellar hydrogen fluxes measured by IBEX-Lo during orbit 23 (spring 2009) using a state-of-the-art kinetic model of the interstellar neutral hydrogen distribution in the heliosphere. This model takes into account the temporal and heliolatitudinal variations of the solar parameters as well as the non-Maxwellian kinetic properties of the hydrogen distribution due to charge exchange in the heliospheric interface. We found that there is a qualitative difference between the IBEX-Lo data and the modeling results obtained with the three-dimensional, time-dependent model. Namely, the model predicts a larger count rate in energy bin 2 (20–41 eV) than in energy bin 1 (11–21 eV), while the data shows the opposite case. We perform study of the model parameter effects on the IBEX-Lo fluxes and the ratio of fluxes in two energy channels. We show that the most important parameter, which has a major influence on the ratio of the fluxes in the two energy bins, is the solar radiation pressure. The parameter fitting procedure shows that the best agreement between the model result and the data occurs when the ratio of the solar radiation pressure to the solar gravitation, μ{sub 0}, is 1.26{sub −0.076}{sup +0.06}, and the total ionization rate of hydrogen at 1 AU is β{sub E,0}=3.7{sub −0.35}{sup +0.39}×10{sup −7} s{sup −1}. We have found that the value of μ{sub 0} is much larger than μ{sub 0} = 0.89, which is the value derived from the integrated solar Lyα flux data for the period of time studied. We discuss possible reasons for the differences.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Drift-kinetic effects of tungsten on plasma response to RMP in ITER

Here, effects of high- Z ( Z is the particle charge number) tungsten impurity ions on the plasma response to the resonant magnetic perturbation (RMP) field are numerically investigated for the ITER 15 MA baseline scenario, where the tungsten contribution to the plasma response is computed with a drift-kinetic model while the bulk thermal particle contributions follow the fluid approximation. The study yields three highlights: (i) the drift-kinetic contribution of the tungsten impurity exerts minor influence on the plasma response compared to that computed by the pure fluid model without tungsten; (ii) a new figure of merit, based on the resonant spectrum perturbation at the plasma boundary surface, results in different optimal coil phasing compared to that previously obtained by maximizing the edge-peeling plasma response; (iii) the optimal $n = 3$ RMP (for edge localized mode (ELM) control, $n$ is the toroidal mode number) is found to induce a large tungsten particle influx near the plasma edge associated with the neoclassical toroidal viscosity. The study thus provides useful data on the compatibility of the full tungsten wall with RMP ELM control in ITER.

ITER↗
DOE OSTI2018

Reduced chemical model for low and high-temperature oxidation of fuel blends relevant to internal combustion engines

A hybrid approach is proposed to develop reduced kinetic models for complex engine-relevant fuels. The reduced mechanism is composed of a core submechanism including C0-C4 chemistry, ethanol chemistry, and NOx chemistry and a fuel-dependent submechanism. The fuel-dependent submechanism consists of three species and ten reactions describing both the low and high-temperature fuel decomposition pathways. Calibrations for these reactions can be made for single component and multicomponent mixtures using experimental targets or results from detailed kinetic mechanisms. In the present study, the methodology is applied to the combustion of gasoline surrogates. Reduced mechanisms are calibrated for primary references fuels and multicomponent mixtures at engine-relevant pressures. The reduced mechanisms capture the low-temperature heat release and negative temperature coefficient (NTC) behavior, which are important to simulations of internal combustion engine performance. The mechanisms are very compact in size and can be easily calibrated for multicomponent mixtures without any detailed knowledge of the chemistry of the fuel components, making them well suited for CFD simulations of new fuel blends.

02 PETROLEUM↗
DOE OSTI2026

Rate constants and third-body collision efficiencies for recombination of Na with OH and O2: Implications for flame inhibition by alkali salts

Flame inhibition by alkali metals has implications for solid fuel combustion, fire safety, and a number of industrial processes. While the mechanism of inhibition is fairly well understood, details related to thermodynamic properties and rate constants are still in question. In the present work, the recombination of Na with OH (R5) and O2 (R7), respectively, was characterized theoretically, and the implications for modeling laminar premixed hydrogen flames doped with sodium species were examined. Third-body collision efficiencies and low-pressure-limit rate constants were obtained using newly fitted ab initio-based potential energy surfaces, classical trajectories, and one-dimensional master equation calculations. The results are consistent with available experimental results and aid in the present modeling study by providing rate information for bath gases and conditions that remain unexplored experimentally. Chemical kinetic modeling of relative Na and absolute H and OH profiles in H2-fueled laminar, premixed flames doped with a sodium salt shows that the most important radical removal cycle is the sequence Na + OH (+M) → NaOH (+M) (R5), NaOH + H → Na + H2O (R12), even under oxidizing conditions. A secondary cycle, Na + O2 (+M) ⇄ NaO2 (+M) (R7), NaO2 + OH → NaOH + O2 (R14), is less important due to the low thermal stability of NaO2. In most flames, reaction R7 is partially equilibrated, and reaction R14 becomes rate-limiting for the second cycle. The flame analysis supports a lower value of k14 than indicated by recent work on KO2 + OH, but more work is required to confirm this.

Jasper, Ahren W.↗
DOE OSTI2016

Intrinsic selectivity and structure sensitivity of Rhodium catalysts for C 2+ oxygenate production [On the intrinsic selectivity and structure sensitivity of Rhodium catalysts for C 2+ oxygenate production]

Synthesis gas (CO + H 2 ) conversion is a promising route to converting coal, natural gas, or biomass into synthetic liquid fuels. Rhodium has long been studied as it is the only elemental catalyst that has demonstrated selectivity to ethanol and other C 2+ oxygenates. However, the fundamentals of syngas conversion over rhodium are still debated. In this work a microkinetic model is developed for conversion of CO and H 2 into methane, ethanol, and acetaldehyde on the Rh (211) and (111) surfaces, chosen to describe steps and close-packed facets on catalyst particles. The model is based on DFT calculations using the BEEF-vdW functional. The mean-field kinetic model includes lateral adsorbate–adsorbate interactions, and the BEEF-vdW error estimation ensemble is used to propagate error from the DFT calculations to the predicted rates. The model shows the Rh(211) surface to be ~6 orders of magnitude more active than the Rh(111) surface, but highly selective toward methane, while the Rh(111) surface is intrinsically selective toward acetaldehyde. A variety of Rh/SiO 2 catalysts are synthesized, tested for catalytic oxygenate production, and characterized using TEM. The experimental results indicate that the Rh(111) surface is intrinsically selective toward acetaldehyde, and a strong inverse correlation between catalytic activity and oxygenate selectivity is observed. Furthermore, iron impurities are shown to play a key role in modulating the selectivity of Rh/SiO 2 catalysts toward ethanol. The experimental observations are consistent with the structure-sensitivity predicted from theory. As a result, this work provides an improved atomic-scale understanding and new insight into the mechanism, active site, and intrinsic selectivity of syngas conversion over rhodium catalysts and may also guide rational design of alloy catalysts made from more abundant elements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Stochasticity and positive feedback enable enzyme kinetics at the membrane to sense reaction size

Here, we present detailed kinetic analyses of a panel of soluble lipid kinases and phosphatases, as well as Ras activating proteins, acting on their respective membrane surface substrates. The results reveal that the mean catalytic rate of such interfacial enzymes can exhibit a strong dependence on the size of the reaction system—in this case membrane area. Experimental measurements and kinetic modeling reveal how stochastic effects stemming from low molecular copy numbers of the enzymes alter reaction kinetics based on mechanistic characteristics of the enzyme, such as positive feedback. For the competitive enzymatic cycles studied here, the final product—consisting of a specific lipid composition or Ras activity state—depends on the size of the reaction system. Furthermore, we demonstrate how these reaction size dependencies can be controlled by engineering feedback mechanisms into the enzymes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Investigation of the decolorization efficiency of two pin-to-plate corona discharge plasma system for industrial wastewater treatment

In this article, a dual pin-to-plate high-voltage corona discharge system is introduced to study experimentally the gap distance, the contact time, the effect of pin and plate materials, the thickness of ground plate and the conductivity on the amount of Acid Blue 25 dye color removal efficiency from polluted water. A study for the optimum air gap distance between dual pin and surface of Acid Blue 25 dye solution is carried out using 3D-EM simulator to find maximum electric field intensity at the tip of both pins. The outcomes display that the best gap for corona discharge is approximately 5 mm for 15-kV source. This separation is constant during the study of other factors. In addition, an investigation of the essential reactive species responsible for oxidation of the dye organic compounds (O{sub 3} in air discharge, O{sub 3} in water, and H{sub 2}O{sub 2}) during the experimental time is conducted. Three various materials such as: stainless steel, copper and aluminum are used for pins and plate. The maximum color removal efficiencies of Acid Blue 25 dyes are 99.03, 82.04, and 90.78% after treatment time 15 min for stainless steel, copper, and aluminum, respectively. Measurement results for the impact of thickness of an aluminum ground plate on color removal competence show color removal efficiencies of 86.3, 90.78, and 98.06% after treatment time 15 min for thicknesses of 2, 0.5, and 0.1 mm, respectively. The increasing of the solution conductivity leads to the reduction of decolorization efficiency. A kinetic model is used to define the performance of corona discharge system. The models of pseudo-zero-order, pseudo-first-order, and pseudo-second-order reaction kinetics are utilized to investigate the decolorization of Acid Blue 25 dye. The rate of degradation of Acid Blue 25 dye follows the pseudo-first-order kinetics in the dye concentration.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Isotherm Modeling and Techno-Economic Analysis of Contactor Technologies for New Tetraamine-Appended MOF for NGCC Applications

Recently, a family of tetraamine-functionalized metal-organic frameworks (MOF) has been reported as promising sorbent materials for capturing CO2 from flue gas conditions relevant to natural gas combined cycle (NGCC) applications. The main advantages of these materials are their two-step cooperative CO2 adsorption, which gives rise to unusual two step-shaped CO2 adsorption profiles and their high thermal stability. This work presents the modelling of the two-transition isotherm of the tetraamine-appended MOF, N,N'-bis(3-aminopropyl)-1,4-diaminobutane (3-4-3)-appended Mg2(dobpdc), and the techno-economic analysis (TEA) of carbon capture processes utilizing this sorbent. Due to the unusual isotherm shapes of the experimental CO2 adsorption data for tetraamine-appended Mg2(dobpdc) and the strong nonlinearity of CO2 loading with respect to temperature and pressure, we tested two different models which use logistic functions for representing the different isotherm behaviors in the different pressure ranges. The first model uses the quadratic isotherm model in the low-pressure region, the Langmuir isotherm model in the middle pressure range, and the dual site Langmuir isotherm model in the high-pressure range. To model the transition between regions we used the arctangent functions independent of temperature and the thermal effect was accounted by using the Clausius-Clapeyron relation]. The second model is an extension of the weighted dual-site Langmuir isotherm model presented by Hughes et al.. In this extended model, the dual-site Langmuir isotherm is employed in the three transition regions, using temperature-dependent logistic functions to activate or deactivate the isotherm model in the low, middle, and high-pressure ranges. Both models fit the experimental data quite well with root mean squared errors (RMSE’s) of 0.41 and 0.17 for model 1 and model 2, respectively. Since model 2 resulted in a lower RMSE, it was leveraged for the development of the gas/solid contactor models used by the TEA. Specifically, two different contactor models, an axial-flow fixed bed and moving bed contactor, were developed as part of this work. These models are dynamic, pressure-driven, and consist of mass, energy, and momentum conservation equations. A kinetic model was also developed by performing parameter estimation using experimental fixed bed breakthrough data. These models are then used to simulate CO2 capture processes from the flue gas generated from a ~600 gross MW NGCC power plant. A cost model was developed which considers the capital cost of the reactors and the significant operating costs such as steam and electricity. Using NETL’s Framework for the Optimization and Quantification of Uncertainty of Uncertainty and Surrogates tool (FOQUS), which has the capability of linking models built using numerous modelling platforms with derivative-free optimization solvers, a techno-economic optimization of the carbon capture processes was performed which minimizes the cost of capture.

Caballero, Daison↗