|
XGCa
|
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
|
XGCa
|
#include <sml.hpp>
Public Member Functions | |
| Simulation (NLReader::NamelistReader &nlr, double transit_time_in=1.0, double psi_norm=1.0) | |
Public Attributes | |
| int | bounce |
| Bounce routine switch 0 for off, 1 for inner boundary, 2 for both boundaries. More... | |
| int | bounce_zero_weight |
| If ==1 and bounce>0, set particle weights to zero after bouncing from the outer boundary. More... | |
| double | bounce_buffer |
| Buffer width between sml_outpsi and where the particle actually bounces (must be >=0) More... | |
| double | dt |
| Time step (s) More... | |
| double | transit_time |
| Torodial transit time of an ion with the characteristic energy. More... | |
| int | sheath_mode |
| Type of sheath (0 is none) More... | |
| bool | sheath_adjust |
| Whether to adjust the sheath. More... | |
| bool | ptb_3db_on |
| Toggle for ptb_3db. More... | |
| bool | drift_on |
| Toggle for using drift. More... | |
| bool | neutrals |
| Toggle for using neutrals. More... | |
| bool | diag_heat_on |
| Use heat diagnostics in sheath. More... | |
| bool | xz_00_up |
| Use upper region of xz when obtaining 00 mode. More... | |
| int | mstep |
| Max number of time steps. More... | |
| int | nrk |
| Order of Runge-Kutta time integration of particles + fields. More... | |
| bool | electron_on |
| Use kinetic electrons. More... | |
| bool | limit_marker_den |
| Whether to limit marker density. More... | |
| int | nthreads |
| Number of OMP threads on the host. More... | |
| bool | field_solver |
| If false, charge deposition and field solve is skipped. More... | |
| bool | helmholtz_spectral |
| Whether to solve Helmholtz-type equations with toroidally spectral solver. More... | |
| bool | thermal_bath_on |
| Switch on thermal-bath and coarse-graining operation, see ccm_param parameters. More... | |
| bool | split_weight_scheme |
| Whether to use the split-weight scheme. More... | |
| bool | dwdt_fix_bg |
| bool | dwdt_exb_only |
| bool | em_mixed_variable |
| Switch for use of mixed-variable formulation. More... | |
| bool | em_control_variate |
| Switch for use of control variate method. More... | |
| int | em_pullback_mode |
| double | em_pullback_dampfac |
| Damping term gamma on -b.grad(phi) in pullback mode 4. More... | |
| bool | em_control_variate_final_cv |
| bool | em_control_variate_flag |
| int | em_control_variate_niter |
| Number of iterations for Ampere solve with control-variate method. More... | |
| bool | em_b_para_eff |
| Effective dB_|| via a modified grad-B drift (Joiner et al., PoP, 2010) More... | |
| bool | em_es_step |
| bool | em_dAsdt_filter_on |
| Switch for applying Fourier filters on RHS of dAs/dt equation (pullback mode 4) More... | |
| bool | no_turb |
| Set all non-axisymmetric field perturbations to zero (electromagnetic version only) More... | |
| bool | em_n0 |
| Include n=0 electromagnetic mode. More... | |
| bool | em_dAsdt_hypvis |
| Use radial hyperviscosity in dA_s/dt (push_As) More... | |
| bool | exclude_private |
| Whether to exclude private region in ion charge deposition. More... | |
| bool | grad_psitheta |
| bool | iter_solver |
| int | poisson_solver_type |
| bool | update_poisson_solver |
| Whether poisson solver is updated. More... | |
| bool | iter_solver_converge |
| bool | positive_phi00_sol |
| bool | heuristic_priv_pot |
| bool | poisson_natural_boundary |
| int | add_pot0 |
| bool | poisson_use_bc |
| (XGCa only) False: 00-bc is phi=0 everywhere; true: 00-boundary More... | |
| bool | dpot_bd_apply |
| Damp (n=0,m>0) potential towards the magnetic axis. More... | |
| double | dpot_bd_width |
| Decay length (normalized flux) for (n=0,m>0) potential towards the magnetic axis. More... | |
| bool | poisson_bias |
| Whether to use a (0,0) bias potential on top of phi_00. More... | |
| bool | poisson_adia_wall |
| EXPERIMENTAL! DON'T USE UNLESS ADVISED BY AN EXPERT. More... | |
| double | dpot_te_limit |
| Max absolute value of dpot/temp in getf0. More... | |
| bool | em_use_dpot_te_limit |
| In EM simulation: whether to force usage of the min-max limiter on the turbulent potential fluctuation. More... | |
| bool | no_fp_in_f |
| double | f0_grid_alpha |
| bool | f0_update_analytic |
| bool | f0_update_analytic_local |
| If .false. –> flux-surface average update, .true. –> local. More... | |
| double | f0_update_analytic_alpha |
| Separate alpha from sml_f0_grid_alpha for Maxwellian contribution. More... | |
| double | f0_update_analytic_damp_width |
| int | f0_grid_alpha_start |
| bool | ignore_f0g |
| bool | symmetric_f |
| Enforce axisymmetry of the total distribution function on the grid. More... | |
| double | update_g_alpha |
| bool | diff_on |
| bool | resamp_on |
| Do resampling. More... | |
| bool | resamp_for_final_restart_write |
| Perform resampling before dumping the final restart file. More... | |
| bool | resamp_restart_read |
| Whether to read a restart file written from a simulation with different grid. More... | |
| double | time |
| Current simulation time. More... | |
| int | gstep |
| Current time step. More... | |
| bool | multirate_timestepping |
| Use multirate timestepping. More... | |
| int | mr_factor |
| double | mr_psi_max |
| maximum normalized psi value of each multirate region More... | |
| int | mstep_rmp_es |
| int | mstep_rmp_em |
| int | rmp_es_to_em_dt_ratio |
| Ratio of EM to ES time step size in RMP penetration calculation. More... | |
| bool | use_em_bounds |
| bool | em_exclude_private |
| double | bd_ext_delta_Ai |
| double | bd_ext_delta_Ao |
| bool | cce_coupling_on |
| Core-edge coupling flag. More... | |
| bool | use_unfused_electron_push_kernel |
| bool | loop_voltage_on |
| Inductive current drive: loop voltage (from Faraday's law curl(E)=-dB/dt) More... | |
| bool | loop_voltage_fsa |
| Make loop voltage and current drive flux-functions if true (default: true) More... | |
| double | loop_voltage |
| double | loop_voltage_psimax |
| Inductive current drive: outer boundary (in pol. flux) of the loop voltage. More... | |
| bool | current_drive_on |
| double | current_drive_p |
| Current drive PID controller: factor for proportional controller ( \(P\)) More... | |
| double | current_drive_i |
| Current drive PID controller: factor for integral controller ( \(I\)) More... | |
| double | current_drive_d |
| Current drive PID controller: factor for derivative controller ( \(D\)) More... | |
| int | mode_select_mode |
| double | mode_select_inpsi |
| double | mode_select_outpsi |
| double | mode_select_bd_width |
Static Public Attributes | |
| static constexpr bool | is_XGCa = true |
| Equivalent to the preprocessor flag for now. More... | |
| static constexpr bool | reduced_deltaf = false |
| Equivalent to the preprocessor flag for now. More... | |
| static constexpr bool | separate_n0 = true |
| Equivalent to the preprocessor flag for now. More... | |
| static constexpr bool | explicit_electromagnetic = false |
| Equivalent to the preprocessor flag for now. More... | |
| static constexpr bool | old_f0_update_analytic = false |
| Equivalent to the preprocessor flag for now. More... | |
| static constexpr int | nhybrid = 2 |
| Number of iterations in electrostatic electron weight evolution scheme. More... | |
| Simulation< Device >::Simulation | ( | NLReader::NamelistReader & | nlr, |
| double | transit_time_in = 1.0, |
||
| double | psi_norm = 1.0 |
||
| ) |
Constructor for simulation
< Fourier filter mode: 0) No filter,
< Whether to solve for fields. Scatter and solve are skipped if set to false
< For sml_poisson_solver_type=0: Split Poisson equation into axisymmetric and
| int Simulation< Device >::add_pot0 |
Set the potential on the solver boundary equal to the flux-surface
!! averaged potential (implies \(\delta\phi_wall=0\)). !! This is supposed to prevent a sheath-like potential drop at the mesh !! boundary and make sure that the logical sheath model takes care of !! the sheath physics. Also, it should prevent \(\delta\phi/\phi\sim1\). Additional electrostatic potential: (0) for off, (1) for reading data from file !! (2) for simple neoclassical axisymmetric field (not operational)
| double Simulation< Device >::bd_ext_delta_Ai |
| double Simulation< Device >::bd_ext_delta_Ao |
| int Simulation< Device >::bounce |
Bounce routine switch 0 for off, 1 for inner boundary, 2 for both boundaries.
| double Simulation< Device >::bounce_buffer |
Buffer width between sml_outpsi and where the particle actually bounces (must be >=0)
| int Simulation< Device >::bounce_zero_weight |
If ==1 and bounce>0, set particle weights to zero after bouncing from the outer boundary.
| bool Simulation< Device >::cce_coupling_on |
Core-edge coupling flag.
| double Simulation< Device >::current_drive_d |
Current drive PID controller: factor for derivative controller ( \(D\))
| double Simulation< Device >::current_drive_i |
Current drive PID controller: factor for integral controller ( \(I\))
| bool Simulation< Device >::current_drive_on |
Whether to use dynamic current drive via adjustments of the loop voltage ! The controller output is implemented as ! \(V_{loop} = 2\pi R \eta_\parallel P \left[\Delta j_\parallel + I \int_0^t \Delta j_\parallel \mathrm{d}\tau + D \frac{\partial \Delta j_\parallel}{\partial t} \right] \) ! Time integral and derivative are with respect to the ion time step instead of time, so \(1/I\) is the integration time scale, ! and \(D\) is the derivative time scale
| double Simulation< Device >::current_drive_p |
Current drive PID controller: factor for proportional controller ( \(P\))
| bool Simulation< Device >::diag_heat_on |
Use heat diagnostics in sheath.
| bool Simulation< Device >::diff_on |
| bool Simulation< Device >::dpot_bd_apply |
Damp (n=0,m>0) potential towards the magnetic axis.
| double Simulation< Device >::dpot_bd_width |
Decay length (normalized flux) for (n=0,m>0) potential towards the magnetic axis.
| double Simulation< Device >::dpot_te_limit |
Max absolute value of dpot/temp in getf0.
| bool Simulation< Device >::drift_on |
Toggle for using drift.
| double Simulation< Device >::dt |
Time step (s)
| bool Simulation< Device >::dwdt_exb_only |
| bool Simulation< Device >::dwdt_fix_bg |
| bool Simulation< Device >::electron_on |
Use kinetic electrons.
| bool Simulation< Device >::em_b_para_eff |
Effective dB_|| via a modified grad-B drift (Joiner et al., PoP, 2010)
| bool Simulation< Device >::em_control_variate |
Switch for use of control variate method.
| bool Simulation< Device >::em_control_variate_final_cv |
| bool Simulation< Device >::em_control_variate_flag |
| int Simulation< Device >::em_control_variate_niter |
Number of iterations for Ampere solve with control-variate method.
| bool Simulation< Device >::em_dAsdt_filter_on |
Switch for applying Fourier filters on RHS of dAs/dt equation (pullback mode 4)
| bool Simulation< Device >::em_dAsdt_hypvis |
Use radial hyperviscosity in dA_s/dt (push_As)
| bool Simulation< Device >::em_es_step |
If true, the ES algorithm will be used within the EM code (details in XGC_core/module.F90)
| bool Simulation< Device >::em_exclude_private |
| bool Simulation< Device >::em_mixed_variable |
Switch for use of mixed-variable formulation.
| bool Simulation< Device >::em_n0 |
Include n=0 electromagnetic mode.
| double Simulation< Device >::em_pullback_dampfac |
Damping term gamma on -b.grad(phi) in pullback mode 4.
| int Simulation< Device >::em_pullback_mode |
(0) ??? (1) ??? (2) ??? (3) mixed-variable pullback with dA_s/dt=0, (4) pullback using Ohm's law dA_s/dt + b.grad(phi) = 0 (5) pullback using dA_s/dt + b.grad(phi) = -eta*j
| bool Simulation< Device >::em_use_dpot_te_limit |
In EM simulation: whether to force usage of the min-max limiter on the turbulent potential fluctuation.
| bool Simulation< Device >::exclude_private |
Whether to exclude private region in ion charge deposition.
|
static |
Equivalent to the preprocessor flag for now.
| double Simulation< Device >::f0_grid_alpha |
Fraction of the particle weight that is transferred to the background distribution function (particle noise control) every sml_f_source_period time steps. As a rule of thumb, sml_f0_grid_alpha=sml_dt*sml_f_source_period can be used.
| int Simulation< Device >::f0_grid_alpha_start |
For delayed onset of total-f particle noise control. (must be >=1) Values >1 cause a linear ramp-up of sml_f0_grid_alpha from 0 to its input value from time step 1 to sml_f0_grid_alpha_start.
| bool Simulation< Device >::f0_update_analytic |
Switch on/off the update of the analytic part of the distribution function
| double Simulation< Device >::f0_update_analytic_alpha |
Separate alpha from sml_f0_grid_alpha for Maxwellian contribution.
| double Simulation< Device >::f0_update_analytic_damp_width |
For width of \(\exp^{-(x/w)^2}\) damping factor for updating density/temperature of the analytical f0
| bool Simulation< Device >::f0_update_analytic_local |
If .false. –> flux-surface average update, .true. –> local.
| bool Simulation< Device >::field_solver |
If false, charge deposition and field solve is skipped.
| bool Simulation< Device >::grad_psitheta |
If true, gradiant operator is calculated as \((\hat{\boldsymbol{\psi}}\cdot\nabla,\hat{\boldsymbol{\theta}}^\ast)\cdot\nabla\) instead of \((R,Z)\) coordinates.
| int Simulation< Device >::gstep |
Current time step.
| bool Simulation< Device >::helmholtz_spectral |
Whether to solve Helmholtz-type equations with toroidally spectral solver.
| bool Simulation< Device >::heuristic_priv_pot |
Override the Poisson solver in the private region and replace the solution by !! the separatrix potential scaled with the initial electron temperature
| bool Simulation< Device >::ignore_f0g |
Ignore f0g in f0 calculation. When true, it skips the update of den_f0_h, instead using the initialized value (which is zero). This is helpful to density conservation of heating when f0g is noisy due to insufficient number of particles.
|
static |
Equivalent to the preprocessor flag for now.
| bool Simulation< Device >::iter_solver |
For sml_poisson_solver_type=0: Split Poisson equation into axisymmetric and !! non-axisymmetric parts. The axisymmetric part has the !! flux-surface averaged potential on the RHS. Solve !! iteratively by lagging the RHS potential by one !! iteration. Stop after sml_iter_solver_niter iterations. !! If .false., split Poisson equation into !! flux-surface averaged part and rest (assuming that the !! polarization term commutes with the flux-surface !! average) and the rest; solve in two steps.
| bool Simulation< Device >::iter_solver_converge |
If .true., keep taking iterations in PETSc until termination criteria is met, !! i.e. until a residual tolerance is reached or it is determined that the iterations !! have diverged or exceeded sml_iter_solver_max_it. If .false., take a fixed number !! of iterations set from sml_iter_solver_niter.
| bool Simulation< Device >::limit_marker_den |
Whether to limit marker density.
| double Simulation< Device >::loop_voltage |
Initial value for the loop-voltage ! Positive loop voltage corresponds to an electron torque opposite to the toroidal ! direction (clockwise looking down on the torus, i.e. in negative Z-direction).
| bool Simulation< Device >::loop_voltage_fsa |
Make loop voltage and current drive flux-functions if true (default: true)
| bool Simulation< Device >::loop_voltage_on |
Inductive current drive: loop voltage (from Faraday's law curl(E)=-dB/dt)
| double Simulation< Device >::loop_voltage_psimax |
Inductive current drive: outer boundary (in pol. flux) of the loop voltage.
| double Simulation< Device >::mode_select_bd_width |
| double Simulation< Device >::mode_select_inpsi |
| int Simulation< Device >::mode_select_mode |
| double Simulation< Device >::mode_select_outpsi |
| int Simulation< Device >::mr_factor |
multirate timestepping for core region ions; mr_ratio (> 1) is the
factor used to accelerate the core region physics.
| double Simulation< Device >::mr_psi_max |
maximum normalized psi value of each multirate region
| int Simulation< Device >::mstep |
Max number of time steps.
| int Simulation< Device >::mstep_rmp_em |
Number of electromagnetic time steps for RMP penetration calculation with damped Newton iteration
| int Simulation< Device >::mstep_rmp_es |
Number of electrostatic time steps for RMP penetration calculation with damped Newton iteration
| bool Simulation< Device >::multirate_timestepping |
Use multirate timestepping.
| bool Simulation< Device >::neutrals |
Toggle for using neutrals.
|
static |
Number of iterations in electrostatic electron weight evolution scheme.
| bool Simulation< Device >::no_fp_in_f |
If .true. the distribution function used for the source routines will not include particle information (only for testing)
| bool Simulation< Device >::no_turb |
Set all non-axisymmetric field perturbations to zero (electromagnetic version only)
| int Simulation< Device >::nrk |
Order of Runge-Kutta time integration of particles + fields.
| int Simulation< Device >::nthreads |
Number of OMP threads on the host.
|
static |
Equivalent to the preprocessor flag for now.
| bool Simulation< Device >::poisson_adia_wall |
EXPERIMENTAL! DON'T USE UNLESS ADVISED BY AN EXPERT.
| bool Simulation< Device >::poisson_bias |
Whether to use a (0,0) bias potential on top of phi_00.
| bool Simulation< Device >::poisson_natural_boundary |
| int Simulation< Device >::poisson_solver_type |
(0) Simple axisymmetric Poisson solver (single equation) !! (–> see sml_iter_solver) !! (1) Linear Poisson equation for the axisymmetric !! potential (as in sml_iter_solver=.true.) but !! with constraint equation for the flux-surface averaged !! potential on the RHS. Solved with Schur complement solver. !! Much faster and accurate than the iterative solver, but !! causes instability in simulations with open field-lines. !! (2) Non-linear Poisson equation
| bool Simulation< Device >::poisson_use_bc |
(XGCa only) False: 00-bc is phi=0 everywhere; true: 00-boundary
| bool Simulation< Device >::positive_phi00_sol |
EXPERIMENTAL! DON'T USE UNLESS ADVISED BY AN EXPERT !! If true, makes sure that \(\langle\phi\rangle > 0\) in SOL and !! private flux region by adding a constant to \(\phi\)
| bool Simulation< Device >::ptb_3db_on |
Toggle for ptb_3db.
|
static |
Equivalent to the preprocessor flag for now.
| bool Simulation< Device >::resamp_for_final_restart_write |
Perform resampling before dumping the final restart file.
| bool Simulation< Device >::resamp_on |
Do resampling.
| bool Simulation< Device >::resamp_restart_read |
Whether to read a restart file written from a simulation with different grid.
| int Simulation< Device >::rmp_es_to_em_dt_ratio |
Ratio of EM to ES time step size in RMP penetration calculation.
|
static |
Equivalent to the preprocessor flag for now.
| bool Simulation< Device >::sheath_adjust |
Whether to adjust the sheath.
| int Simulation< Device >::sheath_mode |
Type of sheath (0 is none)
| bool Simulation< Device >::split_weight_scheme |
Whether to use the split-weight scheme.
| bool Simulation< Device >::symmetric_f |
Enforce axisymmetry of the total distribution function on the grid.
| bool Simulation< Device >::thermal_bath_on |
Switch on thermal-bath and coarse-graining operation, see ccm_param parameters.
| double Simulation< Device >::time |
Current simulation time.
| double Simulation< Device >::transit_time |
Torodial transit time of an ion with the characteristic energy.
| double Simulation< Device >::update_g_alpha |
Fraction of the numerical marker particle density that is transferred to the g2
| bool Simulation< Device >::update_poisson_solver |
Whether poisson solver is updated.
| bool Simulation< Device >::use_em_bounds |
| bool Simulation< Device >::use_unfused_electron_push_kernel |
| bool Simulation< Device >::xz_00_up |
Use upper region of xz when obtaining 00 mode.
1.8.5