OilfieldToolsNet/LBPM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

fix a little bug on linear interpolation of two viscosities according to the local phase field

Browse Source
This commit is contained in:
zherexli 2017-08-30 22:06:01 +10:00
parent 6b7d5c4e1b
commit add17ab75a
2 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
cpu/Color.cpp
View File

@ -1230,7 +1230,7 @@ extern "C" void ScaLBL_D3Q19_ColorCollide_gen( char *ID, double *disteven, doubl
// if phi[n] = -1, tau=tau1
// if phi[n] = 1, tau=tau2
f1=phi[n];
tau=tau1 + 0.5*(1.0+f1)*(tau2-tau1);
tau=tau1 + 0.5*(1.0-f1)*(tau2-tau1);
rlx_setA = 1.f/tau;
rlx_setB = 8.f*(2.f-rlx_setA)/(8.f-rlx_setA);

2
gpu/Color.cu
View File

@ -1231,7 +1231,7 @@ __global__ void dvc_ScaLBL_D3Q19_ColorCollide_gen( char *ID, double *disteven,
// if phi[n] = -1, tau=tau1
// if phi[n] = 1, tau=tau2
f1=phi[n];
tau=tau1 + 0.5*(1.0+f1)*(tau2-tau1);
tau=tau1 + 0.5*(1.0-f1)*(tau2-tau1);
rlx_setA = 1.f/tau;
rlx_setB = 8.f*(2.f-rlx_setA)/(8.f-rlx_setA);