Next Article in Journal
Improvements and Spatial Dependencies in Energy Transition Measures
Next Article in Special Issue
Study on the Convective Heat Transfer and Fluid Flow of Mini-Channel with High Aspect Ratio of Neutron Production Target
Previous Article in Journal
Towards Scalable Economic Photovoltaic Potential Analysis Using Aerial Images and Deep Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Modeling of Sublimation of Ammonium Carbamate Applied to Supply System of NOx Reductant

1
Environment System Research Division, Korea Institute of Machinery and Materials, 156 Gajeong-ro, Yuseong-gu, Daejeon 34103, Korea
2
AI and Mechanical System Center, Institute for Advanced Engineering, 175-28, Goan-ro 51beon-gil, Baegam-myeon, Cheoin-gu, Yongin-si, Yongin 17180, Korea
3
Department of Environment & Energy Mechanical Engineering, University of Science and Technology, 217 Gajeong-ro, Yuseong-gu, Daejeon 34113, Korea
*
Author to whom correspondence should be addressed.
Energies 2021, 14(13), 3795; https://doi.org/10.3390/en14133795
Submission received: 17 May 2021 / Revised: 16 June 2021 / Accepted: 21 June 2021 / Published: 24 June 2021

Abstract

:
Recently, ammonium carbamate (AC) has attracted attention as a substitute for urea, which is a commonly used reductant for NOx emitted from combustion engines. The AC exists as a solid at room temperature, and it is decomposed to NH3 and CO2 gases by heating. Therefore, adequate heat transfer is an essential issue in the design of AC pyrolysis reactor. In this study, a numerical model that describes the sublimation of AC was developed. For modeling, this study considered the three different calculation zones: solid-phase zone, gas-phase zone, and sublimation zone. Additionally, during the sublimation process, collapse of upper solid AC into the hollow space below by the effect of gravity is considered. As a result, it is presented that the modeling shows reasonable information about the AC sublimation in a reactor, such as temperatures in a reactor, pressure of reactor, and flow rate of sublimated gas. However, it is also found that accurate prediction of spatial temperature distribution is challenging because it is related to the accurate prediction of the internal shape of AC and its collapse in a reactor.

1. Introduction

As one of the major sources of air pollutants, NOx emissions have been strictly regulated. [1]. However, approximately 500 million tons of NOx are still emitted annually from vehicles worldwide [2], affecting the quality of human life. Selective catalytic reduction (SCR) has been the most commonly used technology to reduce NOx emissions from vehicles [3,4]. The SCR refers to a catalyst system that converts NO or NO2 into N2 and H2O by reacting with NH3, and it shows an exceptionally high NOx reduction efficiency and high system reliability. However, NH3 is needed as a NOx reducing agent in the SCR. It is a toxic substance and very dangerous in the event of a leak. Additionally, it is difficult to store in a liquid or gaseous state. Therefore, many vehicles currently use urea (CO(NH2)2) water solution to supply NH3. The urea water solution is decomposed into NH3 and CO2 when it is sprayed onto high-temperature exhaust gases. However, effective decomposition of urea solution usually requires a high exhaust gas temperature of 200 °C or more [5]. Due to this requirement, a substantial amount of NOx is emitted during low temperature operating conditions [6]. One of the other disadvantages is that the urea solution can freeze below −12 °C [7]. Additionally, it is not easy to keep the uniform distribution of ammonia in a catalyst during transient engine operating conditions [8].
The AC (ammonium carbamate, H2NCOONH4) has been proposed as an alternative of urea [5]. The AC can be sublimated from solid to gas mixture of NH3 and CO2 at temperatures (approximately 60 °C at 1 atm) lower than those for urea. Another advantage of AC is that the storage capacity of NH3 is about three times larger than that of urea solution. Besides, the NH3 gas produced from AC has better mixing characteristics with exhaust gas than that of urea solution. In general, AC is stored in a reactor. The reactor’s pressure increases when it is heated because some amount of AC is sublimated in the reactor. Then, the NH3 gas is injected into the exhaust gas using the elevated pressure [5]. However, it is difficult to control the sublimation rate of AC because of spatially uneven distribution of AC mass and temperature in the reactor. For example, the AC contacting with the heat sources is sublimated at first, and a gas space is created there. After that, the AC that exists spatially apart from the heat sources is sublimated by the effect of convective heat transfer. Such a temporally and spatially complex sublimation process results in a non-uniform sublimation rate. Therefore, optimal design of the reactor is necessary to control the sublimation rate effectively. Previous researchers have investigated the simple design parameters of AC reactor only by experimental methods [9,10,11]. However, it was not easy to clearly understand the heat and mass transfer phenomena in the reactor. Such difficulties hinder the optimal design of the reactor.
Many numerical models have been developed to solve engineering problems of phase change from solid to liquid [12,13,14] or liquid to gas [15,16]. However, numerical models for phase change from solid to gas have not been proposed well, because it is not a common phenomenon encountered in engineering problems. Kim et al. [17] developed a numerical method to investigate the frosting of water. They tried to predict the variation of heat transfer coefficient owing to the frost. They also simulated the growth of frost by the deposition of H2O molecules in the air on a cold surface. Zhang et al. [18] calculated the amount of heat energy removed by sublimation, evaporation, and freezing that occurred on a porous plate, which was soaked in the water. However, the porous plate’s shape did not vary, and only the heat and mass transfer were analyzed. Zadravec et al. [19] investigated the variation of the interface height by simulating the material’s freeze-drying process in a vial. However, this study considered only the height of the interface. It was not developed to interpret the three-dimensional sublimation phenomenon. As seen above, there have been many simulation studies on the phase change from gases to solid or solid to gases. However, as far as we know, there are no three-dimensional simulation studies about AC sublimation. This seems to be because AC is the newly introduced material for ammonia storage.
This study’s objective is to develop a numerical model, which simulates the sublimation of AC for use in the optimal design of a reactor.

2. Materials and Methods

2.1. Experiment Set-Up

Experiments on AC sublimation were conducted to verify the numerical model developed in this study. As shown in Figure 1, the experimental apparatus is divided into three parts: AC reactor (pyrolysis space), heated water supply device, and heated air supply device. The water heated at 90 °C is circulated outside the reactor to supply heat. Two glass plates are installed on the top ceiling of the reactor to visualize the shape of solid-phase AC. The hot air heated to 120 °C is supplied between the two glass plates to prevent the re-solidification of NH3 and CO2 gases on the surface of the glass plate. In the AC reactor, four fins are installed to enhance heat transfer characteristics. Five thermocouple rods (TRF1, TRF2, TRF3, TRA1, TRA2) are installed in the reactor to measure temperature distribution. TRF1, TRF2, and TRF3 are installed near a fin that are expected to sublimate first. TRA1 and TRA2 are installed in the middle between two fins that are expected to sublimate later. In each of the rods, four thermocouples are implanted at regular intervals. The entire experimental apparatus is covered with insulation material.
The AC exists in the form of powder. However, when it is exposed at a certain condition for a long time, the AC may harden and agglomerate. In this study, the AC was crushed into small powders before the tests.
NH3 and CO2 gas mixtures were exhausted from the bottom of the vessel. The pressure was measured using a pressure sensor (PSH, Sensys, Korea) at the outlet. In the experiment, after 30 min of warm water circulation, the outlet valve was opened, and the flow rate of gas was measured by MFC (mass flow controller). The outlet gas line was heated to 120 °C to prevent clogging, owing to the re-solidification of NH3 and CO2 gases.

2.2. Numerical Method

2.2.1. Sublimation Numerical Analysis Model

Three kinds of calculation zones are considered in the numerical model developed in this study: solid-phase zone, gas-phase zone, and sublimation zone. Only solid AC is present in the solid-phase zone. The gas-phase zone also only contains gases. The sublimation zone is a region where the temperature reaches the critical temperature for sublimation. Both solid and gas exist in the sublimation zone.
The temperature of every zone is monitored at each iteration of the calculation to classify the zone. When the zone’s temperature is higher than the sublimation critical temperature (Tth), the solid zone transforms into a sublimation zone. Here, the Tth is determined by the Figure 2, which is obtained from an experiment of reference [9]. The solid-phase and gas-phase mass ratio in the zones (ρs) is calculated at each iteration of calculation to classify the zone. When one mole of solid AC is pyrolyzed, one mole of CO2 and two moles of NH3 gas are generated. Generated CO2 and NH3 gas are considered by source term of fluid in the sublimation zone. If ρs is decreased to zero, the sublimation zone transforms to a gas-phase zone. In the gas-phase zone, conservation equations of computational fluid dynamics (CFD) are applied to calculate the conserved variables.
Figure 3 shows the concept of a model that explains how to simulate the change in shape of a solid AC where sublimation occurred. Figure 3 is expressed two-dimensionally for ease of understanding, but the real model was calculated three-dimensionally. Step 1 is the initial step at Time i. The distinction between the solid-phase and gas-phase zones is evident, and there is no sublimation zone. Step 2 is the heating stage. The solid-phase zone is transformed to the sublimation zone when the cell temperature (Ts) is greater than the sublimation critical temperature (Tth). The mass ratio of solid-phase and gas-phase (ρs) is also calculated in Step 2. In Step 3, the sublimation zone’s cell transforming to the gas-phase zone is determined based on the ρs. In Step 4, the solid phase zone is shifted in the direction of gravity to fill the voids created by the phase transformation. It was assumed that solid AC was in powder form, as it had weak bonding forces, and it shifted in the direction of gravity without plastic action. Steps 1–4 are repeated until the residue converges below a specific value. Furthermore, the time step is shifted forward from i to i + 1. For the next time step, the zone information (solid, gas, sublimation) calculated in the previous iteration is used as the initial value. In this study, the amount of sublimated AC, the pressure of the pyrolysis vessel, and the amount of NH3 gas generated are predicted by continuing this calculation procedure.
A commercial CFD program ANSYS FLUENT 19.2, was used in this study. The general mass, momentum, and energy conservation equations of the commercial program were used. As a turbulent model, the standard kε model was used. The latent heat of sublimation of AC, Se, was added to the energy equation as follows:
ρ c p T t + ·   ρ c p T = k T + S e ,
where ρ is the density, cp is the specific heat, T is the temperature, k is the thermal conductivity. The generation of NH3 and CO2 gases by AC sublimation is incorporated in the species transport equation as follows:
c k t + ·   c k u = D k c k + S k
where ck is the concentration of each species, u is the velocity, Dk is the diffusion coefficient of each species, and Sk is the source term for each species generation. The NH3 and CO2 gases generated by AC sublimation are calculated in terms of Sk.
Properties of gases and solid material of reactor are listed in Table 1 which are mainly citing NIST (National Institute of Standard and Technology) data. The sublimation heat of solid AC is 159.94 kJ/mol [20].
Commercially available ANSYS Fluent 19.2 software (Canonsburg, Pennsylvania US) was used. Sublimation of solid AC was calculated by a user-defined function. The calculating elements numbered from 1.1 million to 1.7 million were tested for grid independence, and the number of elements was decided to about 1.3 million. The coupled governing equations were discretized using the finite volume method and solved with SIMPLE algorithm. For actual analysis, convergence conditions of solution were taken to be less than 10−5 of relative uncertainty level.

2.2.2. Model Validation

The sublimation model developed in this study was verified by comparing it with experimental results. As shown in Figure 4, the calculating area is identical to the pyrolysis space inside the AC reactor. The inlet boundary conditions of air supplied to the space between the two glass plates are as follows: the temperature is 120 °C, and the velocity of flow is 1 m/s, respectively. The inlet boundary conditions of warm water for supplying heat are a temperature of 90 °C and a flow rate of 5 L/min. All the outlet boundary conditions, such as the heating air, vessel’s outlet, and the warm water outlet, are set as the pressure outlet. All the boundary conditions of the outside wall (except the top of the glass) are set as adiabatic conditions. The glass plates at the top of the experimental apparatus are exposed to ambient air. Therefore, a general convection heat transfer condition is applied to the glass plate surface. The glass plate surface’s heat transfer coefficient is 6 W/m2 °C, and the ambient air temperature is set to 25 °C. All the wall boundary layers exposed to fluid are in the no-slip condition. The pyrolysis vessel’s initial temperature is 37 °C, which was the air temperature of the experimental pyrolysis vessel when that of the warm water attained 90 °C. The model is calculated transiently. Air, NH3, and CO2 are assumed to be ideal gases. The time step of the calculation is 1 min, and the calculation results are summarized every 10 min.

3. Results and Discussion

The variations of the shape of AC in the experiment and simulation are compared in Figure 5. In the experiment, the solid AC area decreases gradually owing to sublimation after heat is supplied. In particular, the solid regions adjacent to the fins and side wall (where heat is supplied directly) are reduced substantially. It means that the fin’s large contact area with AC enhances the sublimation. Hence, it can be mentioned that the increasing the number of fins is an effective way to ensure a consistent pyrolysis performance when designing an AC reactor. The simulation also show similar results of the rapid sublimation of solid AC around the sidewall and the fins.
The reactor pressures and the gas flow rates of experiment and simulation are compared in Figure 6. In the experiment, the reactor pressure increases to 0.1 bar at 30 min because the exhaust line is closed. After the exhaust valve is opened, it decreases to 0.03–0.04 bar. The flow rate of the exhaust from the reactor decreased from 0.27 l/min at 30 min to 0.05 l/min at 150 min in the experiment. The simulation result of pressures and flow rate show similar changes with experiments. From the Figure 6, it can be said that the developed sublimation model well predicted reactor pressure and the ammonia gas flow rate.
The averaged pyrolysis rate was calculated using the above simulation results with 30-min intervals, and is shown in Figure 7. The averaged pyrolysis rate during the first 30 min is 4.32 g/min, and it decreases rapidly to 0.2 g/min. This is because the AC adjacent to the sidewall or fins are sublimated rapidly by conductive heat transfer in the early stage of the experiment. After that, the sublimation rate of the remaining AC is controlled by the heat supplied by convection.
The simulated temperatures are compared with the experimental results in Figure 8. All the temperatures of thermocouples increase linearly by heating initially, and the temperatures are saturated after sublimation has occurred. This means that all energies supplied to the reactor are used as the latent heat of sublimation. The simulation also shows the saturation of temperatures. However, in the experiment, the temperatures are saturated quickly, whereas it takes longer for saturation in the simulation. It seems that this difference comes from the uncertainties of simulation, such as the inexact values of heat capacity or density of solid AC the heat transfer characteristics of porous space in the solid AC may also be affected.
Comparing the temperature distribution of the experiment reveals that TRF3 and TRA2 are the hottest (70–80 °C). It seems that the thermocouples in these regions were rapidly exposed to hot gases because of fast sublimation by a large amount of heat supplied from the sidewall. The temperatures of TRF1 and TRF2 are the second highest (65–72 °C) because these were close to the fin and exposed to the gases after about 90 min. The temperature of TRA1 are the lowest (62–70 °C) because it was farthest from the fin and sidewall. It seems that the thermocouple rod of TRA1 was buried in AC for the entire experiment. Similarly, in the simulation, the temperatures of TRF3 and TRA2 are the highest, and that of TRA1 is the lowest. These are similar with the experiment results.
Generally, the temperature measured in the experiment is high in the order of height. The temperatures are low in the AC or near the AC because the AC absorbs heat during sublimation. However, there is an exceptional temperature distribution in the experiment. For example, the temperature at Point 4 of TRA1 is higher than that at Points 2 and 3 of TRA1 during the entire experiment. Additionally, the temperatures at Point 4 of each TRF1, TRF2, and TRA2 are higher than that at Points 2 and 3 of each TRF1, TRF2, and TRA2 during the first 100 min of the experiment. It seems that the voids in the AC cause an unexpected temperature distribution. The thermocouple exposed to the gases in a void shows a higher temperature than that of the thermocouple buried in AC because of little or no latent heat of AC sublimation. After 100 min, all temperatures of Point 4 were lower than point 2 and point 3. This means that the AC near Point 4 collapsed, and the thermocouples at Point 4 were buried in AC again.
In the simulation, the temperature is high in the order of height without exception. The reason is that the continuous collapse of upper solid AC into the lower space by the effect of gravity was taken into account in the simulation. From the above results, it can be mentioned that it is challenging to accurately predict the spatial temperature distribution in the reactor where AC sublimation occurs because it is necessary to accurately predict the collapse of the solid AC existing in the upper or surrounding area into the lower empty space.

4. Conclusions

In this study, a numerical model for the sublimation of AC is newly developed. To effectively simulate the sublimation of AC in a reactor, the calculation zones are divided into three zones with different characteristics: solid-phase zone, gas-phase zone, and sublimation zone. Depending on the sublimation critical temperature (Tth) and the solid-phase and gas-phase mass ratio (ρs) in the zone, the solid-phase zone changes to the sublimation zone and then to the gas-phase zone. This study also simply modeled the collapse of upper solid AC into the lower empty space, which is created by sublimation.
The simulation results were compared with the experimental results of AC sublimation in a simplified reactor. It is found that the modeling can show very credible information about the phenomena of AC sublimation, such as temperature change with time, the reactor pressure, and sublimated gas flow rate. It is also conjectured that the heat transfer in the reactor is reasonably simulated because the shape variation of solid AC in the experiment is similar to that in the simulation. However, the accurate prediction of spatial temperature distribution is challenging because it is related to the accurate prediction concerning the collapse of the solid AC existing in the upper or surrounding area into the lower empty space, caused by sublimation.

Author Contributions

Model development and simulation, J.-H.N., experiment and writing-original draft, S.-H.W., writing-review and editing, H.R., conceptualization, and supervision, H.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Korea Institute of Machinery & Materials, Republic of Korea, under the project name of “The ultra clean machinery technologies of future electric generation or power systems against PM2.5 fine particles and greenhouse gases”, grant number NK231A.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xu, L.; Tsona, N.T.; You, B.; Zhang, Y.; Wang, S.; Yang, Z.; Xue, L.; Du, L. NOx enhances secondary organic aerosol formation from nighttime γ-terpinene ozonolysis. Atmos. Environ. 2020, 225, 117375. [Google Scholar] [CrossRef]
  2. 2014 National Emissions Inventory Technical Support Document. Available online: https://www.epa.gov/sites/production/files/2018-07/documents/nei2014v2_tsd_05jul2018.pdf (accessed on 23 June 2021).
  3. Gou, X.; Wu, C.; Zhang, K.; Xu, G.; Si, M.; Wang, Y.; Wang, E.; Liu, L.; Wu, J. Low temperature performance of selective catalytic reduction of NO with NH3 under a concentrated Co2 atmosphere. Energies 2015, 8, 12331–12341. [Google Scholar] [CrossRef] [Green Version]
  4. Li, L.; Lin, W.; Zhang, Y. A new dynamic injection system of urea-water solution for a vehicular select catalyst reduction system. Energies 2017, 10, 12. [Google Scholar] [CrossRef] [Green Version]
  5. Fulks, G.; Fisher, G.B.; Rahmoeller, K.; Wu, M.C.; Herde, E. A review of solid materials as alternative ammonia sources for lean NOx reduction with SCR. SAE Int. 2009. [Google Scholar] [CrossRef]
  6. Gramigni, F.; Nasello, N.D.; Selleri, T.; Nova, I.; Tronconi, E.; Dieterich, S.; Weibel, M. Unexpected low-temperature deNOx activity of AdSCR systems for cold start NOx abatement. Emiss. Control Sci. Technol. 2020, 6, 402–409. [Google Scholar] [CrossRef]
  7. Choi, B.C.; Kim, Y.K.; Jhung, W.N.; Lee, C.H.; Hwang, C.Y. Experimental investigation on melting characteristics of frozen urea-water-solutions for a diesel SCR de-Nox-system. Appl. Therm. Eng. 2013, 50, 1235–1245. [Google Scholar] [CrossRef]
  8. Jeong, S.; Kim, H.; Kim, H.; Kwon, O.; Park, E.; Kang, J. Optimization of the urea injection angle and direction: Maximizing the uniformity index of selective catalytic reduction system. Energies 2020, 14, 157. [Google Scholar] [CrossRef]
  9. Kim, Y.; Raza, H.; Lee, S.; Kim, H. Study on the thermal decomposition rate of ammonium carbamate for a diesel NOx reducing agent-generating system. Fuel 2020, 267, 117306. [Google Scholar] [CrossRef]
  10. Kim, H.; Yoon, C.; Lee, J.; Lee, H. A study on the solid ammonium SCR system for control of diesel NOx emissions. SAE Int. 2014. [Google Scholar] [CrossRef]
  11. Lacin, F.; Kotrba, A.; Hayworth, G.; Sullivan, H.; Tatur, M.; Jacques, J.; Tomazic, D.; Cho, H. Solid SCR: Demonstrating an improved approach to NOx reduction via a solid reductant. SAE Int. 2011. [Google Scholar] [CrossRef]
  12. Costa, M.; Oliva, A.; Perez Segarra, C.D.; Alba, R. Numerical simulation of solid-liquid phase change phenomena. Comput. Methods Appl. Mech. Eng. 1991, 91, 1123–1134. [Google Scholar] [CrossRef]
  13. Huo, Y.; Rao, Z. Lattice Boltzmann simulation for solid-liquid phase change phenomenon of phase change material under constant heat flux. Int. J. Heat Mass Transf. 2015, 86, 197–206. [Google Scholar] [CrossRef]
  14. Chabot, C.; Gosselin, L. Solid-liquid phase change around a tube with periodic heating and cooling: Scale analysis, numerical simulations and correlations. Int. J. Therm. Sci. 2017, 112, 345–357. [Google Scholar] [CrossRef]
  15. Gong, S.; Cheng, P. A lattice Boltzmann method for simulation of liquid-vapor phase-change heat transfer. Int. J. Heat Mass Transf. 2012, 55, 4923–4927. [Google Scholar] [CrossRef]
  16. Fang, W.Z.; Tang, Y.Q.; Yang, C.; Tao, W.Q. Numerical simulations of the liquid-vapor phase change dynamic processes in a flat micro heat pipe. Int. J. Heat Mass Transf. 2020, 147, 119022. [Google Scholar] [CrossRef]
  17. Kim, C.; Lee, J.; Lee, K.S. Numerical modeling of frost growth and densification on a cold plate using frost formation resistance. Int. J. Heat Mass Transf. 2017, 115, 1055–1063. [Google Scholar] [CrossRef]
  18. Zhang, C.; Bu, X.; He, J.; Liu, C.; Lin, G.; Miao, J. Simulation of evaporation and sublimation process in porous plate water sublimator based on a reduced CFD model. Int. J. Heat Mass Transf. 2020, 154, 119787. [Google Scholar] [CrossRef]
  19. Zadravec, M.; Casar, Z.; Hribersek, J.R.M. CFD based determination of sublimation mass flux for lyophilization inside a vial. Int. J. Comput. Methods Exp. Meas. 2020, 8, 47–60. [Google Scholar] [CrossRef]
  20. Ahn, C.; Lee, H.; Lee, M.; Chang, Y.; Han, K.; Rhee, C.; Kim, J.; Chun, H.; Park, J. Determination of Ammonium Salt/Ion Speciation in the CO2 Absorption Process Using Ammonia Solution: Modeling and Experimental Approaches. Energy Procedia. 2011, 4, 541–547. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic of experimental set-up for validation of AC sublimation model.
Figure 1. Schematic of experimental set-up for validation of AC sublimation model.
Energies 14 03795 g001
Figure 2. Equilibrium pressure of AC with respect to the temperature measured in a closed chamber [9]. Reprinted with permission from ref. [9]. Copyright 2021 Elsevier.
Figure 2. Equilibrium pressure of AC with respect to the temperature measured in a closed chamber [9]. Reprinted with permission from ref. [9]. Copyright 2021 Elsevier.
Energies 14 03795 g002
Figure 3. Concept of the model: Changes in the shape and location of solid AC where sublimation occurs due to heat transfer.
Figure 3. Concept of the model: Changes in the shape and location of solid AC where sublimation occurs due to heat transfer.
Energies 14 03795 g003
Figure 4. Schematic diagram and boundary definitions of AC sublimation model validation.
Figure 4. Schematic diagram and boundary definitions of AC sublimation model validation.
Energies 14 03795 g004
Figure 5. Variation of the shape of solid ammonium carbamate when the reactor is heated.
Figure 5. Variation of the shape of solid ammonium carbamate when the reactor is heated.
Energies 14 03795 g005
Figure 6. Pyrolysis reactor (a) pressure and (b) exhaust gas flow rate variation.
Figure 6. Pyrolysis reactor (a) pressure and (b) exhaust gas flow rate variation.
Energies 14 03795 g006
Figure 7. Ammonium carbamate pyrolysis rate variation.
Figure 7. Ammonium carbamate pyrolysis rate variation.
Energies 14 03795 g007
Figure 8. Variation in temperature distribution of (a) TRF1, (b) TRF2, (c) TRF3, (d) TRA1, and (e) TRA2.
Figure 8. Variation in temperature distribution of (a) TRF1, (b) TRF2, (c) TRF3, (d) TRA1, and (e) TRA2.
Energies 14 03795 g008
Table 1. Physical properties of working fluid and solid.
Table 1. Physical properties of working fluid and solid.
DescriptionDensity [kg/m3]Viscosity [Pa-s]Specific Heat [J/kg-K]Thermal Conductivity [W/m-K]
FluidAir0.8982.264 × 10−510110.032
Ammonia0.6891.015 × 10−521580.025
Carbon Dioxide0.7301.801 × 10−59080.022
Water965.30.315 × 10−342060.675
SolidSUS 3047930-50216.2
Insulator (Urethane)70-10450.026
Ammonia Carbamate1380-16702.0
Glass2800-7500.7
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Woo, S.-H.; Noh, J.-H.; Raza, H.; Kim, H. Numerical Modeling of Sublimation of Ammonium Carbamate Applied to Supply System of NOx Reductant. Energies 2021, 14, 3795. https://doi.org/10.3390/en14133795

AMA Style

Woo S-H, Noh J-H, Raza H, Kim H. Numerical Modeling of Sublimation of Ammonium Carbamate Applied to Supply System of NOx Reductant. Energies. 2021; 14(13):3795. https://doi.org/10.3390/en14133795

Chicago/Turabian Style

Woo, Sang-Hee, Jung-Hun Noh, Hassan Raza, and Hongsuk Kim. 2021. "Numerical Modeling of Sublimation of Ammonium Carbamate Applied to Supply System of NOx Reductant" Energies 14, no. 13: 3795. https://doi.org/10.3390/en14133795

APA Style

Woo, S. -H., Noh, J. -H., Raza, H., & Kim, H. (2021). Numerical Modeling of Sublimation of Ammonium Carbamate Applied to Supply System of NOx Reductant. Energies, 14(13), 3795. https://doi.org/10.3390/en14133795

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop