Hi all,

I was interested in simulating diffusion in a ternary mixture with palabos, but I’m having problems with numerical (in)stability. My first thought was to try and modify the multiComponent2d example included in palabos source files, however, the system i’m interested in has one component (ethyl glycol) which is much more viscous than the others (about an order of magnitude apart). Therefore when I compute

```
tau=3*nuLB+0.5;
```

I get ratios of TauMax/TauMin ranging from 2 to 10, depending on how I make my physical units dimensionless. however if I try to make this ratio small, the values of the relaxation times come very close to 0.5, which is also numerically unstable.

So far I’ve only been able to simulate a ternary mixture with the Shan-Chen model by choosing unphysical values for tau. Is there a better way to simulate this situation with palabos?