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

May 20, 2021 at 9:10 am #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 3 years, 2 months ago by Marc.