zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

ENH: propellerInfo - protection against Uref = 0

Browse Source
This commit is contained in:
Andrew Heather
2022-10-31 12:20:44 +00:00
parent 913c45afff
commit 3c214e99df
1 changed files with 11 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
src/functionObjects/forces/propellerInfo/propellerInfo.C
View File

@ -729,13 +729,23 @@ void Foam::functionObjects::propellerInfo::writeAxialWake
} }
void Foam::functionObjects::propellerInfo::writeWakeFields(const scalar URef) void Foam::functionObjects::propellerInfo::writeWakeFields(const scalar URef0)
{ {
if (!writeWakeFields_) if (!writeWakeFields_)
{ {
return; return;
} }
scalar URef = URef0;
if (mag(URef) < ROOTSMALL)
{
WarningInFunction
<< "Magnitude of reference velocity should be greater than zero"
<< endl;
URef = ROOTVSMALL;
}
// Normalised velocity // Normalised velocity
const vectorField UDisk(interpolate(U(), vector::uniform(nanValue_))()); const vectorField UDisk(interpolate(U(), vector::uniform(nanValue_))());
const vectorField UrDisk(coordSysPtr_->localVector(UDisk)); const vectorField UrDisk(coordSysPtr_->localVector(UDisk));