Skip to content

Reply To: Output of the root-mean-square(rms) velocity in the example of channel3D

OpenLB – Open Source Lattice Boltzmann Code Forums on OpenLB General Topics Output of the root-mean-square(rms) velocity in the example of channel3D Reply To: Output of the root-mean-square(rms) velocity in the example of channel3D

#5681
Marc
Participant

Dear Xu Yang,

just check the last three entries of the SuperLatticeTimeAveragedF3D in Paraview, should be your RMS Velocity:

for (int iDim = _sData.getDataSize(); iDim < _sData.getDataSize()*2; iDim++)
if (_sDataP2.get(iCloc,input[1],input[2],input[3],(int) iDim-_sDataP2.getDataSize())/_ensembles – _sData.get(iCloc,input[1],input[2],input[3],(int) iDim-_sDataP2.getDataSize())*_sData.get(iCloc,input[1],input[2],input[3],(int) iDim-_sDataP2.getDataSize())/_ensembles/_ensembles<0) {
output[iDim]=0;
}
else {
output[iDim] = sqrt(_sDataP2.get(iCloc,input[1],input[2],input[3],(int) iDim-_sDataP2.getDataSize())/_ensembles – _sData.get(iCloc,input[1],input[2],input[3],(int) iDim-_sDataP2.getDataSize())*_sData.get(iCloc,input[1],input[2],input[3],(int) iDim-_sDataP2.getDataSize())/_ensembles/_ensembles);
}

Best,
Marc

  • This reply was modified 2 years, 11 months ago by Marc.