Abstract
A highly energy-conservative second-order-accurate finite difference method for the cylindrical coordinate system is developed. It is rigorously proved that energy conservation in discretized space is satisfied when appropriate interpolation schemes are used. This argument holds not only for an unequally spaced mesh but also for an equally spaced mesh on cylindrical coordinates but not on Cartesian coordinates. Numerical tests are undertaken for an inviscid flow with various schemes, and it turns out that the proposed scheme offers a superior energy-conservation property and greater stability than the intuitive and previously proposed methods, for both equally spaced and unequally spaced meshes.
| Original language | English |
|---|---|
| Pages (from-to) | 478-498 |
| Number of pages | 21 |
| Journal | Journal of Computational Physics |
| Volume | 181 |
| Issue number | 2 |
| DOIs | |
| Publication status | Published - 2002 Sept 20 |
| Externally published | Yes |
Keywords
- Advection
- Body force
- Cylindrical coordinate system
- Energy conservation
- Finite difference method
- Incompressible flow
- Pipe flow
- Singularity
ASJC Scopus subject areas
- Numerical Analysis
- Modelling and Simulation
- Physics and Astronomy (miscellaneous)
- Physics and Astronomy(all)
- Computer Science Applications
- Computational Mathematics
- Applied Mathematics
