Benzinger, W., Daymo, E., Hettel, M., Maier, L., Antinori, C., Pfeifer, P., Deustchmann, O. Reverse Water Gas Shift (RWGS) over Ni – Spatially-Resolved Measurements and Simulations, Chemical Engineering Journal, 2019. https://doi.org/10.1016/j.cej.2019.01.038
Tonkovich, A., and E. Daymo "Process Intensification", Handbook of Thermal Science and Engineering, Ed. F.A. Kulacki. Springer, 2017. https://doi.org/10.1007/978-3-319-32003-8_34-1
Hettel, M., Daymo, E., Deustchmann, O. 3D Modeling of a CPOX-Reformer including Detailed Chemistry and Radiation Effects with DUO, Computers and Chemical Engineering, 2017. https://doi.org/10.1016/j.compchemeng.2017.11.005
Case Study - 1
Scale-up of Algal Ponds
Algal ponds often show reduced performance at increased production scales. Small reactors are easy to mix and produce minimal dead zones. As the volume increases, considerable recirculation zones can form and mixing deteriorates. Design optimization with the use of Computational Fluid Dynamics (CFD) can specify the number and location of turning vanes along with paddlewheel placement.
Strain selection, growth cycles, and nutrient optimization happen in flasks before validation in small ponds. Results must meet economic model targets before moving to scale-up. The surface of a pond can be observed and validated with the use of CFD. The flow beneath the surface cannot be seen and is analyzed with CFD to identify recirculation zones and optimize design details. Success occurs when the final product is growing well, reproducibly and meeting economic targets.
Case Study - 2
Chemical Kinetics and Reactor Design
Designing and debugging problems with catalytic reactors is aided with the use of reactive Computational Fluid Dynamics (CFD). The methane catalytic partial oxidation reaction was studied with the heterogeneous CFD code, DETCHEM DUO. A complete catalytic mechanism on a Rhodium catalyst is available for this reaction. A conjugate heat transfer analysis with a catalytic elementary-step reaction mechanism quantified composition and temperature within a single channel and a multichannel monolith (figures shown below). Heat losses in the laboratory-scale multichannel monolith were quantified to allow predictive scale-up design. Reactor model predictions matched experimental data.
Often industrial catalytic reactions do not have available detailed chemical reaction mechanisms and kinetics. Further, solving reactor problems is often time sensitive and impractical to quickly develop a detailed set of elementary steps.
The reactive CFD code can be used to fit apparent kinetics from available experimental data as a precursor to understanding the root cause of reactor problems and then vet potential solutions. Apparent kinetics may be a few to about ten steps rather than hundreds of steps found in a typical elementary-step mechanism. Further, apparent rate forms may be power law, “1+”, or others as deduced from the experimental data. Data analysis includes the fitting of activation energies, pre-exponentials, reaction orders, and adsorption terms as necessary.
Practical recommendations are made for the industrial user to resolve the reactor problem. Typical recommendations may include changes in design, laboratory set up or operation, or to the catalyst itself.
Some catalytic oxidation reactions are further challenged by competing (and often unwanted) homogeneous reactions. Homogeneous gas phase kinetic mechanisms and parameters are readily available and are used in conjunction with heterogeneous kinetics. Conjugate heat transfer and reaction models assess performance loss from competitive homogeneous reactions, understand safety considerations, and design complete systems that may require reactant safe mixing in a potentially flammable region.
Liquid phase reactions may occur in the bulk or at a catalytic interface. Mixing and mass transfer is often an important consideration for liquid phase reactions due to slow diffusion and low velocity. Reactive CFD analysis can identify design issues and develop alternatives to improve conversion, reduce side products, and lower costs with compact hardware.
Case Study - 3
Droplet Formation can be predicted with knowledge of fluid properties (density, surface tension, contact angle), dimensions (height, width, length), and flowrate of each phase. Design custom microfluidic hardware. Improved visualization (droplet detachment shown on left using Blender and with standard CFD visualization tools on the right) helps to understand results and communicate with stakeholders.
Case Study - 4
Multicomponent Diffusion with Chemical Reaction
Under some conditions (low flows of highly mobile molecules, like hydrogen), accurate diffusion modeling improves reactive system design. OpenFOAM’s standard solvers assume all gas molecules have the same diffusion coefficient (i.e., Schmidt Number = 1). More realistic diffusion models improve accuracy under reactive laminar flow.
In collaboration with the Karlsruhe Institute of Technology (KIT), several diffusion approaches were implemented in DETCHEM DUO to capture concentration profiles in reactive laminar flow. In order of increasing complexity and accuracy, DUO can perform the species balance assuming Fickian mixture-averaged diffusion, the Hirschfelder-Curtiss approximation for mixture-averaged diffusion, and full multicomponent diffusion (solving the Maxwell-Stefan equations).
Concentration profiles were evaluated in a single channel of a Reverse Water Gas Shift (RWGS) monolithic reactor. At low flow (< 0.1 m/s) inlet conditions, which emphasize forward-diffusion of the reactants and back-diffusion of the products for the rich hydrogen feed composition, equivalent concentration profiles are predicted by both DUO and FLUENT.
Correct model assumptions are vital for engineering cost-effective solutions. OpenFOAM coupled with DUO for reacting flow is an excellent CFD platform with vastly lower licensing costs. It is extendable and transparent, to remove guesses about model parameters or methods for reactor design and scale-up.