marge3d.numeric

TODO : documentation …

Classes

NumericalSolver

3D solver for the Maxey–Riley-Gatignol equation (MaRGE)

Module Contents

class NumericalSolver(x, w, velocity_field, Nt, order, params=None, **kwargs)

3D solver for the Maxey–Riley-Gatignol equation (MaRGE) using Daitche’s Adams–Bashforth methods.

Parameters:

• x (array_like) – Initial particle position (3,)

• w (array_like) – Initial relative velocity (3,)

• velocity_field (object) – Must implement get_velocity, get_gradient, get_dudt

• Nt (int) – Number of time steps

• order (int) – Order of Daitche method (1, 2, or 3)

• params (DaitcheParameters) – Physical parameters of the particle–fluid system

x

w

vel

order

abstractmethod solve(t_v, flag=False)

calculate_G(w1, w2, w3, x, y, z, t)

alpha_jn(j, n)

alpha_v(N)

calc_alpha_mat(N)

Euler(t_v, flag=False)

beta_jn(j, n)

beta_v(N)

calc_beta_mat(N)

AdamBashf2(t_v, flag=False)

gamma_jn(j, n)

gamma_v(N)

calc_gamma_mat(N)

AdamBashf3(t_v, flag=False)