Coupling Navier-Stokes and Gross-Pitaevskii equations for the numerical simulation of two-fluid quantum flows
arXiv:2211.07361 · doi:10.1016/j.jcp.2023.112193
Abstract
Numerical methods for solving the Navier-Stokes equations for classical (or normal) viscous fluids are well established. This is also the case for the Gross-Pitaevskii equation, governing quantum inviscid flows (or superfluids) in the zero temperature limit. In quantum flows, like liquid helium II at intermediate temperatures between zero and 2.17 K, a normal fluid and a superfluid coexist with independent velocity fields. The most advanced existing models for such systems use the Navier-Stokes equations for the normal fluid and a simplified description of the superfluid, based on the dynamics of quantized vortex filaments, with ad hoc reconnection rules. There was a single attempt (C. Coste, The European Physical Journal B - Condensed Matter and Complex Systems, 1998) to couple Navier-Stokes and Gross-Pitaevskii equations in a global model intended to describe the compressible two-fluid liquid helium II. We present in this contribution a new numerical model to couple a Navier-Stokes incompressible fluid with a Gross-Pitaevskii superfluid. Coupling terms in the global system of equations involve new definitions of the following concepts: the regularized superfluid vorticity and velocity fields, the friction force exerted by quantized vortices to the normal fluid, the covariant gradient operator in the Gross-Pitaevskii based on a slip velocity respecting the dynamics of vortex lines in the normal fluid. A numerical algorithm based on pseudo-spectral Fourier methods is presented for solving the coupled system of equations.Finally, we numerically test and validate the new numerical system against well-known benchmarks for the evolution in a normal fluid of different types or arrangements of quantized vortices (vortex crystal, vortex dipole and vortex rings). The new coupling model opens new possibilities to revisit and enrich existing numerical results for complex quantum fluids.
arXiv comment32 pages, 7 figures
References in corpus (15)
- Mathematical theory and numerical methods for Bose-Einstein condensation
- The Magnus force in superfluids and superconductors
- Steady State of Counterflow Quantum Turbulence: Vortex filament Simulation with the Full Biot-Savart Law
- Onsager-Kraichnan Condensation in Decaying Two-Dimensional Quantum Turbulence
- Numerical Studies of Quantum Turbulence
- A superfluid boundary layer
- The sensitivity of the vortex filament method to different reconnection models
- Crossover from interaction to driven regimes in quantum vortex reconnections
- Dual cascade and dissipation mechanisms in helical quantum turbulence
- Derivation of the Biot-Savart equation from the Nonlinear Schrödinger equation
- Quantum turbulence simulations using the Gross-Pitaevskii equation: high-performance computing and new numerical benchmarks
- Friction-enhanced lifetime of bundled quantum vortices
- Critical velocity for vortex nucleation and roton emission in a generalized model for superfluids
- Kolmogorov and Kelvin wave cascades in a generalized model for quantum turbulence
- Coupled dynamics of quantized vortices and normal fluid in superfluid He based on lattice Boltzmann method