MultiComponentReactiveMixtureProject

Project package for the implementation of a model for heat and multi-component, reactive gas phase transport in porous media.

The model is applied to the simulation of a Photo-Thermal catalytic reactor. Concentrated solar irradiation is used to supply heat input into the reactor to drive endothermal chemical reactions. The reactor utilises a catalyst coated porous material which is directly irradiated and through which the reactive gas mixture passses. A 3D (square prismatic) domain of porous material is modelled. Model equations and correlations for effective properties of heat transfer through the porous material are taken from VDI Heat Atlas (chapter D6.3). The simulation utilizes the finite volume method as implemented in the VoronoiFVM.jl package.

Check out the package with git:

git clone https://github.com/DavidBrust/MultiComponentReactiveMixtureProject

After checking out or updating, run from the project root directory:

$ julia --project -e "using Pkg; Pkg.instantiate()"

Computational notebooks based on the Pluto.jl package are used to illustrate this package's functionality, located in the notebooks subfolder.

Notebooks

• LoschmidtTernary.jl: simulation of diffusion in ternary gas mixture (verification of isothermal model)
• Thermodiffusion.jl: demonstration of thermodiffusion (reproduces figures shown in the article)
• Uphill_Diff1D.jl: illustration of transport against a concentration gradient (convection driven) in iso-thermal conditions
• Residual_MoleFrac1D.jl: investigation of grid refinement on residual molar fractions
• PTReactorDemo.jl: application to photo-thermal reactor in 2D and 3D geometry and non-isothermal conditions

Examples shown in the article

Running the following files allows to reproduce the figures shown in the article:

• scripts/PTR3D.jl: script for 3D Photo-thermal reactor simulation with data export for plot (Figure 5, Section 4.1)
• Thermodiffusion.jl: Pluto notebook demonstrating thermodiffusion, export plot (Figure 8, Section 4.2)

To run the script PTR3D.jl, run from the Julia REPL in the project root:

include("scripts/PTR3D.jl")
PTR3D.run()

The script will run the simulation and will export the solution in VTK format to a new folder for plotting with ParaView. Use the provided ParaView state file data/PTR3D.pvsm to reproduce Figure 5 in Section 4.1.