A mixed finite element method for nonlinear diffusion equations

We propose a mixed finite element method for a class of nonlinear diffusion equations, which is based on their interpretation as gradient flows in optimal transportation metrics. We introduce an appropriate linearization of the optimal transport problem, which leads to a mixed symmetric formulation. This formulation preserves the maximum principle in case of the semi-discrete scheme as well as the fully discrete scheme for a certain class of problems. In addition solutions of the mixed formulation maintain exponential convergence in the relative entropy towards the steady state in case of a nonlinear Fokker-Planck equation with uniformly convex potential. We demonstrate the behavior of the proposed scheme with 2D simulations of the porous medium equations and blow-up questions in the Patlak-Keller-Segel model.

Stable stationary states of non-local interaction equations

"Vegeu el resum a l'inici del document del fitxer adjunt."

Equations of hyperelliptic Shimura curves

We describe an algorithm that computes explicit models of hyperelliptic Shimura curves attached to an indefinite quaternion algebra over Q and Atkin-Lehner quotients of them. It exploits Cerednik-Drinfeld’s nonarchimedean uniformisation of Shimura curves, a formula of Gross and Zagier for the endomorphism ring of Heegner points over Artinian rings and the connection between Ribet’s bimodules and the specialization of Heegner points, as introduced in [21]. As an application, we provide a list of equations of Shimura curves and quotients of them obtained by our algorithm that had been conjectured by Kurihara.

Stabilized Petrov-Galerkin methods for the convection-diffusion-reaction and the Helmholtz equations

We present two new stabilized high-resolution numerical methods for the convection–diffusion–reaction (CDR) and the Helmholtz equations respectively. The work embarks upon a priori analysis of some consistency recovery procedures for some stabilization methods belonging to the Petrov–Galerkin framework. It was found that the use of some standard practices (e.g. M-Matrices theory) for the design of essentially non-oscillatory numerical methods is not feasible when consistency recovery methods are employed. Hence, with respect to convective stabilization, such recovery methods are not preferred. Next, we present the design of a high-resolution Petrov–Galerkin (HRPG) method for the 1D CDR problem. The problem is studied from a fresh point of view, including practical implications on the formulation of the maximum principle, M-Matrices theory, monotonicity and total variation diminishing (TVD) finite volume schemes. The current method is next in line to earlier methods that may be viewed as an upwinding plus a discontinuity-capturing operator. Finally, some remarks are made on the extension of the HRPG method to multidimensions. Next, we present a new numerical scheme for the Helmholtz equation resulting in quasi-exact solutions. The focus is on the approximation of the solution to the Helmholtz equation in the interior of the domain using compact stencils. Piecewise linear/bilinear polynomial interpolation are considered on a structured mesh/grid. The only a priori requirement is to provide a mesh/grid resolution of at least eight elements per wavelength. No stabilization parameters are involved in the definition of the scheme. The scheme consists of taking the average of the equation stencils obtained by the standard Galerkin finite element method and the classical finite difference method. Dispersion analysis in 1D and 2D illustrate the quasi-exact properties of this scheme. Finally, some remarks are made on the extension of the scheme to unstructured meshes by designing a method within the Petrov–Galerkin framework.

A subspace correction method for discontinuous Galerkin discretizations of linear elasticity equations

We study preconditioning techniques for discontinuous Galerkin discretizations of isotropic linear elasticity problems in primal (displacement) formulation. We propose subspace correction methods based on a splitting of the vector valued piecewise linear discontinuous finite element space, that are optimal with respect to the mesh size and the Lamé parameters. The pure displacement, the mixed and the traction free problems are discussed in detail. We present a convergence analysis of the proposed preconditioners and include numerical examples that validate the theory and assess the performance of the preconditioners.

Hyperbolic reaction-diffusion equations for a forest fire model

Forest fire models have been widely studied from the context of self-organized criticality and from the ecological properties of the forest and combustion. On the other hand, reaction-diffusion equations have interesting applications in biology and physics. We propose here a model for fire propagation in a forest by using hyperbolic reaction-diffusion equations. The dynamical and thermodynamical aspects of the model are analyzed in detail

Speed of wave-front solutions to hyperbolic reaction-diffusion equations

The asymptotic speed problem of front solutions to hyperbolic reaction-diffusion (HRD) equations is studied in detail. We perform linear and variational analyses to obtain bounds for the speed. In contrast to what has been done in previous work, here we derive upper bounds in addition to lower ones in such a way that we can obtain improved bounds. For some functions it is possible to determine the speed without any uncertainty. This is also achieved for some systems of HRD (i.e., time-delayed Lotka-Volterra) equations that take into account the interaction among different species. An analytical analysis is performed for several systems of biological interest, and we find good agreement with the results of numerical simulations as well as with available observations for a system discussed recently

High order hydrodynamical equations

Hydrodynamical equations act as a link between the local observed magnitudes of galactic motion and the general ones accounting for the behaviour of the Galaxy as a whole. Constraints are set usually in order to use them even in the lower order hierarchy. The authors present in this paper the complete expressions up to their fourth order. These equations will be used in the next future in their general form taking into account both the expected increase of kinematic data that the astrometric mission Hipparcos will provide, and some recent results indicating the possibility to obtain estimates for the momenta gradients.

Influence of the (111) twinning on the formation of diamond cubic/diamond hexagonal heterostructures in Cu-catalyzed Si nanowires

The occurrence of heterostructures of cubic silicon/hexagonal silicon as disks defined along the nanowire (111) growth direction is reviewed in detail for Si nanowires obtained using Cu as catalyst. Detailed measurements on the structural properties of both semiconductor phases and their interface are presented. We observe that during growth, lamellar twinning on the cubic phase along the (111) direction is generated. Consecutive presence of twins along the (111) growth direction was found to be correlated with the origin of the local formation of the hexagonal Si segments along the nanowires, which define quantum wells of hexagonal Si diamond. Finally, we evaluate and comment on the consequences of the twins and wurtzite in the final electronic properties of the wires with the help of the predicted energy band diagram.

Differential equations driven by fractional Brownian motion

A global existence and uniqueness result of the solution for multidimensional, time dependent, stochastic differential equations driven by a fractional Brownian motion with Hurst parameter H> is proved. It is shown, also, that the solution has finite moments. The result is based on a deterministic existence and uniqueness theorem whose proof uses a contraction principle and a priori estimates.

Relaxation times of non-Markovian processes

We consider a general class of non-Markovian processes defined by stochastic differential equations with Ornstein-Uhlenbeck noise. We present a general formalism to evaluate relaxation times associated with correlation functions in the steady state. This formalism is a generalization of a previous approach for Markovian processes. The theoretical results are shown to be in satisfactory agreement both with experimental data for a cubic bistable system and also with a computer simulation of the Stratonovich model. We comment on the dynamical role of the non-Markovianicity in different situations.

Fluctuations in domain growth: Ginzburg-Landau equations with multiplicative noise

Ginzburg-Landau equations with multiplicative noise are considered, to study the effects of fluctuations in domain growth. The equations are derived from a coarse-grained methodology and expressions for the resulting concentration-dependent diffusion coefficients are proposed. The multiplicative noise gives contributions to the Cahn-Hilliard linear-stability analysis. In particular, it introduces a delay in the domain-growth dynamics.

Numerical algorithm for Ginzburg-Landau equations with multiplicative noise: Application to domain growth

We consider stochastic partial differential equations with multiplicative noise. We derive an algorithm for the computer simulation of these equations. The algorithm is applied to study domain growth of a model with a conserved order parameter. The numerical results corroborate previous analytical predictions obtained by linear analysis.

Comment on "Reappraisal of experimental values of third-order elastic constants of some cubic semiconductors and metals"

In a recent paper A. S. Johal and D. J. Dunstan [Phys. Rev. B 73, 024106 (2006)] have applied multivariate linear regression analysis to the published data of the change in ultrasonic velocity with applied stress. The aim is to obtain the best estimates for the third-order elastic constants in cubic materials. From such an analysis they conclude that uniaxial stress data on metals turns out to be nearly useless by itself. The purpose of this comment is to point out that by a proper analysis of uniaxial stress data it is possible to obtain reliable values of third-order elastic constants in cubic metals and alloys. Cu-based shape memory alloys are used as an illustrative example.