Desenvolupament de codis numèrics per a la resolució de fluxos turbulents en ordinadors paral·lels emprant discretitzacions tipus symmetry-preserving

Projecte de recerca elaborat a partir d’una estada a la University of Groningen, Holanda, entre 2007 i 2009. La simulació directa de la turbulència (DNS) és una eina clau dins de la mecànica de fluids computacional. Per una banda permet conèixer millor la física de la turbulència i per l'altra els resultats obtinguts són claus per el desenvolupament dels models de turbulència. No obstant, el DNS no és una tècnica vàlida per a la gran majoria d'aplicacions industrials degut al elevats costos computacionals. Per tant, és necessari cert grau de modelització de la turbulència. En aquest context, s'han introduïts importants millores basades en la modelització del terme convectiu (no lineal) emprant symmetry-preserving regularizations. En tracta de modificar adequadament el terme convectiu a fi de reduir la producció d'escales més i més petites (vortex-stretching) tot mantenint tots els invariants de les equacions originals. Fins ara, aquest models s'han emprat amb èxit per nombres de Rayleigh (Ra) relativament elevats. En aquest punt, disposar de resultats DNS per a configuracions més complexes i nombres de Ra més elevats és clau. En aquest contexte, s'han dut a terme simulacions DNS en el supercomputador MareNostrum d'una Differentially Heated Cavity amb Ra=1e11 i Pr=0.71 durant el primer any dels dos que consta el projecte. A més a més, s'ha adaptat el codi a fi de poder simular el fluxe al voltant d'un cub sobre una pared amb Re=10000. Aquestes simulacions DNS són les més grans fetes fins ara per aquestes configuracions i la seva correcta modelització és un gran repte degut la complexitat dels fluxes. Aquestes noves simulacions DNS estan aportant nous coneixements a la física de la turbulència i aportant resultats indispensables per al progrés de les modelitzacións tipus symmetry-preserving regularization.

Sovereign debt markets in turbulent times: Creditor discrimination and crowding-out effects

In 2007, countries in the Euro periphery were enjoying stable growth, low deficits, and lowspreads. Then the financial crisis erupted and pushed them into deep recessions, raising theirdeficits and debt levels. By 2010, they were facing severe debt problems. Spreads increased and,surprisingly, so did the share of the debt held by domestic creditors. Credit was reallocatedfrom the private to the public sectors, reducing investment and deepening the recessions evenfurther. To account for these facts, we propose a simple model of sovereign risk in which debtcan be traded in secondary markets. The model has two key ingredients: creditor discriminationand crowding-out effects. Creditor discrimination arises because, in turbulent times, sovereigndebt offers a higher expected return to domestic creditors than to foreign ones. This providesincentives for domestic purchases of debt. Crowding-out effects arise because private borrowingis limited by financial frictions. This implies that domestic debt purchases displace productiveinvestment. The model shows that these purchases reduce growth and welfare, and may lead toself-fulfilling crises. It also shows how crowding-out effects can be transmitted to other countriesin the Eurozone, and how they may be addressed by policies at the European level.

Diffusion in spatially and temporarily inhomogeneous media: Effects of turbulent mixing

We consider diffusion of a passive substance C in a phase-separating nonmiscible binary alloy under turbulent mixing. The substance is assumed to have different diffusion coefficients in the pure phases A and B, leading to a spatially and temporarily dependent diffusion ¿coefficient¿ in the diffusion equation plus convective term. In this paper we consider especially the effects of a turbulent flow field coupled to both the Cahn-Hilliard type evolution equation of the medium and the diffusion equation (both, therefore, supplemented by a convective term). It is shown that the formerly observed prolonged anomalous diffusion [H. Lehr, F. Sagués, and J.M. Sancho, Phys. Rev. E 54, 5028 (1996)] is no longer seen if a flow of sufficient intensity is supplied.

Particle dispersion in synthetic turbulent flows

We study particle dispersion advected by a synthetic turbulent flow from a Lagrangian perspective and focus on the two-particle and cluster dispersion by the flow. It has been recently reported that Richardson¿s law for the two-particle dispersion can stem from different dispersion mechanisms, and can be dominated by either diffusive or ballistic events. The nature of the Richardson dispersion depends on the parameters of our flow and is discussed in terms of the values of a persistence parameter expressing the relative importance of the two above-mentioned mechanisms. We support this analysis by studying the distribution of interparticle distances, the relative velocity correlation functions, as well as the relative trajectories.

Langevin approach to generate synthetic turbulent flows

We present an analytical scheme, easily implemented numerically, to generate synthetic Gaussian turbulent flows by using a linear Langevin equation, where the noise term acts as a stochastic stirring force. The characteristic parameters of the velocity field are well introduced, in particular the kinematic viscosity and the spectrum of energy. As an application, the diffusion of a passive scalar is studied for two different energy spectra. Numerical results are compared favorably with analytical calculations.

Front dynamics in turbulent media

A study of a stable front propagating in a turbulent medium is presented. The front is generated through a reaction-diffusion equation, and the turbulent medium is statistically modeled using a Langevin equation. Numerical simulations indicate the presence of two different dynamical regimes. These regimes appear when the turbulent flow either wrinkles a still rather sharp propagating interfase or broadens it. Specific dependences of the propagating velocities on stirring intensities appropriate to each case are found and fitted when possible according to theoretically predicted laws. Different turbulent spectra are considered.

Penrose limits, worldvolume fluxes, and supersymmetry

It is noted that a finite Penrose limit for brane probes with nonzero worldvolume fluxes does not generically exist; this is closely related to the observation by Blau and co-workers that for a brane probe the Penrose limit is equivalent to an infinite-tension limit. It is shown that when the limit exists, however, the number of supersymmetries preserved by the probe does not decrease.

Self-similarity properties of natural images resemble those of turbulent flows

We show that the statistics of an edge type variable in natural images exhibits self-similarity properties which resemble those of local energy dissipation in turbulent flows. Our results show that self-similarity and extended self-similarity hold remarkably for the statistics of the local edge variance, and that the very same models can be used to predict all of the associated exponents. These results suggest using natural images as a laboratory for testing more elaborate scaling models of interest for the statistical description of turbulent flows. The properties we have exhibited are relevant for the modeling of the early visual system: They should be included in models designed for the prediction of receptive fields.

Diffusion in spatially and temporarily inhomogeneous media: Effects of turbulent mixing

We consider diffusion of a passive substance C in a phase-separating nonmiscible binary alloy under turbulent mixing. The substance is assumed to have different diffusion coefficients in the pure phases A and B, leading to a spatially and temporarily dependent diffusion ¿coefficient¿ in the diffusion equation plus convective term. In this paper we consider especially the effects of a turbulent flow field coupled to both the Cahn-Hilliard type evolution equation of the medium and the diffusion equation (both, therefore, supplemented by a convective term). It is shown that the formerly observed prolonged anomalous diffusion [H. Lehr, F. Sagués, and J.M. Sancho, Phys. Rev. E 54, 5028 (1996)] is no longer seen if a flow of sufficient intensity is supplied.

Particle dispersion in synthetic turbulent flows

We study particle dispersion advected by a synthetic turbulent flow from a Lagrangian perspective and focus on the two-particle and cluster dispersion by the flow. It has been recently reported that Richardson¿s law for the two-particle dispersion can stem from different dispersion mechanisms, and can be dominated by either diffusive or ballistic events. The nature of the Richardson dispersion depends on the parameters of our flow and is discussed in terms of the values of a persistence parameter expressing the relative importance of the two above-mentioned mechanisms. We support this analysis by studying the distribution of interparticle distances, the relative velocity correlation functions, as well as the relative trajectories.

Dispersion of passive particles by a quasi-two-dimensional turbulent flow

Using the experimental data of Paret and Tabeling [Phys. Rev. Lett. 79, 4162 (1997)] we consider in detail the dispersion of particle pairs by a two-dimensional turbulent flow and its relation to the kinematic properties of the velocity field. We show that the mean square separation of a pair of particles is governed by rather rare, extreme events and that the majority of initially close pairs are not dispersed by the flow. Another manifestation of the same effect is the fact that the dispersion of an initially dense cluster is not the result of homogeneously spreading the particles within the whole system. Instead it proceeds through a splitting into smaller but also dense clusters. The statistical nature of this effect is discussed.

Langevin approach to generate synthetic turbulent flows

We present an analytical scheme, easily implemented numerically, to generate synthetic Gaussian turbulent flows by using a linear Langevin equation, where the noise term acts as a stochastic stirring force. The characteristic parameters of the velocity field are well introduced, in particular the kinematic viscosity and the spectrum of energy. As an application, the diffusion of a passive scalar is studied for two different energy spectra. Numerical results are compared favorably with analytical calculations.

Front dynamics in turbulent media

A study of a stable front propagating in a turbulent medium is presented. The front is generated through a reaction-diffusion equation, and the turbulent medium is statistically modeled using a Langevin equation. Numerical simulations indicate the presence of two different dynamical regimes. These regimes appear when the turbulent flow either wrinkles a still rather sharp propagating interfase or broadens it. Specific dependences of the propagating velocities on stirring intensities appropriate to each case are found and fitted when possible according to theoretically predicted laws. Different turbulent spectra are considered.