search

ukaea / PROCESS

PROCESS is a systems code at UKAEA that calculates in a self-consistent manner the parameters of a fusion power plant with a specified performance, ensuring that its operating limits are not violated, and with the option to optimise to a given function of these parameters.
https://ukaea.github.io/PROCESS/
MIT License
36 stars 11 forks source link

Error in ST Divertor Model #698

Closed jonmaddock closed 1 year ago

jonmaddock commented 6 years ago

In GitLab by @stuartmuldrew on May 10, 2018, 12:00

Below is an extract of where the divertor area is calculated for spherical tokamaks. The divertor is assumed to be triangular in shape and the area of the three plates is calculated. For the diagonal plate, the area is given as a3=a2/cos(theta) but by my working, and checked by @kellis, this should be a3=a2/(cos(theta))**2.

Additionally it is not entirely clear why r1 and r2 are set to what they are.

!  Thickness of centrepost + first wall at divertor height
r1 = rmajor - rminor*triang - 3.0D0*scrapli
!  Outer radius of divertor region
r2 = rmajor + rminor
!  Angle of diagonal divertor plate from horizontal
if (vgap <= 0.0D0) then
    fdiags(1) = vgap ; call report_error(22)
end if
theta = atan(vgap/(r2-r1))
!  Vertical plate area
a1 = 2.0D0 * pi * r1 * vgap
!  Horizontal plate area
a2 = pi * (r2*r2 - r1*r1)
!  Diagonal plate area
a3 = a2 / cos(theta)
!  Total divertor area (N.B. there are two of them)
areadv = 2.0D0 * (a1+a2+a3)
hldiv = pdivt/areadv
jonmaddock commented 6 years ago

In GitLab by @stuartmuldrew on May 21, 2018, 09:33

mentioned in commit 7353772bac080404d097cdb24fb42dca67532cbd

jonmaddock commented 6 years ago

In GitLab by @jmorris-uk on May 22, 2018, 08:33

mentioned in commit 2812413935393d061263290372d6d90f90012150

jonmaddock commented 6 years ago

In GitLab by @stuartmuldrew on May 22, 2018, 09:15

closed