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

Testing equilibrium function for D3Q7 ADE

Browse Source
This commit is contained in:
James McClure 2014-11-08 10:36:14 -05:00
parent ff16995143
commit 83c6370749
1 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
cpu/Color.cpp
View File

@ -1251,12 +1251,12 @@ extern "C" void MassColorCollideD3Q7(char *ID, double *A_even, double *A_odd, do
A_even[n] = 0.3333333333333333*na;
B_even[n] = 0.3333333333333333*nb;
// Non-Stationary equilibrium distributions
feq[0] = 0.1111111111111111*(1+3*ux);
feq[1] = 0.1111111111111111*(1-3*ux);
feq[2] = 0.1111111111111111*(1+3*uy);
feq[3] = 0.1111111111111111*(1-3*uy);
feq[4] = 0.1111111111111111*(1+3*uz);
feq[5] = 0.1111111111111111*(1-3*uz);
feq[0] = 0.1111111111111111*(1+4.5*ux);
feq[1] = 0.1111111111111111*(1-4.5*ux);
feq[2] = 0.1111111111111111*(1+4.5*uy);
feq[3] = 0.1111111111111111*(1-4.5*uy);
feq[4] = 0.1111111111111111*(1+4.5*uz);
feq[5] = 0.1111111111111111*(1-4.5*uz);
// Construction and streaming for the components
for (idx=0; idx<3; idx++){
//...............................................