Sinan Wang,Yitong Deng,Molin Deng,Hong-Xing Yu,Junwei Zhou,Duowen Chen,Taku Komura,Jiajun Wu,Bo Zhu
Abstract:We present an Eulerian vortex method based on the theory of flow maps to simulate the complex vortical motions of incompressible fluids. Central to our method is the novel incorporation of the flow-map transport equations for line elements, which, in combination with a bi-directional marching scheme for flow maps, enables the high-fidelity Eulerian advection of vorticity variables. The fundamental motivation is that, compared to impulse $\mathbf{m}$, which has been recently bridged with flow maps to encouraging results, vorticity $\boldsymbol{\omega}$ promises to be preferable for its numerical stability and physical interpretability. To realize the full potential of this novel formulation, we develop a new Poisson solving scheme for vorticity-to-velocity reconstruction that is both efficient and able to accurately handle the coupling near solid boundaries. We demonstrate the efficacy of our approach with a range of vortex simulation examples, including leapfrog vortices, vortex collisions, cavity flow, and the formation of complex vortical structures due to solid-fluid interactions.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to develop a new Eulerian Vortex Method to simulate the complex vortex motion in incompressible fluids more efficiently and accurately. Specifically, the author proposes a new vorticity - transport method based on flow maps, which achieves high - fidelity Eulerian advection through the line - element transport equations on long - range flow maps. In addition, the paper also proposes a new Poisson - solving scheme for reconstructing the velocity field from the evolved vorticity, thus avoiding the complexity brought by using the stream function in traditional methods and being able to handle solid - boundary conditions more accurately.
### Main contributions:
1. **Vorticity - transport based on flow maps**:
- The author introduces the line - element transport equations on flow maps and combines a two - way marching scheme to achieve high - fidelity Eulerian advection of vorticity variables.
- Compared with traditional impulse variables, vorticity has obvious advantages in numerical stability and physical interpretability.
2. **New Poisson - solving scheme**:
- A Poisson - solving scheme based on the velocity - vorticity form is proposed, which avoids using the stream function and potential function and directly handles solid - boundary conditions.
- This method uses an efficient matrix - free GPU solver to solve the system equations, ensuring physical fidelity and computational efficiency.
### Problems solved:
- **Numerical stability**: Traditional Eulerian vorticity methods have numerical instabilities when dealing with vorticity - stretching terms and maintaining vorticity structures. The line - element transport method proposed in this paper effectively solves this problem.
- **Physical interpretability**: Vorticity is directly related to the rotational motion of fluids and has a clear physical meaning, while impulse variables are prone to increase rapidly in numerical simulations, leading to instability.
- **Solid - boundary conditions**: Traditional vorticity - stream - function methods need to solve the velocity potential additionally to handle solid - boundary conditions. The method in this paper directly handles these conditions in the Poisson equation, simplifying the solving process.
### Experimental verification:
- The author verifies the effectiveness and robustness of this method through a series of complex vortex - simulation experiments, including vortex jumps, vortex collisions, cavity flows, and complex vortex structures formed by solid - fluid interactions.
In conclusion, by introducing the line - element transport method based on flow maps and the new Poisson - solving scheme, this paper successfully solves the deficiencies of traditional Eulerian vorticity methods in numerical stability and physical interpretability, providing a more efficient and accurate simulation method for complex vortex motions in incompressible fluids.