CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Passive Scalar Diverging in micromixer simulation

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   February 21, 2022, 11:13
Default Passive Scalar Diverging in micromixer simulation
New Member
Alireza Miraliakbar
Join Date: May 2021
Posts: 2
Rep Power: 0
ar.miraliakbar is on a distinguished road
Hello everyone,

I would like to ask for help for my thesis validation simulation where I am reproducing the results of Tsai and Wu. The reference is given below in the link:

An efficient micromixer based on multidirectional vortices due to baffles and channel curvature

I have added the passive scalar as a representation of molar fraction to the simpleFoam solver and ran 2 simulations for Re = 1 and Re = 81 for comparing the results with the experimental ones. Every thing is Ok for the Re=1 but for the Re = 81 the passive scalar is diverging despite the convergence of U and p. Based on my calculation and the data given in the article the Re = 81 is equal to Normal Velocity Inflow velocity of 0.606 m/s.
I did the meshing with Ansys Mesh module with hexahedral elements of 4 million cells.

I would like to know what would be the main cause of the divergence?
how can I diagnose my simulation?

I have attached the related information:
Monitoring Passive Scalar with probe.JPG



Velocity Monitoring with probe.JPG
ar.miraliakbar is offline   Reply With Quote

Old   May 27, 2022, 04:17
New Member
Join Date: Jun 2021
Posts: 1
Rep Power: 0
numb is on a distinguished road
I also calculated this case.
I think your case is diverged, maybe is fvScheme is wrong. I try some divSchemes for div(phi,C), some Schemes is convergent, but some others is diverging.
And I also have several other problems.
1. I'm confused with the inlet velocity. Only the Re number at the inlet is proposed in this paper. The velocity is missing. I also calculated the inlet velocity is 0.606 with the method in this paper. But the max velocity and is the degree of mixing totally different with this paper if u=0.606. And I change some other velocity and try to find the correct inlet velocity. And I find if the velocity at the inlet is 0.81, then the max velocity and the degree of mixing is same with the result in the paper. The result is same as the paper when Re = 9. But when I change the Re to 1, the degree of mixing in my simulation is only about 0.6 rather 0.2 in this paper.
2. And I tried some other divSchemes, but the scalar field can't be limited between 0 to 1 although I have added "bounded" before the Scheme. The scalar field in my simulation is about -0.3 ~ 1.3.
Have others simulated this case? Or someone know how to limited the scalar field?
The first figure is from the paper, and the second one is my simulation result when I set the inlet velocity as 0.81m/s

Last edited by numb; May 27, 2022 at 05:20.
numb is offline   Reply With Quote


microchannel, micromixer, mixing fluid, passive scalar, simplefoam convergence

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Plotting total passive scalar flow Lek yang STAR-CCM+ 0 July 19, 2021 03:40
Passive Scalar Transport, Diffusivity and PimpleFoam Talerbla OpenFOAM Running, Solving & CFD 1 May 6, 2021 09:02
Error in New Surface reaction model (Having multiple reactions) surajkvs OpenFOAM Programming & Development 1 April 12, 2021 05:03
How to add a passive scalar to the simpleFoam solver Uyan OpenFOAM Running, Solving & CFD 6 June 22, 2017 17:36
Passive Scalar Interaction Coefficient ArifArif STAR-CCM+ 0 July 3, 2014 04:27

All times are GMT -4. The time now is 18:47.