fluidsim.solvers.ns2d.strat.solver

NS2D solver (fluidsim.solvers.ns2d.strat.solver)

class fluidsim.solvers.ns2d.strat.solver.Simul(params)

Bases: Simul

Pseudo-spectral solver 2D incompressible Navier-Stokes equations.

InfoSolver

alias of InfoSolverNS2DStrat

static _complete_params_with_default(params)

This static method is used to complete the params container.

tendencies_nonlin(state_spect=None, old=None)

Compute the nonlinear tendencies of the solver ns2d.strat.

Parameters:

state_spectfluidsim.base.setofvariables.SetOfVariables

optional

Array containing the state, i.e. the vorticity, in Fourier space. If state_spect, the variables vorticity and the velocity are computed from it, otherwise, they are taken from the global state of the simulation, self.state.

These two possibilities are used during the Runge-Kutta time-stepping.

Returns:

tendencies_fftfluidsim.base.setofvariables.SetOfVariables

An array containing the tendencies for the vorticity and the buoyancy perturbation.

Notes

The 2D Navier-Stokes equation can be written

\[\partial_t \hat\zeta = \hat N(\zeta) - \sigma(k) \hat \zeta,\]

and

\[\partial_t \hat b = \hat N(b) - \sigma(k) \hat b\]

This function compute the nonlinear terms (“tendencies”) \(N(\zeta) = - \mathbf{u}\cdot \mathbf{\nabla} \zeta + \mathbf{\nabla}\wedge b\mathbf{e_z} = - \mathbf{u}\cdot \mathbf{\nabla} \zeta + \partial_x b\) and \(N(b) = - \mathbf{u}\cdot \mathbf{\nabla} b + N^2u_y\).

check_energy_conservation(rot_fft, b_fft, f_rot_fft, f_b_fft)

Check energy conservation for the inviscid case.

compute_dispersion_relation()

Computes the dispersion relation of internal gravity waves solver ns2d.strat.

Returns:

omega_dispersion_relationarr

Frequency dispersion relation in rad.

Functions

tendencies_nonlin_ns2dstrat(ux, uy, px_rot, ...)

Classes

InfoSolverNS2DStrat([only_root])
Simul(params)	Pseudo-spectral solver 2D incompressible Navier-Stokes equations.