BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Pedro Costa

Pedro Costa contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.flu-dynhep-phComputer Visioneess.IVnucl-thphysics.comp-phMachine LearningComputational Engineering, Finance, and Science

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
20works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

20 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

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
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