|Introduction||Initial and Boundary Conditions|
|Mass Transport Model||Solve the PDE Model|
|Domain||Post-processing and Visualization|
A catalytic converter is an emission control device that reduces certain pollutants in exhaust gas from an internal combustion engine. With the assistance of catalysts, molecules such as nitric oxide, , and carbon monoxide, , can be oxidized and converted to substances that are somewhat less harmful to the environment.
This study is to simulate mass transportation of molecules within a catalytic converter and its evolution in time. The catalytic converter modeled in this example consists of two parallel plates (active surfaces) with a catalyst on them. A gas containing molecules enters the domain from the bottom, flows across the reacting plates, and then exits the converter through the top. During this process the molecules near the plates are continuously oxidized into carbon dioxide:
In the Post-Processing section we will measure the converter performance by calculating the net reduction in the concentration: between the flow inlet and outlet :
For this the mean concentration at the boundary is computed with a boundary integration using NIntegrate:
The symbols and corresponding units used here are summarized in the Nomenclature section.
Please refer to the information provided in "Mass Transport Tutorial" for a more general theoretical background for heat transfer analysis.
The conservative transport equation (1), which is derived from mass conservation, is used to solve for the concentration distribution in a mass transport model:
is the concentration of the transported species ,
is the species diffusivity ,
is the velocity field of possible flow inside of the medium ,
is the mass production/consumption, the volumetric reacting rate of the species .
In this model the reaction of the carbon monoxide, , takes place on the active surfaces only, and will be modeled with a mass flux boundary condition. Since no other reaction occurs within the converter, the volumetric reacting rate is set to zero and the transport equation simplifies to:
The catalytic converter model consists of two parallel reacting plates. If the depth of the plates is reasonably longer than its length and width, then the variation of concentration in the direction can be neglected. Therefore a two dimensional model is sufficient to represent the 3D reactor.
Due to the symmetry along the axis it is effective to only make use of the right half of the reactor as the simulation domain . Here and denote the flow inlet and flow outlet. The active surface and the symmetric boundary are symbolized as and , respectively.
At the top boundary an outflow boundary condition is used to model the flow outlet where molecules are transported out of the domain.
On the active surface the carbon monoxide is absorbed and transformed by the catalyst , and can be modeled by a mass flux boundary condition.
Assuming that each carbon monoxide molecule occupies exactly one catalyst site during the reaction, then the concentration of remaining sites, , is equal to the difference between the total concentration of active sites, , and the number of sites occupied by the absorbed carbon monoxide:
A default symmetric boundary condition is implicitly applied on the symmetric boundary .
See this note about improving the visual quality of the animation.
With the presence of the catalyst , molecules near the right boundary (active surface) are continuously absorbed and converted into carbon dioxide, resulting in a layer-like concentration field. The thickness of this concentration layer grows with time and reaches its steady state around .
Since the concentration has little variation outside this layer, it is known as the "reaction dead space". To minimize the reaction dead space and improve the overall converting effectiveness, the catalyst converter are often made in a shape of honeycombs rather than parallel plates.
To determine the exact time when the system reaches the steady state, we can use WhenEvent to check whether the concentration field has converged within the pre-defined threshold.
To measure the effectivity of the catalytic converter, the net reduction is calculated between the flow inlet and outlet . For this the mean concentration at the boundary is computed with a boundary integration: