Abstract. LISFLOOD-FP 8.0 includes second-order discontinuous Galerkin (DG2) and first-order finite-volume (FV1) solvers of the two-dimensional shallow-water equations for modelling a wide range of flows, including rapidly propagating, supercritical flows, shock waves or flows over very smooth surfaces.
The solvers are parallelised on multi-core CPU and Nvidia GPU architectures and run existing LISFLOOD-FP modelling scenarios without modification.
These new, fully two-dimensional solvers are available alongside the existing local inertia solver (called ACC), which is optimised for multi-core CPUs and integrates with the LISFLOOD-FP sub-grid channel model.
The predictive capabilities and computational scalability of the new DG2 and FV1 solvers are studied for two Environment Agency benchmark tests and a real-world fluvial flood simulation driven by rainfall across a 2500 km2 catchment.
DG2's second-order-accurate, piecewise-planar representation of topography and flow variables enables predictions on coarse grids that are competitive with FV1 and ACC predictions on 2–4 times finer grids, particularly where river channels are wider than half the grid spacing.
Despite the simplified formulation of the local inertia solver, ACC is shown to be spatially second-order-accurate and yields predictions that are close to DG2. The DG2-CPU and FV1-CPU solvers achieve near-optimal scalability up to 16 CPU cores and achieve greater efficiency on grids with fewer than 0.1 million elements. The DG2-GPU and FV1-GPU solvers are most efficient on grids with more than 1 million elements, where the GPU solvers are 2.5–4 times faster than the corresponding 16-core CPU solvers. LISFLOOD-FP 8.0 therefore marks a new step towards operational DG2 flood inundation modelling at the catchment scale.
LISFLOOD-FP 8.0 is freely available under the GPL v3 license, with additional documentation and case studies at https://www.seamlesswave.com/LISFLOOD8.0 (last access: 2 June 2021).
We study the cross-field electron transport in Hall thruster induced by ExB cyclotron drift instability. The investigation tool, consisting of one-dimensional Particle-in-Cell model in the azimuthal drift direction, has been subjected to a convergence study to verify the effects of numerical parameters. The instability evolves keeping the discrete nature of its cyclotron harmonics only for low wavenumbers. A resonance broadening mechanism induced by the electron heating and velocity distribution deformation makes the high-wave numbers disappear in the long non-linear stage. A large wavelength modulation, comparable to the entire azimuthal domain considered is always superimposed. The saturation mechanism is conducted by ions that, due to friction with electrons and trapping in the potential well, heat up and rotate along the electron drift direction. With the best numerical parameters found, the scaling of anomalous mobility with the most important physical quantities (gas and plasma density, magnetic field, accelerating axial electric field, and ion mass) has been obtained. Results show that the electron cross-field mobility has a strong dependence from the plasma density and ion mass: for large plasma density, the system undergoes an abrupt change entering in a mode dominated by fluctuation-induced transport, while lighter ions present larger mobility. The scaling of the dominant wavelength detected is compatible with the first harmonic and no transition toward ion acoustic like instability has been observed.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.