Stochastic generation of homogeneous isotropic turbulence with well-defined-spectra

A precise and simple computational model to generate well-behaved two-dimensional turbulent flows is presented. The whole approach rests on the use of stochastic differential equations and is general enough to reproduce a variety of energy spectra and spatiotemporal correlation functions. Analytical expressions for both the continuous and the discrete versions, together with simulation algorithms, are derived. Results for two relevant spectra, covering distinct ranges of wave numbers, are given.

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.

Phase-field model for Hele-Shaw flows with arabitrary contrast II: numerical approach

We present a phase-field model for the dynamics of the interface between two inmiscible fluids with arbitrary viscosity contrast in a rectangular Hele-Shaw cell. With asymptotic matching techniques we check the model to yield the right Hele-Shaw equations in the sharp-interface limit, and compute the corrections to these equations to first order in the interface thickness. We also compute the effect of such corrections on the linear dispersion relation of the planar interface. We discuss in detail the conditions on the interface thickness to control the accuracy and convergence of the phase-field model to the limiting Hele-Shaw dynamics. In particular, the convergence appears to be slower for high viscosity contrasts.

Interface dynamics in Hele-Shaw flows with centrifugal forces. Preventing cusp singularities with rotation

A class of exact solutions of Hele-Shaw flows without surface tension in a rotating cell is reported. We show that the interplay between injection and rotation modifies the scenario of formation of finite-time cusp singularities. For a subclass of solutions, we show that, for any given initial condition, there exists a critical rotation rate above which cusp formation is suppressed. We also find an exact sufficient condition to avoid cusps simultaneously for all initial conditions within the above subclass.

Viscous fingering as a paradigm of interfacial pattern formation: Recent results and new challenges

We review recent results on dynamical aspects of viscous fingering. The Saffman¿Taylor instability is studied beyond linear stability analysis by means of a weakly nonlinear analysis and the exact determination of the subcritical branch. A series of contributions pursuing the idea of a dynamical solvability scenario associated to surface tension in analogy with the traditional selection theory is put in perspective and discussed in the light of the asymptotic theory of Tanveer and co-workers. The inherently dynamical singular effects of surface tension are clarified. The dynamical role of viscosity contrast is explored numerically. We find that the basin of attraction of the Saffman¿Taylor finger depends on viscosity contrast, and that the sensitivity to this parameter is maximal in the usual limit of high viscosity contrast. The competing attractors are identified as closed bubble solutions. We briefly report on recent results and work in progress concerning rotating Hele-Shaw flows, topological singularities and wetting effects, and also discuss future directions in the context of viscous fingering

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.

Diffusion of passive scalars under stochastic convection

The diffusion of passive scalars convected by turbulent flows is addressed here. A practical procedure to obtain stochastic velocity fields with well¿defined energy spectrum functions is also presented. Analytical results are derived, based on the use of stochastic differential equations, where the basic hypothesis involved refers to a rapidly decaying turbulence. These predictions are favorable compared with direct computer simulations of stochastic differential equations containing multiplicative space¿time correlated noise.

Were Spanish migrants attracted by industrial agglomerations?: an analysis for the interwar years in the light of the new economic geography

In this paper we examine whether access to markets had a significant influence onmigration choices of Spanish internal migrants in the inter-war years. We perform astructural contrast of a New Economic Geography model that focus on the forwardlinkage that links workers location choice with the geography of industrial production,one of the centripetal forces that drive agglomeration in the NEG models. The resultshighlight the presence of this forward linkage in the Spanish economy of the inter-warperiod. That is, we prove the existence of a direct relation between workers¿ localizationdecisions and the market potential of the host regions. In addition, the direct estimationof the values associated with key parameters in the NEG model allows us to simulatethe migratory flows derived from different scenarios of the relative size of regions andthe distances between them. We show that in Spain the power of attraction of theagglomerations grew as they increased in size, but the high elasticity estimated for themigration costs reduced the intensity of the migratory flows. This could help to explainthe apparently low intensity of internal migrations in Spain until its upsurge during the1920s. This also explains the geography of migrations in Spain during this period,which hardly affected the regions furthest from the large industrial agglomerations (i.e.,regions such as Andalusia, Estremadura and Castile-La Mancha) but had an intenseeffect on the provinces nearest to the principal centres of industrial development.

Stochastic generation of homogeneous isotropic turbulence with well-defined-spectra

A precise and simple computational model to generate well-behaved two-dimensional turbulent flows is presented. The whole approach rests on the use of stochastic differential equations and is general enough to reproduce a variety of energy spectra and spatiotemporal correlation functions. Analytical expressions for both the continuous and the discrete versions, together with simulation algorithms, are derived. Results for two relevant spectra, covering distinct ranges of wave numbers, are given.

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.

Velocity fluctuations of fracturelike disruptions of associating polymer solutions

Velocity has been measured as a function of time for propagating crack tips as water is injected into solutions of end-capped associating polymers in a rectanguar Hele-Shaw cell. Measurements were performed for flows with different values of cell gap, channel width, polymer molecular weight, and polymer concentration. The condition for the onset of fracturelike behavior is well described by a Deborah number which uses the shear-thinning shear rate of the polymer solution as a characteristic frequency for network relaxation. At low molecular weight, the onset of fracturelike pattern evolution is accompanied by an abrupt jump in tip velocity, followed by a lower and approximately constant acceleration. At high molecular weight, the transition to fracturelike behavior involves passing through a regime that may be understood in terms of stick-slip dynamics. The crack-tip wanders from side to side and fluctuates (in both speed and velocity along the channel) with a characteristic frequency which depends linearly on the invading fluid injection rate.

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.

Diffusion of passive scalars under stochastic convection

The diffusion of passive scalars convected by turbulent flows is addressed here. A practical procedure to obtain stochastic velocity fields with well¿defined energy spectrum functions is also presented. Analytical results are derived, based on the use of stochastic differential equations, where the basic hypothesis involved refers to a rapidly decaying turbulence. These predictions are favorable compared with direct computer simulations of stochastic differential equations containing multiplicative space¿time correlated noise.

Were Spanish migrants attracted by industrial agglomerations?: an analysis for the interwar years in the light of the new economic geography

In this paper we examine whether access to markets had a significant influence onmigration choices of Spanish internal migrants in the inter-war years. We perform astructural contrast of a New Economic Geography model that focus on the forwardlinkage that links workers location choice with the geography of industrial production,one of the centripetal forces that drive agglomeration in the NEG models. The resultshighlight the presence of this forward linkage in the Spanish economy of the inter-warperiod. That is, we prove the existence of a direct relation between workers¿ localizationdecisions and the market potential of the host regions. In addition, the direct estimationof the values associated with key parameters in the NEG model allows us to simulatethe migratory flows derived from different scenarios of the relative size of regions andthe distances between them. We show that in Spain the power of attraction of theagglomerations grew as they increased in size, but the high elasticity estimated for themigration costs reduced the intensity of the migratory flows. This could help to explainthe apparently low intensity of internal migrations in Spain until its upsurge during the1920s. This also explains the geography of migrations in Spain during this period,which hardly affected the regions furthest from the large industrial agglomerations (i.e.,regions such as Andalusia, Estremadura and Castile-La Mancha) but had an intenseeffect on the provinces nearest to the principal centres of industrial development.

Modeling storm events to investigate the influence of the stream¿catchment interface zone on stream biogeochemistry

We formulate a new mixing model to explore hydrological and chemical conditions under which the interface between the stream and catchment interface (SCI) influences the release of reactive solutes into stream water during storms. Physically, the SCI corresponds to the hyporheic/riparian sediments. In the new model this interface is coupled through a bidirectional water exchange to the conventional two components mixing model. Simulations show that the influence of the SCI on stream solute dynamics during storms is detectable when the runoff event is dominated by the infiltrated groundwater component that flows through the SCI before entering the stream and when the flux of solutes released from SCI sediments is similar to, or higher than, the solute flux carried by the groundwater. Dissolved organic carbon (DOC) and nitrate data from two small Mediterranean streams obtained during storms are compared to results from simulations using the new model to discern the circumstances under which the SCI is likely to control the dynamics of reactive solutes in streams. The simulations and the comparisons with empirical data suggest that the new mixing model may be especially appropriate for streams in which the periodic, or persistent, abrupt changes in the level of riparian groundwater exert hydrologic control on flux of biologically reactive fluxes between the riparian/hyporheic compartment and the stream water.