The paper proposes a numerical technique within the Lagrangian description for propagating the quantum fluid dynamical (QFD) equations in terms of the Madelung field variables R and S, which are connected to the wave function via the transformation psi = exp(R + iS)/[symbol: see text]. The technique rests on the QFD equations depending only on the form, not the magnitude, of the probability density rho = magnitude of psi 2 and on the structure of R = [symbol: see text]/2 ln rho generally being simpler and smoother than rho. The spatially smooth functions R and S are especially suitable for multivariate radial basis function interpolation to enable the implementation of a robust numerical scheme. Examples of two-dimensional model systems show that the method rivals, in both efficiency and accuracy, the split-operator and Chebychev expansion methods. The results on a three-dimensional model system indicates that the present method is superior to the existing ones, especially, for its low storage requirement and its uniform accuracy. The advantage of the new algorithm is expected to increase for higher dimensional systems to provide a practical computational tool.The paper proposes a numerical technique within the Lagrangian description for propagating the quantum fluid dynamical (QFD) equations in terms of the Madelung field variables R and S, which are connected to the wave function via the transformation psi = exp(R + iS)/[symbol: see text]. The technique rests on the QFD equations depending only on the form, not the magnitude, of the probability density rho = magnitude of psi 2 and on the structure of R = [symbol: see text]/2 ln rho generally being simpler and smoother than rho. The spatially smooth functions R and S are especially suitable for multivariate radial basis function interpolation to enable the implementation of a robust numerical scheme. Examples of two-dimensional model systems show that the method rivals, in both efficiency and accuracy, the split-operator and Chebychev expansion methods. The results on a three-dimensional model system indicates that the present method is superior to the existing ones, especially, for its low storage requirement and its uniform accuracy. The advantage of the new algorithm is expected to increase for higher dimensional systems to provide a practical computational tool.

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topicsMay 1, 200030 citations

Physics of Quantum Fluids

Alberto Bramati, Michele Modugno

Springer series in solid-state sciencesJan 1, 201344 citations

Model-independent determination of the shear viscosity of the unitary Fermi gas at high temperature

Marcus Bluhm, Thomas Schäfer

Determinations of the shear viscosity of trapped ultracold gases suffer from systematic, uncontrolled uncertainties related to the treatment of the dilute part of the gas cloud. In this work we present an analysis of expansion experiments based on a new method, anisotropic fluid dynamics, that interpolates between Navier-Stokes fluid dynamics at the center of the cloud and ballistic behavior in the dilute corona. We validate the method using a comparison between anisotropic fluid dynamics and numerical solutions of the Boltzmann equation. We then apply anisotropic fluid dynamics to the expansion data reported by Cao et al. In the high temperature limit we find $\eta=0.282(mT)^{3/2}$ , which agrees within about 5\% with the theoretical prediction $\eta=0.269(mT)^{3/2}$ .Determinations of the shear viscosity of trapped ultracold gases suffer from systematic, uncontrolled uncertainties related to the treatment of the dilute part of the gas cloud. In this work we present an analysis of expansion experiments based on a new method, anisotropic fluid dynamics, that interpolates between Navier-Stokes fluid dynamics at the center of the cloud and ballistic behavior in the dilute corona. We validate the method using a comparison between anisotropic fluid dynamics and numerical solutions of the Boltzmann equation. We then apply anisotropic fluid dynamics to the expansion data reported by Cao et al. In the high temperature limit we find $\eta=0.282(mT)^{3/2}$, which agrees within about 5\% with the theoretical prediction $\eta=0.269(mT)^{3/2}$.

Dec 2, 20150 citations

Physics of quantum fluids

亮五久保, 文彦高野

Jan 1, 197111 citations

Model-Independent Determination of the Shear Viscosity of a Trapped Unitary Fermi gas: Application to High-Temperature Data

Marcus Bluhm, Thomas Schäfer

Determinations of the shear viscosity of trapped ultracold gases suffer from systematic, uncontrolled uncertainties related to the treatment of the dilute part of the gas cloud. In this work we present an analysis of expansion experiments based on a new method, anisotropic fluid dynamics, that interpolates between Navier-Stokes fluid dynamics at the center of the cloud and ballistic behavior in the dilute corona. We validate the method using a comparison between anisotropic fluid dynamics and numerical solutions of the Boltzmann equation. We then apply anisotropic fluid dynamics to the expansion data reported by Cao et al. In the high-temperature limit we find η=0.282(mT)^{3/2}, which agrees within about 5% with the theoretical prediction η=0.269(mT)^{3/2}.Determinations of the shear viscosity of trapped ultracold gases suffer from systematic, uncontrolled uncertainties related to the treatment of the dilute part of the gas cloud. In this work we present an analysis of expansion experiments based on a new method, anisotropic fluid dynamics, that interpolates between Navier-Stokes fluid dynamics at the center of the cloud and ballistic behavior in the dilute corona. We validate the method using a comparison between anisotropic fluid dynamics and numerical solutions of the Boltzmann equation. We then apply anisotropic fluid dynamics to the expansion data reported by Cao et al. In the high-temperature limit we find η=0.282(mT)^{3/2}, which agrees within about 5% with the theoretical prediction η=0.269(mT)^{3/2}.

arXiv:1512.00862Physical Review LettersMar 18, 201638 citations

\eta=0.282(mT)^{3/2}

, which agrees within about 5\% with the theoretical prediction

\eta=0.269(mT)^{3/2}