New and Revamped Chemical Engineering Interfaces

Many exciting features for chemical engineering modeling were introduced in COMSOL Multiphysics version 5.0. In this blog post, I will discuss the most important updates. There are some new modeling interfaces, such as theChemistryinterface, and some that have been revamped and improved, like theReaction Engineeringinterface. To begin with, I will recap the modeling interfaces for reaction engineering and mass transport.

A Recap: Modeling Chemical Reactions and Reacting Flow

If you remember, theReaction Engineeringinterface is the tool that is used to experiment with perfectly mixed chemical reactors or reactors where composition and temperature only vary along one direction. These ideal reactor models are usually described with ordinary differential equations (ODEs) that account for changes in the dependent variables (for example, composition and temperature) along one independent variable (usually time or a single space variable).

When model equations of this type are discretized and converted to numerical model equations, the resulting equation systems have relatively few unknowns (degrees of freedom) compared to an equation system that also includes space discretization. The virtue of this kind of model is that the equations can be solved relatively quickly, little RAM is required, and many “what-if experiments” can be done in a short time.

The results from ideal reactor models can also be validated by comparing them with laboratory measurements, which are usually carried out under perfect mixing. Experimental measurements and ideal models can be used for estimating the parameters that describe, for example, reaction kinetics. These parameters can later be used for space-dependent models.

TheReaction Engineeringinterface sets up small ODE systems of space-independent problems.

In some tubular reactors — for example, plug flow reactors — the composition and temperature can often be decently well described by accounting for variations along the axial direction and neglecting variations in the radial direction. Such models are usually referred to as ideal tubular reactor models. At steady state (no variations in time), these models are described by ODEs with the axial space variable (z) as the only independent variable.

However, we usually need to validate the assumptions done in ideal models by comparing them with space-dependent models and models that describe variations along several independent variables (x,y,z, andt). In addition, scale-up from lab scale to pilot or full scale often introduces spatial and temporal variations in composition, pressure, temperature, and other modeled fields. In order to facilitate this work, COMSOL Multiphysics includes a feature called theGenerate Space-Dependent Modelthat exports the chemical system described in the space-independent model to a space-dependent model.

A perfectly mixed model and a space-dependent model.

With the NewChemistryInterface, Chemical Reactions Are Made Easy

In COMSOL Multiphysics version 5.0, theChemistryinterface shares some similarities with the Materials node in that it harbors information that can be used by many model interfaces in several model components. Whereas the Materials node contains physical properties, theChemistryinterface can contain a library of thermodynamic and kinetic properties as well as reaction mechanisms of a given chemical system.

The functionality is similar to what is available in theReaction Engineeringinterface. However, the primary purpose is not to solve the model equations, but rather to administrate on a global level the chemical reactions and their kinetic and thermodynamic properties for use in space-dependent models. As such, it does not solve any equations in itself.

TheChemistryinterface defines reaction mechanisms and species properties.

The figure below shows an example of model results in which you start with theChemistryinterface, and then continue using theTransport of Diluted SpeciesandDarcy’s lawinterfaces to model the reacting flow and the velocity and pressure fields, respectively, in a 3D packed bed reactor.

A 3D packed bed reactor application with a holed bottom plate. The lines represent the flow streamlines of reacting gas, and the color is the concentration of reactant in moles/m³. Red represents high concentration and blue indicates low concentration.

When you have added theChemistryinterface and defined the mapping to the species in the transport interfaces, the reaction rates appear automatically in theReactionfeature.

With a Chemistry node in the application, the reaction rates are automatically defined on theReactionfeature in the transport interfaces.

A More Efficient Way to Export Reacting Systems from 0D to Space-Dependent Components

When you use theReaction Engineeringinterface, the space-dependent model created by theGenerate Space-Dependent Modelfeature is not linked back to theReaction Engineeringinterface anymore. When the generation is done, aChemistryinterface is instead created separately from theReaction Engineeringinterface. This way, the study settings are simplified. TheChemistryinterface will now include the reaction properties originally provided by theReaction Engineeringinterface.

The introduction of theChemistryinterface makes the modeling steps easier.

The Ability to Handle a Large Number of Reactions

TheSpecies GroupsandReaction Groupsfeatures have been introduced to make it easier to handle large amounts of reactions and species. These features list all reaction properties in table format. Reaction kinetics and thermodynamics included in CHEMKIN files can be imported automatically into these groups.

TheReaction Groupsfeature makes it easier to maintain large numbers of chemical reactions. The example shown above is from theHomogeneous Charge Compression Ignition of Methane model.

Smart Surface Reactions

If a chemical species is designated to be a surface species, the chemical reaction in which it partakes is automatically set up as a surface reaction. Generation of space-dependent models adds these to surfaces in the settings for the boundary conditions.

TheHydrocarbon Dehalogenation in a Tortuous Microreactor modelfrom the Model Gallery demonstrates the new surface reactions capabilities. The surface species can be marked “surface” or named with an “(s)” postfix to create a surface reaction.

New Reactor Model: Generic CSTR

TheReaction Engineeringinterface includes a new reactor type calledGeneric CSTR. This continuous stirred tank reactor has no restrictions of constant mass or constant volume. That means that it can be filled or drained during the simulation. Multiple inlets are possible in this reactor.

A stirred tank reactor with two inlets (top) is modeled with the new Generic CSTR, as shown inthis example from the Model Gallery. Top image by Robert Ashe. Licensed underCC BY 2.5viaWikimedia Commons.

Improved Thermodynamic and Transport Properties for Mixtures

In theReaction EngineeringandChemistryinterfaces, there are improved model properties of mixtures. Viscosity, heat capacity, thermal conductivity, and density are examples of properties that can be calculated for both gas mixtures and liquid mixtures.

An example of automatic mixture properties.

Model Your Kinetics and Equilibrium with Activities

There is now an option to calculate rate expressions and equilibrium relations with the species activities as input rather than concentrations.

Activities can be used in the expressions for equilibrium and kinetics.

In the case of non-ideal mixtures, the concentration can be replaced with the activity. The dimensionless activitya_idepends on the species concentrationc_i, the activity coefficient\gamma_i, and the standard state concentrationc_{0s}, the latter which is oftentimes set to 1 mol/m³but can be tuned by the user:

I hope you will find these new features useful in your simulation work. Happy modeling!