Variables
	MODULE_ROOT = Path(__file__).resolve().parent.parent.parent

	REPO_ROOT = MODULE_ROOT.parent

	m = importlib.import_module(f'radcluster_code.py_utils.{mod_name}')

	I

	V

	i_mobile

	v_mobile

	i_discrete

	v_discrete

	saved = sys.stdout, sys.stderr

	stdout

	stderr

	sim

dict	overrides

	inp = sim.input_data

bool	placed = False

	re = sim.rate_equations

	PM = re.n_mom

	flush

dict	config

	results = sim.run(solver_config=config, save_output=False, timeout_s=60)

	t = results['t']

	y = results['y']

	i_d = re.i_discrete

	i_VAC = i_d + PM * re.I_bin

list	edges = []

int	edge = i_d + 1

int	n_lo = edge

	n_hi = max(int(edge * 1.3), edge + 1)

	mu0 = y[i_d + PM*k, -1]

int	mu1 = 2 else 0

int	n_bar = mu1 / max(abs(mu0), 1e-300)

int	bw = n_hi - n_lo

	S1 = sum(range(n_lo, n_hi))

	S2 = sum(n*n for n in range(n_lo, n_hi))

int	det = bw * S2 - S1 * S1

int	c_lo = 0

tuple	phi0_lo = (S2 - S1 * n_lo) / det

tuple	phi1_lo = (bw * n_lo - S1) / det

tuple	phi0_hi = (S2 - S1 * (n_hi-1)) / det

tuple	phi1_hi = (bw * (n_hi-1) - S1) / det

tuple	c_hi = phi0_hi * mu0 + phi1_hi * mu1

str	flag = ""

	v_d = re.v_discrete

	vac_start = i_VAC + v_d

list	v_edges = []

int	m_lo = edge

	m_hi = max(int(edge * 1.3), edge + 1)

	sia_neg = np.sum(y[:i_d, -1] < -1e-30)

	vac_neg = np.sum(y[i_VAC:i_VAC+v_d, -1] < -1e-30)

	c_boundary = y[i_d - 1, -1]

Detailed Description

Diagnose bin-moment stall: run to t=8000 (just before stall) and inspect state.

Variable Documentation

◆ bw

int diagnose_stall.bw = n_hi - n_lo

Definition at line 106 of file diagnose_stall.py.

◆ c_boundary

diagnose_stall.c_boundary = y[i_d - 1, -1]

Definition at line 183 of file diagnose_stall.py.

◆ c_hi

tuple diagnose_stall.c_hi = phi0_hi * mu0 + phi1_hi * mu1

Definition at line 119 of file diagnose_stall.py.

◆ c_lo

tuple diagnose_stall.c_lo = 0

Definition at line 112 of file diagnose_stall.py.

◆ config

dict diagnose_stall.config

Initial value:

=  {
    't_span': (1e-6, 8e3),
    'n_points': 100,
    'log_time': True,
    'rtol': 1e-5,
    'atol': 1e-20,
    'solver_method': {'linsol': 'dense'},
}

Definition at line 66 of file diagnose_stall.py.

◆ det

int diagnose_stall.det = bw * S2 - S1 * S1

Definition at line 111 of file diagnose_stall.py.

◆ edge

int diagnose_stall.edge = i_d + 1

Definition at line 93 of file diagnose_stall.py.

◆ edges

list diagnose_stall.edges = []

Definition at line 92 of file diagnose_stall.py.

◆ flag

str diagnose_stall.flag = ""

Definition at line 121 of file diagnose_stall.py.

◆ flush

diagnose_stall.flush

Definition at line 63 of file diagnose_stall.py.

◆ I

diagnose_stall.I

Definition at line 20 of file diagnose_stall.py.

◆ i_d

diagnose_stall.i_d = re.i_discrete

Definition at line 86 of file diagnose_stall.py.

◆ i_discrete

diagnose_stall.i_discrete

Definition at line 22 of file diagnose_stall.py.

◆ i_mobile

diagnose_stall.i_mobile

Definition at line 21 of file diagnose_stall.py.

◆ i_VAC

diagnose_stall.i_VAC = i_d + PM * re.I_bin

Definition at line 87 of file diagnose_stall.py.

◆ inp

diagnose_stall.inp = sim.input_data

Definition at line 42 of file diagnose_stall.py.

◆ m

diagnose_stall.m = importlib.import_module(f'radcluster_code.py_utils.{mod_name}')

Definition at line 16 of file diagnose_stall.py.

◆ m_hi

diagnose_stall.m_hi = max(int(edge * 1.3), edge + 1)

Definition at line 141 of file diagnose_stall.py.

◆ m_lo

diagnose_stall.m_lo = edge

Definition at line 140 of file diagnose_stall.py.

◆ MODULE_ROOT

diagnose_stall.MODULE_ROOT = Path(__file__).resolve().parent.parent.parent

Definition at line 7 of file diagnose_stall.py.

◆ mu0

diagnose_stall.mu0 = y[i_d + PM*k, -1]

Definition at line 102 of file diagnose_stall.py.

◆ mu1

int diagnose_stall.mu1 = 2 else 0

Definition at line 103 of file diagnose_stall.py.

◆ n_bar

int diagnose_stall.n_bar = mu1 / max(abs(mu0), 1e-300)

Definition at line 104 of file diagnose_stall.py.

◆ n_hi

diagnose_stall.n_hi = max(int(edge * 1.3), edge + 1)

Definition at line 96 of file diagnose_stall.py.

◆ n_lo

diagnose_stall.n_lo = edge

Definition at line 95 of file diagnose_stall.py.

◆ overrides

dict diagnose_stall.overrides

Initial value:

=  {
        'eta': 0.3, 'f_cl_i': 0.6, 'f_cl_v': 0.15,
        'E_m_1D': 0.4, 'i_mobile': i_mobile, 'L_hat': 71.8,
        'c_C': 1.94e-4, 'E_b_C_SIA': 0.65,
        'rho_d': 1e15, 'Z_i': 1.05, 'Z_ii': 1.2,
        'shape_function': 'linear',
        'i_discrete': i_discrete, 'v_discrete': v_discrete,
        'r_ratio': 1.3,
    }

Definition at line 33 of file diagnose_stall.py.

◆ phi0_hi

tuple diagnose_stall.phi0_hi = (S2 - S1 * (n_hi-1)) / det

Definition at line 117 of file diagnose_stall.py.

◆ phi0_lo

tuple diagnose_stall.phi0_lo = (S2 - S1 * n_lo) / det

Definition at line 114 of file diagnose_stall.py.

◆ phi1_hi

tuple diagnose_stall.phi1_hi = (bw * (n_hi-1) - S1) / det

Definition at line 118 of file diagnose_stall.py.

◆ phi1_lo

tuple diagnose_stall.phi1_lo = (bw * n_lo - S1) / det

Definition at line 115 of file diagnose_stall.py.

◆ placed

bool diagnose_stall.placed = False

Definition at line 44 of file diagnose_stall.py.

◆ PM

diagnose_stall.PM = re.n_mom

Definition at line 61 of file diagnose_stall.py.

◆ re

diagnose_stall.re = sim.rate_equations

Definition at line 60 of file diagnose_stall.py.

◆ REPO_ROOT

diagnose_stall.REPO_ROOT = MODULE_ROOT.parent

Definition at line 8 of file diagnose_stall.py.

◆ results

diagnose_stall.results = sim.run(solver_config=config, save_output=False, timeout_s=60)

Definition at line 76 of file diagnose_stall.py.

◆ S1

diagnose_stall.S1 = sum(range(n_lo, n_hi))

Definition at line 109 of file diagnose_stall.py.

◆ S2

diagnose_stall.S2 = sum(n*n for n in range(n_lo, n_hi))

Definition at line 110 of file diagnose_stall.py.

◆ saved

diagnose_stall.saved = sys.stdout, sys.stderr

Definition at line 24 of file diagnose_stall.py.

◆ sia_neg

diagnose_stall.sia_neg = np.sum(y[:i_d, -1] < -1e-30)

Definition at line 178 of file diagnose_stall.py.

◆ sim

diagnose_stall.sim

Initial value:

=  RadClusterSimulation(
        I=I, V=V, solver_mode='full_system',
        physics_option='bin_moment_CD_fission',
        C_floor=1e-20, he_kinetics='quasi_steady_state',
        i_mobile=i_mobile, v_mobile=v_mobile,
    )

Definition at line 27 of file diagnose_stall.py.

◆ stderr

diagnose_stall.stderr

Definition at line 25 of file diagnose_stall.py.

◆ stdout

diagnose_stall.stdout

Definition at line 25 of file diagnose_stall.py.

◆ t

diagnose_stall.t = results['t']

Definition at line 82 of file diagnose_stall.py.

◆ V

diagnose_stall.V

Definition at line 20 of file diagnose_stall.py.

◆ v_d

diagnose_stall.v_d = re.v_discrete

Definition at line 134 of file diagnose_stall.py.

◆ v_discrete

diagnose_stall.v_discrete

Definition at line 22 of file diagnose_stall.py.

◆ v_edges

list diagnose_stall.v_edges = []

Definition at line 137 of file diagnose_stall.py.

◆ v_mobile

diagnose_stall.v_mobile

Definition at line 21 of file diagnose_stall.py.

◆ vac_neg

diagnose_stall.vac_neg = np.sum(y[i_VAC:i_VAC+v_d, -1] < -1e-30)

Definition at line 179 of file diagnose_stall.py.

◆ vac_start

diagnose_stall.vac_start = i_VAC + v_d

Definition at line 135 of file diagnose_stall.py.

◆ y

diagnose_stall.y = results['y']

Definition at line 83 of file diagnose_stall.py.

Variables

Detailed Description

Variable Documentation

◆ bw

◆ c_boundary

◆ c_hi

◆ c_lo

◆ config

◆ det

◆ edge

◆ edges

◆ flag

◆ flush

◆ I

◆ i_d

◆ i_discrete

◆ i_mobile

◆ i_VAC

◆ inp

◆ m

◆ m_hi

◆ m_lo

◆ MODULE_ROOT

◆ mu0

◆ mu1

◆ n_bar

◆ n_hi

◆ n_lo

◆ overrides

◆ phi0_hi

◆ phi0_lo

◆ phi1_hi

◆ phi1_lo

◆ placed

◆ PM

◆ re

◆ REPO_ROOT

◆ results

◆ S1

◆ S2

◆ saved

◆ sia_neg

◆ sim

◆ stderr

◆ stdout

◆ t

◆ V

◆ v_d

◆ v_discrete

◆ v_edges

◆ v_mobile

◆ vac_neg

◆ vac_start

◆ y