BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Pedro Costa

Pedro Costa appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
hep-phphysics.flu-dynnucl-thComputer Visioneess.IVphysics.comp-phMachine Learningastro-ph.SRComputational Engineering, Finance, and Sciencephysics.ed-ph

Catalog footprint

What is connected

48works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

48 published item(s)

preprint2025arXiv

Mean velocity profile in stably stratified turbulent channel flow

The Monin-Obukhov Similarity Theory (MOST) is a cornerstone of atmospheric science for describing turbulence in stable boundary layers. Extending MOST to stably stratified turbulent channel flows, however, is non-trivial due to confinement by solid walls and the much smaller turbulent length scales involved. In this study, we investigate the applicability of MOST in closed channels and identify where and to what extent the theory remains valid. A key finding is that the ratio of the half-channel height to the Obukhov length serves as a governing parameter for identifying distinct flow regions and determining the scaling of the mean velocity within them. Hence, we propose a closure relation to estimate this ratio directly from the governing input parameters: friction Reynolds and friction Richardson numbers ($Re_τ$ and $Ri_τ$). The framework is tested against a series of direct numerical simulations (DNS) across a range of $Re_τ$ and $Ri_τ$. The reconstructed velocity profiles enable accurate prediction of the skin friction coefficient crucial for quantifying pressure losses in stratified flows in engineering applications.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A FFT-accelerated multi-block finite-difference solver for massively parallel simulations of incompressible flows

We present a multi-block finite-difference solver for massively parallel Direct Numerical Simulations (DNS) of incompressible flows. The algorithm combines the versatility of a multi-block solver with the method of eigenfunctions expansions, to speedup the solution of the pressure Poisson equation. This is achieved by employing FFT-based transforms along one homogeneous direction, which effectively reduce the problem complexity at a low cost. These FFT-based expansions are implemented in a framework that unifies all valid combinations of boundary conditions for this type of method. Subsequently, a geometric multigrid solver is employed to solve the reduced Poisson equation in a multi-block geometry. Particular care was taken here, to guarantee the parallel performance of the multigrid solver when solving the reduced linear systems equations. We have validated the overall numerical algorithm and assessed its performance in a massively parallel setting. The results show that 2- to 8-fold reductions in computational cost may be easily achieved when exploiting FFT acceleration for the solution of the Poisson equation. The solver, SNaC, has been made freely available and open-source under the terms of an MIT license.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

FluTAS: A GPU-accelerated finite difference code for multiphase flows

We present the Fluid Transport Accelerated Solver, FluTAS, a scalable GPU code for multiphase flows with thermal effects. The code solves the incompressible Navier-Stokes equation for two-fluid systems, with a direct FFT-based Poisson solver for the pressure equation. The interface between the two fluids is represented with the Volume of Fluid (VoF) method, which is mass conserving and well suited for complex flows thanks to its capacity of handling topological changes. The energy equation is explicitly solved and coupled with the momentum equation through the Boussinesq approximation. The code is conceived in a modular fashion so that different numerical methods can be used independently, the existing routines can be modified, and new ones can be included in a straightforward and sustainable manner. FluTAS is written in modern Fortran and parallelized using hybrid MPI/OpenMP in the CPU-only version and accelerated with OpenACC directives in the GPU implementation. We present different benchmarks to validate the code, and two large-scale simulations of fundamental interest in turbulent multiphase flows: isothermal emulsions in HIT and two-layer Rayleigh-Bénard convection. FluTAS is distributed through a MIT license and arises from a collaborative effort of several scientists, aiming to become a flexible tool to study complex multiphase flows.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

AutoFITS: Automatic Feature Engineering for Irregular Time Series

A time series represents a set of observations collected over time. Typically, these observations are captured with a uniform sampling frequency (e.g. daily). When data points are observed in uneven time intervals the time series is referred to as irregular or intermittent. In such scenarios, the most common solution is to reconstruct the time series to make it regular, thus removing its intermittency. We hypothesise that, in irregular time series, the time at which each observation is collected may be helpful to summarise the dynamics of the data and improve forecasting performance. We study this idea by developing a novel automatic feature engineering framework, which focuses on extracting information from this point of view, i.e., when each instance is collected. We study how valuable this information is by integrating it in a time series forecasting workflow and investigate how it compares to or complements state-of-the-art methods for regular time series forecasting. In the end, we contribute by providing a novel framework that tackles feature engineering for time series from an angle previously vastly ignored. We show that our approach has the potential to further extract more information about time series that significantly improves forecasting performance.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A volume-of-fluid method for interface-resolved simulations of phase-changing two-fluid flows

We present a numerical method for interface-resolved simulations of evaporating two-fluid flows based on the volume-of-fluid (VoF) method. The method has been implemented in an efficient FFT-based two-fluid Navier-Stokes solver, using an algebraic VoF method for the interface representation, and extended with the transport equations of thermal energy and vaporized liquid mass for the single-component evaporating liquid in an inert gas. The conservation of vaporizing liquid and computation of the interfacial mass flux are performed with the aid of a reconstructed signed-distance field, which enables the use of well-established methods for phase change solvers based on level-set methods. The interface velocity is computed with a novel approach that ensures accurate mass conservation, by constructing a divergence-free extension of the liquid velocity field onto the entire domain. The resulting approach does not depend on the type of interface reconstruction (i.e.\ can be employed in both algebraic and geometrical VoF methods). We extensively verified and validated the overall method against several benchmark cases, and demonstrated its excellent mass conservation and good overall performance for simulating evaporating two-fluid flows in two and three dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Functional renormalization group study of the critical region of the quark-meson model with vector interactions

The critical region of the two flavour quark-meson model with vector interactions is explored using the Functional Renormalization Group technique, a non-perturbative method that takes into account quantum and thermal fluctuations. Special attention is given to the low temperature and high density region of the phase diagram, which is very important for the evolution of proto-neutron stars and to construct the equation of state of compact stars. As in previous studies, without repulsive vector interaction, an unphysical region of negative entropy density is found near the first order chiral phase transition. We explore the connection between this unphysical region and the chiral critical region, especially the first order line and spinodal lines, using also different values for vector interactions. The effect of finite vector interactions in the critical region is explored. We find that the unphysical negative entropy density region appears because the $s=0$ isentropic line, near the critical region, is displaced from its $T=0$ location. For certain values of vector interactions this region is pushed to lower temperatures and high chemical potentials in such way that the negative entropy density region on the phase diagram of the model can even disapear. In the case of finite vector interactions, the location of the critical end point has a non-trivial behaviour in the $T-μ_B$ plane, which differs from that in mean-field calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

GPU acceleration of CaNS for massively-parallel direct numerical simulations of canonical fluid flows

This work presents the GPU acceleration of the open-source code CaNS for very fast massively-parallel simulations of canonical fluid flows. The distinct feature of the many-CPU Navier-Stokes solver in CaNS is its fast direct solver for the second-order finite-difference Poisson equation, based on the method of eigenfunction expansions. The solver implements all the boundary conditions valid for this type of problems in a unified framework. Here, we extend the solver for GPU-accelerated clusters using CUDA Fortran. The porting makes extensive use of CUF kernels and has been greatly simplified by the unified memory feature of CUDA Fortran, which handles the data migration between host (CPU) and device (GPU) without defining new arrays in the source code. The overall implementation has been validated against benchmark data for turbulent channel flow and its performance assessed on a NVIDIA DGX-2 system (16 Tesla V100 32Gb, connected with NVLink via NVSwitch). The wall-clock time per time step of the GPU-accelerated implementation is impressively small when compared to its CPU implementation on state-of-the-art many-CPU clusters, as long as the domain partitioning is sufficiently small that the data resides mostly on the GPUs. The implementation has been made freely available and open-source under the terms of an MIT license.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Improving Lesion Segmentation for Diabetic Retinopathy using Adversarial Learning

Diabetic Retinopathy (DR) is a leading cause of blindness in working age adults. DR lesions can be challenging to identify in fundus images, and automatic DR detection systems can offer strong clinical value. Of the publicly available labeled datasets for DR, the Indian Diabetic Retinopathy Image Dataset (IDRiD) presents retinal fundus images with pixel-level annotations of four distinct lesions: microaneurysms, hemorrhages, soft exudates and hard exudates. We utilize the HEDNet edge detector to solve a semantic segmentation task on this dataset, and then propose an end-to-end system for pixel-level segmentation of DR lesions by incorporating HEDNet into a Conditional Generative Adversarial Network (cGAN). We design a loss function that adds adversarial loss to segmentation loss. Our experiments show that the addition of the adversarial loss improves the lesion segmentation performance over the baseline.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Learned Pre-Processing for Automatic Diabetic Retinopathy Detection on Eye Fundus Images

Diabetic Retinopathy is the leading cause of blindness in the working-age population of the world. The main aim of this paper is to improve the accuracy of Diabetic Retinopathy detection by implementing a shadow removal and color correction step as a preprocessing stage from eye fundus images. For this, we rely on recent findings indicating that application of image dehazing on the inverted intensity domain amounts to illumination compensation. Inspired by this work, we propose a Shadow Removal Layer that allows us to learn the pre-processing function for a particular task. We show that learning the pre-processing function improves the performance of the network on the Diabetic Retinopathy detection task.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Magnetic field dependent 't Hooft determinant extended Nambu--Jona-Lasinio Model

We study the implications of recent lattice QCD results for the magnetic field dependence of the quarks dynamical masses on the 't Hooft determinant extended Nambu--Jona-Lasinio Model in the light and strange quark sectors (\emph{up}, \emph{down} and \emph{strange}). The parameter space is constrained at vanishing magnetic field, using the quarks dynamical masses and the meson spectra, whereas at non-vanishing magnetic field strength the dependence of the dynamical masses of two of the quark flavors is used to fit a magnetic field dependence on the model couplings, both the four-fermion Nambu--Jona-Lasinio interaction and the six-fermion 't Hooft flavor determinant. We found that this procedure reproduces the inverse magnetic catalysis, and the strength of the scalar coupling decreases with the magnetic field, while the strength of the six-fermion 't Hooft flavor determinant increases with the magnetic field.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

O-MedAL: Online Active Deep Learning for Medical Image Analysis

Active Learning methods create an optimized labeled training set from unlabeled data. We introduce a novel Online Active Deep Learning method for Medical Image Analysis. We extend our MedAL active learning framework to present new results in this paper. Our novel sampling method queries the unlabeled examples that maximize the average distance to all training set examples. Our online method enhances performance of its underlying baseline deep network. These novelties contribute significant performance improvements, including improving the model's underlying deep network accuracy by 6.30%, using only 25% of the labeled dataset to achieve baseline accuracy, reducing backpropagated images during training by as much as 67%, and demonstrating robustness to class imbalance in binary and multi-class tasks.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

One-meson-loop NJL model: Effect of collective and noncollective excitations on the quark condensate at finite temperature

We explore the effect of including quantum fluctuations in the two flavor Nambu$-$Jona-Lasinio model at finite temperature. This is accomplished, in a symmetry preserving way, by including collective and noncollective modes in the one-meson-loop gap equation which originate from poles and branch cuts in the complex plane, respectively. The inclusion of a boson cutoff, $Λ_b$, is necessary to regularize the meson-loop momenta. This new parameter is used to study the influence of going beyond the usual mean field approximation in the quark condensate. As the temperature increases, chiral symmetry tends to get restored, the collective modes melt and only noncollective modes contribute to the quark condensate. With the inclusion of such modes, the quark condensate at finite temperature has a different behavior than at mean field level that will be explored.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quark and Polyakov-loop correlations in effective models at zero and non-vanishing density

The aim of this work is to shed light on some lesser known aspects of Polyakov-loop--extended chiral models (namely the Polyakov-loop extended Nambu--Jona-Lasinio and Quark-Meson models), especially on the correlation of the quark sector with the Polyakov loop. We show that the ordering of chiral and Polyakov-loop transitions and their difference in temperature as seen in lattice QCD calculations could be realized with a critical scale of the Polyakov-loop potential that is larger than the one in pure gauge theory. The comparison of the results for the Polyakov-loop susceptibility obtained using the self-consistent medium-dependent quark mass with those obtained while keeping these masses at a fixed value allows to disentangle chiral-symmetry restoration and center-symmetry breaking effects. Furthermore, a confined chirally restored phase is identified by a plateau in the quark contribution to thermodynamics and by sigma and pion spectral functions that coincide but have a small width. We also discuss that, for some large chemical potential values, the explicit center-symmetry breaking is so strong that statistical deconfinement is realized at infinitely small temperatures. Both the missing sensitivity of the Polyakov loop to the quark mass, except at close to the chiral transition, and the Polyakov loop being zero at zero temperature at all chemical potentials, can be interpreted as indications of a missing mechanism which accounts for the quark back-reaction on the Polyakov loop.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Role of the conserved charges in the chiral symmetry restoration phase transition

The effect of conserved baryon, isospin and strangeness charges on the behavior of phase transitions in dense matter is studied. Baryonic matter is described within the three-flavor Polyakov$-$Nambu$-$Jona-Lasinio model and several charge fractions $Y_Q$ are considered. The role of the vector interaction, which can be important to describe dense systems, is discussed. Special attention is given to the case with charge fraction $Y_Q=0.4$, due to its importance in heavy-ion collisions and core-collapse supernova matter. It is shown that the possible formation of chiral-symmetric quark matter in the laboratory will be favored in asymmetric matter. Besides, the inclusion of the vector interaction reinforces the formation of quark matter at lower densities.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The strange critical endpoint and isentropic trajectories in an Extended PNJL Model with Eight Quark Interactions

In this work, we explore the possible existence of several critical endpoints in the phase diagram of strongly interacting matter using an extended PNJL model with 't Hooft determinant and eight quark interactions in the up, down and strange sectors. Besides, we also study the isentropic trajectories crossing both (light and strange) chiral phase transitions and around the critical endpoint in both the crossover and first-order transition regions.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Effects of the finite particle size in turbulent wall-bounded flows of dense suspensions

We use interface-resolved simulations to study finite-size effects in turbulent channel flow of neutrally-buoyant spheres. Two cases with particle sizes differing by a factor of 2, at the same solid volume fraction of 20% and bulk Reynolds number are considered. These are complemented with two reference single-phase flows: the unladen case, and the flow of a Newtonian fluid with the effective suspension viscosity of the same mixture in the laminar regime. As recently highlighted in Costa et al. (PRL 117, 134501), a particle-wall layer is responsible for deviations of the statistics from what is observed in the continuum limit where the suspension is modeled as a Newtonian fluid with an effective viscosity. Here we investigate the fluid and particle dynamics in this layer and in the bulk. In the particle-wall layer, the near wall inhomogeneity has an influence on the suspension micro-structure over a distance proportional to the particle size. In this layer, particles have a significant (apparent) slip velocity that is reflected in the distribution of wall shear stresses. This is characterized by extreme events (both much higher and much lower than the mean). Based on these observations we provide a scaling for the particle-to-fluid apparent slip velocity as a function of the flow parameters. We also extend the flow scaling laws in to second-order Eulerian statistics in the homogeneous suspension region away from the wall. Finite-size effects in the bulk of the channel become important for larger particles, while negligible for lower-order statistics and smaller particles. Finally, we study the particle dynamics along the wall-normal direction. Our results suggest that 1-point dispersion is dominated by particle-turbulence (and not particle-particle) interactions, while differences in 2-point dispersion and collisional dynamics are consistent with a picture of shear-driven interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Experimental investigation of turbulent suspensions of spherical particles in a square duct

We report experimental observations of turbulent flow with spherical particles in a square duct. Three particle sizes namely: $2H/d_{p}$ = 40, 16 and 9 ($2H$ being the duct full height and $d_{p}$ being the particle diameter) are investigated. The particles are nearly neutrally-buoyant with respect to the suspending fluid. Refractive Index Matched - Particle Image Velocimetry is used for fluid velocity measurement even at the highest particle volume fraction (20%) and Particle Tracking Velocimetry for the particle velocity statistics for the flows seeded with particles of the two largest sizes, whereas only pressure measurements are reported for the smallest particles. Settling effects are seen at the lowest bulk Reynolds number $Re_{2H}\approx$ 10000 whereas, at the highest $Re_{2H}\approx$ 27000, particles are in almost full suspension. The friction factor of the suspensions is found to be significantly larger than that of single-phase duct flow at the lower $Re_{2H}$ investigated; however, the difference decreases when increasing the flow rate and the total drag approaches the values of the single phase flow at the higher Reynolds number considered. The pressure drop is found to decrease with the particle diameter for volume fractions lower than $ϕ$ = 10% for nearly all $Re_{2H}$. However, at the highest volume fraction $ϕ$ = 20%, we report a peculiar non-monotonic behavior: the pressure drop first decreases and then increases with increasing particle size. The decrease of the turbulent drag with particle size at the lowest volume fractions is related to an attenuation of the turbulence. The drag increase for the two largest particle sizes at $ϕ$ = 20%, however, occurs despite this large reduction of the turbulent stresses, and it is therefore due to significant particle-induced stresses.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Turbulent channel flow of a dense binary mixture of rigid particles

We study turbulent channel flow of a binary mixture of finite-size neutrally-buoyant rigid particles by means of interface-resolved direct numerical simulations. We fix the bulk Reynolds number and total solid volume fraction, $Re_b = 5600$ and $Φ=20\%$, and vary the relative fraction of small and large particles. The binary mixture consists of particles of two different sizes, $2h/d_l=20$ and $2h/d_s=30$ where $h$ is the half channel height and $d_l$ and $d_s$ the diameter of the large and small particles. While the particulate flow statistics exhibit a significant alteration of the mean velocity profile and turbulent fluctuations with respect to the unladen flow, the differences between the mono-disperse and bi-disperse cases are small. However, we observe a clear segregation of small particles at the wall in binary mixtures, which affects the dynamics of the near wall region and thus the overall drag. This results in a higher drag in suspensions with a larger amount of large particles. As regards bi-disperse effects on the particle dynamics, a non-monotonic variation of the particle dispersion in the spanwise (homogeneous) direction is observed when increasing the percentage of small/large particles. Finally, we note that particles of the same size tend to cluster more at contact whereas the dynamics of the large particles gives highest collision kernels due to a higher approaching speed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Influence of the vector interaction and an external magnetic field on the isentropes near the chiral critical end point

The location of the critical end point (CEP) and the isentropic trajectories in the QCD phase diagram are investigated. We use the (2+1) Nambu$-$Jona-Lasinio model with the Polyakov loop coupling for different scenarios, namely by imposing zero strange quark density, which is the case in the ultra relativistic heavy-ion collisions, and $β$-equilibrium. The influence of strong magnetic fields and of the vector interaction on the isentropic trajectories around the CEP is discussed. It is shown that the vector interaction and the magnetic field, having opposite effects on the first-order transition, affect the isentropic trajectories differently: as the vector interaction increases, the first-order transition becomes weaker and the isentropes become smoother; when a strong magnetic field is considered, the first-order transition is strengthened and the isentropes are pushed to higher temperatures. No focusing of isentropes in region towards the CEP is seen.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Isentropic thermodynamics and scalar-mesons properties near the QCD critical end point

We investigate the QCD phase diagram and the location of the critical end point (CEP) in the SU(2) Polyakov$-$Nambu$-$Jona-Lasinio model with entanglement interaction giving special attention to the $π$ and $σ$-mesons properties, namely the decay widths $σ\rightarrowππ$, for several conditions around the CEP: we focus on the possible $σ\rightarrowππ$ decay along the isentropic trajectories close to the CEP since the hydrodynamical expansion of a heavy-ion collision fireball nearly follows trajectories of constant entropy. It is expected that the type of transition the dense medium goes through as it expands after the thermalization determines the behavior of this decay. It is shown that no pions are produced from the sigma decay in the chirally symmetric phase if the isentropic lines approach the first order line from chemical potentials above it. Near the CEP or above the $σ\rightarrowππ$ decay is possible with a high decay width.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical Study of the Sedimentation of Spheroidal Particles

The gravity-driven motion of rigid particles in a viscous fluid is relevant in many natural and industrial processes, yet this has mainly been investigated for spherical particles. We therefore consider the sedimentation of non-spherical (spheroidal) isolated and particle pairs in a viscous fluid via numerical simulations using the Immersed Boundary Method. The simulations performed here show that the critical Galileo number for the onset of secondary motions decreases as the spheroid aspect ratio departs from 1. Above this critical threshold, oblate particles perform a zigzagging motion whereas prolate particles rotate around the vertical axis while having their broad side facing the falling direction. Instabilities of the vortices in the wake follow when farther increasing the Galileo number. We also study the drafting-kissing-tumbling associated with the settling of particle pairs. We find that the interaction time increases significantly for non-spherical particles and, more interestingly, spheroidal particles are attracted from larger lateral displacements. This has important implications for the estimation of collision kernels and can result in increasing clustering in suspensions of sedimenting spheroids.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Two-solar-mass hybrid stars: a two model description with the Nambu-Jona-Lasinio quark model

Hybrid stars with a quark phase described by the Nambu$-$Jona-Lasinio model are studied. The hadron-quark model used to determine the stellar matter equation of state favors the appearance of quark matter: the coincidence of the deconfinement and chiral transitions and a low vacuum constituent quark mass. These two properties are essential to build equations of state that predict pure quark matter in the center of neutron stars. The effect of vector-isoscalar and vector-isovector terms is discussed, and it is shown that the vector-isoscalar terms are necessary to describe 2$M_\odot$ hybrid stars, and the vector-isovector terms result in larger quark cores and a smaller deconfinement density.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Universal Scaling Laws for Dense Particle Suspensions in Turbulent Wall-Bounded Flows

The macroscopic behavior of dense suspensions of neutrally-buoyant spheres in turbulent plane channel flow is examined. We show that particles larger than the smallest turbulence scales cause the suspension to deviate from the continuum limit in which its dynamics is well described by an effective suspension viscosity. This deviation is caused by the formation of a particle layer close to the wall with significant slip velocity. By assuming two distinct transport mechanisms in the near-wall layer and the turbulence in the bulk, we define an effective wall location such that the flow in the bulk can still be accurately described by an effective suspension viscosity. We thus propose scaling laws for the mean velocity profile of the suspension flow, together with a master equation able to predict the increase in drag as function of the particle size and volume fraction.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Collision model for fully-resolved simulations of flows laden with finite-size particles

We present a collision model for particle-particle and particle-wall interactions in interface-resolved simulations of particle-laden flows. Three types of inter-particle interactions are taken into account: (1) long- and (2) short-range hydrodynamic interactions, and (3) solid-solid contact. Long-range interactions are incorporated through an efficient and second-order accurate immersed boundary method (IBM). Short-range interactions are also partly reproduced by the IBM. However, since the IBM uses a fixed-grid, a lubrication model is needed for an inter-particle gap width smaller than the grid spacing. The lubrication model is based on asymptotic expansions of analytical solutions for canonical lubrication interactions between spheres in the Stokes regime. Roughness effects are incorporated by making the lubrication correction independent of the gap width for gap widths smaller than $\sim 1\%$ of the particle radius. This correction is applied until the particles reach solid-solid contact. To model solid-solid contact we use a variant of a linear soft-sphere collision model capable of stretching the collision time. This choice is computationally attractive because it allows to reduce the number of time steps required for integrating the collision force accurately and is physically realistic, provided that the prescribed collision time is much smaller than the characteristic timescale of particle motion. We verified the numerical implementation of our collision model and validated it against several benchmark cases for immersed head-on particle-wall and particle-particle collisions, and oblique particle-wall collisions. The results show good agreement with experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Influence of the inverse magnetic catalysis and the vector interaction in the location of the critical end point

The effect of a strong magnetic field on the location of the critical end point (CEP) in the QCD phase diagram is discussed under different scenarios. In particular, we consider the contribution of the vector interaction and take into account the inverse magnetic catalysis obtained in lattice QCD calculations at zero chemical potential. The discussion is realized within the (2+1) Polyakov--Nambu--Jona-Lasinio model. It is shown that the vector interaction and the magnetic field have opposite competing effects, and that the winning effect depends strongly on the intensity of the magnetic field. The inverse magnetic catalysis at zero chemical potential has two distinct effects for magnetic fields above $\gtrsim 0.3$ GeV$^2$: it shifts the CEP to lower chemical potentials, hinders the increase of the CEP temperature and prevents a too large increase of the baryonic density at the CEP. For fields $eB<0.1$ GeV$^2$ the competing effects between the vector contribution and the magnetic field can move the CEP to regions of temperature and density in the phase diagram that could be more easily accessible to experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Sensitivity of predictions in an effective model -- application to the chiral critical end point position in the Nambu--Jona-Lasinio model

The measurement of the position of the chiral critical end point (CEP) in the QCD phase diagram is under debate. While it is possible to predict its position by using effective models specifically built to reproduce some of the features of the underlying theory (QCD), the quality of the predictions (\textit{e.g.}, the CEP position) obtained by such effective models, depends on whether solving the model equations constitute a well- or ill-posed inverse problem. Considering these predictions as being inverse problems provides tools to evaluate if the problem is ill-conditioned, meaning that infinitesimal variations of the inputs of the model can cause comparatively large variations of the predictions. If it is ill-conditioned, it has major consequences because of finite variations that could come from experimental and/or theoretical errors. In the following, we shall apply such a reasoning on the predictions of a particular Nambu--Jona-Lasinio model within the mean field + ring approximations, with special attention to the prediction of the chiral CEP position in the $(T-μ)$ plane. We find that the problem is ill-conditioned (\idest very sensitive to input variations) for the $T$-coordinate of the CEP, whereas, it is well-posed for the $μ$-coordinate of the CEP. As a consequence, when the chiral condensate varies in a $10$ MeV range, $μ_{CEP}$ varies far less. As an illustration to understand how problematic this could be, we show that the main consequence when taking into account finite variation of the inputs, is that the existence of the CEP itself cannot be predicted anymore: for a deviation as low as 0.6 \% with respect to vacuum phenomenology (well within the estimation of the first correction to the ring approximation) the CEP may or may not exist.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The QCD phase diagram in the presence of an external magnetic field: the role of the inverse magnetic catalysis

The effect of an external magnetic field in QCD phase diagram, namely, in the the location of the critical end point (CEP) is investigated. Using the 2+1 flavor Nambu--Jona-Lasinio model with Polyakov loop, it is shown that when an external magnetic field is applied its effect on the CEP depends on the strength of the coupling. If the coupling depends on the magnetic field, allowing for inverse magnetic catalysis, the CEP moves to lower chemical potentials eventually disappearing, and the chiral restoration phase transition is always of first order.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Deconfinement, chiral symmetry restoration and thermodynamics of (2+1)--flavor hot QCD matter in an external magnetic field

The entanglement extended Polyakov--Nambu--Jona-Lasinio model at zero chemical potential in the presence of an external magnetic field is studied. The effect of the entanglement parametrization is analyzed, in particular, on the pseudocritical transition temperatures and on the thermodynamical properties of the model. The model predicts that the coincidence or not of both chiral and deconfinement transition temperatures, in the presence of an external magnetic field, depends on the entanglement parametrization chosen.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Inverse magnetic catalysis in the Polyakov-Nambu-Jona-Lasinio and entangled Polyakov-Nambu-Jona-Lasinio models

We investigate the QCD phase diagram at zero chemical potential and finite temperature in the presence of an external magnetic field within the three flavor Polyakov-Nambu-Jona-Lasinio and entangled Polyakov-Nambu-Jona-Lasinio models looking for the inverse magnetic catalysis. Two scenarios for a scalar coupling parameter dependent on the magnetic field intensity are considered. These dependencies of the coupling allow to reproduce qualitatively lattice QCD results for the quark condensates and for the Polyakov loop: due to the magnetic field the quark condensates are enhanced at low and high temperatures and suppressed for temperatures close to the transition temperatures while the Polyakov loop increases with the increasing of the magnetic field.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Phase transition and critical end point driven by an external magnetic field in asymmetric quark matter

The location of the critical end point (CEP) in the QCD phase diagram is determined under different scenarios. The effect of strangeness, isospin/charge asymmetry and an external magnetic field is investigated. The discussion is performed within the 2+1 flavor Nambu--Jona-Lasinio model with Polyakov loop. It is shown that isospin asymmetry shifts the CEP to larger baryonic chemical potentials and smaller temperatures. At large asymmetries the CEP disappears. However, a strong enough magnetic field drives the system into a first order phase transition.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quark matter subject to strong magnetic fields: phase diagram and applications

In the present work we are interested in understanding various properties of quark matter subject to strong magnetic fields described by the Nambu-Jona-Lasinio model with Polyakov loop. We start by analysing the differences arising from two different vector interactions in the Lagrangian densities, at zero temperature, and apply the results to stellar matter. We then investigate the position of the critical end point for different chemical potential and density scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Repulsive Vector Interaction in Three Flavor Magnetized Quark and Stellar Matter

The effect of the vector interaction on three flavor magnetized matter is studied within the SU(3) Nambu--Jona-Lasiono quark model. We have considered cold matter under a static external magnetic field within two different models for the vector interaction in order to investigate how the form of the vector interaction and the intensity of the magnetic field affect the equation of state as well as the strangeness content. It was shown that the flavor independent vector interaction predicts a smaller strangeness content and, therefore, harder equations of state. On the other hand, the flavor dependent vector interaction favors larger strangeness content the larger the vector coupling. We have confirmed that at low densities the magnetic field and the vector interaction have opposite competing effects: the first one softens the equation of state while the second hardens it. Quark stars and hybrid stars subject to an external magnetic field were also studied. Larger star masses are obtained for the flavor independent vector interaction. Hybrid stars may bare a core containing deconfined quarks if neither the vector interaction nor the magnetic field are too strong. Also, the presence of strong magnetic fields seems to disfavor the existence of a quark core in hybrid stars.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Strange quark chiral phase transition in hot 2+1-flavor magnetized quark matter

Using the Polyakov--Nambu--Jona-Lasinio model the strange quark chiral phase transition and the effect of its current mass on a hot magnetized three flavor quark matter at zero chemical potential is investigated. The impact of the 't Hooft mixing term on the restoration of the chiral symmetry on the light and strange sectors is studied and the critical temperature dependence on the magnetic field strength of the chiral strange transition is analyzed. It is shown that the s quark is much less sensitive to the magnetic field than the light quarks. Due to the 't Hooft term, it has a strong influence on the light quarks at all temperatures for small magnetic fields and for temperatures close to the transition temperature at zero magnetic field and above for strong magnetic fields. In particular, the large mass of the s quark makes the chiral transition of the light sector smoother and shifted to larger temperatures. A scalar coupling that weakens when the magnetic field increases will originate an inverse magnetic catalysis also in the s quark and smoothen the signature of the crossover on thermodynamical quantities such as the sound velocities or the specific heat.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The QCD critical end point driven by an external magnetic field in asymmetric quark matter

The effect of the isospin/charge asymmetry and an external magnetic field in the location of the critical end point (CEP) in the QCD phase diagram is investigated. By using the 2+1 flavor Nambu--Jona-Lasinio model with Polyakov loop (PNJL), it is shown that the isospin asymmetry shifts the CEP to larger baryonic chemical potentials and smaller temperatures, and in the presence of a large enough isospin asymmetry the CEP disappears. Nevertheless, a sufficiently high external magnetic field can drive the system into a first order phase transition again.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Deconfinement and chiral restoration within the SU(3) Polyakov--Nambu--Jona-Lasinio and entangled Polyakov--Nambu--Jona-Lasinio models in an external magnetic field

The behavior of the quark condensates at zero chemical potential and finite temperature subject to an external magnetic field is studied within the three flavor Nambu--Jona-Lasinio model with Polyakov loop (PNJL) and its extension, the so-called entangled PNJL model (EPNJL). A comparison with recent lattice QCD data is performed and it is shown that at $T=0$ MeV the light quark condensates are in quantitative agreement. At finite temperature, although there is an overall reasonable agreement with several lattice results, it is shown that in the lattice calculations the effect due to the electric charge quark difference is stronger and the restoration of the $u$ quark chiral symmetry starts at lower temperatures. When considering the EPNJL model with a Polyakov loop scale parameter that depends on the magnetic field, it is possible to obtain an earlier rise of the Polyakov loop with the increase of the magnetic field and due to the entanglement, the inverse magnetic catalysis is found as in the lattice QCD calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effects of entanglement and instanton suppression at finite temperature in a SU(2) EPNJL model with anomaly

We investigate the phase transitions characterized by deconfinement and restoration of chiral and axial symmetries,at finite temperature, in the framework of QCD inspired models. We compare the results obtained in the SU(2) Polyakov-Nambu-Jona-Lasinio model with anomaly and in its extended version, the Entangled Polyakov-Nambu-Jona-Lasinio model. In the last version, four-quark vertices with entanglement between the chiral condensate and the Polyakov loop are considered. The thermodynamics of the phase transitions, the meson spectrum, and in particular the convergence of axial and chiral partners, will be analyzed, as well as the topological susceptibility. We find that an explicit temperature dependence of the coupling vertices is necessary in both models in order to have effective restoration of the U$_A$(1) symmetry.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Restoration of axial symmetry and its possible relation with restoration of chiral symmetry and deconfinement at finite temperature

The phase transitions characterized by deconfinement and restoration of chiral symmetry as well as the restoration of axial symmetry, at finite temperature, are investigated in the framework of SU(2) Polyakov-Nambu-Jona-Lasinio (PNJL) models with the U$_A$(1) anomaly. The thermodynamics of the phase transitions, the topological susceptibility, the meson spectrum, and, in particular, the convergence of axial and chiral partners are analyzed, in the framework of the ordinary PNJL model and in its extension, the entangled Polyakov-Nambu-Jona-Lasinio (EPNJL) model. The latter incorporates entanglement between restoration of chiral symmetry and deconfinement.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Falling chains as variable mass systems: theoretical model and experimental analysis

In the present paper we revisit, theoretical and experimentally, the fall of a folded U-chain and of a pile-chain. The model calculation implies the division of the whole system into two subsystems of variable mass, allowing to explore the role of tensional contact forces at the boundary of the subsystems. This justifies, for instance, that the folded U-chain falls faster than the acceleration due to the gravitational force. This result, which matches quite well with the experimental data independently of the type of chain, implies that the falling chain is well described by energy conservation. We verify that these conclusions are not observed for the pile-chain motion.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Interplay between chiral and axial symmetries in a SU(2) Nambu--Jona-Lasinio Model with the Polyakov loop

We consider a two flavor Polyakov--Nambu--Jona-Lasinio (PNJL) model where the Lagrangian includes an interaction term that explicitly breaks the U$_A(1)$ anomaly. At finite temperature, the restoration of chiral and axial symmetries, signaled by the behavior of several observables, is investigated. We compare the effects of two regularizations at finite temperature, one of them, that allows high momentum quarks states, leading to the full recovery of chiral symmetry. From the analysis of the behavior of the topological susceptibility and of the mesonic masses of the axial partners, it is found in the SU(2) model that, unlike the SU(3) results, the recovery of the axial symmetry is not a consequence of the full recovery of the chiral symmetry. Thus, one needs to use an additional idea, by means of a temperature dependence of the anomaly coefficient, that simulates instanton suppression effects.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The QCD Critical End Point in the Context of the Polyakov--Nambu--Jona-Lasinio Model

We investigate the phase diagram of the so-called Polyakov--Nambu--Jona-Lasinio model at finite temperature and nonzero chemical potential with three quark flavors. Chiral and deconfinement phase transitions are discussed, and the relevant order-like parameters are analyzed. A special attention is payed to the critical end point (CEP): the influence of the strangeness on the location of the CEP is studied; also the strength of the flavor-mixing interaction alters the CEP location, once when it becomes weaker the CEP moves to low temperatures and can even disappear.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Exploring the role of model parameters and regularization procedures in the thermodynamics of the PNJL model

The equation of state and the critical behavior around the critical end point are studied in the context of the Polyakov--Nambu--Jona--Lasinio model. We prove that a convenient choice of the model parameters is crucial to get the correct description of isentropic trajectories. The physical relevance of the effects of the regularization procedure is insured by the agreement with general thermodynamic requirements. The results are compared with simple thermodynamic expectations and lattice data.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

How parameters and regularization affect the PNJL model phase diagram and thermodynamic quantities

We explore the phase diagram and the critical behavior of QCD thermodynamic quantities in the context of the so-called Polyakov--Nambu--Jona-Lasinio model. We show that this improved field theoretical model is a successful candidate for studying the equation of state and the critical behavior around the critical end point. We argue that a convenient choice of the model parameters is crucial to get the correct description of isentropic trajectories. The effects of the regularization procedure in several thermodynamic quantities is also analyzed. The results are compared with simple thermodynamic expectations and lattice data.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

How the Polyakov loop and the regularization affect strangeness and restoration of symmetries at finite T

The effects of the Polyakov loop and of a regularization procedure that allows the presence of high momentum quark states at finite temperature is investigated within the Polyakov-loop extended Nambu--Jona-Lasinio model. The characteristic temperatures, as well as the behavior of observables that signal deconfinement and restoration of chiral and axial symmetries, are analyzed, paying special attention to the behavior of strangeness degrees of freedom. We observe that the cumulative effects of the Polyakov loop and of the regularization procedure contribute to a better description of the thermodynamics, as compared with lattice estimations. We find a faster partial restoration of chiral symmetry and the restoration of the axial symmetry appears as a natural consequence of the full recovering of the chiral symmetry that was dynamically broken. These results show the relevance of the effects of the interplay among the Polyakov loop dynamics, the high momentum quark sates and the restoration of the chiral and axial symmetries at finite temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Phase diagram and critical properties in the Polyakov--Nambu--Jona-Lasinio model

We investigate the phase diagram of the so-called Polyakov--Nambu--Jona-Lasinio model at finite temperature and nonzero chemical potential with three quark flavours. Chiral and deconfinement phase transitions are discussed, and the relevant order-like parameters are analyzed. The results are compared with simple thermodynamic expectations and lattice data. A special attention is payed to the critical end point: as the strength of the flavour-mixing interaction becomes weaker, the critical end point moves to low temperatures and can even disappear.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Phase diagram and critical properties within an effective model of QCD: the Nambu-Jona-Lasinio model coupled to the Polyakov loop

We investigate the phase diagram of the so-called Polyakov--Nambu--Jona-Lasinio model at finite temperature and non-zero chemical potential with three quark flavors. Chiral and deconfinement phase transitions are discussed and the relevant order-like parameters are analyzed. The results are compared with simple thermodynamic expectations and lattice data. We present the phase diagram in the ($T,\,μ_B)$ plane, paying special attention to the critical end point: as the strength of the flavor-mixing interaction becomes weaker, the critical end point moves to low temperatures and can even disappear.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

What Does Low Energy Physics Tell Us About The Zero Momentum Gluon Propagator

The connection between QCD, a nonlocal Nambu--Jona-Lasinio type model and the Landau gauge gluon propagator is explored. This two point function is parameterized by a functional form which is compatible with Dyson-Schwinger and lattice QCD results. Demanding the nonlocal model to reproduce the experimental values for the pion mass, the pion decay constant, $Γ_{π\rightarrow γγ}$ and the light quark condensate we conclude that low energy physics does not distinguish between the so-called decoupling and scaling solutions of the Dyson-Schwinger equations. This result means that, provided that the model parameters are chosen appropriately, one is free to chose any of the above scenarios. Furthermore, the nonlocal Nambu--Jona-Lasinio quark model considered here is chiral invariant and satisfies the GMOR relation at the 1% level of precision.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Behavior of physical observables in the vicinity of the QCD critical end point

Using the SU(3) Nambu-Jona-Lasinio (NJL) model, we study the chiral phase transition at finite $T$ and $μ_B$. Special attention is given to the QCD critical end point (CEP): the study of physical quantities, as the pressure, the entropy, the baryon number susceptibility and the specific heat near the CEP, will provide complementary information concerning the order of the phase transition. We also analyze the information provided by the study of the critical exponents around the CEP.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Effective restoration of the U_A(1) symmetry with temperature and density

We investigate the full U(3)$\otimes$U(3) chiral symmetry restoration, at finite temperature and density, on the basis of a quark model which incorporates the most relevant properties of QCD in this context: explicit and spontaneous breaking of chiral symmetry and axial U$_A$(1) symmetry breaking. A specific lattice-inspired behavior of the topological susceptibility, combined with the convergence of chiral partners, signals the onset of an effective chiral symmetry restoration. The results suggest that the axial part of the symmetry is restored before the possible restoration of the full U(3)$\otimes$U(3) chiral symmetry can occur. This conclusion is valid in the context of both finite temperature and density.

0 citations0 reviewsNo code linkedOpen access