Skip to content
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
21 commits
Select commit Hold shift + click to select a range
8444a8f
🔄 Rename whttf to m_tf_coils_total for clarity and consistency; updat…
chris-ashe Mar 28, 2025
32b83ad
🔄 Rename aiwp to a_tf_coil_wp_turn_insulation for clarity and consist…
chris-ashe Mar 28, 2025
862d6ad
🔄 Rename bmaxtf to b_tf_inboard_peak for clarity and consistency; upd…
chris-ashe Mar 28, 2025
15716a0
🔄 Rename theta_coil to rad_tf_coil_toroidal for clarity and consisten…
chris-ashe Mar 28, 2025
e60e435
🔄 Rename tfareain to a_tf_coil_inboard for clarity and consistency; u…
chris-ashe Mar 28, 2025
b56cd90
🔄 Rename h_cp_top to z_cp_top for clarity and consistency; update ref…
chris-ashe Mar 28, 2025
5f24ac2
🔄 Rename tftort to dx_tf_inboard_out_toroidal for clarity and consist…
chris-ashe Mar 28, 2025
3afa042
🔄 Rename tfc_current to c_tf_coil for clarity and consistency, update…
chris-ashe Mar 28, 2025
cab98cb
🔄 Rename coilshap to tf_coil_shape_inner for clarity and consistency;…
chris-ashe Mar 31, 2025
c4c45cf
🔄 Rename jwptf to j_tf_wp for clarity and consistency; update referen…
chris-ashe Mar 31, 2025
8c3f371
🔄 Rename rbmax to r_b_tf_inboard_peak for clarity and consistency; up…
chris-ashe Mar 31, 2025
993cabd
🔄 Rename tfind to ind_tf_coil for clarity and consistency; update ref…
chris-ashe Mar 31, 2025
6468510
🔄 Rename whtconsh to m_tf_turn_steel_conduit for clarity and consiste…
chris-ashe Mar 31, 2025
db37f48
🔄 Rename tmpcry to temp_tf_cryo for clarity and consistency; update r…
chris-ashe Mar 31, 2025
84b0784
🔄 Rename thkcas to dr_tf_nose_case for clarity and consistency; updat…
chris-ashe Mar 31, 2025
ad3add9
🔄 Rename casths to dx_tf_side_case for clarity and consistency; updat…
chris-ashe Mar 31, 2025
6fa5c2b
🔄 Rename estotft to e_tf_magnetic_stored_total for clarity and consis…
chris-ashe Mar 31, 2025
ae4b1d5
🔄 Update obsolete_vars.py to include renamed variables for clarity an…
chris-ashe Apr 8, 2025
04f1fa4
Update process/io/obsolete_vars.py
chris-ashe Apr 9, 2025
17868a0
🔄 Rename variables for clarity and consistency. update references acr…
chris-ashe Apr 9, 2025
b940ea7
Update for new iteration python file
chris-ashe Apr 10, 2025
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
Original file line number Diff line number Diff line change
Expand Up @@ -87,7 +87,7 @@ The power balance for cryogenics is detailed as in the example below. The calcu
Resistive losses in current leads (MW) (qcl/1.0D6) 2.065E-02 OP
45% allowance for heat loads in transfer lines, storage tanks etc (MW) (qmisc/1.0D6) 3.116E-02 OP
Sum = Total heat removal at cryogenic temperatures (W) (helpow/1.0D6) 1.004E-01 OP
Temperature of cryogenic components (K) (tmpcry) 4.500E+00
Temperature of cryogenic components (K) (temp_tf_cryo) 4.500E+00
Efficiency (figure of merit) of cryogenic plant is 13% of ideal Carnot v 2.028E-03 OP
Electric power for cryogenic plant (MW) (crypmw) 4.952E+01 OP
```
Expand Down
16 changes: 8 additions & 8 deletions documentation/proc-pages/eng-models/tf-coil.md
Original file line number Diff line number Diff line change
Expand Up @@ -58,7 +58,7 @@ The TF coils are assumed to be supporting each other against the net centering f
#### TF coil inboard radial size

<p style='text-align: justify;'>
Following the geometry and its parametrization presented in <em>Figure 1</em>, the TF total thickness <em>dr_tf_inboard</em> \( \left( \Delta R_\mathrm{TF} \right) \) is related with the inner and outer case radial thicknesses (<em>thkcas</em>, \( \Delta R_\mathrm{case}^\mathrm{in} \) and <em>casthi</em>, \( \Delta R_\mathrm{case}^\mathrm{out} \) respectively) and the WP radial thickness <em>dr_tf_wp</em> \(\Delta R_\mathrm{WP}\) by the following equation :
Following the geometry and its parametrization presented in <em>Figure 1</em>, the TF total thickness <em>dr_tf_inboard</em> \( \left( \Delta R_\mathrm{TF} \right) \) is related with the inner and outer case radial thicknesses (<em>dr_tf_nose_case</em>, \( \Delta R_\mathrm{case}^\mathrm{in} \) and <em>casthi</em>, \( \Delta R_\mathrm{case}^\mathrm{out} \) respectively) and the WP radial thickness <em>dr_tf_wp</em> \(\Delta R_\mathrm{WP}\) by the following equation :
</p>

$$
Expand Down Expand Up @@ -89,10 +89,10 @@ $$
Although not physically divided into pieces, three sections of the case can be considered:
</p>
- <p style='text-align: justify;'>
**The nose casing:** this section corresponds to the case separating the WP with the machine center. Due to the presence of net electromechanical centering forces, this case has a major structural purpose and is often much larger than the other sides. The nose case dimension is set by its radial thickness that the user can specify using the `thkcas` input variable (iteration variable 57).
**The nose casing:** this section corresponds to the case separating the WP with the machine center. Due to the presence of net electromechanical centering forces, this case has a major structural purpose and is often much larger than the other sides. The nose case dimension is set by its radial thickness that the user can specify using the `dr_tf_nose_case` input variable (iteration variable 57).
</p>
- <p style='text-align: justify;'>
**Sidewall casing:** this section corresponds to the lateral side of the case, separating the WP with the other vaulted coils. As in the WP geometry is generally squared, the sidewall case thickness may vary with the machine radius. For this reason, the user sets its dimensions though its minimal thickness `casths`. The user can either directly specify `casths` or define it as a fraction of the total coil thickness at the inner radius of the WP (`r_wp_inner`) with the `casths_fraction` input. If `casths_fraction` is set in the input file, the `casths` value will be overwritten.
**Sidewall casing:** this section corresponds to the lateral side of the case, separating the WP with the other vaulted coils. As in the WP geometry is generally squared, the sidewall case thickness may vary with the machine radius. For this reason, the user sets its dimensions though its minimal thickness `dx_tf_side_case`. The user can either directly specify `dx_tf_side_case` or define it as a fraction of the total coil thickness at the inner radius of the WP (`r_wp_inner`) with the `casths_fraction` input. If `casths_fraction` is set in the input file, the `dx_tf_side_case` value will be overwritten.
</p>
- <p style='text-align: justify;'>
**Plasma side casing:** this section corresponds to the case section separating the WP with the plasma. As the geometry of this section is rounded, its thickness is set by its minimal value `casthi` (user input). This parameter can also be defined as a fraction of the total TF coil thickness `dr_tf_inboard` using `casthi_fraction`. If the `casthi_fraction` parametrization is used, the `casthi` value will be overwritten.
Expand Down Expand Up @@ -268,7 +268,7 @@ turns. The number of turns can be parametrized in three different ways :
A much simpler inboard mid-plane geometry is used for resistive TF coils, as shown in <em>Figure 6</em>. The most important difference is the absence of the lateral steel casing structure. Three main sections can be distinguished:
</p>

- **The bucking cylinder:** radial thickness `thkcas` (iteration variable 57), is present to support the centering forces. Its presence is however not mandatory and can be can be removed setting TODO.
- **The bucking cylinder:** radial thickness `dr_tf_nose_case` (iteration variable 57), is present to support the centering forces. Its presence is however not mandatory and can be can be removed setting TODO.
- **The conductor area:** radial thickness `dr_tf_wp` (iteration variable 140). Ground insulation, corresponding to the dark grey area in *Figure 6* is included in this section by convention.
- **The outer cylinder:** radial thickness `casthi`. This cylinder plays no role in the structural models in PROCESS.

Expand Down Expand Up @@ -1180,7 +1180,7 @@ For `i_tf_sc_mat = 4`, important superconductor properties may be input as follo
- Upper critical field at zero temperature and strain: `bcritsc`,
- Critical temperature at zero field and strain: `tcritsc`.

The toroidal field falls off at a rate $1/R$, with the peak value occurring at the outer edge of the inboard portion of the TF coil winding pack (radius `rbmax`).
The toroidal field falls off at a rate $1/R$, with the peak value occurring at the outer edge of the inboard portion of the TF coil winding pack (radius `r_b_tf_inboard_peak`).

Three constraints are relevant to the operating current density $J_{\mbox{op}}$ in the TF coils.

Expand Down Expand Up @@ -1259,9 +1259,9 @@ Another subroutine, `tfspcall` is called outside `stfcoil` to estimate to check
| `dr_tf_inboard` | TF coil maximum radial size <br> calculated if `dr_tf_wp` is used as iteration variable | ixc = 13 | No default | m |
| `tfootfi` | Outboard/inboard TF coil thickness ratio | - | 1 | - |
| `dr_tf_wp` | Winding pack radial thickness <br> calculated if `dr_tf_inboard` is used as iteration variable. Include the ground insulation and the insertion gap. | ixc = 140 | No default | m |
| `thkcas` | Nose/inner case radial thickness | ixc = 57 | 0.3 | m |
| `casths` | Minimal sidewall casing thickness | - | - | m |
| `casths_fraction` | Minimal sidewall casing thickness as a fraction of the TF coil toroidal thickness. Overwites the `casths` input value | - | 0.03 | - |
| `dr_tf_nose_case` | Nose/inner case radial thickness | ixc = 57 | 0.3 | m |
| `dx_tf_side_case` | Minimal sidewall casing thickness | - | - | m |
| `casths_fraction` | Minimal sidewall casing thickness as a fraction of the TF coil toroidal thickness. Overwites the `dx_tf_side_case` input value | - | 0.03 | - |
| `casthi` | Minimal plasma side casing thickness | - | - | m |
| `casthi_fraction` | Minimal plasma side casing thickness as a fraction of the TF thickness (`dr_tf_inboard`). Overwites the `casthi` input value | - | 0.05 | - |
Comment thread
chris-ashe marked this conversation as resolved.
| `i_tf_case_geom` | Plasma side casing geometry option:<br> - 0 : rounder front casing (ITER) <br> - 1 : Straight casing | - | 0 | - |
Expand Down
8 changes: 4 additions & 4 deletions documentation/proc-pages/fusion-devices/stellarator.md
Original file line number Diff line number Diff line change
Expand Up @@ -200,7 +200,7 @@ This file needs to be prepared by hand or can be written automatically by the pr
Alternatively `istell = 1,2,3,4,5` allow for pre-selected stellarator machines.

![alt text](../images/stellartor_windingpack.png "Thingy")
*Figure 3: Differences of the stellarator coil cross section in PROCESS compared to the tokamak description. Note the identical `thkcas` around the cable area.*
*Figure 3: Differences of the stellarator coil cross section in PROCESS compared to the tokamak description. Note the identical `dr_tf_nose_case` around the cable area.*

The stellarator coil model[^6] uses scaling aspects based on a reference calculation of the stellarator configuration, using numerical calculations at a reference point.
Examples for these calculations are inductances, peak field calculations or stellarator forces.
Expand All @@ -209,7 +209,7 @@ The fully three-dimensional shape of the coils is assumed to be fixed, but the s

The stellarator coils are assumed to be superconducting - no resistive coil calculations are performed. The critical field at the superconductor is calculated using circular approximations for the coils in the inductance and field calculations, and the limit is enforced automatically. All superconductor materials that are available for tokamaks are also available for stellarators.

The winding pack cross-section is rectangular for the stellarator coils, rather than the two-step cross-section assumed for tokamaks. The coil thicknesses and most of the dimensions of the materials within the coil cross-section are outputs from the model, instead of being inputs as is the case for tokamaks; see the variable descriptor file for details. In addition, certain iteration variables (`dr_tf_inboard`, no. 13; `thkcas`, no. 57; `cpttf`, no. 60 and `tftort`, no. 77) should not be turned on in the input file as they are calculated self-consistently (`thkcas` is required as input); the code will stop with an error message of this is attempted.
The winding pack cross-section is rectangular for the stellarator coils, rather than the two-step cross-section assumed for tokamaks. The coil thicknesses and most of the dimensions of the materials within the coil cross-section are outputs from the model, instead of being inputs as is the case for tokamaks; see the variable descriptor file for details. In addition, certain iteration variables (`dr_tf_inboard`, no. 13; `dr_tf_nose_case`, no. 57; `cpttf`, no. 60 and `dx_tf_inboard_out_toroidal`, no. 77) should not be turned on in the input file as they are calculated self-consistently (`dr_tf_nose_case` is required as input); the code will stop with an error message of this is attempted.
The conduit insulation thickness (`thicndut`), as well as the steel thickness around each conductor (`thwcndut`) should be given as input parameters together with the dimension of the conductor area (`t_turn_tf`).


Expand All @@ -224,12 +224,12 @@ i_tf_sc_mat = 8 * Switch for superconductor material in tf coils;
sig_tf_wp_max = 4.e8 * Maximal allowable Stress level on Ground insulation for a simple stellarator coil stress module (Pa)
fcutfsu = 0.7 *Copper fraction of cable conductor (TF coils), Schauer: 900 SCU strands, 522 Copper strands. Value for 0.4 Helium
tftmp = 4.75 *Peak helium coolant temperature in TF coils and PF coils (K)
tmpcry = 4.75 * Temperature in TF coils, required for plant efficiency (K)
temp_tf_cryo = 4.75 * Temperature in TF coils, required for plant efficiency (K)
vftf = 0.3 *Coolant fraction of TF coil leg (itfsup=0) this is the same for conductor and strand!
fiooic = 0.78 *Fraction TF coil critical current to operation current (should be iteration variable!)
vdalw = 12.64 * Max voltage across tf coil during quench (kV)
tdmptf = 20 * Dump time (should be iteration variable)
thkcas = 0.1 * Thickness TF Coil case (for stellarators: Also for toroidal direction)
dr_tf_nose_case = 0.1 * Thickness TF Coil case (for stellarators: Also for toroidal direction)
t_turn_tf = 0.048 * Dimension conductor area including steel and insulation. Important parameter.
thicndut = 0.0015 * Conduit insulation thickness (one side) (m)
thwcndut = 0.006 * thickness of steel around each conductor (one side) (m)
Expand Down
2 changes: 1 addition & 1 deletion documentation/proc-pages/io/utilities.md
Original file line number Diff line number Diff line change
Expand Up @@ -309,7 +309,7 @@ This utility plots the output of a PROCESS scan. PROCESS must be run on a scan-e

**Input**: `MFILE.DAT`

**Output** `scan_var1_vs_var2.pdf` (var1 by default is `bmaxtf`, var2 specified by user)
**Output** `scan_var1_vs_var2.pdf` (var1 by default is `b_tf_inboard_peak`, var2 specified by user)

### Usage

Expand Down
Loading
Loading