Oscillations around limit

Please login with a confirmed email address before reporting spam

In my 1D model three variables are solved for. Temperature "T" and a solute concentration "c" are advected through a porous medium. The porous medium reacts with water and produces heat, using conversion variable "X". The reaction front moves as a travelling wave through the domain.

When reacting, the porous medium is converted from X=0 (initial value) to X=1 (fully converted). This means that between the inlet and the wave front, the conversion is 1 and between the wave front and the outlet the conversion is 0. The reaction rate is a function of T, c and X itself. To prevent numerical issues around X=0 and X=1 Heaviside functions are included.

The problem unfortunately is, that throughout the entire domain the values for T and c oscillate heavily, but especially where X=1 (see figure). The mesh is very fine and also the time step is restricted. I have played around with the solver settings but to no avail. Any other possible solutions for this?



0 Replies Last Post Sep 27, 2024, 4:56 a.m. EDT
COMSOL Moderator

Hello Ruben D'Rose

Your Discussion has gone 30 days without a reply. If you still need help with COMSOL and have an on-subscription license, please visit our Support Center for help.

If you do not hold an on-subscription license, you may find an answer in another Discussion or in the Knowledge Base.

Reply

Please read the discussion forum rules before posting.

Please log in to post a reply.

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.