Two new articles by our PhD Students

Two of our students, Eduardo Andres and Luca Chiabò, published recently two journal papers.

Informed scenario-based RRT for aircraft trajectory planning under ensemble forecasting of thunderstorms

Andres, E., González-Arribas, D., Soler, M., Kamgarpour, M., Sanjurjo-Rivo, M., (2021). Transportation Research Part C: Emerging Technologies, 129, 103232.

Thunderstorms represent a major hazard for flights, as they compromise the safety of both the airframe and the passengers. To address trajectory planning under thunderstorms, three variants of the scenario-based rapidly exploring random trees (SB-RRTs) are proposed. During an iterative process, the so-called SB-RRT, the SB-RRT and the Informed SB-RRT find safe trajectories by meeting a user-defined safety threshold. Additionally, the last two techniques converge to solutions of minimum flight length. Through parallelization on graphical processing units the required computational times are reduced substantially to become compatible with near real-time operation. The proposed methods are tested considering a kinematic model of an aircraft flying between two waypoints at constant flight level and airspeed; the test scenario is based on a realistic weather forecast and assumed to be described by an ensemble of equally likely members. Lastly, the influence of the number of scenarios, safety margin and iterations on the results is analyzed. Results show that the SB-RRTs are able to find safe and, in two of the algorithms, close-to-optimum solutions.

Limitations of stationary Vlasov-Poisson solvers in probe theory

Chiabò, L. and Sánchez-Arriaga, G. (2021). Journal of Computational Physics, 438, 110366.

Physical and numerical limitations of stationary Vlasov-Poisson solvers based on backward Liouville methods are investigated with five solvers that combine different meshes, numerical integrators, and electric field interpolation schemes. Since some of the limitations arise when moving from an integrable to a non-integrable configuration, an elliptical Langmuir probe immersed in a Maxwellian plasma was considered and the eccentricity (ep) of its cross-section used as integrability-breaking parameter. In the cylindrical case, ep=0, the energy and angular momentum are both conserved. The trajectories of the charged particles are regular and the boundaries that separate trapped from non-trapped particles in phase space are smooth curves. However, their computation has to be done carefully because, albeit small, the intrinsic numerical errors of some solvers break these conservation laws. It is shown that an optimum exists for the number of loops around the probe that the solvers need to classify a particle trajectory as trapped. For ep≠0, the angular momentum is not conserved and particle dynamics in phase space is a mix of regular and chaotic orbits. The distribution function is filamented and the boundaries that separate trapped from non-trapped particles in phase space have a fractal geometry. The results were used to make a list of recommendations for the practical implementation of stationary Vlasov-Poisson solvers in a wide range of physical scenarios.

Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s