Power balance output and tburn consistency #17

Closed jonmaddock closed 1 year ago

In GitLab by @mkovari on May 13, 2014, 10:04

James has put my changes on the dev_output_mods branch. There is a merge request but I am creating an issue since that is our preferred method.

The power balances don’t work because of the TF coil heating error – which I understand Pete is in the process of correcting. The ITER input files have enormous TF coil heating, because they don’t have a blanket to act as shield. The power balances will change slightly when Chris’s changes are implemented.

9/5/14. Removed constraint 15 (tburn consistency).

This constraint was intended to ensure that the value of tburn converges. I note, however, that constraint 15 is not actually implemented in any of the three pulsed machine test files in the test suite. Moreover, this method means that the value of tburn will change with repeated invocations of 'caller', so will not strictly be a function of the iteration variables. This could pose a problem, especially if the constraints or figure of merit depend on tburn.

To ensure that tburn converges I have added an extra loop in fcnvm1 (evaluators.f90).

I believe that dnbeam, radial build and beta consistency equations are OK as they stand. In each case the two quantities that are constrained to be equal are quite separate and independent. This was not the case with tburn, where tburn0 was reset to equal tburn each time the physics routine was called.

Tested on steady-state and pulsed input files.

9/5/14. Added "Please include the following header and description in any models, presentations and papers based on these results." Increased length of quoted comments to 120 characters.

9/5/14. Added new power balance outputs to plant_power.f90. Also added subroutine olong to output.f90 to allow long lines in the output. (Need to discuss whether 80, 100, or 120 characters is really best.)

( see K:\Power Plant Physics and Technology\PROCESS\Output improvements.docx) Thanks, Michael

In GitLab by @pknight on May 19, 2014, 10:09

Removed tburn consistency equation. Also forced tburn to be non-negative (issue #29) (r278)

In GitLab by @pknight on Jun 12, 2014, 15:45

Okay, new power flow algorithm from Chris Harrington is now in code (r287). Please check new output format (use ipowerflow=1), and the three new figures in the User Guide.

In GitLab by @mkovari on Jun 13, 2014, 08:35

I can’t get this to work. See /home/mkovari/fusion:K/Power Plant Physics and Technology/PROCESS/Solver studies/287_runs Can anyone help?

In GitLab by @mkovari on Jul 17, 2014, 08:40

Once this is all working could you set ipowerflow=1 as default?

In GitLab by @mkovari on Jul 30, 2014, 12:12

Output table for Plasma power balance across separatrix: complete.
Constraint violations in physical units: complete.
Power interception table (ipowerflow=1): complete.
Output tables for Power balance for Reactor, Electrical Power balance and Power balance for power plant: still need to be done I think.
Output_improvements.docx

In GitLab by @mkovari on Sep 5, 2014, 10:51

IMHO having these output tables would help in spotting any possible energy conservation problems (#130)

In GitLab by @pknight on Sep 10, 2014, 09:51

Final output tables are now present if ipowerflow=1.

In GitLab by @pknight on Sep 10, 2014, 10:42

@mkovari comments please

In GitLab by @pknight on Sep 16, 2014, 08:41

Example input file with ipowerflow=1 is attached.

In GitLab by @pknight on Sep 16, 2014, 08:43

...trying yet again with firefox as IE won't attach the file! IN.DAT_ipowerflow_1

In GitLab by @mkovari on Mar 13, 2015, 14:44

Milestone removed

In GitLab by @jmorris-uk on Apr 28, 2015, 13:45

In release r389 (develop) r390 (master)

ukaea / PROCESS

Power balance output and tburn consistency #17