optframework.mcpbe.mcpbe module

Solving multi-component population balance equations for agglomerating systems via MONTE CARLO method. @author: Frank Rhein, frank.rhein@kit.edu, Luzia Gramling, Institute of Mechanical Process Engineering and Mechanics

class optframework.mcpbe.mcpbe.MCPBESolver(dim=2, t_total=601, t_write=10, t_vec=None, verbose=False, load_attr=True, config_path=None, init=True)

Bases: BaseSolver

init_calc(init_Vc=True)

solve_MC_N(N=5, maxiter=100000000.0)

solve_MC(maxiter=100000000.0)

select_two_random()

select_one_random()

select_size()

calc_inter_event_time_agg()

calc_inter_event_time_break()

calc_alpha_ccm(idx1, idx2)

calc_mom_t()

combine_MC(m)

visualize_qQ_t(t_plot=None, t_pause=0.5, close_all=False, scl_a4=1, figsze=[12.8, 9.600000000000001], show_x10=False, show_x50=True, show_x90=False)

visualize_mom_t(i=0, j=0, fig=None, ax=None, clr=array([0., 0.58823529, 0.50980392, 1.]), lbl='MC')

return_distribution(comp='all', t=0, Q3_grid=None)

calc_one_agg()

calc_one_break()

produce_one_frag(particle_to_break)

calc_break_rate()

calc_break_func(num_points=1000)

check_dim_consistency()

optframework.mcpbe.mcpbe.calc_betaarray_jit(COLEVAL, a, G, X, beta0, V, V_c)

Calculate the beta array for Monte Carlo using calc_beta.

Parameters

COLEVALint

Type of beta calculation model.

aint

Total number of particles.

Gfloat

Shear rate or other related parameter.

Xarray

Particle size array.

beta0float

Base collision frequency constant.

Varray

Particle volume array.

Returns

betaarray

Beta array formatted for Monte Carlo.

optframework.mcpbe.mcpbe.calc_b_r_jit_1d(BREAKRVAL, a, G, V, pl_P1, pl_P2)

optframework.mcpbe.mcpbe.calc_alpha_ccm_jit(V, alpha_prim, idx1, idx2)

optframework.mcpbe.mcpbe.calc_inter_event_time_array(Vc, a_tot, betaarray)

optframework.mcpbe.mcpbe.select_size_jit(array)