Skip to content

OLAF FVW stability issue when simulating multiple rotors #3298

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
luwang00 wants to merge 3 commits into
base: rc-5.0.1
Choose a base branch
from
+22 −17
Open

OLAF FVW stability issue when simulating multiple rotors #3298

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
3cfadb3
FVW: Fix notification about wake extent
luwang00 Apr 2, 2026
ebbc2f3
AD: Clarify which rotor is cavitating
luwang00 Apr 2, 2026
f588d9e
Fix FVW instability with mutiple rotors
luwang00 Apr 3, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion modules/aerodyn/src/AeroDyn.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2478,7 +2478,7 @@ subroutine RotCavtCrit(u, p, m, errStat, errMsg)
SigmaCavitCrit = ( p%rotors(iR)%Patm + ( p%rotors(iR)%Gravity * ( abs( u%rotors(iR)%BladeMotion(j)%Position(3,i) + u%rotors(iR)%BladeMotion(j)%TranslationDisp(3,i) ) + p%rotors(iR)%MSL2SWL ) * p%rotors(iR)%airDens ) - p%rotors(iR)%Pvap ) / ( 0.5_ReKi * p%rotors(iR)%airDens * Vreltemp**2 ) ! Critical value of Sigma, cavitation occurs if local cavitation number is greater than this

if ( ( SigmaCavitCrit < SigmaCavit ) .and. ( .not. ( m%rotors(iR)%CavitWarnSet(i,j) ) ) ) then
call WrScr( NewLine//'Cavitation occurred at blade '//trim(num2lstr(j))//' and node '//trim(num2lstr(i))//'.' )
call WrScr( NewLine//'Cavitation occurred at rotor '//trim(num2lstr(iR))//', blade '//trim(num2lstr(j))//', and node '//trim(num2lstr(i))//'.' )
m%rotors(iR)%CavitWarnSet(i,j) = .true.
end if

Expand Down
35 changes: 20 additions & 15 deletions modules/aerodyn/src/FVW.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1489,6 +1489,10 @@ subroutine FVW_CalcOutput(t, u, p, x, xd, z, OtherState, y, m, ErrStat, ErrMsg)
character(*), intent( out) :: ErrMsg !< Error message if ErrStat /= ErrID_None
! Local variables
integer(IntKi) :: nP
integer(IntKi) :: iStepFVW
logical, save :: FirstNotice_FreeNearWake = .true.
logical, save :: FirstNotice_NearWake = .true.
logical, save :: FirstNotice_FullWake = .true.
integer(IntKi) :: ErrStat2
character(ErrMsgLen) :: ErrMsg2
character(*), parameter :: RoutineName = 'FVW_CalcOutput'
Expand All @@ -1514,23 +1518,24 @@ subroutine FVW_CalcOutput(t, u, p, x, xd, z, OtherState, y, m, ErrStat, ErrMsg)
endif

! Write some info to screen when major milestone achieved
if (m%iStep == p%nNWFree .and. p%nNWFree<p%nNWMax) then
nP = CountCPs(p, p%nNWFree, 0)
call WrScr(NewLine//'[INFO] OLAF free near wake is at full extent - '//trim(num2lstr(t))//'s, '//trim(num2lstr(nP))//' points.')
endif
if (m%iStep == p%nNWMax) then
nP = CountCPs(p, p%nNWMax, 0)
if (p%nFWMax==0) then
call WrScr(NewLine//'[INFO] OLAF wake is at full extent - '//trim(num2lstr(t))//'s, '//trim(num2lstr(nP))//' points.')
else
call WrScr(NewLine//'[INFO] OLAF near wake is at full extent - '//trim(num2lstr(t))//'s, '//trim(num2lstr(nP))//' points.')
endif
endif
if (p%nFWMax>0 .and. m%iStep== p%nNWMax+p%nFWMax) then
iStepFVW = floor(t/p%DTfvw+1.0E-6_DbKi)
if (FirstNotice_FullWake .and. iStepFVW == p%nNWMax+p%nFWMax) then
nP = CountCPs(p, p%nNWMax, p%nFWMax)
call WrScr(NewLine//'[INFO] OLAF wake is at full extent - '//trim(num2lstr(t))//'s, '//trim(num2lstr(nP))//' points.')
endif

FirstNotice_FullWake = .false.
FirstNotice_NearWake = .false.
FirstNotice_FreeNearWake = .false.
elseif (FirstNotice_NearWake .and. iStepFVW == p%nNWMax) then ! implies p%nFWMax > 0
nP = CountCPs(p, p%nNWMax, 0_IntKi)
call WrScr(NewLine//'[INFO] OLAF near wake is at full extent - '//trim(num2lstr(t))//'s, '//trim(num2lstr(nP))//' points.')
FirstNotice_NearWake = .false.
FirstNotice_FreeNearWake = .false.
elseif (FirstNotice_FreeNearWake .and. iStepFVW == p%nNWFree) then ! implies p%nNWFree<p%nNWMax
nP = CountCPs(p, p%nNWFree, 0_IntKi)
call WrScr(NewLine//'[INFO] OLAF free near wake is at full extent - '//trim(num2lstr(t))//'s, '//trim(num2lstr(nP))//' points.')
FirstNotice_FreeNearWake = .false.
end if

! --- Export to VTK
if (m%VTKstep==-1) then
! Has never been called, special handling for init
Expand Down
2 changes: 1 addition & 1 deletion modules/openfast-library/src/FAST_Subs.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -718,7 +718,7 @@ SUBROUTINE FAST_InitializeAll( t_initial, m_Glue, p_FAST, y_FAST, m_FAST, ED, SE
if (Failed()) return

! Loop through rotors and add module for each one
do iRot = 1, p_FAST%NRotors
do iRot = p_FAST%NRotors, 1, -1
CALL MV_AddModule(m_Glue%ModData, Module_AD, 'AD', iRot, dt_module, p_FAST%DT, &
Init%OutData_AD%rotors(iRot)%Vars, p_FAST%Linearize, ErrStat2, ErrMsg2, iRotor=iRot)
if (Failed()) return
Expand Down
Loading