Simulate guiding centers using adaptive time-step. More...

#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <omp.h>
#include <math.h>
#include "../ascot5.h"
#include "../endcond.h"
#include "../consts.h"
#include "../physlib.h"
#include "../simulate.h"
#include "../particle.h"
#include "../wall.h"
#include "../diag.h"
#include "../B_field.h"
#include "../E_field.h"
#include "../boozer.h"
#include "../mhd.h"
#include "../rfof.h"
#include "../plasma.h"
#include "simulate_gc_adaptive.h"
#include "step/step_gc_cashkarp.h"
#include "mccc/mccc.h"
#include "mccc/mccc_wiener.h"

Macros
#define	DUMMY_TIMESTEP_VAL 1.0

Functions
real	simulate_gc_adaptive_inidt (sim_data sim, particle_simd_gc p, int i)
	Calculates time step value.

void	simulate_gc_adaptive (particle_queue pq, sim_data sim)
	Simulates guiding centers using adaptive time-step.

void	recalculate_acceleration (Acceleration acc, sim_data sim, particle_simd_gc p, particle_simd_gc p0)

Detailed Description

Simulate guiding centers using adaptive time-step.

Definition in file simulate_gc_adaptive.c.

Macro Definition Documentation

◆ DUMMY_TIMESTEP_VAL

#define DUMMY_TIMESTEP_VAL 1.0

Dummy time step value

Definition at line 33 of file simulate_gc_adaptive.c.

Function Documentation

◆ simulate_gc_adaptive_inidt()

real simulate_gc_adaptive_inidt	(	sim_data *	sim,
		particle_simd_gc *	p,
		int	i )

Calculates time step value.

The returned time step is either directly user-defined, 1/100th of collision frequency or user-defined fraction of gyro-motion.

Parameters

sim	pointer to simulation data struct
p	SIMD array of markers
i	index of marker for which time step is assessed

Returns: Calculated time step

Definition at line 335 of file simulate_gc_adaptive.c.

◆ simulate_gc_adaptive()

void simulate_gc_adaptive	(	particle_queue *	pq,
		sim_data *	sim )

Simulates guiding centers using adaptive time-step.

The simulation includes:

orbit-following with Cash-Karp method
Coulomb collisions with Milstein method

The simulation is carried until all marker have met some end condition or are aborted/rejected. The final state of the markers is stored in the given marker array. Other output is stored in the diagnostic array.

The adaptive time-step is determined by integrator error tolerances as well as user-defined limits for how much marker state can change during a single time-step.

Parameters

pq	particles to be simulated
sim	simulation data

Definition at line 55 of file simulate_gc_adaptive.c.

◆ recalculate_acceleration()

void recalculate_acceleration	(	Acceleration *	acc,
		sim_data *	sim,
		particle_simd_gc *	p,
		particle_simd_gc *	p0 )

Recalculate acceleration factor.

The acceleration is updated when crossing OMP. During the first crossing, the counter for orbit time is started. For the next crossing, we check if OMP was crossed in the same direction as first. If not, the crossing is ignored and for the third crossing we check the direction again. If this is not in the same direction as the first, the counters are nullified, acceleration is set to one, and the process is started again. This way we can account both for passing and banana particles, and for the cases where collisions have changed the orbit topology.

When we have two suitable crossings, the acceleration factor is updated and the counter for the orbit time and crossings are nullified.

Definition at line 385 of file simulate_gc_adaptive.c.

Macros