
By studying the coarsening dynamics of a onedimensional spin1 BoseHubbard
model in a superfluid regime, we analytically find an unconventional universal
dynamical scaling for the growth of the spin correlation length, which is
characterized by the exponential integral unlike the conventional powerlaw or
simple logarithmic behavior, and numerically confirmed with the truncated
Wigner approximation.

We review numerical studies of quantum turbulence. Quantum turbulence is
currently one of the most important problems in low temperature physics and is
actively studied for superfluid helium and atomic BoseEinstein condensates. A
key aspect of quantum turbulence is the dynamics of condensates and quantized
vortices. The dynamics of quantized vortices in superfluid helium are described
by the vortex filament model, while the dynamics of condensates are described
by the GrossPitaevskii model. Both of these models are nonlinear, and the
quantum turbulent states of interest are far from equilibrium. Hence, numerical
studies have been indispensable for studying quantum turbulence. In fact,
numerical studies have contributed in revealing the various problems of quantum
turbulence. This article reviews the recent developments in numerical studies
of quantum turbulence. We start with the motivation and the basics of quantum
turbulence and invite readers to the frontier of this research. Though there
are many important topics in the quantum turbulence of superfluid helium, this
article focuses on inhomogeneous quantum turbulence in a channel, which has
been motivated by recent visualization experiments. Atomic BoseEinstein
condensates are a modern issue in quantum turbulence, and this article reviews
a variety of topics in the quantum turbulence of condensates e.g.
twodimensional quantum turbulence, weak wave turbulence, turbulence in a
spinor condensate, $etc.$, some of which has not been addressed in superfluid
helium and paves the novel way for quantum turbulence researches. Finally we
discuss open problems.

We theoretically and numerically study Bogoliubovwave turbulence in
threedimensional atomic BoseEinstein condensates with the GrossPitaevskii
equation, investigating three spectra for the macroscopic wave function, the
density distribution, and the Bogoliubovwave distribution. In this turbulence,
Bogoliubov waves play an important role in the behavior of these spectra, so
that we call it Bogoliubovwave turbulence. In a previous study [D. Proment
\textit{et al.}, Phys. Rev. A \textbf{80}, 051603(R) (2009)], a $3/2$ power
law in the spectrum for the macroscopic wave function was suggested by using
weak wave turbulence theory, but we find that another $7/2$ power law appears
in both theoretical and numerical calculations. Furthermore, we focus on the
spectrum for the density distribution, which can be observed in experiments,
discussing the possibility of experimental observation. Through these
analytical and numerical calculations, we also demonstrate that the previously
neglected condensate dynamics induced by the Bogoliubov waves is remarkably
important.

We theoretically and numerically study spinwave turbulence in spin1
ferromagnetic spinor BoseEinstein condensates, finding direct and inverse
cascades with powerlaw behavior. To derive these power exponents analytically,
the conventional weak wave turbulence theory is applied to the spin1 spinor
GrossPitaevskii equation. Thus, we obtain the $7/3$ and $5/3$ power laws in
the transverse spin correlation function for the direct and inverse cascades,
respectively. To confirm these power laws, numerical calculations are performed
that obtain results consistent with these power laws.

The logarithmic velocity profile is the most important statistical law of
classical turbulence affected by channel walls. This paper demonstrates
numerically that the logarithmic velocity profile of a superfluid flow appears
in quantum turbulence under pure normal flow in a channel. We investigated the
configuration and dynamics of an inhomogeneous vortex tangle affected by the
walls, and found the characteristic behavior of the loglaw.

Spinsuperflow turbulence (SST) in spin1 ferromagnetic spinor BoseEinstein
condensates is theoretically and numerically studied by using the spinor
GrossPitaevskii (GP) equations. SST is turbulence in which the disturbed spin
and superfluid velocity fields are coupled. Applying the Kolmogorovtype
dimensional scaling analysis to the hydrodynamic equations of spin and velocity
fields, we theoretically find that the 5/3 and 7/3 power laws appear in
spectra of the superflow kinetic and the spindependent interaction energy,
respectively. Our numerical calculation of the GP equations confirms SST with
the coexistence of disturbed spin and superfluid velocity field with two power
laws.

We summarize the recent theoretical and numerical works on spin turbulence
(ST) in spin1 spinor BoseEinstein condensates. When the system is excited
from the ground state, it goes through hy drodynamic instability to ST in
which the spin density vector has various disordered direction. The properties
of ST depend on whether the spindependent interaction is ferromagnetic or
antiferro magnetic. ST has some characteristics different from other kinds of
turbulence in quantum fluids. Firstly, the spectrum of the spindependent
interaction energy exhibits the characteristic power law different from the
usual Kolmogorov 5/3 law. Secondly, ST can show the spinglasslike behavior;
the spin density vectors are spatially random but temporally frozen.

We study numerically the spin turbulence in spin1 ferromagnetic spinor
BoseEinstein condensates (BECs). Spin turbulence (ST) is characterized by a
7/3 power law in the spectrum of the spindependent interaction energy. The
direction of the spin density vector is spatially disordered but temporally
frozen in ST, showing some analogy with the spin glass state. Thus, we
introduce the order parameter of spin glass into ST in spinor BECs. When ST
develops through some instability, the order parameter grows with a 7/3 power
law, thus succeeding in describing ST well.

We theoretically and numerically study spin turbulence (ST) with small spin
magnitude in spin1 spinor BoseEinstein condensates by using the spin1 spinor
GrossPitaevskii (GP)equations. This kind of ST is realized in two cases: (i)
with antiferromagnetic interaction and (ii) with ferromagnetic interaction
under a static magnetic field. The ST with small spin magnitude can exhibit two
characteristic power laws in the spectrum of the spindependent interaction
energy: 1 and 7/3 power laws in the low and highwavenumber regions,
respectively. These power laws are derived from a Kolmogorovtype dimensional
scaling analysis for the equations of motion of the spin vector and nematic
tensor. To confirm these power laws, we perform a numerical calculation of the
spin1 spinor GP equations in a twodimensional uniform system. In case (i),
the 7/3 power law appears in the highwavenumber region, but the spectrum in
the lowwavenumber region deviates from the 1 power law. In contrast, both 1
and 7/3 power laws are found to clearly appear in case (ii).

We numerically study spin turbulence in a twodimensional trapped spin1
spinor BoseEinstein condensate, focusing on the energy spectrum. The spin
turbulence in the trapped system is generated by instability of the helical
structure of the spin density vector in the initial state. Our numerical
calculation finds that in the trapped system the spectrum of the spindependent
interaction energy in the ferromagnetic case exhibits a 7/3 power law, which
was confirmed in a uniform system by our previous study. The relation between
the 7/3 power law and the motion of the spin density vector is discussed by
investigating the orbits of dynamical variables in the spin space.

We theoretically study counterflow instability and turbulence in a spin1
spinor BoseEinstein condensate by the GrossPitaevskii equation and the
Bogoliubovde Gennes equation. Our study considers (i) the dynamics induced by
the counterflow of two components with different magnetic quantum numbers,
which leads to turbulence with spin degrees of freedom, and (ii) the properties
of the turbulence. For (i), the behavior of the condensate induced by the
counterflow strongly depends on whether the spindependent interaction is
ferromagnetic or antiferromagnetic, leading to different behaviors for the
dispersion relation and the spin density vector, $etc$. For (ii), we
numerically calculate the spectrum of the spindependent interaction energy,
which also depends on the spindependent interaction. The spectrum of the
spindependent interaction energy in the ferromagnetic case exhibits a 7/3
power law, whereas that in the antiferromagnetic case does not. The 7/3 power
law can be explained by scaling analysis.

We consider a trapped atomic BoseEinstein condensate penetrated by a
repulsive Gaussian potential and theoretically investigate the dynamics induced
by oscillating the Gaussian potential. Our study is based on the numerical
calculation of the twodimensional GrossPitaevskii equation. Our calculation
reveals the dependence of the characteristic behavior of the condensate on the
amplitude and frequency of the oscillating potential. These dynamics are deeply
related to the nucleation and dynamics of quantized vortices and solitons. When
the potential oscillates with a large amplitude, it nucleates many vortex pairs
that move away from the potential. When the amplitude of the oscillation is
small, it nucleates solitons through annihilation of vortex pairs. We discuss
three issues concerning the nucleation of vortices. The first is the phase
diagram for the nucleation of vortices and solitons near the oscillating
potential. The second is the mechanism and critical velocity of the nucleation.
The critical velocity of the nucleation is an important issue in quantum
fluids, and we propose a new expression for the velocity containing both the
coherence length and the size of the potential. The third is the divergence of
the nucleation time, which is the time it takes for the potential to nucleate
vortices, near the critical parameters for vortex nucleation.

The hydrodynamics of quantized vortices and solitons in an atomic
BoseEinstein condensate excited by an oscillating potential are studied by
numerically solving the twodimensional GrossPitaevskii equation. The
oscillating potential keeps nucleating vortex dipoles, whose impulses
alternatively change their direction synchronously with the oscillation of the
potential. This leads to synergy dynamics of vortices and solitons which have
not been previously reported in quantum fluids.

We study numerically the dynamics of quantized vortices and solitons induced
by an oscillating potential inside a trapped BoseEinstein condensate. The
dynamics of the topological defects is much different from the case for a
linear uniform object; the metamorphosis between vortices and solitons is
characteristic of the dynamics. We discuss how vortices are nucleated by an
oscillating potential.