In this work we develop a new class of high order accurate Arbitrary-Lagrangian–Eulerian (ALE) one-step finite volume schemes for the solution of nonlinear systems of conservative and non-conservative hyperbolic partial differential equations. The numerical algorithm is designed for two and three space dimensions, considering moving unstructured triangular and tetrahedral meshes, respectively. As usual for finite volume schemes, data are represented within each control volume by piecewise constant values that evolve in time, hence implying the use of some strategies to improve the order of accuracy of the algorithm. In our approach high order of accuracy in space is obtained by adopting a WENO reconstruction technique, which produces piecewise polynomials of higher degree starting from the known cell averages. Such spatial high order accurate reconstruction is then employed to achieve high order of accuracy also in time using an element-local space–time finite element predictor, which performs a one-step time discretization. Specifically, we adopt a discontinuous Galerkin predictor which can handle stiff source terms that might produce jumps in the local space–time solution. Since we are dealing with moving meshes the elements deform while the solution is evolving in time, hence making the use of a reference system very convenient. Therefore, within the space–time predictor, the physical element is mapped onto a reference element using a high order isoparametric approach, where the space–time basis and test functions are given by the Lagrange interpolation polynomials passing through a predefined set of space–time nodes. The computational mesh continuously changes its configuration in time, following as closely as possible the flow motion. The entire mesh motion procedure is composed by three main steps, namely the Lagrangian step, the rezoning step and the relaxation step. In order to obtain a continuous mesh configuration at any time level, the mesh motion is evaluated by assigning each node of the computational mesh with a unique velocity vector at each timestep. The nodal solver algorithm preforms the Lagrangian stage, while we rely on a rezoning algorithm to improve the mesh quality when the flow motion becomes very complex, hence producing highly deformed computational elements. A so-called relaxation algorithm is finally employed to partially recover the optimal Lagrangian accuracy where the computational elements are not distorted too much. We underline that our scheme is supposed to be an ALE algorithm, where the local mesh velocity can be chosen independently from the local fluid velocity. Once the vertex velocity and thus the new node location has been determined, the old element configuration at time $$t^n$$ is connected with the new one at time $$t^{n+1}$$ with straight edges to represent the local mesh motion, in order to maintain algorithmic simplicity. The final ALE finite volume scheme is based directly on a space–time conservation formulation of the governing system of hyperbolic balance laws. The nonlinear system is reformulated more compactly using a space–time divergence operator and is then integrated on a moving space–time control volume. We adopt a linear parametrization of the space–time element boundaries and Gaussian quadrature rules of suitable order of accuracy to compute the integrals. We apply the new high order direct ALE finite volume schemes to several hyperbolic systems, namely the multidimensional Euler equations of compressible gas dynamics, the ideal classical magneto-hydrodynamics equations and the non-conservative seven-equation Baer–Nunziato model of compressible multi-phase flows with stiff relaxation source terms. Numerical convergence studies as well as several classical test problems will be shown to assess the accuracy and the robustness of our schemes. Finally we briefly present some variants of the algorithm that aim at improving the overall computational efficiency.