TARANG.lib.force.compute_force_hydro

TARANG.lib.force.compute_force_hydro#

Functions#

compute_force_hydro(para, U, univ)

Compute the hydrodynamic force for the vector field.

compute_random_forcing(para, U, univ)

compute_carati_forcing(para, U)

compute_random_forcing_old(para, U, univ)

compute_rotation(para, U)

Compute the rotational force for the vector field.

craya_to_cartesian(u1, Vkx, Vkz, k_mag, i, j)

craya_to_cartesian_old(u1, Vkx, Vkz, k_mag, i, j)

compute_taylor_green_forcing(para, U)

Computes spectral forcing for Taylor-Green vortex.

compute_toner_tu_forcing(para, U, univ)

Computes Toner-Tu active matter forcing (Local / "Actual" term).

Module Contents#

TARANG.lib.force.compute_force_hydro.compute_force_hydro(para, U, univ)#

Compute the hydrodynamic force for the vector field.

Parameters:

  • U (VectorField) – VectorField object representing the velocity field.

  • univ (UniversalArrays) – UniversalArrays object for temporary storage.

Return type:

None

TARANG.lib.force.compute_force_hydro.compute_random_forcing(para, U, univ)#
TARANG.lib.force.compute_force_hydro.compute_carati_forcing(para, U)#
TARANG.lib.force.compute_force_hydro.compute_random_forcing_old(para, U, univ)#
TARANG.lib.force.compute_force_hydro.compute_rotation(para, U)#

Compute the rotational force for the vector field.

Parameters:

U (VectorField) – VectorField object representing the velocity field.

Return type:

None

TARANG.lib.force.compute_force_hydro.craya_to_cartesian(u1, Vkx, Vkz, k_mag, i, j)#
TARANG.lib.force.compute_force_hydro.craya_to_cartesian_old(u1, Vkx, Vkz, k_mag, i, j)#
TARANG.lib.force.compute_force_hydro.compute_taylor_green_forcing(para, U)#

Computes spectral forcing for Taylor-Green vortex.

3D Forcing: F = F0 * [sin(x)cos(y)cos(z), -cos(x)sin(y)cos(z), 0] 2D Forcing: F = F0 * [sin(x)cos(y), -cos(x)sin(y)]

TARANG.lib.force.compute_force_hydro.compute_toner_tu_forcing(para, U, univ)#

Computes Toner-Tu active matter forcing (Local / “Actual” term).

Equation (Real Space):

F = (alpha - beta * |v|^2) * v

This requires computing the cubic nonlinearity in real space and transforming it to spectral space.

Contents