Natural convection

Mourabit · August 9, 2011, 6:55pm

Hi all,
I’m a user of LB2D_Prime code and I want to study the thermal LBM;
Especially, Natural convection in differentially heated enclosure ;
Unfortunately, I didn’t find good results with this code;
The difficulty is with the Boundary conditions and the initialization of the problem (in latman.c);

case IC_UNIFORM_RHO_IN
```
    {  
```
#if INAMURO_SIGMA_COMPONENT
```
      if( subs==0)  
```
```
      {       
```

        *macro_var_ptr++ = lattice->param.rho_in;

```
      }  a
```
```
      else // subs==1  
```
```
      {  
```

        *macro_var_ptr++ = lattice->param.rho_in;

```
      } 
```

Would you mind helping me and advising me How to initialize the problem and which BCs I have to choose;

Sincerely,