|
[Sponsors] |
November 8, 2023, 14:49 |
Why there is cellV/V_ in radialActuationDisk
|
#1 |
Member
Marķa Rosales
Join Date: Mar 2023
Location: Spain
Posts: 48
Rep Power: 3 |
Hello community,
Why fvOptions radialActuationDiskSource has in its Usource line the contribution of each cell volume divided by total volume PLUS the normalized radial distribution already considered in Tr[i]? Here is the code: https://www.openfoam.com/documentati...8C_source.html HTML Code:
const scalar a = 1.0 - Cp/Ct; const scalar T = 2.0*rhoRef*diskArea_*mag(Uref)*a*(1.0 - a); forAll(cells, i) { const scalar r2 = magSqr(mesh().cellCentres()[cells[i]] - avgCentre); Tr[i] = T *(radialCoeffs_[0] + radialCoeffs_[1]*r2 + radialCoeffs_[2]*sqr(r2)) /intCoeffs; Usource[cells[i]] += ((Vcells[cells[i]]/V_)*Tr[i]*E) & Uref; } |
|
Tags |
actuator disc, actuator disk, fvoptions momentum source |
|
|