pysces.Vortices

class pysces.Vortices(positions=None, strengths=None)
__init__(positions=None, strengths=None)

Methods

__init__([positions, strengths])
append(position, strength)
induced_velocity([x, motion]) Compute the induced velocity at the given point(s)
induced_velocity_single(x, xvort, gam) Compute velocity induced at points x by a single vortex

Attributes

circulation
core_radius
positions
strengths
induced_velocity(x=None, motion=None)

Compute the induced velocity at the given point(s)

induced_velocity_single(x, xvort, gam)

Compute velocity induced at points x by a single vortex

This method returns a vector of velocities at the points x induced by a single vortex of strength gam located at xvort.

Parameters:

x : 2d array

Locations at which to compute induced velocity. Expressed as column vectors (i.e., shape should be (n,2))

xvort : 1d array

Location of vortex (shape should be (2,))

gam : float

Strength of vortex

Notes

Induced velocity is

u_\theta = \frac{\Gamma}{2 \pi r}

where r is the distance between the point and the vortex. If this distance is less than core_radius r_0, the velocity is regularized as solid-body rotation, with

u_\theta = \frac{\Gamma r}{2\pi r_0^2}