fluidsim.base.output.spatiotemporal_spectra

Spatiotemporal Spectra

Provides:

class fluidsim.base.output.spatiotemporal_spectra.SpatioTemporalSpectra3D(output)

Bases: SpecificOutput

Computes the spatiotemporal spectra.

_tag = 'spatiotemporal_spectra'

nb_dim = 3

classmethod _complete_params_with_default(params)

_init_files(arrays_1st_time=None)

_init_new_file(tmin_file=None)

Initializes a new file

_write_to_file(data)

Writes a file with the temporal data

_add_probes_data_to_dict(data, key)

Probes fields in Fourier space and append data to a dict object

_online_save()

Prepares data and writes to file

load_time_series(keys=None, tmin=0, tmax=None, dtype=None)

load time series from files

_compute_spectrum(data)

compute_spectra(tmin=0, tmax=None, dtype=None)

compute spatiotemporal spectra from files

_get_data_probe_from_field(field)

class fluidsim.base.output.spatiotemporal_spectra.SpatioTemporalSpectra2D(output)

Bases: SpatioTemporalSpectra3D

nb_dim = 2

_get_data_probe_from_field(field)

class fluidsim.base.output.spatiotemporal_spectra.SpatioTemporalSpectraNS

Bases: object

_get_path_saved_spectra(tmin, tmax, dtype, save_urud)

_get_path_saved_tspectra(tmin, tmax, dtype, save_urud)

save_spectra_kzkhomega(tmin=0, tmax=None, dtype=None, save_urud=False)

save:

• the spatiotemporal spectra, with a cylindrical average in k-space

• the temporal spectra, with an average on the whole k-space

load_spectra_kzkhomega(tmin=0, tmax=None, dtype=None, save_urud=False)

load kzkhomega spectra from file

compute_omega_emp_vs_kzkh(spectra_kzkhomega, key_spect='spectrum_b')

Compute empirical frequency and fluctuation from the spatiotemporal spectra:

\[ \begin{align}\begin{aligned}\omega_{emp}(k_h, k_z) = \frac{\int ~ \omega ~ S(k_h, k_z, \omega) ~ \mathrm{d}\omega}{\int ~ S(k_h, k_z, \omega) ~ \mathrm{d}\omega},\\\delta \omega_{emp}(k_h, k_z) = \sqrt{\frac{\int ~ (\omega - \omega_{emp})^2 ~ S(k_h, k_z, \omega) ~ \mathrm{d}\omega}{\int ~ S(k_h, k_z, \omega) ~ \mathrm{d}\omega}},\end{aligned}\end{align} \]

where \(\omega_{emp}\) is the empirical frequency and \(\delta \omega_{emp}\) is the empirical frequency fluctuation. \(S(k_h, k_z, \omega)\) is the spectra of key_spect.

plot_kzkhomega(key_field='b', tmin=0, tmax=None, dtype=None, equation=None, xmax=None, ymax=None, cmap=None, vmin=None, vmax=None, plot_omega_emp=False, linscale=False)

plot the spatiotemporal spectra, with a cylindrical average in k-space

equation must start with ‘omega=’, ‘kh=’, ‘kz=’, ‘ikh=’ or ‘ikz=’.

For 3d solvers, key_field can be in State.keys_state_phys + [“Khd”, “Khr”, “Kp”].

compute_spectra_urud(tmin=0, tmax=None, dtype=None)

compute_temporal_spectra(tmin=0, tmax=None, dtype=None, compute_urud=False)

compute the temporal spectra by averaging over Fourier space

_get_default_tmin_periodogram(tmin)

plot_temporal_spectra(key_field=None, tmin=None, tmax=None, xlim=None, ylim=None, dtype=None, xscale='log', coef_compensate=0, plot_resonant_modes=True)

plot the temporal spectra computed from the 4d spectra

load_temporal_spectra(tmin=None, tmax=None, dtype=None, save_urud=False)

load temporal spectra from file

save_temporal_spectra(tmin=0, tmax=None, dtype=None, save_urud=False)

compute temporal spectra from files

_get_default_tmax()

get_spectra(tmin=0, tmax=None, dtype=None)

Functions

filter_tmins_paths(tmin, tmins, paths)
find_index_first_g(arr, value)	find the first index such that arr[index] > value
find_index_first_geq(arr, value)	find the first index such that arr[index] >= value
find_index_first_l(arr, value)	find the first index such that arr[index] < value
get_arange_minmax(times, tmin, tmax)	get a range of index for which tmin <= times[i] <= tmax
sort_files_tmin(paths[, tmins])

Classes

SpatioTemporalSpectra2D(output)
SpatioTemporalSpectra3D(output)	Computes the spatiotemporal spectra.
SpatioTemporalSpectraNS()