-
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.
-
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.
-
aeon: a Python toolkit for learning from time series
Authors:
Matthew Middlehurst,
Ali Ismail-Fawaz,
Antoine Guillaume,
Christopher Holder,
David Guijo Rubio,
Guzal Bulatova,
Leonidas Tsaprounis,
Lukasz Mentel,
Martin Walter,
Patrick Schäfer,
Anthony Bagnall
Abstract:
aeon is a unified Python 3 library for all machine learning tasks involving time series. The package contains modules for time series forecasting, classification, extrinsic regression and clustering, as well as a variety of utilities, transformations and distance measures designed for time series data. aeon also has a number of experimental modules for tasks such as anomaly detection, similarity s…
▽ More
aeon is a unified Python 3 library for all machine learning tasks involving time series. The package contains modules for time series forecasting, classification, extrinsic regression and clustering, as well as a variety of utilities, transformations and distance measures designed for time series data. aeon also has a number of experimental modules for tasks such as anomaly detection, similarity search and segmentation. aeon follows the scikit-learn API as much as possible to help new users and enable easy integration of aeon estimators with useful tools such as model selection and pipelines. It provides a broad library of time series algorithms, including efficient implementations of the very latest advances in research. Using a system of optional dependencies, aeon integrates a wide variety of packages into a single interface while keeping the core framework with minimal dependencies. The package is distributed under the 3-Clause BSD license and is available at https://github.com/ aeon-toolkit/aeon. This version was submitted to the JMLR journal on 02 Nov 2023 for v0.5.0 of aeon. At the time of this preprint aeon has released v0.9.0, and has had substantial changes.
△ Less
Submitted 20 June, 2024;
originally announced June 2024.
-
Fractional domain wall statistics in spin chains with anomalous symmetries
Authors:
Jose Garre Rubio,
Norbert Schuch
Abstract:
We study the statistics of domain wall excitations in quantum spin chains. We focus on systems with finite symmetry groups represented by matrix product unitaries (MPUs), i.e. finite depth quantum circuits. Such symmetries can be anomalous, in which case gapped phases which they support must break the symmetry. The lowest lying excitations of those systems are thus domain wall excitations. We inve…
▽ More
We study the statistics of domain wall excitations in quantum spin chains. We focus on systems with finite symmetry groups represented by matrix product unitaries (MPUs), i.e. finite depth quantum circuits. Such symmetries can be anomalous, in which case gapped phases which they support must break the symmetry. The lowest lying excitations of those systems are thus domain wall excitations. We investigate the behavior of these domain walls under exchange, and find that they can exhibit non-trivial exchange statistics. This statistics is completely determined by the anomaly of the symmetry, and we provide a direct relation between the known classification of MPU symmetry actions on ground states and the domain wall statistics. Already for the simplest case of a $\mathbb Z_2$ symmetry, we obtain that the presence of an anomalous MPU symmetry gives rise to domain wall excitations which behave neither as bosons nor as fermions, but rather exhibit fractional statistics. Finally, we show that the exchange statistics of domain walls is a physically accessible quantity, by devising explicit measurement operators through which it can be determined.
△ Less
Submitted 1 May, 2024;
originally announced May 2024.
-
Classifying symmetric and symmetry-broken spin chain phases with anomalous group actions
Authors:
Jose Garre Rubio,
Andras Molnar,
Yoshiko Ogata
Abstract:
We consider the classification problem of quantum spin chains invariant under local decomposable group actions, covering matrix product unitaries (MPUs), using an operator algebraic approach. We focus on finite group symmetries hosting both symmetric and symmetry broken phases. The local-decomposable group actions we consider have a 3-cocycle class of the symmetry group associated to them. We deri…
▽ More
We consider the classification problem of quantum spin chains invariant under local decomposable group actions, covering matrix product unitaries (MPUs), using an operator algebraic approach. We focus on finite group symmetries hosting both symmetric and symmetry broken phases. The local-decomposable group actions we consider have a 3-cocycle class of the symmetry group associated to them. We derive invariants for our classification that naturally cover one-dimensional symmetry protected topological (SPT) phases. We prove that these invariants coincide with the ones of [J. Garre Rubio et al, Quantum 7, 927 (2023)] using matrix product states (MPSs) techniques, by explicitly working out the GNS representation of MPSs and MPUs, resulting in a useful dictionary between both approaches that could be of independent interest.
△ Less
Submitted 27 March, 2024;
originally announced March 2024.
-
Emergent (2+1)D topological orders from iterative (1+1)D gauging
Authors:
Jose Garre Rubio
Abstract:
Gauging introduces gauge fields in order to localize an existing global symmetry, resulting in a dual global symmetry on the gauge fields that can be gauged again. By iterating the gauging process on spin chains with Abelian group symmetries and arranging the gauge fields in a 2D lattice, the local symmetries become the stabilizer of the $XZZX$-code for any Abelian group. By twisting the gauging m…
▽ More
Gauging introduces gauge fields in order to localize an existing global symmetry, resulting in a dual global symmetry on the gauge fields that can be gauged again. By iterating the gauging process on spin chains with Abelian group symmetries and arranging the gauge fields in a 2D lattice, the local symmetries become the stabilizer of the $XZZX$-code for any Abelian group. By twisting the gauging map we obtain new codes that explicitly confine anyons, which violate an odd number of plaquette terms and whose fusion results in mobile dipole excitations. Our construction naturally realizes any gapped boundary by taking different quantum phases of the initial (1+1)D globally symmetric system. Our method establishes a new route to obtain higher dimensional topological codes from lower ones, to identify their gapped boundaries and their tensor network representations.
△ Less
Submitted 29 May, 2024; v1 submitted 12 March, 2024;
originally announced March 2024.
-
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.
-
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.
-
A closer physico-chemical look to the Layer-by-Layer electrostatic self-assembly of polyelectrolyte multilayers
Authors:
Eduardo Guzman,
Ramon G. Rubio,
Francisco Ortega
Abstract:
The fabrication of polyelectrolyte multilayer films (PEMs) using the Layer-by-Layer (LbL) method is one of the most versatile approaches for manufacturing functional surfaces. This is the result of the possibility to control the assembly process of the LbL films almost at will, by changing the nature of the assembled materials (building blocks), the assembly conditions (pH, ionic strength, tempera…
▽ More
The fabrication of polyelectrolyte multilayer films (PEMs) using the Layer-by-Layer (LbL) method is one of the most versatile approaches for manufacturing functional surfaces. This is the result of the possibility to control the assembly process of the LbL films almost at will, by changing the nature of the assembled materials (building blocks), the assembly conditions (pH, ionic strength, temperature, etc.) or even by changing some other operational parameters which may impact in the structure and physico-chemical properties of the obtained multi-layered films. Therefore, the understanding of the impact of the above mentioned parameters on the assembly process of LbL materials plays a critical role in the potential use of the LbL method for the fabrication of new functional materials with technological interest. This review tries to provide a broad physico-chemical perspective to the study of the fabrication process of PEMs by the LbL method, which allows one to take advantage of the many possibilities offered for this approach on the fabrication of new functional nanomaterials.
△ Less
Submitted 30 January, 2024;
originally announced February 2024.
-
Fabrication of robust capsules by sequential assembly of polyelectrolytes onto charged liposomes
Authors:
Marta Ruano,
Ana Mateos-Maroto,
Francisco Ortega,
Hernan Ritacco,
Jose E. F. Rubio,
Eduardo Guzman,
Ramon G. Rubio
Abstract:
This work presents a simply methodology for coating small unilamellar liposomes bearing different degrees of positive charge with polyelectrolyte multilayers using the sequential Layer-by-Layer (LbL) deposition method. The liposomes were made of mixtures of DOPC (1,2-dioleyl-sn-glycero-3-phosphocoline) and DODAB (dimethyl dioctadecyl ammonium bromide), and coated by alternated layers of the sodium…
▽ More
This work presents a simply methodology for coating small unilamellar liposomes bearing different degrees of positive charge with polyelectrolyte multilayers using the sequential Layer-by-Layer (LbL) deposition method. The liposomes were made of mixtures of DOPC (1,2-dioleyl-sn-glycero-3-phosphocoline) and DODAB (dimethyl dioctadecyl ammonium bromide), and coated by alternated layers of the sodium salt of poly(4-styrenesulfonate), PSS, and poly(allylamine), PAH, as polyanion and polycations, respectively. The results shows that the zeta potential of the liposomes was not very sensitive to the mole fraction of DODAB in the membrane, XD, in the range 0.3< XD<0.8. We were able to coat the liposomes with up to four polymer bilayers. The growth of the capsules size was followed by dynamic light scattering, and in some cases by cryo-TEM, with good agreement between both techniques. The thickness of the layers, measured from the hydrodynamic radius of the coated liposome, depends on the polyelectrolyte used, so that the PSS layers adopt a much more packaged conformation than the PAH layers. An interesting finding is that the PSS amount needed to reach the isoelectric point of the capsules increases linearly with the charge density of the bare liposomes, whereas the amount of PAH does not depend on it. As expected, the preparation of the multilayers has to be done in such a way that when the system is close to the isoelectric point the capsules do not aggregate. For this, we dropped the polyelectrolyte solution quickly, fast stirring and using dilute liposome suspensions. The method is very flexible and not limited to liposomes or polyelectrolyte multilayers, also coatings containing charged nanoparticles can be easily made.
△ Less
Submitted 27 January, 2024;
originally announced January 2024.
-
Physico-chemical study of polymer mixtures formed by a polycation and a zwitterionic copolymer in aqueous solution and upon adsorption onto negatively charged surfaces
Authors:
Laura Fernandez-Pena,
Eduardo Guzman,
Francisco Ortega,
Lionel Bureau,
Fabien Leonforte,
Dandara Velasco,
Ramon G. Rubio,
Gustavo S. Luengo
Abstract:
The adsorption of mixtures of charged polymers onto solid surfaces presents a big interest in different technological and industrial fields, and in particular, in cosmetics. This requires to deepen on the most fundamental physico-chemical bases governing the deposition, which is generally correlated to the interactions occurring in solution. This work explores the interaction in solution of model…
▽ More
The adsorption of mixtures of charged polymers onto solid surfaces presents a big interest in different technological and industrial fields, and in particular, in cosmetics. This requires to deepen on the most fundamental physico-chemical bases governing the deposition, which is generally correlated to the interactions occurring in solution. This work explores the interaction in solution of model polymer mixtures formed by a cationic homopolymer (poly(diallyl-dimethyl-ammonium chloride), PDADMAC) and a zwitterionic copolymer (copolymer of acrylic acid, 3-Trimethylammonium propyl methacrylamide chloride and acrylamide, Merquat 2003), and the adsorption of such mixtures onto negatively charged surfaces. The analysis of the interactions occurring in solution between both polymers performed using dynamic light scattering (DLS), electrophoretic mobility and viscosity measurements, combined with the study of the deposition of the layers of mixtures containing different weight fractions of each polymer using ellipsometry and quartz crystal microbalance with dissipation monitoring (QCM-D) has shown that the interpolymer complexes formed in solution, and their composition, governs the deposition onto the solid surface and the tribological properties of the adsorbed layers as shown the Surface Force Apparatus (SFA) experiments, allowing for a control of the physico-chemical properties and structure of the layers. Furthermore, the use of Self Consistent Mean Fields Calculations (SCF) confirms the picture obtained from the experimental studies of the adsorbed layers, providing a prediction of the distribution of the polymer chains within the adsorbed layers. It is expected that this study can help on the understanding of the correlations existing between the behavior of future associations of innovative and eco-sustainable polymers and their adsorption processes onto solid surface.
△ Less
Submitted 27 January, 2024;
originally announced January 2024.
-
Study of the Liquid/Vapor Interfacial Properties of Concentrated Polyelectrolyte-Surfactant Mixtures Using Surface Tensiometry and Neutron Reflectometry: Equilibrium, Adsorption Kinetics and Dilational Rheology
Authors:
Sara Llamas,
Eduardo Guzman,
Andrew Akanno,
Laura Fernandez-Pena,
Francisco Ortega,
Richard A. Campbell,
Reinhard Miller,
Ramon G. Rubio
Abstract:
The adsorption of mixtures of poly(diallyldimethylammonium chloride) (PDADMAC) and sodium N-lauroyl-N-methyltaurate (SLMT) at the water / vapor interface has been studied using drop profile tensiometry and neutron reflectometry. This study sheds light on the mechanisms involved in the adsorption of polyelectrolyte-oppositely charged surfactants by the characterization of both equilibrium and dynam…
▽ More
The adsorption of mixtures of poly(diallyldimethylammonium chloride) (PDADMAC) and sodium N-lauroyl-N-methyltaurate (SLMT) at the water / vapor interface has been studied using drop profile tensiometry and neutron reflectometry. This study sheds light on the mechanisms involved in the adsorption of polyelectrolyte-oppositely charged surfactants by the characterization of both equilibrium and dynamics features associated with the layer formation at the fluid interface. The results are discussed in terms of an adsorption-equilibration of the interfacial layers as a two-step process: the initial stages involve the adsorption of polyelectrolyte-surfactant complexes formed in the bulk, and a subsequent stage involves reorganization of the interface. This work contributes to the understanding of the physico-chemical features of systems that undergo complex bulk and interfacial interactions with importance in science and technology.
△ Less
Submitted 27 January, 2024;
originally announced January 2024.
-
Equilibration of a Polycation-Anionic Surfactant Mixture at the Water-Vapor Interface
Authors:
Andrew Akanno,
Eduardo Guzman,
Laura Fernandez-Pena,
Sara Llamas,
Francisco Ortega,
Ramon G. Rubio
Abstract:
The adsorption of concentrated polydiallyldimethylammonium chloride, PDADMAC-sodium lauryl ether sulfate (SLES) mixtures at the water-vapor interface has been studied by different surface tension techniques and dilational viscoelasticity measurements. This work tries to shed light on the way in which the formation of polyelectrolyte-surfactant complexes in the bulk affects the interfacial properti…
▽ More
The adsorption of concentrated polydiallyldimethylammonium chloride, PDADMAC-sodium lauryl ether sulfate (SLES) mixtures at the water-vapor interface has been studied by different surface tension techniques and dilational viscoelasticity measurements. This work tries to shed light on the way in which the formation of polyelectrolyte-surfactant complexes in the bulk affects the interfacial properties of mixtures formed by a polycation and an oppositely charged surfactant. The results are discussed in terms of a two-step adsorption-equilibration of PDADMAC-SLES complexes at the interface, with the initial stages involving the diffusion of kinetically trapped aggregates formed in the bulk to the interface followed by the dissociation and spreading of such aggregates at the interface. This latter process becomes the main contribution to the surface tension decrease. This work aids our understanding of the most fundamental basis of the physicochemical behavior of concentrated polyelectrolyte-surfactant mixtures which present complex bulk and interfacial interactions with interest in both basic and applied sciences.
△ Less
Submitted 27 January, 2024;
originally announced January 2024.
-
Effect of molecular structure of eco-friendly glycolipid biosurfactants on the adsorption of hair-care conditioning polymers
Authors:
Laura Fernandez-Pena,
Eduardo Guzman,
Fabien Leonforte,
Ana Serrano-Pueyo,
Krzysztof Regulski,
Lucie Tournier-Couturier,
Francisco Ortega,
Ramon G. Rubio,
Gustavo S. Luengo
Abstract:
Pseudo-binary mixtures of different glycolipids, four different rhamnolipids (RL) and an alkyl polyglucoside (APG), with poly(diallyl-dimethylammonium chloride) (PDADMAC) have been studied in relation to their adsorption onto negatively charged surfaces to shed light on the impact of the molecular structure of biodegradable surfactants from natural sources (instead of synthetic surfactant, such as…
▽ More
Pseudo-binary mixtures of different glycolipids, four different rhamnolipids (RL) and an alkyl polyglucoside (APG), with poly(diallyl-dimethylammonium chloride) (PDADMAC) have been studied in relation to their adsorption onto negatively charged surfaces to shed light on the impact of the molecular structure of biodegradable surfactants from natural sources (instead of synthetic surfactant, such as sodium laureth sulfate) on the adsorption of hair-care formulations. For this purpose, the self-assembly of such mixtures in aqueous solution and their adsorption onto negatively charged surfaces mimicking the negative charge of damaged hair fibres have been studied using experiments and self-consistent field (SCF) calculations. The results show that the specific physico-chemical properties of the surfactants (charge, number of sugar rings present in surfactant structure and length of the hydrocarbon length) play a main role in the control of the adsorption process, with the adsorption efficiency and hydration being improved in relation to conventional sulfate-based systems for mixtures of PDADMAC and glycolipids with the shortest alkyl chains. SCF calculations and Energy Dispersive X-Ray Spectroscopy (EDS) analysis on real hair confirmed such observations. The results allow one to assume that the characteristic of the studied surfactants, especially rhamnolipids, conditions positively the adsorption potential of polyelectrolytes in our model systems. This study provides important insights on the mechanisms underlying the performance of more complex but natural and eco-friendly washing formulations.
△ Less
Submitted 26 January, 2024;
originally announced January 2024.
-
Towards understanding the behavior of polyelectrolyte surfactant mixtures at the water vapor interface closer to technologically relevant conditions
Authors:
Sara Llamas,
Laura Fernandez-Pena,
Andrew Akanno,
Eduardo Guzman,
Victor Ortega,
Francisco Ortega,
Aurelio G. Csaky,
Richard A. Campbell,
Ramon G. Rubio
Abstract:
Polyelectrolyte surfactant mixtures and their interactions with fluid interfaces are an important research field due to their use in technological applications. Most of the existing knowledge on these systems is based on models in which the polyelectrolyte concentration is around 50 times lower than that used in commercial formulations. The present work marks a step to close the gap on the underst…
▽ More
Polyelectrolyte surfactant mixtures and their interactions with fluid interfaces are an important research field due to their use in technological applications. Most of the existing knowledge on these systems is based on models in which the polyelectrolyte concentration is around 50 times lower than that used in commercial formulations. The present work marks a step to close the gap on the understanding of their behavior under more practically relevant conditions. The adsorption of concentrated mixtures of poly(diallyldimethylammonium) chloride and sodium N lauroyl N methyltaurate at the water-vapor interface with a crude mixing protocol has been studied by different surface tension techniques, Brewster angle microscopy, neutron reflectometry, and several bulk characterization techniques. Kinetically trapped aggregates formed during mixing influence the interfacial morphology of mixtures produced in the equilibrium one phase region, yet fluctuations in the surface tension isotherm result depending on the tensiometric technique applied. At low bulk surfactant concentrations, the free surfactant concentration is very low, and the interfacial composition matches the trend of the bulk complexes, which is a behavior that has not been observed in studies on more dilute mixtures. Nevertheless, a transition to synergistic co-adsorption of complexes and free surfactant is observed at the higher bulk surfactant concentrations studied. This transition appears to be a special feature of these more concentrated mixtures, which deserves attention in future studies of systems with additional components.
△ Less
Submitted 26 January, 2024;
originally announced January 2024.
-
DeepCERES: A Deep learning method for cerebellar lobule segmentation using ultra-high resolution multimodal MRI
Authors:
Sergio Morell-Ortega,
Marina Ruiz-Perez,
Marien Gadea,
Roberto Vivo-Hernando,
Gregorio Rubio,
Fernando Aparici,
Maria de la Iglesia-Vaya,
Gwenaelle Catheline,
Pierrick Coupé,
José V. Manjón
Abstract:
This paper introduces a novel multimodal and high-resolution human brain cerebellum lobule segmentation method. Unlike current tools that operate at standard resolution ($1 \text{ mm}^{3}$) or using mono-modal data, the proposed method improves cerebellum lobule segmentation through the use of a multimodal and ultra-high resolution ($0.125 \text{ mm}^{3}$) training dataset. To develop the method,…
▽ More
This paper introduces a novel multimodal and high-resolution human brain cerebellum lobule segmentation method. Unlike current tools that operate at standard resolution ($1 \text{ mm}^{3}$) or using mono-modal data, the proposed method improves cerebellum lobule segmentation through the use of a multimodal and ultra-high resolution ($0.125 \text{ mm}^{3}$) training dataset. To develop the method, first, a database of semi-automatically labelled cerebellum lobules was created to train the proposed method with ultra-high resolution T1 and T2 MR images. Then, an ensemble of deep networks has been designed and developed, allowing the proposed method to excel in the complex cerebellum lobule segmentation task, improving precision while being memory efficient. Notably, our approach deviates from the traditional U-Net model by exploring alternative architectures. We have also integrated deep learning with classical machine learning methods incorporating a priori knowledge from multi-atlas segmentation, which improved precision and robustness. Finally, a new online pipeline, named DeepCERES, has been developed to make available the proposed method to the scientific community requiring as input only a single T1 MR image at standard resolution.
△ Less
Submitted 23 January, 2024; v1 submitted 22 January, 2024;
originally announced January 2024.
-
DeepThalamus: A novel deep learning method for automatic segmentation of brain thalamic nuclei from multimodal ultra-high resolution MRI
Authors:
Marina Ruiz-Perez,
Sergio Morell-Ortega,
Marien Gadea,
Roberto Vivo-Hernando,
Gregorio Rubio,
Fernando Aparici,
Mariam de la Iglesia-Vaya,
Thomas Tourdias,
Pierrick Coupé,
José V. Manjón
Abstract:
The implication of the thalamus in multiple neurological pathologies makes it a structure of interest for volumetric analysis. In the present work, we have designed and implemented a multimodal volumetric deep neural network for the segmentation of thalamic nuclei at ultra-high resolution (0.125 mm3). Current tools either operate at standard resolution (1 mm3) or use monomodal data. To achieve the…
▽ More
The implication of the thalamus in multiple neurological pathologies makes it a structure of interest for volumetric analysis. In the present work, we have designed and implemented a multimodal volumetric deep neural network for the segmentation of thalamic nuclei at ultra-high resolution (0.125 mm3). Current tools either operate at standard resolution (1 mm3) or use monomodal data. To achieve the proposed objective, first, a database of semiautomatically segmented thalamic nuclei was created using ultra-high resolution T1, T2 and White Matter nulled (WMn) images. Then, a novel Deep learning based strategy was designed to obtain the automatic segmentations and trained to improve its robustness and accuaracy using a semisupervised approach. The proposed method was compared with a related state-of-the-art method showing competitive results both in terms of segmentation quality and efficiency. To make the proposed method fully available to the scientific community, a full pipeline able to work with monomodal standard resolution T1 images is also proposed.
△ Less
Submitted 15 January, 2024;
originally announced January 2024.
-
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.
-
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.
-
Numerical and experimental study of open-cell foams for the characterization of heat exchangers
Authors:
Aitor Amatriain,
Corrado Gargiulo,
Gonzalo Rubio
Abstract:
A multiscale model of open-cell foams is developed for the characterization of heat exchangers. The model is applicable to a wide range of materials, cell sizes, and porosities. The microscopic geometry is based on a periodic model that is defined by the porosity and the specific surface area of the foam considered. The representative geometrical scales of the model are validated with microscope i…
▽ More
A multiscale model of open-cell foams is developed for the characterization of heat exchangers. The model is applicable to a wide range of materials, cell sizes, and porosities. The microscopic geometry is based on a periodic model that is defined by the porosity and the specific surface area of the foam considered. The representative geometrical scales of the model are validated with microscope images and computed tomography scans. The outputs of the microscopic model are the coefficients of the parabolic pressure loss curve, the thermal conductivity, and the Nusselt number. These values are used as inputs of the macroscopic model that determines the thermal performance of a macroscopic system. The results given by the models are compared with experimental data obtained from the literature, and from an experimental setup built at CERN. It is concluded that the multiscale model provides accurate results in all open-cell foams considered.
△ Less
Submitted 6 June, 2023;
originally announced June 2023.
-
A flux-differencing formulation with Gauss nodes
Authors:
Andrés Mateo-Gabín,
Andrés M. Rueda-Ramírez,
Eusebio Valero,
Gonzalo Rubio
Abstract:
In this work, we propose an extension of telescopic derivative operators for the DGSEM with Gauss nodes, and we prove that this formulation is equivalent to its usual matrix counterpart. Among other possible applications, this allows extending the stabilization methods already developed for Gauss-Lobatto nodes to Gauss nodes, also ensuring properties such as entropy stability while retaining their…
▽ More
In this work, we propose an extension of telescopic derivative operators for the DGSEM with Gauss nodes, and we prove that this formulation is equivalent to its usual matrix counterpart. Among other possible applications, this allows extending the stabilization methods already developed for Gauss-Lobatto nodes to Gauss nodes, also ensuring properties such as entropy stability while retaining their improved accuracy.
△ Less
Submitted 23 June, 2023; v1 submitted 9 November, 2022;
originally announced November 2022.
-
Planetary nebulae with Wolf-Rayet-type central stars -- IV. NGC 1501 and its mixing layer
Authors:
G. Rubio,
J. Toalá,
H. Todt,
L. Sabin,
E. Santamaría,
G. Ramos-Larios,
M. A. Guerrero
Abstract:
Theory predicts that the temperature of the X-ray-emitting gas ($\sim$10$^{6}$ K) detected from planetary nebulae (PNe) is a consequence of mixing or thermal conduction when in contact with the ionized outer rim ($\sim$10$^{4}$ K). Gas at intermediate temperatures ($\sim$10$^{5}$ K) can be used to study the physics of the production of X-ray-emitting gas, via C IV, N V and O VI ions. Here we model…
▽ More
Theory predicts that the temperature of the X-ray-emitting gas ($\sim$10$^{6}$ K) detected from planetary nebulae (PNe) is a consequence of mixing or thermal conduction when in contact with the ionized outer rim ($\sim$10$^{4}$ K). Gas at intermediate temperatures ($\sim$10$^{5}$ K) can be used to study the physics of the production of X-ray-emitting gas, via C IV, N V and O VI ions. Here we model the stellar atmosphere of the CSPN of NGC 1501 to demonstrate that even this hot H-deficient [WO4]-type star cannot produce these emission lines by photoionization. We use the detection of the C IV lines to assess the physical properties of the mixing region in this PNe in comparison with its X-ray-emitting gas, rendering NGC 1501 only the second PNe with such characterization. We extend our predictions to the hottest [WO1] and cooler [WC5] spectral types and demonstrate that most energetic photons are absorbed in the dense winds of [WR] CSPN and highly ionized species can be used to study the physics behind the production of hot bubbles in PNe. We found that the UV observations of NGC 2452, NGC 6751 and NGC 6905 are consistent with the presence mixing layers and hot bubbles, providing excellent candidates for future X-ray observations.
△ Less
Submitted 17 October, 2022;
originally announced October 2022.
-
Non-relativistic gravity theories in four spacetime dimensions
Authors:
Patrick Concha,
Evelyn Rodríguez,
Gustavo Rubio
Abstract:
In this work we present a non-relativistic gravity theory defined in four spacetime dimensions using the MacDowell-Mansouri geometrical formulation. We obtain a Newtonian gravity action which is constructed from the curvature of a Newton-Hooke version of the so-called Newtonian algebra. We show that the non-relativistic gravity theory presented here contains the Poisson equation in presence of a c…
▽ More
In this work we present a non-relativistic gravity theory defined in four spacetime dimensions using the MacDowell-Mansouri geometrical formulation. We obtain a Newtonian gravity action which is constructed from the curvature of a Newton-Hooke version of the so-called Newtonian algebra. We show that the non-relativistic gravity theory presented here contains the Poisson equation in presence of a cosmological constant. Moreover we make contact with the Modified Newtonian Dynamics (MOND) approach for gravity by considering a particular ansatz for a given gauge field. We extend our results to a generalized non-relativistic MacDowell-Mansouri gravity theory by considering a generalized Newton-Hooke algebra.
△ Less
Submitted 8 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.
-
Gauging quantum states with non-anomalous matrix product operator symmetries
Authors:
José Garre Rubio,
Ilya Kull
Abstract:
Gauging a global symmetry of a system amounts to introducing new degrees of freedom whose transformation rule makes the overall system observe a local symmetry. In quantum systems there can be obstructions to gauging a global symmetry. When this happens the symmetry is dubbed anomalous. Such obstructions are related to the fact that the global symmetry cannot be written as a tensor product of loca…
▽ More
Gauging a global symmetry of a system amounts to introducing new degrees of freedom whose transformation rule makes the overall system observe a local symmetry. In quantum systems there can be obstructions to gauging a global symmetry. When this happens the symmetry is dubbed anomalous. Such obstructions are related to the fact that the global symmetry cannot be written as a tensor product of local operators. In this manuscript we study non-local symmetries that have an additional structure: they take the form of a matrix product operator (MPO). We exploit the tensor network structure of the MPOs to construct local operators from them satisfying the same group relations, that is, we are able to localize even anomalous MPOs. For non-anomalous MPOs, we use these local operators to explicitly gauge the MPO symmetry of a one-dimensional quantum state obtaining non-trivial gauged states. We show that our gauging procedure satisfies all the desired properties as the standard on-site case does. We also show how this procedure is naturally represented in matrix product states protected by MPO symmetries. In the case of anomalous MPOs, we shed light on the obstructions to gauging these symmetries.
△ Less
Submitted 15 February, 2023; v1 submitted 15 September, 2022;
originally announced September 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.
-
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.
-
Three-dimensional Newtonian gravity with cosmological constant and torsion
Authors:
Patrick Concha,
Evelyn Rodríguez,
Gustavo Rubio,
Paola Yañez
Abstract:
In this paper we present an alternative cosmological extension of the three-dimensional extended Newtonian Chern-Simons gravity by switching on the torsion. The theory is obtained as a non-relativistic limit of an enhancement and $U(1)$-enlargement of the so-called teleparallel algebra and can be seen as the teleparallel analogue of the Newtonian gravity theory. The infinite-dimensional extension…
▽ More
In this paper we present an alternative cosmological extension of the three-dimensional extended Newtonian Chern-Simons gravity by switching on the torsion. The theory is obtained as a non-relativistic limit of an enhancement and $U(1)$-enlargement of the so-called teleparallel algebra and can be seen as the teleparallel analogue of the Newtonian gravity theory. The infinite-dimensional extension of our result is also explored through the Lie algebra expansion method. An infinite-dimensional torsional Galilean gravity model is presented which in the vanishing cosmological constant limit reproduces the infinite-dimensional extension of the Galilean gravity theory.
△ Less
Submitted 25 April, 2022;
originally announced April 2022.
-
Common envelope evolution in born-again planetary nebulae -- Shaping the H-deficient ejecta of A30
Authors:
J. B. Rodríguez-González,
E. Santamaría,
J. A. Toalá,
M. A. Guerrero,
B. Montoro-Molina,
G. Rubio,
D. Tafoya,
Y. -H. Chu,
G. Ramos-Larios,
L. Sabin
Abstract:
Born-again planetary nebulae (PNe) are extremely rare cases in the evolution of solar-like stars. It is commonly accepted that their central stars (CSPN) experienced a very late thermal pulse (VLTP), ejecting H-deficient material inside the evolved H-rich PN. Given the short duration of this event and the fast subsequent evolution of the CSPN, details of the mass ejection are unknown. We present t…
▽ More
Born-again planetary nebulae (PNe) are extremely rare cases in the evolution of solar-like stars. It is commonly accepted that their central stars (CSPN) experienced a very late thermal pulse (VLTP), ejecting H-deficient material inside the evolved H-rich PN. Given the short duration of this event and the fast subsequent evolution of the CSPN, details of the mass ejection are unknown. We present the first morpho-kinematic model of the H-deficient material surrounding a born-again PN, namely A30. New San Pedro Mártir observations with the Manchester Echelle Spectrograph were recently obtained to map the inner region of A30 which are interpreted by means of the software SHAPE in conjunction with HST WFC3 images. The SHAPE morpho-kinematic model that best reproduces the observations is composed by a disrupted disk tilted $37^\circ$ with respect to the line of sight and a pair of orthogonal opposite bipolar ejections. We confirm previous suggestions that the structures closer to the CSPN present the highest expansion velocities, that is, the disrupted disk expands faster than the farther bipolar features. We propose that the current physical structure and abundance discrepancy of the H-deficient clumps around the CSPN of A30 can be explained by a common envelope phase following the VLTP event. Our proposed scenario is also compared with other known born-again PNe (A58, A78, HuBi1 and the Sakurai's Object).
△ Less
Submitted 20 June, 2022; v1 submitted 18 April, 2022;
originally announced April 2022.
-
3D mapping of the Wolf-Rayet nebula M1-67: clues for post-common envelope evolution in massive stars
Authors:
S. Zavala,
J. A. Toalá,
E. Santamarí,
G. Ramos-Larios,
L. Sabin,
J. A. Quino-Mendoza,
G. Rubio,
M. A. Guerrero
Abstract:
We present a 3D mapping of the Wolf-Rayet (WR) nebula M1-67 around WR124. We obtained high-resolution San Pedro Mártir (SPM) Manchester Echelle Spectrograph (MES) observations along 17 long-slit positions covering all morphological features in M1-67. We are able to unveil the true morphology of M1-67 and its kinematics by interpreting the SPM MES observations by means of the 3D modelling tool for…
▽ More
We present a 3D mapping of the Wolf-Rayet (WR) nebula M1-67 around WR124. We obtained high-resolution San Pedro Mártir (SPM) Manchester Echelle Spectrograph (MES) observations along 17 long-slit positions covering all morphological features in M1-67. We are able to unveil the true morphology of M1-67 and its kinematics by interpreting the SPM MES observations by means of the 3D modelling tool for Astrophysics SHAPE. Our SHAPE model that best reproduces the SPM MES data includes three concentric bipolar structures composed by a hollow ellipsoidal structure and a torus. In addition, the model requires the presence of expanding jets and broken blisters in order to reproduce specific spectral features. Our results are consistent with the idea that M1-67 and its progenitor star WR124 have formed through a common envelope scenario that occurred 11.8$^{+4.6}_{-0.8}$ kyr ago. Our bipolar model strongly questions previous suggestions of the presence of a bow shock structure surrounding M1-67. We interpret that the bright structures detected in the spectra extracted from the central regions are produced by wind compression at the receding region of the innermost structure in M1-67. Furthermore, WR124 is moving through a low-density region above the Galactic plane that has negligibly affected the formation history of M1-67.
△ Less
Submitted 16 April, 2022;
originally announced April 2022.
-
Rotating Micro-Spheres for Adsorption Monitoring at a Fluid Interface
Authors:
J. Martín-Roca,
M. Jiménez,
F. Ortega,
C. Calero,
Ch. Valeriani,
R. G. Rubio,
F. Martínez-Pedrero
Abstract:
Hypothesis: A broad range of phenomena, such as emulsification and emulsion stability, foam formation or liquid evaporation, are closely related to the dynamics of adsorbing colloidal particles. Elucidation of the mechanisms implied is key to a correct design of many different types of materials.
Experiments: Microspheres forced to rotate near a fluid interface exhibit a roto-translational hydro…
▽ More
Hypothesis: A broad range of phenomena, such as emulsification and emulsion stability, foam formation or liquid evaporation, are closely related to the dynamics of adsorbing colloidal particles. Elucidation of the mechanisms implied is key to a correct design of many different types of materials.
Experiments: Microspheres forced to rotate near a fluid interface exhibit a roto-translational hydrodynamic mechanism that is hindered by capillary torques as soon as the particles protrude the interface. Under these conditions, the time evolution in the ratio of moving spheres provides a direct description of the adsorption kinetics, while microscopy monitoring of particle acceleration\deceleration informs about the adsorption\desorption dynamics. In this work, the proposed strategy is applied at an air/water interface loaded with spherical magnetic particles negatively charged, forced to rotate by the action of a rotating magnetic field.
Findings: The proposed method enables the adsorption/desorption dynamics to be followed during the earliest phase of the process, when desorption of a small fraction of particles is detected, as well as to estimate approximated values of the adsorption/desorption constants. The results obtained show that the addition of a monovalent salt or a cationic (anionic) surfactant promotes (inhibits) both adsorption and formation of permanent bonds between particles.
△ Less
Submitted 19 February, 2022;
originally announced February 2022.
-
A multivariate extension of the Misspecification-Resistant Information Criterion
Authors:
Gery Andrés Díaz Rubio,
Simone Giannerini,
Greta Goracci
Abstract:
The Misspecification-Resistant Information Criterion (MRIC) proposed in [H.-L. Hsu, C.-K. Ing, H. Tong: On model selection from a finite family of possibly misspecified time series models. The Annals of Statistics. 47 (2), 1061--1087 (2019)] is a model selection criterion for univariate parametric time series that enjoys both the property of consistency and asymptotic efficiency. In this article w…
▽ More
The Misspecification-Resistant Information Criterion (MRIC) proposed in [H.-L. Hsu, C.-K. Ing, H. Tong: On model selection from a finite family of possibly misspecified time series models. The Annals of Statistics. 47 (2), 1061--1087 (2019)] is a model selection criterion for univariate parametric time series that enjoys both the property of consistency and asymptotic efficiency. In this article we extend the MRIC to the case where the response is a multivariate time series and the predictor is univariate. The extension requires novel derivations based upon random matrix theory. We obtain an asymptotic expression for the mean squared prediction error matrix, the vectorial MRIC and prove the consistency of its method-of-moments estimator. Moreover, we prove its asymptotic efficiency. Finally, we show with an example that, in presence of misspecification, the vectorial MRIC identifies the best predictive model whereas traditional information criteria like AIC or BIC fail to achieve the task.
△ Less
Submitted 18 February, 2022;
originally announced February 2022.
-
vol2Brain: A new online Pipeline for whole Brain MRI analysis
Authors:
Jose V. Manjon,
Jose E. Romero,
Roberto Vivo-Hernando,
Gregorio Rubio,
Fernando Aparici,
Mariam de la Iglesia-Vaya,
Pierrick Coupe
Abstract:
Automatic and reliable quantitative tools for MR brain image analysis are a very valuable resources for both clinical and research environments. In the last years, this field has experienced many advances with successful techniques based on label fusion and more recently deep learning. However, few of them have been specifically designed to provide a dense anatomical labelling at multiscale level…
▽ More
Automatic and reliable quantitative tools for MR brain image analysis are a very valuable resources for both clinical and research environments. In the last years, this field has experienced many advances with successful techniques based on label fusion and more recently deep learning. However, few of them have been specifically designed to provide a dense anatomical labelling at multiscale level and to deal with brain anatomical alterations such as white matter lesions. In this work, we present a fully automatic pipeline (vol2Brain) for whole brain segmentation and analysis which densely labels (N>100) the brain while being robust to the presence of white matter lesions. This new pipeline is an evolution of our previous volBrain pipeline that extends significantly the number of regions that can be analyzed. Our proposed method is based on a fast multiscale multi-atlas label fusion technology with systematic error correction able to provide accurate volumetric information in few minutes. We have deployed our new pipeline within our platform volBrain (www.volbrain.upv.es) which has been already demonstrated to be an efficient and effective manner to share our technology with users worldwide
△ Less
Submitted 8 February, 2022;
originally announced February 2022.
-
Planetary nebulae with Wolf-Rayet-type central stars -- III. A detailed view of NGC 6905 and its central star
Authors:
V. M. A. Gómez-González,
G. Rubio,
J. A. Toalá,
M. A. Guerrero,
L. Sabin,
H. Todt,
V. Gómez-Llanos,
G. Ramos-Larios,
Y. D. Mayya
Abstract:
We present a multi-wavelength characterisation of the planetary nebula (PN) NGC 6905 and its [Wolf-Rayet]-type ([WR]) central star (CSPN) HD 193949. Our Nordic Optical Telescope (NOT) Alhambra Faint Object Spectrograph and Camera (ALFOSC) spectra and images unveil in unprecedented detail the high-ionization structure of NGC 6905. The high-quality spectra of HD 193949 allowed us to detect more than…
▽ More
We present a multi-wavelength characterisation of the planetary nebula (PN) NGC 6905 and its [Wolf-Rayet]-type ([WR]) central star (CSPN) HD 193949. Our Nordic Optical Telescope (NOT) Alhambra Faint Object Spectrograph and Camera (ALFOSC) spectra and images unveil in unprecedented detail the high-ionization structure of NGC 6905. The high-quality spectra of HD 193949 allowed us to detect more than 20 WR features including the characteristic O-bump, blue bump and red bump, which suggests a spectral type no later than a [WO2]-subtype. Moreover we detect the Ne VII and Ne VIII broad emission lines, rendering HD 193949 yet another CSPN with $T_\mathrm{eff}\lesssim150$ kK exhibiting such stellar emission lines. We studied the physical properties ($T_\mathrm{e}$ and $n_\mathrm{e}$) and chemical abundances of different regions within NGC 6905 including its low-ionization clumps; abundances are found to be homogeneous. We used the PoWR stellar atmosphere code to model the spectrum of HD 193949, which is afterwards used in a photoionization model performed with Cloudy that reproduces the nebular and dust properties for a total mass in the 0.31-0.47 M$_{\odot}$ range and a mass of C-rich dust of $\sim$2 $\times10^{-3}$ M$_{\odot}$. Adopting a current stellar mass of 0.6 M$_{\odot}$, our model suggests an initial mass $\sim$1 M$_\odot$ for HD 193949, consistent with the observations.
△ Less
Submitted 18 October, 2021;
originally announced October 2021.
-
X-ray observations of the nova shell IPHASX J210204.7+471015
Authors:
J. A. Toalá,
G. Rubio,
E. Santamaría,
M. A. Guerrero,
S. Estrada-Dorado,
G. Ramos-Larios,
L. Sabin
Abstract:
We present the analysis of XMM-Newton European Photon Imaging Camera (EPIC) observations of the nova shell IPHASX J210204.7$+$471015. We detect X-ray emission from the progenitor binary star with properties that resemble those of underluminous intermediate polars such as DQ Her: an X-ray-emitting plasma with temperature of $T_\mathrm{X}=(6.4\pm3.1)\times10^{6}$ K, a non-thermal X-ray component, an…
▽ More
We present the analysis of XMM-Newton European Photon Imaging Camera (EPIC) observations of the nova shell IPHASX J210204.7$+$471015. We detect X-ray emission from the progenitor binary star with properties that resemble those of underluminous intermediate polars such as DQ Her: an X-ray-emitting plasma with temperature of $T_\mathrm{X}=(6.4\pm3.1)\times10^{6}$ K, a non-thermal X-ray component, and an estimated X-ray luminosity of $L_\mathrm{X}=10^{30}$ erg s$^{-1}$. Time series analyses unveil the presence of two periods, the dominant with a period of $2.9\pm0.2$ hr, which might be attributed to the spin of the white dwarf, and a secondary of $4.5\pm0.6$ hr that is in line with the orbital period of the binary system derived from optical observations. We do not detect extended X-ray emission as in other nova shells probably due to its relatively old age (130-170 yr) or to its asymmetric disrupted morphology which is suggestive of explosion scenarios different to the symmetric ones assumed in available numerical simulations of nova explosions.
△ Less
Submitted 27 January, 2021;
originally announced January 2021.
-
High--order discontinuous Galerkin approximation for a three--phase incompressible Navier--Stokes/Cahn--Hilliard model
Authors:
Juan Manzanero,
Carlos Redondo,
Miguel Chávez--Módena,
Gonzalo Rubio,
Eusebio Valero,
Susana Gómez--Álvarez,
Ángel Rivero--Jiménez
Abstract:
In this work we introduce the development of a three--phase incompressible Navier--Stokes/Cahn--Hilliard numerical method to simulate three--phase flows, present in many industrial operations. The numerical method is then applied to successfully solve oil transport problems, such as those found in the oil and gas industry. The three--phase model adopted in this work is a Cahn--Hilliard diffuse int…
▽ More
In this work we introduce the development of a three--phase incompressible Navier--Stokes/Cahn--Hilliard numerical method to simulate three--phase flows, present in many industrial operations. The numerical method is then applied to successfully solve oil transport problems, such as those found in the oil and gas industry. The three--phase model adopted in this work is a Cahn--Hilliard diffuse interface model, which was derived by Boyer and Lapuerta et al. 2006. The Cahn--Hilliard model is coupled to the entropy--stable incompressible Navier--Stokes equations model derived by Manzanero et al. 2019. The spatial discretization uses a high--order discontinuous Galerkin spectral element method which yields highly accurate results in arbitrary geometries, while an implicit--explicit (IMEX) method is adopted as temporal discretization. The developed numerical tool is tested for two and three dimensional problems, including a convergence study, a two--dimensional jet, a three--dimensional annular flow, and realistic geometries like T--shaped pipe intersections.
△ Less
Submitted 14 December, 2020;
originally announced December 2020.
-
Unveiling the stellar origin of the Wolf-Rayet nebula NGC6888 through infrared observations
Authors:
G. Rubio,
J. A. Toalá,
P. Jiménez-Hernández,
G. Ramos-Larios,
M. A. Guerrero,
V. M. A. Gómez-González,
E. Santamaría,
J. A. Quino-Mendoza
Abstract:
We present a comprehensive infrared (IR) study of the iconic Wolf-Rayet (WR) wind-blown bubble NGC6888 around WR136. We use Wide-field Infrared Survey Explorer (WISE), Spitzer IRAC and MIPS and Herschel PACS IR images to produce a sharp view of the distribution of dust around WR136. We complement these IR photometric observations with Spitzer IRS spectra in the 5-38 $μ$m wavelength range. The unpr…
▽ More
We present a comprehensive infrared (IR) study of the iconic Wolf-Rayet (WR) wind-blown bubble NGC6888 around WR136. We use Wide-field Infrared Survey Explorer (WISE), Spitzer IRAC and MIPS and Herschel PACS IR images to produce a sharp view of the distribution of dust around WR136. We complement these IR photometric observations with Spitzer IRS spectra in the 5-38 $μ$m wavelength range. The unprecedented high-resolution IR images allowed us to produce a clean spectral energy distribution, free of contamination from material along the line of sight, to model the properties of the dust in NGC6888. We use the spectral synthesis code Cloudy to produce a model for NGC6888 that consistently reproduces its optical and IR properties. Our best model requires a double distribution with the inner shell composed only of gas, whilst the outer shell requires a mix of gas and dust. The dust consists of two populations of grain sizes, one with small sized grains $a_\mathrm{small}$=[0.002-0.008] $μ$m and another one with large sized grains $a_\mathrm{big}$=[0.05-0.5] $μ$m. The population of big grains is similar to that reported for other red supergiants stars and dominates the total dust mass, which leads us to suggest that the current mass of NGC6888 is purely due to material ejected from WR136, with a negligible contribution of swept up interstellar medium. The total mass of this model is 25.5$^{+4.7}_{-2.8}$ M$_{\odot}$, a dust mass of $M_\mathrm{dust}=$0.14$^{+0.03}_{-0.01}$ M$_{\odot}$, for a dust-to-gas ratio of $5.6\times10^{-3}$. Accordingly, we suggest that the initial stellar mass of WR136 was $\lesssim$50 M$_{\odot}$, consistent with current single stellar evolution models.
△ Less
Submitted 11 September, 2020;
originally announced September 2020.
-
Three-dimensional Maxwellian Extended Newtonian gravity and flat limit
Authors:
Patrick Concha,
Lucrezia Ravera,
Evelyn Rodríguez,
Gustavo Rubio
Abstract:
In the present work we find novel Newtonian gravity models in three spacetime dimensions. We first present a Maxwellian version of the extended Newtonian gravity, which is obtained as the non-relativistic limit of a particular $U(1)$-enlargement of an enhanced Maxwell Chern-Simons gravity. We show that the extended Newtonian gravity appears as a particular sub-case. Then, the introduction of a cos…
▽ More
In the present work we find novel Newtonian gravity models in three spacetime dimensions. We first present a Maxwellian version of the extended Newtonian gravity, which is obtained as the non-relativistic limit of a particular $U(1)$-enlargement of an enhanced Maxwell Chern-Simons gravity. We show that the extended Newtonian gravity appears as a particular sub-case. Then, the introduction of a cosmological constant to the Maxwellian extended Newtonian theory is also explored. To this purpose, we consider the non-relativistic limit of an enlarged symmetry. An alternative method to obtain our results is presented by applying the semigroup expansion method to the enhanced Nappi-Witten algebra. The advantages of considering the Lie algebra expansion procedure is also discussed.
△ Less
Submitted 25 September, 2020; v1 submitted 23 June, 2020;
originally announced June 2020.
-
The benefits of a Bayesian analysis for the characterization of magnetic nanoparticles
Authors:
Mathias Bersweiler,
Helena Gavilan Rubio,
Dirk Honecker,
Andreas Michels,
Philipp Bender
Abstract:
Magnetic nanoparticles offer a unique potential for various biomedical applications, but prior to commercial usage a standardized characterization of their structural and magnetic properties is required. For a thorough characterization, the combination of conventional magnetometry and advanced scattering techniques has shown great potential. In the present work, we characterize a powder sample of…
▽ More
Magnetic nanoparticles offer a unique potential for various biomedical applications, but prior to commercial usage a standardized characterization of their structural and magnetic properties is required. For a thorough characterization, the combination of conventional magnetometry and advanced scattering techniques has shown great potential. In the present work, we characterize a powder sample of high-quality iron oxide nanoparticles that are surrounded with a homogeneous thick silica shell by DC magnetometry and magnetic small-angle neutron scattering (SANS). To retrieve the particle parameters such as their size distribution and saturation magnetization from the data, we apply standard model fits of individual data sets as well as global fits of multiple curves, including a combination of the magnetometry and SANS measurements. We show that by combining a standard least-squares fit with a subsequent Bayesian approach for the data refinement, the probability distributions of the model parameters and their cross correlations can be readily extracted, which enables a direct visual feedback regarding the quality of the fit. This prevents an overfitting of data in case of highly correlated parameters and renders the Bayesian method as an ideal component for a standardized data analysis of magnetic nanoparticle samples.
△ Less
Submitted 17 August, 2020; v1 submitted 19 May, 2020;
originally announced May 2020.
-
Spin orbit field in a physically defined p type MOS silicon double quantum dot
Authors:
Marian Marx,
Jun Yoneda,
Ángel Gutiérrez Rubio,
Peter Stano,
Tomohiro Otsuka,
Kenta Takeda,
Sen Li,
Yu Yamaoka,
Takashi Nakajima,
Akito Noiri,
Daniel Loss,
Tetsuo Kodera,
Seigo Tarucha
Abstract:
We experimentally and theoretically investigate the spin orbit (SO) field in a physically defined, p type metal oxide semiconductor double quantum dot in silicon. We measure the magnetic field dependence of the leakage current through the double dot in the Pauli spin blockade. A finite magnetic field lifts the blockade, with the lifting least effective when the external and SO fields are parallel.…
▽ More
We experimentally and theoretically investigate the spin orbit (SO) field in a physically defined, p type metal oxide semiconductor double quantum dot in silicon. We measure the magnetic field dependence of the leakage current through the double dot in the Pauli spin blockade. A finite magnetic field lifts the blockade, with the lifting least effective when the external and SO fields are parallel. In this way, we find that the spin flip of a tunneling hole is due to a SO field pointing perpendicular to the double dot axis and almost fully out of the quantum well plane. We augment the measurements by a derivation of SO terms using group symmetric representations theory. It predicts that without in plane electric fields (a quantum well case), the SO field would be mostly within the plane, dominated by a sum of a Rashba and a Dresselhaus like term. We, therefore, interpret the observed SO field as originated in the electric fields with substantial in plane components.
△ Less
Submitted 17 March, 2020; v1 submitted 16 March, 2020;
originally announced March 2020.
-
Collective Transport of Magnetic Microparticles at a Fluid Interface through Dynamic Self-Assembled Lattices
Authors:
Fernando Martínez-Pedrero,
Francisco Ortega,
Ramón G. Rubio,
Carles Calero
Abstract:
The transport of motile entities across modulated energy landscapes plays an important role in a range of phenomena in biology, colloidal science and solid-state physics. Here, an easily implementable strategy that allows for the collective and monitored transport of microparticles at fluid-fluid interfaces is introduced. Adsorbed magnetic microparticles are carried on time-dependent magnetic pote…
▽ More
The transport of motile entities across modulated energy landscapes plays an important role in a range of phenomena in biology, colloidal science and solid-state physics. Here, an easily implementable strategy that allows for the collective and monitored transport of microparticles at fluid-fluid interfaces is introduced. Adsorbed magnetic microparticles are carried on time-dependent magnetic potentials, generated by a dynamic self-assembled lattice of different-sized magnetic particles. In such binary systems, the sudden reorientation of the applied field triggers the rapid exchange between attractive and repulsive configurations, enabling for the ballistic transfer of the carriers through the lattice. As the number of motile entities increases, the induced current increases, before reaching a maximum, while the loaded interface gradually displays bidirectional transport. The described methodology can be tuned through the applied field and exploited for the monitored guidance of adsorbed molecules on liquid surfaces, the segregation of colloidal mixtures, the induced motion of defects in photonic crystals or the design of new self-assembled microrobots.
△ Less
Submitted 10 March, 2020;
originally announced March 2020.
-
Angular Expansion of Nova Shells
Authors:
E. Santamaría,
M. A. Guerrero,
G. Ramos-Larios,
J. A. Toalá,
L. Sabin,
G. Rubio,
J. A. Quino-Mendoza
Abstract:
Nova shells can provide us with important information on their distance, their interactions with the circumstellar and interstellar media, and the evolution in morphology of the ejecta. We have obtained narrow-band images of a sample of five nova shells, namely DQHer, FHSer, TAur, V476Cyg, and V533Her, with ages in the range from 50 to 130 years. These images have been compared with suitable avail…
▽ More
Nova shells can provide us with important information on their distance, their interactions with the circumstellar and interstellar media, and the evolution in morphology of the ejecta. We have obtained narrow-band images of a sample of five nova shells, namely DQHer, FHSer, TAur, V476Cyg, and V533Her, with ages in the range from 50 to 130 years. These images have been compared with suitable available archival images to derive their angular expansion rates. We find that all the nova shells in our sample are still in the free expansion phase, which can be expected, as the mass of the ejecta is 7-45 times larger than the mass of the swept-up circumstellar medium. The nova shells will keep expanding freely for time periods up to a few hundred years, reducing their time dispersal into the interstellar medium.
△ Less
Submitted 13 March, 2020; v1 submitted 16 February, 2020;
originally announced February 2020.
-
Cosmology from Newton-Chern-Simons gravity
Authors:
S. Lepe,
G. Rubio,
P. Salgado
Abstract:
We study a five-dimensional non-relativistic gravity theory whose action is composed of a gravitational sector and a sector of matter where the gravitational sector is given by the so called Newton--Chern--Simons gravity and where the matter sector is described by a perfect fluid. At time to do cosmology, the obtained field equations shows a close analogy with the projectable version of the Hořava…
▽ More
We study a five-dimensional non-relativistic gravity theory whose action is composed of a gravitational sector and a sector of matter where the gravitational sector is given by the so called Newton--Chern--Simons gravity and where the matter sector is described by a perfect fluid. At time to do cosmology, the obtained field equations shows a close analogy with the projectable version of the Hořava--Lifshitz theory in (3+1)-dimensions. Solutions and their asymptotic limits are found. In particular a phantom solution with a future singularity reminiscent of a Litlle Big Rip future singularity is obtained.
△ Less
Submitted 13 November, 2019;
originally announced November 2019.
-
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.
-
The GAPS Programme with HARPS-N at TNG XIX. Atmospheric Rossiter-McLaughlin effect and improved parameters of KELT-9b
Authors:
F. Borsa,
M. Rainer,
A. S. Bonomo,
D. Barbato,
L. Fossati,
L. Malavolta,
V. Nascimbeni,
A. F. Lanza,
M. Esposito,
L. Affer,
G. Andreuzzi,
S. Benatti,
K. Biazzo,
A. Bignamini,
M. Brogi,
I. Carleo,
R. Claudi,
R. Cosentino,
E. Covino,
M. Damasso,
S. Desidera,
A. Garrido Rubio,
P. Giacobbe,
E. González-Álvarez,
A. Harutyunyan
, et al. (18 additional authors not shown)
Abstract:
In the framework of the GAPS project, we observed the planet-hosting star KELT-9 (A-type star, VsinI$\sim$110 km/s) with the HARPS-N spectrograph at the TNG. In this work we analyse the spectra and the extracted radial velocities (RVs), to constrain the physical parameters of the system and to detect the planetary atmosphere of KELT-9b. We extracted from the high-resolution optical spectra the mea…
▽ More
In the framework of the GAPS project, we observed the planet-hosting star KELT-9 (A-type star, VsinI$\sim$110 km/s) with the HARPS-N spectrograph at the TNG. In this work we analyse the spectra and the extracted radial velocities (RVs), to constrain the physical parameters of the system and to detect the planetary atmosphere of KELT-9b. We extracted from the high-resolution optical spectra the mean stellar line profiles with an analysis based on the Least Square Deconvolution technique. Then, we computed the stellar RVs with a method optimized for fast rotators, by fitting the mean stellar line profile with a purely rotational profile instead of using a Gaussian function. The new spectra and analysis led us to update the orbital and physical parameters of the system, improving in particular the value of the planetary mass to $M_{\rm p}=2.88\pm0.35\,M_{\rm Jup}$. We discovered an anomalous in-transit RV deviation from the theoretical Rossiter-McLaughlin effect solution, calculated from the projected spin-orbit angle $λ=-85.78\pm0.46$ degrees measured with Doppler tomography. We prove that this deviation is caused by the planetary atmosphere of KELT-9b, thus we name this effect Atmospheric Rossiter-McLaughlin effect. By analysing the magnitude of the RV anomaly, we obtained information on the extension of the planetary atmosphere as weighted by the model used to retrieve the stellar mean line profiles, which is up to $1.22\pm0.02\,R_{\rm p}$. The Atmospheric Rossiter-McLaughlin effect will be observable for other exoplanets whose atmosphere has non-negligible correlation with the stellar mask used to retrieve the RVs, in particular ultra-hot Jupiters with iron in their atmosphere. The duration and amplitude of the effect will depend not only on the extension of the atmosphere, but also on the in-transit planetary RVs and on the projected rotational velocity of the parent star.
△ Less
Submitted 26 September, 2019; v1 submitted 23 July, 2019;
originally announced July 2019.
-
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.
-
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.
-
HADES RV program with HARPS-N at TNG. IX. A super-Earth around the M dwarf Gl686
Authors:
L. Affer,
M. Damasso,
G. Micela,
E. Poretti,
G. Scandariato,
J. Maldonado,
A. F. Lanza,
E. Covino,
A. Garrido Rubio,
J. I. Gonzalez Hernandez,
R. Gratton,
G. Leto,
A. Maggio,
M. Perger,
A. Sozzetti,
A. Suarez Mascareno,
A. S. Bonomo,
F. Borsa,
R. Claudi,
R. Cosentino,
S. Desidera,
P. Giacobbe,
E. Molinari,
M. Pedani,
M. Pinamonti
, et al. (3 additional authors not shown)
Abstract:
The HArps-n red Dwarf Exoplanet Survey is providing a major contribution to the widening of the current statistics of low-mass planets, through the in-depth analysis of precise radial velocity measurements in a narrow range of spectral sub-types. As part of that program, we obtained radial velocity measurements of Gl 686, an M1 dwarf at d = 8.2 pc. The analysis of data obtained within an intensive…
▽ More
The HArps-n red Dwarf Exoplanet Survey is providing a major contribution to the widening of the current statistics of low-mass planets, through the in-depth analysis of precise radial velocity measurements in a narrow range of spectral sub-types. As part of that program, we obtained radial velocity measurements of Gl 686, an M1 dwarf at d = 8.2 pc. The analysis of data obtained within an intensive observing campaign demonstrates that the excess dispersion is due to a coherent signal, with a period of 15.53 d. Almost simultaneous photometric observations were carried out within the APACHE and EXORAP programs to characterize the stellar activity and to distinguish periodic variations related to activity from signals due to the presence of planetary companions, complemented also with ASAS photometric data. We took advantage of the available radial velocity measurements for this target from other observing campaigns. The analysis of the radial velocity composite time series from the HIRES, HARPS and HARPS-N spectrographs, consisting of 198 measurements taken over 20 years, enabled us to address the nature of periodic signals and also to characterize stellar physical parameters (mass, temperature, and rotation). We report the discovery of a super-Earth orbiting at a distance of 0.092 AU from the host star Gl 686. Gl 686 b has a minimum mass of 7.1 +/- 0.9 MEarth and an orbital period of 15.532 +/- 0.002 d. The analysis of the activity indexes, correlated noise through a Gaussian process framework and photometry, provides an estimate of the stellar rotation period at 37 d, and highlights the variability of the spot configuration during the long timespan covering 20 yrs. The observed periodicities around 2000 d likely point to the existence of an activity cycle.
△ Less
Submitted 31 January, 2019; v1 submitted 16 January, 2019;
originally announced January 2019.
-
Modified newtonian dynamics and non-relativistic ChSAS gravity
Authors:
G. Rubio,
P. Salgado
Abstract:
In the context of the non-relativistic theories, a generalization of the Chern--Weil-theorem allows us to show that extended Chern--Simons actions for gravity in d=4 invariant under some specific non-relativistic groups lead to modified Poisson equations. In some particular cases, these modified equations have the form of the so-called MOND approach to gravity. The modifications could be understoo…
▽ More
In the context of the non-relativistic theories, a generalization of the Chern--Weil-theorem allows us to show that extended Chern--Simons actions for gravity in d=4 invariant under some specific non-relativistic groups lead to modified Poisson equations. In some particular cases, these modified equations have the form of the so-called MOND approach to gravity. The modifications could be understood as due to the effects of dark matter. This result could leads us to think that dark matter can be interpreted as a non-relativistic limit of dark energy.
△ Less
Submitted 5 November, 2018;
originally announced November 2018.