b2 input – SOLPS-ITER 3.2.0@wg

Specifying the type of electron energy or temperature boundary condition for each segment; makes use of ENEPAR to specify additional information, as indicated

0: default, no boundary condition is applied

default

1: prescribe the value of the electron temperature, ENEPAR(,1) specifies the temperature in eV

const

2: prescribe the gradient of the electron temperature, ENEPAR(,1) specifies the temperature gradient in eV.m^-1

const

4: prescribe the value of the electron temperature, weakly a mixed boundary condition, ENEPAR(,1) specifies the temperature in eV and ENEPAR(,2) specifies the 'strength' of the boundary condition.

const

3: sheath conditions, electron energy transmission, ENEPAR(,1) specifies an additional contribution to the energy transmission coefficient in addition to that of the potential difference [the sound speed used depends on settings of MOMPAR(,ISMAIN,2)]. BCENE = 3 not adapted for WG. Replace with drift-compatible BCENE = 15.

sheath

12: sheath conditions, electron energy transmission coefficient, ENEPAR(,1) specifies an energy transmission factor, delta_e in Q_e = delta_e Γ_e T_e. Obsolete. Replace with drift-compatible BCENE = 15.

sheath

15: sheath boundary condition, recommended when using drifts (see Section C.9.4 of manual for details). Linked to using BCCON=14 and BCMOM=13 for all ion species, BCENI=15, and BCPOT=11.

sheath

drifts

5: prescribe the electron energy flux per unit area, ENEPAR(,1) specifies the energy flux density in W.m^-2.

flux

8: prescribe the total electron heat flux with constant flux density, ENEPAR(,1) specifies the energy flux in W.

flux

13: prescribe the electron energy flux per unit area proportional to temperature, ENEPAR(,1) specifies the energy flux density per temperature in W.m^-2.J^-1 (the temperature here in J).

flux

9: prescribe the decay length for the electron temperature, ENEPAR(,1) specifies the decay length in m (can also use type [19] instead).

decay_len

19: same as [9] but to be used when simultaneously setting BCCON=1 on the same boundary. Not yet available for WG.

decay_len

14: leakage option for electron energy, ENEPAR(,1) specifies the leakage factor, α in Γ_loss = α C_s, collective n_e T_e. The energy loss also includes an electrostatic term as α ENEPAR(,2) C_s,collective e Φ n_e.

leakage

22: radial leakage condition for the electron temperature. ENEPAR(,1) specifies the leakage velocity in units of the electron thermal velocity. A temperature gradient such that the diffusive flux is set to match this leakage is imposed.

leakage

10: should not be used for feedbacks anymore, not converted to WG code.

feedback

16: feedback boundary condition with constant temperature, ENEPAR(,1) specifies the power flux in W across the core boundary, ENEPAR(,2) should be something like 0.1 and specifies the strength of the feedback (if left zero, b2stbc_bc_ref_te will be used). Also see type [17] below. Available if bcene_16_style=0 (default). If bcene_16_style=1, integrated electron heat flux with constant electron temperature, summed over all core boundaries with BCENE=16.

feedback

17: feedback boundary condition with constant shared temperature for both electrons and ions, with ENEPAR(,1) + ENIPAR(,1) giving the total power flux in W across the core boundary, ENEPAR(,2) should be something like 0.1 and specifies the strength of the feedback. Recommended to replace [16] for high densities.

feedback

20: feedback boundary condition with constant temperature, as per type [16] but with a different feedback scheme. ENEPAR(,2) is the strength of the feedback (if left zero, b2stbc_bc_ref_te is used).

feedback

21: constant temperature feedback scaled by temperature. ENEPAR(,1) specifies the desired electron temperature in eV. ENEPAR(,2) is the strength of the feedback (if left zero, b2stbc_bc_ref_te is used).

feedback

6: prescribe the total electron energy flux for a constant electron temperature, ENEPAR(,1) specifies the energy flux in W (no longer supported, use types [16] or [17] instead).

_deprec

7: prescribe the electron temperature profile from the bv_te.dat file (requires b2mndr_boundary_sources.eq.1).

_future

11: not used

_future

18: fractional drop condition. Experimental. Attempts to set the guard cell temperature to be (1 - ENEPAR(,1)) times the boundary cell temperature.

_future

Specifying the type of ion energy or temperature boundary condition for each segment; makes use of ENIPAR to specify additional information, as indicated

0: default, no boundary condition is applied

default

1: prescribe the value of the ion temperature, ENIPAR(,1) specifies the temperature in eV

const

2: prescribe the gradient of the ion temperature, ENIPAR(,1) specifies the temperature gradient in eV.m^-1

const

4: prescribe the value of the ion temperature, weakly a mixed boundary condition, ENIPAR(,1) specifies the temperature in eV and ENIPAR(,2) specifies the 'strength' of the boundary condition.

const

3: sheath conditions, ion energy transmission, ENIPAR(,1) specifies the contribution to the energy transmission coefficient [the sound speed used depends on settings of MOMPAR(,ISMAIN,2)]. BCENI = 3 not adapted for WG. Replace with drift-compatible BCENI = 15.

sheath

11: sheath conditions, ion energy transmission coefficient, ENIPAR(,1) specifies an energy transmission factor, delta_i in Q_i = delta_i T_i sum_a n_a C_s,a. Not yet available for WG.

sheath

12: sheath conditions, ion energy transmission coefficient, ENIPAR(,1) specifies an energy transmission factor, delta_i in Q_i = delta_i T_i sum_a Γ_a. Obsolete. Replace with drift-compatible BCENI = 15.

sheath

15: sheath boundary condition, recommended when using drifts (see Section C.9.5 of manual for details). Linked to using BCCON=14 and BCMOM=13 for all ion species, along with BCENE=15 and BCPOT=11.

sheath

drifts

5: prescribe the ion energy flux per unit area, ENIPAR(,1) specifies the energy flux density in W.m^-2.

flux

8: prescribe the total ion heat flux with constant flux density, ENIPAR(,1) specifies the energy flux in W

flux

13: prescribe the ion energy flux per unit area proportional to temperature, ENIPAR(,1) specifies the energy flux density per temperature in W.m^-2.J^-1 (the temperature here in J).

flux

9: prescribe the decay length for the ion temperature, ENIPAR(,1) specifies the decay length in m (can also use type [19] instead)

decay_len

19: same as [9] but to be used when simultaneously setting BCCON=1 on the same boundary. Not yet available for WG.

decay_len

14: leakage option for ion energy, ENIPAR(,1) specifies the leakage factor, α in Γ_loss = α C_sT_i.

leakage

22: radial leakage condition for the ion temperature. ENIPAR(,1) specifies the leakage velocity in units of the collective ion thermal velocity. A temperature gradient such that the diffusive flux is set to match this leakage is imposed.

leakage

10: should not be used for feedbacks anymore, not converted to WG code.

feedback

16: feedback boundary condition with constant temperature, ENIPAR(,1) specifies the power flux in W across the flux surface with index 'b2stbc_type16_ref' (default=-1), ENIPAR(,2) should be something like 0.1 and specifies the strength of the feedback (if left zero, b2stbc_bc_ref_ti will be used). Also see type [17] below. Available if bceni_16_style=0 (default). If bceni_16_style=1, integrated ion heat flux with constant ion temperature, summed over all core boundaries with BCENI=16. ENIPAR(,1) specifies the power flux in W.

feedback

17: feedback boundary condition with constant shared temperature for both electrons and ions, see BCENE=17 above for description.

feedback

20: feedback boundary condition with constant temperature, as per type [16] but with a different feedback scheme. ENEPAR(,2) is the strength of the feedback (if left zero, b2stbc_bc_ref_ti is used).

feedback

25: constant temperature feedback scaled by temperature on nearby ring. ENIPAR(,1) specifies the desired ion temperature in eV. ENIPAR(,2) is the strength of the feedback (if left zero, b2stbc_bc_ref_ti will be used).

feedback

6: prescribe the total ion energy flux for a constant ion temperature, ENIPAR(,1) specifies the energy flux in W (no longer supported, use types [16] or [17] instead)

_deprec

7: prescribe the ion temperature profile from the bv_ti.dat file (requires b2mndr_boundary_sources.eq.1)

_future

18: fractional drop condition. Experimental. Attempts to set the guard cell temperature to be (1 - ENIPAR(,1)) times the boundary cell temperature.

_future

Specifying the type of electric potential or current boundary condition for each segment; makes use of POTPAR to specify additional information, as indicated

0: default, no boundary condition is applied

default

1: prescribe the value of the potential, POTPAR(,1) specifies the potential in V

const

2: prescribe the gradient of the potential, POTPAR(,1) specifies the potential gradient in V.m^-1

const

4: prescribe the value of the potential weakly a mixed boundary condition, POTPAR(,1) specifies the potential in V and POTPAR(,2) specifies the 'strength' of the boundary condition.

const

3: sheath conditions, POTPAR(,2) used for biasing [see code for details]. BCPOT = 3 not adapted for WG. Replace with drift-compatible BCPOT = 11.

sheath

11: sheath conditions, electron energy transmission, POTPAR(,2) specifies the bias potential in V. Recommended for use in cases with drifts, along with BCENE/I=15, BCCON=14, and BCMOM=13.

sheath

drifts

5: prescribe the current flux density per unit area, POTPAR(,1) specifies the electric current flux density in A.m^-2.

flux

8: prescribe the total electric current with constant flux density, POTPAR(,1) specifies the electric current in A

flux

9: prescribe the decay length for the potential, POTPAR(,1) specifies the decay length in m.

decay_len

16: constant electric potential feedback on imposed total current. The current prescribed is given by the sum of the POTPAR(IB,1) (in A) over all the BCPOT=16 boundaries. Still experimental, will not work for drift cases.

feedback

21: constant potential feedback scaled by potential on the nearby ring. POTPAR(,,1) specifies the desired potential in V. POTPAR(,,2) is the strength of the feedback (if left zero, b2stbc_bc_ref_te will be used).

feedback

12: imposes the currents due to drifts for the core boundary. Must be used in conjunction with istyle_cur_contr_on_S_and_N.eq.2

drifts

13: imposes the currents due to drifts for the private flux and main chamber boundaries. Must be used in conjunction with istyle_cur_contr_on_S_and_N.eq.2

drifts

6: prescribe the total current flux density for a constant potential [not yet available]

_future

7: prescribe the given profile of potential from the bv_po.dat file (requires b2mndr_boundary_sources.eq.1)

_future

10: should not be used for feedbacks anymore, not converted to WG code.

_future

Specifying the type of density boundary condition for each segment and species (fastest varying index is species); makes use of CONPAR to specify additional information, as indicated:

0: default, no boundary condition is applied

default

1: prescribe the value of the density, CONPAR(,,1) specifies the required density in m^-3

const

2: prescribe the gradient of the density, CONPAR(,,1) specifies the required density gradient in m^-4

const

4: prescribe the value of the density, weakly a mixed boundary condition, CONPAR(,,1) specifies the required density in m^-3 and CONPAR(,,2) specifies the 'strength' of the boundary condition.

const

21: prescribe the value of the density and add a poloidal density variation to get a solution which is as close as possible to neoclassical theory. It is recommended to use this boundary condition together with corresponding condition on ion temperature (BCENI=24). The total desired density is summed over all BCCON=21 core boundaries. CONPAR(,,1) specifies the desired density in m^-3.

const

neoclassical

3: sheath conditions, CONPAR(,,1) not used (zero gradient is used)

sheath

5: prescribe the particle flux per unit area, CONPAR(,,1) specifies the required particle flux density in m^-2.s^-1

flux

6: prescribe the total particle flux for a constant density, CONPAR(,,1) specifies the particle flux in s^-1. Not yet available for WG.

flux

8: prescribe the total particle flux with constant flux density, CONPAR(,,1) specifies the particle flux in s^-1

flux

14: sound speed velocity flux, CONPAR(,,1) is a multiplier to the outgoing sound speed C_s. To be used in conjunction with BCENE/I=15, BCPOT=11, and BCMOM=13. Recommended for cases with drifts. If CONPAR(,,2).gt.0, an additional decay-length loss is added, where CONPAR(,,2) is the decay length.

flux

drifts

9: prescribe the decay length for the density, CONPAR(,,1) specifies the gradient length in metres (should use type 15 instead when drifts are turned on)

decay_len

10: leakage option for density, recommended for cases with drifts, CONPAR(,,1) specifies the leakage factor, α in Γ_loss = α C_s n_a, CONPAR(,,2)>0 is the stabilizing factor which may be needed if drifts are activated.

leakage

drifts

15: prescribe a radial leakage velocity, CONPAR(,,1) specifies the leakage velocity in units of the local thermal velocity.

leakage

11: should not be used for feedbacks anymore, not converted to WG code.

feedback

12: should not be used for feedbacks anymore, not converted to WG code.

feedback

13: particle density to achieve specified total flux, CONPAR(,,1) is the specified flux crossing the boundary flux surface, CONPAR(,,2) is the strength of the feedback (if left zero, b2stbc_bc_ref will be used), CONPAR(,,3), when running with Eirene and the 'ionising core' switch is used, is set internally to match the re-entering flux of ionised neutrals that crossed the core boundary (one must then have 'ionising_core'.eq.-IB where IB is the boundary index).

feedback

20: constant density feedback condition, CONPAR(,,1) specifies the desired density in m^-3. Not yet available for WG.

feedback

24: constant density feedback scaled by density on the ring 'bc_type21_ref' away. CONPAR(,,1) specifies the desired density in m^-3. CONPAR(,,2) is the strength of the feedback (if left zero, b2stbc_bc_ref will be used). Not yet available for WG.

feedback

16: particle density to achieve specified total flux, used with ASTRA coupling. The total desired flux is summed over all BCCON=16 core boundaries. CONPAR(,,1) specifies the desired particle flux in s^-1. CONPAR(,,2) is the strength of the feedback (if left zero, b2stbc_bc_ref will be used).

coupling

flux

18: prescribe total main ion particle flux, used with ASTRA coupling. Not yet available for WG.

coupling

const

19: particle flux feedback boundary condition, flux is summed over neutrals and ions, for coupling with ASTRA. The total desired flux is summed over all BCCON=19 core boundaries. This boundary condition type is applied to ions in their highest ionisation stage. CONPAR(,,1) specifies the desired particle flux in s^-1. CONPAR(,,2) is the strength of the feedback (if left zero, b2stbc_bc_ref will be used).

coupling

feedback

7: prescribe the given profile of density from the bv_na.dat file (requires b2mndr_boundary_sources.eq.1)

_future

Specifying the type of parallel momentum or velocity boundary condition for each segment and species (fastest varying index is species); makes use of MOMPAR to specify additional information, as indicated

0: default, no boundary condition is applied

default

1: prescribe the value of the parallel velocity, MOMPAR(,,1) specifies the parallel velocity in m.s^-1

const

2: prescribe the gradient of the parallel velocity, MOMPAR(,,1) specifies the parallel velocity gradient in s^-1

const

4: prescribe the value of the velocity, weakly a mixed boundary condition, MOMPAR(,,1) specifies the parallel velocity in m.s^-1 and MOMPAR(,,2) specifies the 'strength' of the boundary condition.

const

12: condition from b2stbc_spb for the parallel momentum.

const

14: condition from b2stbc_spb for the parallel momentum

const

15: prescribe the value of the parallel velocity, scaled with B_average/B_local. Not yet available for WG.

const

3: sheath conditions, Mach number as input, if MOMPAR(,,2) < 0.5, then the velocity is set to exactly MOMPAR(,,1) * C_s,collective, otherwise the velocity is set to be at least MOMPAR(,,1) * C_s,a, the thermal velocity of that species. BCMOM = 3 was not adapted for WG. Replace with drift-compatible BCMOM = 13.

sheath

13: drift-compatible sheath boundary condition for the parallel momentum. To be used in conjunction with BCENE/I=15, BCPOT=11, and BCCON=14. Recommended for cases with drifts.

sheath

5: prescribe the parallel momentum flux per unit area, MOMPAR(,,1) specifies the parallel momentum flux density in N.m^-2

flux

9: prescribe the total parallel momentum flux with constant flux density, MOMPAR(,,1) specifies the parallel momentum flux in N.

flux

10: prescribe the decay length for the parallel momentum, MOMPAR(,,1) specifies the decay length in m.

decay_len

17: leakage option for parallel momentum, MOMPAR(,,1) specifies the leakage factor, α in Γ_loss = α C_s,a m_a n_a u_a

leakage

8: special : limited shear, imposes zero gradient for the Mach number. [[[Eventually intended to have MOMPAR(,,1) specify the gradient of the Mach number in m^-1]]].

special

11: Rozhansky viscosity condition for the parallel momentum, MOMPAR(,,1) is not used. Not yet available for WG.

special

6: prescribe the total parallel momentum flux for a constant parallel velocity [not yet available]

_future

7: prescribe the given profile of parallel velocity from the bv_ua.dat file (requires b2mndr_boundary_sources.eq.1)

_future

Boundary conditions

BCENE

No condition

Constant value, gradient

Sheath conditions

Flux based

Decay length

Leakage conditions

Feedback for const value

Recommended when drifts are enabled

Neoclassical theory

Coupling with other codes

Special

Deprecated

Not yet implemented

BCENI

No condition

Constant value, gradient

Sheath conditions

Flux based

Decay length

Leakage conditions

Feedback for const value

Recommended when drifts are enabled

Neoclassical theory

Coupling with other codes

Special

Deprecated

Not yet implemented

BCPOT

No condition

Constant value, gradient

Sheath conditions

Flux based

Decay length

Leakage conditions

Feedback for const value

Recommended when drifts are enabled

Neoclassical theory

Coupling with other codes

Special

Deprecated

Not yet implemented

BCCON

No condition

Constant value, gradient

Sheath conditions

Flux based

Decay length

Leakage conditions

Feedback for const value

Recommended when drifts are enabled

Neoclassical theory

Coupling with other codes

Special

Deprecated

Not yet implemented

BCMOM

No condition

Constant value, gradient

Sheath conditions

Flux based

Decay length

Leakage conditions

Feedback for const value

Recommended when drifts are enabled

Neoclassical theory

Coupling with other codes

Special

Deprecated

Not yet implemented