Showing 1–50 of 145 results for author: Ferrer, E

Reinforcement learning for anisotropic p-adaptation and error estimation in high-order solvers

Authors: David Huergo, Martín de Frutos, Eduardo Jané, Oscar A. Marino, Gonzalo Rubio, Esteban Ferrer

Abstract: We present a novel approach to automate and optimize anisotropic p-adaptation in high-order h/p solvers using Reinforcement Learning (RL). The dynamic RL adaptation uses the evolving solution to adjust the high-order polynomials. We develop an offline training approach, decoupled from the main solver, which shows minimal overcost when performing simulations. In addition, we derive a RL-based error… ▽ More We present a novel approach to automate and optimize anisotropic p-adaptation in high-order h/p solvers using Reinforcement Learning (RL). The dynamic RL adaptation uses the evolving solution to adjust the high-order polynomials. We develop an offline training approach, decoupled from the main solver, which shows minimal overcost when performing simulations. In addition, we derive a RL-based error estimation approach that enables the quantification of local discretization errors. The proposed methodology is agnostic to both the computational mesh and the partial differential equation being solved. The application of RL to mesh adaptation offers several benefits. It enables automated, adaptive mesh refinement, reducing the need for manual intervention. It optimizes computational resources by dynamically allocating high-order polynomials where necessary and minimizing refinement in stable regions. This leads to computational cost savings while maintaining solution accuracy. Furthermore, RL allows for the exploration of unconventional mesh adaptations, potentially enhancing the accuracy and robustness of simulations. This work extends our original research, offering a more robust, reproducible, and generalizable approach applicable to complex three-dimensional problems. We provide validation for laminar and turbulent cases: circular cylinders, Taylor Green Vortex and a 10MW wind turbine to illustrate the flexibility of the proposed approach. △ Less

Submitted 26 July, 2024; originally announced July 2024.

Comments: 36 pages, 16 figures, 5 tables

A reinforcement learning strategy to automate and accelerate h/p-multigrid solvers

Authors: David Huergo, Laura Alonso, Saumitra Joshi, Adrian Juanicoteca, Gonzalo Rubio, Esteban Ferrer

Abstract: We explore a reinforcement learning strategy to automate and accelerate h/p-multigrid methods in high-order solvers. Multigrid methods are very efficient but require fine-tuning of numerical parameters, such as the number of smoothing sweeps per level and the correction fraction (i.e., proportion of the corrected solution that is transferred from a coarser grid to a finer grid). The objective of t… ▽ More We explore a reinforcement learning strategy to automate and accelerate h/p-multigrid methods in high-order solvers. Multigrid methods are very efficient but require fine-tuning of numerical parameters, such as the number of smoothing sweeps per level and the correction fraction (i.e., proportion of the corrected solution that is transferred from a coarser grid to a finer grid). The objective of this paper is to use a proximal policy optimization algorithm to automatically tune the multigrid parameters and, by doing so, improve stability and efficiency of the h/p-multigrid strategy. Our findings reveal that the proposed reinforcement learning h/p-multigrid approach significantly accelerates and improves the robustness of steady-state simulations for one dimensional advection-diffusion and nonlinear Burgers' equations, when discretized using high-order h/p methods, on uniform and nonuniform grids. △ Less

Submitted 18 July, 2024; originally announced July 2024.

Comments: 14 pages, 1 figure

Deep Reinforcement Learning for Multi-Objective Optimization: Enhancing Wind Turbine Energy Generation while Mitigating Noise Emissions

Authors: Martín de Frutos, Oscar A. Marino, David Huergo, Esteban Ferrer

Abstract: We develop a torque-pitch control framework using deep reinforcement learning for wind turbines to optimize the generation of wind turbine energy while minimizing operational noise. We employ a double deep Q-learning, coupled to a blade element momentum solver, to enable precise control over wind turbine parameters. In addition to the blade element momentum, we use the wind turbine acoustic model… ▽ More We develop a torque-pitch control framework using deep reinforcement learning for wind turbines to optimize the generation of wind turbine energy while minimizing operational noise. We employ a double deep Q-learning, coupled to a blade element momentum solver, to enable precise control over wind turbine parameters. In addition to the blade element momentum, we use the wind turbine acoustic model of Brooks Pope and Marcolini. Through training with simple winds, the agent learns optimal control policies that allow efficient control for complex turbulent winds. Our experiments demonstrate that the reinforcement learning is able to find optima at the Pareto front, when maximizing energy while minimizing noise. In addition, the adaptability of the reinforcement learning agent to changing turbulent wind conditions, underscores its efficacy for real-world applications. We validate the methodology using a SWT2.3-93 wind turbine with a rated power of 2.3 MW. We compare the reinforcement learning control to classic controls to show that they are comparable when not taking into account noise emissions. When including a maximum limit of 45 dB to the noise produced (100 meters downwind of the turbine), the extracted yearly energy decreases by 22%. The methodology is flexible and allows for easy tuning of the objectives and constraints through the reward definitions, resulting in a flexible multi-objective optimization framework for wind turbine control. Overall, our findings highlight the potential of RL-based control strategies to improve wind turbine efficiency while mitigating noise pollution, thus advancing sustainable energy generation technologies △ Less

Submitted 18 July, 2024; originally announced July 2024.

Modelling Wind Turbines via Actuator Lines in High-Order h/p Solvers

Authors: Oscar A. Marino, Raúl Sanz, Stefano Colombo, Ananth Sivaramakrishnan, Esteban Ferrer

Abstract: This paper compares two actuator line methodologies for modelling wind turbines employing high-order h/p solvers and large-eddy simulations. The methods combine the accuracy of high-order solvers (in this work the maximum order is 6) with the computational efficiency of actuator lines to capture the aerodynamic effects of wind turbine blades. Comparisons with experiments validate the actuator li… ▽ More This paper compares two actuator line methodologies for modelling wind turbines employing high-order h/p solvers and large-eddy simulations. The methods combine the accuracy of high-order solvers (in this work the maximum order is 6) with the computational efficiency of actuator lines to capture the aerodynamic effects of wind turbine blades. Comparisons with experiments validate the actuator line methodologies. We explore the effects of the polynomial order and the smoothing parameter associated with the Gaussian regularization function, and use them to blend the actuator line forcing in the high-order computational mesh, to show that both parameters influence the distribution of forces along the blades and the turbine wake. The greatest impact is obtained when the polynomial order is increased, allowing one to better capture the physics without requiring new meshes. When comparing the actuator line methodologies, we show the advantages of performing weighted sums, over element averages, to compute the blade velocities and forces in high-order solvers. For low-order simulations (low polynomial orders), both methods provide similar results, but as the polynomial order is increased, the weighted sum method shows smoother thrust/torque distributions and better wake resolution. Furthermore, cell averaging introduces nonphysical oscillations when increasing the polynomial order beyond 3 (4th order accuracy). △ Less

Submitted 14 June, 2024; originally announced June 2024.

Low-cost wind turbine aeroacoustic predictions using actuator lines

Authors: Laura Botero-Bolivar, Oscar A Marino, Cornelis H. Venner, Leandro D. de Santana, Esteban Ferrer

Abstract: Aerodynamic noise is a limitation for further exploitation of wind energy resources. As this type of noise is caused by the interaction of turbulent flow with the airframe, a detailed resolution of the flow is necessary to obtain an accurate prediction of the far-field noise. Computational fluid dynamic (CFD) solvers simulate the flow field but only at a high computational cost, which is much incr… ▽ More Aerodynamic noise is a limitation for further exploitation of wind energy resources. As this type of noise is caused by the interaction of turbulent flow with the airframe, a detailed resolution of the flow is necessary to obtain an accurate prediction of the far-field noise. Computational fluid dynamic (CFD) solvers simulate the flow field but only at a high computational cost, which is much increased when the acoustic field is resolved. Therefore, wind turbine noise predictions using numerical approaches remain a challenge. This paper presents a methodology that couples (relatively fast) wind turbine CFD simulations using actuator lines with a fast turn-around noise prediction method. The flow field is simulated using actuator lines and large eddy simulations. The noise prediction method is based on the Amiet-Schlinker's theory for rotatory noise sources, considering leading- and trailing-edge noise as unique noise sources. A 2D panel code (XFOIL) calculates the sectional forces and boundary layer quantities. The resulting methodology for the noise prediction method is of high fidelity since the wind turbine geometry is accounted for in both flow and acoustics predictions. Results are compared with field measurements of a full-scale wind turbine for two operational conditions, validating the results of this research. △ Less

Submitted 8 June, 2024; originally announced June 2024.

Comments: 13 pages, 12 figures

Journal ref: Renewable Energy Volume 227, June 2024, 120476

Implementation of Immersed Boundaries via Volume Penalization in the Industrial Aeronautical Computational Fluid Dynamics Solver CODA

Authors: Jonatan Nunez, David Huergo, Diego Lodares, Suyash Shrestha, Juan Guerra, Juan Florenciano, Esteban Ferrer, Eusebio Valero

Abstract: We present the implementation and validation of an immersed boundary volume penalization method in the computational fluid dynamics solver CODA (from ONERA, DLR, and Airbus). Our goal is to model and simulate turbulent fluid flows in complex 3D aerodynamic configurations through the numerical solution of the Reynolds--averaged Navier--Stokes equations using the Spalart--Allmaras turbulent model. T… ▽ More We present the implementation and validation of an immersed boundary volume penalization method in the computational fluid dynamics solver CODA (from ONERA, DLR, and Airbus). Our goal is to model and simulate turbulent fluid flows in complex 3D aerodynamic configurations through the numerical solution of the Reynolds--averaged Navier--Stokes equations using the Spalart--Allmaras turbulent model. To do that, an immersed boundary method has been implemented in CODA and an efficient preprocessing tool for the construction of unstructured hexahedral meshes with adaptive mesh refinement around immersed geometries has been developed. We report several numerical examples, including subsonic and transonic flow past a NACA0012 airfoil, subsonic flow past an MDA30P30N multi-element airfoil, and subsonic flow around the NASA high-lift CRM aircraft. These simulations have been performed in the CODA solver with a second-order finite volume scheme as spatial discretization and an implicit backward Euler scheme based on the matrix-free GMRES block-Jacobi iterative method. The reported numerical simulations are in very good agreement with their corresponding experimental data. We conclude that the implemented immersed boundary method is efficient, flexible, and accurate and can therefore be used for aeronautical applications in industry. △ Less

Submitted 20 July, 2024; v1 submitted 24 April, 2024; originally announced April 2024.

Comments: 22 pages, 15 figures

Axion-Polaritons in quark stars: a possible solution to the missing pulsar problem

Authors: E. J. Ferrer, V. de la Incera

Abstract: This paper proposes an alternative mechanism to solve the so-called missing pulsar problem, a standing paradox between the theoretical expectations about the number of pulsars that should exist in the galaxy center of the Milky Way and their absence in the observations. The mechanism is based on the transformation of incident $γ$ rays into hybridized modes, known as axion-polaritons, which can exi… ▽ More This paper proposes an alternative mechanism to solve the so-called missing pulsar problem, a standing paradox between the theoretical expectations about the number of pulsars that should exist in the galaxy center of the Milky Way and their absence in the observations. The mechanism is based on the transformation of incident $γ$ rays into hybridized modes, known as axion-polaritons, which can exist inside highly magnetized quark stars with a quark matter phase known as the magnetic dual chiral density wave phase. This phase, which is favored over several other dense matter phases candidates at densities a few times nuclear saturation density, has already passed several important astrophysical tests. In the proposed mechanism, the absence of young magnetars occurs because as electromagnetic waves inside the star can only propagate through the hybridized modes, incident photons coming from a $γ$-ray burst get transformed into massless and massive axion polaritons by the Primakoff effect. Once thermalized, the massive axion-polaritons can self-gravitate up to a situation where their total mass overpasses the Chandrasekhar limit for these bosons, producing a mini blackhole that collapses the star. △ Less

Submitted 24 March, 2024; originally announced March 2024.

Journal ref: Eur.Phys.J.C 84 (2024) 133

Accelerating high order discontinuous Galerkin solvers through a clustering-based viscous/turbulent-inviscid domain decomposition

Authors: Kheir-Eddine Otmani, Andrés Mateo-Gabín, Gonzalo Rubio, Esteban Ferrer

Abstract: We explore the unsupervised clustering technique introduced in [25] to identify viscous/turbulent from inviscid regions in incompressible flows. The separation of regions allows solving the Navier-Stokes equations including Large Eddy Simulation closure models only in the viscous/turbulent ones, while solving the Euler equations in the remaining of the computational domain. By solving different se… ▽ More We explore the unsupervised clustering technique introduced in [25] to identify viscous/turbulent from inviscid regions in incompressible flows. The separation of regions allows solving the Navier-Stokes equations including Large Eddy Simulation closure models only in the viscous/turbulent ones, while solving the Euler equations in the remaining of the computational domain. By solving different sets of equations, the computational cost is significantly reduced. This coupling strategy is implemented within a discontinuous Galerkin numerical framework, which allows discontinuous solutions (i.e., different set of equations) in neighbouring elements that interact through numerical fluxes. The proposed strategy maintains the same accuracy at lower cost, when compared to solving the full Navier-Stokes equations throughout the entire domain. Validation of this approach is conducted across diverse flow regimes, spanning from unsteady laminar flows to unsteady turbulent flows, including an airfoil section at Reynolds numbers Re = 103 and 104 and large angles of attack, and the flow past a wind turbine, modelled using actuator lines. The computational cost is reduced by 25% and 29% for the unsteady turbulent flow around an airfoil section and the flow past the wind turbine, respectively. In addition, to further accelerate the simulations, we combine the proposed decoupling with local P -adaptation, as proposed in [ 30]. When doing so, we reduce the computational cost by 41% and 45% for the flow around the airfoil section and the flow past the wind turbine, respectively △ Less

Submitted 8 March, 2024; originally announced March 2024.

Reinforcement learning to maximise wind turbine energy generation

Authors: Daniel Soler, Oscar Mariño, David Huergo, Martín de Frutos, Esteban Ferrer

Abstract: We propose a reinforcement learning strategy to control wind turbine energy generation by actively changing the rotor speed, the rotor yaw angle and the blade pitch angle. A double deep Q-learning with a prioritized experience replay agent is coupled with a blade element momentum model and is trained to allow control for changing winds. The agent is trained to decide the best control (speed, yaw,… ▽ More We propose a reinforcement learning strategy to control wind turbine energy generation by actively changing the rotor speed, the rotor yaw angle and the blade pitch angle. A double deep Q-learning with a prioritized experience replay agent is coupled with a blade element momentum model and is trained to allow control for changing winds. The agent is trained to decide the best control (speed, yaw, pitch) for simple steady winds and is subsequently challenged with real dynamic turbulent winds, showing good performance. The double deep Q- learning is compared with a classic value iteration reinforcement learning control and both strategies outperform a classic PID control in all environments. Furthermore, the reinforcement learning approach is well suited to changing environments including turbulent/gusty winds, showing great adaptability. Finally, we compare all control strategies with real winds and compute the annual energy production. In this case, the double deep Q-learning algorithm also outperforms classic methodologies. △ Less

Submitted 17 February, 2024; originally announced February 2024.

A comparative study of explicit and implicit Large Eddy Simulations using a high-order discontinuous Galerkin solver: application to a Formula 1 front wing

Authors: Gerasimos Ntoukas, Gonzalo Rubio, Oscar Marino, Alexandra Liosi, Francesco Bottone, Julien Hoessler, Esteban Ferrer

Abstract: This paper explores two Large Eddy Simulation (LES) approaches within the framework of the high-order discontinuous Galerkin solver, Horses3D. The investigation focuses on an Inverted Multi-element Wing in Ground Effect (i.e. 2.5D Imperial Front Wing section) representing a Formula 1 front wing, and compares the strengths and limitations of the two LES methods. The explicit LES formulation relies… ▽ More This paper explores two Large Eddy Simulation (LES) approaches within the framework of the high-order discontinuous Galerkin solver, Horses3D. The investigation focuses on an Inverted Multi-element Wing in Ground Effect (i.e. 2.5D Imperial Front Wing section) representing a Formula 1 front wing, and compares the strengths and limitations of the two LES methods. The explicit LES formulation relies on the Vreman model, that adapts to laminar, transitional and turbulent regimes. The numerical formulation uses nodal basis functions and Gauss points. The implicit LES formulation, does not require explicit turbulence modeling but relies in the discretization scheme. We use the Kennedy-Gruber entropy stable formulation to enhance stability in under resolved simulations, since we recover the continuous properties such as entropy conservation at a discrete level. This formulation employs Gauss-Lobatto points, which downgrades the accuracy of integration but allows for larger time steps in explicit time integration. We compare our results to Nektar++ [1] showing that both LES techniques provide results that agree well with the reference values. The implicit LES shows to better capture transition and allows for larger time steps at a similar cost per iteration. We conclude that this implicit LES formulation is very attractive for complex simulations. △ Less

Submitted 15 January, 2024; originally announced February 2024.

Making Mathematical Research Data FAIR: A Technology Overview

Authors: Tim Conrad, Eloi Ferrer, Daniel Mietchen, Larissa Pusch, Johannes Stegmuller, Moritz Schubotz

Abstract: The sharing and citation of research data is becoming increasingly recognized as an essential building block in scientific research across various fields and disciplines. Sharing research data allows other researchers to reproduce results, replicate findings, and build on them. Ultimately, this will foster faster cycles in knowledge generation. Some disciplines, such as astronomy or bioinformatics… ▽ More The sharing and citation of research data is becoming increasingly recognized as an essential building block in scientific research across various fields and disciplines. Sharing research data allows other researchers to reproduce results, replicate findings, and build on them. Ultimately, this will foster faster cycles in knowledge generation. Some disciplines, such as astronomy or bioinformatics, already have a long history of sharing data; many others do not. The current landscape of so-called research data repositories is diverse. This review aims to perform a technology review on existing data repositories/portals with a focus on mathematical research data. △ Less

Submitted 21 September, 2023; originally announced September 2023.

Bravo MaRDI: A Wikibase Powered Knowledge Graph on Mathematics

Authors: Moritz Schubotz, Eloi Ferrer, Johannes Stegmüller, Daniel Mietchen, Olaf Teschke, Larissa Pusch, Tim OF Conrad

Abstract: Mathematical world knowledge is a fundamental component of Wikidata. However, to date, no expertly curated knowledge graph has focused specifically on contemporary mathematics. Addressing this gap, the Mathematical Research Data Initiative (MaRDI) has developed a comprehensive knowledge graph that links multimodal research data in mathematics. This encompasses traditional research data items like… ▽ More Mathematical world knowledge is a fundamental component of Wikidata. However, to date, no expertly curated knowledge graph has focused specifically on contemporary mathematics. Addressing this gap, the Mathematical Research Data Initiative (MaRDI) has developed a comprehensive knowledge graph that links multimodal research data in mathematics. This encompasses traditional research data items like datasets, software, and publications and includes semantically advanced objects such as mathematical formulas and hypotheses. This paper details the abilities of the MaRDI knowledge graph, which is based on Wikibase, leading up to its inaugural public release, codenamed Bravo, available on https://portal.mardi4nfdi.de. △ Less

Submitted 20 September, 2023; originally announced September 2023.

Comments: Accepted at Wikidata'23: Wikidata workshop at ISWC 2023

On the $Σ$-invariants of Artin groups satisfying the $K(π,1)$-conjecture

Authors: Marcos Escartín Ferrer, Conchita Martínez Pérez

Abstract: We consider $Σ$-invariants of Artin groups that satisfy the $K(π,1)$-conjecture. These invariants determine the cohomological finiteness conditions of subgroups that contain the derived subgroup. We extend a known result for even Artin groups of FC-type, giving a sufficient condition for a character $χ:A_Γ\to\mathbb{R}$ to belong to $Σ^n(A_Γ,\mathbb{Z})$. We also prove some partial converses. As a… ▽ More We consider $Σ$-invariants of Artin groups that satisfy the $K(π,1)$-conjecture. These invariants determine the cohomological finiteness conditions of subgroups that contain the derived subgroup. We extend a known result for even Artin groups of FC-type, giving a sufficient condition for a character $χ:A_Γ\to\mathbb{R}$ to belong to $Σ^n(A_Γ,\mathbb{Z})$. We also prove some partial converses. As applications, we prove that the $Σ^1$-conjecture holds true when there is a prime $p$ that divides $l(e)/2$ for any edge with even label $l(e)>2$, we generalize to Artin groups the homological version of Bestvina-Brady theorem and we compute the $Σ$-invariants of all irreducible spherical and affine Artin groups and triangle Artin groups, which provide a complete classification of the $F_n$ and $FP_n$ properties of their derived subgroup. △ Less

Submitted 11 December, 2023; v1 submitted 6 September, 2023; originally announced September 2023.

Comments: 25 pages. Revised version, some missprints corrected. Accepted for publication in the JLMS

MSC Class: Primary 20J06; 20F36; Secondary 57M07; 55P20

An unsupervised machine-learning-based shock sensor for high-order supersonic flow solvers

Authors: Andrés Mateo-Gabín, Kenza Tlales, Eusebio Valero, Esteban Ferrer, Gonzalo Rubio

Abstract: We present a novel unsupervised machine-learning sock sensor based on Gaussian Mixture Models (GMMs). The proposed GMM sensor demonstrates remarkable accuracy in detecting shocks and is robust across diverse test cases with significantly less parameter tuning than other options. We compare the GMM-based sensor with state-of-the-art alternatives. All methods are integrated into a high-order compres… ▽ More We present a novel unsupervised machine-learning sock sensor based on Gaussian Mixture Models (GMMs). The proposed GMM sensor demonstrates remarkable accuracy in detecting shocks and is robust across diverse test cases with significantly less parameter tuning than other options. We compare the GMM-based sensor with state-of-the-art alternatives. All methods are integrated into a high-order compressible discontinuous Galerkin solver, where two stabilization approaches are coupled to the sensor to provide examples of possible applications. The Sedov blast and double Mach reflection cases demonstrate that our proposed sensor can enhance hybrid sub-cell flux-differencing formulations by providing accurate information of the nodes that require low-order blending. Besides, supersonic test cases including high Reynolds numbers showcase the sensor performance when used to introduce entropy-stable artificial viscosity to capture shocks, demonstrating the same effectiveness as fine-tuned state-of-the-art sensors. The adaptive nature and ability to function without extensive training datasets make this GMM-based sensor suitable for complex geometries and varied flow configurations. Our study reveals the potential of unsupervised machine-learning methods, exemplified by this GMM sensor, to improve the robustness and efficiency of advanced CFD codes. △ Less

Submitted 9 October, 2023; v1 submitted 28 July, 2023; originally announced August 2023.

Thermal phonon fluctuations and stability of the magnetic dual chiral density wave phase in dense QCD

Authors: E. J Ferrer, W. Gyory, V. de la Incera

Abstract: We study the stability against thermal phonon fluctuations of the magnetic dual chiral density wave (MDCDW) phase, an inhomogeneous phase arising in cold dense QCD in a magnetic field. Following a recent study that demonstrated the absence of the Landau-Peierls (LP) instability from this phase, we calculate the (threshold) temperature at which the phonon fluctuations wash out the long-range order… ▽ More We study the stability against thermal phonon fluctuations of the magnetic dual chiral density wave (MDCDW) phase, an inhomogeneous phase arising in cold dense QCD in a magnetic field. Following a recent study that demonstrated the absence of the Landau-Peierls (LP) instability from this phase, we calculate the (threshold) temperature at which the phonon fluctuations wash out the long-range order over a range of magnetic fields and densities relevant to astrophysical applications. Using a high-order Ginzburg-Landau expansion, we find that the threshold temperature is very near the critical temperature for fields of order $10^{18}$ G, and still a sizable fraction of the critical temperature for fields of order $10^{17}$ G. Therefore, at sufficiently large magnetic fields, the long-range order of the MDCDW phase is preserved over most of the parameter space where it is energetically favored; at smaller magnetic fields, the long-range order is still preserved over a considerable region of parameter space relevant to compact stars. We provide general symmetry arguments to explain why a magnetic field alone is not enough to eliminate the LP instability that characterizes single-modulated phases in 3+1 dimensions. △ Less

Submitted 24 March, 2024; v1 submitted 11 July, 2023; originally announced July 2023.

Comments: A new discussion was added on the roles of magnetic fields and condensate symmetry-breaking in eliminating the LP instability

A reinforcement learning strategy for p-adaptation in high order solvers

Authors: David Huergo, Gonzalo Rubio, Esteban Ferrer

Abstract: Reinforcement learning (RL) has emerged as a promising approach to automating decision processes. This paper explores the application of RL techniques to optimise the polynomial order in the computational mesh when using high-order solvers. Mesh adaptation plays a crucial role in improving the efficiency of numerical simulations by improving accuracy while reducing the cost. Here, actor-critic RL… ▽ More Reinforcement learning (RL) has emerged as a promising approach to automating decision processes. This paper explores the application of RL techniques to optimise the polynomial order in the computational mesh when using high-order solvers. Mesh adaptation plays a crucial role in improving the efficiency of numerical simulations by improving accuracy while reducing the cost. Here, actor-critic RL models based on Proximal Policy Optimization offer a data-driven approach for agents to learn optimal mesh modifications based on evolving conditions. The paper provides a strategy for p-adaptation in high-order solvers and includes insights into the main aspects of RL-based mesh adaptation, including the formulation of appropriate reward structures and the interaction between the RL agent and the simulation environment. We discuss the impact of RL-based mesh p-adaptation on computational efficiency and accuracy. We test the RL p-adaptation strategy on a 1D inviscid Burgers' equation to demonstrate the effectiveness of the strategy. The RL strategy reduces the computational cost and improves accuracy over uniform adaptation, while minimising human intervention. △ Less

Submitted 14 June, 2023; originally announced June 2023.

A combined volume penalization / selective frequency damping approach for immersed boundary methods: application to moving geometries

Authors: Jiaqing Kou, Esteban Ferrer

Abstract: This work extends, to moving geometries, the immersed boundary method based on volume penalization and selective frequency damping approach [J. Kou, E. Ferrer, A combined volume penalization/selective frequency damping approach for immersed boundary methods applied to high-order schemes, Journal of Computational Physics (2023)]. To do so, the numerical solution inside the solid is decomposed into… ▽ More This work extends, to moving geometries, the immersed boundary method based on volume penalization and selective frequency damping approach [J. Kou, E. Ferrer, A combined volume penalization/selective frequency damping approach for immersed boundary methods applied to high-order schemes, Journal of Computational Physics (2023)]. To do so, the numerical solution inside the solid is decomposed into a predefined movement and an oscillatory part (spurious waves), where the latter is damped by an SFD approach combined with volume penalization. We challenge the method with two cases. First, a new manufactured solution problem is proposed to show that the method can recover high-order accuracy. Second, we validate the methodology by simulating the laminar flow past a moving cylinder, where improved accuracy of the combined method is reported. △ Less

Submitted 1 June, 2023; originally announced June 2023.

A modified equation analysis for immersed boundary methods based on volume penalization: applications to linear advection-diffusion and high-order discontinuous Galerkin schemes

Authors: Victor J. Llorente, Jiaqing Kou, Eusebio Valero, Esteban Ferrer

Abstract: The Immersed Boundary Method (IBM) is a popular numerical approach to impose boundary conditions without relying on body-fitted grids, thus reducing the costly effort of mesh generation. To obtain enhanced accuracy, IBM can be combined with high-order methods (e.g., discontinuous Galerkin). For this combination to be effective, an analysis of the numerical errors is essential. In this work, we app… ▽ More The Immersed Boundary Method (IBM) is a popular numerical approach to impose boundary conditions without relying on body-fitted grids, thus reducing the costly effort of mesh generation. To obtain enhanced accuracy, IBM can be combined with high-order methods (e.g., discontinuous Galerkin). For this combination to be effective, an analysis of the numerical errors is essential. In this work, we apply, for the first time, a modified equation analysis to the combination of IBM (based on volume penalization) and high-order methods (based on nodal discontinuous Galerkin methods) to analyze a priori numerical errors and obtain practical guidelines on the selection of IBM parameters. The analysis is performed on a linear advection-diffusion equation with Dirichlet boundary conditions. Three ways to penalize the immerse boundary are considered, the first penalizes the solution inside the IBM region (classic approach), whilst the second and third penalize the first and second derivatives of the solution. We find optimal combinations of the penalization parameters, including the first and second penalizing derivatives, resulting in minimum errors. We validate the theoretical analysis with numerical experiments for one- and two-dimensional advection-diffusion equations. △ Less

Submitted 19 December, 2022; originally announced December 2022.

Improving aircraft performance using machine learning: a review

Authors: Soledad Le Clainche, Esteban Ferrer, Sam Gibson, Elisabeth Cross, Alessandro Parente, Ricardo Vinuesa

Abstract: This review covers the new developments in machine learning (ML) that are impacting the multi-disciplinary area of aerospace engineering, including fundamental fluid dynamics (experimental and numerical), aerodynamics, acoustics, combustion and structural health monitoring. We review the state of the art, gathering the advantages and challenges of ML methods across different aerospace disciplines… ▽ More This review covers the new developments in machine learning (ML) that are impacting the multi-disciplinary area of aerospace engineering, including fundamental fluid dynamics (experimental and numerical), aerodynamics, acoustics, combustion and structural health monitoring. We review the state of the art, gathering the advantages and challenges of ML methods across different aerospace disciplines and provide our view on future opportunities. The basic concepts and the most relevant strategies for ML are presented together with the most relevant applications in aerospace engineering, revealing that ML is improving aircraft performance and that these techniques will have a large impact in the near future. △ Less

Submitted 20 October, 2022; originally announced October 2022.

Truncation Error-Based Anisotropic $p$-Adaptation for Unsteady Flows for High-Order Discontinuous Galerkin Methods

Authors: Andrés M. Rueda-Ramírez, Gerasimos Ntoukas, Gonzalo Rubio, Eusebio Valero, Esteban Ferrer

Abstract: In this work, we extend the $τ$-estimation method to unsteady problems and use it to adapt the polynomial degree for high-order discontinuous Galerkin simulations of unsteady flows. The adaptation is local and anisotropic and allows capturing relevant unsteady flow features while enhancing the accuracy of time evolving functionals (e.g., lift, drag). To achieve an efficient and unsteady truncation… ▽ More In this work, we extend the $τ$-estimation method to unsteady problems and use it to adapt the polynomial degree for high-order discontinuous Galerkin simulations of unsteady flows. The adaptation is local and anisotropic and allows capturing relevant unsteady flow features while enhancing the accuracy of time evolving functionals (e.g., lift, drag). To achieve an efficient and unsteady truncation error-based $p$-adaptation scheme, we first revisit the definition of the truncation error, studying the effect of the treatment of the mass matrix arising from the temporal term. Secondly, we extend the $τ$-estimation strategy to unsteady problems. Finally, we present and compare two adaptation strategies for unsteady problems: the dynamic and static $p$-adaptation methods. In the first one (dynamic) the error is measured periodically during a simulation and the polynomial degree is adapted immediately after every estimation procedure. In the second one (static) the error is also measured periodically, but only one $p$-adaptation process is performed after several estimation stages, using a combination of the periodic error measures. The static $p$-adaptation strategy is suitable for time-periodic flows, while the dynamic one can be generalized to any flow evolution. We consider two test cases to evaluate the efficiency of the proposed $p$-adaptation strategies. The first one considers the compressible Euler equations to simulate the advection of a density pulse. The second one solves the compressible Navier-Stokes equations to simulate the flow around a cylinder at Re=100. The local and anisotropic adaptation enables significant reductions in the number of degrees of freedom with respect to uniform refinement, leading to speed-ups of up to $\times4.5$ for the Euler test case and $\times2.2$ for the Navier-Stokes test case. △ Less

Submitted 7 October, 2022; originally announced October 2022.

Aeroacoustic airfoil shape optimization enhanced by autoencoders

Authors: Jiaqing Kou, Laura Botero-Bolívar, Román Ballano, Oscar Marino, Leandro de Santana, Eusebio Valero, Esteban Ferrer

Abstract: We present a framework for airfoil shape optimization to reduce the trailing edge noise for the design of wind turbine blades. Far-field noise is evaluated using Amiet's theory coupled with the TNO-Blake model to calculate the wall pressure spectrum and fast turn-around XFOIL simulations to evaluate the boundary layer parameters. The computational framework is first validated using a NACA0012 airf… ▽ More We present a framework for airfoil shape optimization to reduce the trailing edge noise for the design of wind turbine blades. Far-field noise is evaluated using Amiet's theory coupled with the TNO-Blake model to calculate the wall pressure spectrum and fast turn-around XFOIL simulations to evaluate the boundary layer parameters. The computational framework is first validated using a NACA0012 airfoil at zero angle of attack. Particle swarm optimization is used to find the optimized airfoil configuration. The multi-objective optimization minimizes the A-weighted overall sound pressure level at various angles of attack, while ensuring enough lift and minimum drag. We compare classic parametrization methods to define the airfoil geometry (i.e., CST) to a machine learning method (i.e., a variational autoencoder). We observe that variational autoencoders can represent a wide variety of geometries, with only four encoded variables, leading to efficient optimizations, which result in improved optimal shapes. When compared to the baseline geometry, a NACA0012, the autoencoder-based optimized airfoil reduces by 3% (1.75 dBA) the overall sound pressure level (with decreased noise across the entire frequency range), while maintaining favorable aerodynamic properties in terms of lift and drag. △ Less

Submitted 30 September, 2022; originally announced October 2022.

Machine learning adaptation for laminar and turbulent flows: applications to high order discontinuous Galerkin solvers

Authors: Kenza Tlales, Kheir-Eddine Otmani, Gerasimos Ntoukas, Gonzalo Rubio, Esteban Ferrer

Abstract: We present a machine learning-based mesh refinement technique for steady and unsteady flows. The clustering technique proposed by Otmani et al. arXiv:2207.02929 [physics.flu-dyn] is used to mark the viscous and turbulent regions for the flow past a cylinder at Re=40 (steady laminar flow) and Re=3900 (unsteady turbulent flow). Within this clustered region, we increase the polynomial order to show t… ▽ More We present a machine learning-based mesh refinement technique for steady and unsteady flows. The clustering technique proposed by Otmani et al. arXiv:2207.02929 [physics.flu-dyn] is used to mark the viscous and turbulent regions for the flow past a cylinder at Re=40 (steady laminar flow) and Re=3900 (unsteady turbulent flow). Within this clustered region, we increase the polynomial order to show that it is possible to obtain similar levels of accuracy to a uniformly refined mesh. The method is effective as the clustering successfully identifies the two flow regions, a viscous/turbulent dominated region (including the boundary layer and wake) and an inviscid/irrotational region (a potential flow region). The data used within this framework are generated using a high-order discontinuous Galerkin solver, allowing to locally refine the polynomial order (p-refinement) in each element of the clustered region. For the steady laminar test case we are able to reduce the computational cost up to 32% and for the unsteady turbulent case up to 33%. △ Less

Submitted 6 September, 2022; originally announced September 2022.

Jump penalty stabilisation techniques for under-resolved turbulence in discontinuous Galerkin schemes

Authors: Jiaqing Kou, Oscar A. Marino, Esteban Ferrer

Abstract: Jump penalty stabilisation techniques have been recently proposed for continuous and discontinuous high order Galerkin schemes [1,2,3]. The stabilisation relies on the gradient or solution discontinuity at element interfaces to incorporate localised numerical diffusion in the numerical scheme. This diffusion acts as an implicit subgrid model and stablises under-resolved turbulent simulations. Th… ▽ More Jump penalty stabilisation techniques have been recently proposed for continuous and discontinuous high order Galerkin schemes [1,2,3]. The stabilisation relies on the gradient or solution discontinuity at element interfaces to incorporate localised numerical diffusion in the numerical scheme. This diffusion acts as an implicit subgrid model and stablises under-resolved turbulent simulations. This paper investigates the effect of jump penalty stabilisation methods (penalising gradient or solution) for stabilisation and improvement of high-order discontinuous Galerkin schemes in turbulent regime. We analyse these schemes using an eigensolution analysis, a 1D non-linear Burgers equation (mimicking a turbulent cascade) and 3D turbulent Navier-Stokes simulations (Taylor-Green Vortex problem). We show that the two jump penalty stabilisation techniques can stabilise under-resolved simulations thanks to the improved dispersion-dissipation characteristics (when compared to non-penalised schemes) and provide accurate results for turbulent flows. The numerical results indicate that the proposed jump penalty stabilise under-resolved simulations and improve the simulations, when compared to the original unpenalised scheme and to classic explicit subgrid models (Smagorisnky and Vreman). △ Less

Submitted 24 August, 2022; originally announced August 2022.

Accelerating high order discontinuous Galerkin solvers using neural networks: 3D compressible Navier-Stokes equations

Authors: Fernando Manrique de Lara, Esteban Ferrer

Abstract: We propose to accelerate a high order discontinuous Galerkin solver using neural networks. We include a corrective forcing to a low polynomial order simulation to enhance its accuracy. The forcing is obtained by training a deep fully connected neural network, using a high polynomial order simulation but only for a short time frame. With this corrective forcing, we can run the low polynomial order… ▽ More We propose to accelerate a high order discontinuous Galerkin solver using neural networks. We include a corrective forcing to a low polynomial order simulation to enhance its accuracy. The forcing is obtained by training a deep fully connected neural network, using a high polynomial order simulation but only for a short time frame. With this corrective forcing, we can run the low polynomial order simulation faster (with large time steps and low cost per time step) while improving its accuracy. We explored this idea for a 1D Burgers' equation in (Marique and Ferrer, CAF 2022), and we have extended this work to the 3D Navier-Stokes equations, with and without a Large Eddy Simulation closure model. We test the methodology with the turbulent Taylor Green Vortex case and for various Reynolds numbers (30, 200 and 1600). In addition, the Taylor Green Vortex evolves with time and covers laminar, transitional, and turbulent regimes, as time progresses. The proposed methodology proves to be applicable to a variety of flows and regimes. The results show that the corrective forcing is effective in all Reynolds numbers and time frames (excluding the initial flow development). We can train the corrective forcing with a polynomial order of 8, to increase the accuracy of simulations from a polynomial order 3 to 6, when correcting outside the training time frame. The low order correct solution is 4 to 5 times faster than a simulation with comparable accuracy (polynomial order 6). Additionally, we explore changes in the hyperparameters and use transfer learning to speed up the training. We observe that it is not useful to train a corrective forcing using a different flow condition. However, an already trained corrective forcing can be used to initialise a new training (at the correct flow conditions) to obtain an effective forcing with only a few training iterations. △ Less

Submitted 23 July, 2022; originally announced July 2022.

towards a robust detection of viscous and turbulent flow regions using unsupervised machine learning

Authors: Kheir-Eddine Otmani, Gerasimos Ntoukas, Esteban Ferrer

Abstract: We propose an invariant feature space for the detection of viscous dominated and turbulent regions (i.e., boundary layers and wakes). The developed methodology uses the principal invariants of the strain and rotational rate tensors as input to an unsupervised Machine Learning Gaussian mixture model. The selected feature space is independent of the coordinate frame used to generate the processed da… ▽ More We propose an invariant feature space for the detection of viscous dominated and turbulent regions (i.e., boundary layers and wakes). The developed methodology uses the principal invariants of the strain and rotational rate tensors as input to an unsupervised Machine Learning Gaussian mixture model. The selected feature space is independent of the coordinate frame used to generate the processed data, as it relies on the principal invariants of strain and rotational rate, which are Galilean invariants. This methodology allows us to identify two distinct flow regions: a viscous dominated, rotational region (boundary layer and wake region) and an inviscid, irrotational region (outer flow region). We test the methodology on a laminar and a turbulent (using Large Eddy Simulation) case for flows past a circular cylinder at $Re=40$ and $Re=3900$. The simulations have been conducted using a high-order nodal Discontinuous Galerkin Spectral Element Method (DGSEM). The results obtained are analysed to show that Gaussian mixture clustering provides an effective identification method of viscous dominated and rotational regions in the flow. We also include comparisons with traditional sensors to show that the proposed clustering does not depend on the selection of an arbitrary threshold, as required when using traditional sensors. △ Less

Submitted 6 July, 2022; originally announced July 2022.

HORSES3D: a high-order discontinuous Galerkin solver for flow simulations and multi-physics applications

Authors: E. Ferrer, G. Rubio, G. Ntoukas, W. Laskowski, O. A. Mariño, S. Colombo, A. Mateo-Gabín, F. Manrique de Lara, D. Huergo, J. Manzanero, A. M. Rueda-Ramírez, D. A. Kopriva, E. Valero

Abstract: We present the latest developments of our High-Order Spectral Element Solver (HORSES3D), an open source high-order discontinuous Galerkin framework, capable of solving a variety of flow applications, including compressible flows (with or without shocks), incompressible flows, various RANS and LES turbulence models, particle dynamics, multiphase flows, and aeroacoustics. We provide an overview of t… ▽ More We present the latest developments of our High-Order Spectral Element Solver (HORSES3D), an open source high-order discontinuous Galerkin framework, capable of solving a variety of flow applications, including compressible flows (with or without shocks), incompressible flows, various RANS and LES turbulence models, particle dynamics, multiphase flows, and aeroacoustics. We provide an overview of the high-order spatial discretisation (including energy/entropy stable schemes) and anisotropic p-adaptation capabilities. The solver is parallelised using MPI and OpenMP showing good scalability for up to 1000 processors. Temporal discretisations include explicit, implicit, multigrid, and dual time-stepping schemes with efficient preconditioners. Additionally, we facilitate meshing and simulating complex geometries through a mesh-free immersed boundary technique. We detail the available documentation and the test cases included in the GitHub repository. △ Less

Submitted 20 June, 2022; originally announced June 2022.

Speed of Sound for Hadronic and Quark Phases in a Magnetic Field

Authors: E. J. Ferrer, A. Hackebill

Abstract: In this paper we calculate the speed of sound for three phases that may exist inside a magnetized hybrid neutron star at different density regions: A hadronic phase at low densities, quark-matter in the magnetic dual chiral density wave (MDCDW) phase at intermediate densities and a free-quark phase modeled by the MIT bag model at higher densities. It is found that the speed of sound exhibits a non… ▽ More In this paper we calculate the speed of sound for three phases that may exist inside a magnetized hybrid neutron star at different density regions: A hadronic phase at low densities, quark-matter in the magnetic dual chiral density wave (MDCDW) phase at intermediate densities and a free-quark phase modeled by the MIT bag model at higher densities. It is found that the speed of sound exhibits a non-monotonic behavior, that goes from values smaller than the conformal limit ($c_s^2 < 1/3$) in the hadronic phase, to peak ($c_s^2 > 1/3$) in the MDCDW phase, to finally reach the conformal limit ($c_s^2 \sim 1/3$) at higher densities for quarks in the MIT bag model. Also, the anisotropic speed of sound in the presence of a magnetic field is derived from first principles. This is a consequence of the anisotropy in the system's pressures produced by the breaking of the rotational symmetry in the presence of a magnetic field. The role played by the lowest Landau level contribution in affecting the speed of sound in the magnetized phases is discussed. △ Less

Submitted 17 January, 2023; v1 submitted 30 March, 2022; originally announced March 2022.

Comments: Version to appear in Nucl. Phys. A

Gaka-chu: a self-employed autonomous robot artist

Authors: Eduardo Castelló Ferrer, Ivan Berman, Aleksandr Kapitonov, Vadim Manaenko, Makar Chernyaev, Pavel Tarasov

Abstract: The physical autonomy of robots is well understood both theoretically and practically. By contrast, there is almost no research exploring their potential economic autonomy. In this paper, we present the first economically autonomous robot -- a robot able to produce marketable goods while having full control over the use of its generated income. Gaka-chu ("painter" in Japanese) is a 6-axis robot ar… ▽ More The physical autonomy of robots is well understood both theoretically and practically. By contrast, there is almost no research exploring their potential economic autonomy. In this paper, we present the first economically autonomous robot -- a robot able to produce marketable goods while having full control over the use of its generated income. Gaka-chu ("painter" in Japanese) is a 6-axis robot arm that creates paintings of Japanese characters from an autoselected keyword. By using a blockchain-based smart contract, Gaka-chu can autonomously list a painting it made for sale in an online auction. In this transaction, the robot interacts with the human bidders as a peer not as a tool. Using the blockchain-based smart contract, Gaka-chu can then use its income from selling paintings to replenish its resources by autonomously ordering materials from an online art shop. We built the Gaka-chu prototype with an Ethereum-based smart contract and ran a 6-month long experiment, during which the robot created and sold four paintings, simultaneously using its income to purchase supplies and repay initial investors. In this work, we present the results of the experiments conducted and discuss the implications of economically autonomous robots. △ Less

Submitted 13 March, 2023; v1 submitted 7 March, 2022; originally announced March 2022.

Comments: Accepted for publication in 2023 IEEE International Conference on Robotics and Automation (ICRA)

ACM Class: I.2.9; C.2.4; H.5; K.4

Magnetic Dual Chiral Density Wave: A Candidate Quark Matter Phase for the Interior of Neutron Stars

Authors: E. J. Ferrer, V. de la Incera

Abstract: In this review, we discuss the physical characteristics of the magnetic dual chiral density wave (MDCDW) phase of dense quark matter and argued why it is a promising candidate for the interior matter phase of neutron stars. The MDCDW condensate occurs in the presence of a magnetic field. It is a single-modulated chiral density wave characterized by two dynamically generated parameters: the fermion… ▽ More In this review, we discuss the physical characteristics of the magnetic dual chiral density wave (MDCDW) phase of dense quark matter and argued why it is a promising candidate for the interior matter phase of neutron stars. The MDCDW condensate occurs in the presence of a magnetic field. It is a single-modulated chiral density wave characterized by two dynamically generated parameters: the fermion quasiparticle mass $m$ and the condensate spatial modulation $q$. The lowest Landau level quasiparticle modes in the MDCDW system are asymmetric about the zero energy, a fact that leads to the topological properties and anomalous electric transport exhibited by this phase. The topology makes the MDCDW phase robust against thermal phonon fluctuations, and as such, it does not display the Landau-Peierls instability, a stapled feature of single-modulated inhomogeneous chiral condensates in three dimensions. The topology is also reflected in the presence of the electromagnetic chiral anomaly in the effective action and in the formation of hybridized propagating modes known as an axion-polaritons. Taking into account that one of the axion-polaritons of this quark phase is gapped, we argued how incident $γ$-ray photons can be converted into gapped axion-polaritons in the interior of a magnetar star in the MDCDW phase leading the star to collapse, a phenomenon that can serve to explain the so-called missing pulsar problem in the galactic center. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Journal ref: Universe 7 (2021) 12, 458

A combined volume penalization / selective frequency damping approach for immersed boundary methods applied to high-order schemes

Authors: Jiaqing Kou, Esteban Ferrer

Abstract: There has been an increasing interest in developing efficient immersed boundary method (IBM) based on Cartesian grids, recently in the context of high-order methods. IBM based on volume penalization is a robust and easy to implement method to avoid body-fitted meshes and has been recently adapted to high order discretisations (Kou et al., 2021). This work proposes an improvement over the classic p… ▽ More There has been an increasing interest in developing efficient immersed boundary method (IBM) based on Cartesian grids, recently in the context of high-order methods. IBM based on volume penalization is a robust and easy to implement method to avoid body-fitted meshes and has been recently adapted to high order discretisations (Kou et al., 2021). This work proposes an improvement over the classic penalty formulation for flux reconstruction high order solvers. We include a selective frequency damping (SFD) approach (Aakervik et al., 2006) acting only inside solid body defined through the immersed boundary masking, to damp spurious oscillations. An encapsulated formulation for the SFD method is implemented, which can be used as a wrapper around an existing time-stepping code. The numerical properties have been studied through eigensolution analysis based on the advection equation. These studies not only show the advantages of using the SFD method as an alternative of the traditional volume penalization, but also show the favorable properties of combining both approaches. This new approach is then applied to the Navier-Stokes equation to simulate steady flow past an airfoil and unsteady flow past a circular cylinder. The advantages of the SFD method in providing improved accuracy are reported. △ Less

Submitted 21 July, 2021; originally announced July 2021.

Eigensolution analysis of immersed boundary method based on volume penalization: applications to high-order schemes

Authors: Jiaqing Kou, Aurelio Hurtado-de-Mendoza, Saumitra Joshi, Soledad Le Clainche, Esteban Ferrer

Abstract: This paper presents eigensolution and non-modal analyses for immersed boundary methods (IBMs) based on volume penalization for the linear advection equation. This approach is used to analyze the behavior of flux reconstruction (FR) discretization, including the influence of polynomial order and penalization parameter on numerical errors and stability. Through a semi-discrete analysis, we find that… ▽ More This paper presents eigensolution and non-modal analyses for immersed boundary methods (IBMs) based on volume penalization for the linear advection equation. This approach is used to analyze the behavior of flux reconstruction (FR) discretization, including the influence of polynomial order and penalization parameter on numerical errors and stability. Through a semi-discrete analysis, we find that the inclusion of IBM adds additional dissipation without changing significantly the dispersion of the original numerical discretization. This agrees with the physical intuition that in this type of approach, the solid wall is modelled as a porous medium with vanishing viscosity. From a stability point view, the selection of penalty parameter can be analyzed based on a fully-discrete analysis, which leads to practical guidelines on the selection of penalization parameter. Numerical experiments indicate that the penalization term needs to be increased to damp oscillations inside the solid (i.e. porous region), which leads to undesirable time step restrictions. As an alternative, we propose to include a second-order term in the solid for the no-slip wall boundary condition. Results show that by adding a second-order term we improve the overall accuracy with relaxed time step restriction. This indicates that the optimal value of the penalization parameter and the second-order damping can be carefully chosen to obtain a more accurate scheme. Finally, the approximated relationship between these two parameters is obtained and used as a guideline to select the optimum penalty terms in a Navier-Stokes solver, to simulate flow past a cylinder. △ Less

Submitted 21 July, 2021; originally announced July 2021.

Tidal deformability of strange stars and the GW170817 event

Authors: Odilon Lourenço, César H. Lenzi, Mariana Dutra, Efrain J. Ferrer, Vivian de la Incera, Laura Paulucci, J. E. Horvath

Abstract: In this work we consider strange stars formed by quark matter in the color-flavor-locked (CFL) phase of color superconductivity. The CFL phase is described by a Nambu-Jona-Lasinio model with four-fermion vector and diquark interaction channels. The effect of the color superconducting medium on the gluons are incorporated into the model by including the gluon self-energy in the thermodynamic potent… ▽ More In this work we consider strange stars formed by quark matter in the color-flavor-locked (CFL) phase of color superconductivity. The CFL phase is described by a Nambu-Jona-Lasinio model with four-fermion vector and diquark interaction channels. The effect of the color superconducting medium on the gluons are incorporated into the model by including the gluon self-energy in the thermodynamic potential. We construct parametrizations of the model by varying the vector coupling $G_V$ and comparing the results to the data on tidal deformability from the GW170817 event, the observational data on maximum masses from massive pulsars such as the MSP J0740+6620, and the mass/radius fits to NICER data for PSR J003+0451. Our results points out to windows for the $G_V$ parameter space of the model, with and without gluon effects included, that are compatible with all these astrophysical constraints, namely, $0.21<G_V/G_S<0.4$, and $0.02<G_V/G_S<0.1$, respectively. We also observe a strong correlation between the tidal deformabilites of the GW170817 event and $G_V$. Our results indicate that strange stars cannot be ruled out in collisions of compact binaries from the structural point of view. △ Less

Submitted 15 April, 2021; originally announced April 2021.

Comments: 9 pages, 8 figures. Accepted for publication in Phys. Rev. D

Journal ref: Physical Review D, 2021

Quark matter contribution to the heat capacity of magnetized neutron stars

Authors: E. J. Ferrer, V. de la Incera, P. Sanson

Abstract: In this paper, we find the heat capacity of the magnetic dual chiral density wave (MDCDW) phase of dense quark matter and use it to explore the feasibility of this phase for a neutron star interior. MDCDW is a spatially inhomogeneous phase of quark matter known to be favored at intermediate densities over the chirally symmetric phase and the color-flavor-locked superconducting phase. By comparing… ▽ More In this paper, we find the heat capacity of the magnetic dual chiral density wave (MDCDW) phase of dense quark matter and use it to explore the feasibility of this phase for a neutron star interior. MDCDW is a spatially inhomogeneous phase of quark matter known to be favored at intermediate densities over the chirally symmetric phase and the color-flavor-locked superconducting phase. By comparing our result to the lower limit of the core heat capacity established from observations of transiently accreting neutron stars, we show that the heat capacity of MDCDW quark matter is well above that lower limit and hence cannot be ruled out. This result adds to a wealth of complementary investigations, all of which has served to strengthen the viability of a neutron star interior made of MDCDW quark matter. For completeness, we review the contributions to the heat capacity of the main neutron star ingredients at low, high and intermediate densities, with and without the presence of a magnetic field. △ Less

Submitted 3 June, 2021; v1 submitted 11 January, 2021; originally announced January 2021.

Comments: To appear in PRD

Journal ref: Phys. Rev. D 103, 123013 (2021)

Hadron-Quark Phase Transition at Finite Density in the Presence of a Magnetic Field: Anisotropic Approach

Authors: E. J. Ferrer, A. Hackebill

Abstract: We investigate the hadron-quark phase transition at finite density in the presence of a magnetic field taking into account the anisotropy created by a uniform magnetic field in the system's equations of state. We find a new anisotropic equilibrium condition that will drive the first-order phase transition along the boundary between the two phases. Fixing the magnetic field in the hadronic phase, t… ▽ More We investigate the hadron-quark phase transition at finite density in the presence of a magnetic field taking into account the anisotropy created by a uniform magnetic field in the system's equations of state. We find a new anisotropic equilibrium condition that will drive the first-order phase transition along the boundary between the two phases. Fixing the magnetic field in the hadronic phase, the phase transition is realized by increasing the baryonic chemical potential at zero temperature. It is shown that the magnetic field is mildly boosted after the system transitions from the hadronic to the quark phase. The magnetic-field discontinuity between the two phases is supported by a surface density of magnetic monopoles, which accumulate at the boundary separating the two phases. The mechanism responsible for the monopole charge density generation is discussed. Each phase is found to be paramagnetic with higher magnetic susceptibility in the quark phase. The connection with the physics of neutron stars is highlighted through out the paper. △ Less

Submitted 17 March, 2022; v1 submitted 20 October, 2020; originally announced October 2020.

Comments: 30 pages and 6 figures

Axion-Polaritons in the Magnetic Dual Chiral Density Wave Phase of Dense QCD

Authors: Efrain J. Ferrer, Vivian de la Incera

Abstract: We investigate the propagation of electromagnetic radiation in the magnetic dual chiral density wave (MDCDW) phase of dense quark matter. Considering the theory of low-energy fluctuations in this phase, we show how linearly polarized photons reaching the MDCDW medium couple to the fluctuation field to produce two hybridized modes of propagation that we call in analogy with similar phenomenon in co… ▽ More We investigate the propagation of electromagnetic radiation in the magnetic dual chiral density wave (MDCDW) phase of dense quark matter. Considering the theory of low-energy fluctuations in this phase, we show how linearly polarized photons reaching the MDCDW medium couple to the fluctuation field to produce two hybridized modes of propagation that we call in analogy with similar phenomenon in condensed matter physics axion polaritons, one of them being gapless and the other gapped. The gapped mode's gap is proportional to the background magnetic field and inversely proportional to the amplitude of the inhomogeneous condensate. The generation of axion polaritons can be traced back to the presence of the chiral anomaly in the low-energy theory of the fluctuations. Considering the Primakoff effect in the MDCDW medium, we argued that axion polaritons can be generated inside quark stars bombarded by energetic photons coming from gamma-ray bursts and point out that this mechanism could serve to explain the missing pulsar paradox in the galaxy center. △ Less

Submitted 15 August, 2023; v1 submitted 5 October, 2020; originally announced October 2020.

Comments: Some formulas were changed

Journal ref: Nucl. Phys. B 994 (2023) 116307

Secure Encoded Instruction Graphs for End-to-End Data Validation in Autonomous Robots

Authors: Jorge Peña Queralta, Li Qingqing, Eduardo Castelló Ferrer, Tomi Westerlund

Abstract: As autonomous robots are becoming more widespread, more attention is being paid to the security of robotic operation. Autonomous robots can be seen as cyber-physical systems: they can operate in virtual, physical, and human realms. Therefore, securing the operations of autonomous robots requires not only securing their data (e.g., sensor inputs and mission instructions) but securing their interact… ▽ More As autonomous robots are becoming more widespread, more attention is being paid to the security of robotic operation. Autonomous robots can be seen as cyber-physical systems: they can operate in virtual, physical, and human realms. Therefore, securing the operations of autonomous robots requires not only securing their data (e.g., sensor inputs and mission instructions) but securing their interactions with their environment. There is currently a deficiency of methods that would allow robots to securely ensure their sensors and actuators are operating correctly without external feedback. This paper introduces an encoding method and end-to-end validation framework for the missions of autonomous robots. In particular, we present a proof of concept of a map encoding method, which allows robots to navigate realistic environments and validate operational instructions with almost zero {\it a priori} knowledge. We demonstrate our framework using two different encoded maps in experiments with simulated and real robots. Our encoded maps have the same advantages as typical landmark-based navigation, but with the added benefit of cryptographic hashes that enable end-to-end information validation. Our method is applicable to any aspect of robotic operation in which there is a predefined set of actions or instructions given to the robot. △ Less

Submitted 10 April, 2022; v1 submitted 2 September, 2020; originally announced September 2020.

Comments: To be published in the IEEE Internet of Things Journal

Foundational Oracle Patterns: Connecting Blockchain to the Off-chain World

Authors: Roman Mühlberger, Stefan Bachhofner, Eduardo Castelló Ferrer, Claudio Di Ciccio, Ingo Weber, Maximilian Wöhrer, Uwe Zdun

Abstract: Blockchain has evolved into a platform for decentralized applications, with beneficial properties like high integrity, transparency, and resilience against censorship and tampering. However, blockchains are closed-world systems which do not have access to external state. To overcome this limitation, oracles have been introduced in various forms and for different purposes. However so far common ora… ▽ More Blockchain has evolved into a platform for decentralized applications, with beneficial properties like high integrity, transparency, and resilience against censorship and tampering. However, blockchains are closed-world systems which do not have access to external state. To overcome this limitation, oracles have been introduced in various forms and for different purposes. However so far common oracle best practices have not been dissected, classified, and studied in their fundamental aspects. In this paper, we address this gap by studying foundational blockchain oracle patterns in two foundational dimensions characterising the oracles: (i) the data flow direction, i.e., inbound and outbound data flow, from the viewpoint of the blockchain; and (ii) the initiator of the data flow, i.e., whether it is push or pull-based communication. We provide a structured description of the four patterns in detail, and discuss an implementation of these patterns based on use cases. On this basis we conduct a quantitative analysis, which results in the insight that the four different patterns are characterized by distinct performance and costs profiles. △ Less

Submitted 29 July, 2020; originally announced July 2020.

Lack of Debye and Meissner screening in strongly magnetized quark matter at intermediate densities

Authors: Bo Feng, Efrain J. Ferrer, Israel Portillo

Abstract: We study the static responses of cold quark matter in the intermediate baryonic density region (characterized by a chemical potential $μ$) in the presence of a strong magnetic field. We consider in particular, the so-called Magnetic Dual Chiral Density Wave (MDCDW) phase, which is materialized by an inhomogeneous condensate formed by a particle-hole pair. It is shown, that the MDCDW phase is more… ▽ More We study the static responses of cold quark matter in the intermediate baryonic density region (characterized by a chemical potential $μ$) in the presence of a strong magnetic field. We consider in particular, the so-called Magnetic Dual Chiral Density Wave (MDCDW) phase, which is materialized by an inhomogeneous condensate formed by a particle-hole pair. It is shown, that the MDCDW phase is more stable in the weak-coupling regime than the one considered in the magnetic catalysis of chiral symmetry braking phenomenon (MC$χ$SB) and even than the chiral symmetric phase that was expected to be realized at sufficiently high baryonic chemical potential. The different components of the photon polarization operator of the MDCDW phase in the one-loop approximation are calculated. We found that in the MDCDW phase there is no Debye screening neither Meissner effect in the lowest-Landau-level approximation. The obtained Debye length depends on the amplitude $m$ and modulation $b$ of the inhomogeneous condensate and it is only different from zero if the relation $| μ-b| > m$ holds. But, we found that in the region of interest this inequality is not satisfied. Thus, no Debye screening takes place under those conditions. On the other hand, since the particle-hole condensate is electrically neutral, the U(1) electromagnetic group is not broken by the ground state and consequently there is no Meissner effect. These results can be of interest for the astrophysics of neutron stars. △ Less

Submitted 9 April, 2020; v1 submitted 8 January, 2020; originally announced January 2020.

Comments: 18 pages, 5 figures

Journal ref: Phys. Rev. D 101, 056012 (2020)

A Statically Condensed Discontinuous Galerkin Spectral Element Method on Gauss-Lobatto Nodes for the Compressible Navier-Stokes Equations

Authors: Andrés M. Rueda-Ramírez, Esteban Ferrer, David A. Kopriva, Gonzalo Rubio, Eusebio Valero

Abstract: We present a static-condensation method for time-implicit discretizations of the Discontinuous Galerkin Spectral Element Method on Gauss-Lobatto points (GL-DGSEM). We show that, when solving the compressible Navier-Stokes equations, it is possible to reorganize the linear system that results from the implicit time-integration of the GL-DGSEM as a Schur complement problem, which can be efficiently… ▽ More We present a static-condensation method for time-implicit discretizations of the Discontinuous Galerkin Spectral Element Method on Gauss-Lobatto points (GL-DGSEM). We show that, when solving the compressible Navier-Stokes equations, it is possible to reorganize the linear system that results from the implicit time-integration of the GL-DGSEM as a Schur complement problem, which can be efficiently solved using static condensation. The use of static condensation reduces the linear system size and improves the condition number of the system matrix, which translates into shorter computational times when using direct and iterative solvers. The statically condensed GL-DGSEM presented here can be applied to linear and nonlinear advection-diffusion partial differential equations in conservation form. To test it we solve the compressible Navier-Stokes equations with direct and Krylov subspace solvers, and we show for a selected problem that using the statically condensed GL-DGSEM leads to speed-ups of up to $200$ when compared to the time-explicit GL-DGSEM, and speed-ups of up to three when compared with the time-implicit GL-DGSEM that solves the global system. The GL-DGSEM has gained increasing popularity in recent years because it satisfies the summation-by-parts property, which enables the construction of provably entropy stable schemes, and because it is computationally very efficient. In this paper, we show that the GL-DGSEM has an additional advantage: It can be statically condensed. △ Less

Submitted 12 December, 2019; v1 submitted 3 November, 2019; originally announced November 2019.

Entropy-stable discontinuous Galerkin approximation with summation-by-parts property for the incompressible Navier-Stokes/Cahn-Hilliard system

Authors: Juan Manzanero, Gonzalo Rubio, David A. Kopriva, Esteban Ferrer, Eusebio Valero

Abstract: We develop an entropy stable two-phase incompressible Navier--Stokes/Cahn--Hilliard discontinuous Galerkin (DG) flow solver method. The model poses the Cahn-Hilliard equation as the phase field method, a skew-symmetric form of the momentum equation, and an artificial compressibility method to compute the pressure. We design the model so that it satisfies an entropy law, including free- and no-slip… ▽ More We develop an entropy stable two-phase incompressible Navier--Stokes/Cahn--Hilliard discontinuous Galerkin (DG) flow solver method. The model poses the Cahn-Hilliard equation as the phase field method, a skew-symmetric form of the momentum equation, and an artificial compressibility method to compute the pressure. We design the model so that it satisfies an entropy law, including free- and no-slip wall boundary conditions with non-zero wall contact angle. We then construct a high-order DG approximation of the model that satisfies the SBP-SAT property. With the help of a discrete stability analysis, the scheme has two modes: an entropy conserving approximation with central advective fluxes and the Bassi-Rebay 1 (BR1) method for diffusion, and an entropy stable approximation with an exact Riemann solver for advection and interface stabilization added to the BR1 method. The scheme is applicable to, and the stability proofs hold for, three-dimensional unstructured meshes with curvilinear hexahedral elements. We test the convergence of the schemes on a manufactured solution, and their robustness by solving a flow initialized from random numbers. In the latter, we find that a similar scheme that does not satisfy an entropy inequality had 30% probability to fail, while the entropy stable scheme never does. We also solve the static and rising bubble test problems, and to challenge the solver capabilities we compute a three-dimensional pipe flow in the annular regime. △ Less

Submitted 25 October, 2019; v1 submitted 24 October, 2019; originally announced October 2019.

Magnetic Field Effect in the Fine-Structure Constant and Electron Dynamical Mass

Authors: E. J. Ferrer, A. Sanchez

Abstract: We investigate the effect of an applied constant and uniform magnetic field in the fine-structure constant of massive and massless QED. In massive QED, it is shown that a strong magnetic field removes the so called Landau pole and that the fine-structure constant becomes anisotropic having different values along and transverse to the field direction. Contrary to other results in the literature, we… ▽ More We investigate the effect of an applied constant and uniform magnetic field in the fine-structure constant of massive and massless QED. In massive QED, it is shown that a strong magnetic field removes the so called Landau pole and that the fine-structure constant becomes anisotropic having different values along and transverse to the field direction. Contrary to other results in the literature, we find that the anisotropic fine-structure constant always decreases with the field. We also study the effect of the running of the coupling constant with the magnetic field on the electron mass. We find that in both cases of massive and massless QED, the electron dynamical mass always decreases with the magnetic field, what can be interpreted as an inverse magnetic catalysis effect. △ Less

Submitted 24 October, 2019; originally announced October 2019.

Comments: To appear in Phys. Rev. D

Journal ref: Phys. Rev. D 100, 096006 (2019)

InfraRed Astronomy Satellite Swarm Interferometry (IRASSI): Overview and Study Results

Authors: Hendrik Linz, Divya Bhatia, Luisa Buinhas, Matthias Lezius, Eloi Ferrer, Roger Förstner, Kathrin Frankl, Mathias Philips-Blum, Meiko Steen, Ulf Bestmann, Wolfgang Hänsel, Ronald Holzwarth, Oliver Krause, Thomas Pany

Abstract: The far-infrared (FIR) is one of the few wavelength ranges where no astronomical data with sub-arcsec resolution exist yet. Neither of the medium-term satellite projects like SPICA, Millimetron or OST will resolve this malady. Information at high spatial and spectral resolution in the FIR, taken from atomic fine-structure lines, highly excited CO, and especially from water lines would, however, op… ▽ More The far-infrared (FIR) is one of the few wavelength ranges where no astronomical data with sub-arcsec resolution exist yet. Neither of the medium-term satellite projects like SPICA, Millimetron or OST will resolve this malady. Information at high spatial and spectral resolution in the FIR, taken from atomic fine-structure lines, highly excited CO, and especially from water lines would, however, open the door for transformative science. This calls for interferometric concepts. We present first results of our feasibility study IRASSI (Infrared Astronomy Satellite Swarm Interferometry) for a FIR space interferometer. Extending on the principal concept of the ESPRIT study, it features heterodyne interferometry within a swarm of 5 satellite elements. The satellites can drift in and out within a range of several hundred meters, thereby achieving spatial resolutions of <0.1" over the whole wavelength range of 1-6 THz. Precise knowledge on the baselines will be ensured by metrology methods employing laser-based optical frequency combs, for which preliminary ground-based tests have been designed by us. We first show how the science requirements translate into operational and design parameters. We have put much emphasis on the navigational aspects of such a free-flying satellite swarm operating in relatively close vicinity. We hence present work on the formation geometry, the relative dynamics of the swarm, and aspects of our investigation towards attitude estimation. Furthermore, we discuss issues regarding the real-time capability of the autonomous relative positioning system, which is an important aspect for IRASSI where, due to the large raw data rates expected, the interferometric correlation has to be done onboard. We also address questions regarding the spacecraft architecture and how a thermomechanical model is used to study the effect of thermal perturbations on the spacecraft. (abridged) △ Less

Submitted 18 July, 2019; originally announced July 2019.

Comments: 43 pages (preprint format), 10 figures, 6 tables; in press, to appear in Advances in Space Research, Special Issue on Space-VLBI (eds. Leonid Gurvits et al.)

Entropy-stable discontinuous Galerkin approximation with summation-by-parts property for the incompressible Navier-Stokes equations with variable density and artificial compressibility

Authors: Juan Manzanero, Gonzalo Rubio, David A Kopriva, Esteban Ferrer, Eusebio Valero

Abstract: We present a provably stable discontinuous Galerkin spectral element method for the incompressible Navier-Stokes equations with artificial compressibility and variable density. Stability proofs, which include boundary conditions, that follow a continuous entropy analysis are provided. We define a mathematical entropy function that combines the traditional kinetic energy and an additional energy te… ▽ More We present a provably stable discontinuous Galerkin spectral element method for the incompressible Navier-Stokes equations with artificial compressibility and variable density. Stability proofs, which include boundary conditions, that follow a continuous entropy analysis are provided. We define a mathematical entropy function that combines the traditional kinetic energy and an additional energy term for the artificial compressiblity, and derive its associated entropy conservation law. The latter allows us to construct a provably stable split-form nodal Discontinuous Galerkin (DG) approximation that satisfies the summation-by-parts simultaneous-approximation-term (SBP-SAT) property. The scheme and the stability proof are presented for general curvilinear three-dimensional hexahedral meshes. We use the exact Riemann solver and the Bassi-Rebay 1 (BR1) scheme at the inter-element boundaries for inviscid and viscous fluxes respectively, and an explicit low storage Runge-Kutta RK3 scheme to integrate in time. We assess the accuracy and robustness of the method by solving the Kovasznay flow, the inviscid Taylor-Green vortex, and the Rayleigh-Taylor instability. △ Less

Submitted 16 July, 2019; v1 submitted 12 July, 2019; originally announced July 2019.

Secure and secret cooperation in robotic swarms

Authors: Eduardo Castelló Ferrer, Thomas Hardjono, Alex 'Sandy' Pentland, Marco Dorigo

Abstract: The importance of swarm robotics systems in both academic research and real-world applications is steadily increasing. However, to reach widespread adoption, new models that ensure the secure cooperation of large groups of robots need to be developed. This work introduces a novel method to encapsulate cooperative robotic missions in an authenticated data structure known as Merkle tree. With this m… ▽ More The importance of swarm robotics systems in both academic research and real-world applications is steadily increasing. However, to reach widespread adoption, new models that ensure the secure cooperation of large groups of robots need to be developed. This work introduces a novel method to encapsulate cooperative robotic missions in an authenticated data structure known as Merkle tree. With this method, operators can provide the "blueprint" of the swarm's mission without disclosing its raw data. In other words, data verification can be separated from data itself. We propose a system where robots in a swarm, to cooperate towards mission completion, have to "prove" their integrity to their peers by exchanging cryptographic proofs. We show the implications of this approach for two different swarm robotics missions: foraging and maze formation. In both missions, swarm robots were able to cooperate and carry out sequential operations without having explicit knowledge about the mission's high-level objectives. The results presented in this work demonstrate the feasibility of using Merkle trees as a cooperation mechanism for swarm robotics systems in both simulation and real-robot experiments, which has implications for future decentralized robotics applications where security plays a crucial role such as environmental monitoring, infrastructure surveillance, and disaster management. △ Less

Submitted 1 April, 2021; v1 submitted 19 April, 2019; originally announced April 2019.

Comments: 32 pages, 6 figures, submitted to Science Robotics

Thermodynamics of Neutrons in a Magnetic Field and its Implications for Neutron Stars

Authors: E. J. Ferrer, A. Hackebill

Abstract: We investigate the effects of a magnetic field on the thermodynamics of a neutron system at finite density and temperature. Our main motivation is to deepen the understanding of the physics of a class of neutron stars known as magnetars, which exhibit extremely strong magnetic fields. Taking into account two facts, (i) the existence of a pressure anisotropy in the presence of a magnetic field and… ▽ More We investigate the effects of a magnetic field on the thermodynamics of a neutron system at finite density and temperature. Our main motivation is to deepen the understanding of the physics of a class of neutron stars known as magnetars, which exhibit extremely strong magnetic fields. Taking into account two facts, (i) the existence of a pressure anisotropy in the presence of a magnetic field and (ii) that the quantum field theory contribution to the pressure is non-negligible, we show that the maximum value that the inner magnetic field of a star can reach while being in agreement with the magnetohydrostatic equilibrium between the gravitational and matter pressures becomes $10^{17}$ G, an order of magnitude smaller than the previous value obtained through the scalar virial theorem; that the magnetic field has a negligible effect on the neutron system's equation of state; that the system's magnetic susceptibility increases with the temperature; and that the specific heat $C_V$ does not significantly change with the magnetic field in the range of temperatures characteristic of protoneutron stars. △ Less

Submitted 4 June, 2019; v1 submitted 19 March, 2019; originally announced March 2019.

Journal ref: Phys. Rev. C 99, 065803 (2019)

A free-energy stable nodal discontinuous Galerkin approximation with summation-by-parts property for the Cahn-Hilliard equation

Authors: Juan Manzanero, Gonzalo Rubio, David A. Kopriva, Esteban Ferrer, Eusebio Valero

Abstract: We present a nodal Discontinuous Galerkin (DG) scheme for the Cahn-Hilliard equation that satisfies the summation-by-parts simultaneous-approximation-term (SBP-SAT) property. The latter permits us to show that the discrete free-energy is bounded, and as a result, the scheme is provably stable. The scheme and the stability proof are presented for general curvilinear three-dimensional hexahedral mes… ▽ More We present a nodal Discontinuous Galerkin (DG) scheme for the Cahn-Hilliard equation that satisfies the summation-by-parts simultaneous-approximation-term (SBP-SAT) property. The latter permits us to show that the discrete free-energy is bounded, and as a result, the scheme is provably stable. The scheme and the stability proof are presented for general curvilinear three-dimensional hexahedral meshes. We use the Bassi-Rebay 1 (BR1) scheme to compute interface fluxes, and an IMplicit-EXplicit (IMEX) scheme to integrate in time. Lastly, we test the theoretical findings numerically and present examples for two and three-dimensional problems. △ Less

Submitted 21 February, 2019; originally announced February 2019.

Absence of Landau-Peierls Instability in the Magnetic Dual Chiral Density Wave Phase of Dense QCD

Authors: Efrain J Ferrer, Vivian de la Incera

Abstract: We investigate the stability of the Magnetic Dual Chiral Density Wave (MDCDW) phase of cold and dense QCD against collective low-energy fluctuations of the order parameter. The appearance of additional structures in the system free-energy due to the explicit breaking of the rotational and isospin symmetries by the external magnetic field and the field-induced asymmetry of the lowest Landau level m… ▽ More We investigate the stability of the Magnetic Dual Chiral Density Wave (MDCDW) phase of cold and dense QCD against collective low-energy fluctuations of the order parameter. The appearance of additional structures in the system free-energy due to the explicit breaking of the rotational and isospin symmetries by the external magnetic field and the field-induced asymmetry of the lowest Landau level modes play a crucial role in the analysis. The new structures not only affect the condensate minimum equations, but also the spectrum of the thermal fluctuations, which lacks the transverse soft modes that typically affect single-modulated inhomogeneous phases in the absence of a magnetic field. Consequently, the long-range order of the MDCDW phase is preserved at finite temperature. The lack of Landau-Peierls instabilities in the MDCDW phase makes this inhomogeneous phase of dense quark matter particularly relevant for the physics of neutron stars. △ Less

Submitted 9 May, 2020; v1 submitted 18 February, 2019; originally announced February 2019.

Comments: Typos corrected, new discussions added

Journal ref: Phys. Rev. D 102, 014010 (2020)

Urban Swarms: A new approach for autonomous waste management

Authors: Antonio Luca Alfeo, Eduardo Castelló Ferrer, Yago Lizarribar Carrillo, Arnaud Grignard, Luis Alonso Pastor, Dylan T. Sleeper, Mario G. C. A. Cimino, Bruno Lepri, Gigliola Vaglini, Kent Larson, Marco Dorigo, Alex `Sandy' Pentland

Abstract: Modern cities are growing ecosystems that face new challenges due to the increasing population demands. One of the many problems they face nowadays is waste management, which has become a pressing issue requiring new solutions. Swarm robotics systems have been attracting an increasing amount of attention in the past years and they are expected to become one of the main driving factors for innovati… ▽ More Modern cities are growing ecosystems that face new challenges due to the increasing population demands. One of the many problems they face nowadays is waste management, which has become a pressing issue requiring new solutions. Swarm robotics systems have been attracting an increasing amount of attention in the past years and they are expected to become one of the main driving factors for innovation in the field of robotics. The research presented in this paper explores the feasibility of a swarm robotics system in an urban environment. By using bio-inspired foraging methods such as multi-place foraging and stigmergy-based navigation, a swarm of robots is able to improve the efficiency and autonomy of the urban waste management system in a realistic scenario. To achieve this, a diverse set of simulation experiments was conducted using real-world GIS data and implementing different garbage collection scenarios driven by robot swarms. Results presented in this research show that the proposed system outperforms current approaches. Moreover, results not only show the efficiency of our solution, but also give insights about how to design and customize these systems. △ Less

Submitted 1 March, 2019; v1 submitted 18 October, 2018; originally announced October 2018.

Comments: Manuscript accepted for publication in IEEE ICRA 2019

No Net Charge Separation in Hot QCD in a Magnetic Field

Authors: E. J. Ferrer, V de la Incera

Abstract: We study the realization of axion electrodynamics in QCD in the presence of a background magnetic field at temperatures high enough for the occurrence of topological charge transitions that are reflected in the presence of a $θ$-vacuum term in the action. We show that in this system, the Maxwell equations contain two equal and opposite electric currents that are proportional to the time derivative… ▽ More We study the realization of axion electrodynamics in QCD in the presence of a background magnetic field at temperatures high enough for the occurrence of topological charge transitions that are reflected in the presence of a $θ$-vacuum term in the action. We show that in this system, the Maxwell equations contain two equal and opposite electric currents that are proportional to the time derivative of the axion field $θ$. One of these currents comes directly from the Abelian chiral anomaly term in the action and can be interpreted as a polarization current due to the magnetoelectricity of the system with CP-broken symmetry. The other current is obtained from the regular tadpole diagrams and can be understood as produced by the medium chiral imbalance and the single spin projection of the quarks in the lowest Landau level. Since the two currents cancel out, the net electric charge separation along the magnetic field, a phenomenon known as the Chiral Magnetic Effect, does not take place in hight-T QCD at least in equilibrium, in sharp contrast with many claims in the literature. We discuss the similarities and differences with Weyl semimetals in a magnetic field. △ Less

Submitted 14 July, 2018; originally announced July 2018.

Comments: 17 pages

Journal ref: Phys. Rev. D 98, 074009 (2018)

A p-Multigrid Strategy with Anisotropic p-Adaptation Based on Truncation Errors for High-Order Discontinuous Galerkin Methods

Authors: Andrés M. Rueda-Ramírez, Juan Manzanero, Esteban Ferrer, Gonzalo Rubio, Eusebio Valero

Abstract: High-order DG methods have become a popular technique in computational fluid dynamics because their accuracy increases spectrally in smooth solutions with the order of the approximation. However, their main drawback is that increasing the order also increases the computational cost. Several techniques have been introduced in the past to reduce this cost. On the one hand, local mesh adaptation stra… ▽ More High-order DG methods have become a popular technique in computational fluid dynamics because their accuracy increases spectrally in smooth solutions with the order of the approximation. However, their main drawback is that increasing the order also increases the computational cost. Several techniques have been introduced in the past to reduce this cost. On the one hand, local mesh adaptation strategies based on error estimation have been proposed to reduce the number of degrees of freedom while keeping a similar accuracy. On the other hand, multigrid solvers may accelerate time marching computations for a fixed number of degrees of freedom. In this paper, we combine both methods and present a novel anisotropic p-adaptation multigrid algorithm for steady-state problems that uses the multigrid scheme both as a solver and as an anisotropic error estimator. To achieve this, we show that a recently developed anisotropic truncation error estimator [\textit{Rueda-Ramírez et al., Truncation Error Estimation in the p-Anisotropic DGSEM, J. of Scientific Computing}] is perfectly suited to be performed inside the multigrid cycle with negligible extra cost. Furthermore, we introduce a multi-stage p-adaptation procedure which reduces the computational time when very accurate results are required. The proposed methods are tested for the compressible Navier-Stokes equations, where we investigate two cases: the 2D flow on a flat plate is studied to assess accuracy and computational cost of the algorithm, where a speed-up of 816 is achieved compared to the traditional explicit method; and the 3D flow around a sphere is simulated and used to test the anisotropic properties of the proposed method, where a speed-up of 152 is achieved compared to the explicit method. The proposed multi-stage procedure achieved a speed-up of 2.6 in comparison to the single-stage method in highly accurate simulations. △ Less

Submitted 29 June, 2018; originally announced June 2018.