Stochastic differential equations with random coefficients

In this paper we establish the existence and uniqueness of a solution for different types of stochastic differential equation with random initial conditions and random coefficients. The stochastic integral is interpreted as a generalized Stratonovich integral, and the techniques used to derive these results are mainly based on the path properties of the Brownian motion, and the definition of the Stratonovich integral.

Large deviations for stochastic Volterra equations

This paper is devoted to prove a large-deviation principle for solutions to multidimensional stochastic Volterra equations.

Stochastic delay differential equations driven by fractional Brownian motion with Hurst parameter H 1/2

We consider the Cauchy problem for a stochastic delay differential equation driven by a fractional Brownian motion with Hurst parameter H>¿. We prove an existence and uniqueness result for this problem, when the coefficients are sufficiently regular. Furthermore, if the diffusion coefficient is bounded away from zero and the coefficients are smooth functions with bounded derivatives of all orders, we prove that the law of the solution admits a smooth density with respect to Lebesgue measure on R.

Uniform estimates in the Poincare-Aronszajn theorem on the separation of singularities of analytic functions

We study the possibility of splitting any bounded analytic function $f$ with singularities in a closed set $E\cup F$ as a sum of two bounded analytic functions with singularities in $E$ and $F$ respectively. We obtain some results under geometric restrictions on the sets $E$ and $F$ and we provide some examples showing the sharpness of the positive results.

Generation and nonlinear evolution of shore-oblique/transverse sand bars

The coupling between topography, waves and currents in the surf zone may selforganize to produce the formation of shore-transverse or shore-oblique sand bars on an otherwise alongshore uniform beach. In the absence of shore-parallel bars, this has been shown by previous studies of linear stability analysis, but is now extended to the finite-amplitude regime. To this end, a nonlinear model coupling wave transformation and breaking, a shallow-water equations solver, sediment transport and bed updating is developed. The sediment flux consists of a stirring factor multiplied by the depthaveraged current plus a downslope correction. It is found that the cross-shore profile of the ratio of stirring factor to water depth together with the wave incidence angle primarily determine the shape and the type of bars, either transverse or oblique to the shore. In the latter case, they can open an acute angle against the current (upcurrent oriented) or with the current (down-current oriented). At the initial stages of development, both the intensity of the instability which is responsible for the formation of the bars and the damping due to downslope transport grow at a similar rate with bar amplitude, the former being somewhat stronger. As bars keep on growing, their finite-amplitude shape either enhances downslope transport or weakens the instability mechanism so that an equilibrium between both opposing tendencies occurs, leading to a final saturated amplitude. The overall shape of the saturated bars in plan view is similar to that of the small-amplitude ones. However, the final spacings may be up to a factor of 2 larger and final celerities can also be about a factor of 2 smaller or larger. In the case of alongshore migrating bars, the asymmetry of the longshore sections, the lee being steeper than the stoss, is well reproduced. Complex dynamics with merging and splitting of individual bars sometimes occur. Finally, in the case of shore-normal incidence the rip currents in the troughs between the bars are jet-like while the onshore return flow is wider and weaker as is observed in nature.

On the approximation of the subgrid scale for systems of equations

Postprint (published version)

Kinetic equations for difussion in the presence of entropic barriers

We use the mesoscopic nonequilibrium thermodynamics theory to derive the general kinetic equation of a system in the presence of potential barriers. The result is applied to a description of the evolution of systems whose dynamics is influenced by entropic barriers. We analyze in detail the case of diffusion in a domain of irregular geometry in which the presence of the boundaries induces an entropy barrier when approaching the exact dynamics by a coarsening of the description. The corresponding kinetic equation, named the Fick-Jacobs equation, is obtained, and its validity is generalized through the formulation of a scaling law for the diffusion coefficient which depends on the shape of the boundaries. The method we propose can be useful to analyze the dynamics of systems at the nanoscale where the presence of entropy barriers is a common feature.

Integro-difference equations for interacting species and the Neolithic transition

We introduce a set of sequential integro-difference equations to analyze the dynamics of two interacting species. Firstly, we derive the speed of the fronts when a species invades a space previously occupied by a second species, and check its validity by means of numerical random-walk simulations. As an example, we consider the Neolithic transition: the predictions of the model are consistent with the archaeological data for the front speed, provided that the interaction parameter is low enough. Secondly, an equation for the coexistence time between the invasive and the invaded populations is obtained for the first time. It agrees well with the simulations, is consistent with observations of the Neolithic transition, and makes it possible to estimate the value of the interaction parameter between the incoming and the indigenous populations

Fronts from integrodifference equations and persistence effects on the Neolithic transition

We extend a previous model of the Neolithic transition in Europe [J. Fort and V. Méndez, Phys. Rev. Lett. 82, 867 (1999)] by taking two effects into account: (i) we do not use the diffusion approximation (which corresponds to second-order Taylor expansions), and (ii) we take proper care of the fact that parents do not migrate away from their children (we refer to this as a time-order effect, in the sense that it implies that children grow up with their parents, before they become adults and can survive and migrate). We also derive a time-ordered, second-order equation, which we call the sequential reaction-diffusion equation, and use it to show that effect (ii) is the most important one, and that both of them should in general be taken into account to derive accurate results. As an example, we consider the Neolithic transition: the model predictions agree with the observed front speed, and the corrections relative to previous models are important (up to 70%)

JPEG standard uniform quantization error modeling with applications to sequential and progressive operation modes

In this paper we propose a method for computing JPEG quantization matrices for a given mean square error or PSNR. Then, we employ our method to compute JPEG standard progressive operation mode definition scripts using a quantization approach. Therefore, it is no longer necessary to use a trial and error procedure to obtain a desired PSNR and/or definition script, reducing cost. Firstly, we establish a relationship between a Laplacian source and its uniform quantization error. We apply this model to the coefficients obtained in the discrete cosine transform stage of the JPEG standard. Then, an image may be compressed using the JPEG standard under a global MSE (or PSNR) constraint and a set of local constraints determined by the JPEG standard and visual criteria. Secondly, we study the JPEG standard progressive operation mode from a quantization based approach. A relationship between the measured image quality at a given stage of the coding process and a quantization matrix is found. Thus, the definition script construction problem can be reduced to a quantization problem. Simulations show that our method generates better quantization matrices than the classical method based on scaling the JPEG default quantization matrix. The estimation of PSNR has usually an error smaller than 1 dB. This figure decreases for high PSNR values. Definition scripts may be generated avoiding an excessive number of stages and removing small stages that do not contribute during the decoding process with a noticeable image quality improvement.

Kumar's algorithm for solving Toeplitz systems of equations over finite fields

Adaptació de l'algorisme de Kumar per resoldre sistemes d'equacions amb matrius de Toeplitz sobre els reals a cossos finits en un temps 0 (n log n).

An expression of mixed animal model equations to account for different means and variances in the base population

This paper presents a general expression to predict breeding values using animal models when the base population is selected, i.e. the means and variances of breeding values in the base generation differ among individuals. Rules for forming the mixed model equations are also presented. A numerical example illustrates the procedure.

Weierstrass integrability of differential equations

The integrability problem consists in finding the class of functions a first integral of a given planar polynomial differential system must belong to. We recall the characterization of systems which admit an elementary or Liouvillian first integral. We define {\it Weierstrass integrability} and we determine which Weierstrass integrable systems are Liouvillian integrable. Inside this new class of integrable systems there are non--Liouvillian integrable systems.

Flow resistance equations for mountain rivers

Three models of flow resistance (a Keulegan-type logarithmic law and two models developed for large-scale roughness conditions: the full logarithmic law and a model based on an inflectional velocity profile) were calibrated, validated and compared using an extensive database (N = 1,533) from rivers and flumes, representative of a wide hydraulic and geomorphologic range in the field of gravel-bed and mountain channels. It is preferable to apply the model based on an inflectional velocity profile in the relative submergence (y/d90) interval between 0.5 and 15, while the full logarithmic law is preferable for values below 0.5. For high relative submergence, above 15, either the logarithmic law or the full logarithmic law can be applied. The models fitted to the coarser percentiles are preferable to those fitted to the median diameter, owing to the higher explanatory power achieved by setting a model, the smaller difference in the goodness-of-fit between the different models and the lower influence of the origin of the data (river or flume).

Predicción de la resistencia al flujo en ríos de montaña

Se han calibrado, validado y comparado tres modelos de resistencia al flujo de contorno granular: un modelo potencial y otros dos modelos desarrollados para condiciones de alta rugosidad relativa (uno basado en una modificación de la ley logarítmica de Prandtl-von Karman y otro fundamentado en un perfil de velocidad configurado en dos zonas: una uniforme en las proximidades de los elementos de rugosidad y otra superior que sigue una distribución logarítmica). Se ha empleado para ello un numeroso conjunto de 1.533 datos tomados en ríos y en canales de laboratorio, representativo de un amplio intervalo hidráulico y geomorfológico en el ámbito de ríos de grava y de montaña. Han resultado preferibles las ecuaciones ajustadas con los percentiles granulométricos mayores (d90 o d84) que las ajustadas con el diámetro mediano (d50), debido a la mayor capacidad explicativa alcanzada dado un modelo, la menor diferencia en la bondad de ajuste entre los diferentes modelos y la menor influencia del origen de los datos (río o canal de laboratorio). Las ecuaciones ajustadas de acuerdo con los modelos en donde se contemplan condiciones de alta rugosidad relativa presentan predicciones similares, exceptuando el intervalo macrorrugoso (y/d90 < 1), en el que es preferible la correspondiente al modelo fundamentado en el perfil de velocidad configurado en dos zonas. Se recomienda restringir la aplicación de la ecuación ajustada con arreglo a la ley potencial al intervalo de y/d90 comprendido entre uno y veinte, puesto que fuera de dicho intervalo tiende a infraestimar notablemente la resistencia al flujo.