Suite à une suggestion d'E. Lévêque, il pourrait être intéressant d'étudier des quantités purement réelles (et non pas complexes comme les modes des shells). Un bon candidat est le flux d'énergie entre le shell (k) et le shell (k+1).
Ce flux est défini dans l'article de Kadanov et al "Scaling and dissipation in the GOY shell model", Phys. of Fluids 7 (1995)
(un papier très intéressant par ailleurs), dans l'équation (14).
Le code permet maintenant d'avoir la trace temporelle de ce flux (dans le fichier "flux.dat") si l'option DO_FLUX du fichier "parameters.h" est activée. Le fichier produit contient N-1 valeurs, mais la 1ere et la dernière (k=0 et k=N-2) soont "incommplètes" (il n'est pas possible de calculer le flux sur les "bords" du domaine en k, ou plutôt, celui-ci n'y est pas proprement défini).
Le filtrage passe-bas FIR peut éventuellement être appliqué à cette nouvelle quantité (avec le même réglage qu'auparavant: option DO_FIR dans "parameters.h").
Suite à une suggestion d'E. Lévêque, il pourrait être intéressant d'étudier des quantités purement réelles (et non pas complexes comme les modes des shells). Un bon candidat est le flux d'énergie entre le shell (k) et le shell (k+1).
Ce flux est défini dans l'article de Kadanov et al "Scaling and dissipation in the GOY shell model", Phys. of Fluids 7 (1995)
(un papier très intéressant par ailleurs), dans l'équation (14).
Le code permet maintenant d'avoir la trace temporelle de ce flux (dans le fichier "flux.dat") si l'option DO_FLUX du fichier "parameters.h" est activée. Le fichier produit contient N-1 valeurs, mais la 1ere et la dernière (k=0 et k=N-2) soont "incommplètes" (il n'est pas possible de calculer le flux sur les "bords" du domaine en k, ou plutôt, celui-ci n'y est pas proprement défini).
Le filtrage passe-bas FIR peut éventuellement être appliqué à cette nouvelle quantité (avec le même réglage qu'auparavant: option DO_FIR dans "parameters.h").